Local Well-Posedness of the Cahn-Hilliard-Biot System
Pith reviewed 2026-05-18 21:36 UTC · model grok-4.3
The pith
The Cahn-Hilliard-Biot system admits unique short-time solutions for fluid flow through a two-phase deformable porous medium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The Cahn-Hilliard-Biot system is locally well-posed in time. For suitable initial data the coupled equations possess a unique solution on a positive but possibly small time interval, with continuous dependence on the data. The result covers both the purely elastic case, treated via semigroup methods on Hilbert spaces, and the viscoelastic case, treated via Banach scales to keep spatial regularity assumptions minimal. Phase-field dependent coefficients appear in the Biot part without destroying the contraction property.
What carries the argument
Reduction of the nonlinear coupled system to a fixed-point equation solved by contraction mapping, using maximal regularity theory (semigroup approach on Hilbert spaces without viscosity, Banach scales with viscosity).
If this is right
- Unique local-in-time solutions exist for initial data of sufficient regularity.
- The solution depends continuously on the initial data in the chosen function spaces.
- The result holds uniformly for both the viscoelastic and non-viscoelastic versions of the model.
- Material coefficients that vary with the phase field remain compatible with the short-time existence theory.
Where Pith is reading between the lines
- The local existence theory supplies a rigorous basis for starting numerical simulations of poroelastic phase separation before possible singularities appear.
- Similar fixed-point reductions could be attempted for other diffuse-interface models that couple phase fields to linear elasticity or Darcy flow.
- Global-in-time results might follow if additional structural assumptions, such as small initial energy or monotonicity, are imposed on top of the local theory.
Load-bearing premise
The nonlinear coupling between the phase field and the poroelastic stresses can be recast as a contraction mapping problem whose fixed point yields the solution via maximal regularity estimates.
What would settle it
Constructing smooth initial data for which the fixed-point map fails to be contractive, or for which no solution exists on any positive time interval, would show the local well-posedness claim is false.
read the original abstract
We show short-time well-posedness of a diffuse interface model describing the flow of a fluid through a deformable porous medium consisting of two phases. The system non-linearly couples Biot's equations for poroelasticity, including phase-field dependent material properties, with the Cahn-Hilliard equation to model the evolution of the solid, where we further distinguish between the absence and presence of a visco-elastic term of Kelvin-Voigt type. While both problems will be reduced to a fixed-point equation that can be solved using maximal regularity theory along with a contraction argument, the first case relies on a semigroup approach over suitable Hilbert spaces, whereas treating the second case under minimal assumptions with respect to spatial regularity necessitates the application of Banach scales.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proves short-time local well-posedness for a nonlinearly coupled Cahn-Hilliard-Biot system modeling two-phase fluid flow through a deformable porous medium. The system is reduced to a fixed-point equation solved by contraction mapping; the non-viscoelastic case uses semigroup theory on Hilbert spaces while the viscoelastic (Kelvin-Voigt) case employs maximal regularity on Banach scales under minimal spatial regularity assumptions on the coefficients.
Significance. If the contraction estimates close, the result supplies a rigorous existence theory for a diffuse-interface poroelasticity model with phase-dependent material parameters, which is relevant to applications in biomechanics and geomechanics. The technical distinction between the two cases and the use of Banach-scale maximal regularity to accommodate low-regularity coefficients constitute a clear advance over existing well-posedness results for Biot-type systems.
major comments (1)
- [Abstract and the fixed-point / maximal-regularity sections (viscoelastic case)] The central contraction argument (outlined in the abstract and presumably detailed in the fixed-point sections) applies maximal regularity to the Biot operator whose coefficients (elasticity tensor, permeability, etc.) depend nonlinearly on the phase field φ. Standard maximal-regularity theorems for parabolic systems with variable coefficients require the coefficients to lie in spaces such as W^{1,∞} or suitable multiplier spaces with controlled time derivatives. The Cahn-Hilliard component yields φ with regularity no better than L^∞(0,T;H¹) ∩ L²(0,T;H²) locally in time. The manuscript must explicitly verify that these coefficients satisfy the multiplier conditions (or supply a separate regularization/approximation step) so that the linear theory applies directly to the coupled system; without this verification the contraction estimate does not close.
minor comments (2)
- [Notation and setup] Define the precise Banach-scale norms and the precise dependence of the Biot coefficients on φ at the beginning of the viscoelastic-case analysis.
- [Introduction] Add a short remark comparing the obtained regularity to existing results for the pure Biot system or the pure Cahn-Hilliard equation.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive feedback on our manuscript concerning the local well-posedness of the Cahn-Hilliard-Biot system. The positive assessment of the significance is appreciated. We address the single major comment below and will incorporate clarifications to strengthen the presentation.
read point-by-point responses
-
Referee: The central contraction argument (outlined in the abstract and presumably detailed in the fixed-point sections) applies maximal regularity to the Biot operator whose coefficients (elasticity tensor, permeability, etc.) depend nonlinearly on the phase field φ. Standard maximal-regularity theorems for parabolic systems with variable coefficients require the coefficients to lie in spaces such as W^{1,∞} or suitable multiplier spaces with controlled time derivatives. The Cahn-Hilliard component yields φ with regularity no better than L^∞(0,T;H¹) ∩ L²(0,T;H²) locally in time. The manuscript must explicitly verify that these coefficients satisfy the multiplier conditions (or supply a separate regularization/approximation step) so that the linear theory applies directly to the coupled system; without this verification the contraction estimate does not close.
Authors: We thank the referee for this precise observation. For the viscoelastic case we deliberately invoke maximal regularity on Banach scales (as stated in the abstract and detailed in the fixed-point section) precisely because this framework accommodates coefficients with the limited regularity inherited from φ ∈ L^∞(0,T;H¹) ∩ L²(0,T;H²). The material coefficients are smooth (in fact analytic) functions of φ; under the cited abstract theorems for Banach-scale maximal regularity, such compositions remain admissible multipliers when φ satisfies the stated Sobolev regularity. Nevertheless, we agree that an explicit verification of the multiplier conditions for the elasticity tensor, permeability, and Biot coupling terms was not written out in sufficient detail. In the revised manuscript we will insert a short paragraph (or subsection) immediately after the statement of the linear theory, confirming that the coefficient maps satisfy the required time-integrability and multiplier estimates. This addition will make the contraction-mapping argument fully rigorous without introducing regularization or approximation steps. revision: yes
Circularity Check
No circularity: standard fixed-point reduction via external maximal regularity theory
full rationale
The derivation reduces the coupled Biot-Cahn-Hilliard system to a fixed-point equation solved by contraction mapping, invoking maximal regularity theory (or semigroups on Hilbert spaces) as an external tool. These are independent, pre-existing results from functional analysis that do not presuppose the target well-posedness statement. No self-definitional loops, fitted inputs renamed as predictions, or load-bearing self-citations appear in the abstract or described chain; the argument remains self-contained once the coefficient regularity conditions induced by the phase field are checked against the hypotheses of the cited linear theory.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Maximal regularity theory applies to the linearized operators arising after the fixed-point reduction
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
While both problems will be reduced to a fixed-point equation that can be solved using maximal regularity theory along with a contraction argument... the first case relies on a semigroup approach over suitable Hilbert spaces, whereas treating the second case under minimal assumptions... necessitates the application of Banach scales.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the aggregate corresponds to a dissipative operator on an appropriate Hilbert space and deduce that it generates an analytic semigroup
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
-
[2]
H. Abels and J. Weber , Local well-posedness of a quasi-incompressible two-phase flow , J. Evol. Equ., 21 (2021), pp. 3477–3502
work page 2021
- [3]
-
[4]
H. Amann, Linear and Quasilinear Parabolic Problems: Volume I: Abstract Linear Theory , Monographs in Mathe- matics, Springer International Publishing, 1995
work page 1995
-
[5]
H. Amann and J. Escher , Analysis I, Grundstudium Mathematik, Birkh¨ auser Basel, 2013
work page 2013
- [6]
-
[7]
Auriault, Dynamic behaviour of a porous medium saturated by a Newtonian fluid , Int
J. Auriault, Dynamic behaviour of a porous medium saturated by a Newtonian fluid , Int. J. Eng. Sci., 18 (1980), pp. 775–785
work page 1980
- [8]
-
[9]
A. Behzadan and M. Holst , Multiplication in Sobolev spaces, revisited, Ark. Mat., 59 (2021), pp. 275–306
work page 2021
-
[10]
M. A. Biot, General theory of three-dimensional consolidation , J. Appl Phys., 12 (1941), pp. 155–164
work page 1941
-
[11]
, Theory of deformation of a porous viscoelastic anisotropic solid , J. Appl. Phys., 27 (1956), pp. 459–467
work page 1956
-
[12]
M. A. Biot and D. G. Willis , The elastic coefficients of the theory of consolidation , J. Appl. Mech., (1957)
work page 1957
-
[13]
J. F. Blowey and C. M. Elliott , The Cahn–Hilliard gradient theory for phase separation with non-smooth free energy Part I: Mathematical analysis , European J. Appl. Math., 2 (1991), pp. 233–280
work page 1991
- [14]
- [15]
-
[16]
E. Bonetti, P. Colli, W. Dreyer, G. Gilardi, G. Schimperna, and J. Sprekels , On a model for phase separation in binary alloys driven by mechanical effects , Phys. D, 165 (2002), pp. 48–65
work page 2002
-
[17]
J. W. Both, I. S. Pop, and I. Yotov, Global existence of weak solutions to unsaturated poroelasticity, ESAIM: M2AN, 55 (2021), pp. 2849–2897
work page 2021
-
[18]
A. Brunk and M. Fritz , Structure-preserving approximation of the Cahn-Hilliard-Biot system , Numer. Methods Partial Differential Equations, 41 (2025), pp. Paper No. e23159, 14
work page 2025
-
[19]
S.-S. Byun and L. Wang , Gradient estimates for elliptic systems in non-smooth domains , Math. Ann., 341 (2008), pp. 629–650
work page 2008
-
[20]
J. Cahn and F. Larch ´e, The effect of self-stress on diffusion in solids , Acta Metall., 30 (1982), pp. 1835–1845
work page 1982
-
[21]
J. W. Cahn, C. M. Elliott, and A. Novick-Cohen , The Cahn–Hilliard equation with a concentration dependent mobility: Motion by minus the laplacian of the mean curvature , European J. Appl. Math., 7 (1996), pp. 287–301
work page 1996
-
[22]
J. W. Cahn and J. E. Hilliard , Free energy of a nonuniform system. I. Interfacial free energy , J. Chem. Phys., 28 (1958), pp. 258–267
work page 1958
-
[23]
Y. Cao, S. Chen, and A. J. Meir , Analysis and numerical approximations of equations of nonlinear poroelasticity , Discrete Contin. Dyn. Syst. Ser. B, 18 (2013), pp. 1253–1273
work page 2013
-
[24]
M. Carrive, A. Miranville, and A. Pi ´etrus, The Cahn–Hilliard equation for deformable elastic continua , Adv. Math. Sci. Appl., 10 (2000), pp. 539–569
work page 2000
-
[25]
G. Cavalleri, A phase field model of Cahn-Hilliard type for tumour growth with mechanical effects and damage , J. Math. Anal. Appl., 550 (2025), pp. Paper No. 129627, 40. 44 HELMUT ABELS AND JONAS HASELB ¨OCK
work page 2025
-
[26]
P. Cl´ement and J. Pr ¨uss, An operator-valued transference principle and maximal regularity on vector-valued Lp- spaces, in Evolution equations and their applications in physical and life sciences (Bad Herrenalb, 1998), vol. 215 of Lecture Notes in Pure and Appl. Math., Dekker, New York, 2001, pp. 67–87
work page 1998
- [27]
- [28]
-
[29]
L. de Simon, Un’applicazione della teoria degli integrali singolari allo studio delle equazioni differenziali lineari astratte del primo ordine , Rend. Sem. Mat. Univ. Padova, 34 (1964), pp. 205–223
work page 1964
-
[30]
R. Denk, G. Dore, M. Hieber, J. Pr ¨uss, and A. Venni , New thoughts on old results of R. T. Seeley , Math. Ann., 328 (2004), pp. 545–583
work page 2004
-
[31]
R. Denk, M. Hieber, and J. Pr ¨uss, R-boundedness, Fourier multipliers and problems of elliptic and parabolic type , Mem. Amer. Math. Soc., 166 (2003), pp. viii+114
work page 2003
-
[32]
K. Disser, H.-C. Kaiser, and J. Rehberg , Optimal Sobolev regularity for linear second-order divergence elliptic operators occurring in real-world problems, SIAM J. Math. Anal., 47 (2015), pp. 1719–1746
work page 2015
-
[33]
G. Dore, Lp regularity for abstract differential equations, in Functional Analysis and Related Topics, 1991, H. Komatsu, ed., Berlin, Heidelberg, 1993, Springer Berlin Heidelberg, pp. 25–38
work page 1991
-
[34]
M. Ebenbeck and H. Garcke , Analysis of a Cahn–Hilliard–Brinkman model for tumour growth with chemotaxis , J. Diff. Equ., 266 (2019), pp. 5998–6036
work page 2019
-
[35]
C. Elliott and H. Garcke, On the Cahn–Hilliard equation with degenerate mobility, SIAM J. Math. Anal., 27 (1996), pp. 404–423
work page 1996
-
[36]
J. Elschner, J. Rehberg, and G. Schmidt , Optimal regularity for elliptic transmission problems including C1 interfaces, Interfaces Free Bound., 9 (2007), pp. 233–252
work page 2007
-
[37]
K. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution Equations, Graduate Texts in Mathematics, Springer New York, 1999
work page 1999
-
[38]
M. Fritz, On the well-posedness of the Cahn-Hilliard-Biot model and its applications to tumor growth, Discrete Contin. Dyn. Syst. Ser. S, 17 (2024), pp. 3533–3563
work page 2024
-
[39]
Garcke, On Cahn–Hilliard systems with elasticity , Proc
H. Garcke, On Cahn–Hilliard systems with elasticity , Proc. Roy. Soc. Edinburgh Sect. A, 133 (2003), p. 307–331
work page 2003
-
[40]
Garcke, On a Cahn–Hilliard model for phase separation with elastic misfit , Ann
H. Garcke, On a Cahn–Hilliard model for phase separation with elastic misfit , Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire, 22 (2005), pp. 165–185
work page 2005
- [41]
- [42]
-
[43]
H. Garcke and K. F. Lam , Global weak solutions and asymptotic limits of a Cahn–Hilliard–Darcy system modelling tumour growth, AIMS Mathematics, 1 (2016), pp. 318–360
work page 2016
-
[44]
J. Geng, W 1,p estimates for elliptic problems with Neumann boundary conditions in Lipschitz domains , Adv. Math., 229 (2012), pp. 2427–2448
work page 2012
-
[45]
M. Giaquinta and L. Martinazzi, An introduction to the regularity theory for elliptic systems, harmonic maps and minimal graphs, vol. 11 of Appunti. Scuola Normale Superiore di Pisa (Nuova Serie) [Lecture Notes. Scuola Normale Superiore di Pisa (New Series)], Edizioni della Normale, Pisa, second ed., 2012
work page 2012
-
[46]
D. Gilbarg and N. Trudinger , Elliptic Partial Differential Equations of Second Order , Classics in Mathematics, Springer Berlin Heidelberg, 2001
work page 2001
-
[47]
K. Gr ¨oger, A W 1,p-estimate for solutions to mixed boundary value problems for second order elliptic differential equations, Mathematische Annalen, 283 (1989), pp. 679–687
work page 1989
-
[48]
R. Haller-Dintelmann, A. Jonsson, D. Knees, and J. Rehberg , Elliptic and parabolic regularity for second-order divergence operators with mixed boundary conditions , Math. Methods Appl. Sci., 39 (2016), pp. 5007–5026
work page 2016
-
[49]
E. Holland and R. E. Showalter , Poro-visco-elastic compaction in sedimentary basins , SIAM J. Math. Anal., 50 (2018), pp. 2295–2316
work page 2018
- [50]
-
[51]
A. Hosseinkhan and R. E. Showalter, Semilinear degenerate Biot-Signorini system, SIAM J. Math. Anal., 55 (2023), pp. 5643–5665
work page 2023
-
[52]
J. Lowengrub, E. Titi, and K. Zhao , Analysis of a mixture model of tumor growth , European J.Appl. Math., 24 (2013), pp. 691–734
work page 2013
-
[53]
A. McIntosh , Operators which have an H∞ functional calculus, in Miniconference on operator theory and partial differential equations (North Ryde, 1986), vol. 14 of Proc. Centre Math. Anal. Austral. Nat. Univ., Austral. Nat. Univ., Canberra, 1986, pp. 210–231
work page 1986
-
[54]
W. McLean, Strongly elliptic systems and boundary integral equations, Cambridge University Press, Cambridge, 2000
work page 2000
-
[55]
C. B. Morrey, Jr. , Multiple integrals in the calculus of variations , vol. Band 130 of Die Grundlehren der mathema- tischen Wissenschaften, Springer-Verlag New York, Inc., New York, 1966
work page 1966
-
[56]
V. C. Mow, S. Kuei, W. M. Lai, and C. G. Armstrong , Biphasic creep and stress relaxation of articular cartilage in compression: Theory and experiments , ASME J. Biomech. Eng., (1980)
work page 1980
-
[57]
S. Nicaise, About the Lam´ e system in a polygonal or a polyhedral domain and a coupled problem between the Lam´ e system and the plate equation. I. Regularity of the solutions , Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 19 (1992), pp. 327–361
work page 1992
-
[58]
Onuki, Ginzburg–Landau approach to elastic effects in the phase separation of solids , J
A. Onuki, Ginzburg–Landau approach to elastic effects in the phase separation of solids , J. Phys. Soc. Jap., 58 (1989), pp. 3065–3068
work page 1989
-
[59]
Y. Oono and S. Puri , Study of phase-separation dynamics by use of cell dynamical systems. I. Modeling , Phys. Rev. A, 38 (1988), pp. 434–453. LOCAL WELL-POSEDNESS OF THE CAHN–HILLIARD–BIOT SYSTEM 45
work page 1988
-
[60]
A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations , Applied Mathematical Sciences, Springer New York, 2012
work page 2012
-
[61]
J. Pr¨uss and G. Simonett , Moving interfaces and quasilinear parabolic evolution equations , vol. 105 of Monographs in Mathematics, Birkh¨ auser/Springer, [Cham], 2016
work page 2016
-
[62]
M. Renardy and R. C. Rogers , An introduction to partial differential equations , vol. 13 of Texts in Applied Mathe- matics, Springer-Verlag, New York, second ed., 2003
work page 2003
-
[63]
C. Riethm ¨uller, E. Storvik, J. W. Both, and F. A. Radu , Well-posedness analysis of the Cahn-Hilliard-Biot model, Nonlinear Anal. Real World Appl., 84 (2025), pp. Paper No. 104271, 27
work page 2025
-
[64]
T. Runst and W. Sickel, Sobolev Spaces of Fractional Order, Nemytskij Operators, and Nonlinear Partial Differential Equations, De Gruyter Series in Nonlinear Analysis and Applications, De Gruyter, 2011
work page 2011
- [65]
-
[66]
Savar´e, Regularity and perturbation results for mixed second order elliptic problems , Comm
G. Savar´e, Regularity and perturbation results for mixed second order elliptic problems , Comm. Partial Differential Equations, 22 (1997), pp. 869–899
work page 1997
-
[67]
Seeley, Interpolation in lp with boundary conditions, Studia Mathematica, 44 (1972), pp
R. Seeley, Interpolation in lp with boundary conditions, Studia Mathematica, 44 (1972), pp. 47–60
work page 1972
-
[68]
Shamir, Regularization of mixed second-order elliptic problems , Isr
E. Shamir, Regularization of mixed second-order elliptic problems , Isr. J. Math., 6 (1968), pp. 150–168
work page 1968
-
[69]
Showalter, Diffusion in poro-elastic media , J
R. Showalter, Diffusion in poro-elastic media , J. Math. Anal., 251 (2000), pp. 310–340
work page 2000
-
[70]
R. E. Showalter, Diffusion in poro-elastic media , J. Math. Anal. Appl., 251 (2000), pp. 310–340
work page 2000
-
[71]
R. E. Showalter and U. Stefanelli, Diffusion in poro-plastic media, Math. Methods Appl. Sci., 27 (2004), pp. 2131– 2151
work page 2004
-
[72]
R. E. Showalter and N. Su, Partially saturated flow in a poroelastic medium , Discrete Contin. Dyn. Syst. Ser. B, 1 (2001), pp. 403–420
work page 2001
-
[73]
E. Storvik, J. W. Both, J. M. Nordbotten, and F. A. Radu , A Cahn–Hilliard–Biot system and its generalized gradient flow structure , Appl. Math. Lett., 126 (2022), p. 107799
work page 2022
-
[74]
E. Storvik and C. Bringedal , Sharp-interface limit of the Cahn–Hilliard–Biot equations , Appl. Math. Lett., 166 (2025), p. Paper No. 109522
work page 2025
-
[75]
E. Storvik, C. Riethm ¨uller, J. W. Both, and F. A. Radu , Sequential solution strategies for the Cahn-Hilliard- Biot model, in Numerical mathematics and advanced applications—ENUMATH 2023. Vol. 2, vol. 154 of Lect. Notes Comput. Sci. Eng., Springer, Cham, 2025, pp. 369–378
work page 2023
-
[76]
Valent, Boundary value problems of finite elasticity , vol
T. Valent, Boundary value problems of finite elasticity , vol. 31 of Springer Tracts in Natural Philosophy, Springer- Verlag, New York, 1988. Local theorems on existence, uniqueness, and analytic dependence on data
work page 1988
-
[77]
L. Weis, The H ∞-holomorphic functional calculus for sectorial operators – a survey , in Partial Differential Equations and Functional Analysis: The Philippe Cl´ ement Festschrift, E. Koelink, J. van Neerven, B. de Pagter, G. Sweers, A. Luger, and H. Woracek, eds., Birkh¨ auser Basel, Basel, 2006, pp. 263–294
work page 2006
-
[78]
Werner, Funktionalanalysis, Springer-Lehrbuch, Springer Berlin Heidelberg, 2018
D. Werner, Funktionalanalysis, Springer-Lehrbuch, Springer Berlin Heidelberg, 2018
work page 2018
-
[79]
A. ˇZen´ıˇsek, The existence and uniqueness theorem in Biot’s consolidation theory , Aplikace matematiky, 29 (1984), pp. 194–211. Appendix A. W2,p-regularity for elliptic systems The aim of this section is to provide a proof of the regularity result for elliptic systems with mixed boundary conditions as claimed in Theorem 3.2. While there are many classic...
work page 1984
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.