Parabolic Inertia Lamé Operators

• Parabolic inertia Lamé systems are defined by coupling the classical Lamé elasticity operator with a parabolic time derivative and an inertial transport term, applicable in both linear and nonlinear settings.
• They facilitate the study of fractional operator formulations, unique continuation properties, and ill-posed mixed boundary-value problems in compressible and incompressible flow models.
• These models inform penalization approaches for incompressible Navier–Stokes equations and extend to solving complex eigenvalue problems in thin-shell vibration analyses.

“Parabolic inertia Lam” most naturally refers, in the cited mathematical literature, to time-dependent Lamé-type systems in which the spatial Lamé operator is combined with a parabolic time derivative and, in one recent usage, with the nonlinear inertial transport term $(u\cdot\nabla)u$. Across the cited literature, the expression does not denote a single universally fixed operator. Instead, it spans at least three closely related constructions: the linear parabolic Lamé operator $\mathbb H=\partial_t-L$ used to define fractional powers and extension problems (Banerjee et al., 2022); the parabolic Lamé-type operator arising from linearization of the compressible Navier–Stokes equations (Puzyrev et al., 2019); and the “parabolic inertia Lamé equations” $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$, introduced as a penalized auxiliary model for incompressible Navier–Stokes (Liu, 24 Jul 2025). A separate, older usage appears in thin-shell elasticity, where “Lamé system” refers to the three-dimensional elastic eigenproblem and “parabolic” refers to shell geometry rather than PDE type (Chaussade-Beaudouin et al., 2016).

1. Terminology and principal usages

The term “Lamé” refers to the isotropic elasticity operator built from the Lamé parameters $\mu$ and $\lambda$. In the linear parabolic setting studied by Banerjee and Senapati, one defines

$$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$

with $\mu\ge\delta_0$ and $2\mu+\lambda\ge\delta_0>0$, so that $L$ is strongly elliptic (Banerjee et al., 2022). In the mixed-boundary-value setting of Puzyrev and Shlapunov, the operator is written

$L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$

where $\mathbb H=\partial_t-L$0 is the Lamé operator in divergence form and $\mathbb H=\partial_t-L$1 collects lower-order terms coming from linearization of compressible Navier–Stokes (Puzyrev et al., 2019). In the nonlinear penalization framework of the 2025 paper, the “parabolic inertia Lamé” system is

$\mathbb H=\partial_t-L$2

on $\mathbb H=\partial_t-L$3 (Liu, 24 Jul 2025).

Usage	Model	Source
Fractional parabolic Lamé operator	$\mathbb H=\partial_t-L$4	(Banerjee et al., 2022)
Parabolic Lamé-type operator	$\mathbb H=\partial_t-L$5	(Puzyrev et al., 2019)
Parabolic inertia Lamé equations	$\mathbb H=\partial_t-L$6	(Liu, 24 Jul 2025)
3D Lamé system for parabolic shells	Elastic eigenproblem on thin shells	(Chaussade-Beaudouin et al., 2016)

This multiplicity of usage matters technically. In (Banerjee et al., 2022) and (Puzyrev et al., 2019), the operator is linear and parabolic in the standard PDE sense. In (Liu, 24 Jul 2025), the same Lamé diffusion is supplemented by the Navier–Stokes-type inertial nonlinearity. In (Chaussade-Beaudouin et al., 2016), by contrast, “parabolic” describes developable shell geometry, namely cylinders and truncated cones, while the Lamé system is the three-dimensional elasticity eigenproblem rather than a time-evolution equation.

2. Core operator structure

In isotropic elasticity, the spatial Lamé operator is

$\mathbb H=\partial_t-L$7

acting on vector fields $\mathbb H=\partial_t-L$8. The corresponding parabolic operator is

$\mathbb H=\partial_t-L$9

Its full Fourier transform satisfies

$$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$0

and the matrix symbol $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$1 is self-adjoint and positive-definite for $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$2 real (Banerjee et al., 2022). This yields a direct definition of fractional powers,

$$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$3

For the linearized compressible-flow setting, the elliptic part is written

$$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$4

and the lower-order operator is

$$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$5

The coefficients satisfy

$$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$6

on $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$7 (Puzyrev et al., 2019). This formulation makes explicit the relation between the parabolic Lamé-type operator and the linearization of the compressible Navier–Stokes momentum balance.

The nonlinear “parabolic inertia Lamé” equation retains the same Lamé diffusion but adds the inertial transport $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$8: $$\partial_tu-\mu\Delta u-(\lambda+\mu)\,\mathrm{grad}\,\mathrm{div}\,u+(u\cdot\nabla)u=0$$9 The intended interpretation is penalization of incompressibility by the large coefficient $\mu$0 in front of $\mu$1, with $\mu$2 fixed and $\mu$3 (Liu, 24 Jul 2025).

3. Fractional powers, extension problems, and unique continuation

A central development for the linear operator $\mu$4 is the analysis of its fractional powers $\mu$5, defined both by Fourier symbol and by semigroup subordination. The Bochner–Balakrishnan representation is

$\mu$6

where $\mu$7 and $\mu$8 (Banerjee et al., 2022).

The associated extension problem is posed on the half-space $\mu$9 with $\lambda$0: $\lambda$1 An explicit Poisson-kernel representation is obtained through the heat kernel $\lambda$2 of $\lambda$3 and the kernel $\lambda$4, leading to

$\lambda$5

The construction is accompanied by energy bounds,

$\lambda$6

(Banerjee et al., 2022).

The same work reduces the vector-valued extension to a degenerate system for

$\lambda$7

proves boundary $\lambda$8 and Hölder estimates, and establishes a space-like strong unique continuation theorem for

$\lambda$9

when $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$0 and $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$1 is bounded. In that regime, infinite-order vanishing at a point implies $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$2 (Banerjee et al., 2022). This places fractional parabolic Lamé operators within the modern extension-problem and frequency-function framework previously developed for scalar nonlocal parabolic equations.

4. Boundary-value theory and ill-posed mixed Cauchy problems

For the parabolic Lamé-type operator

$$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$3

in the finite cylinder $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$4, the mixed problem considered in (Puzyrev et al., 2019) prescribes both displacement and boundary stress on an open boundary portion $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$5: $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$6 The stress tensor is

$$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$7

The main conclusions are sharply dichotomous. First, uniqueness holds: if $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$8 contains at least one interior point of $$L u=\mu\Delta u+(\mu+\lambda)\nabla\!\operatorname{div}u,\qquad \mathbb H:=\partial_t-L,$$9 and

$\mu\ge\delta_0$0

then $\mu\ge\delta_0$1 in $\mu\ge\delta_0$2 (Puzyrev et al., 2019). Second, the problem is ill-posed in the natural spaces of smooth functions and in the corresponding Hölder spaces: there is no continuous dependence of $\mu\ge\delta_0$3 on the Cauchy data, and adding initial data $\mu\ge\delta_0$4 does not restore well-posedness.

The solvability criterion is expressed through layer potentials built from a fundamental solution $\mu\ge\delta_0$5. If one can find an extension $\mu\ge\delta_0$6 satisfying

$\mu\ge\delta_0$7

and matching the boundary sum of the volume, single-layer, and double-layer potentials of the prescribed data, then the solution is represented by

$\mu\ge\delta_0$8

This gives necessary and sufficient solvability conditions and shows that the parabolic Lamé-type Cauchy problem behaves much like classical elliptic Cauchy problems: uniqueness without stability (Puzyrev et al., 2019).

5. The nonlinear “parabolic inertia Lamé equations” and the Navier–Stokes penalization limit

The 2025 paper introduces the initial-value problem

$\mu\ge\delta_0$9

for $2\mu+\lambda\ge\delta_0>0$0, $2\mu+\lambda\ge\delta_0>0$1, and $2\mu+\lambda\ge\delta_0>0$2, and refers to it as the “parabolic inertia Lamé” system (Liu, 24 Jul 2025). The operator

$2\mu+\lambda\ge\delta_0>0$3

generates an analytic semigroup $2\mu+\lambda\ge\delta_0>0$4 on $2\mu+\lambda\ge\delta_0>0$5 with smoothing estimates uniform in $2\mu+\lambda\ge\delta_0>0$6, including

$2\mu+\lambda\ge\delta_0>0$7

The nonlinear problem is rewritten in Duhamel form,

$2\mu+\lambda\ge\delta_0>0$8

and the paper states that for any $2\mu+\lambda\ge\delta_0>0$9, $L$0, the problem admits a unique global smooth solution

$L$1

for each $L$2 (Liu, 24 Jul 2025).

The key connection to incompressible Navier–Stokes is the formal limit $L$3. The incompressible system is

$L$4

The paper argues that large $L$5 enforces incompressibility by penalization, proving uniform-in-$L$6 estimates such as

$L$7

and weak convergence

$L$8

A subsequence $L$9 is then extracted so that

$L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$0

and the limit $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$1 is asserted to satisfy the incompressible Navier–Stokes equations globally and smoothly in $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$2 (Liu, 24 Jul 2025).

The paper explicitly states that no smallness assumption on $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$3 is imposed and that the result is global-in-time in $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$4. It also notes that, in the classical literature, global smoothness in 3D is a famous open problem, and that the paper’s result, if correct, would resolve the global regularity question for 3D Navier–Stokes in the affirmative under Schwartz-class initial data (Liu, 24 Jul 2025). This makes the work mathematically significant but also places it in a context where independent verification is indispensable.

6. Distinct elastic-shell usage: parabolic shells and the 3D Lamé vibration system

A different use of “parabolic” and “Lamé” appears in thin-shell elasticity. There, the three-dimensional Lamé system is the elastic eigenvalue problem on a shell of thickness $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$5, with laterally clamped boundary conditions, and “parabolic” refers to the shell midsurface satisfying $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$6, hence $L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$7; the two principal examples are the circular cylinder ($L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$8) and the truncated cone ($L_{n+1}u=\partial_tu-L_nu-A(x,t)u=f(x,t),$9) (Chaussade-Beaudouin et al., 2016).

The 3D Lamé operator is written

$\mathbb H=\partial_t-L$00

and the first eigenvalue $\mathbb H=\partial_t-L$01 gives the square of the lowest vibration frequency. Via angular Fourier decomposition and reduction from the Koiter shell model to a scalar normal-component operator, the parabolic case leads to

$\mathbb H=\partial_t-L$02

with $\mathbb H=\partial_t-L$03 a positive fourth-order operator and

$\mathbb H=\partial_t-L$04

Optimization over the angular frequency $\mathbb H=\partial_t-L$05 gives the power law

$\mathbb H=\partial_t-L$06

so that

$\mathbb H=\partial_t-L$07

The first mode is therefore non-axisymmetric and oscillates strongly in the angular variable as $\mathbb H=\partial_t-L$08 (Chaussade-Beaudouin et al., 2016).

This shell-theoretic usage is not a parabolic evolution equation, but it is directly relevant to the phrase “parabolic Lamé” because it couples the Lamé elastic system with a parabolic geometric class and an inertia interpretation through free vibrations. A plausible implication is that searches for “parabolic inertia Lamé” can conflate two distinct literatures: one on parabolic Lamé operators in PDE analysis, and another on inertial vibrations of parabolic shells in elasticity.

7. Mathematical significance and conceptual boundaries

Taken together, the cited works place parabolic Lamé-type models at the intersection of elasticity, parabolic PDE theory, inverse problems, and fluid mechanics. The linear theory supports spectral, nonlocal, and boundary-value analyses: fractional powers, explicit extension formulas, weighted energy estimates, and unique continuation on one side (Banerjee et al., 2022), and integral representation, uniqueness, ill-posedness, and solvability criteria for mixed Cauchy data on the other (Puzyrev et al., 2019). The nonlinear theory, in the 2025 formulation, uses the same Lamé diffusion as a penalized, uniformly parabolic surrogate for incompressible Navier–Stokes (Liu, 24 Jul 2025).

At the same time, the phrase has clear conceptual boundaries. In the shell-vibration literature, the Lamé system is elliptic in space and enters an eigenvalue problem, while “parabolic” describes the geometry of the shell midsurface (Chaussade-Beaudouin et al., 2016). By contrast, in the operator-theoretic PDE literature, “parabolic” refers to the presence of $\mathbb H=\partial_t-L$09, and “inertia” may either be absent, encoded in nonlinear transport $\mathbb H=\partial_t-L$10, or appear only through a physical interpretation of the model.

For technical work, the most precise practice is therefore to specify the exact model: the fractional parabolic Lamé operator $\mathbb H=\partial_t-L$11 (Banerjee et al., 2022), the parabolic Lamé-type operator $\mathbb H=\partial_t-L$12 from compressible-flow linearization (Puzyrev et al., 2019), the nonlinear parabolic inertia Lamé penalization of incompressible Navier–Stokes (Liu, 24 Jul 2025), or the 3D Lamé shell-vibration system in the parabolic geometric class (Chaussade-Beaudouin et al., 2016).