Dynamic stress response kernels for dislocations and cracks: unified anisotropic Lagrangian formulation
Pith reviewed 2026-05-16 17:49 UTC · model grok-4.3
The pith
Dynamic stress kernels for anisotropic dislocations and cracks depend only on the prelogarithmic Lagrangian factor L(v) in Fourier space.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The stress-response convolution kernels for planar defects in anisotropic elastodynamics have Fourier representations that rest exclusively on the prelogarithmic Lagrangian factor L(v), while their space-time forms involve the prelogarithmic impulsion function p(v) = L'(v). The theory is obtained via the Stroh formalism, pays explicit attention to causality, and therefore covers subsonic, intersonic, and supersonic regimes of defect motion without modification. A direct consequence is that the stress appearing in the Weertman model of steadily moving dislocations can be rewritten solely in terms of L(v). The resulting kernels are immediately usable in phase-field Fourier-based codes for two-
What carries the argument
The prelogarithmic Lagrangian factor L(v) together with its derivative p(v) = L'(v), which alone determine the Fourier and space-time kernels.
If this is right
- The Weertman steady-motion dislocation stress can be expressed directly through L(v) rather than through separate wave integrals.
- The same kernel expression applies without alteration to subsonic, intersonic, and supersonic defect speeds.
- The kernels plug immediately into existing Fourier-based phase-field codes for planar defect systems in anisotropic media.
- Cracks and flat-core dislocations are treated by the identical kernel structure, differing only in the applied boundary condition.
Where Pith is reading between the lines
- Numerical schemes could precompute L(v) once per material and velocity range, replacing repeated full anisotropic wave solves at each time step.
- The same Lagrangian reduction may apply to three-dimensional or curved defects once an appropriate generalization of the Stroh representation is available.
- Energy-release-rate calculations in dynamic fracture could be simplified by substituting the p(v) form of the kernel into the dissipation integral.
- Recovery of the known isotropic generalized-function kernels in the appropriate limit would serve as a basic consistency check.
Load-bearing premise
The Stroh formalism, developed for static or steady-state problems, extends without change to the full dynamic convolution kernels while preserving causality in every velocity regime.
What would settle it
Direct numerical evaluation of the derived space-time kernel for an orthotropic crystal at a chosen supersonic velocity, followed by comparison with an independent elastodynamic Green's-function simulation of the same configuration.
read the original abstract
Elastodynamic cohesive-zone models for defects such as cracks or dislocations (such as the Geubelle-Rice model for cracks, or the Dynamic Peierls Equation for flat-core dislocations), feature the same stress-response convolution kernel in space and time. It accounts for in-plane elastic wave propagation, while its associated instantaneous radiative term accounts for radiative losses in the surrounding medium. These objects are well-known for isotropic elasticity, with their space-time representations involving generalized functions. For anisotropic elasticity they were unknown. The paper presents a derivation using the Stroh formalism. Their Fourier representation rests exclusively on the so-called prelogarithmic Lagrangian factor $L(v)$, while their space-time form involves its derivative $p(v)=L'(v)$, the prelogarithmic impulsion function. A straightforward consequence is the reformulation of the stress in the Weertman model of steadily-moving dislocations in terms of $L(v)$. Special care being paid to the causality constraint, the theory covers indifferently subsonic, intersonic and supersonic regimes of motion. The theory proposed is suitable to phase-field-type Fourier-based numerical codes for planar systems of defects in anisotropic elastodynamics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives dynamic stress-response convolution kernels for dislocations and cracks in anisotropic elasticity via the Stroh formalism. Fourier-space kernels depend only on the prelogarithmic Lagrangian factor L(v) from the velocity-dependent Stroh eigenvalue problem; space-time kernels follow from its derivative p(v)=L'(v). The formulation unifies subsonic, intersonic and supersonic regimes with explicit attention to causality and is positioned for phase-field Fourier-based numerical codes; a byproduct is a reformulation of the Weertman steady-motion model in terms of L(v).
Significance. If the central derivation holds, the result supplies a compact, parameter-free extension of the known isotropic elastodynamic kernels to general anisotropy. It enables direct use of existing Stroh solvers for dynamic defect modeling across all velocity regimes and directly supports numerical implementations in materials science. The Lagrangian framing and causality handling are strengths that could make the kernels a standard reference object.
major comments (2)
- [§3] §3 (velocity-dependent Stroh matrix and contour integration): the assertion that residues depend only on L(v) with no supplementary continuous-spectrum or radiation-condition contributions must be shown explicitly for intersonic and supersonic regimes, where the PDE changes type. The abstract states special care for causality, but the load-bearing step is whether the static Stroh procedure adapts without additional branch-cut adjustments; an explicit residue calculation for v > c_s is required to confirm the kernels remain causal.
- [§5] §5 (space-time kernels via p(v)=L'(v)): the radiative-loss term must be verified to reduce exactly to the known isotropic expression when anisotropy vanishes. Without this limit check, the claim that the formulation covers all regimes uniformly rests on an untested extrapolation.
minor comments (2)
- [Introduction] Notation: the distinction between L(v) (prelogarithmic Lagrangian) and the conventional energy factor should be stated once at first use to avoid confusion with prior literature.
- [Figure 2] Figure 2 (kernel plots): axis labels and velocity normalization should be made explicit so that subsonic/intersonic transitions are immediately visible.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We address each major point below, providing clarifications from the derivation and indicating where the manuscript will be revised to strengthen the presentation.
read point-by-point responses
-
Referee: [§3] §3 (velocity-dependent Stroh matrix and contour integration): the assertion that residues depend only on L(v) with no supplementary continuous-spectrum or radiation-condition contributions must be shown explicitly for intersonic and supersonic regimes, where the PDE changes type. The abstract states special care for causality, but the load-bearing step is whether the static Stroh procedure adapts without additional branch-cut adjustments; an explicit residue calculation for v > c_s is required to confirm the kernels remain causal.
Authors: The contour integration in §3 is performed in the complex plane with the iε prescription to enforce causality, enclosing only the poles associated with the Stroh eigenvalues whose residues yield the L(v)-dependent kernels; the radiation condition at infinity eliminates continuous-spectrum contributions even when the PDE changes type for v > c_s. We agree that an explicit residue evaluation for a supersonic example would make this transparent and will add it to the revised §3. revision: yes
-
Referee: [§5] §5 (space-time kernels via p(v)=L'(v)): the radiative-loss term must be verified to reduce exactly to the known isotropic expression when anisotropy vanishes. Without this limit check, the claim that the formulation covers all regimes uniformly rests on an untested extrapolation.
Authors: When the elastic tensor is taken isotropic, L(v) reduces to the standard isotropic prelogarithmic factor and p(v) = L'(v) reproduces the known radiative-loss term (including the correct coefficient involving the Rayleigh speed). We will insert this explicit reduction as a short verification in the revised §5 (or as a new appendix) to confirm the uniform coverage of regimes. revision: yes
Circularity Check
No circularity: kernels derived from established Stroh formalism as functions of pre-existing L(v)
full rationale
The derivation applies the standard Stroh eigenvalue problem (external to the paper) to obtain dynamic kernels whose Fourier form depends only on the known prelogarithmic factor L(v) and whose space-time form follows by differentiation to p(v)=L'(v). This is an independent reorganization and extension of the formalism to time-dependent convolution kernels across velocity regimes, with no step that redefines L(v) in terms of the kernels, fits parameters to data, or relies on self-citation chains for the central result. Causality is asserted via the construction but does not reduce the claim to its inputs by definition.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Stroh formalism applies to dynamic elastodynamic kernels for planar defects
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.