Continuum Mechanics and Field Theory

Intuition [Beginner]

So far we have treated objects as point particles or rigid bodies: collections of discrete masses with fixed relationships. The real world is messier. A block of rubber deforms under load. Water flows around an obstacle. Sound waves propagate through air. In every case, matter is not a collection of isolated points but a continuum -- a medium where every small volume has a density, a velocity, and internal forces.

The continuum hypothesis is the assumption that matter is continuously distributed. We do not track individual atoms; instead we define fields -- functions of position and time -- that describe bulk properties: mass density $\rho (\mathbf{x},t)$, velocity $\mathbf{v}(\mathbf{x},t)$, pressure $p(\mathbf{x},t)$. This is an approximation. It breaks down when the length scale of interest approaches the mean free path between molecules (rarefied gases, shock fronts). But for the vast majority of engineering and physics problems, it works spectacularly well.

Two perspectives exist for describing motion in a continuum. The Lagrangian description tracks individual material elements: you label a piece of fluid by its initial position $\mathbf{X}$ and follow it. The Eulerian description fixes a point in space $\mathbf{x}$ and watches what flows through it. Imagine watching a river. Lagrangian: you ride a leaf downstream. Eulerian: you stand on a bridge and measure the current at a fixed point. Both carry the same physics; they differ in bookkeeping.

The internal forces in a continuum are captured by the stress tensor ${\sigma }_{ij}$. When you squeeze a material, the force on a surface element depends not just on the magnitude of the squeeze but on the orientation of the surface. A single number (pressure) is not enough; you need a $3\times 3$ tensor. For a fluid at rest, stress is isotropic (pure pressure: ${\sigma }_{ij}=-p{\delta }_{ij}$). For a deforming solid or viscous fluid, off-diagonal components appear -- these are shear stresses.

This unit bridges particle mechanics to field theory. Once you can describe a continuum by fields obeying partial differential equations, you are ready for electromagnetism (the EM field is a continuum too, just not made of matter) and general relativity (spacetime itself is a continuum with geometry governed by the stress-energy tensor).

Visual [Beginner]

LAGRANGIAN vs EULERIAN DESCRIPTION
===================================

Lagrangian: follow the material          Eulerian: watch the point
                                          
  t=0        t=1        t=2              t=0   t=1   t=2
                                         
  A--B       A--B        A--B            x1: T  S    R
  |  |  -->  |  /   -->   \  |           x2: U  T    S
  C--D       C--D          C--D          x3: V  U    T
                                         
 (material labels          (fixed spatial
  deform with body)         grid, fluid
                            flows through)

STRESS TENSOR COMPONENTS
=========================

         sigma_xx (normal)     
            |                  
            v                  
   +--------+--------+        
   |                  |  --> sigma_xy (shear)
   |   material       |        
   |   element        |  --> sigma_xz (shear)
   |                  |        
   +--------+--------+        
                             

   Full tensor:  [ sigma_xx  sigma_xy  sigma_xz ]
                 [ sigma_yx  sigma_yy  sigma_yz ]
                 [ sigma_zx  sigma_zy  sigma_zz ]

   For fluid at rest: sigma_ij = -p * delta_ij
     (only diagonal, all equal, negative = compressive)

Concept	Particle mechanics	Continuum mechanics
Degrees of freedom	Finite ($3N$)	Infinite (field values)
Configuration	Positions ${\mathbf{r}}_i$	Fields $\varphi (\mathbf{x},t)$
Equations of motion	ODEs (Newton, Lagrange)	PDEs (Euler, Navier-Stokes)
Internal forces	Discrete contact/constraint	Stress tensor ${\sigma }_{ij}$
Mass	Point mass $m_i$	Density $\rho (\mathbf{x},t)$

Worked example [Beginner]

Problem: Incompressible fluid of density $\rho$ flows steadily through a pipe that narrows from radius $R_1$ to radius $R_2$. If the entrance velocity is $v_1$, find the exit velocity $v_2$ and the pressure drop.

Solution:

Step 1: Conservation of mass (continuity equation). For incompressible flow, the volume flux is constant:

$$A_1{v}_1=A_2{v}_2$$

where $A_i=\pi R_i^2$.

Step 2: Solve for $v_2$:

$$v_2=v_1\frac{A_1}{A_2}=v_1\frac{R_1^2}{R_2^2}$$

Step 3: Apply Bernoulli's equation (steady, incompressible, inviscid, along a streamline):

$$p_1+\frac{1}{2}\rho v_1^2=p_2+\frac{1}{2}\rho v_2^2$$

Step 4: Solve for pressure drop:

$$p_1-p_2=\frac{1}{2}\rho (v_2^2-v_1^2)=\frac{1}{2}\rho v_1^2\left(\frac{R_1^4}{R_2^4}-1\right)$$

Since $R_1>R_2$, the exit velocity is higher and the pressure is lower. This is the Bernoulli effect: where the pipe narrows, the fluid speeds up and the pressure drops. This principle drives carburetors, aspirators, and the lift on an airplane wing.

Check your understanding [Beginner]

Formal definition [Intermediate+]

Continuum hypothesis. Let $\Delta V$ be a volume element with characteristic length $\ell$. If $\ell$ is much larger than the molecular mean free path $\lambda$ (Knudsen number $\mathrm{Kn}=\lambda /\ell \ll 1$), then matter in $\Delta V$ can be treated as continuously distributed. We define macroscopic fields as averages over $\Delta V$:

$$\rho (\mathbf{x},t)=\frac{\Delta m}{\Delta V},\mathbf{v}(\mathbf{x},t)=\frac{1}{\Delta m}\sum _{\alpha \in \Delta V}{m}_{\alpha }{\mathbf{v}}_{\alpha }$$

Kinematics. The motion of a continuum is described by the deformation map $\mathbf{\chi }(\mathbf{X},t)$ which maps material coordinates $\mathbf{X}$ (labeling a material element in its reference configuration) to spatial coordinates $\mathbf{x}$:

$$\mathbf{x}=\mathbf{\chi }(\mathbf{X},t)$$

The deformation gradient is:

$$F_{iJ}=\frac{\partial {\chi }_i}{\partial X_J}$$

Strain. The Green-Lagrange strain tensor measures deformation from the reference configuration:

$$E_{IJ}=\frac{1}{2}\left(F_{kI}{F}_{kJ}-{\delta }_{IJ}\right)$$

For small deformations this reduces to the Cauchy infinitesimal strain:

$${ϵ}_{ij}=\frac{1}{2}\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right)$$

where $\mathbf{u}=\mathbf{x}-\mathbf{X}$ is the displacement.

Stress. The Cauchy stress tensor ${\sigma }_{ij}$ relates the traction vector $\mathbf{t}$ on a surface with normal $\mathbf{n}$:

$$t_i={\sigma }_{ij}{n}_j$$

Conservation laws. The fundamental equations of continuum mechanics derive from conservation of mass, momentum, and energy applied to an arbitrary material volume.

Mass conservation (continuity equation):

$$\frac{\partial \rho }{\partial t}+\frac{\partial (\rho v_j)}{\partial x_j}=0$$

Momentum conservation (Cauchy momentum equation):

$$\rho \frac{D{v}_i}{Dt}=\frac{\partial {\sigma }_{ji}}{\partial x_j}+\rho b_i$$

where $D/Dt=\partial /\partial t+v_j\partial /\partial x_j$ is the material derivative and $b_i$ is the body force per unit mass.

Euler equations (inviscid fluid). For an ideal (inviscid, incompressible) fluid, ${\sigma }_{ij}=-p{\delta }_{ij}$:

$$\rho \frac{D{v}_i}{Dt}=-\frac{\partial p}{\partial x_i}+\rho b_i$$

Navier-Stokes equations (Newtonian viscous fluid). With a linear constitutive relation:

$${\sigma }_{ij}=-p{\delta }_{ij}+\mu \left(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\frac{\partial v_k}{\partial x_k}\right)+\zeta {\delta }_{ij}\frac{\partial v_k}{\partial x_k}$$

the momentum equation becomes:

$$\rho \frac{D{v}_i}{Dt}=-\frac{\partial p}{\partial x_i}+\mu \frac{{\partial }^2{v}_i}{\partial x_j\partial x_j}+\left(\zeta +\frac{\mu }{3}\right)\frac{\partial }{\partial x_i}\left(\frac{\partial v_k}{\partial x_k}\right)+\rho b_i$$

For incompressible flow ($\nabla \cdot \mathbf{v}=0$):

$$\rho \frac{D\mathbf{v}}{Dt}=-\nabla p+\mu {\nabla }^2\mathbf{v}+\rho \mathbf{b}$$

Elastic solid (Hookean). For a linear isotropic elastic solid, the constitutive law is:

$${\sigma }_{ij}=\lambda {ϵ}_{kk}{\delta }_{ij}+2\mu {ϵ}_{ij}$$

where $\lambda$ and $\mu$ are the Lame parameters. Combined with the momentum equation and the strain-displacement relation, this yields the Navier equation of elasticity:

$$\rho \frac{{\partial }^2{u}_i}{\partial t^2}=(\lambda +\mu )\frac{{\partial }^2{u}_k}{\partial x_i\partial x_k}+\mu \frac{{\partial }^2{u}_i}{\partial x_j\partial x_j}+\rho b_i$$

Lagrangian field theory. The transition from discrete to continuous systems replaces the Lagrangian $L={\sum }_i{T}_i-V_i$ with a Lagrangian density $\mathcal{L}(\varphi ,{\partial }_{\mu }\varphi ,x^{\mu })$:

$$S=\int \mathcal{L}{d}^{4}x$$

The Euler-Lagrange equations for fields are:

$$\frac{\partial \mathcal{L}}{\partial \varphi }-\frac{\partial }{\partial x^{\mu }}\left(\frac{\partial \mathcal{L}}{\partial ({\partial }_{\mu }\varphi )}\right)=0$$

This formalism applies to elasticity, fluid dynamics, electromagnetism, and all fundamental field theories.

Key results [Intermediate+]

1. Reynolds transport theorem. For any intensive property $\Psi$ per unit mass in a material volume $V(t)$:

$$\frac{d}{dt}{\int }_{V(t)}\rho \Psi dV={\int }_V\rho \frac{D\Psi }{Dt}dV$$

This is the bridge between Lagrangian and Eulerian formulations.

2. Helmholtz decomposition. Any sufficiently smooth vector field can be decomposed into an irrotational (curl-free) part and a solenoidal (divergence-free) part:

$$\mathbf{v}=\nabla \varphi +\nabla \times \mathbf{A}$$

This underlies potential flow theory and the decomposition of elastic waves into longitudinal (P) and transverse (S) waves.

Advanced treatment [Master]

Rational continuum mechanics. The modern, axiomatically rigorous treatment of continuum mechanics (Truesdell, Noll) begins with balance laws (conservation of mass, balance of linear and angular momentum, balance of energy, entropy inequality) imposed as axioms on a body manifold $\mathcal{B}$ mapped into physical space $\mathcal{S}$. Constitutive equations are restricted by the principle of material frame indifference (objectivity) and the entropy inequality (Clausius-Duhem).

The general constitutive equation for a simple material is:

$$\mathbf{\sigma }(t)=\hat{\mathbf{\sigma }}\left[{\mathbf{F}}^t(s),s\ge 0\right]$$

where ${\mathbf{F}}^t(s)=\mathbf{F}(t-s)$ is the history of the deformation gradient. Special cases include elastic (depends only on current $\mathbf{F}$), viscoelastic (depends on history), and hyperelastic (derivable from a stored energy function $W(\mathbf{F})$).

Non-Newtonian fluids. The Navier-Stokes constitutive relation is linear in the strain rate. Many real fluids (polymer solutions, blood, mud, glacier ice) are non-Newtonian: the effective viscosity depends on the strain rate. The power-law model ${\sigma }_{ij}^{\text{(visc)}}\propto ({\dot{ϵ}}_{kl}{\dot{ϵ}}_{kl}{)}^{(n-1)/2}{\dot{ϵ}}_{ij}$ is one parametrization. Rheology is the systematic study of these constitutive relations.

Existence and regularity of Navier-Stokes solutions. Whether the 3D incompressible Navier-Stokes equations always possess smooth solutions for smooth initial data is one of the Clay Millennium Prize problems. In 2D, existence and uniqueness are settled. In 3D, weak (Leray) solutions exist, but whether they are smooth (regular) remains open. The difficulty is the nonlinear convective term $(\mathbf{v}\cdot \nabla )\mathbf{v}$, which can transfer energy to ever-smaller scales (turbulent cascade) faster than viscosity can dissipate it.

Elastic waves. The Navier equation supports two families of plane waves in an infinite isotropic solid:

• Longitudinal (P-wave): $c_P=\sqrt{(\lambda +2\mu )/\rho }$
• Transverse (S-wave): $c_S=\sqrt{\mu /\rho }$

The ratio $c_P/c_S=\sqrt{(\lambda +2\mu )/\mu }>1$. Seismology exploits this: P-waves arrive before S-waves at a seismometer, and the time difference constrains the distance to the hypocenter.

From continuum mechanics to field theory. The Lagrangian density formalism of continuum mechanics is the direct precursor to classical field theory in particle physics. The displacement field $\mathbf{u}(\mathbf{x},t)$ of elasticity, the velocity potential $\varphi (\mathbf{x},t)$ of irrotational flow, and the scalar field $\varphi (x^{\mu })$ of a relativistic theory all obey the same Euler-Lagrange structure. The key difference is the metric: non-relativistic theories use Galilean-invariant Lagrangians ($t$ and $\mathbf{x}$ are distinct), while relativistic theories use Lorentz-invariant ones ($x^{\mu }$ is a 4-vector on Minkowski spacetime). General relativity promotes this further: the metric itself becomes a dynamical field.

Connection to electromagnetism. The EM field is a continuum field with its own stress tensor (the Maxwell stress tensor):

$$T_{ij}^{\text{EM}}={ϵ}_0\left(E_i{E}_j-\frac{1}{2}{\delta }_{ij}{E}_k{E}_k\right)+\frac{1}{{\mu }_0}\left(B_i{B}_j-\frac{1}{2}{\delta }_{ij}{B}_k{B}_k\right)$$

This enters the momentum balance of a charged continuum as a body force via the Lorentz force density $\mathbf{f}={\rho }_e\mathbf{E}+\mathbf{J}\times \mathbf{B}$.

Connections [Master]

• 09.04.01 -- Lagrangian mechanics generalizes directly to the Lagrangian density for fields.
• 09.06.01 -- Rigid body dynamics is the special case where deformation is forbidden; continuum mechanics permits it.
• 10.01.01 -- Electromagnetism is the first fundamental field theory; Maxwell's equations are PDEs for the $\mathbf{E}$ and $\mathbf{B}$ fields, analogous to the Navier-Stokes equations for $\mathbf{v}$.
• 13.01.01 -- General relativity treats the metric $g_{\mu \nu }$ as a continuum field governed by the Einstein field equations, sourced by the stress-energy tensor $T_{\mu \nu }$.
• 13.05.01 -- The stress-energy tensor in GR is the relativistic generalization of the Cauchy stress tensor plus energy-momentum density.
• 12.12.01 -- Quantum field theory applies the canonical quantization procedure to fields, beginning with the same Lagrangian density formalism.

Bibliography [Master]

1. Landau, L. D. and Lifshitz, E. M. Fluid Mechanics, 2nd ed., Butterworth-Heinemann, 1987. The standard reference for theoretical fluid mechanics, Chapters 1-3 on the Euler and Navier-Stokes equations.

2. Landau, L. D. and Lifshitz, E. M. Theory of Elasticity, 3rd ed., Butterworth-Heinemann, 1986. Comprehensive treatment of elasticity theory, stress-strain relations, and elastic waves.

3. Truesdell, C. and Noll, W. The Non-Linear Field Theories of Mechanics, 3rd ed., Springer, 2004. The definitive axiomatic treatment of continuum mechanics and constitutive theory.

4. Batchelor, G. K. An Introduction to Fluid Dynamics, Cambridge University Press, 2000. A thorough and readable treatment of viscous and inviscid flows.

5. Segel, L. A. Mathematics Applied to Continuum Mechanics, SIAM, 2007. Covers the mathematical foundations, including tensor analysis and variational methods for continua.

6. Soper, D. E. Classical Field Theory, Dover, 2008. Bridges continuum mechanics to relativistic field theory, emphasizing the Lagrangian density formalism.