A hyperbolic model for two-layer thin film flow with a perfectly soluble anti-surfactant
Pith reviewed 2026-05-23 02:27 UTC · model grok-4.3
The pith
Neglecting capillarity and diffusion yields a strictly hyperbolic system for two-layer anti-surfactant thin films that admits entropy pairs and explicit Riemann solutions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Under the assumptions of negligible capillarity and diffusion and perfect solubility of the solute, the lubrication-derived equations reduce to a conservative first-order system in film heights and concentration gradients. This reduced system is strictly hyperbolic for a certain set of states and admits an entire class of entropy/entropy-flux pairs, including a strictly convex entropy, which implies well-posedness of the Cauchy problem. The system is almost Temple-class, enabling explicit computation of solutions to the Riemann problem.
What carries the argument
The conservative first-order hyperbolic system in film heights and concentration gradients obtained after the lubrication limit and the stated reductions.
If this is right
- The Cauchy problem for the system is well-posed.
- Explicit solutions to the Riemann problem can be computed from the near-Temple-class structure.
- A Godunov-type finite volume scheme can be constructed using the exact Riemann solver.
- The model admits multiple entropy/entropy-flux pairs.
Where Pith is reading between the lines
- The same reduction steps could produce hyperbolic models for other thin-film problems that involve surface-active agents.
- The explicit Riemann solutions offer a benchmark for testing numerical methods on related multi-layer conservation laws.
- The identified hyperbolicity region may help select parameter ranges for laboratory validation of the reduced model.
Load-bearing premise
Capillarity and diffusion effects must be negligible and the solute must be perfectly soluble to reach the conservative first-order hyperbolic system.
What would settle it
A direct numerical simulation of the full lubrication equations or a physical experiment that produces flow behavior differing measurably from the hyperbolic model's predictions when capillarity remains small but nonzero.
Figures
read the original abstract
We consider the motion of a two-layer thin film that consists of two immiscible viscous fluids and is endowed with an anti-surfactant solute. The presence of such solute particles induces variations of the surface tension and interfacial stress driving a Marangoni-type flow. We first analyze a lubrication limit and derive one-dimensional evolution equations for film heights and solute concentrations. Then, under the assumption that the capillarity and diffusion effects are negligible and the solute is perfectly soluble, we obtain a conservative first-order system in terms of film heights and concentration gradients. This reduced system is found to be strictly hyperbolic for a certain set of states and to admit an entire class of entropy/entropy-flux pairs. We also provide a strictly convex entropy for the hyperbolic system. Thus, the well-posedness for the Cauchy problem is given. Moreover, the system is almost a Temple-class system, which allows to compute explicit solutions of the Riemann problem. The paper concludes with numerical experiments using a Godunov-type finite volume method, which relies on the exact Riemann solver.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives one-dimensional lubrication equations for the heights and solute concentrations of a two-layer thin film with an anti-surfactant. Under the explicit assumptions that capillarity and diffusion are negligible and the solute is perfectly soluble, the system reduces to a conservative first-order hyperbolic system in film heights and concentration gradients. The reduced system is shown to be strictly hyperbolic on a nonempty open set of states, to admit an entire family of entropy/entropy-flux pairs together with a strictly convex entropy (hence well-posedness of the Cauchy problem), and to be almost Temple-class, permitting explicit Riemann solutions. The paper closes with Godunov-type finite-volume numerics that employ the exact Riemann solver.
Significance. If the reduction is justified, the work supplies a new, mathematically structured model for Marangoni-driven two-layer films whose almost-Temple property yields explicit Riemann solutions and whose strictly convex entropy directly implies local well-posedness. These analytic features are genuine strengths that distinguish the contribution from purely numerical thin-film studies.
major comments (2)
- [§3 (reduction step)] The passage from the lubrication equations to the claimed conservative first-order system (abstract and §3) rests on dropping capillarity, diffusion, and solubility-related terms without quantitative error estimates or a priori bounds on the neglected contributions. Because this reduction is the load-bearing step that produces the hyperbolic system whose hyperbolicity, entropy, and Temple properties are then proved, the absence of such estimates leaves the domain of validity of the subsequent analysis unclear.
- [§4 (hyperbolicity)] The statement that the system is “strictly hyperbolic for a certain set of states” (abstract) is central to all later claims. The manuscript must explicitly identify this open set in the (h1,h2,c) variables, compute the eigenvalues of the flux Jacobian, and verify that they remain real and distinct throughout the set; without these explicit expressions the hyperbolicity claim cannot be checked.
minor comments (2)
- [§2–3] Notation for the two film heights and the concentration gradient should be introduced once and used consistently; several passages mix dimensional and nondimensional variables without explicit warning.
- [§6] The numerical section would benefit from a brief statement of the CFL condition employed with the exact Riemann solver and from a short convergence table (L1 error versus mesh size) for at least one Riemann problem whose exact solution is known.
Simulated Author's Rebuttal
Thank you for the referee's careful reading and constructive comments on our manuscript. We address each major comment point by point below, providing the strongest honest defense of the work while noting where revisions will strengthen the presentation.
read point-by-point responses
-
Referee: [§3 (reduction step)] The passage from the lubrication equations to the claimed conservative first-order system (abstract and §3) rests on dropping capillarity, diffusion, and solubility-related terms without quantitative error estimates or a priori bounds on the neglected contributions. Because this reduction is the load-bearing step that produces the hyperbolic system whose hyperbolicity, entropy, and Temple properties are then proved, the absence of such estimates leaves the domain of validity of the subsequent analysis unclear.
Authors: The reduction is a standard formal asymptotic step in lubrication theory, where capillarity, diffusion, and related terms are neglected to obtain the leading-order hyperbolic model under the assumptions of thin films and dominant Marangoni effects. The manuscript's primary contribution lies in analyzing the resulting conservative system rather than deriving rigorous approximation bounds. We will add a clarifying paragraph in §3 on the physical regime of validity (small capillary number, negligible diffusion) to better delineate applicability, but quantitative error estimates lie outside the paper's scope. revision: partial
-
Referee: [§4 (hyperbolicity)] The statement that the system is “strictly hyperbolic for a certain set of states” (abstract) is central to all later claims. The manuscript must explicitly identify this open set in the (h1,h2,c) variables, compute the eigenvalues of the flux Jacobian, and verify that they remain real and distinct throughout the set; without these explicit expressions the hyperbolicity claim cannot be checked.
Authors: We agree that explicit verification strengthens the claim. In the revised manuscript we will define the precise open set in (h1, h2, c) variables on which the system is strictly hyperbolic, compute the eigenvalues of the flux Jacobian explicitly, and state the conditions under which they are real and distinct. These details will be inserted into §4. revision: yes
Circularity Check
No significant circularity; derivation is self-contained
full rationale
The paper begins from standard lubrication theory for two-layer films, applies explicit physical assumptions (negligible capillarity/diffusion, perfect solubility) to drop higher-order terms, and obtains a first-order conservative system. Hyperbolicity, existence of entropy pairs, strict convexity of an entropy, and almost-Temple structure are then proved directly from the resulting PDE system using standard techniques for hyperbolic conservation laws. No fitted parameters, self-definitional relations, or load-bearing self-citations appear in the derivation chain. The central claims rest on independent mathematical analysis of the reduced equations rather than on any input being renamed as output.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Lubrication approximation is valid for the two-layer thin film geometry
- ad hoc to paper Capillarity and diffusion effects are negligible and the solute is perfectly soluble
Forward citations
Cited by 1 Pith paper
-
Global Existence for a Class of Keyfitz--Kranzer Systems with Application to Thin-Film Flows
Global weak entropy solutions exist for a class of non-symmetric Keyfitz-Kranzer systems that includes thin-film flow models.
Reference graph
Works this paper leans on
-
[1]
Anamika, R. Barthwal, and T. Raja Sekhar. Construction of solutions of the Riemann problem for a two-dimensional Keyfitz-Kranzer type model governing a thin film flow. Accepted for publication in Applied Mathematics and Computation
-
[2]
B. Andreianov, C. Donadello, S. S. Ghoshal, and U. Razafison. On the attainable set for a class of triangular systems of conservation laws. Journal of Evolution Equations , 15:503–532, 2015
work page 2015
-
[3]
R. Barthwal and T. Raja Sekhar. Two-dimensional non-self-similar Riemann solutions for a thin film model of a perfectly soluble anti-surfactant solution. Quarterly of Applied Mathematics, 80(4):717–738, 2022. 27
work page 2022
-
[4]
R. Barthwal, T. Raja Sekhar, and G. Raja Sekhar. Construction of solutions of a two- dimensional Riemann problem for a thin film model of a perfectly soluble antisurfactant solution. Mathematical Methods in the Applied Sciences , 46(6):7413–7434, 2023
work page 2023
-
[5]
A. Bertozzi, A. M¨ unch, and M. Shearer. Undercompressive shocks in thin film flows. Physica D: Nonlinear Phenomena , 134(4):431–464, 1999
work page 1999
-
[6]
J. Billingham. Surface tension-driven flow in a slender wedge. SIAM Journal on Applied Mathematics, 66(6):1949–1977, 2006
work page 1949
-
[7]
G. Br¨ ull. Modeling and analysis of a two-phase thin film model with insoluble surfactant. Nonlinear Analysis: Real World Applications , 27:124–145, 2016
work page 2016
-
[8]
C. Chalons, P. Engel, and C. Rohde. A conservative and convergent scheme for under- compressive shock waves. SIAM Journal on Numerical Analysis , 52(1):554–579, 2014
work page 2014
-
[9]
J. Conn, B. Duffy, D. Pritchard, S. Wilson, and K. Sefiane. Simple waves and shocks in a thin film of a perfectly soluble anti-surfactant solution. Journal of Engineering Mathematics, 107(1):167–178, 2017
work page 2017
-
[10]
J. Conn, B. R. Duffy, D. Pritchard, S. K. Wilson, P. J. Halling, and K. Sefiane. Fluid- dynamical model for antisurfactants. Physical Review E, 93(4):043121, 2016
work page 2016
-
[11]
J. J. Conn. Stability and dynamics of anti-surfactant solutions . PhD thesis, University of Strathclyde, 2017
work page 2017
-
[12]
B. Cook, A. L. Bertozzi, and A. Hosoi. Shock solutions for particle-laden thin films. SIAM Journal on Applied Mathematics , 68(3):760–783, 2008
work page 2008
-
[13]
R. V. Craster and O. K. Matar. Dynamics and stability of thin liquid films. Reviews of Modern Physics, 81(3):1131–1198, 2009
work page 2009
-
[14]
C. M. Dafermos. Hyperbolic conservation laws in continuum physics , volume 325 of Grundlehren der mathematischen Wissenschaften . Springer-Verlag, Berlin, 2000
work page 2000
-
[15]
D. P. Gaver and J. B. Grotberg. The dynamics of a localized surfactant on a thin film. Journal of Fluid Mechanics , 213:127–148, 1990
work page 1990
-
[16]
O. Jensen and J. Grotberg. Insoluble surfactant spreading on a thin viscous film: shock evolution and film rupture. Journal of Fluid Mechanics , 240:259–288, 1992
work page 1992
-
[17]
A. Kalogirou and M. Blyth. The role of soluble surfactants in the linear stability of two-layer flow in a channel. Journal of Fluid Mechanics , 873:18–48, 2019
work page 2019
-
[18]
T. Kato. The Cauchy problem for quasi-linear symmetric hyperbolic systems. Archive for Rational Mechanics and Analysis , 58(3):181–205, 1975
work page 1975
-
[19]
J. B. Keller and M. J. Miksis. Surface tension driven flows. SIAM Journal on Applied Mathematics, 43(2):268–277, 1983. 28
work page 1983
- [20]
-
[21]
A. Majda. Compressible fluid flow and systems of conservation laws in several space variables, volume 53. Springer Science & Business Media, 2012
work page 2012
-
[22]
O. Matar and S. Kumar. Rupture of a surfactant-covered thin liquid film on a flexible wall. SIAM Journal on Applied Mathematics , 64(6):2144–2166, 2004
work page 2004
-
[23]
Minhajul, T. Raja Sekhar, and G. P. Raja Sekhar. Stability of solutions to the Riemann problem for a thin film model of a perfectly soluble anti-surfactant solution. Commu- nications on Pure & Applied Analysis , 18(6):3367–3386, 2019
work page 2019
-
[24]
T. Myers. Thin films with high surface tension. SIAM Review, 40(3):441–462, 1998
work page 1998
-
[25]
S. B. G. O’Brien and L. W. Schwartz. Theory and modeling of thin film flows. Ency- clopedia of surface and colloid science , 1:5283–5297, 2002
work page 2002
-
[26]
J. Picardo, T. Radhakrishna, A. B. Vir, S. Ramji, and S. Pushpavanam. Modelling extraction in microchannels with stratified flow: channel geometry, flow configuration and Marangoni stresses. Indian Chemical Engineer, 57(3-4):322–358, 2015
work page 2015
-
[27]
J. R. Picardo, T. Radhakrishna, and S. Pushpavanam. Solutal Marangoni instability in layered two-phase flows. Journal of Fluid Mechanics , 793:280–315, 2016
work page 2016
- [28]
-
[29]
D. Serre. Systems of conservation laws 1: hyperbolicity, entropies, shock waves . Cam- bridge University Press, 1999
work page 1999
-
[30]
B. Temple. Systems of conservation laws with invariant submanifolds. Transactions of the American Mathematical Society, 280(2):781–795, 1983. 29
work page 1983
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.