The time derivatives



The most important difference between those two methods that describes the fluid motion comes in the time rate of change. So in Lagrange notation the total derivative Df/Dt of some property f of the fluid is the time rate of change of f in a portion of the fluid, as it moves though the point x at time t. In contrast, in the Euler notation the partial derivative ∂f/∂t is the change of f at the fixed point x, as the fluid streams past. Thus ∂f/∂t includes not only the time rate of change of f in the moving fluid, but also the change in f at x because the fluid is moving past, if f differs from point to point in the fluid.
Both methods have some advantages. The main advantage of Lagrange description lies in the time derives. So total derivative dρ/dt measures fluid expansion or compression and du/dt measures true acceleration of a portion of the fluid, and thus is proportional to the net force acting on the fluid at x,t. Euler description of the fluid motion has many other advantages so it will be useful to find the measure of the difference between the two time derivatives.
To find the total derivative of some property of the fluid, in terms of the partial derivatives at t,x, we compare the value of f at point x at time t with its value at time t+dt at the point x+udt to which the fluid has moved in time dt. Accoridng to the definitions of derivatives, this is:
$$\begin{align} & \frac{df}{dt}dt=f(x+udt,t+dt)-f(x,t) \\ & \frac{df}{dt}dt=f(x,t)+\frac{\partial f}{\partial x}udt+\frac{\partial f}{\partial t}dt-f(x,t) \\ & \frac{df}{dt}=\frac{\partial f}{\partial t}+u\frac{\partial f}{\partial x} \\ \end{align}$$
 
The actual change in f, as the fluid moves past x with velocity u, is the change in f with time at the fixed point x, plus the spatial change in f times its velocity u of flow past x. Thus the true acceleration of the fluid at x at time t is given, in terms of the Eulerian velocity function u(x,t), by
$$\frac{du}{dt}=\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}$$
We see that the correction term u(∂u/∂x) is second order in the fluid velocity, and therefore, if u is small (as it often is for sound waves), du/dt = ∂u/∂t to the first order of approximation. The distinction cannot be forgotten, however, because the sound energy and intensity are second-order expressions, and the correction term may have to be included.

Fluid displacement and velocity



Our next task is to work out an approximate equation of motion for small-amplitude sound waves in a fluid. But before we do this, we must develop a notation which will adequately describe the motion. As a matter of fact, we shall develop two notations, one of which will be generally employed; but the other will be occasionally useful. The most straightforward way, the extension of the notation we have used for the chain of masses and springs, is to label each portion of the fluid  by giving its position at t = 0, and then to give the subsequent position of each portion as a function of time and of its initial position. In this notation, called the Lagrange description, we follow the fluid in its motion.


In one-dimensional motion, the portion of fluid initially at xo has position X at time t, where X is a function of t and of xo. All other quantities similarly follow the motion of the fluid. The density ρ(xo,t) is the density of the portion of the fluid which was initially at Xo and which is at X at time t. This notation gives us a good insight into what is going on inside the fluid as it moves about. For example, if the time rate of change dρ/dt of this density is less than zero, we can be sure the fluid is expanding.

The derivative du/dt, where u = dX/dt is the fluid velocity, is the acceleration of the element of fluid which was originally at xo, and is thus equal to the force acting on the element, divided by its mass. However, the notation also has its disadvantages, for the coordinate system moves with the fluid, and it is not easy to determine what the fluid is doing at some specified point in space at a given time.
So the alternative notation, the Euler description, sets up a coordinate system (x, in one dimension) fixed in space and describes the properties of the fluid which happens to be at a given point at a given time. Therefore, in this notation, the various quantities are functions of x and of t, for one dimension. These are quantities, such as density, velocity, etc., measured at a fixed point in space at a given time, and therefore corresponding to different portions of the fluid, as the fluid streams past the point.
This is the description we shall be using most of the time, with p and u given as explicit functions of x and t; by it we can tell what is happening where, without having to go back to find which part of the fluid happens to be there at time t. In particular, the Euler description enables us to compute easily the spatial variation of a given quantity at time t, the partial derivative ∂p/∂x being the rate of change of density at x, as might be measured from an instantaneous photograph of the density at time t.

Acoustic energy and momentum



Consider a long pipe of uniform cross section A, and having rigid walls. Fluid is confined to the interior of a long pipe. 

Figure 1 - Pressure wave generated by motion of piston D in pipe of
area A


So at one end is a piston of negligible mass, which constrains the fluid to stay to its right. The piston was at rest at point 0 held by the force F=AP counteracting the equilibrium pressure P of the fluid. Suppose that at time t=0, the piston is started into motion to the right with constant velocity u. This action will push the piston to the right increasing the pressure to P+p. If we want to get a piston into motion we need to increase the force from AP to A(P+p), where p is related to the change in density of the fluid by the equation δ=κρp where coefficient of compressibility can be isothermal or adiabatic.
Originally δ is related to the change in volume of the compressed fluid, and this is related to the amount of fluid which has been set in motion at any specified time t. Since fluid has it’s compressibility and inertia it takes time to get it into motion. So when the piston is moved from equilibrium position the starting wave front will travel away from the piston reaching the point B at time t. Let’s say that in time t the wave has traveled the distance ct where c is the velocity of the wave. Then all the fluid to the right of this front, at B is in its original equilibrium condition at rest with density ρ. The fluid in the region D to B in motion with velocity u of the piston compressed to a density and thus occupies a volume ρ/( ρ +δ) times its original volume Al=Act. The difference between this volume and the original volume Act is of course the volume Aut, swept out by the piston in time t so that
$$\begin{align} & Aut=Act-\frac{A\rho ct}{\rho +\delta } \\ & \frac{u}{c}=1-\frac{\rho }{\rho +\delta }=\frac{\rho +\delta -\rho }{\rho +\delta }=\frac{\delta }{\rho } \\ \end{align}$$

Where the last formula is valid only as long as the density change δ caused by the motion is small compared with the equilibrium density.  In traveling the distance ut the piston has done an amount of work A(P+p)ut which must have gone into additional energy of the fluid. Part of this energy will be kinetic energy one-half of the mass of the fluid in motion, 0.5*Aρct times the square of its velocity u². The other part is potential energy of compression the work done in compressing the fluid from volume Act to volume Act[ρ/(ρ+δ)]. The work is the integral of the pressure P+p times the change in volume dV, integrated between these limits. Thus the potential energy is:
$$\int{\left( P+p \right)V\kappa dp=PV}\kappa p+\frac{1}{2}V\kappa {{p}^{2}}\simeq Aut\left( P+\frac{u}{2\kappa c} \right)$$
Thus the equation for energy balance becomes:

Force x distance = potential + kinetic energy

$$Aut\left( P+p \right)=AutP+Aut\frac{u}{2\kappa c}+Aut\frac{2}{\rho cu}$$
The first terms on both sides cancel. In fact, the work APut should not be charged to energy of the sound wave, since it is the work done pushing against the equilibrium pressure; if there had been fluid pressure P on the lef hand surface of the piston, it would have provided the AP part of the force and only force Ap would have been required to move the piston.

Potential energy density = 0.5*κp²

Kinetic energy density = 0.5*ρ*u²

For the energies per unit volume of the sound, correct ot the second order in the small quantities p and u.

What we have shown by this manipulation is that the mass of the fluid, represented by the equilibrium density ρ, and its elasticity, represented by κ, so combine that it takes time for the motion of one part of the fluid to be transmitted to another part, the speed of transmission being equal to c=Sqrt(1/ρκ) as long as the fluid motion is small enough so that the excess pressure p is small compared with the equilibrium pressure P. The impulse-momentum relationship could also be used to compute c. The excess force Ap contributes an impulse Apt to the fluid in time t, which must equal the momentum contributed to the fluid newly set into motion in time t. The acoustic momentum density of the fluid is up, so that
$$Apt=\left( u\rho \right)\left( Act \right)\simeq A{{c}^{2}}\delta t$$
It is obvious that the wave velocity c will be independent of the amplitude and shape of the sound wave only as long as the sound pressure p is small compared with the equilibrium pressure P.
For fluids like water, where K is small enough so that δ/ρ is small even though piP is not small, another limitation enters. For example, if the piston were pushed to the left, a wave of rarefaction would be sent along the tube, and if u were large enough, the p computed from the formulas might be equal to or larger than P; the piston would be withdrawn rapidly enough to produce a vacuum behind it. When this happens the fluid loses contact with the piston, leaving a large bubble of vapor in the interspace, and it is obvious that our equations no longer hold. Long before this happens the value of K begins to change, so that K cannot be considered independent of
p unless Abs[p] is small compared with P.

Some general properties of sound



Fluids such as water or air have mass density and volume elasticity thus they can be approximated with multi degree of freedom systems for example chain of masses and springs. The elasticity or volume elasticity helps the fluid to resist being compressed and when the external forces act on fluid, a fluid tends to return to its original state.
At equilibrium fluid has density ρ  [kg/m^2], has a uniform temperature T [K] and uniform pressure  P [Pa, atm, N/m^2] . These three quantities are connected by an equation of state, which can be given in two forms. These forms are explicit equation of state which connects these three quantities in explicit form  or  in terms of partial derivatives, assuming density as function of T and P.
            The equation of state is usually expressed in terms of volume V that occupied by n moles of the fluid. If nM kg of the fluid occupies volume V the density of the fluid is obviously ρ=nM/V. The two partial derivatives which are usually employed to measure the equation of state of the fluid near the equilibrium state are:
$$\begin{align} & {{\kappa }_{T}}=-\frac{1}{V}{{\left( \frac{\partial V}{\partial P} \right)}_{T}}=\frac{1}{\rho }{{\left( \frac{\partial \rho }{\partial P} \right)}_{T}} \\ & \beta =\frac{1}{V}{{\left( \frac{\partial V}{\partial T} \right)}_{P}}=-\frac{1}{\rho }{{\left( \frac{\partial \rho }{\partial T} \right)}_{P}} \\ \end{align}$$

Quantity κt the fraction rate of change of volume or  density with pressure at constant temperature is called isothermal compressibility of the fluid.
Quantity β is the fractional change in volume with temperature at constant pressure is called the coefficient of thermal expansion of the fluid.
Both κt and β are functions of P and T though usually their rates of change are not very large. From κt and β we can compute the other partial derivatives of ρ, P and T near the equilibrium state of the fluid.
If the P, β and κt are expressed as functions of ρ and T, then
$${{\left( \frac{\partial P}{\partial T} \right)}_{\rho }}=\frac{\beta }{{{\kappa }_{T}}}$$
            Kinetic theory indicates that the pressure of a perfect gas is caused by its molecular motion and that P=1/3 * ρ where is the mean square molecular velocity and ro is the gas density. The compressibility of a perfect gas is thus inversely proportional to the mean kinetic energy of the gas molecules. The velocity propagation of the spring-mass system is:
$$v=\sqrt{\frac{1}{\kappa \varepsilon }}$$
Where k is the compressibility of the spring, and ε is the mass per unit length of the chain. In acoustics the velocity of propagation of sound waves in a fluid can be determined from:
$$c=\sqrt{\frac{1}{\kappa \rho }}$$
Where k is the compressibility of the air and ρ is density. For perfect or ideal gases
$${{\kappa }_{T}}=\frac{1}{P}=\frac{3}{\rho \left\langle {{v}^{2}} \right\rangle }$$
And the speed of sound is:
$$c=\sqrt{\frac{1}{{{\kappa }_{T}}\rho }}\simeq \sqrt{\frac{1}{\frac{3}{\rho \left\langle {{v}^{2}} \right\rangle }\rho }}=\sqrt{\frac{\left\langle {{v}^{2}} \right\rangle }{3}}$$
So we can see from previous equation that the speed of sound in perfect or ideal gas is approximately equal to the root mean square of the molecular speed.
In liquids and solid the speed of sound is much larger than the root mean square speed of vibration of an atom about its equilibrium position. The c or speed of sound in solids is roughly equal to the mean speed a molecule would have if its amplitude of vibration were equal to the mean distance between molecules. In reality the actual amplitudes are much smaller than the mean distance between molecules so the speed of sound is not equal to a mean speed of a molecule. In case that the mean speed of molecule is roughly equal to the speed of sound the solids would melt or vaporize.
Presence of a sound wave causes the change in density, pressure and temperature in the fluid and each change is being proportional to the amplitude of the wave. The pressure changes are most easily measurable, though some sound detectors measure the motion of the fluid, caused by the sound. You can also measure the changes in temperature of density but the results are not particularly accurate; so they are usually computed from the measured pressure change.

In fluids with very large thermal conductivity the temperature of the fluid is practically unchanged by the passage of the sound wave. In this case the isothermal compressibility kt is directly applicable. If the pressure change caused by the sound p and the equilibrium pressure is P, so that the total pressure in the presence of sound is P+p then the density of the fluid is:
$$\begin{align} & \rho +\delta =\rho +{{\left( \frac{\partial \rho }{\partial P} \right)}_{T}}p=\rho +\rho {{\kappa }_{T}}p \\ & \rho =\rho {{\kappa }_{T}}p \\ \end{align}$$
Where ρ is the equilibrium pressure and δ is the small increase in density produced by the sound.
If the frequency of sound is smaller than 10^9 cps or so in gases it is better approximation to assume that the compression is adiabatic than the entropy content of the gas remains constant during compression. In the case of adiabatic compression we need to formulate differential relations between density and pressure and between temperature and pressure and they are:
$$\begin{align} & {{\left( \frac{\partial \rho }{\partial P} \right)}_{s}}=\rho {{\kappa }_{s}}=\frac{1}{\gamma }{{\left( \frac{\partial \rho }{\partial P} \right)}_{T}} \\ & {{\kappa }_{s}}=\frac{{{\kappa }_{T}}}{\gamma } \\ & {{\left( \frac{\partial T}{\partial P} \right)}_{s}}=\frac{\gamma -1}{\gamma }{{\left( \frac{\partial T}{\partial P} \right)}_{\rho }}=\left( \gamma -1 \right)\frac{{{\kappa }_{s}}}{\beta } \\ & \rho +\delta =\rho +\rho {{\kappa }_{s}}p \\ & T+\tau =T+\left( \gamma -1 \right)\frac{{{\kappa }_{s}}}{\beta }p \\ \end{align}$$
γ is the ration between specific heat at constant pressure of the gas and the specific heat at constant volume. τ is the small change in temperature caused by the sound wave. For a perfect:
-          Monatomic gas (Helium) γ=5/3,
-          Diatomic gas (hydrogen or air) γ=7/5,
-          Polyatomic gases γ=4/3.
Just for remainder for the perfect gas the state equation is MP=RTρ.
$$\begin{align} & {{\kappa }_{s}}=\frac{1}{\gamma P} \\ & {{\left( \frac{\partial \rho }{\partial P} \right)}_{s}}=\frac{\rho }{\gamma P}=\frac{M}{\gamma RT} \\ & {{\left( \frac{\partial T}{\partial P} \right)}_{s}}=\frac{\gamma -1}{\gamma }\frac{T}{P}=\left( \gamma -1 \right)\frac{M}{\gamma R\rho } \\ \end{align}$$
If the sound is present in fluid than as we know the pressure is P+p, the density and the temperature are:
$$\begin{align} & \rho +\delta =\rho +\left( \frac{\rho }{\gamma P} \right)p \\ & \delta =\rho {{\kappa }_{s}}p=\frac{M}{\gamma RT}p \\ & T+\tau =T+\frac{\left( \gamma -1 \right)T}{\gamma P}p=T+\frac{\left( \gamma -1 \right)M}{\gamma R\rho }p \\ & \tau =\frac{\left( \gamma -1 \right)M}{\gamma R\rho }p \\ \end{align}$$
for adiabatic compression, for a perfect gas.
We can see that the ratio between the density change and the acoustic pressure p is independent of the equilibrium pressure but are inversely proportional to the equilibrium temperature of the gas. Also the ratio between the change in temperature and sound pressure is independent of temperature but inversely proportional to the equilibrium density. All these relationships are valid only to the first order in the small quantity p/P