Open problems for the coupled interfacial evolution equations of inertialess stratified three-layer channel flows: fourth-order Majda-Pego stability, Feigenbaum universality, and large-zeta asymptotics

Papaefthymiou, Papageorgiou and Pavliotis (J. Fluid Mech. 734 (2013) 114-143) derived fully nonlinear coupled evolution equations for the two interfaces of an inertialess three-layer channel flow, including a symmetric system whose dissipation is purely fourth order (surface tension) and a kinematically modified coupled Kuramoto-Sivashinsky system. They left open a three-part program for these systems, which remains open as of 2026: (i) a rigorous extension of the Majda-Pego stable-viscosity-matrix theory from second-order to fourth-order dissipation matrices, so that the stability criterion their numerics support becomes a theorem; (ii) verification (or refutation) of Feigenbaum universal scaling for the reverse (period-halving) cascades that regularize chaotic dynamics into traveling waves as the kinematic parameter grows; and (iii) asymptotic and numerical analysis of the large-kinematic-parameter limit and of the general unequal-density parameter regimes.

Background

When three immiscible viscous liquids flow in layers through a channel, the two internal interfaces can become unstable even at vanishingly small Reynolds number, where the classical inertial route to instability is absent. In this inertialess (Stokes-flow) limit, long-wavelength interfacial instabilities survive: Yih-type viscosity-stratified instability (the long-wave instability of a sheared flow whose layers have different viscosities, identified by Yih), and, specific to three or more layers, kinetic instabilities arising from the interaction of the two interfaces rather than from inertia. Kliakhandler and Sivashinsky (1995) identified a kinetic alpha effect - an instability mechanism of exactly this coupling-driven type - in viscosity-stratified creeping three-layer Poiseuille flow (pressure-driven flow between parallel plates), and Kliakhandler (1999) showed that long interfacial waves in multilayer thin films are governed by coupled Kuramoto-Sivashinsky systems. The scalar Kuramoto-Sivashinsky equation, $u_t+uu_x+u_{xx}+u_{xxxx}=0$ , is the canonical long-wave thin-film model combining a destabilizing second-order term (negative effective diffusion, driving long-wave growth) with stabilizing fourth-order surface-tension damping of short waves; coupling several such equations through the shared fluid layers produces the coupled Kuramoto-Sivashinsky (cKS) systems studied here.

Papaefthymiou, Papageorgiou and Pavliotis (2013) used systematic long-wave asymptotics for three immiscible layers in an inclined channel driven by gravity and a pressure gradient to derive fully nonlinear coupled evolution equations for the two interfacial positions in the inertialess limit. Removing an artificial small parameter through a distinguished asymptotic detuning of the physical parameter ratios (viscosities, densities, undisturbed layer thicknesses, surface tensions) away from a symmetric state at which the unregularized flux Jacobian has an umbilic degeneracy (coincident real eigenvalues that can turn complex, i.e. the flux is of mixed hyperbolic-elliptic type), they obtained two canonical systems: a symmetric inertialess $2\times2$ system whose only dissipation is fourth order, entering through a surface-tension dissipation matrix $d_{ij}$ multiplying interfacial curvature terms, and a kinematically modified cKS system in which a kinematic parameter $\zeta$ (the magnitude of the parameter detuning) enhances the kinematic instability terms relative to the inertial ones.

The relevant mathematical precedent is the theory of Majda and Pego (1985), who studied strictly hyperbolic $2\times2$ systems of conservation laws regularized by second-order viscosity, $\mathbf{u}_t+A(\mathbf{u})\mathbf{u}_x=(D(\mathbf{u})\mathbf{u}_x)_x$ . When the flux Jacobian $A$ is symmetric and the viscosity matrix $D$ is positive definite, a constant state is stable precisely when $l_i(\bar{\mathbf{u}})D(\bar{\mathbf{u}})r_i(\bar{\mathbf{u}})>0$ , where $l_i$ and $r_i$ are the left and right eigenvectors of $A$ ; violation produces the counterintuitive Majda-Pego instability, in which added dissipation destabilizes a hyperbolic state. Papaefthymiou et al. (2013) found the same mechanism operating through their fourth-order surface-tension dissipation: their computations support the analogous condition on the fourth-order dissipation matrix, but the 1985 analysis is proved only for second-order viscosity, and they state that the higher-order analysis remains open to their knowledge.

Numerically, the 2013 paper discovered a novel regularization phenomenon: as the kinematic parameter $\zeta$ increases, chaotic interfacial dynamics are regularized into steady traveling waves through a reverse Feigenbaum cascade, i.e. a sequence of period-halving bifurcations. Feigenbaum universality is the statement that period-doubling cascades in dissipative systems approach universal limiting constants, independent of the specific system: the bifurcation-parameter ratios converge to $\delta\approx4.6692$ and the orbit-amplitude scaling factors to $\alpha\approx2.5029$ . Universality has been verified numerically for scalar dissipative PDEs in small domains, e.g. the Nikolaevskii equation, a scalar model PDE combining a long-wave instability with stabilizing diffusion whose chaotic dynamics near onset is a standard test case for such cascades (Tanaka 2005), but not for the coupled three-layer systems studied here. The only experiments bearing on three-layer inertialess instability, by Jiang, Helenbrook, Lin and Weinstein (2005) on aqueous gelatin systems, confirmed the linear instability and predate the nonlinear predictions. The authors note that the specific four-parameter detunings used in their illustrative computations may be difficult to access experimentally, although analogous dynamics were found for simpler depth-only detunings.

Problem Statement

The problem is the future-work program of Papaefthymiou, Papageorgiou and Pavliotis (2013) for the coupled interfacial evolution equations they derived for inertialess stratified three-layer channel flows - the symmetric inertialess $2\times2$ system whose dissipation is fourth order through the surface-tension matrix, and the kinematically modified coupled Kuramoto-Sivashinsky system with kinematic parameter $\zeta$ . The common objective is to complete the mathematical theory of these same systems; three linked questions constitute that completion:

(i) Rigorous fourth-order stability theory. Establish or refute a generalization of the Majda-Pego stable-viscosity-matrix theory to systems whose regularization is fourth order (dissipation matrices multiplying fourth spatial derivatives, as produced by surface tension), covering the derived symmetric inertialess system. In particular: prove, or disprove, that the Majda-Pego-type condition - the left-eigenvector/dissipation-matrix/right-eigenvector products being positive at a constant state, evaluated with the fourth-order dissipation matrix - is a valid stability criterion, as the numerical results of 2013 (and later computations) indicate but do not establish. The fourth-order Majda-Pego-type criterion is settled for the derived symmetric inertialess system by a proof or by a counterexample in that system.

(ii) Feigenbaum universality of the period-halving cascades. Determine, by sufficiently resolved long-time computations and/or analysis of the kinematically modified coupled Kuramoto-Sivashinsky system, whether the reverse Feigenbaum cascades (sequences of period-halving bifurcations through which chaotic dynamics are regularized into steady traveling waves as $\zeta$ increases) satisfy Feigenbaum universal scaling, i.e. whether the bifurcation-parameter ratios and amplitude scalings converge to the universal Feigenbaum constants as the cascades accumulate, as the 2013 numerical findings are 'consistent with' but did not establish.

(iii) Large- $\zeta$ asymptotics and unequal densities. Analyze the large- $\zeta$ limit of the kinematically modified system, in which the kinematic terms dominate the inertial ones and chaotic or time-periodic dynamics are regularized into steady traveling waves, by asymptotic and numerical methods (the nature of the limiting traveling waves, the mechanism and rate of regularization, and the structure of the alternating periodic/chaotic windows); and extend the analyses of the derived systems from the illustrative equal-density slices to the general unequal-density parameter regimes, which the authors state can be analysed analogously but did not pursue. Part (iii) comprises both the large- $\zeta$ limit analysis, covering the three named aspects, and the extension to the general unequal-density regimes.

References

E. S. Papaefthymiou, D. T. Papageorgiou, G. A. Pavliotis, 'Nonlinear interfacial dynamics in stratified multilayer channel flows', Journal of Fluid Mechanics 734 (2013) 114-143. DOI: 10.1017/jfm.2013.443, https://doi.org/10.1017/jfm.2013.443
A. Majda, R. L. Pego, 'Stable viscosity matrices for systems of conservation laws', Journal of Differential Equations 56 (1985) 229-262. DOI: 10.1016/0022-0396(85)90107-X, https://doi.org/10.1016/0022-0396(85)90107-X
I. Kliakhandler, G. Sivashinsky, 'Kinetic alpha effect in viscosity stratified creeping flows', Physics of Fluids 7 (1995) 1866-1871. DOI: 10.1063/1.868501, https://doi.org/10.1063/1.868501
I. L. Kliakhandler, 'Long interfacial waves in multilayer thin films and coupled Kuramoto-Sivashinsky equations', Journal of Fluid Mechanics 391 (1999) 45-65. DOI: 10.1017/S0022112099005297, https://doi.org/10.1017/S0022112099005297
W. Y. Jiang, B. T. Helenbrook, S. P. Lin, S. J. Weinstein, 'Low-Reynolds-number instabilities in three-layer flow down an inclined wall', Journal of Fluid Mechanics 539 (2005) 387-416. DOI: 10.1017/S0022112005005781, https://doi.org/10.1017/S0022112005005781
D. Tanaka, 'Bifurcation scenario to Nikolaevskii turbulence in small systems', Journal of the Physical Society of Japan 74 (2005) 2223-2225. DOI: 10.1143/JPSJ.74.2223, https://doi.org/10.1143/JPSJ.74.2223 (also arXiv:nlin/0504031, https://arxiv.org/abs/nlin/0504031)
E. S. Papaefthymiou, D. T. Papageorgiou, 'Vanishing viscosity limits of mixed hyperbolic-elliptic systems arising in multilayer channel flows', Nonlinearity 28 (2015) 1607-1631. DOI: 10.1088/0951-7715/28/6/1607, https://doi.org/10.1088/0951-7715/28/6/1607
E. S. Papaefthymiou, D. T. Papageorgiou, 'Nonlinear stability in three-layer channel flows', Journal of Fluid Mechanics 829 (2017) R2. DOI: 10.1017/jfm.2017.605, https://doi.org/10.1017/jfm.2017.605
S. N. Gomes, D. T. Papageorgiou, G. A. Pavliotis, 'Stabilizing non-trivial solutions of the generalized Kuramoto-Sivashinsky equation using feedback and optimal control', IMA Journal of Applied Mathematics 82 (2017) 158-194. DOI: 10.1093/imamat/hxw011, https://doi.org/10.1093/imamat/hxw011
M. Vécsei, M. Dietzel, S. Hardt, 'Interfacial instability of liquid films coating the walls of a parallel-plate channel and sheared by a gas flow', Microfluidics and Nanofluidics 22 (2018), article 91. DOI: 10.1007/s10404-018-2111-z, https://doi.org/10.1007/s10404-018-2111-z (preprint arXiv:1604.05975, https://arxiv.org/abs/1604.05975)
J. P. Alexander, D. T. Papageorgiou, 'Ordered and disordered dynamics in inertialess stratified three-layer shear flows', Physical Review Fluids 7 (2022) 014804. DOI: 10.1103/PhysRevFluids.7.014804, https://doi.org/10.1103/PhysRevFluids.7.014804
V. Park, 'Analytical Study of Inertialess Instability in Three-Layer Film Flow', Advances in Applied Mathematics and Mechanics (2026) 1-32. DOI: 10.4208/aamm.oa-2025-0022, https://doi.org/10.4208/aamm.oa-2025-0022
S. Chu, I. Kliakhandler, A. Kurganov, 'On the Gelfand Problem and Viscosity Matrices for Two-Dimensional Hyperbolic Systems of Conservation Laws' (2024). arXiv:2405.04214, https://arxiv.org/abs/2405.04214
E. S. Papaefthymiou, 'Systems of nonlinear PDEs arising in multilayer channel flows', PhD thesis, Imperial College London (2014). Available at https://spiral.imperial.ac.uk/handle/10044/1/24974 (also https://hdl.handle.net/10044/1/24974)