Bifurcation Analysis Reveals Solution Structures of Phase Field Models

doi:10.1007/s42967-022-00221-1

Abstract

Abstract: The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems. Here, we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models. Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions. To elucidate the idea, we apply this analytical approach to three representative phase field equations: the Allen-Cahn equation, the Cahn-Hilliard equation, and the Allen-Cahn-Ohta-Kawasaki system. The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

Key words: Phase field modeling, Bifurcations, Multiple solutions

Xinyue Evelyn Zhao, Long-Qing Chen, Wenrui Hao, Yanxiang Zhao. Bifurcation Analysis Reveals Solution Structures of Phase Field Models[J]. Communications on Applied Mathematics and Computation, 2024, 6(1): 64-89.

TrendMD

References

[1] Akrivis, G., Crouzeix, M., Makridakis, C.:Implicit-explicit multistep finite element methods for nonlinear parabolic problems. Math. Comput. Am. Math. Soc. 67, 457 (1998)
[2] Anderson, D.M., McFadden, G.B., Wheeler, A.A.:Diffuse-interface methods in fluid mechanics. Annu. Rev. Fluid Mech. 30, 139 (1998)
[3] Bats, F.S., Fredrickson, G.H.:Block copolymers-designer soft materials. Phys. Today 52(2), 32 (1999)
[4] Bauer, L., Keller, H.B., Reiss, E.:Multiple eigenvalues lead to secondary bifurcation. SIAM Rev. 17(1), 101-122 (1975)
[5] Benesova, B., Melcher, C., Suli, E.:An implicit midpoint spectral approximation of nonlocal Cahn-Hilliard equations. SIAM J. Numer. Anal. 52, 1466 (2014)
[6] Boettinger, W.J., Warren, J.A., Beckermann, C., Karma, A.:Phase-field simulation of solidification. Annu. Rev. Mater. Res. 32, 163-194 (2002)
[7] Borden, M.J., Verhoosel, C.V., Scott, M.A., Hughes, T.J., Landis, C.M.:A phase-field description of dynamic brittle fracture. Comput. Methods Appl. Mech. Eng. 217, 77 (2012)
[8] Brower, R., Kessler, D., Koplik, J., Levine, H.:Geometrical models of interface evolution. Phy. Rev. A 29, 1335 (1984)
[9] Cahn, J., Allen, S.:A microscopic theory for domain wall motion and its experimental verification in Fe-Al alloy domain growth kinetics. J. de Physique 38, C7-51 (1977)
[10] Cahn, J., Hilliard, J.:Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys. 28, 258 (1958)
[11] Cahn, J.W.:On spinodal decomposition. Acta Metallurgica 9(9), 795-801 (1961)
[12] Camley, B., Zhao, Y., Li, B., Levine, H., Rappel, W.-J.:Crawling and turning in a minimal reaction-diffusion cell motility model:coupling cell shape and biochemistry. Phys. Rev. E 95, 012401 (2017)
[13] Chan, H., Nejad, M., Wei, J.:Lamellar phase solutions for diblock copolymers with nonlocal diffusions. Phys. D 388, 22-32 (2019)
[14] Chen, L.Q.:Phase-field models for microstructure evolution. Annu. Rev. Mater. Res. 32, 113 (2002)
[15] Chen, L.Q., Shen, J.:Applications of semi-implicit Fourier-spectral method to phase field equations. Comput. Phys. Commun. 108, 147 (1998)
[16] Chen, L.Q., Zhao, Y.H.:From classical thermodynamics to phase-field method. Progress in Materials Science 124, 10086 (2022)
[17] Cheng, W., Yang, X., Shen, J.:Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model. J. Comput. Phys. 341, 44 (2017)
[18] Choi, H., Zhao, Y.:Second-order stabilized semi-implicit energy stable schemes for bubble assemblies in binary and ternary systems. DCDS-B 27(8), 4649-4683 (2021)
[19] Choksi, R.:Scaling laws in microphase separation of diblock copolymers. J. Nonlinear Sci. 11, 223-236 (2011)
[20] Choksi, R.:On global minimizers for a variational problem with long-range interactions. Quart. Appl. Math. 70, 517-537 (2012)
[21] Cox, S.M., Matthews, P.C.:Exponential time differencing for stiff systems. J. Comput. Phys. 176, 430 (2002)
[22] Crandall, M.G., Rabinowitz, P.H.:Bifurcation from simple eigenvalues. J. Funct. Anal. 8(2), 321-340 (1971)
[23] Du, Q., Feng, X.:The phase field method for geometric moving interfaces and their numerical approximations. Handbook of Numerical Analysis 21, 425 (2020)
[24] Du, Q., Ju, L., Li, X., Qiao, Z.:Maximum bound principles for a class of semilinear parabolic equations and exponential time-differencing schemes. SIAM Rev. 63, 317 (2021)
[25] Du, Q., Nicolaides, R.A.:Numerical analysis of a continuum model of phase transition. SIAM J. Numer. Anal. 28, 1310 (1991)
[26] Eyre, D. J.:Unconditionally gradient stable time marching the Cahn-Hilliard equation. In:Computational and Mathematical Models of Microstructural Evolution (San Francisco, CA, 1998), Mater. Res. Soc. Sympos. Proc., vol. 529, pp. 39-46. MRS, Warrendale (1998)
[27] Fix, G.J.:Phase field model for free boundary problems. In:Fasano, A., Primicerio, M. (eds) Free Boundary Problems:Theory and Applications, vol. 2. Pitman Publishing Inc., Massachusetts, USA (1983)
[28] Friedman, A., Hu, B.:Bifurcation from stability to instability for a free boundary problem arising in a tumor model. Arch. Ration. Mech. Anal. 180(2), 293-330 (2006)
[29] Friedman, A., Hu, B.:Bifurcation for a free boundary problem modeling tumor growth by Stokes equation. SIAM J. Math. Anal. 39(1), 174-194 (2007)
[30] Friedman, A., Reitich, F.:Symmetry-breaking bifurcation of analytic solutions to free boundary problems:an application to a model of tumor growth. Trans. Amer. Math. Soc. 353(4), 1587-1634 (2001)
[31] Gennip, Y., Peletier, M.:Copolymer-homopolymer blends:global energy minimisation and global energy bounds. Calc. Var. Partial Differ. Equ. 33, 75-111 (2008)
[32] Ginzburg, V.Q., Landau, L.E.:On the theory of superconductivity. Soviet Physics-JETP 20(12), 1064-1082 (1950)
[33] Goldstein, G.R., Miranville, A., Schimperna, G.:A Cahn-Hilliard model in a domain with non-permeable walls. Phys. D 240(8), 754-766 (2011)
[34] Haber, R., Unbehauen, H.:Structure identification of nonlinear dynamic systems survey on input/output approaches. Automatica 26(4), 651-677 (1990)
[35] Hao, W.:An adaptive homotopy tracking algorithm for solving nonlinear parametric systems with applications in nonlinear ODEs. Appl. Math. Lett. 125, 107767 (2022)
[36] Hao, W., Crouser, E., Friedman, A.:Mathematical model of sarcoidosis. Proc. Natl. Acad. Sci. USA 111(45), 16065-16070 (2014)
[37] Hao, W., Friedman, A.:The LDL-HDL profile determines the risk of atherosclerosis:a mathematical model. PloS One 9(3), e90497 (2014)
[38] Hao, W., Hauenstein, J., Hu, B., Liu, Y., Sommese, A., Zhang, Y.-T.:Multiple stable steady states of a reaction-diffusion model on zebrafish dorsal-ventral patterning. Discrete Contin. Dyn. Syst. Ser. S 4(6), 1413-1428 (2011)
[39] Hao, W., Hauenstein, J., Hu, B., Liu, Y., Sommese, A.J., Zhang, Y.-T.:Bifurcation for a free boundary problem modeling the growth of a tumor with a necrotic core. Nonlinear Analysis:Real World Applications 13(2), 694-709 (2012)
[40] Hao, W., Hauenstein, J., Hu, B., Sommese, A.:A three-dimensional steady-state tumor system. Appl. Math. Comput. 218(6), 2661-2669 (2011)
[41] Hao, W., Hauenstein, J., Shu, C.-W., Sommese, A., Xu, Z., Zhang, Y.-T.:A homotopy method based on WENO schemes for solving steady state problems of hyperbolic conservation laws. J. Comput. Phys. 250, 332-346 (2013)
[42] Hao, W., Nepomechie, R., Sommese, A.:Completeness of solutions of Bethe's equations. Phys. Rev. E 88(5), 052113 (2013)
[43] Hao, W., Nepomechie, R., Sommese, A.:Singular solutions, repeated roots and completeness for higher-spin chains. J. Stat. Mech. Theory Exp. 2014(3), P03024 (2014)
[44] Hao, W., Xue, C.:Spatial pattern formation in reaction-diffusion models:a computational approach. J. Math. Biol. 80, 521-543 (2020)
[45] Hao, W., Zheng, C.:An adaptive homotopy method for computing bifurcations of nonlinear parametric systems. J. Sci. Comput. 82(3), 53 (2020)
[46] Hao, W., Zheng, C.:Bifurcation analysis of a free boundary model of the atherosclerotic plaque formation associated with the cholesterol ratio. Chaos 30(9), 093113 (2020)
[47] Hao, W., Zheng, C.:A stochastic homotopy tracking algorithm for parametric systems of nonlinear equations. J. Sci. Comput. 87(3), 87 (2021)
[48] Joo, S., Xu, X., Zhao, Y.:Analysis and computation for Allen-Cahn-Ohta-Nakazawa model in ternary system. Interfaces Free Bound. 23, 535-559 (2021)
[49] Kessler, D., Koplik, J., Levine, H.:Geometrical models of interface evolution. II. Numerical simulation. Phys. Rev. A 30, 3161 (1984)
[50] Kessler, D., Koplik, J., Levine, H.:Geometrical models of interface evolution. III. Theory of dendritic growth. Phys. Rev. A 31, 1712 (1985)
[51] Knopf, P., Lam, K.F., Liu, C., Metzger, S.:Phase-field dynamics with transfer of materials:the Cahn-Hilliard equation with reaction rate dependent dynamic boundary conditions. ESAIM:Mathematical Modelling and Numerical Analysis 55(1), 229-282 (2021)
[52] Liu, C., Shen, J.:A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method. Phys. D 179(3/4), 211-228 (2003)
[53] Liu, C., Wu, H.:An energetic variational approach for the Cahn-Hilliard equation with dynamic boundary condition:model derivation and mathematical analysis. Arch. Ration. Mech. Anal. 233(1), 167-247 (2019)
[54] Morgan, A., Sommese, A.:Computing all solutions to polynomial systems using homotopy continuation. Appl. Math. Comput. 24(2), 115-138 (1987)
[55] Nakazawa, H., Ohta, T.:Microphase separation of ABC-type triblock copolymers. Macromolecules 26, 5503-5511 (1993)
[56] Ohta, T., Kawasaki, K.:Equilibrium morphology of block copolymer melts. Macromolecules 19, 2621-2632 (1986)
[57] Ren, X., Truskinovsky, L.:Finite scale microstructures in nonlocal elasticity. J. Elasticity 59, 319-355 (2000)
[58] Ren, X., Wei, J.:On the multiplicity of solutions of two nonlocal variational problems. SIAM J. Math. Anal. 4, 909-924 (2000)
[59] Ren, X., Wei, J.:Triblock copolymer theory:ordered ABC lamellar phase. J. Nonlinear Sci. 13, 175-208 (2003)
[60] Ren, X., Wei, J.:A double bubble in a ternary system with inhibitory long range interaction. Arch. Ration. Mech. Anal. 208, 201-253 (2013)
[61] Ren, X., Wei, J.:A double bubble assembly as a new phase of a ternary inhibitory system. Arch. Ration. Mech. Anal. 215, 967-1034 (2015)
[62] Rheinboldt, W.:Numerical methods for a class of finite dimensional bifurcation problems. SIAM J. Numer. Anal. 15(1), 1-11 (1978)
[63] Rheinboldt, W.:Numerical analysis of continuation methods for nonlinear structural problems. Comput. & Struct. 13(1), 103-113 (1981)
[64] Rowlinson, J.S.:Translation of J. D. van der Waals' "The Thermodynamic Theory of Capillarity Under the Hypothesis of a Continuous Variation of Density". J. Stat. Phys. 20(2), 197-200 (1979)
[65] Shen, J., Xu, J., Yang, J.:The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys. 353, 407 (2018)
[66] Song, H., Shu, C.-W.:Unconditional energy stability analysis of a second order implicit-explicit local discontinuous Galerkin method for the Cahn-Hilliard equation. J. Sci. Comput. 73, 1178-1203 (2017)
[67] Steinbach, I.:Phase-field model for microstructure evolution at the mesoscopic scale. Annu. Rev. Mater. Res. 43, 89-107 (2013)
[68] Stephen, M.J., Suhl, H.:Weak time dependence in pure superconductors. Phys. Rev. Lett. 13(26), 797-800 (1964)
[69] Van der Waals, J.D.:Theorie thermodynamique de la capillarite, dans l'hypothese d'une variation continue de la densite. Archives Neerlandaises des Sciences Exactes et Naturelles XXVIII, 121-209 (1979)
[70] Wang, C., Ren, X., Zhao, Y.:Bubble assemblies in ternary systems with long range interaction. Comm. Math. Sci. 17, 2309-2324 (2019)
[71] Xu, J., Li, Y., Wu, S., Bousquetd, A.:On the stability and accuracy of partially and fully implicit schemes for phase field modeling. Comput. Methods Appl. Mech. Eng. 345, 826 (2019)
[72] Xu, X., Zhao, Y.:Energy stable semi-implicit schemes for Allen-Cahn-Ohta-Kawasaki model in binary system. J. Sci. Comput. 80, 1656-1680 (2019)
[73] Xu, X., Zhao, Y.:Maximum principle preserving schemes for binary systems with long-range interactions. J. Sci. Comput. 84, 33 (2020)
[74] Xu, Z., Du, Q.:On the ternary Ohta-Kawasaki free energy and its one-dimensional global minimizers. arXiv:2111.09877 (2021)
[75] Yang, X.:Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends. J. Comput. Phys. 327, 294 (2016)
[76] Zhao, X.Y., Hu, B.:Symmetry-breaking bifurcation for a free-boundary tumor model with time delay. J. Differential Equations 269(3), 1829-1862 (2020)
[77] Zhao, X.Y., Hu, B.:Bifurcation for a free boundary problem modeling a small arterial plaque. J. Differential Equations 288, 250-287 (2021)
[78] Ziebert, F., Aranson, I.S.:Computational approaches to substrate-based cell motility. npj Computational Materials, 2(1), 103-118 (2016)