The last term is 1-k, because a lot of people find it confusing to see
P(choke)=P(up)(2/(k+1)^(k/(k-1))
P(up)=P(choke)/[(2/(k+1)^(k/(k-1))]
P(up)=P(choke)[(2/(k+1)^-(k/(k-1))]
So I took the minus sign into the expression (^-k/(k-1)=^+k/(1-k)
When I do the calc with air (k=1.4) at sea level (P(atm)=14.73 psia) I get P(up)=27.88 psia = 13.1 psig. Different gases and different elevations give you different numbers.
We all accept the term "choked flow" without asking "Why?". I did it for 20 years until a @$%@$% student asked me "Why is that?" and I had to think about it. Thinking about it for the last year has made my head ache, but I think I have some insight.
For approximately compressible flow (think velocity less than 0.6 Mach), gas sort of acts like an avalanche--the upstream gas "falls" into the downstream lower pressure. As the velocity approaches Mach 1.0, that stops happening as simply. Something limits the velocity to the speed of sound (unless you do something special that I'll get to in a minute). That "something" is a standing wave. The standing wave is a largely impermeable barrier to flow. Gas at Mach 1.0 can escape, but anything slower is trapped and builds up pressure. In that model, pressure upstream of the standing wave is at P(choked) and pressure downstream is at some lower pressure. So if you have a 10,000 psig PSV discharging into 50 ft of 2-inch pipe,
[ul]
[li]then critical pressure (right at the PSV) is 5276 psia. [/li]
[li]That is a bunch of pressure, so it will set up a second standing wave in the pipe with a critical pressure at 2780 psig. [/li]
[li]Still too much so there would be a third wave at 1462 psig.[/li]
[li]If I'm not at the end of the pipe yet, then I'll have a 4th wave maintaining 765 psig. [/li]
[li]Next would be 397 psig, [/li]
[li]followed by 203 psig, [/li]
[li]followed by 100 psig, [/li]
[li]then 46 psig, [/li]
[li]then 17 psig which would probably tend to move back up the pipe a bit and actually discharge to atmosphere. [/li]
[/ul]
The velocity at every point in the pipe is equal to every other point because they're all at sonic velocity. The theoretical mass flow rate is significantly lower at each standing wave, and you'll eventually reach a pseudo-steady-state point where the mass flow rate (which actually must be the same at every point in the pipe) will reach an equilibrium somewhere around v(sonic)*A(pipe)*rho(last wave).
We all learned that in incompressible flow, pressure is directly related to the square of velocity. Sonic flow is not incompressible, and as cross sectional area increases the standing wave is spread thinner and velocity is able to increase above Mach 1.0 (velocity increases to the point that the standing wave for that cross section is adequate to hold back the new velocity > 1.0 Mach). That is the only way I know of to get greater than the speed of sound in a pipe. Contentions that the traffic jam of standing waves I described above result in super sonic flow aren't supported in the data.
David Simpson, PE
MuleShoe Engineering
Please see FAQ731-376 for tips on how to make the best use of Eng-Tips Fora.
"It is always a poor idea to ask your Bridge Club for medical advice or a collection of geek engineers for legal advice"