Abstract
In this paper, we provide a solution to the Cahn–Hilliard equation using the q-homotopy analysis method (q-HAM). The q-HAM is a more general, simple and widely used method for solving stiff nonlinear partial differential equations. The Cahn–Hilliard equation is a classical model in material sciences that describe spinodal decomposition and phase separation in two-phase flows. Using the q-HAM, the effect of various parameters of physical interest such as diffusive parameter, thickness parameter, advection and reaction terms on concentration is studied. The comparison of the computed solution with the exact solution is presented for some fixed parameter values to validate the solution obtained using the q-HAM.
1. Introduction
Many natural and physical phenomena can be described by reaction, diffusion and advection to describe a variety of physical and chemical processes arising in various fields of sciences and engineering, such as the flow of heat, wave propagation on shallow water, movement of electricity in conductors, the flow of fluids, plasma physics, chemical reactions and quantum mechanical systems. The mathematical modelling of such problems leads to nonlinear partial differential equations (PDEs). These models play a vital role in studying the physical behaviour of such types of natural phenomena. Numerous mathematical models and their solutions have been discussed in the literature, such as Wazwaz–Benjamin–Bona–Mahony equation, Schrödinger equation, diffusion equation, geophysical Korteweg–de Vries equation, Regularized Long Wave equations, Ablowitz-Kaup-Newell-Segur water wave dynamical equation, etc. [Citation1–12]. These models introduce new ideas about the relations of diffusion, reaction, advection and nonlinearity.
The Cahn–Hilliard (CH) equation was introduced by American mathematicians and scientists Cahn and Hilliard [Citation13]. It is an important mathematical model of mathematical physics that describes phase separation processes like spinodal decomposition in multi-phase systems. The nonlinear nature of the CH equation makes it difficult to find its exact solution. It becomes more complicated when it contains diffusion parameters, advection parameters, interface thickness parameters and reaction terms. An essential and important characteristic of the CH equation is the interface thickness between two phases that has a finite thickness. In the literature, researchers and scientists explore and analyse various forms of the CH equation with its applications in different fields of engineering and sciences [Citation14–22]. Bouhassoun and Cherif [Citation23] analyse the fractional form of the CH equation using the homotopy perturbation method. Tripathi et al. [Citation24] discuss the solution of the time-fractional CH equation with reaction term using the homotopy analysis method. Zhu et al. [Citation25] reported the solution of the CH equation with variable mobility. Barrett et al. [Citation26] studied the solution behaviour of the CH equation with degenerate mobility. Berti and Bochiccchio [Citation27] discussed generalized CH equation including thermal and mixing effects. Colby et al. [Citation28] used the adaptive neural networks approach for solving the Allen Cahn and the CH equations numerically. A. Shah and A. A Siddiqui [Citation29] used the variational iteration method to solve the viscous CH equation. S. Hussain and A. Shah [Citation30] applied the homotopy perturbation method and variational iteration method to solve the CH equation. It is known that the exact analytical solution of the CH equation is difficult to obtain due to the sharp interface width and jump discontinuity at the interface. The literature review reveals that the CH equation with parameters of physical interest is not studied extensively due to its complicated nature. Therefore, we propose a series solution to the CH equation using the q-HAM. We also explore the effect of different parameters on concentration. To validate the computed results, we compare the solution with the existing one [Citation31] for some fixed value of the parameters. The graphical illustrations of the results are also given.
2. Mathematical model
The generalized form of the CH equation [Citation25–27] is (1) (1) In Equation (1), is the diffusion parameter, is the thickness of the transition region, is the coefficient of the advection term and is the coefficient of the reaction term.
3. The q-HAM
In 2012, Tawil and Huseen [Citation32] introduced the modified version of the Homotopy Analysis Method known as q-HAM for solving several nonlinear PDEs. They proved that HAM is a special case of q-HAM [Citation33]. However, q-HAM is fast convergent with a large convergence region including other advantages over HAM. It also provides a more appropriate approach to tackle the region of convergence.
3.1. Basic idea of q-HAM
To understand the key idea of q-HAM, consider the nonlinear equation of the following form (2) (2) where is a nonlinear operator, “” is a known function and “” is an unknown function with x and t as independent variables. Let us make the zeroth-order deformation equation as follows (3) (3) where is a non-zero auxiliary parameter, is the embedding parameter, is the auxiliary function, is the initial solution of and is a linear operator. When and and , Equation (3) gives, and , respectively.
Now expanding with respect to using the Taylor series, we have (4) (4) where (5) (5) If one chooses the linear operator , auxiliary parameter , auxiliary function and initial approximation properly, the series in Equation (4) converges, then we get (6) (6) Define the vector (7) (7) then order deformation equation obtained using Equation (3) is given by (8) (8) where (9) (9) and (10) (10) Applying on both sides of Equation (8) and after simplification, we get (11) (11) Using Equation (11) one can approximate the analytical solution of the following form (12) (12) Due to the presence of a factor , q-HAM converges faster than HAM.
4. Numerical experiments
In this section, we solve the CH equation with the given initial conditions.
4.1. Example
Consider the CH equation of the following form [Citation24] (13) (13) with the initial condition (14) (14) By considering the initial condition as an initial guess, we have (15) (15) Using Equations (13) and (14) in Equations (9) and (11), we get the first approximation as (16) (16)
Similarly, the second and third iterations are (17) (17) (18) (18) Continuing in the same line, we can find up to the required number of iterations.
4.1.1. Effect of diffusion on concentration
In , the effect of diffusion parameter “” on concentration has been studied for different values. It is clear from the figure that the concentration increases with the increase of diffusion.
4.1.2. Effect of thickness on concentration
In , the effect of the thickness parameter “” on concentration is depicted by taking various values of the thickness parameter. From , one can conclude that the thickness parameter is directly proportional to the concentration.
4.1.3. Effect of advection on concentration
In , the effect of the advection parameter “” on concentration is represented by taking various values of the advection parameter. From we can observe that the concentration upsurges with the increase of advection.
4.1.4. Effect of reaction term on concentration
In , the effect of the reaction parameter “” on concentration is validated by taking various values of “k”. It is observed that the concentration enhances in the presence of the source term i.e. . Whereas, the concentration decreases in the presence of sink term i.e. , while a very small increase in concentration is observed in the absence of source and sink term i.e. .
4.1.5. Convergence analysis
In order to calculate the effective region of “h”, the h-curves are plotted for different values of independent variables “x” and “t”. From Figures , it can be observed that q-HAM provides a large convergence region as compared to the classical HAM, as shown in .
From Figures , it can be observed that the convergence region of parameter “h” increases with the increase of “n”. states the convergence region of parameter “h” with the increase of “n”.
4.2. Numerical comparison
In this section, we compare computed results with the results obtained in [Citation31] through absolute error analysis taking in Equation (13). (19) (19) The exact solution of Equation (19) is given as [Citation31] (20) (20) To implement q-HAM, consider as an initial guess (21) (21) Using Equations (20) and (21) in Equations (9) and (11), we get the first approximation as (22) (22) Similarly, the second and third iterations are (23) (23) (24) (24) Proceeding in the same manner we can calculate .
4.2.1. Convergence analysis
The valid region for convergence parameter “h” is studied for different values of “n”. It is observed that the convergence region increases with the increase in “n”. The results are elaborated in .
4.2.2. Error analysis
In , the absolute error is calculated for different values of “x” and “t”. The table illustrates the convergence of the approximate solution obtained using the q-HAM.
5. Conclusion
In this work, the effect of various parameters of practical interest on the concentration of the multi-phase system described by the CH equation is investigated. The effect of each parameter is illustrated in Figures . The q-HAM is used to obtain the series solution of the CH equation. The convergence of the h-curve is also studied through Figures and for increasing values of “n”. gives the range of solution controlling parameter “h” for higher “n”. It is clear from Tables and that the q-HAM provides a large range of “h” as compared to HAM. It also suggests that HAM is a special case of q-HAM when n = 1. The comparison of exact and approximate solutions is also provided using the absolute error. validates the application of q-HAM for solving the CH equation involving various parameters of physical importance. The following key observations are made in the light of obtained solutions.
It is observed that concentration increases with the increase of diffusion.
Concentration increases the increase of thickness parameter.
Concentration is directly proportional to the advection parameter.
The results of the reaction term reveal that concentration upsurges in the presence of source term i.e. . Concentration decreases in the presence of sink term i.e. . While a very small increase in concentration is observed in the absence of source and sink term i.e. .
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Akram U, Aly RS, Rizvi STR, et al. Traveling wave solutions for the fractional Wazwaz–Benjamin–Bona–Mahony model in arising shallow water waves. Results Phys. 2021;20:103725. doi:10.1016/j.rinp.2020.103725.
- Rizvi ST, Seadawy AR, Ashraf F, et al. Lump and interaction solutions of a geophysical Korteweg–de Vries equation. Results Phys. 2020;19:Article 103661. doi:10.1016/j.rinp.2020.103661.
- Younas U, Younis M, Aly RS, et al. Diverse exact solutions for modified nonlinear Schrödinger equation with conformable fractional derivative. Results Phys. 2021;20(8):Article 103766. doi:10.1016/j.rinp.2020.103766.
- Younis M, Ali S, Rizvi STR, et al. Investigation of solitons and mixed lump wave solutions with (3 + 1)-dimensional potential-YTSF equation. Commun Nonlinear Sci Numer Simul. 2021;94:Article 105544. doi:10.1016/j.cnsns.2020.105544.
- Ali A, Seadawy AR, Lu D. Soliton solutions of the nonlinear Schrödinger equation with the dual power law nonlinearity and resonant nonlinear schrödinger equation and their modulation instability analysis. Optik (Stuttg). 2017;145:79–88. doi:10.1016/j.ijleo.2017.07.016.
- Aly RS, Younis M, Baber MZ, et al. Diverse acoustic wave propagation to confirmable time–space fractional KP equation arising in dusty plasma. Commun Theor Phys. 2021;73(11):Article 115004.
- Younis M, Nadia C, Mahmood AS, et al. A variety of exact solutions to (2 + 1)-dimensional schrödinger equation. Waves Random Complex Media. 2018;30:1–10. doi:10.1080/17455030.2018.1532131.
- Asghar A, Seadawy AR, Dianchen L. Computational methods and traveling wave solutions for the fourth-order nonlinear Ablowitz-Kaup-Newell-Segur water wave dynamical equation via two methods and its applications. Open Physics. 2018;16(1):219–226. DOI: 10.1515/phys-2018-0032.
- Dianchen L, Seadawy AR, Asghar A. Applications of exact traveling wave solutions of modified Liouville and the symmetric regularized long wave equations via two new techniques. Results Phys. 2018;9:1403–1410. doi:10.1016/j.rinp.2018.04.039.
- Asghar A, Seadawy AR, Dianchen L. Dispersive solitary wave soliton solutions of (2 + 1)-dimensional Boussinesq dynamical equation via extended simple equation method. J King Saud Univ Sci. 2019;31(4):653–658. doi:10.1016/j.jksus.2017.12.015.
- Alruwaili AD, Seadawy AR, Asghar A, et al. Novel analytical approach for the space-time fractional (2 + 1)-dimensional breaking soliton equation via mathematical methods. Mathematics. 2021;9:3253. doi:10.3390/math9243253.
- Alruwaili AD, Seadawy AR, Ali A, et al. Soliton solutions of Calogero– Degasperis–Fokas dynamical equation via modified mathematical methods. Open Physics. 2022;20:174–187. DOI: 10.1515/phys-2022-0016.
- Cahn JW, Hilliard J. Free energy of a nonuniform system. I. interfacial free energy. J Chem Phys. 1958;28:258–267.
- Bertozzi AL, Esedoglu S, Gillette A. Inpainting of binary images using the Cahn-Hilliard equation. IEEE Trans Image Process. 2007;16(1):285–291. doi:10.1109/TIP.2006.887728.
- Junseok K, Seunggyu L, Yongho C, et al. Basic principles and practical applications of the Cahn–Hilliard equation. Math Probl Eng. 2016: 11. Article ID 9532608. doi:10.1155/2016/9532608.
- Miranville A, Piétrus A, Rakotoson J. Dynamical aspect of a generalized Cahn–Hilliard equation based on a microforce balance. Asymptotic Anal. 1998;16:315–345. Corpus ID:118719560.
- Nicolaenko B, Scheurer B, Temam R. Some global dynamical properties of a class of pattern formation equations. Commun Partial Differ Equ. 1989;14(2):245–297. doi:10.1080/03605308908820597.
- Novick-Cohen A, Segel LA. Nonlinear aspects of the Cahn-Hilliard equation. Physica D. 1984;10:277–298. doi:10.1016/0167-2789(84)90180-5.
- Gurtin M. Generalised Ginzburg-Landau and Cahn-Hilliard equations based on a microforce balance. Physica D. 1996;92:178–192. doi:10.1016/0167-2789(95)00173-5.
- Choksi R, Peletier MA, Williams JF. On the phase diagram for microphase separation of diblock copolymers: an approach via a nonlocal Cahn–Hilliard functional. SIAM J Appl Math. 2009;69:1712–1738. Doi:10.1137/080728809.
- Tang P, Qiu F, Zhang H, et al. Phase separation patterns for diblock copolymers on spherical surfaces: a finite volume method. Phys Rev E Stat Nonlin Soft Matter Phys. 2005;72:016710. doi:10.1103/PhysRevE.72.016710.
- Hu SY, Chen LQ. A phase-field model for evolving microstructures with strong elastic inhomogeneity. Acta Mater. 2001;49:1879–1890. doi:10.1016/S1359-6454(01)00118-5.
- Bouhassoun A, Cherif H. Homotopy perturbation method For solving The fractional Cahn-Hilliard equation. J Interdiscip Math. 2015;18:513–524. doi:10.1080/10288457.2013.867627.
- Tripathi NK, Das S, Ong SH, et al. Solution of time-fractional Cahn–Hilliard equation with reaction term using homotopy analysis method. Adv Mech Eng. 2017;9:168781401774077. doi:10.1177/1687814017740773.
- Zhu J, Long-Qing C, Shen J, et al. Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1999;60:3564–3572. doi:10.1103/PhysRevE.60.3564.
- Barrett J, Blowey J, Garcke H. Finite element approximation of the Cahn-Hilliard equation with degenerate mobility. SIAM J Numer Anal. 2000;37. doi:10.1137/S0036142997331669.
- Berti A, Bochicchio I. A mathematical model for phase separation: a generalized Cahn- Hilliard equation. Math Method Appl Sci. 2011;34:1193–1203. doi:10.1002/mma.1432.
- Wight CL, Zhao J. Solving Allen-Cahn and Cahn-Hilliard equations using the adaptive physics informed neural networks. Commun Comput Phys. 2021;29:930–954. doi:10.4208/cicp.OA-2020-0086.
- Shah A, Siddiqui AA. Variational iteration method for the solution of viscous Cahn-Hilliard equation. World Appl Sci J. 2010;11(7):813–818. Corpus ID: 16027042.
- Hussain S, Shah A. An analysis of two iterative techniques for solution of Cahn-Hilliard equation. Int J Nonlinear Sci. 2011;12(1):42–47.
- Ugurlu Y, Kaya D. Solutions of the Cahn–Hilliard equation. Comput Math Appl. 2008;56:3038–3045. doi:10.1016/j.camwa.2008.07.007.
- El-Tawil MA, Huseen S. The q-homotopy analysis method (q-HAM). Int J Appl Math Mech. 2012;8:51–75. Corpus ID:124341343.
- Liao SJ. Beyond perturbation: introduction to the homotopy analysis method. Boca Raton (FL): CRC Press, Chapman & Hall; 2003.