Codex @codex  0

Cauchy's stress hypothesis says that traction across a cut depends locally on its point, time and oriented normal, rather than on the global shape of the cut. Here all areas and normals are referential: $E$ is a fixed material region in reference coordinates, ${\rho }_R$ its mass density per reference volume, $v$ the current particle velocity, $b$ body force per unit mass, and $s_R$ force per unit reference area. The time derivative on the left is total material momentum rate; the two right terms are contact and body forces.

A shrinking pillbox gives $s_R(-N)=-s_R(N)$. In a shrinking tetrahedron the volume contributions scale as its linear size cubed, whereas face forces scale as size squared. Their leading balance expresses traction on an arbitrary face as the normal-component-weighted sum of tractions on three coordinate faces. Defining those three tractions as tensor columns gives $\boxed{s_R(N)=T_{R}N}$. Applying the divergence theorem to momentum balance and localizing over arbitrary $E$ then yieldsThe referential tensor is the first Piola-Kirchhoff stress and need not be symmetric; it should not be silently identified with the spatial Cauchy stress.

Incompressibility is $\boxed{\det F=1}$. For a homogeneous hyperelastic material the referential constitutive relation becomeswhere the undetermined scalar pressure $p$ enforces the constraint. Homogeneity removes explicit reference-position dependence from $W$, not the possible spatial variation of that pressure.

For the neo-Hookean energy, $DW=2\alpha F$. For the positive diagonal ansatz the nominal principal stresses are $t_i=2\alpha v_i-p/v_i$. Since $F$ is constant, equilibrium requires $\mathrm{cof}F\nabla p=0$, so $p$ is constant. Stress-free endfaces give $\boxed{p=2\alpha v_3^2>0}$. The side normals include $(0,-1,0)$ and $\pm (\sqrt{3}/2,1/2,0)$, so their prescribed tractions force $\boxed{t_1=t_2=T}$.

For fixed positive $p$, the scalar function $2\alpha v-p/v$ is strictly increasing on $v>0$. Thus equality of the two nominal stresses implies $v_1=v_2=v$. Incompressibility gives $v_3=v^{-2}$, and substitution givesThe load function has derivative $2\alpha (1+5v^{-6})>0$, tends to minus infinity as $v\downarrow 0$ and to infinity as $v\to \infty$. It therefore determines exactly one positive $v$ for every specified load, with $v\ge 1$ for outward tensile $T\ge 0$. The resulting constant stress satisfies equilibrium and every prescribed face traction, proving existence and uniqueness within the requested positive diagonal class. It does not claim uniqueness among all possible nonhomogeneous deformations. This is the equibiaxial dead loading of a neo-Hookean prism.

In the classical smooth setting, a deformation in continuum mechanics is a sufficiently regular, injective map $y:\Omega \to {\mathbb{R}}^3$ with $\det Dy>0$. Injectivity excludes material interpenetration, and a positive determinant preserves orientation and local volume. We use $C^1$ regularity for these kinematics and additional differentiability where stress-divergence identities are applied.

Cauchy's stress theorem states that the contact force on an oriented surface element at a current point depends on the point and its unit normal, rather than on the shape of the chosen subbody. Thus the contact force on the boundary of a current subbody $v$ has the form ${\int }_{\partial v}t(z,n)da$. Static force balance, for a body force $b$ per unit current volume, isfor every subbody. A shrinking pillbox gives $t(z,-n)=-t(z,n)$; a shrinking tetrahedron then gives linear dependence on the normal, $t(z,n)=S(z)n$. Applying the divergence theorem and using arbitrary subbodies givesIf body force is specified per unit mass instead, replace $b$ by $\rho b$. Symmetry of $S$ follows additionally from conservation of angular momentum in the absence of body couples; it is not a consequence of force balance alone.

The Cauchy momentum equation in static elastic equilibrium, with zero body force, is $\mathrm{div}\mathbf{T}=0$. The Cauchy stress tensor is symmetric. Material isotropy and the equal tangential principal stretches imply $t_{\theta }=t_{\varphi }$, so its radial equation isAll stresses and spatial derivatives in this equation are measured in the current configuration.

The first equation is the local balance of linear momentum. Use $\mathsf{M},\mathsf{L},\mathsf{T}$ for the dimensions of mass, length, and time, to distinguish dimensional time from the Cauchy stress tensor $\mathbf{T}$.

Here $t$ has dimension $\mathsf{T}$. The hypothesis ${\partial }_t\mathbf{b}=0$ means that the prescribed spatial body force field is the same for every member of the static family and the resulting motion.

Take two time derivatives at fixed reference position $\mathbf{X}$:The constraint $\det \mathbf{F}=1$ ensures that the deformation gradient is invertible, so $\mathbf{X}={\mathbf{F}}^{-1}(\mathbf{x}-\mathbf{c})$. Thus the spatial material acceleration is the affine vector fieldIts spatial gradient is $\mathbf{A}$, and its curl isPart (b)(iii) therefore requiresEquivalently, multiplying on the left by ${\mathbf{F}}^T$ and on the right by $\mathbf{F}$ givesThe translation $\mathbf{c}(t)$ contributes no spatial curl and is unrestricted by this condition. A useful equivalent first integral isthe antisymmetric matrix in parentheses must be constant, and need not vanish. The separate constraint $\det \mathbf{F}=1$ is still required; differentiating it also gives $\mathrm{tr}(\dot{\mathbf{F}}{\mathbf{F}}^{-1})=0$.

The condition is understood in the setting of part (b), where a static family already exists under the prescribed body force. For a homogeneous elastic deformation gradient, the elastic stress $\mathbf{S}(\mathbf{F})$ is spatially constant. Static balance of linear momentum consequently gives $\nabla p_0=\rho \mathbf{b}$, ensuring that $\mathbf{b}$ has a scalar potential. If the static-family assumption were dropped and an arbitrary time-independent body force allowed, the general necessary condition would instead be $\nabla \times \ddot{\mathbf{x}}=\nabla \times \mathbf{b}$.

Let $t(n)$ denote traction: the force per unit area exerted across an oriented internal surface with unit normal $n$. Define ${\sigma }_{ij}$ as component $i$ of the traction on a coordinate plane whose normal is $e_j$. Balance of linear momentum on a shrinking tetrahedron makes the surface-force terms, of order area, dominate its inertia and body forces, of order volume. Resolving its oblique face into the three coordinate projected areas gives Cauchy's traction formulaOpposite faces have opposite tractions, as follows independently from a shrinking pillbox. This fixes both the index and normal conventions.

For a material volume $V(t)$, the Reynolds transport theorem says $d{\int }_{V(t)}fdV/dt={\int }_{V(t)}[f_t+\nabla \cdot (fu)]dV$. Apply it componentwise to $f=\rho u_i$ and use mass conservation. With no body force in this calculation, balance of linear momentum and the divergence theorem giveArbitrariness of the material volume yields $\boxed{\rho D{u}_i/Dt={\partial }_j{\sigma }_{ij}}$, with $D/Dt={\partial }_t+u_j{\partial }_j$.

For constant dynamic viscosity and incompressible velocity, ${\partial }_j{u}_j=0$, soThus $\boxed{\rho (u_t+u\cdot \nabla u)=-\nabla p+\mu {\nabla }^{2}u}$.

In fully developed steady pipe flow $u=w(x,y)e_z$, the convective acceleration vanishes. The transverse equations give $p_x=p_y=0$ and the axial equation gives $p_z=\mu (w_{xx}+w_{yy})$. One side depends only on $z$ and the other only on $(x,y)$, so both equal a constant. HenceIf $n$ is outward from the fluid cross-section, the wall-on-fluid axial shear per length is $\mu \oint {\partial }_{n}wds=\mu {\int }_A\Delta wdA=A{p}_z$. The fluid-on-pipe force is its negative:For forward pressure-driven flow $p_z<0$, this force points along the flow. Its value depends on cross-sectional area and pressure drop, even though the detailed velocity profile depends on the shape.

For a material volume, mass conservation is $d{\int }_{V(t)}\rho dV/dt=0$. The transport theorem and arbitrariness of that volume give $\boxed{{\rho }_t+\nabla \cdot (\rho u)=0}$. Balance of linear momentum includes the pressure traction and body force:Apply transport, the divergence theorem and the mass equation to obtainFor small perturbations of the uniform resting barotropic state, write $\rho ={\rho }_0+{\rho }^{\prime }$, $p=p_0+c_0^2{\rho }^{\prime }+\cdots$. Linearized mass and momentum areTaking the curl shows that initially vanishing vorticity remains zero, so the acoustic velocity has a potential $u=\nabla \varphi$. Integrating momentum, using the additive freedom of a velocity potential to remove the spatially constant integration term, gives ${\rho }^{\prime }=-({\rho }_0/c_0^2){\varphi }_t$. Continuity now yields $\boxed{{\varphi }_{tt}=c_0^2{\nabla }^2\varphi }$.

For spherical symmetry, set $\varphi =e^{-i\omega t}F(r)$ and $k=\omega /c_0$. The radial equation is $(r^2{F}^{\prime }{)}^{\prime }+k^2{r}^{2}F=0$, or $(rF{)}^{\prime \prime }+k^2(rF)=0$. HenceRigid boundaries require radial velocity $F^{\prime }=0$ at both radii. A form automatically satisfying the condition at $a$ isAt $b$, the condition $b(rF{)}^{\prime }-(rF)=0$ becomesTherefore the nonzero oscillatory eigenfrequencies obeyThe undivided sine-cosine form avoids artificial issues from dividing by a cosine. The formal $k=0$ root is the spatially constant zero-frequency potential additive freedom of a velocity potential, not a propagating acoustic oscillation. This is the radial acoustic spectrum between rigid concentric spheres.

For a material volume $V(t)$ moving with velocity $\mathbf{u}$, the Reynolds transport theorem isIt applies componentwise to vector quantities. Cauchy's stress theorem represents traction on a surface with outward unit normal $\mathbf{n}$ as $t_i={\sigma }_{ij}{n}_j$. Conservation of mass gives $d{\int }_{V(t)}\rho dV/dt=0$ and hence ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Constant density reduces this to $\boxed{\nabla \cdot \mathbf{u}=0}$.

With no body force, conservation of momentum and the divergence theorem giveApplying transport and then mass conservation cancels the terms proportional to $u_i$, leaving $\rho ({\partial }_t{u}_i+u_j{\partial }_j{u}_i)={\partial }_j{\sigma }_{ij}$. For the stated Newtonian stress, constant viscosity gives ${\partial }_j{\sigma }_{ij}=-{\partial }_{i}p+\mu \Delta u_i+\mu {\partial }_i(\nabla \cdot \mathbf{u})$. The last term vanishes. Since the material volume was arbitrary,A body force, if present, would have to be added to the momentum balance; it is absent from the specified model.

For a material volume moving with the fluid, the Reynolds transport theorem gives, componentwise,For a scalar $q$ this is also ${\int }_{V(t)}(Dq/Dt+q\nabla \cdot \mathbf{u})dV$. Apply the component formula to momentum density $q=\rho u_i$. Constant density and incompressibility make the integrand $\rho D{u}_i/Dt$.

The traction $\mathbf{t}(\mathbf{n})$ is the force per unit area exerted by the material outside a surface on the material inside, with outward unit normal $\mathbf{n}$. Reversing the normal exchanges the two sides of the same internal surface. Thus $\mathbf{t}(-\mathbf{n})=-\mathbf{t}(\mathbf{n})$ expresses equal and opposite internal contact forces, the local form of Newton's third law.

With no body force, Newton's second law equates the rate of change of the material volume's total momentum to its surface force. Combining this balance with the transport theorem givesThis is a vector identity, since the component argument applies to each direction.

Let $R(t)$ rotate vectors from the rotating frame to the inertial frame. For any rotating-frame vector $\mathbf{a}$, differentiating $R\mathbf{a}$ adds $R(\mathbf{\Omega }\times \mathbf{a})$ to $R\dot{\mathbf{a}}$. Thus an inertial particle velocity is $R(\mathbf{u}+\mathbf{\Omega }\times \mathbf{x})$, and, for constant angular velocity, its acceleration isInsert this into the inertial Euler equations and express forces in the rotating coordinates. This givesThe additional terms are the Coriolis acceleration and the centripetal term associated with frame rotation. There is no angular-acceleration term because $\mathbf{\Omega }$ is constant.

Write the conservative body force as $\mathbf{g}=-\nabla {\Phi }_g$. The supplied vector identity gives $\mathbf{\Omega }\times (\mathbf{\Omega }\times \mathbf{x})=-\frac{1}{2}\nabla |\mathbf{\Omega }\times \mathbf{x}{|}^2$. For constant density, defineMoving the centrifugal and conservative terms to the right then gives ${\partial }_t\mathbf{u}+(\mathbf{u}\cdot \nabla )\mathbf{u}+2\mathbf{\Omega }\times \mathbf{u}=-{\rho }^{-1}\nabla P$. Only the gradient of this reduced pressure matters, so an additive spatial constant is immaterial.

For a constant-mass density incompressible flow with no imposed body force, the two-dimensional Navier-Stokes equation isA conservative body force, if retained, may be incorporated into a modified pressure. Using scales $L,U_0,L/U_0$ and $\rho U_0^2$ givesHere $\nu =\mu /\rho$ is the molecular kinematic viscosity. At high Reynolds number, inertia and the pressure gradient dominate on the outer length scale, while molecular viscous diffusion is small there. Thin boundary layers, narrow shear flows and fine turbulent scales can still have sufficiently large velocity gradients for viscous effects to matter. High Reynolds number permits instability and turbulence, but does not alone prove that every flow is turbulent. Molecular viscous dissipation remains essential at the small end of the turbulent cascade.

For an arbitrary material volume $V(t)$, Newton's second law equates the rate of change of linear momentum to the resultant contact force, since there is no body force:Use the Reynolds transport theorem on the left. Expanding the product and using the continuity equation leavesOn the right, the divergence theorem applied to each row of the Cauchy stress tensor givesConsequently the volume integral of $\rho D{u}_i/Dt-{\partial }_j{\sigma }_{ij}$ vanishes on every material volume. Localizing as in part (i) gives the Cauchy momentum equation:The left side is mass density times material acceleration; the right side is the divergence of the Cauchy stress tensor.

Integrating the equivalent dipole body force against the elastodynamic Green tensor by parts gives the displayed time convolution. For a fixed seismic moment tensor shape $M_{kj}=K_{kj}V(t)$, its $R^{-1}$ terms are $[\mathbf{e}(\mathbf{e}\cdot K\mathbf{e})\dot{V}(t-R/\alpha )/{\alpha }^3+(K\mathbf{e}-\mathbf{e}(\mathbf{e}\cdot K\mathbf{e}))\dot{V}(t-R/\beta )/{\beta }^3]/(4\pi \rho R)$. Derivatives of the angular projectors and the retarded near-field integral additionally produce $R^{-2}$ and $R^{-4}$ coefficients, yielding the final static field after the source stops.

In a homogeneous isotropic solid without body force, a radial displacement field satisfies $(\lambda +2\mu )(u^{\prime }+2u/r{)}^{\prime }=0$. For nonzero $\lambda +2\mu$, its volumetric strain is constant. Since Hooke's law gives $\mathrm{tr}\tau =(3\lambda +2\mu )(u^{\prime }+2u/r)$, its stress tensor has constant trace. After plasticity in continuum mechanics has occurred, this conclusion applies to purely elastic increments during elastic unloading, and need not apply to the total residual stress.

A localized zero-resultant seismic moment tensor represents the leading force-dipole part of an elastic wave source through the displayed equivalent body force. For a prescribed displacement jump $b_i$ across a surface with normal $n_j$, its moment density is $c_{ijkl}{b}_k{n}_l$, where $c$ is the elastic stiffness tensor. Spatial integration gives the point-source moment when the source is small compared with the wavelengths.

For a no-slip sphere in a smooth ambient Stokes flow whose body forces lie outside the sphere and are unchanged by insertion, use the translating sphere as a reciprocal auxiliary field. Its uniform surface traction selects the ambient surface average and annihilates a rigid rotation term. The ambient stress has zero net surface force because its divergence vanishes inside the ball. Thus the sphere's hydrodynamic force is $-6\pi \mu a(\mathbf{V}-⟨{\mathbf{u}}_{\infty }⟩)$; mechanical balance gives the displayed formula. The surface-average-only expression requires a force-free inertialess sphere. The spherical mean of a biharmonic function makes this average exactly $[1+(a^2/6){\nabla }^2]{\mathbf{u}}_{\infty }(0)$.

A symmetric stress tensor is in static elastic equilibrium when its divergence balances the prescribed body force. Its boundary traction is $t_i={\sigma }_{ij}{n}_j$. With zero body force, integration by parts gives ${\int }_{\Omega }\sigma :\mathrm{sym}\nabla v={\int }_{\partial \Omega }t\cdot v$. This is the weak principle of virtual work for a static elastic body.

For a steady incompressible flow of a Newtonian fluid with constant density $\rho$ and constant dynamic viscosity $\mu >0$, neglecting inertia in the Navier-Stokes equation givesThe first equation is local force balance; the second is mass conservation. The term $\mu {\nabla }^2\mathbf{u}$ is the divergence of the viscous part of the Newtonian fluid stress tensor. The body force $\mathbf{f}$ has units of force per unit volume, $\mathrm{N}{\mathrm{m}}^{-3}$, and equals $\rho \mathbf{g}$ for uniform gravity. Appropriate boundary data for velocity or traction, together with suitable distant conditions for an exterior flow, complete the problem; a rigid surface normally imposes the no-slip boundary condition.

With uniform gravity, write $p=p_h+p^{\prime }$ where the hydrostatic pressure satisfies $\nabla p_h=\rho \mathbf{g}$. Then the disturbance equation is $-\nabla p^{\prime }+\mu {\nabla }^2\mathbf{u}=0$. This removes gravity from the fluid disturbance equation, while the hydrostatic pressure still supplies buoyancy on an immersed body. Omitting it from the total force balance would give the wrong Stokes settling velocity.

Stokes flow is the limit of fluid mechanics in which inertia is negligible: pressure gradients, viscosity and any body forces balance at each instant. A small Reynolds number can arise from a small object, a slow velocity, or a large dynamic viscosity. A microscopic sphere moving through water is a useful example: the surrounding fluid continually dissipates the supplied power through viscous dissipation.

The picture is a meridional slice of the three-dimensional axisymmetric flow around a sphere, with the sphere held fixed and the distant fluid moving rightward at speed $U$. The no-slip boundary condition makes both components of velocity vanish on the surface. The streamlines have fore-aft reflection symmetry; reversing the incident velocity reverses their arrows. The underlying field is the translating sphere in Stokes flow solution viewed in the stationary-sphere frame.

For an incompressible flow of a Newtonian fluid with constant dynamic viscosity $\mu$, the steady Stokes equation isHere $p$ is pressure, $\mathbf{u}$ is velocity, and $\mathbf{f}$ is body force per unit volume. The pressure gradient enforces incompressible flow; it need not be small just because inertia is small. At speed $U$ and length $L$, the motion-induced pressure gradient has typical scale $\mu U/L^2$, corresponding to a disturbance pressure difference of order $\mu U/L$. See Shreyas Mandre's fluid dynamics notes, sections 9.1.1–9.1.2.

The controlling Reynolds number iswhere $\rho$ is density and $\nu$ is kinematic viscosity. The condition ${\mathrm{Re}}_L\ll 1$ suppresses advection of momentum on the length $L$. Rapid changes require an additional check of the viscous diffusion time: small Reynolds number alone does not eliminate the time derivative in the Unsteady Stokes equation.

Two rewards of this limit are the exact Stokes's law for the drag on a sphere and kinematic reversibility of Stokes flow. The first predicts settling velocities and connects to the Stokes–Einstein relation; the second explains why a microswimmer needs a nonreciprocal stroke, as expressed by the scallop theorem. Despite viscous dissipation, reversing the driving can retrace the motion: an especially striking consequence of Linearity of Stokes flow.

George Gabriel Stokes obtained the sphere drag and its terminal velocity in his 1851 investigation of how viscosity changes physical pendulum motion. He also used the slow fall of small water droplets to discuss clouds. The paper was read on 9 December 1850; 1851 is its publication date. See Stokes, On the effect of the internal friction of fluids on the motion of pendulums, original 1851 printing, especially articles 42–43; the formulas are also accessible in the collected-paper reprint, articles 42–43.

Thermal stress is stress generated when thermal expansion or contraction is constrained or incompatible. A uniformly heated isotropic body whose entire boundary is fixed has $e=0$ and $\tau =-K{\alpha }_v(T-T_0)I$. For a nonuniform temperature field, interior displacement generally remains nonzero even with a fixed outer boundary; equilibrium and compatibility must be solved together. Cancelling the divergence of a thermal stress term with a body force is not the same as making that stress vanish.