CFD Online Logo CFD Online URL
Home > Wiki > Gresho vortex

Gresho vortex

From CFD-Wiki

Revision as of 09:56, 30 August 2013 by HannesGS (Talk | contribs)
Jump to: navigation, search


Classical Definition

The Gresho vortex was originally designed for incompressible flow, but here we consider the case of inviscid flow (Euler equations).

The Gresho problem [1] is a rotating vortex problem independent of time. Angular velocity u_{\phi} depends only on radius and the centrifugal force is balanced by the gradient of the pressure p.

\left( u_{\phi}(r), p(r) \right) = 
  \begin{cases} \left( 5r  , 5 + \frac{25}{2} r^2 \right), & 0   \leq r < 0.2, \\
                \left( 2-5r, 9 - 4 \ln 0.2 + \frac{25}{2} r^2 
                             - 20 r + 4 \ln r \right),     & 0.2 \leq r < 0.4, \\
                \left( 0   , 3 + 4 \ln 2 \right), & 0.4 \leq r.

The radial velocity is zero and the density is one everywhere.

Furthermore, we can modify this setup and let the vortex move with constant velocity. Customarily, we use periodic boundary conditions in each direction.

Dependence on Mach number

Miczek [2] modified the setup by introducing a reference Mach number. Pressure is scaled such that the rotation acts on this Mach number.

p_0 = \frac{\rho}{\gamma {\rm Ma}^2}, \left( u_{\phi}(r), p(r) \right) = 
  \begin{cases} \left( 5r  , p_0 + \frac{25}{2} r^2 \right),      & 0   \leq r < 0.2, \\
                \left( 2-5r, p_0 + \frac{25}{2} r^2 
                             + 4 (1 - 5r -\ln 0.2 +\ln r \right), & 0.2 \leq r < 0.4, \\
                \left( 0   , p_0 -2 + 4 \ln 2 \right), & 0.4 \leq r.

Now, we can use the Gresho vortex as test case how well the numerical scheme performs in dependence of the Mach number.

Initial condition for the azimuthal velocity in the modified setup with {\rm Ma}=0.1. Here, the computational domain is [-1,1] \times [-1,1] and the Gresho vortex was stretched by a factor 2.

Properties and Applications

In the original formulation, the Gresho vortex is a test for the following properties [1,3]:

  • Accuracy of numerical scheme: the empirical order of accuracy is around 1.4 for this problem.
  • Preservation of symmetry: since the analytical solution is pure rotation, the vortex should not change its shape.
  • Conservation of angular momentum: the smoothing of the azimuthal velocity profile is a measure how well the code preserves angular momentum.

Furthermore, [3] for the case of Voronoi meshes and [4] for curvilinear grids use the Gresho vortex to check the errors arising from the chosen numerical grid.

The modified version from [2] can be used to check the ability of the numerical scheme to handle low Mach number flows [2,4,5]. For small Mach numbers, the numerical dissipation of many finite difference schemes will lead to strong deformations of the vortex.


[1] Liska R., Wendroff B. Comparison of Several Difference Schemes on 1D and 2D Test Problems for the Euler Equations. SIAM J. Sci. Comput., 25(3), 995–1017 (23 pages).

[2] Miczek F. Simulation of low Mach number astrophysical flows. München, Technische Universität München, Diss., 2013.

[3] Springel V. E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh. Mon. Not. R. Astron. Soc. 401, 791–851 (2010).

[4] Grimm-Strele H., Kupka F., Muthsam H. Curvilinear Grids for WENO Methods in Astrophysical Simulations. Submitted to Computer Physics Communications, 2013.

[5] Happenhofer N., Grimm-Strele H., Kupka F., Löw-Baselli B., Muthsam H. A low Mach number solver: Enhancing applicability. Journal of Computational Physics Volume 236, 1 March 2013, Pages 96-118.

My wiki