Fluid Kinematics

Learning Objectives

  • Classify flows using time dependence, spatial variation, dimensionality, compressibility, and rotationality.
  • Distinguish Eulerian field descriptions from Lagrangian particle descriptions and common flow-visualization lines.
  • Compute local, convective, and total material acceleration from a velocity field.
  • Apply conservation of mass to fixed control volumes, streamtubes, junctions, and variable-storage systems.
  • Interpret fluid-element strain, rotation, vorticity, and circulation from velocity gradients.
  • Use stream functions and velocity potentials only under their appropriate dimensional and physical assumptions.
  • Interpret two-dimensional potential-flow and seepage flow nets while respecting boundary and isotropy assumptions.

Lagrangian description

A description that follows identified fluid particles and records their positions, velocities, accelerations, and other properties as functions of time.

Eulerian description

A description that specifies fluid properties as fields at fixed spatial locations, such as V(x,y,z,t)\mathbf{V}(x,y,z,t), and observes particles as they move through those fields.

Eulerian and Lagrangian Viewpoints

Most hydraulic calculations use an Eulerian description because pressure, velocity, discharge, and acceleration are evaluated at sections or locations in space. Lagrangian tracking is valuable when the history of particular particles matters, such as sediment paths, pollutant trajectories, or particle-image velocimetry.

Neither description changes the physics; they are different ways of representing the same motion.

Eulerian Velocity Field

Represents the three Cartesian velocity components as functions of position and time.

V(x,y,z,t)=u i+v j+w k\mathbf{V}(x,y,z,t)=u\,\mathbf{i}+v\,\mathbf{j}+w\,\mathbf{k}

Variables

SymbolDescriptionUnit
u,v,wu,v,wVelocity components in the x, y, and z directionsm/s
ttTimes

Steady flow

Flow for which a property observed at a fixed location does not change with time; for that property, its partial derivative with respect to time is zero.

Uniform flow

Flow for which the selected property is spatially constant along the direction or region being considered at a specified instant.

Common Flow Classifications

  • Steady versus unsteady: Concerns variation with time at a fixed location.
  • Uniform versus nonuniform: Concerns variation with spatial position.
  • One-dimensional model: Uses cross-sectional mean quantities that vary mainly along one coordinate.
  • Two- or three-dimensional model: Retains two or three spatial components and gradients where needed.
  • Incompressible versus compressible: Concerns whether density changes are negligible for the analysis.
  • Rotational versus irrotational: Concerns local fluid-element spin, measured by vorticity.
  • Laminar, transitional, turbulent: Concerns flow stability and the relative roles of viscous and inertial effects.

Uniformity does not by itself require constant conduit area. Geometry and continuity may imply an area relation only after the flow model, density behavior, and discharge constraints are stated.

Steady Flow Can Still Accelerate

Steady flow removes explicit time dependence at a fixed point, but particles can accelerate while moving through spatial velocity gradients. A steady nozzle flow is a standard example: local acceleration is zero while convective acceleration is nonzero.

Material derivative

The rate of change of a field quantity experienced by a moving fluid particle.

Material Derivative of a Scalar

Converts an Eulerian scalar field into the rate experienced along a particle trajectory.

DϕDt=∂ϕ∂t+u∂ϕ∂x+v∂ϕ∂y+w∂ϕ∂z\frac{D\phi}{Dt} = \frac{\partial \phi}{\partial t} +u\frac{\partial \phi}{\partial x} +v\frac{\partial \phi}{\partial y} +w\frac{\partial \phi}{\partial z}

Variables

SymbolDescriptionUnit
ϕ\phiScalar flow property such as pressure, density, or temperature-
D/DtD/DtMaterial or substantial derivative following a particle-

Fluid-Particle Acceleration

Separates local and convective acceleration of the velocity field.

a=DVDt=∂V∂t+(V⋅∇)V\mathbf{a} = \frac{D\mathbf{V}}{Dt} = \frac{\partial \mathbf{V}}{\partial t} + (\mathbf{V}\cdot\nabla)\mathbf{V}ax=∂u∂t+u∂u∂x+v∂u∂y+w∂u∂za_x= \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z}

Variables

SymbolDescriptionUnit
a\mathbf{a}Material acceleration vectorm/s2m/s^2
axa_xAcceleration component in the x directionm/s2m/s^2
∇\nablaSpatial gradient operator1/m
Local and convective accelerationParticle acceleration combines time and spatial velocity-field changes.local changespatial gradientacceleration

Local and convective acceleration

Particle acceleration combines time and spatial velocity-field changes.

Steady One-Dimensional Convective Acceleration

Applies when the mean speed varies mainly along a streamline coordinate.

as=VdVdsa_s=V\frac{dV}{ds}

Variables

SymbolDescriptionUnit
asa_sAcceleration along streamline coordinate sm/s2m/s^2
VVMean speedm/s
ssDistance along the streamlinem

Streamline

A curve everywhere tangent to the instantaneous local velocity vector. At that instant, the normal velocity component across the streamline is zero.

Pathline

The actual trajectory traced by one identified fluid particle through time.

Streakline

The instantaneous locus of all particles that previously passed through a specified fixed point, such as a continuous dye-injection location.

Timeline

A line formed by particles marked simultaneously; later distortion of the line reveals spatial velocity gradients and deformation.

Streamline Differential Equation

Defines the local direction of a streamline in a three-dimensional velocity field.

dxu=dyv=dzw\frac{dx}{u}=\frac{dy}{v}=\frac{dz}{w}

Variables

SymbolDescriptionUnit
dx,dy,dzdx,dy,dzDifferential coordinate changes along a streamlinem
u,v,wu,v,wLocal Cartesian velocity componentsm/s

When Streamlines, Pathlines, and Streaklines Coincide

For a sufficiently smooth steady velocity field, streamlines, pathlines, and streaklines coincide. In unsteady flow they generally differ, so the visualization method must be identified before interpreting a measured curve.

Streamline, pathline, and streaklineThree common flow-description lines and their geometric meaning.streamlinepathlinestreakline

Streamline, pathline, and streakline

Three common flow-description lines and their geometric meaning.

Interactive 3D Fluid Kinematics & Flow Visualizer

Explore 3D velocity vector fields, Lagrangian particle pathlines, instantaneous streamlines, and dye streaklines. Toggle between steady stagnation flow, Rankine vortices, and shear flows, and observe how fluid elements translate, rotate, and undergo linear and angular strain.

Fluid Kinematics 3D Vector Studio

Spatial velocity fields, 3D streamline integration, Lagrangian pathline tracing, vorticity vectors, and fluid element strain deformation.

Loading 3D Fluid Kinematics Simulation…

Kinematic Flow Controls

Flow Field Strength (kk)1.20
Particle Integration Speed1.0x
Divergence (Continuity) ∇⋅V\nabla \cdot \mathbf{V}∂u/∂x + ∂v/∂y = k - k = 0
Vorticity (Spin) ω=∇×V\mathbf{\omega} = \nabla \times \mathbf{V}ω_z = ∂v/∂x - ∂u/∂y = 0
Strain Rate Tensor ε˙\dot{\mathbf{\varepsilon}}ε_xx = k, ε_yy = -k
Flow Dynamics ClassStagnation flow with orthogonal streamlines and constant elongational strain.

Volumetric discharge

The volume flow rate crossing a specified surface, obtained by integrating the velocity component normal to that surface.

Volumetric Discharge

Integrates the normal velocity component over a surface.

Q=∫AV⋅n dAQ=\int_A \mathbf{V}\cdot\mathbf{n}\,dAQ=VavgAQ=V_{\text{avg}}A

Variables

SymbolDescriptionUnit
QQVolumetric dischargem3/sm^3/s
n\mathbf{n}Chosen unit normal to the surface-
VavgV_{\text{avg}}Area-average normal velocitym/s

Point Velocity Is Not Mean Velocity

A Pitot tube, acoustic probe, or numerical sample may provide a local velocity. Discharge calculations require an area integral or an appropriate section-average velocity; energy and momentum calculations may additionally require correction coefficients when the profile is strongly nonuniform.

Conservation of mass

Mass can neither be created nor destroyed within the continuum model; accumulation inside a control volume is balanced by net mass flux across its boundary.

Integral Continuity for a Fixed Control Volume

General mass balance for a control volume fixed in space.

ddt∫CVρ dV+∫CSρ V⋅n dA=0\frac{d}{dt}\int_{CV}\rho\,dV + \int_{CS}\rho\,\mathbf{V}\cdot\mathbf{n}\,dA =0

Variables

SymbolDescriptionUnit
CVCVControl volume fixed in space-
CSCSControl surface bounding the fixed control volume-
ρ\rhoFluid densitykg/m3kg/m^3
Continuity through a control volumeMass inflow, storage, and outflow are balanced.inflowstorageoutflow

Continuity through a control volume

Mass inflow, storage, and outflow are balanced.

Moving Control Volumes Require Relative Velocity

The fixed-control-volume continuity equation above uses the fluid velocity V\mathbf{V} through a stationary boundary. If the control surface itself moves or deforms, the mass-flux term must use the fluid velocity relative to that control surface. Do not label the fixed-boundary equation as a general moving-control-volume relation.

Differential Continuity Equation

Local conservation of mass for a continuum.

∂ρ∂t+∇⋅(ρV)=0\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{V})=0

Variables

SymbolDescriptionUnit
ρ\rhoFluid densitykg/m3kg/m^3
V\mathbf{V}Velocity fieldm/s

Constant-Density Incompressible Continuity

Velocity-field compatibility condition for constant density.

∇⋅V=∂u∂x+∂v∂y+∂w∂z=0\nabla\cdot\mathbf{V} = \frac{\partial u}{\partial x} +\frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} =0

Variables

SymbolDescriptionUnit
u,v,wu,v,wCartesian velocity componentsm/s

Steady One-Dimensional Mass Continuity

Relates mass flow at discrete sections of a streamtube.

m˙=ρAV=constant\dot m=\rho AV=\text{constant}A1V1=A2V2=Qfor constant-density flowA_1V_1=A_2V_2=Q \qquad \text{for constant-density flow}

Variables

SymbolDescriptionUnit
m˙\dot mMass flow ratekg/s
AACross-sectional area normal to mean flowm2m^2
VVSection-average velocitym/s

Junctions and Storage

At a junction with negligible storage, the algebraic sum of incoming and outgoing mass flows is zero. In a tank or reservoir with changing storage, inflow and outflow need not be equal instantaneously; their difference changes stored volume or mass.

For an incompressible liquid in a vertical tank of plan area AtA_t,

Atdhdt=Qin−QoutA_t\frac{dh}{dt}=Q_{\text{in}}-Q_{\text{out}}

provided the plan area is constant over the level range considered.

Normal strain rate

The instantaneous extension or compression rate of a fluid element along a coordinate direction.

Shear strain rate

The instantaneous angular deformation rate associated with cross-gradients of velocity components.

Cartesian Strain Rates

Selected normal and engineering shear strain rates for a velocity field.

ε˙x=∂u∂x,ε˙y=∂v∂y,ε˙z=∂w∂z\dot\varepsilon_x=\frac{\partial u}{\partial x}, \qquad \dot\varepsilon_y=\frac{\partial v}{\partial y}, \qquad \dot\varepsilon_z=\frac{\partial w}{\partial z}γ˙xy=∂u∂y+∂v∂x\dot\gamma_{xy} = \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}

Variables

SymbolDescriptionUnit
ε˙x,ε˙y,ε˙z\dot\varepsilon_x,\dot\varepsilon_y,\dot\varepsilon_zNormal strain rates1/s
γ˙xy\dot\gamma_{xy}Engineering shear strain rate in the xy plane1/s

Vorticity

The curl of the velocity field. The vorticity vector equals twice the local angular-velocity vector of an infinitesimal fluid element.

Vorticity and Fluid-Element Rotation

Relates velocity curl to local average angular rotation.

ω=∇×V\boldsymbol{\omega}=\nabla\times\mathbf{V}Ω=12ω\boldsymbol{\Omega}=\frac{1}{2}\boldsymbol{\omega}

Variables

SymbolDescriptionUnit
ω\boldsymbol{\omega}Vorticity vector1/s
Ω\boldsymbol{\Omega}Average angular-velocity vector of a fluid elementrad/s

Curved Motion Does Not Prove Rotational Flow

A particle can follow a curved path while the fluid element has zero average spin. Rotationality is determined from ∇×V\nabla\times\mathbf{V}, not from visible streamline curvature alone.

Circulation

The closed-line integral of the tangential velocity component around a specified contour.

Circulation

Measures integrated tangential motion around a closed curve.

Γ=∮CV⋅ds\Gamma=\oint_C\mathbf{V}\cdot d\mathbf{s}

Variables

SymbolDescriptionUnit
Γ\GammaCirculationm2/sm^2/s
CCClosed integration contour-

Circulation and Vorticity

For a sufficiently smooth field, Stokes' theorem relates circulation around a closed contour to the flux of vorticity through a spanning surface. This connection is useful for checking whether a visibly curved flow field is actually rotational.

Stream function

A scalar function used for two-dimensional incompressible flow so that the velocity components are generated from cross-derivatives and continuity is satisfied identically.

Two-Dimensional Stream Function

Defines velocity components for two-dimensional incompressible flow.

u=∂ψ∂y,v=−∂ψ∂xu=\frac{\partial\psi}{\partial y}, \qquad v=-\frac{\partial\psi}{\partial x}

Variables

SymbolDescriptionUnit
ψ\psiStream functionm2/sm^2/s
u,vu,vVelocity components in the x and y directionsm/s

Discharge Between Streamlines

Gives discharge per unit thickness between two streamlines in a two-dimensional incompressible flow.

q′=ψ2−ψ1q'=\psi_2-\psi_1

Variables

SymbolDescriptionUnit
q′q'Discharge per unit thickness normal to the flow planem2/sm^2/s
ψ1,ψ2\psi_1,\psi_2Stream-function values on the bounding streamlinesm2/sm^2/s

Stream-Function Scope

The Cartesian relations above are specifically for a two-dimensional incompressible flow representation. A stream function is not a universal scalar description for arbitrary three-dimensional flow.

Velocity potential

A scalar field whose gradient equals the velocity in an irrotational region.

Velocity Potential

Represents an irrotational velocity field locally as a gradient.

V=∇ϕ\mathbf{V}=\nabla\phi

Variables

SymbolDescriptionUnit
ϕ\phiVelocity potentialm2/sm^2/s
V\mathbf{V}Velocity fieldm/s

Local and Global Potential Requirements

Zero vorticity is the local compatibility condition for a smooth velocity potential. A single-valued potential throughout an entire region also depends on domain topology; simply connected regions avoid circulation around holes that can prevent a global single-valued potential.

Laplace Equation for Incompressible Potential Flow

Governs a velocity potential when the flow is both incompressible and irrotational.

∇2ϕ=0\nabla^2\phi=0

Variables

SymbolDescriptionUnit
∇2\nabla^2Laplacian operator1/m2whenactingonadimensionlessfield1/m^2 when acting on a dimensionless field
ϕ\phiVelocity potentialm2/sm^2/s

Two-Dimensional Potential-Flow Nets

For a regular two-dimensional incompressible irrotational flow, streamlines and equipotential lines intersect orthogonally. A useful graphical net has noncrossing streamlines and closely approximates curvilinear squares where the scale is appropriate.

The same mathematical structure appears in saturated, steady, two-dimensional seepage through homogeneous isotropic soil, although hydraulic head replaces velocity potential in the seepage formulation.

Orthogonal potential-flow netStreamlines and equipotentials form an orthogonal flow net.streamlinesequipotentialsflow net

Orthogonal potential-flow net

Streamlines and equipotentials form an orthogonal flow net.

Seepage Discharge from an Isotropic Flow Net

Estimates steady two-dimensional seepage per unit thickness through homogeneous isotropic soil.

q=kHNfNdq=kH\frac{N_f}{N_d}

Variables

SymbolDescriptionUnit
qqSeepage discharge per unit thicknessm2/sm^2/s
kkHydraulic conductivitym/s
HHTotal hydraulic-head differencem
NfN_fNumber of flow channels-
NdN_dNumber of equipotential drops-

Anisotropic Seepage Requires Transformation

The curvilinear-square construction applies directly to homogeneous isotropic media. In anisotropic soil, coordinates or conductivity must be transformed appropriately before the standard graphical relation is used.

Kinematics Analysis Workflow

  1. Define coordinates and distinguish local point velocities from section-average quantities.
  2. Classify the flow by time dependence, spatial variation, dimensionality, density behavior, and rotationality.
  3. Use the material derivative whenever the rate experienced by a moving particle is required.
  4. Apply mass conservation before energy or momentum relations.
  5. For a fixed control volume, use the actual velocity through the stationary control surface; for moving boundaries, use relative velocity.
  6. Check divergence for incompressibility and curl for rotationality when a velocity field is given.
  7. Use stream functions, potentials, and flow nets only when their dimensional, incompressibility, irrotationality, and boundary assumptions are satisfied.
Key Takeaways
  • Steady and uniform are different classifications: steady is temporal, while uniform is spatial.
  • The material derivative combines local and convective changes experienced by a moving fluid particle.
  • Fixed-control-volume continuity uses fluid velocity through a stationary boundary; moving control surfaces require relative velocity.
  • Constant-density incompressible flow requires a divergence-free velocity field.
  • Vorticity is the curl of velocity and equals twice the local fluid-element angular velocity.
  • Streamlines, pathlines, and streaklines coincide in a smooth steady flow but generally differ in unsteady flow.
  • A two-dimensional stream function enforces incompressible continuity; a velocity potential requires irrotational flow.
  • Potential-flow and seepage flow nets are powerful only when their boundary, dimensionality, and isotropy assumptions are respected.