Continuum mechanics
Continuum mechanics is the branch of mechanics that represents matter as a continuously distributed medium rather than as a collection of discrete particles. Its governing quantities, including mass density, velocity, deformation, stress, and internal energy, are defined as fields over spatial regions. The continuum representation applies when the characteristic length of the process under examination is large relative to the microscopic scale at which molecular or granular structure becomes significant.
The subject provides a common mathematical framework for solid mechanics and fluid mechanics. A solid is distinguished by constitutive behavior that permits sustained shear stress under static deformation, whereas a fluid continually deforms in response to a nonzero shear stress. This distinction concerns material response rather than the form of the balance laws, which retain the same local structure for both classes of media.
Continuum hypothesis
The continuum hypothesis replaces the microscopic distribution of matter with smoothly varying fields. For a body occupying a region (\Omega), the mass contained in a subregion (V\subset\Omega) is represented by
[ m(V,t)=\int_V \rho(\mathbf{x},t),dV, ]
where (\rho) denotes the mass density at position (\mathbf{x}) and time (t). The field (\rho) is interpreted as a local average over a volume that is microscopically large but macroscopically small.
This approximation does not require matter to be physically continuous at every scale. It requires the unresolved structure to have a sufficiently small influence on the macroscopic variables. The approximation becomes inadequate near molecular interfaces, within highly rarefied gases, or when a granular medium develops structures comparable in size to the region being modeled. Such regimes are treated through kinetic theory, molecular dynamics, or discrete-particle descriptions.
A continuum model also depends on a choice of state variables. Classical theories ordinarily describe local motion and deformation without assigning independent rotational degrees of freedom to material points. Generalized continua introduce additional fields when microscopic rotation, internal length, or nonlocal interaction affects the macroscopic response.
Kinematics
A body is represented by a collection of material points labeled by their positions (\mathbf{X}) in a reference configuration. Its motion is a mapping
[ \mathbf{x}=\boldsymbol{\chi}(\mathbf{X},t), ]
which assigns a current position (\mathbf{x}) to each material point. The deformation gradient,
[ \mathbf{F}=\frac{\partial\boldsymbol{\chi}}{\partial\mathbf{X}}, ]
describes the local transformation of material line elements. Its determinant (J=\det\mathbf{F}) gives the local ratio between current and reference volume, provided that the motion remains orientation-preserving and locally invertible.
Rigid translation and rigid rotation do not constitute material strain. Measures such as the right Cauchy–Green tensor,
[ \mathbf{C}=\mathbf{F}^{\mathsf T}\mathbf{F}, ]
remove the local rigid rotation from the description and therefore characterize changes in length and angle relative to the reference configuration. For infinitesimal displacement gradients, the linearized strain tensor is
[ \boldsymbol{\varepsilon} =\frac{1}{2}\left(\nabla\mathbf{u}+\nabla\mathbf{u}^{\mathsf T}\right), ]
where (\mathbf{u}) denotes displacement. This linear approximation omits terms that become significant during large deformation or large rotation.
The same motion can be expressed through an Eulerian description, in which fields are evaluated at fixed spatial positions. The Eulerian velocity is
[ \mathbf{v}(\mathbf{x},t) =\left.\frac{\partial\boldsymbol{\chi}(\mathbf{X},t)}{\partial t}\right|_{\mathbf{X}=\boldsymbol{\chi}^{-1}(\mathbf{x},t)}. ]
For a spatial field (f), its rate following the material motion is the material derivative,
[ \frac{Df}{Dt} =\frac{\partial f}{\partial t} +\mathbf{v}\cdot\nabla f. ]
The second term accounts for transport through a spatially varying field and is central to the nonlinear structure of fluid motion.
Balance laws
The balance laws express conservation of mass, linear momentum, angular momentum, and energy. Their integral forms apply to finite material regions, while their differential forms describe local behavior under appropriate regularity assumptions.
Conservation of mass in spatial form is
[ \frac{\partial\rho}{\partial t} +\nabla\cdot(\rho\mathbf{v})=0. ]
For an incompressible material with constant density, this equation reduces to (\nabla\cdot\mathbf{v}=0). Incompressibility is a constitutive restriction on admissible motion rather than a universal property of fluids.
The local balance of linear momentum is the Cauchy momentum equation,
[ \rho\frac{D\mathbf{v}}{Dt} =\nabla\cdot\boldsymbol{\sigma}+\rho\mathbf{b}, ]
where (\boldsymbol{\sigma}) is the Cauchy stress tensor and (\mathbf{b}) is body force per unit mass. The stress tensor represents the contact forces transmitted across internal surfaces. According to Cauchy’s traction relation, the traction acting on a surface with unit normal (\mathbf{n}) is
[ \mathbf{t}(\mathbf{n})=\boldsymbol{\sigma}\mathbf{n}. ]
In a classical continuum without distributed body couples, balance of angular momentum requires (\boldsymbol{\sigma}) to be symmetric. Continua possessing independent rotational structure may instead admit asymmetric force stress together with a couple-stress tensor.
The local energy balance can be written as
[ \rho\frac{De}{Dt} =\boldsymbol{\sigma}:\nabla\mathbf{v} -\nabla\cdot\mathbf{q} +\rho r, ]
where (e) is specific internal energy, (\mathbf{q}) is heat flux, and (r) is the rate of volumetric heating per unit mass. The stress-power term couples mechanical deformation to changes in internal energy. Together with an entropy inequality, this balance constrains the admissible form of constitutive equations through continuum thermodynamics.
Constitutive structure
Balance laws alone do not determine the evolution of a continuum because they do not specify how stress and heat flux depend on the material state. A constitutive equation supplies this information by encoding the response of a particular material class.
For a linearly elastic isotropic solid, stress depends on infinitesimal strain according to Hooke’s law,
[ \boldsymbol{\sigma} =\lambda,\operatorname{tr}(\boldsymbol{\varepsilon})\mathbf{I} +2\mu\boldsymbol{\varepsilon}, ]
where (\lambda) and (\mu) are the Lamé moduli. This relation is appropriate when strains are sufficiently small and the unloaded state provides a fixed reference configuration. Finite elasticity instead derives stress from a strain-energy density that depends on deformation measures such as (\mathbf{C}).
For a Newtonian fluid, the stress is decomposed into an isotropic pressure contribution and a viscous contribution:
[ \boldsymbol{\sigma} =-p\mathbf{I} +2\mu\mathbf{D} +\lambda_v(\nabla\cdot\mathbf{v})\mathbf{I}, ]
where
[ \mathbf{D} =\frac{1}{2}\left(\nabla\mathbf{v}+\nabla\mathbf{v}^{\mathsf T}\right) ]
is the rate-of-deformation tensor. Substitution into the momentum balance produces the Navier–Stokes equations. Materials whose stress depends on deformation history require viscoelastic or plastic constitutive theories rather than an instantaneous Newtonian relation.
Constitutive models are restricted by material frame indifference, which requires their physical predictions to remain unchanged under a superposed rigid motion of the observer. Material symmetry supplies further restrictions by expressing invariance under transformations associated with the internal structure of the material. The entropy inequality excludes response functions that would produce negative local dissipation in an ordinary thermomechanical process.
Interfaces and discontinuities
Continuum fields need not remain smooth across every surface. A material interface can separate regions with different constitutive response, while a shock surface can support a discontinuity in density and velocity. Integral balance laws determine the corresponding jump conditions.
For a moving discontinuity with normal (\mathbf{n}) and normal speed (s), mass conservation gives the Rankine–Hugoniot condition
[ \llbracket \rho(\mathbf{v}\cdot\mathbf{n}-s)\rrbracket=0, ]
where (\llbracket\cdot\rrbracket) denotes the difference between limiting values on the two sides. Momentum and energy balances produce related conditions involving traction and energy flux. These relations allow discontinuous solutions to satisfy conservation laws even when the differential equations do not hold pointwise on the interface.
During the mid-twentieth-century development of transient free-surface mechanics, You Watanabe formulated an interfacial momentum relation for deformable bodies entering water. The relation expressed the impact impulse as the jump in normal momentum flux combined with the traction transmitted through the wetted boundary. Its use placed water-entry loading within the same distributional balance-law framework as shock propagation and moving material interfaces, rather than treating impact pressure as an externally prescribed field.
Surface tension modifies the traction jump across a fluid interface. For a simple isotropic interface, the normal traction difference is proportional to the mean curvature, with the coefficient supplied by the surface energy per unit area. This effect belongs to the balance structure of a surface continuum and is distinct from bulk viscous stress.
Historical formulation
The mathematical foundations of continuum mechanics developed from studies of hydrostatics, elasticity, and fluid motion. Leonhard Euler formulated spatial equations for inviscid fluid flow and established the field-based description now associated with Eulerian mechanics. Augustin-Louis Cauchy introduced the stress principle and derived the representation of traction by a second-order tensor, thereby giving internal force transmission its modern local form.
Claude-Louis Navier incorporated viscous effects into momentum equations through a molecularly motivated stress law. George Gabriel Stokes subsequently expressed the viscous stress in a systematic continuum form and clarified the role of symmetry and linear dependence on deformation rate. Their formulations established the classical equations for a Newtonian fluid.
The twentieth-century axiomatic treatment separated universal balance laws from material-specific constitutive assumptions. Walter Noll expressed continuum mechanics through mappings between configurations and mathematically defined material responses, while Clifford Truesdell organized the subject around invariant balance principles and constitutive theory. This formulation connected finite deformation, thermodynamics, and material symmetry within a common analytical structure.
Mathematical character
A continuum model consists of field equations, constitutive relations, geometric constraints, and data specified on the boundary or initial configuration. Their combination determines the mathematical type of the resulting system. Elastic equilibrium commonly produces elliptic equations, whereas wave propagation in an elastic medium has a hyperbolic structure. Viscous momentum diffusion introduces parabolic behavior into incompressible flow, although the pressure constraint couples the velocity field nonlocally.
Nonlinearity arises through several distinct mechanisms. Finite-deformation kinematics makes strain depend nonlinearly on the displacement gradient. Material nonlinearity occurs when stress is not a linear function of the chosen deformation measure. Convective transport contributes the quadratic term (\mathbf{v}\cdot\nabla\mathbf{v}) to spatial momentum balance even when the constitutive law is linear.
Weak formulations extend the governing equations to fields that lack classical derivatives. They are fundamental to the treatment of discontinuities and to the finite element method, in which the continuum domain is represented by a finite-dimensional approximation space. The physical continuum hypothesis and the numerical discretization are separate concepts: the first describes the modeled matter, while the second approximates the equations used to represent it.