A Plumber’s Guide to Starships Part 2: Gases in Pipes
posted by Michel Lamontagne on April 15, 2014
In this section, we will mainly be looking at gases in pipes from the heat transfer point of view. Although there are plenty of other reasons to move gases around in a starship, starting with fuel feed, for example, heat transfer covers the most complicated cases. The information will be applicable to simpler problems.
Gases are trickier than liquids. They compress and change volume. They turn back into liquids under pressure. They leak. But kg for kg, they have higher specific heat than most liquids, so for a mass driven machine like a starship, they are an interesting choice of heat transfer fluid.
The equations for gases and liquids are remarkably similar from the point of view of a starship plumber. We will basically run through the same steps as for “Part 1: Liquids in Pipes”. Pumps will change for compressors, but the rules for these two types of equipment are the same.
The main differences are geometric, stemming from the large difference in density between the two fluids at the same temperature. Liquid Beryllium, for example, has a density of 1,615 kg/m^{3} at 2200 K. Helium, at this temperature, even if it is compressed up to 3,500 kPa (500 psi) has a density of only 0.75 kg/m^{3}; almost 2,000 times less! Therefore the pipes will be much larger. However, gases can flow much faster than liquids, since they are less viscous. So the cycle time will be shorter than for liquid metals, giving us a much lower coolant weight. But the pipework will be larger, possible leading to heavier radiator weight.
A bit of, this, a bit of that, it’s the classic iterative design process.
Another important parameter is the low thermal conductivity of gases compared to liquids. Taking Beryllium again as our reference, its thermal conductivity is 200 W/mK, while the thermal conductivity of helium is only 0.6 W/mK at the temperatures our system will be using. This means that because of the internal resistance of the gas to heat transfer, its low thermal conductivity, the radiator walls will be significantly cooler than the gas. This results in larger radiators, and more weight for our starships.
The high pressure and low density of the coolant will affect the radiation shield design as well. The coolant will probably run through pipes built into the shield, to limit pipe diameter and reduce hoop stress to manageable levels. And, to avoid catastrophic failures, the system should probably be divided into a multitude of independent circuits that can be maintained individually.
Again, we are interested in limiting the power required to make the gas circulate through the circuit, and in keeping the pipes as small as possible, so the pressure will create less strain in the walls.
About Gases
Before we start our design work, we should get a little more familiar with the behavior of gases and the equations that govern this behavior. Our starting point is the Ideal Gas Law:
P*V = n*R*T (1)
Delightfully simple, isn’t it? But a bit abstract, so we will change right away to a more useful version:
ρ=Pw/RT (2)
Where:
Ρ = Gas density (kg/m^{3})
P = Gas Pressure (kPa)
W = Molecular Weight (g/mole)
R = Ideal Gas Constant = 8.314 J/kg mole
T = Absolute Gas Temperature (K)
The interesting thing about the molecular weight expressed in g/mole is that it is exactly the same number as the atomic weight (or relative atomic mass) of the substance, that can be found in the periodic tables. We will round off the fractions, because real gases aren’t perfect anyway, so who needs so much precision?
Gas 
molecular weight 
Hydrogen 
2 
Deuterium 
4 
Helium 
4 
Nitrogen 
28 
Oxygen 
32 
Fluorine 
38 
Neon 
20 
Table 1, gas molecular weight
In this form we can easily find the density of a gas at various temperatures and pressures, and from the density determine the weight of the coolant. The density will also be required later for our pressure drop calculations.
Another property we will want to determine about the coolant gas is its viscosity, since this changes with the temperature of the gas and we need to know it to determine the pressure drop, and therefore the size and power of the compressors. Mr. Sutherland did the work for us here and his findings are compiled in the following equation and table:
μ = μ_{0} x ((T_{0}+C)/(T+C)) x (T/T_{o})^{3/2} (3)
Where:
T = actual gas temperature (K)
T_{0} = reference gas temperature (K)
μ = actual gas viscosity (Ns/m^{2})
μ_{0} = reference gas viscosity (Ns/m^{2})
and C and T_{o} can be found in the following table:
Gas 
Sutherland Constant 
Reference Temperature 
Reference Viscosity 
Viscosity at System Temperature 
Specific Heat Capacity 

C 
T_{0} 
μ_{0} 
μ at 2200 K 
C_{p} 


[K] 
[μNs/m^{2}] 
Ns/m2 
kJ/kgK 
Air 
120 
291.15 
18.27 
0,000067 

Nitrogen 
111 
300.55 
17.81 
0,000063 
1.3 
Oxygen 
127 
292.25 
20.18 
0,000075 

carbon dioxide 
240 
293.15 
14.8 
0,000066 

carbon monoxide 
118 
288.15 
17.2 
0,000064 

Hydrogen 
72 
293.85 
8.76 
0,000029 
14 
Deuterium 
? 
? 
? 
? 
5.1 
Ammonia 
370 
293.15 
9.82 
0,000052 

sulfur dioxide 
416 
293.65 
12.54 
0,000070 

Helium 
79.4 
273 
19 
0.000067 
5.4 
Beryllium 
n/a 
n/a 
n/a 
0,0011 
1,92 
Table 2, gas properties and factors for viscosity calculations, beryllium is included for comparison purposes
<If anyone can find the viscosity properties of deuterium, it would be much appreciated!>
The Parameters
There are many possible arrangements for this system of a hot gas circulating in pipes. The one illustrated here is simply one of the simplest to calculate, but probably not the one that will be eventually chosen.
Here are the cooling process parameters_{(1)}:
The power to be handled and radiated out is 100 GW total, 50 GW per side.
The coolant is helium gas at 3,500 kPa (500 psi).
It is heated, and then cooled, from 1,800K to 2,600K.
The mass flow is 24,000 kg/s (24 tonnes per second), or 32,000 m^{3}/s in volume flow.
The average density of the gas is 0.75 kg/m^{3}.
The Geometry
The system is divided into 100 sections. Each section is made of two elements: a 30m diameter heat exchanger, about 2m thick, which holds the radiation absorbing material; and a long pipe loop, stretching 300m out and returning, for a total length of 600m, that serves as the radiator. There are 50 loops of pipe per side, 1 GW per loop, and 240 kg/s, or 320 m^{3}/s of helium flowing in each loop.
The Calculations
The equation linking compression power to gas flow:
W=(Q*dp)/n (4)
Where:
W = compressor power (kW)
Dp = pressure change (kPa)
Q = flow = 32,000 (m^{3}/s)
n = compressor efficiency (usually between 0.6 and 0.8)
The pressure drop in the pipe is calculated by the Darcy–Weisbach equation:
dp = f*(L/D)*((p*v^{2})/2) (5)
Where:
dp = pressure drop (Pa)
f = friction factor
L = length of pipe (m)
D = Diameter of pipe (m)
ρ = density = 0.75 (kg/m^{3})
v = average gas velocity (m/s)
To choose the velocity, we can apply the upper limit of industrial design velocities for gases of about 100 m/s. We find that the cycle time will be 600/100 = 6 seconds in the pipes, plus about 2 seconds in the radiator. 8 seconds x 24 tonnes per second = 192 tonnes. That’s very good (not too heavy), especially if we compare this to the 1,425 tonnes we calculated for liquid metals.
With the fluid velocity known, we can determine the pipe size required for each of the 100 cooling circuits. First the area of each pipe:
Q = v*A (6)
Q = Flow = 320 m^{3}/s
v = velocity = 100 m/s
A = Area = 3.2 m^{2}
and then the diameter:
D = 2*(A/π)^{1/2} (m) à D = 2*((3.2)/π)^{1/2} = 2m in diameter (7)
All we are missing now to determine the pressure drop is (f), the friction factor. The Haaland friction equation is used to find it:
f = (1.8*log(((e/D)/3.7)^{1.11} + (6.9/Re)))^{2} (8)
Where:
e = pipe roughness factor (m)
R_{e} = Reynolds number (dimensionless)
D = pipe diameter (m)
Roughness factors (e) is about 0.0000015 for the smooth pipe we should be using.
And the Reynolds number (R_{e}):
R_{e} = ρvD/μ (9)
Where:
ρ = density = 0.75 (kg/m^{3})
v = average fluid velocity = 100 (m/s)
D = Diameter of pipe = 2 (m)
μ = dynamic viscosity = 0,000067 (kg/m*s)(Pa*s) , From equation 3, and table 2.
Putting it all together, for helium at about 2200K
9) R_{e}= 0.75 * 100 * 2/0,000067 = 2,238,800
8) f = (1.8*log(((0.0000015/1.1)/3.7)^{1.11} + (6.9/2 238 800)))^{2} = 0.0102
5) dp = 0.0102*(600/2)*((0.75*100^{2})/2) = 11,400 Pa = 11.4 kPa
4) W = (320*11.4)/0.7 = 5,240 kW = 5.2 MW
For 100 circuits, the power is 520 MW.
We can also do the hoop stress right away:
σ = Pr/t (10)
Where:
σ = Hoop stress (Pa)
P = Internal pressure 3,500,000 (Pa) + pressure from the compressor 11,400 (Pa)
r = Interior pipe radius (m)
t = Pipe wall thickness (m)
And the relationship with the material of the pipe is:
S/SF = σ (11)
S = Material tensile strength (Pa)
SF = Safety Factor
σ = the allowable hoop stress (Pa)
The pipes are very large; if we chose the same middle of the range carboncarbon pipe as we had for liquid metals, at 360 MPa, the radiators would weight 14,000 tonnes! So we will move up to stronger, more speculative stuff; carbon nanotube reinforced ceramics, with a tensile strength of 3,600 MPa.
S/SF = σ = 3,600 MPa / 3 = 1,200 MPa or 1,200,000 kPa
For our radiator, 2 m in diameter with a 1 m radius, the wall thickness will be:
t = Pr/σ = 3511*1 /1.200000 = 0.003m, or about 3 mm thick
So each pipe wall will have a volume of 0.003 x 600 x 2 x pi = 11 m^{3}
With a density of 1500 kg/m^{3}, each pipe will weigh 16 tonnes.
For 100 circuits, the total weight will be 1,600 tonnes. Adding the 192 tonnes of coolant we calculated for Equation #5, we find a total weight of about 1,800 tonnes.
The Analysis
The pressure drop for gases in this arrangement is very low, compared to liquid metals. However the volume flow rate is huge. The coolant races through the system in just a few seconds and everything is moving very fast. There is very little energy accumulated in the system, so it has very little buffer. If the flow stops for as little as a minute, the radiation shield will melt, probably with explosive results. So we must plan for the drive to be able to shut down quickly.
The poor heat transfer characteristics of gases tend to negate their excellent heat capacity and low weight. This leaves them pretty much equal to liquid metals. The high pressures required are hard on the pipes, however. The large diameter of the radiator pipes forces us to move to more speculative materials to keep the weight reasonable.
We haven’t worked out the details of the heat exchanger yet. That’s for a future essay. The pressure loss in a heat exchanger is usually the main pressure loss of a piping system, just below the piping itself. So we can guess that the total pressure drop of the system might be twice what we have calculated, for a total power of 540 MW x 2 = 1080 MW. So the power require to circulate the gas may be up to twice the power required to circulate an equivalent liquid metal, at least for the present model.
It’s pretty obvious that we haven’t optimized the circuit. We should look into reducing the pipe diameter, either by increasing the gas velocity over standard values (as we did for liquid metals) or by increasing the radiator size and increasing the number of circuits. If we can determine materials that will resist corrosion from hot gases, we might switch to hydrogen cooling for its much higher specific heat. Deuterium does not offer this advantage, since its specific heat is lower than the one of helium, so its use as a substitute for helium is problematic.
Alternatively, a compound gas might be interesting; a mixture of nitrogen and helium, or of neon and helium, might provide better overall performances.
This guide provides the equations. The work is up to the starship designer.
Are we done yet?_{(2)}
Not quite. We’ve mentioned several times that the surface temperature of a radiator will be lower than the cooling fluid’s temperature, but we have not given the details yet. Here they are.
Surface Temperature of Radiators
In our example for radiators with liquid metal coolants, we assumed the temperature of the pipe walls, and so of the radiator, would be the same as the temperature of the liquid metal. In reality, this is not the case, but the error is not very large since liquid metals are very conductive and the pipes are small. For gases, we really need to take the thermal conductivity of the gas into account, the convection rate, as well as the size of the pipes and the other physical characteristics of the gas. Generations of experimenters have labored tirelessly to produce the following equation, linking all these elements together, the Gnielinski correlation_{(3)}:
h = ((f/8) x (R_{e} – 1000) x Pr) / (1+(12.7 x (f/8)^{1/2 }x (P_{r} ^{2/3} 1))) x k/D (12)
Where:
h = convective heat transfer rate for a gas in a pipe (W/m2K)
f = friction factor
R_{e} = Reynolds number
P_{r} = Prandtl number
D = Pipe diameter (m)
k = Thermal conductivity of the gas
Admittedly, it’s a lot less elegant than the Ideal Gas Law.
The Reynolds number, the friction factor (f) and the pipe diameter we already know.
P_{r} is the Prandtl number, a dimensionless number, standing in for the following:
P_{r} = c_{p}μ/k (13)
Where:
c_{p }=specific heat = 5,400 (J/kgK)
μ = dynamic viscosity = 0.00067 (Ns/m^{2}) from Equation #3
k = thermal conductivity (W/mK) from Equation #4
We need to determine the thermal conductivity (k) for both Equations #12 and #13. It varies with the temperature and the pressure. Again we turn to experimental results and to the equations that were worked out during the great age of nuclear experimentation, in the fifties and sixties_{(4)}:
14a) Helium: k=2.684×103*T0.71
14b) Deuterium: k= 0.7073+0.02368x103T+0.1048x106T2 (mW/cmK)
14c) Hydrogen: k=0.7897+0.03623x103T+0.01809x106T2 (mW/cmK)
14d) Nitrogen: k=0.051+0.7438x103T0.1573x106T2 (mW/cmK)
14e) Oxygen: k=0.07979+0.6671x103T0.05479x106T2 (mW/cmK)
Putting the numbers into equations 13) and 14) we find:
14) k = 2.684×103*22000.71= 0.6
13) P_{r} = 5,400*0.000067/0.6 = 0.57
12) h= ((0.0102/8)x(2,238,8001000)x0.57) / (1+(12.7 x (0.0102/8)1/2 x(0.57 2/3 1)))x0.6/2 = 600
As our final step, we want to determine the wall surface temperature. To solve this we will use the fact that the surface temperature for the radiation equation is the same number as the wall temperature for the convection in the pipes, and that the power that radiates out from the radiator has the same value as the power that moves out from the gas and to the radiator walls.
Q_{r} = eBA(T_{wall}^{4}T_{space}^{4}). (15)
Q_{t} = hA(T_{fluid}T_{wall}) (16)
We can combine these two equations into a single one, since Q_{r} = Q_{t} and find our last equation:
eB(T_{wall}^{4}T_{space}^{4}) – h(T_{fluid}T_{wall}) = 0 (17)
Where
e = emissivity
B = StephenBoltzman constant
T_{space} = 10K
T_{fluid} = 2200K
h = thermal conductivity = 600 (W/m^{2}K)
We know all the numbers here, except T_{wall}. To find this temperature, you can either embark on a long and complex search for analytical solutions, or use the brute force approach available in most spreadsheet programs, using the solve function in Excel, for example.
In this case, the solution is 570 K of difference. So for an average gas temperature of 2200 K, the average radiator surface temperature will be 1630 K.
It is left as an exercise to the reader to determine the increase in radiator area this represents.
Tools of the Trade
You will find the tables and equations mentioned in this article in the following spreadsheet:
Plumbers guide to starshipsGases
Notes
(1) These are the parameters used in the Starship radiator post, and are repeated in the Tools of the Trade spreadsheet.
(2) Kitty was found here on the Internet: http://bit.ly/1kNw1aJ
(3) The Gnielinski correlation is actually expressed as a function of the adimensional Nusselt number, not of convective heat transfer. But since its only use for us here is to combine the equation with the heat transfer equation h=Nu D / k, I did it right away and saved us all a step. The Dittus Boelter correlation is also often used for heat transfer in pipes, but it goes way off in the ranges of temperature we are working in, so we will not use it here. And since we have the Haaland equation to determine ‘f’ directly, there is no longer a problem with using the Gnielinski correlation for simple calculations.
((4) The actual equation for helium is more complex and takes pressure into account
k = (2.682×103)(1 + 1.23×103*P)*T0.71(10.0002*P) it can be found here: http://1.usa.gov/1lBsDgu.
For the other gases, the equations were taken from ‘Transport properties of gases and gaseous mixtures at high temperatures, S.C.Saxena, 1970’. The curve fitting was only tested by the authors up to 1100 K, so the results of these equations should be taken with a grain of salt.
Excellent work! Michel
Thanks! This is a work in progress; if you think of anythink that might improve it, i’ll be glad to make the changes.
Pingback: Plumbers guide to starships Part 2  Mostly drawings
very good and interesting series of posts.
an interesting possibility is to analyse the tradeoff between splitting the gas pipe between multiple ones, increasing the net gaspipe contact area, and partially compensate the relatively low thermal conductivity of the gases. Although I’ve seen this applied only for small hydraulics, so I’m not sure how that would scale to a configuration with pipes running hundreds of meters. You made a mention about multiple circuits on the post, so it is clear that you already knew this though
Hi Charles,
I did the simulations: for a 2m pipe the temperture difference is 570K, as mentionned in the paper. for a 10cm pipe, the temperature difference is 435K and for a small 1cm pipe the temperature difference is 325K. This corresponds to increases in efficiency of 37% and 75%. So it is a very effective solution, from the thermal point of view. The physical arrangement is probably heavier though. The numbers are here: http://bit.ly/RoG7mz.
Regards,
Michel Lamontagne
yes Indeed a great set of articles,
my math skills are poor however i still struggle on with my very small project,
http://yellowdragonblog.com/2014/03/LH2%20volume%20~1,700%20ft3%20LO2%20volume%20~600%20ft3%20The%20Centaur%20tank%20sidewall%20is%20~25%20ft%20long%20and%2010%20ft%20in%20diameter.%20You%20can%20calculate%20the%20area%20from%20this.%20We%20fly%20SOFI%20on%20the%20sidewall%20on%20all%20missions%20for%20insulation.%20Some%20missions%20also%20have%20MLI.26/hybridsmallfissionreactorpressurizedargonheliumcapsuleasionpoweredsepstage/
A ULA engineer who is interested in my idea provided this data for the Centaur in space stage;
begin quote
“LH2 volume ~1,700 ft3
LO2 volume ~600 ft3
The Centaur tank sidewall is ~25 ft long and 10 ft in diameter. You can calculate the area from this.
We fly SOFI on the sidewall on all missions for insulation. Some missions also have MLI.”
end quote
The Engineer did not mention the alloy that makes up the Centaur,s tanks
So I assume a small fission reactor radiator in the spent LO2 with gasses Xenon at 600 ft3 and the Effective radiator temperature 475 K
Radiator size required 5.035 m2
Radiator length required 240 cm
This data is from the decadel survey small fission reactor study
I assume gas law and the surface area losses through the fuel tank walls and friction in the Stirling engines
because the SFR radiator has displaced LO2 in the tank needed for the chemical burn we need to increase the size of the O2 tank thus altering the tank volume needed to compute the gas law and surface area equation!
I was inspired by the ISS radiator NASA paper as it is the largest in space radiator to date!
if you would like to see this paper I have posted it on my blog with my take on how to recover waste heat from the ISS radiator.
would there be some starship designs where radiators could be re purposed?
http://yellowdragonblog.com/2014/04/24/theisseeatcsradiatororucentaurcogenerationexperiment/
Pingback: Project Persphone: British Scientists Building 'Living Space Ark' To Save Humanity  Really Upbeat
Pingback: Project Persphone: British Scientists Building 'Living Space Ark' To Save Humanity  Worth a Giggle
Pingback: Project Persphone: British Scientists Building ‘Living Space Ark’ To Save Humanity 
Would supercritical CO2 be a middle ground of sorts? Ease of gas pumping, but thermal properties more like a liquid?
It might be interesting. Its molecular weight is 44, or 10 times denser than helium. However its specific heat is only 0,84 compared to 5 for helium, so my guess is that you would need a lot more weight, but that the overall volume would be smaller, so lighter pipes; so not a clear advantage.
Oups! I have updated table 1, there was a mistake in the molecular weight table of the non inert gases. I will also be providing a more exact equation for compressor power in my upcoming pumps and compressor post.
Regarding deuterium viscosity:
Assael M.J.; Mixafendi S.; Wakeham W.A. The Viscosity of Normal Deuterium in the Limit of Zero Density J. Phys. Chem. Ref. Data 1987 Vol. 2, No. 2
A pdf reprint is available on the NIST website http://www.nist.gov/data/PDFfiles/jpcrd318.pdf
Thanks! That fills a big hole in my database.
Michel, I admire your work on this. Have you considered fractal metamaterials for heat transfer between the shield and the radiator? Link: http://nextbigfuture.com/2013/05/metamaterialsremoveheatandcontrol.html#links
This is a near field effect so it won’t replace the need for radiators. I realize that the fractals are a bit speculative and one can’t easily put numbers to it, but what are your thoughts on this?
Hi Devon,
Thanks for the link.
Frankly, it looks like a series of cavities transmitting infra red radiation, added to closed convection cells. It might be a good replacement for pipe walls, as long as the overall product was gas tight. I’m a bit skeptical about ‘special’ electromagnetic waves, whatever they are.
Yes, the word “special” is not too far from “magic”, and a poor descriptor of evanescent wave nearfield effects. But it looks like the rate of heat transfer exceeds StefanBoltzmann by orders of magnitude, certainly much higher than convection and conduction also.
Correct me if I’m wrong, but I don’t think there’s necessarily a gas medium? Seems like a fine vacuum effect unless I’m missing something.
The only bottleneck may be conduction of the heat out of the neutron shield into the metamaterial, and then out of the metamaterial into the radiator wings, but this should be no worse than traditional conductive transfer. The advantage would be mass savings, no moving parts, no power or pumps, etc with regard to the plumbing.
Hi Devon,
Do you have a link to some more substantial material than the Bigfuture page? It’s a new field for me and I need to study a bit more to understand and give useful comments.
Okay, so the company making the claims doesn’t disclose very much about how it works, other than basic description like this (only slightly better than that Bigfuture page):
http://www.businesswire.com/news/home/20110120006450/en/FractalsOfferCool#.U_7S8oJXuY
For example, they claim that the metamaterial structure itself is responsible for the passive and directed flow of IR evanescent waves.
Although I can’t find much else on fractal metamaterials being used for radiative heat transfer, there’s quite a bit out there on the general theory and practice of nearfield radiative heat transfer.
Basics about evanescent waves:
http://en.wikipedia.org/wiki/Evanescent_wave
Here is a more serious intro to radiative heat transfer. This one’s pretty good. If you have a big library go request this book, or just scan the pages here. Scroll down to CH 24 for the nanoscale bit. Even if this is over your head, he has a good list of reference papers at the end of the chapter. (Intermediate E&M course background would help you here) http://books.google.com/books?id=J2KZq0e4lCIC&pg=PA805&lpg=PA805&dq=evanescent+radiative+heat+transfer&source=bl&ots=_z84t7SWb1&sig=LO9BUSBJuTkTKuJh0Y7RGCQ6RxE&hl=en&sa=X&ei=8crU97_MufTiwLFzYHQBQ&ved=0CD4Q6AEwAw#v=onepage&q=evanescent%20radiative%20heat%20transfer&f=false
Follow this link, scroll to the bottom for a list of references to theory papers, some of which seem like good intros from the abstracts. I don’t have subscriptions to the journals, maybe you’ll have more luck? http://www.ichmt.org/rad13/images/abstracts/rad13%20poster%2019.pdf
Here’s a paper free for download, but it might be too indepth if you’re new to the field. The gist of it is that the radiative heat transfer between materials is strongly dependent upon material properties (no kidding): http://download.springer.com/static/pdf/805/art%253A10.1007%252Fs1143401145869.pdf?auth66=1409381024_fb947863b9a8e6a11050f224357d1141&ext=.pdf
Hopefully at least some of those links are useful…
Interesting links, thanks! So it seems to be a kind of heat transfer using infra red waves transmitted in the cavities of the nano materials? And for nano structures that are smaller, or in the range of the wavelength, the transfer rate of the infrared waves is higher than the one predicted by the radiation equation, with little attenuation? It would be a bit like accelerating the phonons of classic thermal conductivity, creating, more or less, (there seem to be some losses) thermal superconductors? In the end we still have to radiate through classic radiation, but I think I see your point that the heat transfer mechanism from the source to the radiators would be much simpler this way. Very exciting! I will try to find out more and perhaps work it into a future article. In a sense, it would make the whole concept of space plumbing moot 😉 and replace it with something better. So it goes with many professions!