Heat transfer

Heat transfer is a very complicated process involving many parameters and conditions. This page plays the role of introduction to heat transfer and a documentation on the methods used to calculate the heat transfers applied to our vehicle.

Heat transfers

A heat transfer is basically an energy transfer that can produce work or a change of temperature difference between the interacting parts. It can take three forms: conductive, convective and radiative. In the case of a heat transfer due to the atmosphere around an aircraft, the three mechanisms are effective. This page currently focuses on convective heat transfer, which implies conductive heat transfer. Radiative transfer is probably negligible before the two others in our particular case study.

Introduction to heat convection and conduction

Heat convection occurs when there is a fluid flowing around a control volume at a temperature different than the control volume's. The flow can be either natural, when the fluid is flowing due to density gradients (buoyancy force), or forced, when the fluid is flowing because an external process force it to flow or make the control volume move through the fluid.

When air is the considered fluid, convection is always associated to conduction heat transfer. The reason is that a film is formed on the surface of the object, called the boundary layer, and it is partly steady and at a temperature closer than the temperature of the object than the temperature of the fluid. In that case, conduction applies.

An insulation layer's effect can be easily calculated as its thermal transmittance (U-value) or resistance (R-value), the amount of heat that it allows to be transferred through it. This is directly obtained from the material's thermal conductivity k (unit: W/m.K) and the thickness of the insulation L. R = L/k and U = k/L. Unit of U is W/m^2.K. The transmitter heat is then Φ = A × U × (T1 - T2), in Watt (= Joule/s), where A is the external area of the insulation layer, T1 and T2 are the internal and external temperature. Examples for the transmittance of insulation layers can be found here [1],

Heat transfer coefficient h

The rate of heat loss of a body by convection is proportional to the difference in temperatures between the body and its surroundings, as stated by Newton's law of cooling: , where h is the heat transfer coefficient. This h depends on many parameters (flow rate, surface roughness, fluid properties, and others) and is very hard to calculate accurately. Approximations exist for some conditions and determining them is still an active research topic for some conditions. It will be our main problem in the case of air to aircraft heat transfer.

Approximations and conditions

h should be calculated from Computational fluid dynamics analysis, which requires a good expertise in the domain and complex software to be used. Since we don't have this capability, we explore existing approximations of the value. These approximations are only valid in some very specific conditions, each condition having a different approximation. Assumptions are also made, for example the temperature of the fluid and the body over the area of the heat transfer are assumed constant, as well as the flow rate. Approximation are often based on a fluid temperature equal to the arithmetic mean between the wall and the free stream. This is called the mean film temperature.

The Nusselt number Nu

The Nusselt number Nu is introduced when equating Newton's equation to the conduction heat:

= Nu

The Nusselt number is then the ratio of the temperature gradient at the surface to the reference temperature gradient, meaning that its value indicates the shape of the temperature gradient.

The Prandtl number Pr

The Prandtl number depends only on the fluid and its state, not on a characteristic length.

The Prandtl number controls the relative thickness of the momentum and thermal boundary layers. When Pr is small, it means that the heat diffuses very quickly compared to the velocity (momentum). This means that for liquid metals the thickness of the thermal boundary layer is much bigger than the velocity boundary layer. See page 223 of [2] for other interpretations of the values.

The Reynolds number Re

The Reynolds number gives a measure of the ratio of inertial forces to viscous forces.

It is also useful because its value can indicate if the flow in the boundary layer is laminar (value < 350,000), turbulent (value > 500,000) or transitional between the two, in which case it depends on other factors such as surface roughness.

The Rayleigh number Ra

The Rayleigh number is an indicator for natural convection heat transfer. When it is below the critical value for a fluid, heat transfer is primarily in the form of conduction; when it exceeds the critical value, heat transfer is primarily in the form of convection.

α is the thermal diffusivity  ; β is the thermal expansion coefficient, for an isobaric process it can be approximated to β = 1/T. x is the distance from the leading edge.

Cases of application

We currently use heat transfer to estimate the rate of vaporization of the cryogenic propellants on the ground and during ascent to rocket ignition altitude. The latter could however be easily used to calculate the drag force of the aircraft fuselage on its tank part. In both case we can assume the fluids are incompressible and steady and that the temperatures are constant over the body and in the free stream fluid for a given altitude.

Natural convection for horizontal cryogenic tank

On the ground, the cryogenic propellant tank undergoes a large temperature difference and since it doesn't move it's the natural convection that is at work, if the wind is neglected. An approximation exists for the Nusselt number in this condition: horizontal cylinder in external flow under natural convection.

D is the diameter of the approximated cylinder, in our case we take 0.266m, given by a first evaluation of the rocket mass program.

Air density is derived from the atmospheric model at sea level, corrected to the actual ground temperature. For a 25°C or 298.15 K ground temperature:

• the density (rho) is 1.18391 kg/m^3

and at the mean film temperature of -78.07°C or 195.075 K, given the propellant and tank temperature of 92 K:

• the viscosity µ is 1.30015e-05 Pa.s (calculated using Sutherland's law)
• the kinematic viscosity nu is µ / rho = 1.0982e-5 m^2/s
• the specific heat Cp is 1007.68 J/kg.K
• the Prusselt number is 0.874462
• the thermal conductivity k is 0.0177375 W/m.K
• α is thus k/(rho.Cp) = 1.4868e-5 m^2/s
• β is 1/195.075 = 5.126e-3
• Ra_D is ((9.80665 * 5.126e-3) / (1.0982e-5 * 1.4868e-5)) * (92 - 298.15) * 0.266^3 = 6.3467e10 * 0.266^3 = 1.1945e9.
• finally, h is 0.0177375/0.266 * (0.6 + (0.387 * 1.1945e9^1/6) / (1 + (0.559/0.874462)^9/16)^8/27)^2 = 0.024382 W/m^2.K

Forced convection on aircraft fuselage during flight

Section 5.8 in [2], Heat transfer in turbulent flow over a flat plate, is the closest known answer to this problem. The tank's cylinder can be considered as a flat plate since there is no pressure change in the local y direction. The Colburn analogy is the approximation that applies in this case. It states that the local Nusselt number, assuming constant wall temperature, is Nu_x = 0.0292 Re_x^4/5 Pr^1/3. This equation is however only valid for Re_c < Re_x < 10⁷ and 0.7 < Pr < 100.

From what we calculated with the ISA model, Pr drops from 0.87 at sea level to 0.013 at 30km altitude, and the Reynolds number drops from 4.36714e+07 to 668455. If Re is nearly in the allowed range, Pr is far from it. We can't apply the Colburn analogy for altitudes above 2.5km, where Pr gets lower than 0.7.

References

1. Principles of Heat Transfer: SI Edition, Frank Kreith, Raj M. Manglik, Mark S. Bohn, 2010. Google books

2. Textbook Of Heat Transfer (4th Edition), S.P. Sukhatme, 2006. Google books