On a thin film model with insoluble surfactant

Gabriele Bruell; Rafael Granero-Belinch\'on

REVIEW 2 major objections 4 minor 45 references

On a thin film model with insoluble surfactant

T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read For medium-sized initial data, the full surfactant thin-film system has global weak solutions that decay exponentially to the flat state.

desk verdict First global existence result for the full surfactant system, but Theorem 1's L²_t H²_x regularity does not follow from the estimates and likely fails; referee it after a correction. read the letter →

arxiv 1908.06406 v2 pith:D2P4FYP5 submitted 2019-08-18 math.AP

Gabriele Bruell , Rafael Granero-Belinchón This is my paper · ORCID

classification math.AP MSC 35D3035B4035K5235K6576A20

keywords thinfilmequationsinsolublesurfactantquasilinearparabolicsystemsdegenerateglobalweaksolutionsexponentialdecayratesWienerspaces

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper treats the lubrication model (1) for a viscous film whose free surface carries an insoluble surfactant, with gravity, capillarity, van der Waals forces, and surface diffusion all present. It aims to prove that for initial data that are medium-sized in the Wiener norm — measured by the $\ell^1$ sum of Fourier coefficients, not by Sobolev norms — the coupled parabolic system has a global weak solution, and that this solution decays exponentially to the flat film with an explicitly computable rate. Earlier analytical results covered gravity-only or capillary-only special cases, while the present size condition also admits highly oscillatory initial data for which Sobolev smallness fails. The proof controls the Wiener energy of the perturbations $(f,\Theta)$ around the flat steady state, and the same a priori estimates give conditional uniqueness for slightly more regular weak solutions.

What carries the argument

The machinery is the scale of Wiener spaces $\dot A^s(\mathbb{T})$, the spaces of $2\pi$-periodic integrable functions with $\sum_k |k|^s |\hat u(k)|$ finite; $\dot A^0$ is a Banach algebra, so products can be estimated without losing derivatives. The paper measures the perturbation $(f,\Theta)=(h-h_\sharp,\Gamma-\Gamma_\sharp)$ through $E^0_0(f,\Theta)=\|f\|_{\dot A^0}+\|\Theta\|_{\dot A^0}$. On the interval where $E^0_0$ stays below $\min\{h_\sharp,\Gamma_\sharp\}$, the geometric series expansions for $1/(1+f/h_\sharp)$ and its reciprocal powers converge, the algebra and interpolation inequalities bound every nonlinear term by constants $\Lambda_i$ times the energy, and the positive coefficients $C_i$ absorb them, giving a differential inequality in which the energy can only decrease. This keeps Galerkin approximations bounded, compactness turns them into a global weak solution, and the Poincare-type inequality on the Wiener scale converts the derivative bound into the exponential decay.

What would settle it

Run the Galerkin scheme for the even initial data $h_0=1+\mu\sin(1000x)$, $\Gamma_0=1/2+\mu\cos(1000x)$ with $\mu<1/356$, $A=0$, $G=D=1$, and evaluate $\partial_x h(t,0)$, $\partial_x\Gamma(t,0)$, and, when $S>0$, $S\partial_x^3 f(t,0)$ at several positive times; any nonzero value would show the even-reflection bridge to the interval no-flux problem has failed, so the theorems would apply only on the torus.

Watch

Extended reading notes

Core claim

The central claim is that, for both the gravity-driven flow ($S=0$) and the capillary-driven flow ($S>0$), if the initial perturbation energy $E^0_0(f_0,\Theta_0)=\|f_0\|_{\dot A^0}+\|\Theta_0\|_{\dot A^0}$ is below $\min\{h_\sharp,\Gamma_\sharp\}$, where $h_\sharp$ and $\Gamma_\sharp$ are the spatial means (the flat steady state), and if the structural constants $C_1,C_2,C_3$ in (20)--(24) are positive and dominate the initial-energy-weighted constants $\Lambda_i$, then (10) has at least one global weak solution in the sense of Definition 1. The solution satisfies exponential decay $\|f(t)\|_{L^\infty}+\|\Theta(t)\|_{L^\infty}\le E^0_0(f_0,\Theta_0)e^{-\delta t}$, with $\delta=\min\{\gamma_1,\gamma_2,\gamma_3\}$ for $S>0$ and $\delta=\min\{\gamma_1,\gamma_2\}$ for $S=0$, where each $\gamma_i$ is an explicit positive number computed from the parameters and the initial data. Uniqueness holds whenever the weak solution has the additional regularity $L^1(0,T;\dot A^4)\times L^1(0,T;\dot A^2)$ in the capillary case, and $L^1(0,T;\dot A^2)^2$ in the gravity case.

Load-bearing premise

The load-bearing premise is that solving on the torus is equivalent to solving on the interval with no-flux boundary conditions: the paper asserts that even reflection of the initial data preserves $\partial_x h=\partial_x\Gamma=0$ (and $S\partial_x^3 f=0$) under the nonlinear evolution, but gives no proof of that preservation.

Editorial extensions

If this is right

For both $S=0$ and $S>0$, medium-sized initial data in $\dot A^0$ produce global weak solutions, so under the stated size conditions neither the film height nor the surfactant concentration can develop a finite-time singularity.
The positivity information $\|f(t)\|_{L^\infty}<h_\sharp$ and $\|\Theta(t)\|_{L^\infty}<\Gamma_\sharp$ is preserved, so $h$ and $\Gamma$ remain strictly positive for all times.
The decay rate $\delta$ is explicit and computable from $G,S,A,D$ and the means $h_\sharp,\Gamma_\sharp$, so the result gives a quantitative prediction for how quickly the film relaxes to flat.
The hypotheses admit highly oscillatory initial data such as $h_0=1+\mu\sin(1000x)$, $\Gamma_0=1/2+\mu\cos(1000x)$ with $\mu<1/356$ in the example $A=0,G=D=1$; this is medium-sized in the Wiener norm even though the $H^2$ norm is of order $10^3$.
Weak solutions with one additional integrability level, $L^1(0,T;\dot A^2)^2$ for $S=0$ and $L^1(0,T;\dot A^4)\times L^1(0,T;\dot A^2)$ for $S>0$, are unique.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the even-reflection step can be justified, the theorems would apply to the physical interval no-flux problem; until then the proved statements should be read as torus-periodic.
Because the argument uses only the Banach algebra and interpolation properties of $\dot A^0$, one could replace the linear surface-tension law $\sigma(s)=1-s$ by any constitutive relation admitting a convergent expansion around the steady state, with the constants (20)--(24) recomputed accordingly.
The explicit decay rate $\delta$ is a quantitative prediction that numerical experiments could test directly: a simulated $L^\infty$ distance to the flat state staying below $E^0_0(f_0,\Theta_0)e^{-\delta t}$ would support the energy estimate, while persistent oscillations would point to missing damping.
One can probe sharpness of the medium-size condition by increasing the initial energy until some $\gamma_i$ becomes nonpositive and checking numerically whether the flat state loses stability; an earlier instability would show the condition is sufficient but not necessary.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The genuinely new thing here is the first global weak existence and exponential decay result for the full thin-film surfactant system with all forces G, S, A, D positive. The Wiener-space energy method is borrowed from the authors' own Muskat work, but applying it to this mixed-order system is a real step forward. The a priori estimate leading to (30) and the decay inequality are worked out in detail and look correct. The smallness conditions are explicit and the paper honestly maps the prior special cases.

The main problem is the regularity stated in Theorem 1. For S = 0, the estimates give L∞_t A0 ∩ L¹_t A2 with time derivatives in L¹_t A0. That implies L²_t H¹ by interpolation with the L∞ bound, but it does not imply L²_t H². The scalar heat equation with Fourier coefficients ε k^{-3/2} satisfies exactly these bounds, and its H² norm squared diverges in time, so the theorem's claim is not just unproven—it is false for admissible data. Section 4.2 only claims L²_t H¹, so the theorem statement is internally inconsistent. For S > 0, the analogous f∈L²_t H² can be recovered from L∞_t A0 ∩ L¹_t A4 via interpolation, and Θ∈L²_t H¹ from L¹_t A2, so Theorem 2 may survive, but Section 5.2's statement that L¹(A4) ⊂ L²(H²) is not true as written.

The other soft spot is the even-reflection reduction in Section 3.1. It is asserted, not proved. However, it is a standard symmetry argument: the system is reflection-invariant, so even initial data produce even solutions, and the no-flux boundary conditions follow from odd derivatives vanishing. This should be easy to fix.

The conditional uniqueness proofs are only sketches, but they are conditional on extra regularity, so that seems acceptable.

Who is this for? Anyone working on degenerate parabolic systems or thin-film surfactant models. The method is useful and the full-force result matters. The paper deserves a serious referee, but not acceptance as is. A major revision that corrects Theorem 1's regularity—either lowering it to L²_t H¹ or adding assumptions that actually imply H²—would make the paper solid.

Referee Report

2 major / 4 minor

Summary. The paper studies a one-dimensional thin film equation with an insoluble surfactant layer, including gravitational, capillary, van der Waals, and surface-diffusion effects. The authors reformulate the interval problem with no-flux boundary conditions as a periodic problem on the torus, introduce a Wiener-algebra framework, and prove two global existence theorems: one for the gravity-driven case (S = 0) and one for the capillary-driven case (S > 0). The solutions are obtained via Galerkin approximations and compactness arguments, they satisfy explicit exponential decay toward the flat equilibrium, and conditional uniqueness is claimed under additional regularity. The paper emphasizes that the smallness condition on the initial data is explicitly computable and that the initial data may be highly oscillatory.

Significance. The paper is a solid contribution to the analytical theory of thin film equations with surfactant. Its main novelty is the treatment of the full system with all four physical effects simultaneously, combined with a low-regularity Wiener-space setting that allows medium-sized, highly oscillatory initial data. The energy estimates in Sections 4.1 and 5.1 are detailed and appear correct, and the decay rates are given by explicit constants depending only on the parameters and the initial data. The example in Remark 2 is instructive and demonstrates that the hypotheses are not vacuous. If the issues raised below are addressed, the paper would be a useful and citable contribution to the field.

major comments (2)

[Theorem 1 and §4.2] The regularity assertion (f,Θ) ∈ (L²(0,T;H²))² in Theorem 1 is not supported by the a priori estimates. The uniform bounds (35)–(36) give boundedness in (L∞(0,T;A0) ∩ L¹(0,T;A2))² for the solutions and (L¹(0,T;A0))² for the time derivatives. As the text itself notes in §4.2, these imply at most uniform boundedness in (L²(0,T;H¹))². In general, the estimate ∫₀ᵀ ||f||_{A2} dt ≤ C together with ||f||_{L∞(A0)} ≤ C does not imply ∫₀ᵀ ||f||_{H²}² dt < ∞. For example, with Fourier coefficients c_k = ε k^{-3/2}, which are admissible in A0, the Galerkin sequence for a linear heat equation satisfies the counterparts of (35)–(36), yet ∫₀ᵀ ||f_M||_{H²}² dt behaves like ε² log M as M→∞. Thus the claimed H² regularity in the S=0 case does not follow from the proof and fails for admissible initial data in general. The statement should be weakened to L²(0,T;H¹) unless an additional estimate is supplied.
[Section 3.1] The reduction from the interval problem with no-flux boundary conditions (8) to the periodic torus is asserted but not proven. The sentence 'the evenness of initial data is preserved' requires justification, because the global weak solutions are constructed by compactness and no uniqueness in the class of weak solutions is available to conclude that the limit of even Galerkin approximants is even. While this is likely fixable by observing that the Galerkin system is equivariant under the reflection x ↦ -x, so the approximants remain even, the equivalence is load-bearing for the physical interpretation of Theorems 1 and 2, and it should be either proved or stated as an assumption.

minor comments (4)

[Sections 4.4 and 5.4] The conditional uniqueness statements are only sketched. The key differential inequality for the difference of two solutions is asserted without derivation, and it is not made explicit which terms require the smallness conditions from the theorems and the additional regularity assumptions. Since uniqueness is a secondary result, this is less serious, but the proof should be expanded to a level that allows verification of the constants and the integrability of the Gronwall factor.
[Section 4.3] The exponential decay inequality is stated for all t ≥ 0, but the proof via lower semicontinuity from the Galerkin approximations yields the inequality for almost every t (or in a time-averaged sense). To claim pointwise validity for all t, one would need continuity in time of the L∞ norm, which is not established for the weak solutions. The statement should be qualified as holding for almost every t ≥ 0, or an additional argument should be provided.
[Theorem 1 and Theorem 2] The notation L^{2/r}(0,T;W^{r,∞}) (and L^{4/s}(0,T;W^{s,∞}) in Theorem 2) is typeset ambiguously in the manuscript. It would be clearer to write the exponent explicitly as 2/r or 4/s, respectively, so that the reader does not confuse it with an L² space with a weight r.
[Section 4.2] In the paragraph following (37), the sentence 'From the previous fact we can infer that actually (f,Θ) ∈ (L^{2/r}(0,T;W^{r,∞}))²' is abrupt. Since this is an important step in the regularity argument, a short justification via the interpolation inequality in Wiener spaces would improve readability.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: existence and decay follow from explicit a priori energy estimates; the H2-regularity gap and interval-to-torus reduction are correctness concerns, not circular reasoning.

full rationale

The derivation is self-contained. The main theorems are proved by Galerkin approximation: local existence via Picard–Lindelöf, uniform Wiener-space energy estimates (equations (30), (34), and (41)) with explicit constants C_i, Λ_i, and γ_i computed from the parameters and initial data, and passage to the limit via standard compactness arguments. The exponential decay is derived from the differential inequality d/dt E_0^0 ≤ −δ E_0^0, with δ = min{γ1,γ2,γ3} defined from the verified positivity conditions; it is not an assumed conclusion or a fitted output. The only citation to the authors' companion paper [10] concerns the choice of method (Wiener-space framework and an analogous compactness passage), while the actual compactness steps are displayed in (35)–(37) and rely on the standard result [43]; no uniqueness theorem, ansatz, or conclusion is imported from a self-citation. Two concerns raised by the text are real but are not circularity: the Section 3.1 assertion that even reflection preserves the no-flux boundary conditions is not proven, and Theorem 1's stated L2(0,T;H2) regularity for S=0 does not follow from the displayed L∞(A0) ∩ L1(A2) plus time-derivative bounds. These are correctness or rigor gaps, not reductions of the conclusion to the hypotheses, so they do not raise the circularity score.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

No free parameters are fitted; the constants C1, C2, C3, Λ1, Λ2, Λ3 and δ are explicitly computed from the physical parameters and initial means. All additional premises are either standard mathematical tools or model assumptions inherent to the thin film surfactant problem.

assumptions (6)

standard math Wiener algebra (A^s(T), ||·||_{A^s}) is a Banach algebra and satisfies the interpolation inequalities (4)-(6).
Used throughout Sections 4 and 5 to bound nonlinear terms in the a priori estimates.
domain assumption The lubrication approximation model (1) with linear equation of state σ(s)=1-s is the correct description of the thin film with insoluble surfactant.
The entire analysis is performed on this model; the paper states this simplification is common in applications and numerics.
domain assumption Periodic extension via even reflection preserves the no-flux boundary conditions (8).
Section 3.1, used to reduce the interval problem to the torus; not proven in detail.
domain assumption Initial data belong to A^0(T) and satisfy the positivity condition E^0_0(f0,Θ0) < min{h♯,Γ♯}.
A hypothesis of Theorems 1 and 2, ensuring h and Γ remain positive throughout the evolution.
standard math Picard-Lindelöf theorem ensures local existence of the Galerkin ODE system.
Section 4.2, used to construct Galerkin approximations.
standard math Simon's compactness theorem [43, Corollary 4] provides strong convergence of the Galerkin sequence.
Section 4.2, used to pass to the limit and obtain weak solutions.

how reviews work

0 comments

Cite this review

Pith. "Pith review of On a thin film model with insoluble surfactant." pith.science (2026). https://pith.science/paper/D2P4FYP5

@misc{pith2026190806406,
  author       = {Pith},
  title        = {Pith review of: On a thin film model with insoluble surfactant},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/D2P4FYP5}},
  note         = {Machine review of arXiv:1908.06406}
}

read the original abstract

This paper studies the existence and asymptotic behavior of global weak solutions for a thin film equation with insoluble surfactant under the influence of gravitational, capillary and van der Waals forces. We prove the existence of global weak solutions for \emph{medium sized} initial data in \emph{large function spaces}. Moreover, exponential decay towards the flat equilibrium state is established, where an estimate on the decay rate can be computed explicitly.

Figures

Figures reproduced from arXiv: 1908.06406 by the authors.

**Figure 1.** Scheme of a thin film flow with insoluble surfactant Here, ΩT := (0, T) × Ω denotes the time-space domain for the unknown functions h and Γ, with Ω ⊂ R being an open, bounded interval. The system (1) is supplemented with initial conditions h(0, x) = h0(x) and Γ(0, x) = Γ0(x) for all x ∈ Ω, where h0 and Γ0 are given functions and boundary conditions ∂xh = ∂xΓ = 0, S∂ 3 x f = 0 for all x ∈ ∂Ω. The appearing parameters… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

45 extracted references · 44 canonical work pages

[1]

J. W. Barrett, H. Garcke, and R. N¨ urnberg. Finite elemen t approximation of surfactant spreading on a thin ﬁlm. SIAM J. Numer. Anal. , 41(4):1427–1464, 2003

work page 2003
[2]

J. W. Barrett and R. N¨ urnberg. Convergence of a ﬁnite-el ement approximation of surfactant spreading on a thin ﬁlm in the presence of van der Waals forces. IMA J. Numer. Anal. , 24(2):323–363, 2004

work page 2004
[3]

Beretta, M

E. Beretta, M. Bertsch, and R. Dal Passo. Nonnegative sol utions of a fourth-order nonlinear degenerate parabolic equation. Archive for Rational Mechanics and Analysis , 129(2):175–200, 1995

work page 1995
[4]

Bernis and A

F. Bernis and A. Friedman. Higher order nonlinear degene rate parabolic equations. Journal of Diﬀerential Equations, 83(1):179–206, 1990

work page 1990
[5]

A. L. Bertozzi and M. C. Pugh. The lubrication approximat ion for thin viscous ﬁlms : regularity and long-time behavior of weak solutions. Comm. Pure Appl. Math. , 49(2):85–123, 1996

work page 1996
[6]

M. S. Borgas and J. B. Grotberg. Monolayer ﬂow on a thin ﬁlm . J. Fluid Mech. , 193:151–170, 1988

work page 1988
[7]

G. Bruell. Modeling and analysis of a two-phase thin ﬁlm m odel with insoluble surfactant. Nonlinear Anal. Real World Appl. , 27:124–145, 2016

work page 2016
[8]

G. Bruell. Weak solutions to a two-phase thin ﬁlm model wi th insoluble surfactant driven by capillary eﬀects. Journal of Evolution Equations , 17(4):1341–1379, 2017. 19

work page 2017

Show all 45 references

[9]

G. Bruell. Well-posedness and stability for a mixed orde r system arising in thin ﬁlm equations with surfactant. to appear in Mathematische Nachrichten (2019)

2019
[10]

Bruell and R

G. Bruell and R. Granero-Belinch´ on. On the thin ﬁlm Mus kat and the thin ﬁlm Stokes equations. Journal of Mathematical Fluid Mechanics , 21(2):33, 2019

2019
[11]

Burczak and R

J. Burczak and R. Granero-Belinch´ on. On a generalized doubly parabolic Keller–Segel system in one spatial dimension. Mathematical Models and Methods in Applied Sciences , 26(01):111–160, 2016

2016
[12]

Chugunova and R

M. Chugunova and R. M. Taranets. Nonnegative weak solut ions for a degenerate system modeling the spread- ing of surfactant on thin ﬁlms. Appl. Math. Res. Express. AMRX , (1):102–126, 2013

2013
[13]

Constantin, T

P. Constantin, T. Dupont, R. Goldstein, L. Kadanoﬀ, M. S helley, and S. Zhou. Droplet breakup in a model of the Hele-Shaw cell. Physical Review E , 47(6):4169, 1993

1993
[14]

Constantin, D

P. Constantin, D. Crdoba, F. Gancedo, L. Rodrguez-Piaz za and R. Strain. On the Muskat problem: global in time results in 2D and 3D. American Journal of Mathematics , 138(6):1455–1494, 2016

2016
[15]

Constantin, D

P. Constantin, D. Crdoba, F. Gancedo and R. Strain. On th e global existence for the Muskat problem. J. Eur. Math. Soc. , 15(1), 201–227, 2013

2013
[16]

B. D. Edmonstone, R. V. Craster, and O. K. Matar. Surfact ant-induced ﬁngering phenomena beyond the critical micelle concentration. J. Fluid Mech. , 564:105–138, 2006

2006
[17]

Escher, M

J. Escher, M. Hillairet, Ph. Lauren¸ cot, and Ch. Walker . Global weak solutions for a degenerate parabolic system modeling the spreading of insoluble surfactant. Indiana Univ. Math. J. , 60(6):1975–2019, 2011

1975
[18]

Escher, M

J. Escher, M. Hillairet, Ph. Lauren¸ cot, and Ch. Walker . Thin ﬁlm equations with soluble surfactant and gravity: modeling and stability of steady states. Math. Nachr. , 285(2-3):210–222, 2012

2012
[19]

Escher, M

J. Escher, M. Hillairet, Ph. Lauren¸ cot, and Ch. Walker . Weak solutions to a thin ﬁlm model with capillary eﬀects and insoluble surfactant. Nonlinearity, 25(9):2423–2441, 2012

2012
[20]

Escher, M

J. Escher, M. Hillairet, Ph. Lauren¸ cot, and Ch. Walker . Traveling waves for a thin ﬁlm with gravity and insoluble surfactant. SIAM J. Appl. Dyn. Syst. , 14(4):1991–2012, 2015

1991
[21]

Escher, Ph

J. Escher, Ph. Lauren¸ cot, and B.-V. Matioc. Existence and stability of weak solutions for a degenerate parabolic system modelling two-phase ﬂows in porous media. Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire, 28(4):583–598, 2011

2011
[22]

Escher and Ch

J. Escher and Ch. Lienstromberg. Travelling waves in di latant non-Newtonian thin ﬁlms. J. Diﬀerential Equations, 264(3):2113–2132, 2018

2018
[23]

Escher, A.-V

J. Escher, A.-V. Matioc, and B.-V. Matioc. Modelling an d analysis of the Muskat problem for thin ﬂuid layers. Journal of Mathematical Fluid Mechanics , 14(2):267–277, 2012

2012
[24]

Escher, A.-V

J. Escher, A.-V. Matioc, and B.-V. Matioc. Thin-ﬁlm app roximations of the two-phase Stokes problem. Nonlinear Analysis: Theory, Methods & Applications , 76:1–13, 2013

2013
[25]

Escher and B.-V

J. Escher and B.-V. Matioc. Existence and stability of s olutions for a strongly coupled system modelling thin ﬂuid ﬁlms. NoDEA: Nonlinear Diﬀerential Equations and Applications , pages 1–17, 2013

2013
[26]

Escher and B.-V

J. Escher and B.-V. Matioc. Non-negative global weak so lutions for a degenerated parabolic system approx- imating the two-phase Stokes problem. Journal of Diﬀerential Equations , 256(8):2659–2676, 2014

2014
[27]

Gancedo, E

F. Gancedo, E. Garcia-Juarez, N. Patel, and R. Strain. O n the Muskat problem with viscosity jump: Global in time results. Advances in Mathematics , 345:552–597, 2019

2019
[28]

Garcke and S

H. Garcke and S. Wieland. Surfactant spreading on thin v iscous ﬁlms: nonnegative solutions of a coupled degenerate system. SIAM J. Math. Anal. , 37(6):2025–2048, 2006

2025
[29]

D. P. Gaver and J. B. Grotberg. The dynamics of a localize d surfactant on a thin ﬁlm. J. Fluid Mech. , 214:127–148, 1990

1990
[30]

Granero-Belinch´ on and M

R. Granero-Belinch´ on and M. Magliocca. Global existe nce and decay to equilibrium for some crystal surface models. Discrete & Continuous Dynamical Systems-A , 39(4):2101–2131, 2019

2019
[31]

Granero-Belinch´ on and S

R. Granero-Belinch´ on and S. Scrobogna. On an asymptot ic model for free boundary darcy ﬂow in porous media. arXiv preprint arXiv:1810.11798, 2018

2018 arXiv
[32]

Greenspan

H. Greenspan. On the motion of a small viscous droplet th at wets a surface. Journal of Fluid Mechanics , 84(1):125–143, 1978

1978
[33]

J. B. Grotberg and O. E. Jensen. The spreading of heat or s oluble surfactant along a thin liquid ﬁlm. Physics of Fluids A , 5(58), 1993

1993
[34]

O. E. Jensen and J. B. Grotberg. Insoluble surfactant sp reading on a thin viscous ﬁlm: shock evolution and ﬁlm rupture. J. Fluid Mech. , 240:259–288, 1992

1992
[35]

Lauren¸ cot and B.-V

Ph. Lauren¸ cot and B.-V. Matioc. A gradient ﬂow approac h to a thin ﬁlm approximation of the Muskat problem. Calculus of Variations and Partial Diﬀerential Equations , 47(1-2):319–341, 2013

2013
[36]

Lauren¸ cot and B.-V

Ph. Lauren¸ cot and B.-V. Matioc. A thin ﬁlm approximati on of the Muskat problem with gravity and capillary forces. Journal of the Mathematical Society of Japan , 66(4):1043–1071, 2014

2014
[37]

Lauren¸ cot and B.-V

Ph. Lauren¸ cot and B.-V. Matioc. Finite speed of propagation and waiting time for a thin-ﬁlm Muskat problem. Proceedings of the Royal Society of Edinburgh Section A: Mat hematics, 147(4):813–830, 2017

2017
[38]

Lauren¸ cot and B.-V

Ph. Lauren¸ cot and B.-V. Matioc. Self-similarity in a t hin ﬁlm Muskat problem. SIAM Journal on Mathe- matical Analysis, 49(4):2790–2842, 2017. 20 G. BRUELL AND R. GRANERO-BELINCH ´ON

2017
[39]

B.-V. Matioc. Non-negative global weak solutions for a degenerate parabolic system modelling thin ﬁlms driven by capillarity. Proceedings of the Royal Society of Edinburgh Section A: Mat hematics, 142(5):1071– 1085, 2012

2012
[40]

Pernas-Casta˜ no and J

T. Pernas-Casta˜ no and J. Vel´ azquez. Analysis of a thin ﬁlm approximation for two-ﬂuid Taylor-Couette ﬂows. arXiv preprint arXiv:1905.13606, 2019

1905 arXiv
[41]

M. Renardy. On an equation describing the spreading of s urfactants on thin ﬁlms. Nonlinear Anal. , 26(7):1207–1219, 1996

1996
[42]

Sheludko

A. Sheludko. Thin liquid ﬁlms. Advances in Colloid and Interface Science , 1(4):391–464, 1967

1967
[43]

J. Simon. Compact sets in the space Lp(0, T; B). Ann. Mat. Pura Appl. (4) , 146:65–96, 1987

1987
[44]

M. R. E. Warner, R. V. Craster, and O. K. Matar. Fingering phenomena associated with insoluble surfactant spreading on thin liquid ﬁlms. J. Fluid Mech. , 510:169–200, 2004

2004
[45]

S. G. Yiantsios and B. G. Higgins. A mechanism of Marango ni instability in evaporating thin liquid ﬁlms due to soluble surfactant. Physics of Fluids , 22(022102), 2010. Institute for Analysis, Karlsruher Institute of Technolog y (KIT), D-76128 Karlsruhe, Germany E-mail addres...

2010

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

J. W. Barrett, H. Garcke, and R. N¨ urnberg. Finite elemen t approximation of surfactant spreading on a thin ﬁlm. SIAM J. Numer. Anal. , 41(4):1427–1464, 2003

work page 2003

[2] [2]

J. W. Barrett and R. N¨ urnberg. Convergence of a ﬁnite-el ement approximation of surfactant spreading on a thin ﬁlm in the presence of van der Waals forces. IMA J. Numer. Anal. , 24(2):323–363, 2004

work page 2004

[3] [3]

Beretta, M

E. Beretta, M. Bertsch, and R. Dal Passo. Nonnegative sol utions of a fourth-order nonlinear degenerate parabolic equation. Archive for Rational Mechanics and Analysis , 129(2):175–200, 1995

work page 1995

[4] [4]

Bernis and A

F. Bernis and A. Friedman. Higher order nonlinear degene rate parabolic equations. Journal of Diﬀerential Equations, 83(1):179–206, 1990

work page 1990

[5] [5]

A. L. Bertozzi and M. C. Pugh. The lubrication approximat ion for thin viscous ﬁlms : regularity and long-time behavior of weak solutions. Comm. Pure Appl. Math. , 49(2):85–123, 1996

work page 1996

[6] [6]

M. S. Borgas and J. B. Grotberg. Monolayer ﬂow on a thin ﬁlm . J. Fluid Mech. , 193:151–170, 1988

work page 1988

[7] [7]

G. Bruell. Modeling and analysis of a two-phase thin ﬁlm m odel with insoluble surfactant. Nonlinear Anal. Real World Appl. , 27:124–145, 2016

work page 2016

[8] [8]

G. Bruell. Weak solutions to a two-phase thin ﬁlm model wi th insoluble surfactant driven by capillary eﬀects. Journal of Evolution Equations , 17(4):1341–1379, 2017. 19

work page 2017