PWR Pressurizer Surge Line Thermal Stratification Analysis

Fundamentals of Thermal Stratification in PWR Surge Lines

The pressurizer surge line connects the Pressurizer Surge Line to the primary reactor coolant loop hot leg in a Pressurized Water Reactor (PWR). Its main job is to transfer water back and forth as the primary coolant expands or contracts. When flow velocities inside this pipe drop to low speeds on the order of centimeters per second (cm/s), fluid mixing stops. Thermal Stratification occurs when fluids of different temperatures and densities flow into horizontal piping without significant mixing. The hotter, less dense water floats on top while the colder, denser water settles at the bottom, establishing a clear cross-sectional temperature gradient across the pipe.

This dynamic is not limited to emergency situations. Thermal stratification occurs during normal operating conditions such as low flow, heat-up, and cool-down, not exclusively during abnormal or emergency operational transients. During heat-up and cool-down transients in 900 MW PWRs, the maximum temperature difference across surge lines can reach up to 110°C. In operating nuclear plants, direct measurements of fluid temperature differences between the hot and cold water layers typically range from 120-150°C, while surface temperature measuring points on the pipe outer wall have recorded maximum cross-sectional differences reaching 130-140°C. The design temperature difference for thermal stratification in the Korean Standard Nuclear Power Plant (KSNP) surge line is set even higher at 178°C.

Because the top half of the pipe expands more than the cooler bottom half, differential thermal expansion induces severe structural stress. Stratification produces significant circumferential temperature gradients, causing local through-wall axial stresses and global bending stresses throughout the piping system, rather than pure axial stress alone. Over time, these cyclic forces lead to Thermal Fatigue, causing physical deformation, vertical movements of up to 4 inches on 14-28 inch diameter piping configurations, and structural cracking.

Fluid Dynamics and Mixing Layer Behavior During Out-surge and In-surge

Fluid movement through the surge line changes direction based on primary loop pressure changes. An Insurge occurs when cooler water enters the pressurizer from the hot leg during primary coolant volume increases. An Outsurge occurs when hotter water exits the pressurizer into the hot leg during primary system depressurization. Because outsurges push hot fluid over cold, and insurges push cold fluid underneath hot, both transients continuously alter the height and shape of the fluid interface.

To determine whether buoyancy will force the fluid into stratified layers or whether inertial forces will mix it, engineers evaluate the Richardson Number (Ri). The dimensionless ratio is defined as:

$Ri = g · β · Δ T · L/U²$

Here, g is gravitational acceleration, β is the volumetric thermal expansion coefficient, Δ T is the temperature difference between fluid layers, L is the characteristic length (typically the inner pipe diameter), and U is the average flow velocity. A higher Richardson number indicates that natural convection buoyancy forces dominate forced convection inertial forces, driving the fluid to separate into distinct thermal layers.

Computational Fluid Dynamics (CFD) models this transient, dynamic phenomenon where temperature profiles, boundary layer mixing, and spatial interfaces shift continuously over time. CFD codes solve the general transport equation across spatial domains:

$partial(ρ Phi)/partial t + div(ρ U Phi) = div(Gamma grad Phi) + S_Phi$

Where ρ is fluid density, Phi represents mass, momentum, or energy variables, Gamma is the diffusion coefficient, and S_Phi is the source term. At the interface between hot and cold layers, high fluid velocity gradients can induce localized instability called Thermal Striping. Increasing fluid flow velocity does not always eliminate thermal fatigue risks, as higher velocities can aggravate thermal striping loads and exacerbate high-cycle fatigue.

Conjugate Heat Transfer and Pipe Wall Temperature Gradients

The fluid temperature profile does not instantly transfer into the metal structure. Conjugate Heat Transfer (CHT) governs the coupled thermal dynamics between fluid convection inside the pipe boundary layer and heat conduction within the solid pipe wall itself. Thermal stratification significantly impacts pipe wall surface temperatures, direct wall gradients, and long-term fatigue life, rather than affecting fluid temperature alone.

To evaluate the localized structural response, you can compute a simplified baseline inner thermal stress using the classic linear thermoelastic formula:

$σ_thermal = E · α · Δ T_wall$

Where E is Young's Modulus, α is the coefficient of thermal expansion, and Δ T_wall is the temperature differential across or within the wall section.

Spatial geometry strongly influences these heat transfer rates. For example, in marine reactor surge lines, maximum fluid temperature differences have been measured at 59.6 K in horizontal sections and 23.9 K in vertical sections. Consequently, thermal stratification in marine reactor surge lines causes maximum stress fluctuations of 36 MPa in horizontal sections, compared to only 5 MPa in vertical sections.

Structural Stress Distribution and Thermal Fatigue Calculation Methods

Calculating the true structural stress caused by thermal stratification requires looking at both local and global effects simultaneously. When the top half of a surge line is significantly hotter than the bottom half, the upper portion tries to expand longitudinally while the lower portion stays fixed. Because the pipe is constrained at its end connections, namely the pressurizer nozzle and the hot leg nozzle, this differential expansion forces the pipe to arch upward. This structural distortion, known as pipe bananaing, induces large global bending moments across the entire piping system alongside localized through-wall axial stresses.

To compute the cumulative damage from these dynamic loads, plant engineers use specialized fatigue usage factor methodologies, such as two-dimensional to one-dimensional (2D-1D) stress computation frameworks developed by FRAMATOME. These frameworks take transient temperature data from boundary layers and map it directly into finite element structural models. The structural analysis evaluates localized wall stresses caused by stratification-induced bending moments and local non-linear wall temperature profiles. By summing the transient stress ranges over thousands of operating cycles, engineers calculate the Cumulative Usage Factor (CUF). If the total CUF approaches 1.0, the fatigue life of the material is depleted, indicating a high risk of micro-cracking and structural failure.

Mitigation Strategies and In-Service Monitoring Data Integration

Managing thermal stratification in nuclear power plants requires a combination of smart mechanical design and continuous real-time monitoring. Mitigation and management techniques include adjusting piping layout and inclination angle, alongside installing online thermocouple monitoring systems to calculate fatigue usage factors continuously. By adding a continuous slight slope to horizontal pipe runs, buoyancy forces encourage natural mixing rather than sharp thermal separation. Advanced reactor designs, such as the Westinghouse AP1000, modify surge line geometry and layout specifically to minimize buoyancy-driven stratification risks by eliminating long, uninclined horizontal pipe sections near low-flow zones.

For operating units built with legacy layouts, utilities rely on skin temperature monitoring. Arrays of thermocouples are clamped around the outer circumference of the surge line at multiple axial stations. The real-time surface temperature data feeds directly into monitoring software, converting temperature differences into ongoing fatigue usage calculations. This data integration allows plant operators to adjust operational procedures, such as controlling heat-up rates or adjusting spray valve operation during insurge and outsurge events, extending the operational life of the primary coolant system.

Questions students ask

Does thermal stratification only occur during emergency plant trips?

No. Thermal stratification occurs during normal operating conditions such as low flow, heat-up, and cool-down transients, and is not exclusive to emergency operational events.

What fluid flow velocities are typically associated with thermal stratification in surge lines?

Thermal stratification typically occurs when fluid flow velocities inside the surge line drop to low speeds on the order of centimeters per second (cm/s).

How large can temperature differences get across operating surge lines?

Fluid temperature differences between hot and cold layers reach 120-150°C, while maximum surface temperature differences across the pipe outer wall reach 130-140°C in operating plants. Design values can reach 178°C.

Can increasing fluid flow velocity completely eliminate thermal fatigue risks?

No. Increasing fluid flow velocity does not always eliminate fatigue risks because higher velocities can intensify thermal striping loads at the fluid interface, aggravating high-cycle thermal fatigue.

How does thermal stratification affect global piping displacement?

Differential thermal expansion causes structural bending (bananaing), which can produce global vertical displacements of up to 4 inches on 14-28 inch diameter piping systems.