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.

What is Drag in Physics

In fluid dynamics, drag is a force acting opposite to the relative motion of any moving object. The force a flowing fluid exerts on a body in the flow direction. Unlike other resistive forces, such as dry friction, which are nearly independent of velocity, drag forces depend on velocity. Drag force is proportional to the velocity for a laminar flow and the squared velocity for a turbulent flow. Drag is generally caused by two phenomena:

  • Drag - skin friction and form drag

    Source: wikipedia.org License: CC BY-SA 3.0

    Skin Friction.  In general, when a fluid flows over a stationary surface, e.g. the flat plate, the bed of a river, or the wall of a pipe, the fluid touching the surface is brought to rest by the shear stress to at the wall. The region in which flow adjusts from zero velocity at the wall to a maximum in the main stream of the flow is termed the boundary layer. Therefore, a moving fluid exerts tangential shear forces on the surface because of the no-slip condition caused by viscous effects.  This type of drag force, depends especially on the geometry, the roughness of the solid surface and on the type of fluid flow.

  • Form Drag. Form drag known also as pressure drag arises because of the shape and size of the object. This type of drag force is an interesting consequence the Bernoulli’s effect. According to Bernoulli’s principle, faster moving air exerts less pressure. This causes, that there can be a pressure difference between surfaces of the object. The general size and shape of the body are the most important factors in form drag. In general, bodies with a larger presented geometric cross-section will have a higher drag than thinner bodies.

Both of these forces, in general, have components in the direction of flow, and thus the resulting drag force is due to the combined effects of pressure and skin friction forces in the flow direction.

When the friction and pressure drag coefficients are available, the total drag coefficient is determined by simply adding them:

skin friction - form drag - coefficients

At low Reynolds numbers, most drag is due to friction drag. This is especially the case for highly streamlined bodies such as airfoils. On the other hand, at high Reynolds number, the pressure drop is significant, which increases form drag.

The components of the pressure and skin friction forces in the normal direction to flow tend to move the body in that direction, and their sum is called lift.

In aeronautics, the lift is an upward-acting force on an aircraft wing or airfoil. Bernoulli’s principle requires airfoil to be of an asymmetrical shape.

Drag Force – Drag Equation

The drag force, FD,depends on the density of the fluid, the upstream velocity, and the size, shape, and orientation of the body, among other things. One way to express this is by means of the drag equation. The drag equation is a formula used to calculate the drag force experienced by an object due to movement through a fluid.

drag force - drag equation - formula

The reference area, A, is defined as the area of the orthographic projection of the object on a plane perpendicular to the direction of motion. For hollow objects, the reference area may be significantly larger than the cross sectional area, but for non-hollow objects, it is exactly the same as a cross sectional area.

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: