Physics:Nuclear reactor heat removal

The removal of heat from nuclear reactors is an essential step in the generation of energy from nuclear reactions. In nuclear engineering there are a number of empirical or semi-empirical relations used for quantifying the process of removing heat from a nuclear reactor core so that the reactor operates in the projected temperature interval that depends on the materials used in the construction of the reactor. The effectiveness of removal of heat from the reactor core depends on many factors, including the cooling agents used and the type of reactor. Common liquid coolants for nuclear reactors include: deionized water (with boric acid as a chemical shim during early burnup), heavy water, the lighter alkaline metals (such as sodium and lithium), lead or lead-based eutectic alloys like lead-bismuth, and NaK, a eutectic alloy of sodium and potassium. Gas cooled reactors operate with coolants like carbon dioxide, helium or nitrogen but some very low powered research reactors have even been air-cooled with Chicago Pile 1 relying on natural convection of the surrounding air to remove the negligible thermal power output. There is ongoing research into using supercritical fluids as reactor coolants but thus far neither the supercritical water reactor nor a reactor cooled with supercritical Carbon Dioxide nor any other kind of supercritical-fluid-cooled reactor has ever been built.

Theoretical framework

The thermal energy produced in nuclear fuel comes mainly from the kinetic energy of fission fragments. Therefore, the heat generated per volume unit is proportional to the fraction of nuclear fissionable fuel burned in the unit of time:

[math]\displaystyle{ Q=a\phi\Nu_0\sigma_f }[/math]

where [math]\displaystyle{ N_0 }[/math] represents the number of atoms in a cubic meter of fuel, a is the amount of energy released in the fuel in each fission reaction (~181 MeV),  [math]\displaystyle{ \phi }[/math] is the neutronic flux, and [math]\displaystyle{ \sigma_f }[/math] is the effective section of the fission.

The total heat produced in the nuclear reactor is:

[math]\displaystyle{ Q=a\overline{\phi}\Nu_0\sigma_fV }[/math]

where [math]\displaystyle{ \overline{\phi} }[/math] is the mean neutronic flux and V is the fuel volume (normally measured in [math]\displaystyle{ m^3 }[/math]).

Recovery of this amount of heat is achieved by using cooling fluids whose temperature at the entrance to the reactor channel [math]\displaystyle{ T_{f_0} }[/math]  will increase with the distance traveled in the channel. The thermal balance of the channel is expressed by the relationship:

[math]\displaystyle{ Gc_pdT_f=q(z)n\pi r_0^2dz }[/math]

where [math]\displaystyle{ G }[/math] is the flow rate of the cooling agent, [math]\displaystyle{ c_p }[/math] is the specific heat at constant pressure, [math]\displaystyle{ dT_f }[/math] is the increase in the temperature of the fluid after passing a distance [math]\displaystyle{ dz }[/math] in the channel, [math]\displaystyle{ q(z) }[/math] is the heat generated per unit volume of the fuel, [math]\displaystyle{ r_0 }[/math] is the fuel cell radius and [math]\displaystyle{ n }[/math] is the number of channel bars.

Under these conditions, the temperature of the cooling agent at distance z travelled into the cooling channel inside nuclear reactor is obtained by integrating the previous equation:

[math]\displaystyle{ T_f(z)=T_{f_0}+n(\pi r_0^2)(Gc_p)^{-1} \textstyle \int\limits_{0}^{z} \displaystyle q(z)dz }[/math]

The difference between the temperature of the outer surface of the tube-channel [math]\displaystyle{ T }[/math] and the temperature of the fluid is obtained from the relationship:

[math]\displaystyle{ \varphi=h(T_s-T_f) }[/math]

where  [math]\displaystyle{ \varphi }[/math] is the local heat flow on the casing - cooler contact surface unit and [math]\displaystyle{ h }[/math] is the heat transfer agent casing-cooling agent.

The heat discharge from the PWR and PHWR reactors is made by pressurized water under forced convection. The general expression for determining the transfer coefficient is given by the Dittus - Boelter equation:

[math]\displaystyle{ \aleph u=aRe^{0.8}Pr^{0.4} }[/math]

where [math]\displaystyle{ \aleph u }[/math] is Nusselt's number ( [math]\displaystyle{ \aleph u=\tfrac{hd_n}{\lambda} }[/math], [math]\displaystyle{ h }[/math] is the heat transfer coefficient, [math]\displaystyle{ d_n }[/math] is the equivalent diameter, [math]\displaystyle{ \lambda }[/math] is the thermal conductivity of the fluid); [math]\displaystyle{ a }[/math] is a constant ([math]\displaystyle{ a }[/math]=0.023);  [math]\displaystyle{ Re }[/math] is the number of Reynolds ( [math]\displaystyle{ Re=V\tfrac{d_n}{\mu/\rho} }[/math]) [math]\displaystyle{ V }[/math] is the average velocity of the fluid in the section considered, [math]\displaystyle{ \rho }[/math] is the density of the fluid and [math]\displaystyle{ \mu }[/math] is its dynamic viscosity); [math]\displaystyle{ Pr }[/math] is the number of Prandtl ([math]\displaystyle{ Pr=\tfrac{\mu c_p}{\lambda} }[/math]).

If the flow of the fluid is made under conditions of a great difference between its temperature and the contact surface, the transfer coefficient is determined from the relationship:

[math]\displaystyle{ \tfrac{hd_n}{\lambda}=0.023 Re^{0.8}Pr^{1/3}(\mu/\mu_s)^{0.14} }[/math]

where [math]\displaystyle{ \mu_s }[/math] is the dynamic viscosity of the coolant at the temperature of the adhering fluid film at the surface of the casing. The relation presented above is valid in the case of a long channel with [math]\displaystyle{ L/d_n\gt 60 }[/math], where [math]\displaystyle{ L }[/math] is the length of the channel.

The transfer coefficient for cooling the pipes by natural convection is obtained from:

[math]\displaystyle{ \aleph u=0.5(Gr\cdot Pr)^{0.25}(\tfrac{\mu}{\mu_s})^{0.25} }[/math]

where [math]\displaystyle{ Gr }[/math] is the Grashof number given by the expression:

[math]\displaystyle{ Gr=\beta\Delta T\tfrac{gd_n^3\rho^2}{\mu^2} }[/math]

We use the notation  [math]\displaystyle{ \beta }[/math] for the volume expansion coefficient of the fluid, [math]\displaystyle{ g }[/math] is the gravitational acceleration and [math]\displaystyle{ \Delta T }[/math] is the difference between the average wall temperatures of the casing and the cooling agent.

In boiling water cooled reactors (BWR) and partly in pressure water cooled reactors (PWR and PHWR) the heat transfer is made with a vapor phase in the cooling medium, which is why this type of heat transfer is called heat transfer in a biphasic system. This allows obtaining much higher transfer coefficients than the one-phase heat transfer described in the Dittus-Boelter equation.

Coolant flow

Increasing the flow of heat, reducing the agent flow and lowering the pressure can lead to increased temperature of the cooled surface. If the temperature of the fluid in the channel section that we consider is lower than the boiling temperature under local pressure conditions, the vaporization is limited to the immediate vicinity of the surface and in this case the boiling is called submerged boiling. There is no proportionality between the heat flow and the difference between the surface temperature and the coolant temperature that allows the definition of a heat transfer coefficient similar to the one-phase case. In this situation we can use the equation of Jens and Lottes, which establishes a connection between the difference [math]\displaystyle{ \Delta T_s }[/math] between the surface temperature and the boiling temperature of the cooling agent under local pressure conditions [math]\displaystyle{ P }[/math] below the thermal flux [math]\displaystyle{ \varphi }[/math]:      

[math]\displaystyle{ \Delta T_s=a\varphi^{1/4}exp(-P/P_0) }[/math]  

where [math]\displaystyle{ a=7.9^\circ C cm^{0.5}W^{-0.25} }[/math] and [math]\displaystyle{ P_0=6.43 MPa }[/math]

If the temperature of the fluid in the channel section considered is slightly higher than the boiling temperature under local pressure conditions, the heat transfer is by boiling with nucleation, forming vapor bubbles trained by the cooling agent (that becomes biphasic throughout its entire volume). However, the vapor content is relatively small and the continuous phase remains the liquid phase. The vapor content of the PHW-CANDU reactor is about 0.03-0.04 kg steam / kg of agent, thus increasing the amount of heat transported by the unit mass of agent by over 10%. If the cooled surface temperature far exceeds the boiling temperature of the cooling agent in the channel section, the vapor content of the agent increases considerably, the continuous phase becoming the vapor phase and the liquid phase becoming only a suspension between vapors. The cooled surface remains covered with a liquid film which still provides a very high heat transfer coefficient, [math]\displaystyle{ \thicksim6\cdot10^4 W/m^2\cdot K }[/math] at BWR compared to [math]\displaystyle{ \thicksim3\cdot10^4 W/m^2\cdot K }[/math] at PWR. The film of liquid is continuously fed with drops from the agent suspension.

A further increase in surface temperature leads to a temporary interruption of continuity of the liquid film adhering to the cooled surface. Watering of the surface continues, however, by the drops of liquid in the suspension that are present in the cooling agent as long as the heat flow remains below a value that depends on local conditions (value that is called critical flux). Over this flux there is a thermal transfer crisis characterized by a sudden decrease in the transfer coefficient due to the presence of only one-phase transfer. The heat transfer coefficient in the pre-crisis period can be determined from the relationship:

[math]\displaystyle{ \tfrac{\aleph u}{Re^{0.8}Pr^{0.4}}=k_1B_0+k_2(1/\chi)^n }[/math]

where [math]\displaystyle{ \chi=((\Delta p)_l/\Delta p_g)_g)^{1/2} }[/math]

In these formulas the following notations were made: [math]\displaystyle{ dp }[/math] is pressure losses for the two phases (water and vapors), [math]\displaystyle{ B_0=\varphi/d_n H_{gl} }[/math] ( [math]\displaystyle{ \varphi }[/math] - the thermal flux, [math]\displaystyle{ H_{gl} }[/math]- the enthalpy of the biphasic liquid-gaseous mixture). The heat transfer coefficient during the crisis is related to the critical heat flow [math]\displaystyle{ \varphi_{cr}(W/m^2) }[/math] through a linear relationship, of the equation type [math]\displaystyle{ \varphi=h(T_s-T_f) }[/math] that was presented before:

[math]\displaystyle{ \varphi_{cr}=h_{cr}(T_s-T_{st}) }[/math]

Where [math]\displaystyle{ T_s }[/math] is the temperature of the surface in thermal transfer crisis, and [math]\displaystyle{ T_{st} }[/math] is the temperature of the vapor at saturation.

The critical flow is obtained by using the Kutateladze's formula:

[math]\displaystyle{ \rho_{cr}=0.14r\rho_v^{1/2}[g\sigma(\rho_l-\rho_v)]^{1/4} }[/math]

where [math]\displaystyle{ r }[/math] (J/kg)is the latent heat of vaporization, [math]\displaystyle{ \rho_l }[/math] and [math]\displaystyle{ \rho_v }[/math] are density of the liquid and saturation vapor, [math]\displaystyle{ \sigma }[/math] is the superficial tension in N / m and [math]\displaystyle{ g }[/math] is the gravitational acceleration. The heat transfer to the gas-cooled reactors is carried out by forced convection. For a gaseous thermal agent, the heat transfer coefficient can be deduced from a relation of the type Dittus-Boelter, but taking into account, for the intervening sizes, the values corresponding to the average temperature of the fluid film denoted by the index [math]\displaystyle{ m }[/math]:

[math]\displaystyle{ \aleph u_m=0.0205Re_m^{0.8}Pr_m^{0.4} }[/math]

which differs in the use of water by a slightly lower value of the coefficient a.

Forced flow relationships established for fluids are also not valid for liquid metals. The coefficient of heat transfer for circular pipelines with constant heat flux, where the heat evacuation is achieved by the turbulent flow of the molten metals, can be estimated with a relation of the type:

[math]\displaystyle{ \aleph u=0.625Pe^{0.4} }[/math]

where [math]\displaystyle{ Pe }[/math] is the number of Peclet ([math]\displaystyle{ Pe=Vd_n\rho c_p/\lambda }[/math]).

Examples of heat evacuation hydrodynamic parameters

For exemplification of the above formulas the hydrodynamic parameters of some types of reactors can be found in the following table:

Reactor	Cooling agent	Pressure [Mpa]	Speed of the fluid [m/s]	Re	Pr
G1	[math]\ce{ Acr }[/math]	0.1	20	26,000	0.73
EL-4	[math]\ce{ CO2 }[/math]	5.9	50	820,000	0.7
MTR	[math]\ce{ H2O }[/math]	0.1	1	9,000	4.57
VVER-440	[math]\ce{ H2O }[/math]	12.5	3.7	330,000	0.90
FBR	[math]\ce{ Na }[/math]	-	5	93,000	0.36

G1 and EL-4 are reactors that were built in France, while VVER-440 is a reactor that has been constructed in the Soviet Union.