Laminar CFD Solver programmed in C++

The wake of an incompressible flow over a square cylinder is numerically studied for Re = 50, 100 and s* = 0.5, 1. Distinct modes are observed for each individual flow configuration. A steady wake is observed for Re = 50 and s* = 1. However, on decreasing s* to 0.5 without altering Re, the wake sets into a periodically inphase oscillation. While increasing Re to 100 and keeping s* constant at 1, an inantiphase mode is observed where the flow structures even though being unsteady remain symmetrical about the centerline splitter plate. Therefore, the splitter plate experiences minimal lift but in a stochastic manner after 400 nondimensional time steps. [2]


References:
[1] Patankar, S. V. and Spalding, D.B., "A calculation procedure for heat, mass and momentum transfer in threedimensional parabolic flows", Int. J. of Heat and Mass Transfer, Volume 15, Issue 10, October 1972, Pages 17871806.
[2] S. Anantharamu and Gaonkar, V. A. and Kadoli, R., Effect of Centerline Splitter Plate on Flow over a Pair of SideBySide Square Cylinders: A Numerical Study, 4th International Conference on Fluid Mechanics and Fluid Power, Indian Institute of Technology Kanpur, India, 2014.
[1] Patankar, S. V. and Spalding, D.B., "A calculation procedure for heat, mass and momentum transfer in threedimensional parabolic flows", Int. J. of Heat and Mass Transfer, Volume 15, Issue 10, October 1972, Pages 17871806.
[2] S. Anantharamu and Gaonkar, V. A. and Kadoli, R., Effect of Centerline Splitter Plate on Flow over a Pair of SideBySide Square Cylinders: A Numerical Study, 4th International Conference on Fluid Mechanics and Fluid Power, Indian Institute of Technology Kanpur, India, 2014.
Turbulent CFD solver programmed on MATLAB
Wind Farm, Actuator Disc, ZeroOrder Turbulence Model
While working as a student research assistant with Stuttgart Wind Energy (SWE), I developed a transient two dimensional Navier Stokes solver. The solver was mainly used to calculate the flow field of a wind farm in which the position and size of the wind turbines were given as a user input. A structured mesh was generated as a collection of geometric progressions with simple but powerful controls given to the user. This ensured that the gradients near the turbine could be resolved at the accuracy necessary. This would especially come in handy because the actuator disc theory was used in this solver to calculate the forces exerted by the wind turbines and it is in the proximity of the wind turbine that the gradients of velocity are largest. The upwind and central differencing scheme were hybridised for the momentum fluxes to keep the solver stable but still predict accurate values. The transient solver was solved using a quasi Linear Parameter Varying model mainly for speed reasons. Soon after, a zero equation turbulence model was added to the solver based on the mixing length model.
References:
[1] Versteeg, H. K. and Malalasekera, W., An introduction to computational fluid dynamics, Pearson Education Ltd. Essex England, 1995
References:
[1] Versteeg, H. K. and Malalasekera, W., An introduction to computational fluid dynamics, Pearson Education Ltd. Essex England, 1995