REVIEW 2 major objections 4 minor 32 references
Bounded solutions and their asymptotics for a doubly nonlinear Cahn-Hilliard system
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper proves that a doubly nonlinear Cahn–Hilliard system has a unique bounded solution even when the bulk potential is singular and the viscosity nonlinearity is arbitrary, and that both regularization limits converge.
desk verdict The existence theorem is a real advance and the maximum-principle argument is the genuine novelty; the ε→0 limit has one explicitly formal, load-bearing estimate in Section 4.3 that should be fixed before that theorem is relied on. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument is carried by a time-discretization maximum principle. The regularized system is rewritten as a doubly nonlinear inclusion of the form $\partial\Xi((v^k-v^{k-1})/\tau)+A(v^k)\ni f^k$, and each discrete step is tested against the positive part $(v^k-b'_0)_+$ and $(a'_0-v^k)_+$. The sign of the potential term is controlled by the fast growth of $\psi'$ near the edges of $(a,b)$ and by properties of the resolvent $(I+\lambda\gamma)^{-1}$, forcing every discrete value to remain between $a'_0$ and $b'_0$; passing to the limit in the discretization and in the standard smoothings of the monotone graphs gives the uniform bound on $u$ that makes the singular nonlinearities well defined. The same boundedness is then used as an a priori estimate to let $\varepsilon$ or $\delta$ tend to zero, with compactness results and monotone graph closure arguments identifying the limit selection $\xi\in\beta(\partial_t u)$.
What would settle it
Carry out the omitted smoothing-regularization computation behind Section 4.3. If the calculation produces boundary terms on $\partial\Omega$ or terms involving $\beta'_\lambda(\partial_t u_\lambda)$ that cannot be bounded independently of the smoothing parameter $\lambda$, then the bound (4.5)–(4.6) is not available and Theorem 2.4's $L^\infty$ conclusion would need weaker alternatives. A concrete numerical test of the same phenomenon: take $\beta$ as the sign graph and $\psi$ with $|\psi''(r)|\le M(1+|r|^5)$, and check whether $\sup_{\varepsilon>0}\|u_\varepsilon\|_{L^\infty(Q)}$ is finite as $\varepsilon\to 0$ for an initial datum obeying (2.20)–(2.22).
Extended reading notes
Core claim
On its own terms, the paper establishes three statements. First, under the hypotheses (2.2)–(2.10), Theorem 2.2 gives a unique triplet $(u,\mu,\xi)$ with $u\in W^{1,\infty}(0,T;H)\cap H^1(0,T;V)\cap L^\infty(0,T;W_n)$, with $a'_0\le u\le b'_0$ almost everywhere for a compact subinterval $[a'_0,b'_0]\subset(a,b)$, with $\xi\in\beta(\partial_t u)$ almost everywhere, and with (2.16)–(2.18) holding. This removes the growth restrictions on $\beta$ and $\psi$ that earlier existence proofs required. Second, Theorem 2.4 says that with $\delta>0$ fixed and data satisfying (2.20)–(2.23) or (2.36), the family of solutions computed with $\varepsilon>0$ has a subsequence converging to a solution of the $\varepsilon=0$ limit problem, including the selection $\xi\in\beta(\partial_t u)$ in the limit. Third, Theorem 2.6 says that with $\varepsilon>0$ fixed and data satisfying (2.37)–(2.42), the solutions converge as $\delta\downarrow 0$ to the $\delta=0$ limit problem, and the difference obeys the rate estimate (2.50): $\|\mu_\delta-\mu\|_{L^2(0,T;V_0)}+\|u_\delta-u\|_{H^1(0,T;H)}\le M(\delta^{1/4}+\|u_{0\delta}-u_0\|_H+\|g_\delta-g\|_{L^2(0,T;H)})$.
Load-bearing premise
The load-bearing step is the third a priori estimate in Section 4.3, which the authors derive by formally testing the chemical-potential inclusion against $-\delta\Delta u_t+\partial_t\gamma(u)$ and describe as formal, deferring the rigorous justification to a smoothing regularization of $\beta$ that is only sketched in words; if that estimate cannot be made rigorous, the uniform $L^\infty$ bounds on $u$ and $\xi$ that the $\varepsilon\to 0$ theorem needs would collapse.
Editorial extensions
If this is right
- For fixed $\varepsilon,\delta>0$, the system is well posed for physically natural singular potentials such as the logarithmic double-well potential, with the solution confined to a compact subinterval of the potential's domain.
- With $\delta$ fixed, sending $\varepsilon$ to zero shows that the diffusive term alone can handle a noncoercive, multivalued viscosity law, so the vanishing-viscosity limit has a meaningful solution.
- With $\varepsilon$ fixed, sending $\delta$ to zero converges the whole family to the hysteresis-type limit equation, and the error is controlled by $\delta^{1/4}$ plus the approximation errors in the initial data and forcing term.
- In both limits, $\psi'(u)$ and the selection $\xi$ stay in bounded spaces, so the nonlinearities do not leave the framework even though they were singular or multivalued.
Reading between the lines
- The same positive-part test on time-discretization steps should extend to other doubly nonlinear gradient-flow systems with two competing regularizations, for instance with nonconstant mobility or different boundary conditions, because the sign argument only uses monotonicity and the steepness of $\psi$ outside an interval.
- The rate in (2.50) likely has room to move: the proof yields $\delta^{1/4}$ from balancing $\sqrt\delta$ and $\delta^{3/4}$ estimates, and a sharper comparison of the limit problem might show whether $\delta^{1/2}$ or a similar improved rate is attainable.
- If the formal Section 4.3 estimate were to fail, the alternative assumption (2.36) in Theorem 2.4 suggests a natural fallback: retain well-posedness but replace the $L^\infty$ regularity assertions with the corresponding $L^2$ statements rather than abandoning the $\varepsilon\to0$ limit.
- The boundedness argument might be adapted to prove existence for logarithmic potentials in related phase-separation models with hysteresis, since the only point where the potential enters is via $\psi'$ exceeding a fixed threshold outside the interval $[a'_0,b'_0]$.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the doubly nonlinear Cahn–Hilliard system (1.1)–(1.4), in which the chemical potential contains both a viscous term ε∂tu and a diffusive term −δΔu, with a maximal monotone nonlinearity β acting on ∂tu and a possibly singular potential ψ. The main results are: Theorem 2.2, existence and uniqueness of a bounded solution with values in a compact subinterval of the effective domain of ψ, without polynomial growth restrictions on β or ψ; Theorem 2.3, continuous dependence; Theorem 2.4, convergence as ε → 0 with δ > 0 fixed; and Theorem 2.6, convergence as δ → 0 with ε > 0 fixed, including the rate estimate (2.50). The existence proof combines a Yosida-type approximation with a nonstandard maximum-principle argument for time-discretized problems, and the asymptotic proofs use energy estimates, compactness, and monotonicity. The paper explicitly labels the third estimate of Section 4.3 as formal and defers its rigorous justification to a Yosida regularization that is described only in words.
Significance. If the formal estimates are supplied, the results are a genuine advance: the no-growth existence and uniqueness theorem extends previous work in [5] to singular ψ and β, and the two asymptotic limits answer natural regularization questions for this system. The maximum-principle argument for time-discretizations is nonstandard and is the main technical novelty of the existence proof. The δ → 0 rate estimate (2.50) is a concrete quantitative result that goes beyond mere compactness. The paper is also transparent about the formal step, which makes the gap identifiable rather than hidden. However, the proof of Theorem 2.4 is incomplete as written because the uniform bounds (4.5)–(4.6) rest on an unjustified formal computation; a parallel gap appears in the δ → 0 proof in Section 5.3. These gaps are fillable, and no internal inconsistency or circularity is apparent.
major comments (2)
- [Section 4.3, Eqs. (4.5)–(4.6)] The third a priori estimate, which yields the uniform bounds (4.5)–(4.6), is explicitly formal. At this point the available regularity is ∂tuε ∈ L2(0,T;V) and ξε ∈ L∞(0,T;H), so the test function −δΔ∂tuε + ∂tγ(uε) is not legitimate: −Δ∂tuε is only a distribution, and the displayed identity contains δ∫∇ξε·∇∂tuε and ∫γ′(uε)ξε∂tuε, neither of which is meaningful at this regularity. The authors state that a rigorous version can be obtained by replacing β with its Yosida approximation and that the resulting estimate is independent of λ, but no proof is supplied. These bounds are load-bearing: Section 4.5 uses (4.5)–(4.6) to pass to the limit and to identify ξ ∈ β(∂tu). Theorem 2.4 is therefore not fully proved as written. This is a fillable gap, but it must be closed by a complete regularization argument with constants independent of both ε and the regularization parameter.
- [Section 5.3, Eqs. (5.15)–(5.16)] The estimate leading to (5.15)–(5.16) is obtained by 'formally testing' (2.17) by −δ1/2Δ∂tuδ. At this stage ξδ is only known to lie in L∞(0,T;H), so the term δ1/2∫∇ξδ·∇∂tuδ is not a well-defined quantity, and no Yosida regularization or time-discretization argument is given for this computation. Since (5.15) and (5.16) are used in Section 5.4 to obtain compactness and the rate estimate (2.50), the proof of Theorem 2.6 has the same type of gap. If the regularization argument for Section 4.3 is carried out, the same argument should be adapted and written out for Section 5.3.
minor comments (4)
- [Theorem 2.6, Eq. (2.48)] Equation (2.48) contains a sign error: the limit equation should contain +g(t), not −g(t), consistent with (2.17) and with the computations in Sections 5.1 and 5.3. Please verify and correct this typo.
- [Section 5.4, convergence list] In the convergence list before (5.17), the line 'δ1/2u → 0 in H1(0,T;V)' should read 'δ1/2uδ → 0'; the subscript is missing.
- [Section 3.1, Proposition 3.1 application] The functional Φ is defined as Φ(v) := 1/2∫(|∇v|² + λ|v2|); this should be λ|v|².
- [Section 3.1, maximum principle] The lower-bound half of the maximum-principle argument is omitted with the note 'for brevity'. Since the upper-bound argument uses specific sign choices for b′0, the lower-bound analogue should be written out or the sign conventions should be explained, for completeness.
Circularity Check
No circularity: the load-bearing estimates are derived from the PDE and independent compactness/monotonicity arguments, and the self-citations are technical borrowings rather than assumed conclusions.
full rationale
Walking the claimed derivation chain: Theorem 2.2 is proved by a Yosida/truncation approximation in Section 3.1. Existence of the smooth approximating pair (u_lambda, mu_lambda) is imported from the authors' earlier paper [5], but [5] is a published, parameter-free existence theorem for the same regularized system under growth assumptions, and it does not contain the target no-growth boundedness result; the new maximum-principle argument using Proposition 3.1 and the Rothe discretization supplies the L-infinity control independently. The continuous-dependence estimate in Theorem 2.3 is obtained by a direct Gronwall estimate on differences, not by assuming the result. For Theorem 2.4, the epsilon-to-0 limit is based on the a priori estimates (4.1)-(4.6) derived from testing the equations, the Young inequality, the compactness inequality (2.1), and Lemma 4.1; the identification xi in beta(u_t) is a standard monotonicity and weak-strong closure argument. The delta-to-0 limit and the rate estimate (2.50) in Theorem 2.6 are likewise obtained from direct estimates, the maximum principle with constants depending only on ||mu_delta - g_delta||_{L^infty(Q)}, and the difference estimate (5.17). The only flagged weakness is that the third estimate in Section 4.3 is explicitly called formal and deferred to a Yosida regularization described only in words; this is a rigor gap that, if unfillable, would threaten Theorem 2.4, but it is not a circularity because the estimate is not assumed as the conclusion and the authors do not invoke the target theorem to justify it. There are no fitted inputs, no definitional coincidences, and no self-citation chain that forces the stated theorems. Hence the circularity score is 0.
Assumptions & free parameters
assumptions (7)
- standard math Proposition 3.1 (from Roubicek [28, Prop. 11.6]): existence and uniform estimates for time-discretized doubly nonlinear inclusions.
- standard math Lemma 2.1 (variation of Lions compactness lemma [19, Lem. 5.1]): interpolation inequality (2.1).
- standard math Aubin-Lions-Simon compactness lemma (Simon [31, Cor. 4]).
- standard math Convex analysis and maximal monotone operator theory, including Yosida approximations and the Young inequality for convex conjugates (Barbu [3]).
- standard math Elliptic regularity theory on smooth bounded domains.
- domain assumption Assumptions (2.2)-(2.8): psi in C2(a,b), psi >= 0, psi' blows up at the endpoints of (a,b), psi'' > -K, and beta = d beta-hat for a proper convex lower semicontinuous beta-hat with beta-hat(0) = 0.
- domain assumption The boundary value problem (1.1)-(1.4) with normal derivative of u equal to zero and mu equal to zero in a smooth bounded domain Omega in R3.
Cite this review
Pith. "Pith review of Bounded solutions and their asymptotics for a doubly nonlinear Cahn-Hilliard system." pith.science (2026). https://pith.science/paper/64JRHK46
@misc{pith2026190802079,
author = {Pith},
title = {Pith review of: Bounded solutions and their asymptotics for a doubly nonlinear Cahn-Hilliard system},
year = {2026},
howpublished = {\url{https://pith.science/paper/64JRHK46}},
note = {Machine review of arXiv:1908.02079}
}
read the original abstract
In this paper we deal with a doubly nonlinear Cahn-Hilliard system, where both an internal constraint on the time derivative of the concentration and a potential for the concentration are introduced. The definition of the chemical potential includes two regularizations: a viscosity and a diffusive term. First of all, we prove existence and uniqueness of a bounded solution to the system using a nonstandard maximum-principle argument for time-discretizations of doubly nonlinear equations. Possibly including singular potentials, this novel result brings improvements over previous approaches to this problem. Secondly, under suitable assumptions on the data, we show the convergence of solutions to the respective limit problems once either of the two regularization parameters vanishes.
Reference graph
Works this paper leans on
-
[5]
E. Bonetti, P. Colli, L. Scarpa and G. Tomassetti. A doubl y nonlinear Cahn-Hilliard system with nonlinear viscosity. Commun. Pure Appl. Anal. 17 (2018), 1001–1022
work page 2018
-
[1]
A. Agosti, P.F. Antonietti, P. Ciarletta, M. Grasselli, and M. Verani. A Cahn-Hilliard- type equation with application to tumor growth dynamics. Math. Methods Appl. Sci. 40 (2017), 7598–7626
work page 2017
-
[2]
F. Bagagiolo and A. Visintin. Hysteresis in filtration th rough porous media. Z. Anal. Anwendungen 19 (2000), 977–997
work page 2000
-
[3]
Nonlinear differential equations of monotone types in Banach spaces
V. Barbu. “Nonlinear differential equations of monotone types in Banach spaces”. Springer Monographs in Mathematics. Springer, New York (2010). Asymptotics for a doubly nonlinear CH system 25
work page 2010
-
[4]
E. Bonetti, P. Colli and G. Tomassetti. A non-smooth regu larization of a forward- backward parabolic equation. Math. Models Methods Appl. Sci. 27 (2017), 641–661
work page 2017
-
[6]
A.L. Bertozzi, S.Esedoglu, and A. Gillette. Inpainting of binary images using the Cahn–Hilliard equation. IEEE Trans. Image Process. 16 (2007), 285–291
work page 2007
- [7]
-
[8]
J.W. Cahn. On spinodal decomposition. Acta Metall. 9 (1961), 795–801
work page 1961
Show all 32 references
-
[9]
Cahn and J.E
J.W. Cahn and J.E. Hilliard. Free energy of a nonuniform s ystem. I. Interfacial free energy. J. Chem. Phys. 28 (1958), 258–267
1958
-
[10]
Colli, G
P. Colli, G. Gilardi, P. Podio-Guidugli and J. Sprekels . Well-posedness and long-time behavior for a nonstandard viscous Cahn–Hilliard system. SIAM J. Appl. Math. 71 (2011), 1849–1870
2011
-
[11]
Colli, G
P. Colli, G. Gilardi, P. Podio-Guidugli and J. Sprekels . Global existence and unique- ness for a singular/degenerate Cahn–Hilliard system with v iscosity. J. Differential Equations 254 (2013), 4217–4244
2013
-
[12]
Colli, G
P. Colli, G. Gilardi, E. Rocca and J. Sprekels. Vanishin g viscosities and error estimate for a Cahn–Hilliard type phase field system related to tumor g rowth. Nonlinear Anal. Real World Appl. 26 (2015), 93–108
2015
-
[13]
Colli, G
P. Colli, G. Gilardi, E. Rocca and J. Sprekels. Asymptot ic analyses and error esti- mates for a Cah–Hilliard type phase field system modelling tu mor growth. Discrete Contin. Dyn. Syst. Ser. S 10 (2017), 37–54
2017
-
[14]
Colli and L
P. Colli and L. Scarpa. From the viscous Cahn–Hilliard e quation to a regularized forward-backward parabolic equation. Asymptot. Anal. 99 (2016), 183–205
2016
-
[15]
P.C. Fife. Models for phase separation and their mathem atics. Electron. J. Differen- tial Equations 48 (2000), 26 pp
2000
-
[16]
C.G. Gal, M. Grasselli and A. Miranville. Cahn-Hilliar d-Navier-Stokes systems with moving contact lines. Calc. Var. Partial Differential Equations 55 (2016), Art. 50, 47 pp
2016
-
[17]
M.E. Gurtin. Generalized Ginzburg-Landau and Cahn-Hi lliard equations based on a microforce balance. Phys. D 92 (1996), 178–192
1996
-
[18]
Latroche
M. Latroche. Structural and thermodynamic properties of metallic hydrides used for energy storage. J. Phys. Chem. Solids 65 (2004) 517–522
2004
-
[19]
Quelques méthodes de résolution des probl èmes aux limites non linéaires
J. L. Lions. “Quelques méthodes de résolution des probl èmes aux limites non linéaires”. Dunod; Gauthier-Villars, Paris (1969). 26 Bonetti – Colli – Scarpa – Tomassetti
1969
-
[20]
Q.X. Liu, M. Rietkerk, P.M.J. Herman, T. Piersma, J.M. F ryxell, and J. van de Kop- pel. Phase separation driven by density-dependent movemen t: A novel mechanism for ecological patterns. Phys. Life Rev. , 19 (2016) 107–121
2016
-
[21]
Miranville
A. Miranville. Some generalizations of the Cahn–Hilli ard equation. Asymptot. Anal. 22 (2000) 235–259
2000
-
[22]
Miranville and G
A. Miranville and G. Schimperna. On a doubly nonlinear C ahn-Hilliard-Gurtin sys- tem. Discrete Contin. Dyn. Syst. Ser. B 14 (2010), 675–697
2010
-
[23]
Miranville and S
A. Miranville and S. Zelik. Doubly nonlinear Cahn-Hill iard-Gurtin equations. Hokkaido Math. J. 38 (2009), 315–360
2009
-
[24]
Miranville
A. Miranville. The Cahn–Hilliard equation and some of i ts variants. AIMS Math. 2 (2017), 479–544
2017
-
[25]
Material instabilities in continuum mechanics (Edinburgh, 1985–1986)
A. Novick-Cohen. On the viscous Cahn-Hilliard equatio n, in “Material instabilities in continuum mechanics (Edinburgh, 1985–1986)”. Oxford Sci. Publ. , Oxford Univ. Press, New York, (1988), 329–342
1988
-
[26]
Novick-Cohen and R
A. Novick-Cohen and R. L. Pego. Stable patterns in a visc ous diffusion equation. Trans. Amer. Math. Soc. 324 (1991), 331–351
1991
-
[27]
Podio-Guidugli
P. Podio-Guidugli. Models of phase segregation and diff usion of atomic species on a lattice. Ric. Mat. 55 (2006), 105–118
2006
-
[28]
Nonlinear partial differential equation s with applications
T. Roubíček. “Nonlinear partial differential equation s with applications”. Birkhäuser Verlag, Basel (2005)
2005
-
[29]
L. Scarpa. Existence and uniqueness of solutions to sin gular Cahn–Hilliard equations with nonlinear viscosity terms and dynamic boundary condit ions. J. Math. Anal. Appl. 469 (2019), 730–764
2019
-
[30]
Schweizer
B. Schweizer. Hysteresis in porous media: modelling an d analysis. Interfaces Free Bound. 19 (2017), 417–447
2017
-
[31]
J. Simon. Compact sets in the space Lp(0,T ;B). Ann. Mat. Pura Appl. (4) 146 (1987), 65–96
1987
-
[32]
Tomassetti
G. Tomassetti. Smooth and non-smooth regularizations of the nonlinear diffusion equation. Discrete Contin. Dyn. Syst. Ser. S 10 (2017) 1519-1537
2017
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.