Abstract
The scalar auxiliary variable (SAV) approach of Shen et al. (2018), which presents a novel way to discretize a large class of gradient flows, has been extended and improved by many authors for general dissipative systems. In this work we consider a Cahn–Hilliard system with mass source that, for image processing and biological applications, may not admit a dissipative structure involving the Ginzburg–Landau energy. Hence, compared to previous works, the stability of SAV-discrete solutions for such systems is not immediate. We establish, with a bounded mass source, stability and convergence of time discrete solutions for a first-order relaxed SAV scheme in the sense of Jiang et al. (2022), and apply our ideas to Cahn–Hilliard systems with mass source appearing in diblock co-polymer phase separation, tumor growth, image inpainting and segmentation.
References
[1] A.C. Aristotelous, O.A. Karakashian and S.M. Wise, Adaptive, second-order in time, primitive-variable discontinuous Galerkin schemes for a Cahn–Hilliard equation with a mass source, IMA J. Numer. Anal. 35 (2015), 1167–1198.10.1093/imanum/dru035Search in Google Scholar
[2] A. Bertozzi, S. Esdoḡlu and A. Gillette, Analysis of a two-scale Cahn–Hilliard model for binary image inpainting, Multiscale Model. Simul. 6 (2007), 913–936.10.1137/060660631Search in Google Scholar
[3] A. Bertozzi, S. Esdoḡlu and A. Gillette, Inpainting of binary images using the Cahn–Hilliard equation, IEEE Trans. Image Process. 16 (2007), 285–291.10.1109/TIP.2006.887728Search in Google Scholar PubMed
[4] H. Chen, J. Mao and J. Shen, Optimal error estimates for the scalar auxiliary variable finite-element schemes for gradient flows, Numer. Math. 145 (2020), 167–19610.1007/s00211-020-01112-4Search in Google Scholar
[5] Q. Cheng, C. Liu and J. Shen, A new Lagrange multiplier approach for gradient flows, Comput. Method Appl. Math. Engrg. 367 (2020), 113070.10.1016/j.cma.2020.113070Search in Google Scholar
[6] Q. Cheng, X. Yang and J. Shen, Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model, J. Comput. Phys. 341 (2017), 44–60.10.1016/j.jcp.2017.04.010Search in Google Scholar
[7] L. Cherfils, H. Fakih and A. Miranville, Finite-dimensional attractors for the Bertozzi-Esedoglu-Gillette-Cahn-Hilliard equation in image inpainting, Inv. Probl. Imag. 9 (2015, 105–125.10.3934/ipi.2015.9.105Search in Google Scholar
[8] P. Clément, Approximation by finite element functions using local regularization, ESAIM: M2NA 9 (1975), 77–84.10.1051/m2an/197509R200771Search in Google Scholar
[9] V. Cristini and J. Lowengrub, Multiscale Modeling of Cancer: An Integrated Experimental and Mathematical Modeling Approach, Cambridge University Press, 2010.10.1017/CBO9780511781452Search in Google Scholar
[10] C.M. Elliott and A.M. Stuart, The global dynamics of discrete semilinear parabolic equations, SIAM J. Numer. Anal. 30 (1993), 1622–1663.10.1137/0730084Search in Google Scholar
[11] D.J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation, Mater. Res. Soc. Symp. Proc. 529 (1998), 39–46.10.1557/PROC-529-39Search in Google Scholar
[12] H. Fakih, R. Mghames and N. Nasreddine, On the Cahn–Hilliard equation with mass source for biological applications, Commun. Pure Appl. Anal. 20 (2021), 495–510.10.3934/cpaa.2020277Search in Google Scholar
[13] F. Huang and J. Shen, A new class of implicit-explicit BDFk SAV schemes for general dissipative systems and their error analysis, Comput. Methods Appl. Mech. Engrg. 392 (2022), 114718.10.1016/j.cma.2022.114718Search in Google Scholar
[14] F. Huang, J. Shen and Z. Yang, A highly efficient and accurate new scalar auxiliary variable approach for gradient flows, SIAM J. Sci. Comput. 42 (2020), A2514–A2536.10.1137/19M1298627Search in Google Scholar
[15] M. Jiang, Z. Zhang and J. Zhao, Improving the accuracy and consistency of the scalar auxiliary variable (sav) method with relaxation, J. Comput. Phys. (2022), 110954.10.1016/j.jcp.2022.110954Search in Google Scholar
[16] H. Garcke and K.F. Lam, Well-posedness of a Cahn–Hilliard system modelling tumour growth with chemotaxis and active transport, Eur. J. Appl. Math. 28 (2017), 284–316.10.1017/S0956792516000292Search in Google Scholar
[17] H. Garcke, K.F. Lam, E. Sitka and V. Styles, A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport, Math. Models Methods Appl. Sci. 26 (2016), 1095–1148.10.1142/S0218202516500263Search in Google Scholar
[18] H. Garcke, and D. Trautwein, Numerical analysis for a Cahn–Hilliard system modelling tumour growth with chemotaxis and active transport, J. Numer. Math. 30 (2022), 295–324.10.1515/jnma-2021-0094Search in Google Scholar
[19] A. Giorgini, M. Grasselli and A. Miranville, The Cahn–Hilliard–Oono equation with singular potential, Math. Models Methods Appl. Sci. 27 (2017), 2485–2510.10.1142/S0218202517500506Search in Google Scholar
[20] P. Grisvard, Elliptic Problems in Nonsmooth Domains, SIAM, Philadelphia, 2011.10.1137/1.9781611972030Search in Google Scholar
[21] F. Guillén-González and G. Tierra, On linear schemes for a Cahn–Hilliard diffuse interface model, J. Comput. Phys. 234 (2013), 140–171.10.1016/j.jcp.2012.09.020Search in Google Scholar
[22] A. Hawkins-Daarud, K.G. van der Zee and J.T. Oden, Numerical simulation of a thermodynamically consistent four-species tumor growth model, Int. J. Numer. Method Biomed. Eng. 28 (2012), 3–24.10.1002/cnm.1467Search in Google Scholar PubMed
[23] D. Hou, M. Azaiex and C. Xu, A variant of scalar auxiliary variable approaches for gradient flows, J. Comput. Phys. 395 (2019), 307–332.10.1016/j.jcp.2019.05.037Search in Google Scholar
[24] J. Kim, S. Lee, Y. Choi, S.-M. Lee and D. Jeong, Basic Principles and Practical Applications of the Cahn–Hilliard Equation, Math. Probl. Eng. 2016 (2016), Article ID 9532608, 11 pages.10.1155/2016/9532608Search in Google Scholar
[25] P. Knopf, K.F. Lam, C. Liu and S. Metzger, Phase-field dynamics with transfer of materials: The Cahn–Hilliard equation with reaction rate dependent dynamical boundary conditions, ESAIM: M2AN 55 (2021), 229–282.10.1051/m2an/2020090Search in Google Scholar
[26] X. Li and J. Shen, On a SAV-MAC scheme for the Cahn–Hilliard–Navier–Stokes phase field model and its error analysis for the corresponding Cahn–Hilliard–Stokes case, Math. Models Methods Appl. Sci. 30 (2020), 2263–2297.10.1142/S0218202520500438Search in Google Scholar
[27] Z. Liu and X. Li, The exponential scalar auxiliary variable (E-SAV) approach for phase field models and its explicit computing, SIAM J. Sci. Comput. 42 (2020), pp. B630–B655.10.1137/19M1305914Search in Google Scholar
[28] S. Metzger, A convergent SAV scheme for Cahn–Hilliard equations with dynamic boundary conditions, IMA J. Numer. Anal., (2023), doi:10.1093/imanum/drac078.10.1093/imanum/drac078Search in Google Scholar
[29] L. Li, A. Miranville and R. Guillevin, Cahn–Hilliard models for glial cells, Appl. Math. Optim. 84 (2021), 1821–1842.10.1007/s00245-020-09696-xSearch in Google Scholar
[30] A. Miranville, E. Rocca and G. Schimperna, On the long time behavior of a tumor growth model, J. Differential Equations 267 (2019), 2616–2642.10.1016/j.jde.2019.03.028Search in Google Scholar
[31] T. Ohta and K. Kawasaki, Equilibrium morphology of block copolymer melts, Marcomolecules 19 (1986), 2621– 2632.10.1021/ma00164a028Search in Google Scholar
[32] Y. Oono and S. Puri, Computationally efficient modeling of ordering of quenched phases, Phys. Rev. Lett. 58 (1987), 836–839.10.1103/PhysRevLett.58.836Search in Google Scholar PubMed
[33] J. Shen and J. Xu, Convergence and error analysis of the scalar auxiliary variable (SAV) schemes to gradient flows, SIAM J. Numer. Anal. 56 (2018), 2895–2912.10.1137/17M1159968Search in Google Scholar
[34] J. Shen, J. Xu and J. Yang, The scalar auxiliary variable (SAV) approach for gradient flows, J. Comput. Phys. 353 (2018), 407–416.10.1016/j.jcp.2017.10.021Search in Google Scholar
[35] J. Shen, J. Xu and J. Yang, A new class of efficient and robust energy stable schemes for gradient flows, SIAM Review 61 (2019), 474–506.10.1137/17M1150153Search in Google Scholar
[36] X. Shen, L. Wu, J. Wen and J. Zhang, SAV Fourier-spectral method for diffuse-interface tumor-growth model, Com-put. Math. with Appl. (2022) DOI: 10.1016/j.camwa.2022.09.031.10.1016/j.camwa.2022.09.031Search in Google Scholar
[37] J. Shen and X. Yang, Numerical approximations of Allen-Cahn and Cahn-Hilliard equations, Discrete Contin. Dyn. Syst. 28 (2010), 1669–1691.10.3934/dcds.2010.28.1669Search in Google Scholar
[38] J. Shen and X. Yang, The IEQ and SAV approaches and their extensions for a class of highly nonlinear gradient flow systems, Contemp. Math. 754 (2020), 217–245.10.1090/conm/754/15147Search in Google Scholar
[39] J. Simon, Compact Sets in the Space Lp(0 T;B), Ann. di Mat. Pura ed Appl. (IV) 146 (1987), 65–96.10.1007/BF01762360Search in Google Scholar
[40] S.M. Wise, J.S. Lowengrub, H.B. Frieboes and V. Cristini, Three-dimensional multispecies nonlinear tumor growth - I: Model and numerical method, J. Theor. Biol. 253 (2008), 524–543.10.1016/j.jtbi.2008.03.027Search in Google Scholar PubMed PubMed Central
[41] W. Yang, Z. Huang and W. Zhu, Image segmentation using the Cahn–Hilliard equation, J. Sci. Comput. 79 (2019), 1057–1077.10.1007/s10915-018-00899-7Search in Google Scholar
[42] Z. Yang and S. Dong, A roadmap for discretely energy-stable schemes for dissipative systems based on a generalized auxiliary variable with guaranteed positivity, J. Comput. Phys. 404 (2020), 109121.10.1016/j.jcp.2019.109121Search in Google Scholar
[43] Y. Zhang and J. Shen, A generalized SAV approach with relaxation for dissipative systems, J. Comput. Phys. 464 (2022), 111311.10.1016/j.jcp.2022.111311Search in Google Scholar
[44] N. Zheng and X. Li, Error analysis of the SAV Fourier-spectral method for the Cahn-Hilliard-Hele-Shaw system, Adv. Comput. Math. 47 (2021), 71.10.1007/s10444-021-09897-0Search in Google Scholar
© 2023 Walter de Gruyter GmbH, Berlin/Boston
