Reactor Cooling Rate

The formula Q = ṁ c_p ΔT uses mass flow rate because heat capacity is mass-specific (J/kg·K). Using volumetric flow rate would require multiplying by density, which varies significantly with temperature and pressure in the primary loop. For example, water density at 300°C and 15 MPa is about 750 kg/m³, whereas at 280°C it's about 780 kg/m³—a 4% difference. This variation would introduce errors if a constant density were assumed. In nuclear plant calculations, mass flow is the conserved quantity, so using ṁ ensures accuracy regardless of temperature-dependent density changes.

For a typical 3000 MWth PWR, the primary coolant mass flow rate is about 15,000-18,000 kg/s (about 15-18 tonnes per second). The coolant enters the core at about 280°C and leaves at about 320-330°C, giving a ΔT of 40-50 K. Using water's specific heat at these conditions (~5.5 kJ/kg·K), Q = 16,000 × 5,500 × 45 ≈ 3,960,000,000 W = 3.96 GWth, which matches the reactor power. The product ṁ × ΔT is essentially fixed by the thermal power; if the flow rate increases, the ΔT decreases proportionally, and vice versa.

The formula Q = ṁ c_p ΔT shows that if Q remains constant (reactor power unchanged) and ṁ decreases, then ΔT must increase to remove the same heat. For example, if flow drops by 20%, ΔT must increase by 25% to maintain the same Q. This raises the coolant outlet temperature, which could approach the saturation temperature and increase the risk of boiling or departure from nucleate boiling (DNB). This is precisely why loss-of-flow accidents are a major design basis event—the formula quantifies the thermal response and sets the required trip setpoints.

In a PWR operating at ~15.5 MPa, the specific heat of water increases from about 5.2 kJ/kg·K at 280°C to about 5.8 kJ/kg·K at 330°C. A 10% variation over the core's temperature range means that using a constant c_p would introduce a few percent error in the heat removal calculation. For accurate thermal-hydraulic analysis, the integral form Q = ṁ ∫_{T_in}^{T_out} c_p(T) dT is used, or the enthalpy difference is computed directly from steam tables. Many safety codes use property tables to capture this non-linearity, especially at high temperatures near the pseudo-critical point.

The cooling rate formula gives the heat removed, but the flow rate ṁ required to achieve that heat removal also determines the pressure drop through the core (ΔP = f × (L/D) × (ρv²/2)). The pump must overcome this pressure drop to circulate the coolant. Higher flow rates increase pressure drop, requiring more pump power. The design optimization balances the need for a high enough flow to achieve adequate cooling (as per Q = ṁ c_p ΔT) with the pump's hydraulic capability and power consumption. The pump sizing is directly linked to the required ṁ from the cooling rate calculation.

For a multi-channel core, the total heat removal is the sum over all channels: Q_total = Σ ṁ_i c_p ΔT_i, where i indexes each fuel assembly or subchannel. The flow distribution ṁ_i and ΔT_i vary because of power peaking and flow resistance differences. Hot channels with higher power require higher ΔT_i for the same flow, or if flow is reduced, they become limiting. Thermal-hydraulic codes use this summation to evaluate the margin to DNB and to ensure that the hottest channel does not exceed the critical heat flux. The sum of all channel heat removals must equal the reactor thermal power.

After a reactor trip, the fission power drops to decay heat levels (initially about 6-7% of full power, then decreasing). The formula Q_decay = ṁ c_p ΔT is used to determine the required flow to remove this decay heat. At early times, the decay heat is still high (e.g., 100 MW for a 3000 MWth plant), so the flow must be sufficient to keep the ΔT within limits. If natural circulation is used (no pumps), the flow is driven by buoyancy, and the formula is used to check that the resultant ΔT is within safe limits. This is critical for accident analysis and for sizing emergency cooling systems.

For a PWR, the primary coolant remains subcooled (no boiling) under normal operation, so a single-phase formula Q = ṁ c_p ΔT is valid. For a BWR, boiling occurs in the core, so the coolant is two-phase; the formula becomes Q = ṁ (h_out - h_in), where h is the enthalpy, including the latent heat of vaporization. The two-phase flow requires a more complex model because the specific heat is not constant and the quality changes along the channel. In BWRs, the cooling rate is calculated using the enthalpy rise, which accounts for both sensible and latent heat, and is often expressed as Q = ṁ [ (1 - x) c_p ΔT + x h_fg ], where x is the quality.

The local heat flux q'' (kW/m²) on a fuel rod is related to the linear heat rate (q') and the rod surface area. The cooling rate for a single channel is Q_channel = ṁ c_p ΔT, and the total power in that channel is the sum of linear heat rates of all rods in the channel. If the total Q_channel is fixed by the channel power, the ΔT is determined. The maximum ΔT that can be tolerated without exceeding the DNB limit defines the maximum allowable linear heat rate. The formula links the global cooling rate to the local thermal limits, ensuring that the hottest rod stays within the DNB ratio (DNBR) limits.

If a coolant has a higher c_p (or higher volumetric heat capacity), then for the same Q and ΔT, the required mass flow rate ṁ = Q / (c_p ΔT) would be lower. For example, sodium has a c_p of about 1.3 kJ/kg·K, which is lower than water's 5.5 kJ/kg·K, so sodium requires a much higher mass flow for the same ΔT. However, sodium's density is lower, and it operates at higher temperatures, so the volumetric flow and pump design differ. The formula is universal: a lower c_p means a higher flow rate to remove the same heat, impacting the pump size and the system's hydraulic design. This is why liquid metal coolants often require high-flow pumps.

The formula Q = ṁ c_p ΔT gives the total heat removed by the coolant, which includes heat transferred from the fuel rods plus any direct heating of the coolant by radiation (gamma and fast neutron absorption). In a typical LWR, about 2-3% of the total heat is deposited directly in the coolant by radiation. This contribution is usually included in the effective heat source for the coolant. In detailed calculations, the heat generation in the coolant is added as a volume source term, and the total ΔT includes both the heat from the fuel and the direct radiation heating. The formula itself is a heat balance; it doesn't distinguish the source, but the input Q must include all heat sources in the coolant.

In a LOCA, when the core becomes partially uncovered, the coolant is no longer a single-phase liquid, and the flow is two-phase or even steam-only. The simple single-phase formula is replaced by a two-phase heat transfer model that accounts for film boiling, droplet entrainment, and reflooding. The heat removal is then Q = ṁ_gas (h_g_out - h_g_in) + ṁ_liquid (h_l_out - h_l_in), with complex heat transfer coefficients. The temperature rise ΔT for the steam is much larger than for liquid because steam's c_p is lower. The formula is adapted to the local flow regime and is used in LOCA analysis codes (e.g., TRACE, RELAP) to predict cladding temperature evolution.

For a fixed mass flow and ΔT, the cooling rate Q = ṁ c_p ΔT is independent of pressure because c_p is a function of temperature (and weakly of pressure). However, if the pressure changes enough to change the state of the coolant (e.g., approaching saturation), the specific heat can change. For subcooled water, a 10% pressure change at constant temperature changes c_p by less than 1%, so the effect is negligible. But if the pressure drops below the saturation pressure, boiling begins, and the two-phase enthalpy rise must be used instead of the single-phase formula. Thus, pressure does not directly appear in the single-phase formula, but it determines the phase state, which dictates which formula applies.

The total required flow rate ṁ_total is calculated from Q_rector = ṁ_total c_p ΔT, where Q_reactor is the core thermal power. This total flow is then divided among the coolant loops (typically 2-4 loops in a PWR). Each loop's flow determines the size of the primary coolant pumps and the steam generator heat transfer area. The ΔT is chosen to balance the steam generator efficiency and the pump power. For example, a 3000 MWth PWR with 4 loops might have ṁ_per_loop ≈ 4,000 kg/s. The steam generator is then designed to transfer heat from that flow rate at the specified primary inlet/outlet temperatures. The formula is the starting point for the entire primary system sizing.

The reactor thermal power is often inferred from Q = ṁ c_p ΔT, using calibrated flow meters and RTDs. The relative uncertainty in Q is the sum of the uncertainties in ṁ, c_p, and ΔT. Since ΔT is only about 40-50 K, a 0.5 K error in ΔT (typical for high-quality RTDs) translates to a 1% error in power. For a 3000 MWth plant, 1% is 30 MW, which is significant for safety analysis and for verifying the reactor's thermal limits. This is why the primary coolant temperature measurement system is highly redundant and calibrated, and why the heat balance is a key part of the plant's startup and periodic testing program.