Jerry,
Apologies for the delay but as usual time is short.
The calculation we use in our work includes;
- an orifice plate calculation (Crane) for the nozzle
- adiabatic expansion in the bottle over a short time period
- heat transfer from ambient to bottle to gas over the same time period and an adjustment to temperature
- Joulle Thompson cooling across the orifice
Several key inputs are estimated with the most important for this case being the Overall Heat Transfer Coefficient
- typical values air, free convection 3 - 35 W/m2.°C, forced convection air 30 - 850 W/m2.°C.
- I tried a few values then settled on 100 W/m2.°C
I started with everything at 25°C
The next most important is the time to blowdown, determined by the orifice size. You did not state whether you had a valve at the bottle outlet or how long it took to blow down.
- I tried 2, 3, 5, 10, 15 mm orifices.
- Resulting blowdown times 0.10, 0.20, 1.02, 2.55, 6.12 secs
- Resulting gas in bottle temp 10.7, -0.12, -26.3, -57.7, -71.5 °C.
- Resulting bottle temp 17.9, 19.6, 21.4, 24.2, 24.6 °C
As you can see with a slower blowdown the air temp has dropped only 14.3°C and the bottle 7.1°C which is in line with what you indicated.
I still have serious doubts about the quick blowdown results as I have experienced air blowdown from vessels and never felt anything like the -70°C condition, then again they did not happen in 0.1 sec.
I would love to see this case run as an undergraduate study of first theory and then for real in the lab.
While our method is far from perfect it provides a useful tool for our design and with time we will find ways to improve it.
Regards
Dennis K-B
Dennis Kirk Engineering