The Allen-Cahn (AC) model is a mathematical equation that represents the phase separation process. The AC equation has numerous applications in various disciplines, such as image processing, physics, and biology. It models phase transitions, such as solidification and grain growth in materials, pattern formation in chemical reactions, and domain coarsening in biological systems like lipid membranes. Numerical methods are crucial for solving the AC equation due to its complexity and nonlinear nature. Analytical solutions are often extremely difficult to obtain. Therefore, the development of efficient numerical techniques is indispensable for approximating solutions and studying phase transitions, material behavior, and pattern formation accurately. We investigate the stability of an explicit finite difference method (FDM) used to numerically solve the two-dimensional (2D) AC model with a high-order polynomial potential, which was recently proposed to preserve a more intricate structure of interfaces. To demonstrate the precision and optimal estimate of our stability constraints, we conduct various computational tests using the derived time step formulas that ensure the maximum principle.
Citation: Jaeyong Choi, Seokjun Ham, Soobin Kwak, Youngjin Hwang, Junseok Kim. Stability analysis of an explicit numerical scheme for the Allen-Cahn equation with high-order polynomial potentials[J]. AIMS Mathematics, 2024, 9(7): 19332-19344. doi: 10.3934/math.2024941
The Allen-Cahn (AC) model is a mathematical equation that represents the phase separation process. The AC equation has numerous applications in various disciplines, such as image processing, physics, and biology. It models phase transitions, such as solidification and grain growth in materials, pattern formation in chemical reactions, and domain coarsening in biological systems like lipid membranes. Numerical methods are crucial for solving the AC equation due to its complexity and nonlinear nature. Analytical solutions are often extremely difficult to obtain. Therefore, the development of efficient numerical techniques is indispensable for approximating solutions and studying phase transitions, material behavior, and pattern formation accurately. We investigate the stability of an explicit finite difference method (FDM) used to numerically solve the two-dimensional (2D) AC model with a high-order polynomial potential, which was recently proposed to preserve a more intricate structure of interfaces. To demonstrate the precision and optimal estimate of our stability constraints, we conduct various computational tests using the derived time step formulas that ensure the maximum principle.
[1] | C. Lee, H. Kim, S. Yoon, S. Kim, D. Lee, J. Park, et al., An unconditionally stable scheme for the Allen-Cahn equation with high-order polynomial free energy, Commun. Nonlinear Sci., 95 (2021), 105658. https://doi.org/10.1016/j.cnsns.2020.105658 doi: 10.1016/j.cnsns.2020.105658 |
[2] | J. Shin, J. Yang, C. Lee, J. Kim, The Navier-Stokes-Cahn-Hilliard model with a high-order polynomial free energy, Acta Mech., 231 (2020), 2425–2437. https://doi.org/10.1007/s00707-020-02666-y doi: 10.1007/s00707-020-02666-y |
[3] | S. Lee, S. Yoon, J. Kim, A linear convex splitting scheme for the Cahn-Hilliard equation with a high-order polynomial free energy, Int. J. Numer. Meth. Eng., 124 (2023), 3586–3602. https://doi.org/10.1002/nme.7288 doi: 10.1002/nme.7288 |
[4] | S. M. Allen, J. W. Cahn, A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metall., 27 (1979), 1085–1095. https://doi.org/10.1016/0001-6160(79)90196-2 doi: 10.1016/0001-6160(79)90196-2 |
[5] | Y. Li, S. Lan, X. Liu, B. Lu, L. Wang, An efficient volume repairing method by using a modified Allen-Cahn equation, Pattern Recogn., 107 (2020), 107478. https://doi.org/10.1016/j.patcog.2020.107478 doi: 10.1016/j.patcog.2020.107478 |
[6] | T. K. Akinfe, A. C. Loyinmi, An improved differential transform scheme implementation on the generalized Allen-Cahn equation governing oil pollution dynamics in oceanography, Part. Differ. Eq. Appl. Math., 6 (2022), 100416. https://doi.org/10.1016/j.padiff.2022.100416 doi: 10.1016/j.padiff.2022.100416 |
[7] | H. Kim, C. Lee, S. Kwak, Y. Hwang, S. Kim, Y. Choi, et al., Three-dimensional volume reconstruction from multi-slice data using a shape transformation, Comput. Math. Appl., 113 (2022), 52–58. https://doi.org/10.1016/j.camwa.2022.03.018 doi: 10.1016/j.camwa.2022.03.018 |
[8] | Y. Li, X. Song, S. Kwak, J. Kim, Weighted 3D volume reconstruction from series of slice data using a modified Allen-Cahn equation, Pattern Recogn., 132 (2022), 108914. https://doi.org/10.1016/j.patcog.2022.108914 doi: 10.1016/j.patcog.2022.108914 |
[9] | Z. Han, H. Xu, J. Wang, A simple shape transformation method based on phase-field model, Comput. Math. Appl., 147 (2023), 121–129. https://doi.org/10.1016/j.camwa.2023.07.020 doi: 10.1016/j.camwa.2023.07.020 |
[10] | B. Xia, R. Yu, X. Song, X. Zhang, J. Kim, An efficient data assimilation algorithm using the Allen-Cahn equation, Eng. Anal. Bound. Elem., 155 (2023), 511–517. https://doi.org/10.1016/j.enganabound.2023.06.029 doi: 10.1016/j.enganabound.2023.06.029 |
[11] | C. Liu, Z. Qiao, Q. Zhang, Multi-phase image segmentation by the Allen-Cahn Chan-Vese model, Comput. Math. Appl., 141 (2023), 207–220. https://doi.org/10.1016/j.camwa.2022.12.020 doi: 10.1016/j.camwa.2022.12.020 |
[12] | J. Gao, B. Song, Z. Mao, Structural topology optimization through implicit boundary evolution based on the Allen-Cahn equation, Eng. Optimiz., 53 (2021), 125–144. https://doi.org/10.1080/0305215X.2019.1705288 doi: 10.1080/0305215X.2019.1705288 |
[13] | X. Xiao, X. Feng, A second-order maximum bound principle preserving operator splitting method for the Allen-Cahn equation with applications in multi-phase systems, Math. Comput. Simulat., 202 (2022), 36–58. https://doi.org/10.1016/j.matcom.2022.05.024 doi: 10.1016/j.matcom.2022.05.024 |
[14] | D. Lee, The numerical solutions for the energy-dissipative and mass-conservative Allen-Cahn equation, Comput. Math. Appl., 80 (2020), 263–284. https://doi.org/10.1016/j.camwa.2020.04.007 doi: 10.1016/j.camwa.2020.04.007 |
[15] | H. Zhang, J. Yan, X. Qian, X. Chen, S. Song, Explicit third-order unconditionally structure-preserving schemes for conservative Allen-Cahn equations, J. Sci. Comput., 90 (2022), 1–29. |
[16] | H. Zhang, J. Yan, X. Qian, S. Song, Up to fourth-order unconditionally structure-preserving parametric single-step methods for semilinear parabolic equations, Comput. Method. Appl. M., 393 (2022), 114817. |
[17] | J. Feng, Y. Zhou, T. Hou, A maximum-principle preserving and unconditionally energy-stable linear second-order finite difference scheme for Allen-Cahn equations, Appl. Math. Lett., 118 (2021), 107179. https://doi.org/10.1016/j.aml.2021.107179 doi: 10.1016/j.aml.2021.107179 |
[18] | H. Zhang, J. Yan, X. Qian, S. Song, Numerical analysis and applications of explicit high order maximum principle preserving integrating factor Runge-Kutta schemes for Allen-Cahn equation, Appl. Numer. Math., 161 (2021), 372–390. https://doi.org/10.1016/j.apnum.2020.11.022 doi: 10.1016/j.apnum.2020.11.022 |
[19] | D. Lee, Computing the area-minimizing surface by the Allen-Cahn equation with the fixed boundary, AIMS Math., 8 (2023), 23352–23371. https://doi.org/10.3934/math.20231187 doi: 10.3934/math.20231187 |
[20] | Y. Deng, Z. Weng, Barycentric interpolation collocation method based on Crank-Nicolson scheme for the Allen-Cahn equation, AIMS Math., 6 (2021), 3857–3873. https://doi.org/10.3934/math.2021229 doi: 10.3934/math.2021229 |
[21] | J. Park, C. Lee, Y. Choi, H. G. Lee, S. Kwak, Y. Hwang, et al., An unconditionally stable splitting method for the Allen-Cahn equation with logarithmic free energy, J. Eng. Math., 132 (2022), 18. https://doi.org/10.1007/s10665-021-10203-6 doi: 10.1007/s10665-021-10203-6 |
[22] | H. Zhang, X. Qian, J. Xia, S. Song, Unconditionally maximum-principle-preserving parametric integrating factor two-step Runge-Kutta schemes for parabolic Sine-Gordon equations, CSIAM T. Appl. Math., 4 (2023), 177–224. |
[23] | M. Liao, D. Wang, C. Zhang, H. Jia, The error analysis for the Cahn-Hilliard phase field model of two-phase incompressible flows with variable density, AIMS Math., 8 (2023), 31158–31185. https://doi.org/10.3934/math.20231595 doi: 10.3934/math.20231595 |
[24] | C. Lee, S. Kim, S. Kwak, Y. Hwang, S. Ham, S. Kang, et al., Semi-automatic fingerprint image restoration algorithm using a partial differential equation, AIMS Math., 8 (2023), 27528–27541. https://doi.org/10.3934/math.20231408 doi: 10.3934/math.20231408 |
[25] | J. Wang, Z. Han, W. Jiang, J. Kim, A fast, efficient, and explicit phase-field model for 3D mesh denoising, Appl. Math. Comput., 458 (2023), 128239. https://doi.org/10.1016/j.amc.2023.128239 doi: 10.1016/j.amc.2023.128239 |
[26] | H. G. Lee, S. Ham, J. Kim, Isotropic finite difference discrettization of Laplacian operator, Appl. Comput. Math., 22 (2023), 259–274. https://doi.org/10.30546/1683-6154.22.2.2023.259 doi: 10.30546/1683-6154.22.2.2023.259 |
[27] | S. Ham, J. Kim, Stability analysis for a maximum principle preserving explicit scheme of the Allen-Cahn equation, Math. Comput. Simulat., 207 (2023), 453–465. https://doi.org/10.1016/j.matcom.2023.01.016 doi: 10.1016/j.matcom.2023.01.016 |
[28] | H. Zhang, X. Qian, J. Xia, S. Song, Efficient inequality-preserving integrators for differential equations satisfying forward Euler conditions, ESAIM-Math. Model. Num., 57 (2023), 1619–1655. https://doi.org/10.1051/m2an/2023029 doi: 10.1051/m2an/2023029 |