Helmholtz separation of elastic waves Created 2026-09-24 Updated 2026-11-05
For a homogeneous isotropic solid with constant mass density and Lamé parameters , the Navier-Cauchy equation is
Taking the divergence, obeys with . Taking the curl, obeys with , because the curl of a gradient vanishes. These are the longitudinal P wave and transverse S wave speeds.
A suitable Helmholtz decomposition into a scalar-potential gradient and a divergence-free part separates the corresponding motions. In two-dimensional in-plane propagation, the P-SV displacement potentials implement this separation with two scalar wave equations. The constant-coefficient and isotropy assumptions matter: interfaces can couple P and SV modes, and arbitrary heterogeneous coefficients cannot be passed through the divergence and curl as if they were constants.
For a twice differentiable vector potential , commuting partial derivatives in gives the necessary condition
The potential is not unique: for any sufficiently smooth scalar , has the same curl, since the curl of a gradient is zero. This is the usual gradient freedom of a vector potential. A divergence-free condition alone need not suffice globally on every domain; here only its necessity is requested.
Stokes theorem states that for an oriented surface with compatible boundary orientation,
For the integrand in this problem, a convenient globally smooth vector potential is
The boundary of the upper ellipsoid is the ellipse . The outward normal on the upper surface induces counterclockwise orientation viewed from above, so parameterize it as , . Along this curve the first component of is zero and , leaving
Since the last angular integral is ,
The result is independent of the vertical semiaxis . As a sign check, the field is divergence-free, so closing the cap with the downward-oriented base disc makes its flux the positive integral of over that disc, which has the same value.
Integrate with respect to , obtaining . Matching gives , so . Finally, must equal the third component, so . Thus a scalar potential is
The fundamental theorem for line integrals now gives for either path, confirming both direct computations. Since the curl of a gradient is zero for this smooth function, it also confirms the result in part (b). The potential is defined throughout , so no domain obstruction to path independence occurs.
Each component depends only on its own coordinate, so all cross partial derivatives vanish. Therefore the curl is
This also follows from the identity that the curl of a gradient vanishes for a twice continuously differentiable potential energy.
The curl of a gradient vanishes because the mixed second partial derivatives of the smooth function agree. Therefore everywhere, and
Unlike part (i), this integrand does not involve the problematic auxiliary field.
The curl of the first vector field is
which is not identically zero. Thus is not conservative on and cannot equal a negative gradient there.
For the second field, direct differentiation shows
Therefore is conservative, and the requested negative-gradient potential is
Its curl vanishes, consistently with the curl of a gradient identity.
Taking each partial derivative of the scalar potential gives
The factor in the middle component comes from the chain rule. The potential is smooth, so the curl of a gradient identity gives .
Direct calculation of the curl gives
Since , the first vector field cannot be a gradient field: the curl of a gradient is zero. The second is a smooth polynomial field on all of the simply connected space , so its vanishing curl makes it a conservative vector field. A potential is also constructed explicitly in the next part. Thus is a conservative vector field and is not. A vanishing divergence alone, as for , is a different property and does not establish conservativity.
For a twice continuously differentiable scalar potential, the curl of its gradient is
Equality of mixed partial derivatives makes every component zero. Thus
This is the curl of a gradient identity.
Away from the origin, . Hence the curl of a gradient gives . Its divergence is the Laplacian of the radial coordinate, which in spherical coordinates is
Equivalently, differentiating in Cartesian coordinates gives .
In Cartesian coordinates, use the Levi-Civita symbol and the Einstein summation convention to expand
The contraction of the two symbols uses the Kronecker delta, and the final step is the product rule. Thus the convective acceleration identity is
Put . The curl of a gradient vanishes, and the curl of a cross product identity gives
Since the divergence of a curl is zero, taking the curl of the first identity yields
Twice continuous differentiability suffices for the mixed-derivative identities used here.
The divergence theorem states that a continuously differentiable vector field on a volume with piecewise smooth closed boundary satisfies
where is the outward unit normal. Take . Since the curl of a gradient is zero for a twice continuously differentiable scalar field, part (i) gives
Integrating this identity and applying the divergence theorem proves
Since the electric field is irrotational, introduce an electrostatic potential with . Gauss's law becomes the Poisson equation
The curl equation is automatic because the curl of a gradient vanishes. Independence of makes constant, and its required decay forces it to be zero, so a -independent potential suffices.
Velocity potential Created 2026-09-24 Updated 2026-10-13
Use a velocity potential to reconstruct an irrotational flow from one scalar field. A velocity potential satisfies
It has units and is unchanged as a description of the velocity field when any function of time alone is added. For twice continuously differentiable , the identity curl of a gradient proves zero vorticity. In incompressible flow, taking the divergence also gives , the Laplace equation. This turns a three-component velocity field into a single harmonic function with appropriate boundary conditions. [ref]
For uniform flow , take . With , the change from to is . By the fundamental theorem for line integrals, the line integral of velocity between these points equals that change along any path in the domain; a closed curve gives zero circulation.
For the potential line vortex, choose a local angle branch in polar coordinates and set . Then and . One counterclockwise turn changes the continued velocity potential by , although the velocity field returns to its original value. With , this change is , exactly the circulation. A full annulus therefore admits local angle branches and a multivalued velocity potential, but no global single-valued velocity potential when . Cutting the annulus along a radial line allows a single-valued angle branch on the remaining domain. The obstruction is path independence, not a failure of the local gradient formula.