Reactor Heat Removal Rate

The core is not uniform—different fuel assemblies have different power levels, coolant flow rates, and inlet temperature distributions due to radial and axial power peaking. Summing over each channel (or subchannel) accounts for these local variations, which is critical for identifying hot spots. A single average value would underestimate the peak temperatures and might miss a local boiling crisis. Thermal-hydraulic codes divide the core into hundreds or thousands of channels to capture the detailed thermal margins.

In a typical PWR, the coolant enters the core at about 280°C and exits at 320-330°C, giving a total core average ΔT of about 40-50 K across the core. However, individual channels can have higher ΔT due to lower flow or higher power. For example, a hot channel near the center might see a ΔT of 60-70 K, while peripheral channels with higher flow might see only 30 K. The formula uses the local ΔT for each channel, and the sum over all channels gives the total heat removal.

The specific heat of water (or any coolant) varies with temperature and pressure. At typical PWR conditions (≈15 MPa, 300°C), cp is about 5.5-5.8 kJ/(kg·K) for liquid water. In BWRs at lower pressure, cp is slightly different. Using a constant average cp can introduce errors, especially when ΔT is large or near the saturation point. Accurate calculations use cp as a function of temperature (and pressure) from steam tables. For safety analyses, designers use the minimum or maximum cp in the expected range to bound the heat removal uncertainty.

In a steady-state operation, the heat removal rate must exactly equal the thermal power generated in the core (minus a small amount of heat stored in the fuel and structure). If the calculated Q_rem is greater than the core power, it indicates either the flow is higher than expected, the temperature rise is overestimated, or the power measurement is incorrect. If Q_rem is less than core power, the core would heat up. During transients, the heat removal rate may be less than power, causing temperature rise; during cooldown, it may exceed power. The formula is used to verify the energy balance.

The flow per channel is determined by the core inlet plenum pressure, the pressure drop through each fuel assembly (which depends on the assembly geometry, spacer grids, and the local flow conditions), and the core outlet pressure. Flow variations arise from differences in: (1) assembly hydraulic resistance (due to manufacturing tolerances or debris), (2) local power peaking causing different boiling (in BWRs), (3) control rod insertion changing the flow area, (4) thermal-hydraulic feedback (density changes), and (5) pump characteristics. In design, flow distribution is calculated using computational fluid dynamics (CFD) and validated with in-core flow measurements. The formula sums over the actual flow distribution, not just the average flow.

In forced convection, the coolant is driven by primary pumps, ensuring high flow and large heat removal capacity. The formula Q_rem = Σ (flow × cp × ΔT) still applies, but the flow is set by pump speed. In natural convection (e.g., after a loss of forced cooling), the flow is driven by density gradients: hot coolant rises, cold coolant sinks. The flow is much lower and the ΔT is larger. The heat removal rate is limited by the buoyancy-driven flow and is used in accident analysis. The same formula applies, but the flow and ΔT are determined by the buoyancy balance, not pumps.

The formula shows that if the inlet temperature (and thus ΔT) changes, the heat removal rate changes if flow is constant. A higher inlet temperature reduces ΔT for a given outlet temperature, lowering heat removal, which may require adjusting control rods or flow to maintain power. In load-following, the turbine demand changes, and the reactor control system adjusts both the average coolant temperature and flow to match the desired heat removal. The formula is fundamental to designing the control algorithms that keep the reactor balanced with the turbine load.

During a loss-of-coolant accident (LOCA), the ECCS must inject coolant to remove decay heat. The formula Q_rem = Σ (flow × cp × ΔT) is applied to the transient conditions: the flow from the ECCS pumps, the temperature of the injected water, and the ΔT of the coolant as it passes through the core. Safety analysis uses this to ensure that the heat removal rate always exceeds the decay heat generation, preventing core melt. The sum over channels ensures that even the hottest channels receive sufficient cooling. The formula is used in the safety analysis report (SAR) to demonstrate compliance with regulatory criteria.

Yes, but with care. In supercritical water reactors (SCWRs), the coolant operates above the critical point (≈22.1 MPa, 374°C). The specific heat cp peaks sharply near the pseudo-critical temperature (the temperature of maximum cp at a given pressure). This makes the integral of cp over ΔT nonlinear; using a constant cp would be inaccurate. The formula is applied in an integrated form: Q_rem = Σ flow × ∫_{T_in}^{T_out} cp(T) dT, effectively calculating the enthalpy rise. Many thermal-hydraulic codes for SCWRs use enthalpy-based formulations rather than simple cp × ΔT.

The heat removal rate per channel determines the local heat flux (power per fuel rod surface area). The DNB ratio (or CHF ratio) is the margin between the actual heat flux and the CHF limit at which a vapor film forms on the cladding, leading to a sharp temperature rise. The formula Q_rem = flow × cp × ΔT gives the total heat removed, but it must be combined with the heat flux distribution to check that each channel's local heat flux does not exceed the CHF. Thermal-hydraulic analysis uses the ΔT and flow to calculate the channel's enthalpy rise and then evaluates CHF using empirical correlations. A low DNB ratio indicates a high risk of dryout.

In high-fidelity CFD models, the core is not divided into discrete channels but is represented as a continuous porous medium with distributed flow and heat sources. The summation is replaced by an integral over the core volume: Q_rem = ∫∫∫ (ρ v · ∇h) dV, where ρ is density, v is velocity, and h is enthalpy. This approach captures 3D flow and temperature distributions, including cross-flow between assemblies, with no need to manually sum channels. However, for licensing and simpler safety analyses, the discrete channel summation is still used because it's more transparent and easier to validate with measured data.

In a BWR, boiling occurs in the core, so the coolant is a two-phase mixture. The formula using cp × ΔT for liquid water is not directly applicable because the enthalpy change includes the latent heat of vaporization. Instead, the heat removal rate is calculated as Q_rem = Σ (flow × (h_out - h_in)), where h is the specific enthalpy (including sensible and latent heat). The mass flow rate is the total flow, and the exit quality is determined by the heat input. The sum over channels accounts for different void fractions and flow distributions. This enthalpy-based form is more general and is used for two-phase systems.

The hot channel factor is the ratio of the power in the hottest fuel assembly (or rod) to the core average. In the heat removal sum, the hot channel will have a higher heat flux, requiring either higher flow or a larger ΔT to remove that heat. The formula is used to check that the hot channel's ΔT does not exceed the limit for DNB or cladding temperature. The hot channel factor is a key input to the thermal-hydraulic design: it determines the required margin in the cooling system. The sum over all channels uses the actual power distribution, not an average.

Yes. The total heat removal rate Q_rem must equal the reactor's thermal power. This determines the required total coolant flow rate: total flow = Q_rem / (cp × ΔT_avg). This total flow sets the pump capacity. The steam generators must be able to transfer this heat from the primary to the secondary side, so their heat transfer area is sized to handle Q_rem at the design primary and secondary temperatures. The formula is fundamental to matching the core, the pumps, and the steam generators in the balance of plant design. During plant design, the sum over channels is used to calculate the required flow distribution that meets all thermal limits.

In a LOFA, the primary pump flow decreases, reducing the flow in each channel. The formula shows that if flow decreases while power remains constant (until the reactor trips), the ΔT in each channel must increase to remove the same heat. If ΔT exceeds the limit (e.g., approaching saturation or causing DNB), the cladding may overheat. The sum over channels shows which channels reach their limits first—usually the hot channels with the lowest flow-to-power ratio. The safety analysis uses the formula to predict the time to fuel damage and the effectiveness of the reactor trip and emergency cooling system.