Moving expansion waves: Difference between revisions

From Flowpedia
Jump to navigation Jump to search
No edit summary
No edit summary
 
(One intermediate revision by the same user not shown)
Line 1: Line 1:
[[Category:Compressible flow]]
<!--
[[Category:Unsteady waves]]
-->[[Category:Compressible flow]]<!--
[[Category:Inviscid flow]]
-->[[Category:Unsteady waves]]<!--
<noinclude>
-->[[Category:Inviscid flow]]<!--
[[Category:Compressible flow:Topic]]
--><noinclude><!--
</noinclude>
-->[[Category:Compressible flow:Topic]]<!--
--></noinclude><!--


__TOC__
--><nomobile><!--
-->__TOC__<!--
--></nomobile><!--


--><noinclude><!--
-->{{#vardefine:secno|6}}<!--
-->{{#vardefine:eqno|134}}<!--
--></noinclude><!--
-->
==== Moving Expansion Waves ====
==== Moving Expansion Waves ====


Line 24: Line 32:
The expansion is propagating into stagnant fluid in region four (the driver section), which means that the flow properties ahead of the expansion wave are constant.
The expansion is propagating into stagnant fluid in region four (the driver section), which means that the flow properties ahead of the expansion wave are constant.


<math display="block">
{{NumEqn|<math>
J^+_a=J^+_b
J^+_a=J^+_b
</math>
</math>}}


<math>J^+</math> invariants constant along <math>C^+</math> characteristics
<math>J^+</math> invariants constant along <math>C^+</math> characteristics


<math display="block">
{{NumEqn|<math>
J^+_a=J^+_c=J^+_e
J^+_a=J^+_c=J^+_e
</math>
</math>}}


<math display="block">
{{NumEqn|<math>
J^+_b=J^+_d=J^+_f
J^+_b=J^+_d=J^+_f
</math>
</math>}}


Since <math>J^+_a=J^+_b</math> this also implies <math>J^+_e=J^+_f</math>. In fact, since the flow properties ahead of the expansion are constant, all <math>C^+</math> lines will have the same <math>J^+</math> value.
Since <math>J^+_a=J^+_b</math> this also implies <math>J^+_e=J^+_f</math>. In fact, since the flow properties ahead of the expansion are constant, all <math>C^+</math> lines will have the same <math>J^+</math> value.
Line 42: Line 50:
<math>J^-</math> invariants constant along <math>C^-</math> characteristics
<math>J^-</math> invariants constant along <math>C^-</math> characteristics


<math display="block">
{{NumEqn|<math>
J^-_c=J^-_d
J^-_c=J^-_d
</math>
</math>}}


<math display="block">
{{NumEqn|<math>
J^-_e=J^-_f
J^-_e=J^-_f
</math>
</math>}}


<math display="block">
{{NumEqn|<math>
\left.
\left.
\begin{aligned}
\begin{aligned}
Line 59: Line 67:
\end{aligned}
\end{aligned}
\right\}\Rightarrow u_e=u_f \Rightarrow a_e=a_f
\right\}\Rightarrow u_e=u_f \Rightarrow a_e=a_f
</math>
</math>}}


Due to the fact the <math>J^+</math> is constant in the entire expansion region, <math>u</math> and <math>a</math> will be constant along each <math>C^-</math> line.
Due to the fact the <math>J^+</math> is constant in the entire expansion region, <math>u</math> and <math>a</math> will be constant along each <math>C^-</math> line.
Line 65: Line 73:
The constant <math>J^+</math> value can be used to obtain relations for the variation of flow properties through the expansion region. Evaluation of the <math>J^+</math> invariant at any position within the expansion region should give the same value as in region 4.
The constant <math>J^+</math> value can be used to obtain relations for the variation of flow properties through the expansion region. Evaluation of the <math>J^+</math> invariant at any position within the expansion region should give the same value as in region 4.


<math display="block">
{{NumEqn|<math>
u+\frac{2a}{\gamma-1}=u_4+\frac{2a_4}{\gamma-1}=0+\frac{2a_4}{\gamma-1}
u+\frac{2a}{\gamma-1}=u_4+\frac{2a_4}{\gamma-1}=0+\frac{2a_4}{\gamma-1}
</math>
</math>}}


and thus
and thus


<math display="block">
{{NumEqn|<math>
\frac{a}{a_4}=1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)
\frac{a}{a_4}=1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)
</math>
</math>}}


Eqn. \ref{eq:expansion:a} and <math>a=\sqrt{\gamma RT}</math> gives
Eqn. \ref{eq:expansion:a} and <math>a=\sqrt{\gamma RT}</math> gives


<math display="block">
{{NumEqn|<math>
\frac{T}{T_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^2
\frac{T}{T_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^2
</math>
</math>}}


Using isentropic relations, we can get pressure ratio and density ratio
Using isentropic relations, we can get pressure ratio and density ratio


<math display="block">
{{NumEqn|<math>
\frac{p}{p_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2\gamma/(\gamma-1)}
\frac{p}{p_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2\gamma/(\gamma-1)}
</math>
</math>}}


<math display="block">
{{NumEqn|<math>
\frac{\rho}{\rho_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2/(\gamma-1)}
\frac{\rho}{\rho_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2/(\gamma-1)}
</math>
</math>}}

Latest revision as of 13:36, 1 April 2026

Moving Expansion Waves

The expansion wave propagation into the driver section in a shock tube can be described using characteristic lines.


The expansion is propagating into stagnant fluid in region four (the driver section), which means that the flow properties ahead of the expansion wave are constant.

Ja+=Jb+(Eq. 6.135)

J+ invariants constant along C+ characteristics

Ja+=Jc+=Je+(Eq. 6.136)
Jb+=Jd+=Jf+(Eq. 6.137)

Since Ja+=Jb+ this also implies Je+=Jf+. In fact, since the flow properties ahead of the expansion are constant, all C+ lines will have the same J+ value.

J invariants constant along C characteristics

Jc=Jd(Eq. 6.138)
Je=Jf(Eq. 6.139)
ue=12(Je++Je)uf=12(Jf++Jf)Je=JfJe+=Jf+}ue=ufae=af(Eq. 6.140)

Due to the fact the J+ is constant in the entire expansion region, u and a will be constant along each C line.

The constant J+ value can be used to obtain relations for the variation of flow properties through the expansion region. Evaluation of the J+ invariant at any position within the expansion region should give the same value as in region 4.

u+2aγ1=u4+2a4γ1=0+2a4γ1(Eq. 6.141)

and thus

aa4=1γ12(ua4)(Eq. 6.142)

Eqn. \ref{eq:expansion:a} and a=γRT gives

TT4=[1γ12(ua4)]2(Eq. 6.143)

Using isentropic relations, we can get pressure ratio and density ratio

pp4=[1γ12(ua4)]2γ/(γ1)(Eq. 6.144)
ρρ4=[1γ12(ua4)]2/(γ1)(Eq. 6.145)