Interfacial Jump Conditions

• Interfacial jump conditions are clear definitions of local discontinuities in field variables and their normal derivatives across sharply defined interfaces in multiphase systems.
• They ensure physical fidelity by encoding conservation of mass, momentum, and energy through precise flux and transmission conditions at material discontinuities or singular sources.
• Numerical methods such as immersed interface and ghost fluid methods leverage these conditions to achieve high-order accuracy and stability in simulations with sharp property transitions.

Interfacial Jump Conditions

Interfacial jump conditions rigorously specify the local discontinuities of field variables and their normal derivatives across a sharply defined interface—a surface separating media that admit distinct constitutive parameters or support singular sources. They arise naturally in models of multiphase flows, composite materials, turbulence closures, chemical reactions, and conjugate heat/mass transfer, and are essential to guarantee physical fidelity and numerical accuracy in the presence of sharp transitions. Jump conditions encode the conservation or exchange of measures such as mass, momentum, energy, or scalar concentrations at the differential level and dictate the coupling of bulk PDEs with interface laws. Precise formulation and numerical enforcement of these conditions are central to high-order and robust algorithms for interface-resolved PDEs.

1. Mathematical Formulation of Interfacial Jump Conditions

Interfacial jump conditions relate the limiting values of solution fields and their (co-)normal fluxes on either side of the interface. For a smooth codimension-one interface $\Gamma$ splitting a domain $\Omega$ into subdomains $\Omega^-$ and $\Omega^+$, the canonical expressions are: $$[u]_\Gamma \equiv u^+|_{\Gamma} - u^-|_{\Gamma} = g_D(x), \quad [\mathsf{F} \cdot n]_\Gamma \equiv \mathsf{F}^+ \cdot n - \mathsf{F}^- \cdot n = g_N(x)$$ where $u$ is the solution (scalar or vector), $\mathsf{F}$ the flux (often $\mathsf{F} = \mathcal{A} \nabla u$), $n$ the unit normal pointing from $\Omega^-$ to $\Omega^+$, and $g_D$, $g_N$ are prescribed jumps, possibly nonlinear ($g(x,u^+,u^-)$).

In elliptic, parabolic, or hyperbolic settings, jump conditions result from pillbox (distributional) integration of the PDE across $\Gamma$, enforcing the correct transmission or creation of the conserved quantity. In multi-physics contexts—e.g., two-phase turbulent flow (Elsayed et al., 2023), vapor-liquid phase change (Chen, 2022, Kalempa et al., 2024, Chen, 28 Jan 2025)—they encode Rankine-Hugoniot type balances, including molecular accommodation and kinetic-theory corrections.

For fluid-solid boundaries or walls, jump conditions typically reduce to velocity slip and temperature jump relations, given, for instance, by the Maxwell and Smoluchowski formulas in kinetic theory (Shu et al., 2017, Shu et al., 2016).

2. Physical Origin and Classification

Interfacial jumps arise through several mechanisms:

(i) Material Discontinuity: Contrasts in physical properties (density, viscosity, diffusivity, thermal conductivity) naturally induce flux discontinuities, often with solution continuity (Adjerid et al., 2020, Ji et al., 17 May 2025, Bochkov et al., 2019).

(ii) Singular Source/Sink on the Interface: Localized chemical reactions, surface tension, phase change, or deposition ablate or inject the conservative quantity only at $\Gamma$, directly generating $\delta$-function sources, e.g., $[\partial_n u] = g(x,u)$ (Bhaskar et al., 2013).

(iii) Rankine-Hugoniot or Momentum/Stress Balance: In multiphase flows, continuity of stress or mass-flux with possible interfacial force or transfer term, as in the Stokes equations with singular interfacial force (Gera et al., 2013), or RANS turbulence closures with discontinuous viscosity/density (Elsayed et al., 2023).

(iv) Molecular-level Kinetics: Kinetic-theory jump conditions (velocity/temperature/density) at fluid-solid and vapor-liquid interfaces, incorporating accommodation coefficients, Knudsen layers, and phonon mismatch (Shu et al., 2017, Shu et al., 2016, Kalempa et al., 2024).

(v) Porous Media and Reservoir Simulation: Sharp saturation or capillarity-induced "flood front" leading to jumps in phase velocities/pressures, with only bulk combination (e.g., total pressure/velocity) continuous (Peng et al., 2016, Peng et al., 2016).

A concise taxonomy appears in the following table:

Jump Type	Mathematical Form	Prototypical System
Dirichlet/solution jump	$[u] = g_D$	Composite media
Normal flux jump	$[A \partial_n u] = g_N$	Heat/solute, anisotropic elliptic PDE
Stress/momentum jump	$[n \cdot \sigma n] = f_\Gamma$	Multiphase Stokes/Navier-Stokes
Velocity/temperature jump	$[u] = b \partial_n u$	Kinetic theory (Maxwell/Smoluchowski)
Dissipation/turbulence	$[\omega] \propto [\nu]$	Two-phase RANS, turbulence across interface

3. Analytical and Computational Enforcement

The numerical imposition of interfacial jump conditions requires algorithmic strategies to preserve consistency, order, and stability irrespective of interface location with respect to the computational grid. Approaches include:

(a) Fitted Discretizations: Mesh aligns with $\Gamma$ (e.g., standard finite element methods), giving rise to exact subdomain coupling. Rarely practical for evolving/complex geometry.

(b) Embedded/Immersed Methods:

• Immersed Interface Method (IIM): Discrete correction terms constructed via Taylor expansion and local interface geometry to account for discontinuities (including up to the second derivatives for Stokes flows) (Gera et al., 2013).
• Ghost Fluid Method (GFM): Ghost values extrapolated across $\Gamma$ to ensure the prescribed jump, enabling standard stencils (Cho et al., 2020).
• Jump Splicing/Symmetric Correction: Finite-difference stencils are "spliced" with jump corrections based on normal Taylor expansion, preserving symmetry and high-order accuracy (Preskill et al., 2016).
• Immersed Finite Element (IFE): Local basis functions or enrichment constructed to locally satisfy interface jumps on cut elements. Recent developments prove unisolvence, mesh-independent conditioning, and optimal convergence for high-contrast, anisotropic, and nonhomogeneous jumps (Adjerid et al., 2020, Ji et al., 17 May 2025).
• Volume Penalization (VPM): Introduce a source term in a unified Cartesian discretization representing integrated effect of the local jump; the diffuse interface formulation ensures correct net jump imposition without spurious leakage (Liu et al., 15 Jan 2026).

(c) Correction-Function and Least-Squares Approaches: Explicitly construct a correction (singular) function matching the jump data and enforce the condition weakly via least squares over suitable manifolds when only an implicit (level-set) interface description is available (Marques et al., 2017).

(d) Model-Specific Treatments: Highly specialized closures and coupling, as in turbulence models where only certain composite quantities remain continuous (e.g., "inverse turbulence area" in RANS $k\!-\!\omega$ models, (Elsayed et al., 2023)).

(e) Variational and A Posteriori Analysis: Weak formulations and residual-based estimators account for Robin-type or nonlinear (flux) jump conditions, with stability and error control ensured via interface-adapted decompositions (Bhaskar et al., 2013, Lee, 2023).

4. Selected Examples Across Physical Models

A. Scalar Elliptic/Parabolic PDEs: Standard jump conditions are

$$[u] = g_D(x), \qquad [\mathcal{A} \nabla u \cdot n] = g_N(x).$$

Both IFE (Ji et al., 17 May 2025) and sharp finite-volume (Bochkov et al., 2019) methods demonstrate that combination of basis modification and local Taylor expansion yields $O(h^2)$ accuracy for solution and mesh-independent conditioning even for large coefficient jumps and arbitrary interface geometry.

B. Fluid-Solid Velocity Slip/Temperature Jump: Maxwell slip and Smoluchowski temperature jump conditions are fundamental for rarefied gas and microfluidic applications: $$u_{\rm slip} = b\,\partial_n u|_{n=0},\qquad \Delta T = b_T\,\partial_n T|_{n=0},$$ with $b = ({2-\alpha})/{\alpha}\,\lambda$, $b_T=({2-\beta})/{\beta}\,{2}/{(\gamma+1)}\,\lambda$; extensions to liquids involve molecular accommodation, surface diffusion, and phonon-mismatch induced jumps (Kapitza resistance) (Shu et al., 2017, Shu et al., 2016).

C. Multiphase Darcy Flow and Flood Fronts: In water-oil displacement, only the total Darcy velocity $v_t$ and (Antoncev-Chavent) total pressure $P_{\rm tot}$ are continuous. Individual phase velocities and pressures exhibit jumps corresponding to saturation and capillary pressure transitions (Peng et al., 2016, Peng et al., 2016).

D. Turbulence/Two-Equation RANS Models: Sharp interfaces in air-water systems induce jumps in dissipation rate $\omega$ and viscous dissipation $\epsilon$, strictly proportional to kinematic viscosity jump and velocities' gradient continuity (Elsayed et al., 2023).

E. Phase Change and Liquid-Vapor Interfaces: At a vapor-liquid interface, detailed kinetic-theory analysis yields a trio of coupled jump relations for mass flux, heat flux, and the discontinuities in temperature, pressure, and density, all involving accommodation coefficients and matching at the Knudsen layer boundary. Table of jump coefficients for various models and accommodation regimes confirm weak dependence on collision model (Chen, 2022, Kalempa et al., 2024, Chen, 28 Jan 2025).

5. Theoretical and Practical Implications

Accurate specification and imposition of interfacial jump conditions ensure:

• Physical Consistency: Conservation laws and interface physics (mass, energy, momentum, chemical reaction rate) are respected in the sharp interface limit.
• Numerical Accuracy: High-order convergence and stability in finite difference, finite volume, and finite element methods even with nonconforming or moving interfaces.
• Robustness Across Parameter Regimes: Methods based on local algebraic elimination or basis-enrichment maintain conditioning and convergence independently of property contrast, geometry, or interface orientation (Bochkov et al., 2019, Ji et al., 17 May 2025).
• Model Closedness/Well-Posedness: For multiphysics coupling (e.g., phase-change, chemical reaction, multi-field turbulence models), the number and structure of jump conditions is essential for system solvability.
• Experimental and Simulation Validation: Quantitative agreement with experiment (e.g., measured temperature/velocity jumps, turbulence decay, pressure at flood front) is only achievable when jump conditions reflect correct microphysics—e.g., proper capillarity corrections, molecular accommodation (Shu et al., 2017, Chen, 2022).

6. Extensions and Current Research Challenges

Research continues to expand both the rigorous derivation and efficient enforcement of jump conditions:

• Generalization to Nonlinear and Nonlocal Jumps: Inclusion of nonlinear interface kinetics (chemical reactions, strongly coupled thermomechanics), higher-order derivatives, or nonlocal effects, as emergent in advanced composite, reactive, or soft-matter systems (Bhaskar et al., 2013).
• Adaptive and Multiscale Methods: Interface-resolving adaptivity, a posteriori error control, and explicit linkages between continuum-scale PDEs and underlying kinetic or atomistic jump laws (Lee, 2023, Ji et al., 17 May 2025, Liu et al., 15 Jan 2026).
• Mobile/Evolving Interfaces and Topological Change: Robust representation and enforcement of jump conditions on time-dependent manifolds, including level-set and phase-field approaches, and handling of topologically complex or multi-valued interfaces (Cho et al., 2020, Marques et al., 2017).
• Unification of Continuum-Kinetic (Gas-Liquid-Solid) Interface Laws: Toward a consistent, parameterizable boundary condition model for all regimes, mapping molecular-scale processes (adsorption, desorption, phonon scattering) to continuum jump parameters (Shu et al., 2016, Shu et al., 2017).
• Physical Origin of "Anomalous" Jumps: Negative temperature jumps (Kapitza resistance sign change), control of phase-velocity jumps in capillary-limited/reservoir flows, and complex turbulent sink/source behavior are active areas of investigation (Cao et al., 2015, Elsayed et al., 2023).

Interfacial jump conditions remain both a fundamental theoretical construct and a practical computational necessity across the broad landscape of multiphase, multi-physics, and heterogeneous systems.