Physics:Simple wave

A simple wave is a flow in a region adjacent to a region of constant state. In the language of Riemann invariant, the simple wave can also be defined as the zone where one of the Riemann invariant is constant in the region of interest, and consequently, a simple wave zone is covered by arcs of characteristics that are straight lines.

Simple waves occur quite often in nature. There is a theorem (see Courant and Friedrichs) that states that a non-constant state of flow adjacent to a constant value is always a simple wave. All expansion fans including Prandtl–Meyer expansion fan are simple waves. Compressive waves until shock wave forms are also simple waves. Weak shocks (including sound waves) are also simple waves up to second-order approximation in the shock strength.

Simple waves are also defined by the behavior that all the characteristics under hodograph transformation collapses into a single curve. This means that the Jacobian involved in the hodographic transformation is zero.

Unsteady one-dimensional simple waves

Let [math]\displaystyle{ \rho }[/math] be the gas density, [math]\displaystyle{ u }[/math] the velocity, [math]\displaystyle{ p }[/math] the pressure and [math]\displaystyle{ c=\sqrt{(\partial p/\partial \rho)_s} }[/math] the speed of sound. In isentropic flows, entropy [math]\displaystyle{ s }[/math] is constant and if the initial state of the gas is homogenous, then entropy is a constant everywhere at all times and therefore the pressure is a function only of [math]\displaystyle{ \rho }[/math], i.e., [math]\displaystyle{ p=p(\rho) }[/math] In simple waves, all dependent variables are just function of any one of the dependent variables (this is certainly the case in one-dimensional sound waves) and therefore we can assume the velocity to be also a function only of [math]\displaystyle{ \rho }[/math]. i.e., [math]\displaystyle{ u=u(\rho). }[/math] This latter property is the cause of origin of the name simple wave, although the wave is nonlinear.

From the one-dimensional Euler equations, we have

[math]\displaystyle{ \frac{\partial \rho}{\partial t} + \frac{\partial (\rho u)}{\partial x} =0 }[/math]

[math]\displaystyle{ \frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + \frac{1}{\rho}\frac{\partial p}{\partial x}=0 }[/math]

which, because [math]\displaystyle{ u=u(\rho) }[/math], can be written as

[math]\displaystyle{ \frac{\partial \rho}{\partial t} + \frac{d(\rho u)}{d\rho}\frac{\partial \rho }{\partial x} =0 }[/math]

[math]\displaystyle{ \frac{\partial u}{\partial t} + \left(u+\frac{1}{\rho}\frac{dp}{du}\right) \frac{\partial u}{\partial x} =0. }[/math]

Further, since (remember that the time derivative of a function [math]\displaystyle{ f(x,t) }[/math] integrated along a curve [math]\displaystyle{ x=\varphi(t) }[/math] is given by [math]\displaystyle{ (df/dt)_\varphi = \partial f/\partial t + (dx/dt)_\varphi\partial f/\partial x }[/math])

[math]\displaystyle{ \frac{\partial \rho/\partial t}{\partial \rho/\partial x} = - \left(\frac{\partial x}{\partial t}\right)_\rho, \quad \frac{\partial u/\partial t}{\partial u/\partial x} = - \left(\frac{\partial x}{\partial t}\right)_u, }[/math]

the two equations lead to

[math]\displaystyle{ \left(\frac{\partial x}{\partial t}\right)_\rho= \frac{d(\rho u)}{d\rho}= u + \rho \frac{du}{d\rho}, \quad \left(\frac{\partial x}{\partial t}\right)_u = u + \frac{1}{\rho}\frac{dp}{du}. }[/math]

However, since [math]\displaystyle{ \rho }[/math] determines [math]\displaystyle{ u }[/math] and therefore the above derivatives must be equal so that [math]\displaystyle{ \rho du/d\rho=(1/\rho)dp/du=(c^2/\rho)d\rho/du }[/math]. Thus, we obtain [math]\displaystyle{ du/d\rho=\pm c/\rho }[/math], whence

[math]\displaystyle{ u = \pm \int \frac{c}{\rho}d\rho = \pm \int \frac{dp}{\rho c}. }[/math]

This equation provides the required relation [math]\displaystyle{ u=u(\rho) }[/math] or, [math]\displaystyle{ c=c(u) }[/math] or, [math]\displaystyle{ u=u(p) }[/math] etc. The above equation is just a statement that either the [math]\displaystyle{ J_+ }[/math] or the [math]\displaystyle{ J_- }[/math] Riemann invariant is constant.

Thus, we obtain

[math]\displaystyle{ \left(\frac{\partial x}{\partial t}\right)_u = u \pm c(u) }[/math],

which upon integration becomes

[math]\displaystyle{ x = t[u \pm c(u)] + f(u) }[/math]

where [math]\displaystyle{ f(u) }[/math] is an arbitrary function. This equation indicates that the characteristics in the [math]\displaystyle{ x }[/math]-[math]\displaystyle{ t }[/math] plane are just straight lines. When [math]\displaystyle{ f(u)=0 }[/math] and when consequently length scale and time scale associated with the initial function disappears, the problem is self-similar and the solution depends only on the ratio [math]\displaystyle{ x/t }[/math]. This particular case is referred as the centred simple wave.

Steady two-dimensional simple waves

Similar to the unsteady one-dimensional waves, simple waves in steady two-dimensional system cab be derived. In this case, the solution is given by

[math]\displaystyle{ y = xf_1(p) + f_2(p) }[/math]

where [math]\displaystyle{ f_1(p) = (\partial y/\partial x)_p }[/math] and [math]\displaystyle{ f_2(p) }[/math] is an arbitrary function of pressure. The characteristics in the [math]\displaystyle{ x }[/math]-[math]\displaystyle{ y }[/math] plane are straight lines. Similarly, the case corresponding to [math]\displaystyle{ f_2(p)=0 }[/math] is referred as the centred simple wave; the Prandtl–Meyer expansion fan is a special case of this centred wave.