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

(Difference between revisions)
Jump to: navigation, search
(Fixed mistake in pressure as pointed out by Kevin Schaal.)
 
(3 intermediate revisions not shown)
Line 1: Line 1:
-
The test case involves convection of an isentropic vortex in inviscid flow.
+
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  
-
The free-stream conditions are  
+
:<math>
:<math>
\begin{matrix}
\begin{matrix}
\rho &=& 1 \\
\rho &=& 1 \\
-
u &=& 0.5\\
+
u &=& u_\infty \\
-
v &=& 0\\
+
v &=& v_\infty \\
-
p &=& 1/\gamma
+
p &=& 1
\end{matrix}
\end{matrix}
</math>
</math>
-
Perturbations are added to the free-stream in such a way that there is no
+
Perturbations are added to the [[free-stream]] in such a way that there is no
-
entropy gradient in the flow-field. The perturbations are given by
+
[[entropy]] gradient in the [[flow-field]]. The perturbations in the velocity are given by:
-
:<math>
+
<math>
 +
(\delta u, \delta v)  = \frac{\beta}{2\pi} \exp\left( \frac{1-r^2}{2}
 +
\right) [ -(y-y_o), (x-x_o) ]
 +
</math>
 +
 
 +
while the temperature, density and pressure are given by
 +
 
 +
<math>
\begin{matrix}
\begin{matrix}
-
(\delta u, \delta v) &=& \frac{\beta}{2\pi} \exp\left( \frac{1-r^2}{2}
+
T &=& 1 - \frac{ (\gamma-1)\beta^2}{8\gamma\pi^2} \exp\left(
-
\right) [ -(y-y_o), (x-x_o) ] \\
+
1-r^2\right)\\
-
\rho &=& \left[ 1 - \frac{ (\gamma-1)\beta^2}{8\gamma\pi} \exp\left(
+
\rho &=& T^{\frac{1}{\gamma-1}} \\
-
1-r^2\right) \right]^{\frac{1}{\gamma-1}} \\
+
p &=& \rho T = \rho^\gamma
-
p &=& \frac{ \rho^\gamma }{\gamma}
+
\end{matrix}
\end{matrix}
</math>
</math>
Line 30: Line 35:
</math>
</math>
-
is distance from the vortex center <math>(x_o, y_o)</math>. One choice for the domain
+
is distance from the [[vortex]] center <math>(x_o, y_o)</math>.  
-
and parameters are
+
 
 +
One choice for the domain and parameters is:
:<math>
:<math>
Line 39: Line 45:
</math>
</math>
-
As a result of isentropy, the exact solution corresponds to a pure advection
+
As a result of [[isentropy]], the exact solution corresponds to a pure [[advection]]
-
of the vortex at the free-stream velocity. Further details can be found in Yee et al. (1999).
+
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==
==References==
*{{reference-paper | author=Yee, H-C., Sandham, N. and Djomehri, M., | year=1999 | title=Low dissipative high order shock-capturing methods using characteristic-based filters| rest=JCP, Vol. 150}}
*{{reference-paper | author=Yee, H-C., Sandham, N. and Djomehri, M., | year=1999 | title=Low dissipative high order shock-capturing methods using characteristic-based filters| rest=JCP, Vol. 150}}
 +
 +
 +
{{stub}}

Latest revision as of 04:00, 21 October 2014

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.


This article is a stub, a short article which needs to be improved. You can help by expanding it.
Retrieved from "https://openfoam.cfd-online.com/Wiki/2-D_vortex_in_isentropic_flow"
My wiki