# Gresho vortex

(Difference between revisions)
 Revision as of 07:33, 29 August 2013 (view source)HannesGS (Talk | contribs) (Created page with "== Classical Definition == The Gresho vortex was originally designed for incompressible flow, but here we consider the case of Euler equations. The Gresho problem [1] i...")← Older edit Revision as of 09:39, 30 August 2013 (view source)HannesGS (Talk | contribs) mNewer edit → Line 14: Line 14: The radial velocity is zero and the density is one everywhere. 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. Line 31: Line 33: Now, we can use the Gresho vortex as test case how well the numerical scheme performs in dependence of the Mach number. Now, we can use the Gresho vortex as test case how well the numerical scheme performs in dependence of the Mach number. + + + == 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]. == References == == References == Line 40: Line 55: [http://www.dart-europe.eu/full.php?id=749060 Simulation of low Mach number astrophysical flows]. München, Technische Universität München, Diss., 2013. [http://www.dart-europe.eu/full.php?id=749060 Simulation of low Mach number astrophysical flows]. München, Technische Universität München, Diss., 2013. - [3] Happenhofer N., Grimm-Strele H., Kupka F., Löw-Baselli B., Muthsam H. + [3] Springel V. - [http://dx.doi.org/10.1016/j.jcp.2012.11.002 A low Mach number solver: Enhancing applicability]. Journal of Computational Physics Volume 236, 1 March 2013, Pages 96-118. + [10.1111/j.1365-2966.2009.15715.x 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. [4] Grimm-Strele H., Kupka F., Muthsam H. [http://arxiv.org/abs/1308.3066 Curvilinear Grids for WENO Methods in Astrophysical Simulations]. Submitted to Computer Physics Communications, 2013. [http://arxiv.org/abs/1308.3066 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. + [http://dx.doi.org/10.1016/j.jcp.2012.11.002 A low Mach number solver: Enhancing applicability]. Journal of Computational Physics Volume 236, 1 March 2013, Pages 96-118. {{Stub}} {{Stub}}

## Classical Definition

The Gresho vortex was originally designed for incompressible flow, but here we consider the case of 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. \end{cases}$

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. \end{cases}$

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

## 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].

## References

[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. [10.1111/j.1365-2966.2009.15715.x 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.