VORTEX FLOWS
   THERMAL FLOWS
   NUMERICAL METHODS
   ACOUSTICS
   AERODYNAMICS
   WAVES
    GRANULAR
   ROTATING FLOWS
NUMERICAL METHODS
 


A control volume formulation and an iterative semi-implicit scheme on uniform or non-uniform grids, is used to solve the time varying Navier Stokes equations in two and three dimensional domains. The solver developed at our Lab can be coupled with the energy equation to perform natural convection simulations.

 2D flows over arbitrary bodies
 
Governing equations are solved with our code in primitive variables using the finite volume formulation
SIMPLER algorithm. Discretization of the physical domain incorporates uniform and staggered grids where the number of grid points is a choice that depends on the spatial resolution, looking for both grid independent solutions and enough resolution for the eventual unstable wavy behavior of developing instabilities .

The use of the block-off method to emulate solid boundaries inside the numerical domain is a simple way to solve the flow equations with arbitrary bluff bodies like in the figure above, where we show the axial flow velocity colormap around our lab logo.
 
 Coupled Benard von Karman Instability
 
 

The BvK instability can be generated in the wake of a single circular cylinder (diameter d and length L) or with an array of them. Our goal here is to understand the effects of an open loop control stratregy on the vortex shedding for three side-by-side cylinders. Vortex shedding will start when the Reynolds number is greater than the critical value creating a wake flow of three coupled von Karman streets.

For the case of two side-by-side cylinders, at low Reynolds numbers, both a symmetric and an asymmetric shedding mode of two parallel vortex streets can be observed. If we include a third cylinder in between two peripheral cylinders, we can fix the starting phase of all of three shedding modes using a numerical version of the wake control method by suction of the stagnation point at each cylinder.

A proper control of the shedding phase of will produce an array of vortices forming a kind of lattice where the number of shed vortices (at the Strohual frequency fs ) is determined with the control method and the lattice size with the wavelength l = U/fs and cylinder separation. The three-body situation allow us to increase the number of vortices over a transverse length scale, and to get a wider lattice-like vortex distribution. It also helps to understand the role of overall flow symmetries associated to local broken flow symmetries for each wake.

  Universidad de Chile, Facultad de Ciencas Físicas y Matemáticas