This is supposed to help with what might happen in atmospheric uplift, as in a hurricane. The pressure reduces and the gas expands in the same way. If you focus on the moving air, in fluid mechanics that is the Lagrangian view.
It's better for intuition - the boundaries of the volume travel with the fluid, and by assumption in the atmosphere, heat and other species like water vapor travel with it too, and remain in the boundaries.. That is because diffusion over the large distances is relatively small, although thermal radiation could be an issue.
But there are problems when you want to think about the processes in relation to other things that are not moving with that fluid. Fixed boundaries are one issue - not so common in the atmosphere. But in the ongoing argument about Anastassia Makarieva's contraction theories (TaV1, JC and TaV2, Lucia1, Lucia2), AM prefers the Eulerian frame, a fixed volume in space.
The good thing there is that you can assume steady conditions. At every point the velocity, temperature etc will remain steady, even though as you follow it in the flow, it changes. The bad thing is that the boundaries are no longer material boundaries, and you have to account for whatever is crossing it. Not so intuitive.
Scientists tend to talk Lagrangian to students, but Eulerian to their computers.
Here I'll take the Eulerian view of expansion and condensation.
A 1D (vertical) problem
In accounting for boundaries, it helps to have only two (upper and lower). I'll look at a 1 m thick slice, assuming zero horizontal gradients. The notation is generally as in the previous post. We'll say z goes from 0 to 1, and I'll suffix accordingly (but if suffix is missing, assume zero). So T,V,P, ρ are assumed known at z=0. The incoming air is saturated.
I'll also assume, as is conventional, hydrostatic pressure. So P1 =P0 -ρ g
Then each incoming packet of air changes according to the cylinder analysis, with this pressure increment:
T1 =T0 -ρ g dT/dPwhere dT/dP = (T/P) DT/DP = (T/P) (1+c+c/b)/(3.5+b+b*c), b=L/(R*T), c = b*nw/n
In the previous post, n, nw were moles in a volume V. Here I'll take V=1, so they are molar densities.
nw1 =nw0 -ρ g dnw/dPwhere dnw/dP = (nw/P) Dnw/DP = (nw/P) (1+c+c/b)/(3.5+b+b*c)
Density is calculated from dV. The mass change in the unit volume element Mw*dnw, ( Mw = molar mass of water) so the density change is
d ρ /dP = Mw*dnw /dP + ρ/P DV/DP= (Mw*nw /P *(b-3.5)+ ρ/P (2.5+b*c+2.5*c/b-c))/(3.5+c+b*c)
Velocity and precipitation reconciliationThe molar mass flux of dry air through z0 is v0*Nd
The number of moles is conserved; in the notation of the previous post, Nd varies inversely as V, ie DNd/DP=-DV/DP.
So with v=1 m/s,
dv/dP=(v/P)Dv/DP = (v/P)DV/DP=(v/P)a4/a3
Velocity v increases with z.
The precipitation rate r moles/m3/s is
r=v dnw/dz= vρ g (nw/P) (1+c+c/b)/(3.5+b+b*c)
ResultsAgain the meaningful test is to compare saturated air with unsaturated air under the same conditions. With z in km:
InterpretationIt is a steady flow. But also expanding. The upward velocity would have increased just because of the hydrostatic density gradient, but it increases more because of the condensation. This is due to the release of latent heat, which dominates the contraction due to mass loss.
The same expansion is seen in the density gradient (reducing faster). and the warming is seen in the well-known reduction of lapse rate with condensation.
Further noteI mentioned in the previous post that the formula for DT/DP could be found, with slight modification, in the Phys Met notes of Dr Caballero. I realised then that it is just the formula for the saturated adiabatic lapse rate, which has a much longer history.
As explained in that earlier post, my form isn't quite the same. There is a discrepancy due to the mole concentration I made in the Ideal Gas Law. I've kept with it, because it does seem to be the irght thing to do.
T=298.15 V=v=1 P=98131 R=8.314 N=P*V/(R*T) Nw=0.03128*N Cv=20.85 Mw=0.018 Ma=0.02896 L=2304900*Mw g=9.8 rho = N*Ma b=L/R/T c=b*Nw/N a1=1+c+c/b a2=b-3.5 a3=b*a1-a2 a4= -(2.5*a1+c*a2) dp=c(0,0,0,0) dp=dTdP=(T/P)*(a1/a3) dp=dNwdP=(Nw/P)*(a2/a3) dp=drhodP = (Mw*(Nw/P)*a2 - (rho/P)*a4)/a3 dp=(v/P)*a4/a3 dz=-dp*rho*g*1000 # print(dp); print(c(a1,a2,a3,a4)) print(dz)