Finite element approximation of coupled Cahn-Hilliard equations with a logarithmic potential and nondegenerate mobility

Document Type : Research Article

Authors

1 Department of Mathematics, College of Sciences, University of Basrah, Basrah, Iraq & Department of Mathematics, College of Education for Pure Sciences, University of Thi-Qar, Thi-Qar, Iraq

2 Department of Mathematics, College of Sciences, University of Basrah, Basrah, Iraq

Abstract

This research presents a numerical analysis conducted on a system of coupled Cahn-Hilliard equations featuring a logarithmic potential, nondegenerate mobility, and homogeneous Neumann boundary conditions. These equations are derived from a model describing phase separation in a thin film of binary liquid mixture. The study proposes semi-discrete and fully-discrete piecewise linear finite element approximations to the continuous problem. Existence, uniqueness, and various stability estimates for the approximate solutions are established. Fully-discrete error bounds are derived, and optimal time discretisation error is demonstrated. An iterative method is introduced for solving the resulting nonlinear algebraic system, and linear stability analysis in one space dimension is investigated. The research concludes with numerical experiments, providing illustrations of some of the theoretical findings, conducted in both one and two space dimensions.

Keywords

Main Subjects

References

[1] G.A. Al-Juaifri, A.J. Harfash, Finite element analysis of nonlinear reaction–diffusion system of Fitzhugh–Nagumo type with Robin boundary conditions, Math. Comput. Simul. 203 (2023) 486–517.
[2] J.W. Barrett, J.F. Blowey, An optimal error bound for a finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy, ESAIM: Math. Model. Numer. Anal. 33 (1999) 971–987.
[3] J.W. Barrett, J.F. Blowey, Finite element approximation of an Allen–Cahn/Cahn–Hilliard system, IMA J. Numer. Anal. 22 (2002) 11–71.
[4] J.W. Cahn, J.E. Hilliard, Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys. 28 (1958) 258–267.
[5] T. Cazenave, Semilinear Schrödinger Equations, American Mathematical Society, 2003.
[6] J. Ciavaldini, Analyse numérique d'un problème de Stefan à deux phases par une méthode d'éléments finis, SIAM J. Numer. Anal. 12 (1975) 464–487.
[7] C.M. Elliott, H. Garcke, On the Cahn–Hilliard equation with degenerate mobility, SIAM J. Math. Anal. 27 (1996) 404–423.
[8] C.M. Elliott, S. Luckhaus, A generalised diffusion equation for phase separation of a multi-component mixture with interfacial free energy, MA, University of Minnesota, preprint 887, 1991.
[9] H. Gomez, V.M. Calo, Y. Bazilevs, T.J. Hughes, Isogeometric analysis of the Cahn–Hilliard phase-field model, Comput. Methods Appl. Mech. Eng. 197 (2008) 4333–4352.
[10] J.D. Gunton, R. Toral, A. Chakrabarti, Numerical studies of phase separation in models of binary alloys and polymer blends, Phys. Scripta 33 (1990) 12–19.
[11] N. Habibi, A. Mesforush, Semi-algebraic mode analysis for multigrid method on regular rectangular and triangular grids, J. Math. Model. 11 (2023) 547–572.
[12] S.M. Hassan, A.J. Harfash, Finite element analysis of a two-species chemotaxis system with two chemicals, Appl. Numer. Math. 182 (2022) 148–175.
[13] S.M. Hassan, A.J. Harfash, Finite element approximation of a Keller–Segel model with additional self- and cross-diffusion terms and a logistic source, Commun. Nonlinear Sci. Numer. Simul. 104 (2022) 106063.
[14] Y. He, Y. Liu, T. Tang, On large time-stepping methods for the Cahn–Hilliard equation, Appl. Numer. Math. 57 (2007) 616–628.
[15] M. Izadi, M. Afshar, Solving the Basset equation via Chebyshev collocation and LDG methods, J. Math. Model. 9 (2021) 61–79.
[16] M. Jalili, R. Salehi, The approximate solution of one dimensional stochastic evolution equations by meshless methods, J. Math. Model. 9 (2021) 599–609.
[17] J. Kim, K. Kang, A numerical method for the ternary Cahn–Hilliard system with a degenerate mobility, Appl. Numer. Math. 59 (2009) 1029–1042.
[18] J. Langer, M. Bar-On, H.D. Miller, New computational method in the theory of spinodal decomposition, Phys. Rev. A 11 (1975) 1417.
[19] A. Mesforush, S. Larsson, A posteriori error analysis for the Cahn–Hilliard equation, J. Math. Model. 10 (2022) 437–452.
[20] L. Modica, The gradient theory of phase transitions and the minimal interface criterion, Arch. Rational Mech. Anal. 98 (1987) 123–142.
[21] M. Mohammed, T. Mouhcine, A finite element approximation of a current-induced magnetization dynamics model, J. Math. Model. 10 (2022) 53–69.
[22] D.M. Saylor, C.-S. Kim, D.V. Patwardhan, J.A. Warren, Diffuse-interface theory for structure formation and release behavior in controlled drug release systems, Acta Biomater. 3 (2007) 851–864.
[23] S. Tremaine, On the origin of irregular structure in Saturn's rings, Astrophys. J. 125 (2003) 894.
Journal of Mathematical Modeling
Volume 13, Issue 1
March 2025
Pages 49-68

Files

Share

How to cite

Statistics