Lamé Transmission Problem in Elasticity

• Lamé transmission problem is a fundamental concept in linear elasticity, enforcing displacement and traction continuity across interfaces with differing Lamé parameters.
• It involves boundary integral formulations and jump relations using single- and double-layer potentials for an accurate representation of elastic fields.
• The analysis extends to inverse boundary determination, dynamic modeling, and high-contrast scattering, highlighting its importance in theoretical and applied mechanics.

The Lamé transmission problem arises in the analysis of elastic media composed of domains with differing Lamé parameters, where the displacement field and its associated tractions must satisfy continuity conditions across interfaces. This foundational problem in linear elasticity underpins advances in inverse problems, resonance theory, and numerical modeling, with central importance for both theoretical and applied mechanics in heterogeneous and composite materials.

1. Mathematical Formulation and Transmission Conditions

Consider a bounded domain $\Omega \subset \mathbb{R}^m$ (typically $m=2$ or $3$), partitioned into a finite number of subdomains $\Omega_j$ with piecewise-constant or piecewise-smooth Lamé parameters $(\lambda_j, \mu_j)$. The displacement field $u : \Omega \to \mathbb{R}^m$ satisfies the system: $$L u := \mu \Delta u + (\lambda + \mu) \nabla(\nabla \cdot u) = 0$$ within each subdomain, subject to boundary or initial conditions, and interface/transmission conditions on each surface $$\Sigma_{jk} = \partial \Omega_j \cap \partial \Omega_k$$. Across each interface, the continuity requirements are: $$[u] := u|_{\Omega_j} - u|_{\Omega_k} = 0, \qquad [\sigma(u) n] = 0,$$ where $$\sigma(u) = \lambda (\operatorname{div} u) I_m + 2\mu \varepsilon(u)$$ and $\varepsilon(u) = (\nabla u + (\nabla u)^\top)/2$. In time-dependent (elastodynamic) formulations, similar conditions are imposed for both displacement and traction traces across jump interfaces at each time $t>0$ (Beretta et al., 2013, Caday et al., 2019).

2. Jump Relations, Boundary Integral Formulations, and Clifford Analysis

A comprehensive treatment of the jump problem for generalized Lamé-Navier systems employs Clifford algebra and the Euclidean Dirac operator, offering a first-order factorization of the Lamé operator. Solutions can be represented as

$$u(x) = D[\psi](x) - S[\varphi](x), \quad x \in \Omega_{+} \cup \Omega_{-}$$

where $S[\varphi]$ and $D[\psi]$ are single-layer and double-layer potentials constructed from the fundamental solution $\Gamma(x)$. The Plemelj-Sokhotski jump relations yield that the double layer operator introduces a jump in displacement, while the single layer is continuous but its traction has a jump: $$D[\psi]_{\pm}(x) = (\pm \tfrac{1}{2} I + K)[\psi](x), \quad \sigma_n^+ S[\varphi](x) - \sigma_n^- S[\varphi](x) = -\varphi(x)$$ Solving the transmission problem then reduces to a system of boundary integral equations for the unknown densities, valid even for admissible domains with non-smooth or fractal boundaries (Santiesteban et al., 7 Nov 2025).

3. Uniqueness, Stability, and Inverse Boundary Determination

Uniqueness and Lipschitz stability for the determination of piecewise constant Lamé parameters $(\lambda, \mu)$ from (local or global) Dirichlet-to-Neumann (DtN) maps have been established. For subdomains with Lipschitz interfaces, if two tensors yield the same (local) DtN map over a patch $E$, then their Lamé parameters coincide in all subdomains. The inverse mapping from the DtN operator to the parameter vector $L$ is globally Lipschitz—quantitatively,

$$\| \lambda - \lambda' \|_\infty + \| \mu - \mu' \|_\infty \leq C\| \Lambda_{\lambda, \mu} - \Lambda_{\lambda', \mu'} \|_{\mathcal{L}(H^{1/2}, H^{-1/2})}$$

with the constant depending on convexity bounds for the parameters and geometric features of the partition (Beretta et al., 2013). Key proof elements include Alessandrini's identity, quantitative three-sphere inequalities, and finite-dimensional inverse function theorems.

4. Dynamic Boundary Data and Microlocal Layer-Stripping

When Lamé parameters are piecewise smooth with jumps across a finite set of smooth hypersurfaces, dynamic (time-dependent) boundary measurements determine the wave speeds $c_P = \sqrt{(\lambda + 2\mu)/\rho}$ and $c_S = \sqrt{\mu/\rho}$ via the dynamic Dirichlet-to-Neumann operator. Under a convex foliation condition—requiring strictly geodesically convex level sets containing all interfaces—the exterior measurement operator determines $(\lambda, \mu)$ up to uniqueness. A microlocal layer-stripping procedure proceeds by reconstructing the parameters on successive convex "leaves," utilizing the lens relation, local boundary rigidity, and the reflection operator's principal symbol. This yields a constructive and stable algorithm for full recovery of the Lamé parameters from partial dynamic data (Caday et al., 2019).

5. High-Contrast Transmission and Scattering Resonances

In the high-contrast regime—where the interior-to-exterior Lamé moduli and densities scale as $1/\tau$ with $\tau \to 0$—the transmission problem exhibits a rich structure of scattering resonances. Near each nonzero Neumann eigenvalue of the interior Lamé operator, there are clusters of resonances; generically, the lifetime (imaginary part) is of order $\tau$, but in a codimension-one admissible set, resonances of order $\tau^2$ (long-lived subwavelength resonances) occur. Near $\omega=0$, the leading-order behaviour transitions between monopole and dipole types, determined by explicit criteria on the effective matrices derived from the interior Neumann data. Precise asymptotic expansions for the resolvent and resonance locations enable a quantitative understanding of resonant elastic scattering in composites and metamaterials (Li et al., 15 Jan 2026).

6. Transmission Problems with Damping and Numerical Analysis

One-dimensional Lamé transmission problems incorporating localized fractional time-derivative damping (in the Caputo sense) admit a semigroup framework showing strong stability and optimal polynomial energy decay. With parameters $\alpha\in(0,1)$, the decay rate is $E(t)\lesssim (1+t)^{-1}$ for $\eta=0$, or $E(t)\lesssim (1+t)^{-2/(1-\alpha)}$ for $\eta>0$. Numerical simulation via a finite-volume discretization in space, Newmark-$\beta$ and Crank-Nicolson schemes in time, and augmented memory variables for the fractional derivatives accurately reproduces the theoretical predictions for both wave propagation and long-time decay. Transmission conditions at the interface ensure continuity of displacement and flux, and correct reflection/transmission phenomena are confirmed numerically (Ammari et al., 5 May 2025).

7. Generalizations, Regularity, and Domains

All functional analytic results for the layer potentials, Teodorescu transform, and transmission boundary operators remain valid for Lipschitz as well as $C^{1,\alpha}$ boundaries. Further, the framework admits extension to fractal boundaries via $d$-summable sets in $\mathbb{R}^m$, preserving continuity and boundedness properties of the relevant operators and yielding well-posedness for generalized transmission/jump problems in highly irregular geometries (Santiesteban et al., 7 Nov 2025).