Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
158 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
45 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

A novel high-order linearly implicit and energy-stable additive Runge-Kutta methods for gradient flow models (2307.03905v1)

Published 8 Jul 2023 in math.NA and cs.NA

Abstract: This paper introduces a novel paradigm for constructing linearly implicit and high-order unconditionally energy-stable schemes for general gradient flows, utilizing the scalar auxiliary variable (SAV) approach and the additive Runge-Kutta (ARK) methods. We provide a rigorous proof of energy stability, unique solvability, and convergence. The proposed schemes generalizes some recently developed high-order, energy-stable schemes and address their shortcomings. On the one other hand, the proposed schemes can incorporate existing SAV-RK type methods after judiciously selecting the Butcher tables of ARK methods \cite{sav_li,sav_nlsw}. The order of a SAV-RKPC method can thus be confirmed theoretically by the order conditions of the corresponding ARK method. Several new schemes are constructed based on our framework, which perform to be more stable than existing SAV-RK type methods. On the other hand, the proposed schemes do not limit to a specific form of the nonlinear part of the free energy and can achieve high order with fewer intermediate stages compared to the convex splitting ARK methods \cite{csrk}. Numerical experiments demonstrate stability and efficiency of proposed schemes.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (63)
  1. Energy-decaying extrapolated RK-SAV methods for the Allen–Cahn and Cahn–Hilliard equations. SIAM J. Sci. Comput., 41(6):A3703–A3727, 2019.
  2. A review on phase-field models of brittle fracture and a new fast hybrid formulation. Comput. Mech., 55:383–405, 2015.
  3. K. Burrage. Efficiently implementable algebraically stable Runge–Kutta methods. SIAM J. Numer. Anal., 19(2):245–258, 1982.
  4. K. Burrage and J. C. Butcher. Stability criteria for implicit Runge–Kutta methods. SIAM J. Numer. Anal., 16(1):46–57, 1979.
  5. A microscopic theory for domain wall motion and its experimental verification in Fe-Al alloy domain growth kinetics. J. Phys. Colloq, 38:C7–51–C7–54, 1977.
  6. Free energy of a nonuniform system. I. Interficial free energy. J. Chem. Phys., 28:258–267, 1958.
  7. A novel second-order scheme for the molecular beam epitaxy model with slope selection. Commun. Comput. Phys., 25(4):1024–1044, 2019.
  8. A third order exponential time differencing numerical scheme for No-Slope-Selection epitaxial thin film model with energy stability. J. Sci. Comput., 81:154–185, 2019.
  9. A new Lagrange multiplier approach for gradient flows. Comput. Methods Appl. Mech. Engrg., 367:113030, 2020.
  10. Generalized SAV approaches for gradient systems. J. Comput. Appl. Math., 394:113532, 2021.
  11. Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model. J. Comput. Phys., 341:44–60, 2017.
  12. Fast and stable explicit operator splitting methods for phase-field models. J. Comput. Phys., 303:45–65, 2015.
  13. S. Clarke and D. Vvedensky. Origin of reflection high-energy electron-diffraction intensity oscillations during molecular-beam epitaxy: A computational modeling approach . Phys. Rev. Lett, 58:2235–2238, 1987.
  14. Maximum principle preserving exponential time differencing schemes for the nonlocal Allen-Cahn equation. SIAM J. Numer. Anal., 57:875–898, 2019.
  15. Maximum bound principles for a class of semilinear parabolic equations and exponential time-differencing schemes. SIAM Rev., 63:317–359, 2021.
  16. Analysis of fully discrete approximations for dissipative systems and application to time-dependent nonlocal diffusion problems. J. Sci. Comput., 78(3):1438–1466, 2019.
  17. C. Elliot and A. Stuart. The global dynamics of discrete semilinear parabolic equations. SIAM J. Numer. Anal., 30:1622–1663, 1993.
  18. M. Elsey and B. Wirth. A simple and efficient scheme for phase field crystal simulation. ESIAM: M2AN, 47:1413–1432, 2013.
  19. D. J. Eyre. Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation. Mater. Res. Soc. Sympos. Proc., 529:39–46, 1998.
  20. Stabilized crank-nicolson/adams-bashforth schemes for phase field models. East Asian J. Appl. Math., 3(1):59–80, 2013.
  21. D. Furihata. A stable and conservative finite difference scheme for the Cahn-Hilliard equation, 2001.
  22. D. Furihata and T. Matsuo. a. Chapman and Hall/CRC, 1st edition, 2010.
  23. Supplementary variable method for thermodynamically consistent partial differential equations. Comput. Methods Appl. Mech. Engrg., 381:113746, 2021.
  24. A conservative Fourier pseudo-spectral method for the nonlinear Schrödinger equation. J. Comput. Phys., 328:354–370, 2017.
  25. Arbitrarily high-order linear energy stable schemes for gradient flow models. J. Comput. Phys., 419:109610, 2020.
  26. F. Guillén-González and G. Tierra. On linear schemes for a Cahn–Hilliard diffuse interface model. J. Comput. Phys., 234:140–171, 2013.
  27. F. Guo and W. Dai. Arbitrarily high-order accurate and energy-stable schemes for solving the conservative Allen–Cahn equation. Numer. Methods Partial Differential Eq., 39:187–212, 2022.
  28. Z. Guo and P. Lin. A thermodynamically consistent phase-field model for two-phase flows with thermocapillary effects. J. Fluid Mech., 766:226–271, 2015.
  29. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations. Springer-Verlag, Berlin, 2nd edition, 2006.
  30. A new second-order maximum-principle preserving finite difference scheme for Allen-Cahn equations with periodic boundary conditions. Appl. Math. Lett., 104:106256, 2020.
  31. F. Huang and J. Shen. Stability and error analysis of a class of high-order IMEX scheme for Navier–Stokes equations with periodic boundary conditions. SIAM J. Numer. Anal., 59:2926–2954, 2021.
  32. J. Huang. Energy stable schemes for gradient flows based on the DVD method. arXiv:2210.11960v1, 2022.
  33. T. Ide. Some energy preserving finite element schemes based on the discrete variational derivative method. Appl. Math. Comput., 175:277–296, 2006.
  34. Improving the accuracy and consistency of the scalar auxiliary variable (SAV) method with relaxation. J. Comput. Phys., 456:110954, 2022.
  35. Stabilized Exponential-SAV cchemes preserving energy dissipation law and maximum bound principle for the Allen–Cahn type equations. J Sci Comput, 92(2):66, 2022.
  36. Maximum bound principle preserving integrating factor Runge-Kutta methods for semilinear parabolic equations. J. Comput. Phys., 439(110405):18, 2021.
  37. Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection. Math. Comp., 87(312):1859–1885, 2017.
  38. B. Li and J.-G. Liu. Thin film epitaxy with or without slope selection. Eur. J. Appl. Math., 14:713–743, 2003.
  39. B. Li and J.-G. Liu. Stability analysis of large time-stepping methods for epitaxial growth models. SIAM J. Numer. Anal., 44:1759–1779, 2006.
  40. D. Li and W. Sun. Linearly implicit and high-order energy-conserving schemes for nonlinear wave equations. J. Sci. Comput., 83:17, 2020.
  41. Linear high-order energy-preserving schemes for the nonlinear Schrödinger equation with wave operator using the scalar auxiliary variable approach. J. Sci. Comput., 88:25, 2021.
  42. New SAV-pressure correction methods for the Navier-Stokes equations: stability and error analysis. Math. Comp., 92(340):141–167, 2022.
  43. On energy stable, maximum-principle preserving, second-order BDF scheme with variable steps for the Allen-Cahn equation. SIAM J. Numer. Anal., 58:2294–2314, 2020.
  44. Z. Liu and X. Li. The Exponential Scalar Auxiliary Variable (E-SAV) Approach for Phase Field Models and Its Explicit Computing. SIAM J. Sci. Comput., 42:B630–B655, 2020.
  45. Z. Liu and X. Li. A novel lagrange multiplier approach with relaxation for gradient flows. arXiv:2210.02723v1 [math.NA], 2022.
  46. Local structure-preserving algorithms for the molecular beam epitaxy model with slope selection. Am. J. Math, 26:4745–4765, 2021.
  47. Margination of white blood cells: a computational approach by a hydrodynamic phase field model. J. Fluid Mech., 790:389–406, 2016.
  48. Error analysis of a mixed finite element method for the molecular beam epitaxy model. SIAM J. Numer. Anal., 53:184–205, 2015.
  49. A. Sandu and M. Günther. A generalized–structure approach to additive Runge–Kutta methods. SIAM J. Numer. Anal., 53(1):17–42, 2015.
  50. The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys., 353:407–416, 2018.
  51. A new class of efficient and robust energy stable schemes for gradient flows. SIAM Rev., 61:474–506, 2019.
  52. Unconditionally stable methods for gradient flow using convex splitting Runge-Kutta scheme. J. Comput. Phys., 347:367–381, 2017.
  53. Z. Tan and H. Tang. A general class of linear unconditionally energy stable schemes for the gradient flows. J. Comput. Phys., 464(111372):32, 2022.
  54. T. Tang and J. Yang. Implicit-explicit scheme for the Allen-Cahn equation preserves the maximum principle. J. Comput. Math., 34:471–481, 2016.
  55. Simulating binary fluid-surfactant dynamics by a phase field model . Discrete Contin. Dyn. Syst. Ser. B, 4(17):1289–1307, 2012.
  56. Phase-field model for isothermal phase transitions in binary alloys. Phys. Rev. A, 45:7424–7438, 1992.
  57. Efficient IMEX and consistently energy-stable methods of diffuse-interface models for incompressible three-component flows. Comput. Phys. Commun., 282, 2023.
  58. X. Yang. X. Yang, Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends. J. Comput. Phys., 327:294–316, 2016.
  59. X. Yang and L. Ju. Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model. Comput. Methods Appl. Mech. Eng., 135:691–712, 2017.
  60. Numerical approximations for a three components Cahn–Hilliard phase-field model based on the invariant energy quadratization method. Math. Models Methods Appl. Sci., 27(11):1993–2030, 2017.
  61. On the preserving of the maximum principle and energy stability of high-order implicit-explicit Runge-Kutta schemes for the space-fractional Allen-Cahn equation. Numer. Algorithms, 88(3):1309–1336, 2021.
  62. Y. Zhang and J. Shen. A generalized SAV approach with relaxation for dissipative systems. J. Comput. Phys., 464:111311, 2022.
  63. N. Zheng and X. Li. Error analysis of the SAV Fourier-spectral method for the Cahn-Hilliard-Hele-Shaw system. Adv. Comput. Math., 47(71), 2021.

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now