# 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 Latest revision as of 09:25, 4 February 2014 (view source)HannesGS (Talk | contribs) (9 intermediate revisions not shown) Line 1: Line 1: == Classical Definition == == 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 for the case of [[inviscid flow]] (Euler equations). Angular velocity $u_{\phi}$ depends only on radius and the centrifugal force is balanced by the gradient of the pressure $p$. - + - 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) = [itex]\left( u_{\phi}(r), p(r) \right) = Line 14: Line 12: 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. == Dependence on Mach number == == 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. + Miczek [2] modified the setup by introducing a reference [[Mach number]]. Pressure is scaled such that the rotation acts on this Mach number. [itex]p_0 = \frac{\rho}{\gamma {\rm Ma}^2},$ $p_0 = \frac{\rho}{\gamma {\rm Ma}^2},$ Line 31: Line 31: 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. + [[File:Gresho.png|600px|thumb|left|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. == References == == References == Line 40: Line 53: [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. + [http://dx.doi.org/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://dx.doi.org/10.1016/j.cpc.2013.11.005 Curvilinear Grids for WENO Methods in Astrophysical Simulations]. Computer Physics Communications 184 (2014) 764-776. - + [5] Happenhofer N., Grimm-Strele H., Kupka F., Löw-Baselli B., Muthsam H. - {{Stub}} + [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.

## Classical Definition

The Gresho problem [1] is a rotating vortex problem independent of time for the case of inviscid flow (Euler equations). 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.

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.

## 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. 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. Computer Physics Communications 184 (2014) 764-776.

[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.