Parabolic-Hyperbolic Skeleton Equation

• The parabolic-hyperbolic skeleton equation is a structural framework coupling diffusion and wave dynamics, defined by different behaviors across interfaces or via singular limits.
• Its formulations span canonical mixed-type PDEs, fractional diffusion–hyperbolic systems with nonlocal gluing, and controlled hydrodynamic equations in large-deviation theory.
• Explicit interfacial reduction, energy estimates, and asymptotic analysis ensure well-posedness and unique decay rates, offering practical insights into variational and dynamic regimes.

The parabolic-hyperbolic skeleton equation denotes a class of canonical or limiting equations in which parabolic and hyperbolic dynamics are coupled either across a type-change interface, through a nonlocal gluing law, or via a singular-perturbation or large-deviation reduction. In the mixed-type boundary-value literature, the canonical model is the piecewise equation

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0, \end{cases}$$

posed on a domain partitioned by the interface $\Gamma=\{y=0\}$; in the non-equilibrium large-deviations literature, the skeleton equation is

$$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$

and in the weakly dissipative Kirchhoff setting, the parabolic limit

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$

acts as the skeleton for a hyperbolic equation with small inertia (Shishkina et al., 2023, Fehrman et al., 2019, Ghisi et al., 2012). This suggests that the expression is contextual rather than universal.

1. Terminological scope and model classes

Within the supplied literature, the term “skeleton equation” is used explicitly in the large-deviation analysis of conservative interacting particle systems, where it names a degenerate parabolic-hyperbolic PDE with irregular drift on $\mathbb T^d\times(0,T)$ (Fehrman et al., 2019). By contrast, the 2023 boundary-value study does not use the term and instead treats a canonical mixed parabolic–hyperbolic PDE in a bounded planar domain split by $y=0$ (Shishkina et al., 2023). A fractional analogue replaces the parabolic part by a Caputo fractional diffusion equation and couples it to a hyperbolic equation by an integral gluing condition on the interface segment $AB=\{(x,0):0<x<1\}$ (Karimov et al., 2013). In a different but related usage, the weakly dissipative Kirchhoff literature identifies the parabolic limit of a hyperbolic equation as the “skeleton” or leading-order dynamics in the singular limit $\varepsilon\to0$ (Ghisi et al., 2012).

These usages share a structural theme rather than a single governing formula. Each formulation isolates a reduced dynamics that mediates between diffusion-dominated and wave- or transport-dominated behavior. In the planar mixed-type problems, the reduction is geometric and interfacial: the PDE changes type across a prescribed curve. In the large-deviation setting, the reduction is variational and control-theoretic: the skeleton equation is the controlled hydrodynamic constraint entering the rate function. In the Kirchhoff setting, the reduction is asymptotic: the skeleton is the parabolic equation obtained by formally suppressing the inertial term.

2. Canonical mixed-type equation and geometric classification

The canonical mixed-type model is

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$

For $y>0$, the second-order part has $\Gamma=\{y=0\}$0, $\Gamma=\{y=0\}$1, $\Gamma=\{y=0\}$2, hence $\Gamma=\{y=0\}$3, so the equation is parabolic. For $\Gamma=\{y=0\}$4, the second-order part is $\Gamma=\{y=0\}$5, with $\Gamma=\{y=0\}$6, $\Gamma=\{y=0\}$7, $\Gamma=\{y=0\}$8, hence $\Gamma=\{y=0\}$9, so the equation is hyperbolic (Shishkina et al., 2023).

The domain $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$0 is a simply connected bounded subset of the $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$1-plane. In the upper part $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$2 it is bounded by the vertical segments on $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$3 and $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$4 and by the line $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$5. In the lower part $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$6 it is bounded by the characteristics

$$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$7

issuing from $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$8 and $$\partial_t \rho=\Delta\Phi(\rho)-\nabla\cdot\big(\Phi^{1/2}(\rho)\,g\big);$$9. The type-change curve is the horizontal line

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$0

The change of type is therefore imposed by switching the governing operator across $b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$1, not by a single operator whose discriminant vanishes on the interface (Shishkina et al., 2023).

The corresponding boundary data are also split by type. In the parabolic region,

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$2

In the hyperbolic region, the data are prescribed as a linear combination of characteristic traces:

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$3

with real constants $b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$4 satisfying $b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$5 (Shishkina et al., 2023). This choice is characteristic of mixed-type formulations in which the upper and lower subdomains require different boundary mechanisms.

A fractional Tricomi-type variant is posed on

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$6

with

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$7

and equations

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$8

$b(t)u'(t)+m(|A^{1/2}u(t)|^2)Au(t)=0$9

Here

$\mathbb T^d\times(0,T)$0

so the upper region is a time-fractional diffusion-type regime and the lower region remains hyperbolic (Karimov et al., 2013).

3. Interface traces, transmission, and gluing laws

For the canonical mixed problem, regularity imposes continuity across $\mathbb T^d\times(0,T)$1:

$\mathbb T^d\times(0,T)$2

where

$\mathbb T^d\times(0,T)$3

The parabolic equation above the interface implies the relation

$\mathbb T^d\times(0,T)$4

The hyperbolic region is then treated as a Cauchy problem with data

$\mathbb T^d\times(0,T)$5

so the interface is the locus where the two subproblems are glued through continuity and the PDE-implied identity $\mathbb T^d\times(0,T)$6 (Shishkina et al., 2023).

The boundary data must satisfy a compatibility condition obtained by evaluating the characteristic relation at the endpoints:

$\mathbb T^d\times(0,T)$7

In the theorem of the model problem, the data regularity is

$\mathbb T^d\times(0,T)$8

together with $\mathbb T^d\times(0,T)$9 (Shishkina et al., 2023).

The fractional problem replaces local gluing by a nonlocal integral condition on $y=0$0:

$y=0$1

With the trace notation

$y=0$2

the interface system becomes

$y=0$3

$y=0$4

The hyperbolic side yields

$y=0$5

This gluing condition generalizes the standard continuity/flux match and reduces to a local classical gluing when $y=0$6 (Karimov et al., 2013).

A recurring misconception is that transmission for mixed-type problems must be encoded by an independent jump or flux law. In the canonical constant-coefficient model, no separate jump or flux conditions are imposed beyond continuity across $y=0$7 and the PDE-implied relation $y=0$8 (Shishkina et al., 2023). In the fractional model, by contrast, the interface law is explicitly nonlocal and weighted by $y=0$9 on the parabolic side (Karimov et al., 2013).

4. Interfacial reduction, explicit representations, and well-posedness

For the canonical model, the hyperbolic Cauchy problem in $AB=\{(x,0):0<x<1\}$0 is solved by the d’Alembert-type formula

$AB=\{(x,0):0<x<1\}$1

Substituting this into the characteristic boundary condition and differentiating gives

$AB=\{(x,0):0<x<1\}$2

Combining this with $AB=\{(x,0):0<x<1\}$3 produces the linear ODE

$AB=\{(x,0):0<x<1\}$4

subject to

$AB=\{(x,0):0<x<1\}$5

Multiplication by the integrating factor

$AB=\{(x,0):0<x<1\}$6

reduces the problem to a first-order equation for $AB=\{(x,0):0<x<1\}$7, and this yields a unique $AB=\{(x,0):0<x<1\}$8 under the stated hypotheses (Shishkina et al., 2023).

Once $AB=\{(x,0):0<x<1\}$9 is determined, the solution is reconstructed explicitly in each subdomain. In $\varepsilon\to0$0, the d’Alembert formula applies. In $\varepsilon\to0$1, the parabolic component is represented by

$\varepsilon\to0$2

where the Dirichlet heat kernel on $\varepsilon\to0$3 is

$\varepsilon\to0$4

These formulas define a function

$\varepsilon\to0$5

satisfying the PDE, all boundary data, and the transmission conditions (Shishkina et al., 2023).

Uniqueness in the homogeneous case is based on the energy-type identity

$\varepsilon\to0$6

With $\varepsilon\to0$7, this forces $\varepsilon\to0$8, hence $\varepsilon\to0$9, then $$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$0 in both subdomains. The resulting theorem states that if

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$1

and the compatibility condition holds, then the problem has a unique regular solution in $$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$2. No Lopatinskii–Shapiro or smallness assumptions are invoked; existence and uniqueness are proved by explicit interfacial reduction to an ODE for $$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$3 combined with explicit solution operators in each subdomain (Shishkina et al., 2023).

The fractional Tricomi problem follows a parallel but more nonlocal route. The interface equations reduce the problem to a second-order ODE for $$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$4 and then to a Fredholm equation of the second kind,

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$5

whose unique solution is written via the resolvent kernel $$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$6:

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$7

The uniqueness theorem assumes

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$8

while the existence theorem requires

$$0=\begin{cases} u_{xx}-u_y, & y>0,\ u_{yy}-u_{xx}, & y<0. \end{cases}$$9

The solution in $y>0$0 is represented by a Green function for the fractional diffusion equation involving a Wright-type kernel

$y>0$1

and the hyperbolic region is reconstructed by the Cauchy formula (Karimov et al., 2013).

5. Energy-critical skeleton PDE in non-equilibrium large deviations

In the large-deviation theory of conservative interacting particle systems, the skeleton equation is the controlled hydrodynamic PDE

$y>0$2

with nonnegative increasing $y>0$3 and irregular conservative drift

$y>0$4

The entropy functional is

$y>0$5

and the initial data belong to

$y>0$6

For porous medium nonlinearities $y>0$7, the equation is energy-critical for $y>0$8 and $y>0$9; it is supercritical for $\Gamma=\{y=0\}$00 initial data with $\Gamma=\{y=0\}$01 (Fehrman et al., 2019).

The analytic framework is built on renormalized kinetic solutions. The kinetic function is

$\Gamma=\{y=0\}$02

and the parabolic defect measure is

$\Gamma=\{y=0\}$03

The kinetic equation takes the conservative form

$\Gamma=\{y=0\}$04

A renormalized kinetic solution is a nonnegative $\Gamma=\{y=0\}$05 such that

$\Gamma=\{y=0\}$06

and the renormalized kinetic identity holds against compactly supported test functions. Classical weak solutions are defined by the same optimal regularity together with the integrated identity obtained from

$\Gamma=\{y=0\}$07

Under the structural hypotheses denoted as_unique and as_equiv, renormalized kinetic solutions coincide with classical weak solutions (Fehrman et al., 2019).

The key a priori estimate is the entropy-energy inequality

$\Gamma=\{y=0\}$08

Under as_unique, solutions satisfy the $\Gamma=\{y=0\}$09 contraction

$\Gamma=\{y=0\}$10

Under as_unique and as_compact, existence holds for every $\Gamma=\{y=0\}$11 and $\Gamma=\{y=0\}$12, and every renormalized kinetic solution admits an $\Gamma=\{y=0\}$13-continuous representative in time. Stability is strong in $\Gamma=\{y=0\}$14 under weak convergence of controls and initial data, provided uniform entropy bounds hold (Fehrman et al., 2019).

This PDE enters the rate function

$\Gamma=\{y=0\}$15

and the well-posedness theory supports several applications: identification of the lower semicontinuous envelope of restricted rate functions for zero range processes, zero-noise large deviations for conservative singular SPDE, and $\Gamma=\{y=0\}$16-convergence of finite-mode rate functions. The paper states that the identification of the lower semicontinuous envelope solves a longstanding open problem in the large deviations for zero range processes (Fehrman et al., 2019).

6. Parabolic regime as skeleton limit for weakly dissipative hyperbolic equations

A distinct notion of parabolic-hyperbolic skeleton equation appears in the coercive Kirchhoff problem

$\Gamma=\{y=0\}$17

posed in a real Hilbert space $\Gamma=\{y=0\}$18 with a self-adjoint nonnegative coercive operator $\Gamma=\{y=0\}$19 satisfying

$\Gamma=\{y=0\}$20

and a nondegenerate locally Lipschitz Kirchhoff nonlinearity $\Gamma=\{y=0\}$21 with

$\Gamma=\{y=0\}$22

Formally setting $\Gamma=\{y=0\}$23 yields the parabolic limit

$\Gamma=\{y=0\}$24

which serves as the skeleton or leading-order dynamics (Ghisi et al., 2012).

The dissipation threshold is determined by

$\Gamma=\{y=0\}$25

For $\Gamma=\{y=0\}$26 this is the parabolic regime; for $\Gamma=\{y=0\}$27 it is the hyperbolic regime. In the coercive case, the parabolic problem has a unique global solution and satisfies

$\Gamma=\{y=0\}$28

For the hyperbolic equation with $\Gamma=\{y=0\}$29 and small $\Gamma=\{y=0\}$30, the functional

$\Gamma=\{y=0\}$31

obeys the optimal decay bounds

$\Gamma=\{y=0\}$32

$\Gamma=\{y=0\}$33

$\Gamma=\{y=0\}$34

The paper emphasizes a “quite surprising fact”: in the coercive case the analogy between parabolic equations and dissipative hyperbolic equations is weaker than in the noncoercive case (Ghisi et al., 2012).

The singular perturbation is sharpened by a boundary-layer corrector $\Gamma=\{y=0\}$35, defined by

$\Gamma=\{y=0\}$36

Writing

$\Gamma=\{y=0\}$37

the remainder satisfies decay-error estimates of order $\Gamma=\{y=0\}$38 with the same decay profiles as the hyperbolic problem. For example, if $\Gamma=\{y=0\}$39,

$\Gamma=\{y=0\}$40

$\Gamma=\{y=0\}$41

$\Gamma=\{y=0\}$42

The optimality theorem shows that no faster decay can hold for nontrivial hyperbolic solutions, so the parabolic limit is the skeleton only after the corrector and the sharp decay mismatch are taken into account (Ghisi et al., 2012).

Taken together, these formulations show that “parabolic-hyperbolic skeleton equation” is best understood as a structural designation. It may refer to a piecewise mixed-type PDE with explicit interfacial reconstruction, a fractional diffusion–wave system with nonlocal gluing, an energy-critical controlled hydrodynamic equation underpinning non-equilibrium large deviations, or the parabolic limit of a weakly dissipative hyperbolic evolution. The common mathematical content is the extraction of a reduced equation that organizes well-posedness, interface coupling, decay, or variational rate-function structure across a parabolic-hyperbolic divide.