Fluid Kinematics

Learning Objectives

  • Classify flows as steady or unsteady, uniform or nonuniform, one-, two-, or three-dimensional, and laminar or turbulent.
  • Distinguish Lagrangian particle tracking from the Eulerian field description used in most engineering analyses.
  • Calculate local, convective, and total material acceleration from a velocity field.
  • Apply integral and differential conservation of mass to incompressible and compressible flow.
  • Interpret translation, deformation, angular rotation, circulation, and vorticity of fluid elements.
  • Construct and use streamlines, pathlines, streaklines, stream functions, velocity potentials, and flow nets.

Fluid kinematics describes how fluid moves without first asking what forces produce that motion. The subject supplies the velocity, acceleration, deformation, and continuity relationships later used in energy, momentum, pipe-flow, seepage, and open-channel analyses.

Eulerian and Lagrangian Descriptions

  • Lagrangian description: Follow a particular fluid particle and record its position, velocity, and acceleration as functions of time.
  • Eulerian description: Specify flow properties at fixed locations in space, such as V(x,y,z,t)\mathbf{V}(x,y,z,t), and observe particles as they pass through the field.

Laboratory instruments, numerical models, and most civil-engineering calculations use the Eulerian description, while particle tracking is useful for pollutant transport, sediment trajectories, and visualization.

Eulerian Velocity Field

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

V(x,y,z,t)=ui+vj+wk\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

Common Flow Classifications

  • Steady: At a fixed point, flow properties do not change with time, so ()/t=0\partial()/\partial t=0.
  • Unsteady: At least one property at a fixed point changes with time.
  • Uniform: A selected property does not change along the direction under consideration at an instant.
  • Nonuniform: Spatial gradients are present.
  • One-dimensional model: Cross-sectional variation is represented by average quantities that change mainly along one coordinate.
  • Two- or three-dimensional flow: Two or three spatial velocity components and gradients are required.
  • Laminar, transitional, turbulent: Classified by the relative importance of viscous and inertial effects and by flow instability.
  • Rotational or irrotational: Classified by whether fluid elements have nonzero average angular rotation.

Steady Does Not Mean Zero Acceleration

A steady flow can have substantial acceleration when velocity changes with position. Water moving through a contraction is steady if the field is time-invariant, yet each particle accelerates as it enters the smaller area. This is convective acceleration.

Material Derivative

The rate of change experienced by a moving fluid particle combines the explicit time change at a fixed location with the change caused by the particle moving through spatial gradients.

Material Derivative of a Scalar

Converts an Eulerian field derivative into the rate experienced by a moving particle.

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
ϕ\phiAny scalar flow property, such as pressure, density, or temperature-
D/DtD/DtMaterial or substantial derivative following a particle1/s when applied to a dimensionless scalar

Fluid-Particle Acceleration

Separates local and convective acceleration of the velocity field.

a=DVDt=Vtlocal+(V)Vconvective\mathbf{a} = \frac{D\mathbf{V}}{Dt} = \underbrace{\frac{\partial \mathbf{V}}{\partial t}}_{\text{local}} + \underbrace{(\mathbf{V}\cdot\nabla)\mathbf{V}}_{\text{convective}}

For the x component,

ax=ut+uux+vuy+wuza_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
axa_xParticle acceleration in the x directionm/s2m/s^2
\nablaSpatial gradient operator1/m

One-Dimensional Convective Acceleration

For steady flow whose mean velocity varies mainly along a streamline coordinate ss,

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

This relation explains acceleration through nozzles, contractions, bridge openings, and varying open-channel sections.

Streamline

A line tangent everywhere to the instantaneous velocity vector. No fluid crosses a streamline because its normal velocity component is zero.

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}

Pathline

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

Streakline

The locus at an instant 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. Its later distortion reveals deformation and velocity gradients.

When the Flow Lines Coincide

In steady flow, streamlines, pathlines, and streaklines coincide. In unsteady flow they can have very different shapes, so the term used must match the measurement or visualization method.

Discharge and Mean Velocity

The volumetric flow rate through a surface is the integral of the velocity component normal to that surface. A single average velocity is defined so that it carries the same discharge through the area.

Volumetric Discharge

Integrates normal velocity over a cross-sectional area.

Q=AVndA=VavgAQ=\int_A \mathbf{V}\cdot\mathbf{n}\,dA =V_{avg}A

Variables

SymbolDescriptionUnit
QQVolumetric flow ratem3/sm^3/s
n\mathbf{n}Unit normal vector directed outward from the surface-
VavgV_{avg}Area-average normal velocitym/s

Point Velocity Is Not Automatically Mean Velocity

Pitot tubes, acoustic probes, and numerical output may provide local velocity. Energy, momentum, and discharge calculations usually require cross-sectional averages and may also require kinetic-energy or momentum correction coefficients when the profile is nonuniform.

Conservation of Mass

For a control volume, accumulation of mass inside plus net mass outflow through its boundary must equal zero.

Integral Continuity Equation

General conservation of mass for a fixed or moving control volume.

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

Variables

SymbolDescriptionUnit
CVCVControl volume-
CSCSControl surface bounding the control volume-
ρ\rhoFluid densitykg/m3kg/m^3

Differential Continuity Equation

Local mass-conservation equation for a continuum.

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

For constant-density incompressible flow,

V=ux+vy+wz=0\nabla\cdot\mathbf{V} = \frac{\partial u}{\partial x} +\frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} =0

Steady One-Dimensional Continuity

Relates mass flow between sections of a streamtube.

m˙=ρAV=constant\dot m=\rho AV=\text{constant}

For incompressible flow,

A1V1=A2V2=QA_1V_1=A_2V_2=Q

Variables

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

Translation, Deformation, and Rotation

A small fluid element can translate, change volume, change shape, and rotate. Velocity gradients control these motions.

  • Normal strain rates describe extension or compression along coordinate directions.
  • Shear strain rates describe angular distortion.
  • Average angular velocity is one half of the vorticity vector.

Normal and Shear Strain Rates

Describes instantaneous deformation of a fluid element in Cartesian coordinates.

ε˙x=ux,ε˙y=vy,ε˙z=wz\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=uy+vx\dot\gamma_{xy}=\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}

Vorticity

The curl of the velocity field. It measures twice the local angular velocity of an infinitesimal fluid element.

Vorticity and Fluid-Element Rotation

Relates the curl of velocity to local 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}Angular velocity of a fluid elementrad/s

Curved Streamlines Do Not Necessarily Mean Rotational Flow

A particle can travel along a curved path while the fluid element itself has zero average spin. Rotationality is determined by ×V\nabla\times\mathbf{V}, not simply by the visible curvature of streamlines.

Circulation

The line integral of tangential velocity around a closed curve.

Circulation

Measures the integrated rotational tendency around a closed contour.

Γ=CVds\Gamma=\oint_C \mathbf{V}\cdot d\mathbf{s}

By Stokes' theorem,

Γ=A(×V)ndA\Gamma=\int_A (\nabla\times\mathbf{V})\cdot\mathbf{n}\,dA

Stream Function for Two-Dimensional Incompressible Flow

A stream function ψ(x,y)\psi(x,y) can be defined so that continuity is satisfied automatically:

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

Curves of constant ψ\psi are streamlines, and the discharge per unit depth between two streamlines is Δψ\Delta\psi.

Stream-Function Discharge

Gives two-dimensional discharge per unit thickness between streamlines.

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

Variables

SymbolDescriptionUnit
qq'Discharge per unit thickness normal to the flow planem2/sm^2/s

Velocity Potential and Irrotational Flow

A velocity potential ϕ\phi exists in a simply connected irrotational region when

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

so that u=ϕ/xu=\partial\phi/\partial x, v=ϕ/yv=\partial\phi/\partial y, and w=ϕ/zw=\partial\phi/\partial z. Combining irrotationality with incompressibility gives Laplace's equation.

Laplace Equation for Potential Flow

Governs incompressible irrotational potential flow.

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

For a two-dimensional stream function in an irrotational region,

2ψ=0\nabla^2\psi=0

Flow Nets

A flow net is an orthogonal network of constant-potential lines and streamlines for two-dimensional, incompressible, irrotational flow. In geotechnical and hydraulic engineering, analogous flow nets represent seepage through isotropic soil.

A properly drawn flow net has approximately curvilinear square fields, no crossing streamlines, and boundaries that satisfy the physical no-flow or constant-head conditions.

Seepage Discharge from a Flow Net

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

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

Variables

SymbolDescriptionUnit
qqSeepage discharge per unit thicknessm2/sm^2/s
kkHydraulic conductivitym/s
HHTotal head loss across the flow regionm
NfN_fNumber of flow channels-
NdN_dNumber of equipotential drops-

Anisotropic Seepage

For anisotropic soil, physical coordinates must be transformed or an equivalent conductivity treatment must be used before applying the curvilinear-square rule. A flow net drawn directly in untransformed anisotropic coordinates is generally not orthogonal.

Kinematics Analysis Workflow

  1. State the coordinate system and write the velocity field or section-average velocities.
  2. Classify time dependence, dimensionality, and compressibility assumptions.
  3. Use the material derivative when particle acceleration or property change is required.
  4. Apply continuity before energy or momentum equations.
  5. Distinguish local point values from cross-sectional averages.
  6. Check velocity-field compatibility using divergence and, where relevant, curl.
  7. Use stream functions, potentials, or flow nets only when their incompressibility and irrotationality assumptions are satisfied.
Key Takeaways
  • The Eulerian field description becomes a particle rate through the material derivative.
  • Steady flow can accelerate through spatial velocity gradients.
  • Conservation of mass is expressed generally by the integral continuity equation and locally by ρ/t+(ρV)=0\partial\rho/\partial t+\nabla\cdot(\rho\mathbf{V})=0.
  • For incompressible flow, the velocity field is divergence-free.
  • Vorticity is the curl of velocity and equals twice the local fluid-element angular velocity.
  • A stream function satisfies two-dimensional incompressible continuity automatically; a velocity potential requires irrotational flow.
  • Flow nets are powerful only when boundary conditions, isotropy, and geometric construction rules are respected.