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.
Funding statement: The authors gratefully acknowledge the support by the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. HKBU 14302319) and Hong Kong Baptist University (Project No. RC-OFSGT2/20-21/SCI/006).
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–196.10.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. Engrg. 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. Engrg. 2016 (2016), Art. ID 9532608.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), 10.1093/imanum/drac078.Search 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, Comput. Math. Appl. (2022) 10.1016/j.camwa.2022.09.031.Search 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
© 2024 Walter de Gruyter GmbH, Berlin/Boston