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
m (Provided interlinks between this page and other pages (many may require creation))
(Fixed mistake in pressure as pointed out by Kevin Schaal.)
 
(2 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 40: Line 46:
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