# Cfd simulation of vortex shedding

### From CFD-Wiki

(improved writing) |
(improved linking) |
||

Line 1: | Line 1: | ||

- | To numerically simulate [[ | + | To numerically simulate [[vortex shedding]], [[CFD]] is used to calculate the unsteady flow that arises from a fluid moving past an obstruction. [[Downwind]] of the obstruction are regions of lower pressure which causes the fluid which initially was deflected around the obstruction to get sucked into these regions, begin to circulate, and form [[vortices]]. |

- | For a [[Newtonian | + | For a [[Newtonian fluid]], it is sufficient to use the [[Navier-Stokes equations]] as the governing equations, and choose either the [[Finite element|finite element method]], '''FEM''', or [[Finite difference|finite differencing]], '''FD''', to do the simulation. Vortex shedding is more likely in [[Incompressible flow|incompressible fluids]] than [[Compressible flow|compressible fluids]], so the former is used when capturing this phenomenon. Although both types of fluids are represented by the [[Navier-Stokes equations]], their methods of solution differ significantly. For instance, when using '''FD''' for incompressible fluids, it is common to have an explicit expression for the velocity as you integrate the equations in time. In [[Compressible flow|compressible fluids]], the governing equations are often written in terms of the [[conservative form]] of [[Euler's equation]], with velocity and pressure extracted from the conserved quantities(e.g. mass, momentum, energy). Also, compressible flows usually involve the calculation of additional quantities such as density and temperature. |

==Example of 2 Dimensional Problem== | ==Example of 2 Dimensional Problem== | ||

Line 49: | Line 49: | ||

==Pressure Corrector== | ==Pressure Corrector== | ||

- | The incompressibility of the fluid is enforced through the [[continuity equation]] using pressure as a corrective quantity whose gradient in either the x or y direction is added to the velocity in the respective direction. Calculating the pressure correction involves solving the [[Discrete Poisson equation|discrete Poisson equation]] for pressure. | + | The incompressibility of the fluid is enforced through the [[Continuity_Equation|continuity equation]] using pressure as a corrective quantity whose gradient in either the x or y direction is added to the velocity in the respective direction. Calculating the pressure correction involves solving the [[Discrete Poisson equation|discrete Poisson equation]] for pressure. |

==Staggered grid== | ==Staggered grid== |

## Revision as of 15:56, 24 May 2007

To numerically simulate vortex shedding, CFD is used to calculate the unsteady flow that arises from a fluid moving past an obstruction. Downwind of the obstruction are regions of lower pressure which causes the fluid which initially was deflected around the obstruction to get sucked into these regions, begin to circulate, and form vortices.

For a Newtonian fluid, it is sufficient to use the Navier-Stokes equations as the governing equations, and choose either the finite element method, **FEM**, or finite differencing, **FD**, to do the simulation. Vortex shedding is more likely in incompressible fluids than compressible fluids, so the former is used when capturing this phenomenon. Although both types of fluids are represented by the Navier-Stokes equations, their methods of solution differ significantly. For instance, when using **FD** for incompressible fluids, it is common to have an explicit expression for the velocity as you integrate the equations in time. In compressible fluids, the governing equations are often written in terms of the conservative form of Euler's equation, with velocity and pressure extracted from the conserved quantities(e.g. mass, momentum, energy). Also, compressible flows usually involve the calculation of additional quantities such as density and temperature.

## Contents |

## Example of 2 Dimensional Problem

The images below were created for the case of 2-D flow past a square obstruction with a Reynolds number of 250 and freestream boundary conditions on the top, bottom, and right boundaries. The quantities of the Navier-Stokes equations were non-dimensionalized and the QUICKEST finite differencing algorithm by B.P. Leonard was used. The grid was uniform with 204 points in the x direction and 90 points in y.

## Pressure Corrector

The incompressibility of the fluid is enforced through the continuity equation using pressure as a corrective quantity whose gradient in either the x or y direction is added to the velocity in the respective direction. Calculating the pressure correction involves solving the discrete Poisson equation for pressure.

## Staggered grid

Use of a staggered grid, where pressure and velocities in the x and y directions are calculated at alternating grid points, is commonly used with upwinding to avoid oscillatory behavior.

## See also

## References

- Davis, R. W., and E. F. Moore,
*A numerical study of vortex shedding from rectangles*, Journal of Fluid Mechanics, 116, 475 (1982).

- Davis, R. W., E. F. Moore, and L. P. Purtell,
*A numerical-experimental study of confined flow around rectangular cylinders*, Physics of Fluids, 27, 46 (1984).

- Hoffman, Joe D.,
*Numerical Methods for Engineers and Scientist, fourth edition*, McGraw-Hill Inc., New York, 1992.

- Leonard, B. P.,
*A stable and accurate convective modelling procedure*, Computer Methods in Applied Mechanics and Engineering, 19, pages 59-98 (1979).