Drag Force in Nuclear Engineering

Analysis of the hydraulic lift force is one of the most important analysis in designing a fuel assembly and analyzing the hydraulic compatibility of mixed cores. The vertical forces are induced by upward high-velocity flow through the reactor core. The flow path for the reactor coolant through the reactor vessel would be:

  • Continuity Equation - Flow Rates through Reactor

    Example of flow rates in a reactor. It is an illustrative example, data do not represent any reactor design.

    The coolant enters the reactor vessel at the inlet nozzle and hits against the core barrel.

  • The core barrel forces the water to flow downward in the space between the reactor vessel wall and the core barrel, this space is usually known as the downcomer.
  • From the bottom of the pressure vessel, the flow is reversed up through the core in order to pass through the fuel assemblies, where the coolant temperature increases as it passes through the fuel rods.
  • Finally the hotter reactor coolant enters the upper internals region, where it is routed out the outlet nozzle into the hot legs of primary circuit and goes on to the steam generators.

Fuel assemblies are held by the upper guide structure assembly, which defines the top of the core. This assembly is made of stainless steel and has many purposes. The upper guide structure assembly exerts an axial force on fuel assemblies (through springs in the top nozzle), thus defines the exact position of the fuel assembly in the core. The upper guide structure assembly flange is held in place and preloaded by the RPV closure head flange. The upper guide structure assembly also guides and protects control rod assemblies and in-core instrumentation.

Required downforce of the upper guide structure assembly on fuel assemblies must very carefully calculated. Insufficient downforce can result in the lift of the fuel assembly, on the other hand, an excessive downforce can result in bowing of fuel assembly, which is also unacceptable.

Example: Drag Force – Drag Coefficient – Fuel Bundle

Hydraulic Diameter - Fuel ChannelCalculate the friction drag of a single fuel rod inside a reactor core at normal operation (design flow rate). Assume that this fuel rod is part of a fuel bundle with the rectangular fuel lattice and this fuel bundle does not contain spacing grids. Its height is h = 4m and the core flow velocity is constant and equal to Vcore = 5 m/s.

Assume that:

  • the outer diameter of the cladding is: d = 2 x rZr,1 = 9,3 mm
  • the pitch of fuel pins is: p = 13 mm
  • the relative roughness is ε/D = 5×10-4
  • the fluid density is: ρ = 714 kg/m3
  • the core flow velocity is constant and equal to Vcore = 5 m/s
  • the average temperature of reactor coolant is: Tbulk = 296°C

Calculation of the Reynolds number

To calculate the Reynolds number, we have to know:

  • the outer diameter of the cladding is: d = 2 x rZr,1 = 9,3 mm (to calculate the hydraulic diameter)
  • the pitch of fuel pins is: p = 13 mm  (to calculate the hydraulic diameter)
  • the dynamic viscosity of saturated water at 300°C is: μ = 0.0000859 N.s/m2
  • the fluid density is: ρ = 714 kg/m3

The hydraulic diameter, Dh, is a commonly used term when handling flow in non-circular tubes and channels. The hydraulic diameter of the fuel channel, Dh, is equal to 13,85 mm.

See also: Hydraulic Diameter

The Reynolds number inside the fuel channel is then equal to:

reynolds number - example

This fully satisfies the turbulent conditions.

Calculation of the Skin Friction Coefficient

The friction factor for turbulent flow depends strongly on the relative roughness. It is determined by the Colebrook equation or can be determined using the Moody chart. The Moody chart for Re = 575 600 and ε/D = 5 x 10-4 returns following values:

Therefore the skin friction coefficient is equal to:

skin friction coefficient - example

Calculation of the Drag Force

To calculate the drag force, we have to know:

  • the skin friction coefficient, which is: CD,friction = 0.00425
  • the area of pin surface, which is: A = π.d.h = 0.1169 m2
  • the fluid density, which is: ρ = 714 kg/m3
  • the core flow velocity, which is constant and equal to Vcore = 5 m/s

From the skin friction coefficient, which is equal to the Fanning friction factor we can calculate the frictional component of the drag force. The drag force is given by:

drag force - example

Assuming that a fuel assembly can have, for example, 289 fuel pins (17×17 fuel assembly), the frictional component of the drag force is then of the order of kilonewtons. Moreover, this drag force originates purely from the skin friction on fuel bundle. But typical PWR fuel assembly contains another components, which influence the fuel assembly hydraulics:

  • Fuel rods. Fuel rods contain the fuel and burnable poisons.
  • Top nozzle. Provides the mechanical support for the fuel assembly structure.
  • Bottom nozzle. Provides the mechanical support for the fuel assembly structure.
  • Spacing grid. Ensures an exact guiding of the fuel rods.
  • Guide thimble tube. Vacant tube for control rods or in-core instrumentation.

As was written, the second component of the drag force is the form drag. Form drag known also as pressure drag arises because of the shape and size of the object. The pressure drag is proportional to the difference between the pressures acting on the front and back of the immersed body, and the frontal area.

Pressure Drop – Fuel Assembly

In general, total fuel assembly pressure drop is formed by fuel bundle frictional drop (dependent on relative roughness of fuel rods, reynolds number, hydraulic diameter etc.) and other pressure drops of structural elements (top and bottom nozzle, spacing grids or mixing grids).

In general, it is not so simple to calculate pressure drops in fuel assemblies (especially the spacing grids) and it belongs to key know-how of certain fuel manufacturer. Mostly, pressure drops are measured in experimental hydraulic loops, rather than calculated.

Engineers use the pressure loss coefficient, PLC. It is noted K or ξ  (pronounced “xi”). This coefficient characterizes pressure loss of a certain hydraulic system or of a part of a hydraulic system. It can be easily measured in hydraulic loops. The pressure loss coefficient can be defined or measured for both straight pipes and especially for local (minor) losses.

PLC - Pressure loss coefficient - equations

Using data from above mentioned example, the pressure loss coefficient (only frictional from straight pipe) is equal to ξ = fDL/DH = 4.9. But the overall pressure loss coefficient (including spacing grids, top and bottom nozzles etc.) is usually about three times higher. This PLC (ξ = 4.9) causes that the pressure drop is of the order of (using the previous inputs) Δpfriction =  4.9 x 714 x 52/ 2 = 43.7 kPa (without spacing grids, top and bottom nozzles). About three times higher real PLC means about three times higher Δpfuel will be.

The overall reactor pressure loss, Δpreactor, must include:

  • downcomer and reactor bottom
  • lower support plate
  • fuel assembly including spacing grids, top and bottom nozzles and other structural components – Δpfuel
  • upper guide structure assembly

In result the overall reactor pressure loss – Δpreactor is usually of the order of hundreds kPa (let say 300 – 400 kPa) for design parameters.

Reactor Physics and Thermal Hydraulics:

  1. J. R. Lamarsh, Introduction to Nuclear Reactor Theory, 2nd ed., Addison-Wesley, Reading, MA (1983).
  2. J. R. Lamarsh, A. J. Baratta, Introduction to Nuclear Engineering, 3d ed., Prentice-Hall, 2001, ISBN: 0-201-82498-1.
  3. W. M. Stacey, Nuclear Reactor Physics, John Wiley & Sons, 2001, ISBN: 0- 471-39127-1.
  4. Glasstone, Sesonske. Nuclear Reactor Engineering: Reactor Systems Engineering, Springer; 4th edition, 1994, ISBN: 978-0412985317
  5. Todreas Neil E., Kazimi Mujid S. Nuclear Systems Volume I: Thermal Hydraulic Fundamentals, Second Edition. CRC Press; 2 edition, 2012, ISBN: 978-0415802871
  6. Zohuri B., McDaniel P. Thermodynamics in Nuclear Power Plant Systems. Springer; 2015, ISBN: 978-3-319-13419-2
  7. Moran Michal J., Shapiro Howard N. Fundamentals of Engineering Thermodynamics, Fifth Edition, John Wiley & Sons, 2006, ISBN: 978-0-470-03037-0
  8. Kleinstreuer C. Modern Fluid Dynamics. Springer, 2010, ISBN 978-1-4020-8670-0.
  9. U.S. Department of Energy, THERMODYNAMICS, HEAT TRANSFER, AND FLUID FLOW. DOE Fundamentals Handbook, Volume 1, 2 and 3. June 1992.
  10. White Frank M., Fluid Mechanics, McGraw-Hill Education, 7th edition, February, 2010, ISBN: 978-0077422417

See above: