Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
194 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

Unique solvability and error analysis of the Lagrange multiplier approach for gradient flows (2405.03415v1)

Published 6 May 2024 in math.NA, cs.NA, and math.AP

Abstract: The unique solvability and error analysis of the original Lagrange multiplier approach proposed in [8] for gradient flows is studied in this paper. We identify a necessary and sufficient condition that must be satisfied for the nonlinear algebraic equation arising from the original Lagrange multiplier approach to admit a unique solution in the neighborhood of its exact solution, and propose a modified Lagrange multiplier approach so that the computation can continue even if the aforementioned condition is not satisfied. Using Cahn-Hilliard equation as an example, we prove rigorously the unique solvability and establish optimal error estimates of a second-order Lagrange multiplier scheme assuming this condition and that the time step is sufficient small. We also present numerical results to demonstrate that the modified Lagrange multiplier approach is much more robust and can use much larger time step than the original Lagrange multiplier approach.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (30)
  1. Scalar auxiliary variable/Lagrange multiplier based pseudospectral schemes for the dynamics of nonlinear Schrödinger/Gross-Pitaevskii equations. J. Comput. Phys., 437:Paper No. 110328, 19, 2021.
  2. A mixed discontinuous Galerkin, convex splitting scheme for a modified Cahn-Hilliard equation and an efficient nonlinear multigrid solver. Discrete Contin. Dyn. Sys. B, 18:2211–2238, 2013.
  3. W. Bao and Q. Du. Computing the ground state solution of bose–einstein condensates by a normalized gradient flow. SIAM J. Sci. Comput., 25(5):1674–1697, March 2004.
  4. Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential. J. Comput. Phys.: X, 3:100031, 2019.
  5. An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation. J. Comput. Appl. Math., 362:574–595, 2019.
  6. A third order accurate in time, BDF-type energy stable scheme for the Cahn-Hilliard equation. Numer. Math. Theor. Meth. Appl., 15(2):279–303, 2022.
  7. A second-order, weakly energy-stable pseudo-spectral scheme for the Cahn-Hilliard equation and its solution by the homogeneous linear iteration method. J. Sci. Comput., 69:1083–1114, 2016.
  8. A new lagrange multiplier approach for gradient flows. Comput. Methods Appl. Mech. Engrg., 367:13070, 2020.
  9. Q. Cheng and J. Shen. Global constraints preserving scalar auxiliary variable schemes for gradient flows. SIAM J. Sci. Comput., 42:A2514–A2536, 2020.
  10. Stability and convergence of a second order mixed finite element method for the Cahn-Hilliard equation. IMA J. Numer. Anal., 36:1867–1897, 2016.
  11. A positivity-preserving, energy stable scheme for a ternary Cahn-Hilliard system with the singular interfacial parameters. J. Comput. Phys., 442:110451, 2021.
  12. Optimal rate convergence analysis of a numerical scheme for the ternary Cahn-Hilliard system with a Flory-Huggins-deGennes energy potential. J. Comput. Appl. Math., 406:114474, 2022.
  13. A positivity-preserving, energy stable and convergent numerical scheme for the Cahn-Hilliard equation with a Flory-Huggins-deGennes energy. Commun. Math. Sci., 17:921–939, 2019.
  14. Q. Du and F. Lin. Numerical approximations of a norm-preserving gradient flow and applications to an optimal partition problem. Nonlinearity, 22(1):67–83, December 2008.
  15. Simulating the deformation of vesicle membranes under elastic bending energy in three dimensions. J. Comput. Phys., 212(2):757–777, March 2006.
  16. The global dynamics of discrete semilinear parabolic equations. SIAM J. Numer. Anal., 30:1622–1663, 1993.
  17. D. Eyre. Unconditionally gradient stable time marching the Cahn-Hilliard equation. In J. W. Bullard, R. Kalia, M. Stoneham, and L.Q. Chen, editors, Computational and Mathematical Models of Microstructural Evolution, volume 53, pages 1686–1712, Warrendale, PA, USA, 1998. Materials Research Society.
  18. An H2superscript𝐻2H^{2}italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT convergence of a second-order convex-splitting, finite difference scheme for the three-dimensional Cahn-Hilliard equation. Commu. Math. Sci., 14:489–515, 2016.
  19. An improved error analysis for a second-order numerical scheme for the Cahn-Hilliard equation. J. Comput. Appl. Math., 388:113300, 2021.
  20. Generalized SAV-exponential integrator schemes for Allen-Cahn type gradient flows. SIAM J. Numer. Anal., 60(4):1905–1931, 2022.
  21. Stabilized exponential-SAV schemes preserving energy dissipation law and maximum bound principle for the Allen-Cahn type equations. J. Sci. Comput., 92(2):Paper No. 66, 34, 2022.
  22. D. Li and Z. Qiao. On second order semi-implicit Fourier spectral methods for 2D Cahn-Hilliard equations. J. Sci. Comput., 70:301–341, 2017.
  23. D. Li and Z. Qiao. On the stabilization size of semi-implicit Fourier-spectral methods for 3D Cahn-Hilliard equations. Commun. Math. Sci., 15:1489–1506, 2017.
  24. Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations. SIAM J. Numer. Anal., 54:1653–1681, 2016.
  25. The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys., 353:407–416, 2018.
  26. A new class of efficient and robust energy stable schemes for gradient flows. SIAM Review, 61(3):474–506, 2019.
  27. J. Shen and X. Yang. Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Sys. A, 28:1669–1691, 2010.
  28. A second-order energy stable BDF numerical scheme for the Cahn-Hilliard equation. Commun. Comput. Phys., 23:572–602, 2018.
  29. Numerical approximations for a three components Cahn-Hilliard phase-field model based on the invariant energy quadratization method. M3AS: Mathematical Models and Methods in Applied Sciences, 27(11), August 2017.
  30. An energy stable finite element scheme for the three-component Cahn-Hilliard-type model for macromolecular microsphere composite hydrogels. J. Sci. Comput., 87:78, 2021.
Citations (1)

Summary

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

X Twitter Logo Streamline Icon: https://streamlinehq.com