CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > 2-D vortex in isentropic flow

2-D vortex in isentropic flow

From CFD-Wiki

Revision as of 04:00, 21 October 2014 by Praveen (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

This is one of the few exact solutions for the compressible Euler equations. The test case involves convection of an isentropic vortex in inviscid flow. It is used to show the ability of numerical schemes to accurately capture vortical flows. The free-stream conditions are


\begin{matrix}
\rho &=& 1 \\
u &=& u_\infty \\
v &=& v_\infty \\
p &=& 1
\end{matrix}

Perturbations are added to the free-stream in such a way that there is no entropy gradient in the flow-field. The perturbations in the velocity are given by:


(\delta u, \delta v)  = \frac{\beta}{2\pi} \exp\left( \frac{1-r^2}{2}
\right) [ -(y-y_o), (x-x_o) ]

while the temperature, density and pressure are given by


\begin{matrix}
T &=& 1 - \frac{ (\gamma-1)\beta^2}{8\gamma\pi^2} \exp\left(
1-r^2\right)\\
\rho &=& T^{\frac{1}{\gamma-1}} \\
p &=& \rho T = \rho^\gamma
\end{matrix}

where


r = [ (x-x_o)^2 + (y-y_o)^2 ]^{1/2}

is distance from the vortex center (x_o, y_o).

One choice for the domain and parameters is:


\Omega = [0,10] \times [-5,5], \quad
(x_o, y_o) = (5,0), \quad
\beta = 5

As a result of isentropy, the exact solution corresponds to a pure advection of the vortex at the free-stream velocity. This problem is usually solved with periodic boundary conditions. Further details can be found in Yee et al. (1999).

References

  • Yee, H-C., Sandham, N. and Djomehri, M., (1999), "Low dissipative high order shock-capturing methods using characteristic-based filters", JCP, Vol. 150.


My wiki