## Conservation of Mass

The mass can neither be created nor destroyed.

Continuity Equation – Definition

This principle is generally known as the conservation of matter principle and states that the mass of an object or collection of objects never changes over time, no matter how the constituent parts rearrange themselves. This principle can be use in the analysis of flowing fluids. Conservation of mass in fluid dynamics states that all mass flow rates into a control volume are equal to all mass flow rates out of the control volume plus the rate of change of mass within the control volume. This principle is expressed mathematically by following equation:

in = ṁout +∆m∆t

Mass entering per unit time = Mass leaving per unit time + Increase of mass in the control volume per unit time

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

This equation describes nonsteady-state flow. Nonsteady-state flow refers to the condition where the fluid properties at any single point in the system may change over time. Steady-state flow refers to the condition where the fluid properties (temperature, pressure, and velocity) at any single point in the system do not change over time. But one of the most significant properties that is constant in a steady-state flow system is the system mass flow rate. This means that there is no accumulation of mass within any component in the system.

## Continuity Equation

The continuity equation is simply a mathematical expression of the principle of conservation of mass. For a control volume that has a single inlet and a single outlet, the principle of conservation of mass states that, for steady-state flow, the mass flow rate into the volume must equal the mass flow rate out.

in = ṁout

Mass entering per unit time = Mass leaving per unit time

This equation is called the continuity equation for steady one-dimensional flow. For a steady flow through a control volume with many inlets and outlets, the net mass flow must be zero, where inflows are negative and outflows are positive.

This principle can be applied to a streamtube such as that shown above. No fluid flows across the boundary made by the streamlines so mass only enters and leaves through the two ends of this streamtube section.

When a fluid is in motion, it must move in such a way that mass is conserved. To see how mass conservation places restrictions on the velocity field, consider the steady flow of fluid through a duct (that is, the inlet and outlet flows do not vary with time).

## Differential Form of Continuity Equation

A general continuity equation can also be written in a differential form:

∂⍴∂t + ∇ . (⍴ ͞v) = σ

where

• ∇ . is divergence,
• ρ is the density of quantity q,
• ⍴ ͞v is the flux of quantity q,
• σ is the generation of q per unit volume per unit time. Terms that generate (σ > 0) or remove (σ < 0) q are referred to as a “sources” and “sinks” respectively. If q is a conserved quantity (such as energy), σ is equal to 0.

## Continuity Equation – Multiple Inlets and Outlets

For a control volume with multiple inlets and outlets, the principle of conservation of mass requires that the sum of the mass flow rates into the control volume equal the sum of the mass flow rates out of the control volume. The continuity equation for this more general situation is expressed by following equation:

∑ṁin = ∑ṁout

Sum of mass flow rates entering per unit time = Sum of mass flow rates leaving per unit time

## Continuity Equation – Examples

### Flow rate through a reactor core

In this example, we will calculate the flow rate through a reactor core from continuity equation. It is an illustrative example, following data do not represent any reactor design.

in = ṁout

(ρAv)in = (ρAv)out

____________________________

Density as a function of temperature of water

Pressurized water reactors are cooled and moderated by high-pressure liquid water (e.g. 16MPa). At this pressure water boils at approximately 350°C (662°F).  Inlet temperature of the water is about 290°C (⍴ ~ 720 kg/m3). The water (coolant) is heated in the reactor core to approximately 325°C (⍴ ~ 654 kg/m3) as the water flows through the core.
The primary circuit of typical PWR is divided into 4 independent loops (piping diameter ~ 700mm), each loop comprises a steam generator and one main coolant pump. Inside the reactor pressure vessel (RPV), the coolant first flows down outside the reactor core (through the downcomer). From the bottom of the pressure vessel, the flow is reversed up through the core, where the coolant temperature increases as it passes through the fuel rods and the assemblies formed by them.

### Calculate:

• the primary piping volumetric flow rate (m3/s),
• the primary piping flow velocity (m/s),
• the core inlet flow velocity (m/s),
• the core outlet flow velocity (m/s)

when

• the mass flow rate in the hot leg of primary piping is equal to 4648 kg/s,
• Reactor core flow cross-section is equal to 5m2,
• Primary piping flow cross-section (single loop) is equal to 0.38 m2

### Results:

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

Cold leg volumetric flow rate:

Qcold = ṁ / ⍴ = 4648 / 720 = 6.46 m3/s = 23240 m3/hod

Cold leg flow velocity:

A1 = π.d2 / 4

vcold = Qcold / A1 = 6.46 / (3.14 x 0.72 / 4) = 6.46 / 0.38 = 17 m/s

Hot leg volumetric flow rate:

Qhot = ṁ / ⍴ = 4648 / 654 = 7.11 m3/s = 25585 m3/hod

Hot leg flow velocity:

A = π.d2 / 4

vhot = Qhot / A1 = 7.11 / (3.14 x 0.72 / 4) = 7.11 / 0.38 = 18,7 m/s

or according to the continuity equation:

1 . A1 . v1 = ⍴2 . A2 . v2

vhot =  vcold . ⍴cold / ⍴hot = 17 x 720 / 654 = 18.7 m/s

Core inlet flow velocity:

Acore = 5m2

Apiping = 4 x A1 = 4 x 0.38 = 1.52 m2

inlet = ⍴cold

according to the continuity equation:

inlet . Acore . vinlet = ⍴cold . Apiping . vcold

vinlet =  vcold . Apiping / Acore = 17 x 1.52 / 5 = 5.17 m/s

Core outlet flow velocity:

inlet = ⍴cold

outlet = ⍴hot

according to the continuity equation:

outlet . Acore . voutlet = ⍴inlet . Acore . vinlet
voutlet =  vinlet . ⍴inlet / ⍴outlet = 5.17 x 720 / 654 = 5.69 m/s

#### References:

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.