Nonlinear Elastic Interfaces

• Nonlinear elastic interfaces are lower-dimensional subsets in bulk materials where a nonlinear constitutive law governs complex phenomena like bending, buckling, and surface tension effects.
• The topic details variational formulations and surface-polyconvex models that couple bulk and interface mechanics, addressing anisotropy, higher-order derivatives, and contact nonlinearity.
• Computational approaches such as finite element, isogeometric, and FE-BE coupling are employed for accurate simulation and experimental parameter identification of interface behaviors.

A nonlinear elastic interface is a lower-dimensional subset (typically a surface or curve) embedded in a bulk elastic material across which the mechanical response is governed by a generally nonlinear constitutive law. These interfaces capture a wide class of phenomena including surface energy effects in soft solids, strongly nonlinear transmission or contact conditions, and singular geometric responses such as bending or buckling. The mathematical and computational treatment of such interfaces involves variational principles, advanced constitutive models, and often specialized numerical discretizations to resolve the coupling of bulk and interfacial mechanics, including anisotropy, inhomogeneity, and higher-order geometric terms.

1. Variational Formulations and Surface Energy Models

The rigorous treatment of nonlinear elastic interfaces begins with a variational framework. Let $\Omega \subset \mathbb{R}^3$ be the reference domain and $\Gamma \subset \partial\Omega$ the interface. A deformation $\varphi: \Omega \to \mathbb{R}^3$ describes the bulk configuration, with bulk and surface deformation gradients $F = \nabla\varphi$ in $\Omega$ and $F_s = \nabla_\Gamma\varphi = P F$ on $\Gamma$, where $P = I - N \otimes N$ is the tangential projector and $N$ is the unit normal.

The total stored energy with a nonlinear elastic interface is

$$E[\varphi] = \int_\Omega W_b(F) \, dV + \int_{\Gamma} W_s(F_s) \, dA$$

where $W_b$ is the bulk energy density (often polyconvex in $F$), and $W_s$ is a surface energy density encoding the interfacial mechanics. Surface-polyconvexity extends the classical notion of polyconvexity to $W_s$ by requiring existence of a convex, lower semi-continuous function $\Psi_s$ in suitable minors (e.g., $F_s$, $\mathrm{Cof} \ F_s$, $J_s = \det F_s$) such that $W_s(F_s) = \Psi_s(F_s, \mathrm{Cof}\ F_s, J_s)$, with appropriate growth and coercivity conditions for well-posedness (Horák et al., 31 Mar 2025).

These models yield coupled Euler–Lagrange equations, enforcing in the bulk $\mathrm{Div}\ S = 0$ and on the interface $S N - \mathrm{Div}_\Gamma S_s = 0$, where $S = \partial W_b/\partial F$ and $S_s = \partial W_s/\partial F_s$ are the (bulk and surface) first Piola–Kirchhoff stresses.

2. Constitutive Structure of Interface Laws

Nonlinear elastic interfaces exhibit a rich variety of constitutive behaviors:

• Isotropic and anisotropic surface tension or elasticity: For example, the Gurtin–Murdoch model posits $W_s(F_s) = \gamma J_s$, penalizing only area stretch, while a surface neo-Hookean law may take $$W_s(F_s) = \mu_s ( \tfrac12 \operatorname{Tr} C_s - 1 - \ln J_s ) + \frac{\lambda_s}{2}(J_s - 1)^2$$, with $C_s = F_s^T F_s$. Orthotropic elasticity further involves specific principal stretches $\lambda_1, \lambda_2$ along preferred surface directions (Horák et al., 31 Mar 2025).
• Fractional and memory effects: For interfaces coupled to a viscous bulk, fractional-derivative dynamics can arise, as in the nonlinear Lucassen–Kappler equation for viscous-coupled surface waves,

$$\left[ K^{(0)} + K^{(2)} \left(\bar{a}(1 + \partial_x U) - a_0 \right)^2 \right] \partial_x^2 U(x,t) = \mu D_t^{(3/2)} U(x,t)$$

with $D_t^{(3/2)}$ the Caputo $3/2$-order derivative, and $K_{2D}(a)$ is a strongly nonlinear area modulus reflecting phase transitions in experimental monolayers (Kappler et al., 2017).

• Contact and friction: For interfaces modeling frictional or adhesive contact, strongly nonlinear, possibly set-valued, interface laws involving Coulomb friction or general hemivariational inequalities are used, as in transmission problems coupling nonlinear p-Laplacian-type models in the interior to Laplacian or other fields in the exterior, with intricate interface laws (0911.1433, Gwinner, 2022, Gwinner et al., 2021).

3. Geometric and Higher-Order Effects

Many nonlinear elastic interfaces involve geometric functionals with higher-order derivatives or curvature dependence.

• Fourth/fifth-order geometric evolution: For example, the Helfrich bending energy for a 2D membrane gives $$E_H[\Gamma] = \frac12 \nu_0 \int_\Gamma \kappa^2 ds$$ with $\kappa$ the signed curvature, leading to interface conditions involving $\kappa_{ss}$ (where subscript $s$ is arclength), rendering governing equations of fourth order (Lu et al., 2019). In two-phase Muskat-type problems with elastic sheets, the bending energy generates a fifth-order quasilinear parabolic PDE in the interface graph $\eta$ (Wan et al., 4 Jan 2026).
• Curvature-driven instabilities: Nonlinear geometric terms enable the quantitative description of complex phenomena such as Rayleigh–Plateau instability, necking, capillarity-driven breakup, and tip-splitting morphogenesis in both engineered and biological materials (Horák et al., 31 Mar 2025, Lu et al., 2019).

4. Computational Methods

Nonlinear elastic interfaces require tailored discretization and solution techniques:

• Finite Element & Isogeometric Analysis: Variational formulations with surface–polyconvex densities are discretized using standard volume finite elements in the bulk, coupled with surface elements on interfaces or boundaries. The full energy is assembled across bulk and surface contributions and differentiated to yield residuals and tangent operators, ensuring correct coupling of surface and bulk mechanics. Newton-type solvers and mesh adaptivity (including gradient-recovery and residual estimators) are used to handle strong solution gradients near interfaces (Horák et al., 31 Mar 2025, 0911.1433).
• FE–BE Coupling and Hemivariational Formulations: Contact and transmission problems (especially with unbounded exterior domains) are handled by coupling finite elements in the interior to boundary element methods on the interface, employing boundary-integral representations (e.g., Steklov–Poincaré operators) and reducing the problem to finite-dimensional, possibly non-smooth (hemivariational) inequalities (Gwinner et al., 2021, Gwinner, 2022).
• Specialized Time Integration/Nonlocal Solvers: For high-order or fractional interface models, semi-implicit, spectral, or quadrature-based discretizations are implemented to efficiently handle stiff terms or memory kernels (Kappler et al., 2017, Lu et al., 2019).

5. Physical Phenomena and Applications

Nonlinear elastic interfaces underlie a diverse class of physical behaviors:

• Surface-dominated instabilities: The fidelity of surface-polyconvex models is critical for replicating surface-dominated phenomena such as the liquid bridge instability, Rayleigh–Plateau instability, and tip-splitting in soft gels, as demonstrated by quantitative agreement of computed critical aspect ratios and fastest-growing mode predictions with classical theory and experiment (Horák et al., 31 Mar 2025, Lu et al., 2019).
• Switching and thresholding: Nonlinearities in interfacial compressibility and modulus underpin amplitude-dependent propagation and threshold “switching,” as observed in phospholipid monolayers where wave propagation lengths increase sharply above a threshold amplitude and wave speed non-monotonically depends on compression (Kappler et al., 2017).
• Contact nonlinearity and harmonic generation: Nonlinear stick–slip and contact laws cause the transmission and reflection of higher harmonics in ultrasonic experiments, enabling quantitative inference of frictional and dissipative interface properties through universal scaling laws (Meziane et al., 2011, Gao et al., 2018).
• Fluid-structure and biological mechanics: Nonlinear interface models are foundational in biointerface problems (e.g., traction force microscopy, tumor-host interface evolution), as well as in fluid-structure interaction, where global existence and exponential decay results have been established for coupled nonlinear elasticity–fluid systems with boundary dissipation (Sarnighausen et al., 2024, Lu et al., 2019, Qin et al., 2018).

6. Analytical Well-Posedness and Optimal Control

The mathematical well-posedness of nonlinear elastic interface models depends on monotonicity, convexity, and coercivity properties at both the bulk and interface levels:

• Existence and uniqueness: Polyconvex and surface-polyconvex energy densities provide lower-semicontinuity and coercivity, guaranteeing the existence of minimizers for the variational problems associated with nonlinear interfaces (Horák et al., 31 Mar 2025, Sarnighausen et al., 2024). For transmission/contact problems, strong monotonicity and smallness conditions (relative to nonsmooth, frictional, or set-valued terms) assure unique solvability even for hemivariational inequalities with nonmonotone friction (0911.1433, Gwinner, 2022, Gwinner et al., 2021).
• Regularization and convergence: Non-differentiable contact/friction laws are regularized via smooth approximants, with rigorous error and convergence analyses as the regularization parameter vanishes, ensuring numerical tractability (Gwinner et al., 2021).
• Optimal control: Well-posedness and stability under parameter variation enable optimal control frameworks for nonlinear elastic interface problems, including distributed, boundary, and obstacle controls, with Mosco-convergence arguments ensuring stability and convergence of minimizers (Gwinner, 2022).

7. Experimental and Inverse Problems

Emerging techniques in soft matter and biomechanics increasingly probe interfacial nonlinearities via direct experiments or inverse problems:

• Parameter identification: Inverse approaches to traction force microscopy in nonlinear hyperelastic substrates leverage the well-posedness and differentiability of parameter-to-state maps for reconstructing traction fields from observed displacements, with algorithms based on Tikhonov regularization and finite-element discretizations yielding quantitative reconstructions under realistic noise (Sarnighausen et al., 2024).
• Nondestructive evaluation: The sensitivity of solitary wave reflection and harmonic generation to interface properties enables advanced diagnostic tools for defect detection, friction quantification, and material characterization in engineering and geophysical systems (Yang et al., 2010, Meziane et al., 2011, Gao et al., 2018).

---