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


[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 Schrodinger Equations, American Mathematical Society, 2003.
[6] J. Ciavaldini, Analyse numerique dun probleme de stefan adeux phases par une methode delements
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, Uniersity 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 and A. Mesforush, Semi-algebraic mode analysis for multigrid method on regular rect-
angular 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 decompo-
sition, 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.