Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential
HTML articles powered by AMS MathViewer
- by Yunzhuo Guo, Cheng Wang, Steven M. Wise and Zhengru Zhang;
- Math. Comp. 93 (2024), 2185-2214
- DOI: https://doi.org/10.1090/mcom/3916
- Published electronically: November 20, 2023
- HTML | PDF | Request permission
Abstract:
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment for the singular logarithmic terms and the surface diffusion term, and an explicit update for the expansive concave term. The convective term for the phase variable, as well as the coupled term in the Stokes equation, is approximated in a semi-implicit manner. In the spatial discretization, the marker and cell difference method is applied, which evaluates the velocity components, the pressure and the phase variable at different cell locations. Such an approach ensures the divergence-free feature of the discrete velocity, and this property plays an important role in the analysis. The positivity-preserving property and the unique solvability of the proposed numerical scheme are theoretically justified, utilizing the singular nature of the logarithmic term as the phase variable approaches the singular limit values. An unconditional energy stability analysis is standard, as an outcome of the convex-concave decomposition technique. A convergence analysis with accompanying error estimate is provided for the proposed numerical scheme. In particular, a higher order consistency analysis, accomplished by supplementary functions, is performed to ensure the separation properties of numerical solution. In turn, using the approach of rough and refined error estimates, we are able to derive an optimal rate convergence. To conclude, several numerical experiments are presented to validate the theoretical analysis.References
- John W. Barrett and James F. Blowey, Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility, Math. Comp. 68 (1999), no. 226, 487–517. MR 1609678, DOI 10.1090/S0025-5718-99-01015-7
- Wenbin Chen, Wenqiang Feng, Yuan Liu, Cheng Wang, and Steven M. Wise, A second order energy stable scheme for the Cahn-Hilliard-Hele-Shaw equations, Discrete Contin. Dyn. Syst. Ser. B 24 (2019), no. 1, 149–182. MR 3932721, DOI 10.3934/dcdsb.2018090
- Wenbin Chen, Shufen Wang, Yichao Zhang, Daozhi Han, Cheng Wang, and Xiaoming Wang, Error estimate of a decoupled numerical scheme for the Cahn-Hilliard-Stokes-Darcy system, IMA J. Numer. Anal. 42 (2022), no. 3, 2621–2655. MR 4454932, DOI 10.1093/imanum/drab046
- W. Chen, J. Jing, Q. Liu, C. Wang, and X. Wang, A second order accurate, positivity-preserving numerical scheme of the Cahn-Hilliard-Navier-Stokes system with Flory-Huggins potential, Commun. Comput. Phys., 2023. accepted and in press.
- Wenbin Chen, Jianyu Jing, Cheng Wang, and Xiaoming Wang, A positivity preserving, energy stable finite difference scheme for the Flory-Huggins-Cahn-Hilliard-Navier-Stokes system, J. Sci. Comput. 92 (2022), no. 2, Paper No. 31, 24. MR 4443516, DOI 10.1007/s10915-022-01872-1
- Wenbin Chen, Jianyu Jing, Cheng Wang, Xiaoming Wang, and Steven M. Wise, A modified Crank-Nicolson numerical scheme for the Flory-Huggins Cahn-Hilliard model, Commun. Comput. Phys. 31 (2022), no. 1, 60–93. MR 4350698, DOI 10.4208/cicp.oa-2021-0074
- Wenbin Chen, Yuan Liu, Cheng Wang, and Steven M. Wise, Convergence analysis of a fully discrete finite difference scheme for the Cahn-Hilliard-Hele-Shaw equation, Math. Comp. 85 (2016), no. 301, 2231–2257. MR 3511281, DOI 10.1090/mcom3052
- Wenbin Chen, Cheng Wang, Xiaoming Wang, and Steven M. Wise, Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential, J. Comput. Phys. X 3 (2019), 100031, 29. MR 4116082, DOI 10.1016/j.jcpx.2019.100031
- Craig Collins, Jie Shen, and Steven M. Wise, An efficient, energy stable scheme for the Cahn-Hilliard-Brinkman system, Commun. Comput. Phys. 13 (2013), no. 4, 929–957. MR 2982944, DOI 10.4208/cicp.171211.130412a
- Francesco Della Porta, Andrea Giorgini, and Maurizio Grasselli, The nonlocal Cahn-Hilliard-Hele-Shaw system with logarithmic potential, Nonlinearity 31 (2018), no. 10, 4851–4881. MR 3856121, DOI 10.1088/1361-6544/aad52a
- Amanda E. Diegel, Cheng Wang, and Steven M. Wise, Stability and convergence of a second-order mixed finite element method for the Cahn-Hilliard equation, IMA J. Numer. Anal. 36 (2016), no. 4, 1867–1897. MR 3556407, DOI 10.1093/imanum/drv065
- Amanda E. Diegel, Cheng Wang, Xiaoming Wang, and Steven M. Wise, Convergence analysis and error estimates for a second order accurate finite element method for the Cahn-Hilliard-Navier-Stokes system, Numer. Math. 137 (2017), no. 3, 495–534. MR 3712284, DOI 10.1007/s00211-017-0887-5
- Lixiu Dong, Cheng Wang, Steven M. Wise, and Zhengru Zhang, A positivity-preserving, energy stable scheme for a ternary Cahn-Hilliard system with the singular interfacial parameters, J. Comput. Phys. 442 (2021), Paper No. 110451, 29. MR 4270095, DOI 10.1016/j.jcp.2021.110451
- Lixiu Dong, Cheng Wang, Hui Zhang, and Zhengru Zhang, A positivity-preserving, energy stable and convergent numerical scheme for the Cahn-Hilliard equation with a Flory-Huggins-deGennes energy, Commun. Math. Sci. 17 (2019), no. 4, 921–939. MR 4030506, DOI 10.4310/CMS.2019.v17.n4.a3
- Lixiu Dong, Cheng Wang, Hui Zhang, and Zhengru Zhang, A positivity-preserving second-order BDF scheme for the Cahn-Hilliard equation with variable interfacial parameters, Commun. Comput. Phys. 28 (2020), no. 3, 967–998. MR 4126520, DOI 10.4208/cicp.oa-2019-0037
- Chenghua Duan, Wenbin Chen, Chun Liu, Cheng Wang, and Xingye Yue, A second order accurate, energy stable numerical scheme for the one-dimensional porous medium equation by an energetic variational approach, Commun. Math. Sci. 20 (2022), no. 4, 987–1024. MR 4407510, DOI 10.4310/CMS.2022.v20.n4.a3
- Chenghua Duan, Wenbin Chen, Chun Liu, Cheng Wang, and Shenggao Zhou, Convergence analysis of structure-preserving numerical methods for nonlinear Fokker-Planck equations with nonlocal interactions, Math. Methods Appl. Sci. 45 (2022), no. 7, 3764–3781. MR 4407950, DOI 10.1002/mma.8015
- Chenghua Duan, Chun Liu, Cheng Wang, and Xingye Yue, Convergence analysis of a numerical scheme for the porous medium equation by an energetic variational approach, Numer. Math. Theory Methods Appl. 13 (2020), no. 1, 63–80. MR 4044415, DOI 10.4208/nmtma.oa-2019-0073
- D. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation, In J. W. Bullard, R. Kalia, M. Stoneham, and L.Q. Chen, editors, Computational and Mathematical Models of Microstructural Evolution, volume 53, pages 1686–1712, Warrendale, PA, USA, 1998. Materials Research Society.
- Xiaobing Feng and Steven Wise, Analysis of a Darcy-Cahn-Hilliard diffuse interface model for the Hele-Shaw flow and its fully discrete finite element approximation, SIAM J. Numer. Anal. 50 (2012), no. 3, 1320–1343. MR 2970745, DOI 10.1137/110827119
- Zhen Guan, John Lowengrub, and Cheng Wang, Convergence analysis for second-order accurate schemes for the periodic nonlocal Allen-Cahn and Cahn-Hilliard equations, Math. Methods Appl. Sci. 40 (2017), no. 18, 6836–6863. MR 3742100, DOI 10.1002/mma.4497
- Zhen Guan, Cheng Wang, and Steven M. Wise, A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation, Numer. Math. 128 (2014), no. 2, 377–406. MR 3259494, DOI 10.1007/s00211-014-0608-2
- Daozhi Han and Xiaoming Wang, A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation, J. Comput. Phys. 290 (2015), 139–156. MR 3324579, DOI 10.1016/j.jcp.2015.02.046
- Z. Hu, S. M. Wise, C. Wang, and J. S. Lowengrub, Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation, J. Comput. Phys. 228 (2009), no. 15, 5323–5339. MR 2541456, DOI 10.1016/j.jcp.2009.04.020
- Hyeong-Gi Lee, J. S. Lowengrub, and J. Goodman, Modeling pinchoff and reconnection in a Hele-Shaw cell. I. The models and their calibration, Phys. Fluids 14 (2002), no. 2, 492–513. MR 1878351, DOI 10.1063/1.1425843
- X. Li, Z. Qiao, and C. Wang, Double stabilizations and convergence analysis of a second- order linear numerical scheme for the nonlocal Cahn-Hilliard equation, Sci. China Math., 2022. Accepted and in press.
- Xiao Li, Zhonghua Qiao, and Cheng Wang, Stabilization parameter analysis of a second-order linear numerical scheme for the nonlocal Cahn-Hilliard equation, IMA J. Numer. Anal. 43 (2023), no. 2, 1089–1114. MR 4568441, DOI 10.1093/imanum/drab109
- Chun Liu, Cheng Wang, and Yiwei Wang, A structure-preserving, operator splitting scheme for reaction-diffusion equations with detailed balance, J. Comput. Phys. 436 (2021), Paper No. 110253, 22. MR 4234219, DOI 10.1016/j.jcp.2021.110253
- Chun Liu, Cheng Wang, Steven M. Wise, Xingye Yue, and Shenggao Zhou, A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system, Math. Comp. 90 (2021), no. 331, 2071–2106. MR 4280293, DOI 10.1090/mcom/3642
- Chun Liu, Cheng Wang, Steven M. Wise, Xingye Yue, and Shenggao Zhou, An iteration solver for the Poisson-Nernst-Planck system and its convergence analysis, J. Comput. Appl. Math. 406 (2022), Paper No. 114017, 13. MR 4356837, DOI 10.1016/j.cam.2021.114017
- Yuan Liu, Wenbin Chen, Cheng Wang, and Steven M. Wise, Error analysis of a mixed finite element method for a Cahn-Hilliard-Hele-Shaw system, Numer. Math. 135 (2017), no. 3, 679–709. MR 3606459, DOI 10.1007/s00211-016-0813-2
- Yiran Qian, Cheng Wang, and Shenggao Zhou, A positive and energy stable numerical scheme for the Poisson-Nernst-Planck-Cahn-Hilliard equations with steric interactions, J. Comput. Phys. 426 (2021), Paper No. 109908, 17. MR 4195851, DOI 10.1016/j.jcp.2020.109908
- Jie Shen, Cheng Wang, Xiaoming Wang, and Steven M. Wise, Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: application to thin film epitaxy, SIAM J. Numer. Anal. 50 (2012), no. 1, 105–125. MR 2888306, DOI 10.1137/110822839
- Jie Shen and Xiaofeng Yang, Decoupled, energy stable schemes for phase-field models of two-phase incompressible flows, SIAM J. Numer. Anal. 53 (2015), no. 1, 279–296. MR 3303685, DOI 10.1137/140971154
- C. Wang and S. M. Wise, An energy stable and convergent finite-difference scheme for the modified phase field crystal equation, SIAM J. Numer. Anal. 49 (2011), no. 3, 945–969. MR 2802554, DOI 10.1137/090752675
- S. M. Wise, Unconditionally stable finite difference, nonlinear multigrid simulation of the Cahn-Hilliard-Hele-Shaw system of equations, J. Sci. Comput. 44 (2010), no. 1, 38–68. MR 2647498, DOI 10.1007/s10915-010-9363-4
- S. M. Wise, C. Wang, and J. S. Lowengrub, An energy-stable and convergent finite-difference scheme for the phase field crystal equation, SIAM J. Numer. Anal. 47 (2009), no. 3, 2269–2288. MR 2519603, DOI 10.1137/080738143
- Xiaofeng Yang and Jia Zhao, On linear and unconditionally energy stable algorithms for variable mobility Cahn-Hilliard type equation with logarithmic Flory-Huggins potential, Commun. Comput. Phys. 25 (2019), no. 3, 703–728. MR 3878519, DOI 10.4208/cicp.oa-2017-0259
- Maoqin Yuan, Wenbin Chen, Cheng Wang, Steven M. Wise, and Zhengru Zhang, An energy stable finite element scheme for the three-component Cahn-Hilliard-type model for macromolecular microsphere composite hydrogels, J. Sci. Comput. 87 (2021), no. 3, Paper No. 78, 30. MR 4251690, DOI 10.1007/s10915-021-01508-w
- Juan Zhang, Cheng Wang, Steven M. Wise, and Zhengru Zhang, Structure-preserving, energy stable numerical schemes for a liquid thin film coarsening model, SIAM J. Sci. Comput. 43 (2021), no. 2, A1248–A1272. MR 4241519, DOI 10.1137/20M1375656
Bibliographic Information
- Yunzhuo Guo
- Affiliation: School of Mathematical Sciences, Beijing Normal University, Beijing 100875, People’s Republic of China
- ORCID: 0009-0003-0502-6832
- Email: yunzguo@mail.bnu.edu.cn
- Cheng Wang
- Affiliation: Department of Mathematics, The University of Massachusetts, North Dartmouth, Massachusetts 02747
- MR Author ID: 652762
- Email: cwang1@umassd.edu
- Steven M. Wise
- Affiliation: Department of Mathematics, The University of Tennessee, Knoxville, Tennessee 37996
- MR Author ID: 615795
- ORCID: 0000-0003-3824-2075
- Email: swise1@utk.edu
- Zhengru Zhang
- Affiliation: Laboratory of Mathematics and Complex Systems, Beijing Normal University, Beijing 100875, People’s Republic of China
- ORCID: 0000-0001-6753-0653
- Email: zrzhang@bnu.edu.cn
- Received by editor(s): March 20, 2023
- Received by editor(s) in revised form: August 15, 2023, and September 17, 2023
- Published electronically: November 20, 2023
- Additional Notes: The second author was partially supported by the NSF DMS-2012269, DMS-2309548. The third author was partially supported by the NSF DMS-2012634, DMS-2309547. The fourth author was partially supported by the NSFC No. 11871105 and Science Challenge Project No. TZ2018002.
The second author is the corresponding author. - © Copyright 2023 American Mathematical Society
- Journal: Math. Comp. 93 (2024), 2185-2214
- MSC (2020): Primary 35K35, 35K55, 49J40, 65M06, 65M12
- DOI: https://doi.org/10.1090/mcom/3916
- MathSciNet review: 4759373