The fundamental governing equations are the incompressible continuity equation and the Navier–Stokes momentum equations. For an isothermal, incompressible Newtonian fluid with constant density (\rho) and dynamic viscosity (\mu), they are

[ \nabla \cdot \mathbf{v}=0 ]

and

[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} +\mathbf{v}\cdot\nabla\mathbf{v}\right)

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

Together, these equations determine the three velocity components and pressure field, subject to the appropriate equipment geometry, inlet, outlet, wall, and initial conditions.

Core takeaway: Isothermal, incompressible Newtonian flow in pilot-plant equipment is modeled primarily by mass conservation and momentum conservation. Energy, species, reaction, or multiphase equations are added only when the unit operation includes effects that cannot be represented by these two equations alone.

The Governing Equation Set

Mass conservation: the continuity equation

The general mass conservation equation is

[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0. ]

For an incompressible fluid with constant density,

[ \boxed{\nabla\cdot\mathbf{v}=0} ]

This condition states that a fluid element does not expand or contract as it moves. In Cartesian coordinates,

[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0, ]

where (u), (v), and (w) are the velocity components in the (x), (y), and (z) directions.

Momentum conservation: the Navier–Stokes equations

The linear momentum balance is

[ \boxed{ \rho\frac{D\mathbf{v}}{Dt}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g} } ]

where the material derivative is

[ \frac{D\mathbf{v}}{Dt}

\frac{\partial\mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}. ]

Thus, the momentum equation represents a balance among:

Fluid acceleration: (\rho D\mathbf{v}/Dt)
Pressure forces: (-\nabla p)
Viscous forces: (\mu\nabla^2\mathbf{v})
Body forces: (\rho\mathbf{g}), usually gravity

The convective term ((\mathbf{v}\cdot\nabla)\mathbf{v}) is important in bends, contractions, expansions, mixers, and other regions where velocity varies spatially.

Cartesian component form

For velocity (\mathbf{v}=(u,v,w)), the momentum equations are

[ \rho\left( \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z} \right)

-\frac{\partial p}{\partial x} +\mu\nabla^2u+\rho g_x, ]

[ \rho\left( \frac{\partial v}{\partial t} +u\frac{\partial v}{\partial x} +v\frac{\partial v}{\partial y} +w\frac{\partial v}{\partial z} \right)

-\frac{\partial p}{\partial y} +\mu\nabla^2v+\rho g_y, ]

[ \rho\left( \frac{\partial w}{\partial t} +u\frac{\partial w}{\partial x} +v\frac{\partial w}{\partial y} +w\frac{\partial w}{\partial z} \right)

-\frac{\partial p}{\partial z} +\mu\nabla^2w+\rho g_z. ]

These three equations, combined with continuity, give four equations for four primary fields: (u), (v), (w), and (p).

Why Newtonian Behavior Matters

Newton’s law of viscosity

A Newtonian fluid has viscous stresses that vary linearly with the rate of deformation. In tensor form,

[ \boldsymbol{\tau}

2\mu\mathbf{D}, ]

for incompressible flow, where

[ \mathbf{D}

\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right] ]

is the rate-of-deformation tensor.

For constant viscosity and incompressible flow, the divergence of the viscous stress reduces to

[ \nabla\cdot\boldsymbol{\tau}

\mu\nabla^2\mathbf{v}. ]

That reduction is what produces the viscous term in the standard Navier–Stokes equation.

Physical interpretation

The viscous term describes momentum diffusion. It accounts for the transfer of momentum caused by velocity gradients, such as those between a stationary pipe wall and the moving fluid in the center of the pipe.

Conditions Needed to Close the Model

No-slip wall condition

At a stationary solid wall,

[ \mathbf{v}=0. ]

This condition is central to predicting wall shear stress, pressure loss, and velocity profiles in pipes, channels, heat exchangers, and packed equipment.

For a moving wall, the fluid velocity at the wall equals the wall velocity.

Inlet conditions

An inlet may be specified by a velocity profile,

[ \mathbf{v}=\mathbf{v}_{\text{in}}, ]

or by a mass flow rate,

[ \dot{m}

\int_A \rho,\mathbf{v}\cdot\mathbf{n},dA. ]

A uniform inlet profile is often used for simplified studies, while developed or measured profiles are more appropriate for detailed pilot-plant analysis.

Outlet conditions

Typical outlet specifications include a known pressure,

[ p=p_{\text{out}}, ]

or a suitable outflow condition when the downstream flow is sufficiently developed. The selected outlet condition can materially affect predicted pressure fields and recirculation.

Initial condition

For transient operation, the starting velocity field must be defined:

[ \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}). ]

A pressure reference is also required because incompressible pressure is generally determined only up to an arbitrary constant.

Common Simplifications in Pilot-Plant Modules

Steady flow

If operating conditions do not vary significantly with time,

[ \frac{\partial\mathbf{v}}{\partial t}=0. ]

The momentum equation becomes

[ \rho(\mathbf{v}\cdot\nabla)\mathbf{v}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

This approximation is commonly used for steady pipe, channel, and heat-exchanger analyses.

Fully developed internal flow

For fully developed flow in a straight channel or pipe, the velocity profile does not change in the axial direction. The governing equations can then reduce to a balance between pressure gradient and viscous resistance.

For example, for one-dimensional flow between parallel plates,

[ 0=-\frac{dp}{dx}+\mu\frac{d^2u}{dy^2}. ]

The resulting velocity profile is parabolic when the pressure gradient is constant.

Negligible inertia

At sufficiently low Reynolds number, inertial effects are small compared with viscous effects. The momentum equation reduces to the Stokes-flow form,

[ 0=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

This approximation is useful for highly viscous fluids, narrow channels, and low-flow-rate micro or laboratory modules.

Gravity and hydrostatic pressure

If velocity is negligible, the momentum equation reduces to

[ \nabla p=\rho\mathbf{g}. ]

This describes hydrostatic pressure variation. In flowing systems, gravity may remain important in vertical piping, packed columns, tanks, and natural-convection-related configurations.

How the Equations Apply to Unit Operations

Piping networks and pumps

The equations describe velocity distributions and pressure gradients throughout piping. Integrated or reduced forms are then used to estimate pressure losses and determine pump requirements.

For a pilot plant, the full field equations may be solved computationally, while simpler pipe sections are often represented through friction-factor correlations derived from these balances.

Heat exchangers

For strictly isothermal fluid-flow modeling, continuity and momentum determine the hydrodynamics, including velocity distribution, wall shear, and pressure drop.

If the objective includes temperature change or heat-transfer duty, an energy equation must also be added. A common stream-level balance is

[ Q=\dot{m}C_p\Delta T, ]

while exchanger performance is often related through

[ Q=UA\Delta T_{\mathrm{lm}}F. ]

These are additional thermal relations, not replacements for the flow equations.

Reactors and mixing equipment

In stirred tanks and flow reactors, the Navier–Stokes equations predict circulation, recirculation, and mixing-related velocity fields.

Reaction, species transport, and residence-time behavior require additional conservation equations. The fluid-flow equations alone do not predict conversion or concentration.

Packed beds and columns

In packed equipment, the pore-scale Navier–Stokes equations can be used for detailed modeling, but engineering-scale models commonly replace them with porous-media pressure-drop relations.

The underlying principles remain mass and momentum conservation, with the packing represented through additional resistance terms.

Two-fluid or multiphase modules

For two immiscible fluids, each phase requires its own conservation equations, along with interface conditions. At a liquid–liquid interface, typical conditions include

[ \mathbf{v}_1=\mathbf{v}_2 ]

for tangential velocity continuity and

[ \boldsymbol{\tau}_1\cdot\mathbf{n}

\boldsymbol{\tau}_2\cdot\mathbf{n} ]

for viscous stress continuity, when interfacial effects are treated appropriately.

The single-fluid incompressible equations are therefore insufficient when phase boundaries, interphase momentum transfer, or phase change are important.

Understanding the Trade-offs

Isothermal does not mean energy is always irrelevant

An isothermal model assumes that temperature is fixed or that thermal effects do not influence the flow solution. It does not provide temperature distributions, heat-transfer coefficients, or thermal duty.

For heat exchangers, heated reactors, or systems with strong viscosity variation, an energy equation may be necessary even if the initial flow analysis uses an isothermal approximation.

Incompressibility has a defined range of validity

Liquids are often modeled as incompressible because density changes are small. Gases may also be treated as incompressible when velocity is sufficiently below the speed of sound and pressure-induced density changes are negligible.

The approximation should be reconsidered for high-speed gas flow, large pressure changes, strong temperature variation, or systems involving compressibility-driven effects.

Constant viscosity may be a poor approximation

The standard equation assumes constant (\mu). This is reasonable for many liquids over modest temperature and shear-rate ranges, but not for strongly temperature-dependent fluids or non-Newtonian materials.

If viscosity varies spatially, the viscous term should generally be retained in stress-divergence form rather than written simply as (\mu\nabla^2\mathbf{v}).

Laminar assumptions are not required by Navier–Stokes

The Navier–Stokes equations apply to Newtonian flow whether it is laminar or turbulent. However, direct resolution of turbulent scales is computationally demanding, so turbulence models or experimentally validated correlations are often required for pilot-scale equipment.

“Viscous dissipation” belongs to the energy balance

Viscous stresses affect the momentum equation directly. The heat generated by viscous dissipation appears in the energy equation and is omitted from a strictly isothermal formulation.

This distinction prevents the common mistake of treating viscous dissipation as an additional term in the momentum equation.

Making the Right Choice for Your Goal

The correct governing model depends on what the pilot-plant analysis must predict.

If your primary focus is velocity profiles and pressure loss: Use incompressible continuity plus the Navier–Stokes momentum equations, with density and viscosity appropriate to the operating conditions.
If your primary focus is heat-transfer performance: Add an energy equation and use the flow solution to determine convection, pressure drop, and heat-transfer behavior.
If your primary focus is reactor conversion or species distribution: Add species or reaction balances to the hydrodynamic equations.
If your primary focus is liquid–liquid or gas–liquid operation: Use phase-specific conservation equations and include interfacial velocity, stress, and mass-transfer conditions.
If your primary focus is rapid equipment-scale design: Use validated reduced models, such as friction-factor, porous-media, or residence-time correlations, while recognizing their operating limits.

With the assumptions clearly stated and the boundary conditions correctly specified, these conservation equations provide a reliable foundation for analyzing pilot-plant fluid-flow behavior.

Summary Table:

Model element	Governing relation	Primary purpose
Mass conservation	∇ · v = 0	Ensures incompressible flow continuity
Momentum conservation	ρ(∂v/∂t + v · ∇v) = −∇p + μ∇²v + ρg	Predicts velocity, pressure, viscous, and body-force effects
Newtonian viscosity	τ = 2μD	Relates viscous stress to deformation rate
Boundary conditions	Inlet, outlet, wall, and initial conditions	Closes the mathematical model
Additional physics	Energy, species, reaction, or multiphase equations	Extends the model for heat transfer, reactions, or multiple phases

Need reliable pilot-plant equipment for fluid-flow education, process development, or scale-up? Contact LABPARK to discuss your requirements. LABPARK provides Educational and Vocational Unit Operations Pilot Plants for chemical engineering, bioprocess and biotechnology, and environmental and water-treatment applications. Our modular systems help universities, research institutes, and enterprises study pressure drop, velocity profiles, mixing, heat transfer, separation, reaction, and process control with practical, application-focused equipment. Tell us your target process, operating conditions, teaching or research objectives, and required scale so our team can recommend a suitable configuration and support your pilot-plant planning.