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

An exact divergence-free spectral method for incompressible and resistive magneto-hydrodynamic equations in two and three dimensions (2312.12218v1)

Published 19 Dec 2023 in math.NA and cs.NA

Abstract: In this paper, we present exact divergence-free spectral method for solving the incompressible and resistive magneto-hydrodynamic (MHD) equations in two and three dimensions, as well as the efficient solution algorithm and unconditionally energy-stable fully-discretized numerical schemes. We introduce new ideas of constructing two families of exact divergence-free vectorial spectral basis functions on domains diffeomorphic to squares or cubes. These bases are obtained with the help of orthogonality and derivative relation of generalised Jacobi polynomials, several de Rham complexes, as well as the property of contravariant Piola transformation. They are well-suited for discretizing the velocity and magnetic fields, respectively, thereby ensuring point-wise preservation of the incompressibility condition and the magnetic Gauss's law. With the aid of these bases, we propose a family of exact divergence-free implicit-explicit $k$-step backward differentiation formula (DF-BDF-$k$) fully-discretized schemes for the MHD system. These schemes naturally decouple the pressure field from the velocity field. Consequently, the stability of the space-time fully-discretized numerical schemes based on these bases are significantly enhanced. These schemes exhibit unconditional stability for $k=1,2$, and demonstrate exceptional stability and accuracy for $k=3,4$, verified with extensive numerical results for long time simulations using large time step sizes. Furthermore, we present efficient solution algorithms for these two decoupled equations for the velocity and magnetic fields, respectively, by exploiting the sparsity and structure of the resultant linear algebraic systems. Ample numerical examples in two and three dimensions are provided to demonstrate the distinctive accuracy, efficiency and stability of our proposed method.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (39)
  1. R. Adams and J. Fournier. Sobolev Spaces. Elsevier, 2003.
  2. Finite element exterior calculus: from Hodge theory to numerical stability. Bull. Am. Math. Soc., 47(2):281–354, 2010.
  3. Mixed Finite Element Methods And Applications, volume 44. Springer, 2013.
  4. J. Boland and R. Nicolaides. Stability of finite elements under divergence constraints. SIAM J. Numer. Anal., 20(4):722–731, 1983.
  5. The effect of nonzero ∇⋅𝑩⋅∇𝑩\nabla\cdot\bm{B}∇ ⋅ bold_italic_B on the numerical solution of the magnetohydrodynamic equations. Journal of Computational Physics, 35(3):426–430, 1980.
  6. Divergence-free 𝑯⁢(div)𝑯div\bm{H}({\rm div})bold_italic_H ( roman_div )-conforming hierarchical bases for magnetohydrodynamics (MHD). Commun. Math. Stat., 1(1):19–35, 2013.
  7. H. Choi and J. Shen. Efficient splitting schemes for magneto-hydrodynamic equations. Sci. China Math., 59(8):1495–1510, 2016.
  8. P. Ciarlet. Three-Dimensional Elasticity. Elsevier, 1988.
  9. P. Concus and G. Golub. Use of fast direct methods for the efficient numerical solution of nonseparable elliptic equations. SIAM J. Numer. Anal., 10(6):1103–1120, 1973.
  10. A new approximate block factorization preconditioner for two-dimensional incompressible (reduced) resistive MHD. SIAM J. Sci. Comput., 35(3):B701–B730, 2013.
  11. C. Evans and J. Hawley. Simulation of magnetohydrodynamic flows: a constrained transport method. Astrophys. J., 332:659–677, 1988.
  12. R. Falk and M. Neilan. Stokes complexes and the construction of stable finite elements with pointwise mass conservation. SIAM J. Numer. Anal., 51(2):1308–1326, 2013.
  13. E. Gawlik and F. Gay-Balmaz. A finite element method for MHD that preserves energy, cross-helicity, magnetic helicity, incompressibility, and div𝑩𝑩\bm{B}bold_italic_B= 0. J. Comput. Phys., 450:110847, 2022.
  14. V. Girault and P. Raviart. Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms, volume 5. Springer Science & Business Media, 2012.
  15. G. Golub and C. Van Loan. Matrix Computations. JHU Press, 2013.
  16. B. Guo and Y. Jiao. Spectral method for Navier–Stokes equations with non-slip boundary conditions by using divergence-free base functions. J. Sci. Comput., 66(3):1077–1101, 2016.
  17. Generalized Jacobi polynomials/functions and their applications. Appl. Numer. Math., 59(5):1011–1028, 2009.
  18. J. Guzmán and M. Neilan. Conforming and divergence-free Stokes elements on general triangular meshes. Math. Comput., 83(285):15–36, 2014.
  19. A fully divergence-free finite element method for magnetohydrodynamic equations. Math. Models Methods Appl. Sci., 28(4):659–695, 2018.
  20. Stable finite element methods preserving ∇⋅𝑩=0⋅∇𝑩0\nabla\cdot\bm{B}=0∇ ⋅ bold_italic_B = 0 exactly for MHD models. Numer. Math., 135(2):371–396, 2017.
  21. K. Hu and J. Xu. Structure-preserving finite element methods for stationary MHD models. Math. Comput., 88(316):553–581, 2019.
  22. S. Jardin. Computational Methods in Plasma Physics. CRC Press, 2010.
  23. C. Kelley. Iterative Methods for Linear and Nonlinear Equations. SIAM, 1995.
  24. Structure-preserving and helicity-conserving finite element approximations and preconditioning for the Hall MHD equations. J. Comput. Phys., 492:112410, 2023.
  25. F. Li and C. Shu. Locally divergence-free discontinuous Galerkin methods for MHD equations. J. Sci. Comput., 22(1-3):413–442, 2005.
  26. A robust finite element method for Darcy–Stokes flow. SIAM J. Numer. Anal., 40(5):1605–1631, 2002.
  27. P. Monk. Finite Element Methods for Maxwell’s Equations. Oxford University Press, 2003.
  28. A discrete divergence free weak Galerkin finite element method for the Stokes equations. Appl. Numer. Math., 125(1):172–182, 2018.
  29. M. Neilan. Discrete and conforming smooth de Rham complexes in three dimensions. Math. Comput., 84(295):2059–2081, 2015.
  30. A highly efficient and accurate divergence-free spectral method for curl-curl equation in two and three dimensions. arXiv:2308.12865, 2023.
  31. Y. Saad. Iterative Methods for Sparse Linear Systems. SIAM, 2003.
  32. Towards a scalable fully-implicit fully-coupled resistive mhd formulation with stabilized fe methods. J. Comput. Phys., 229(20):7649–7671, 2010.
  33. J. Shen. Efficient spectral-Galerkin method I. direct solvers for second- and fourth-order equations by using Legendre polynomials. SIAM J. Sci. Comput., 15(6):1489–1505, 1994.
  34. Spectral Methods: Algorithms, Analysis and Applications, volume 41. Springer Science & Business Media, 2011.
  35. J. Stone and M. Norman. ZEUS-2D: a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. II. the magnetohydrodynamic algorithms and tests. Astrophys. J., Suppl. Ser., 80:791, 1992.
  36. X. Tai and R. Winther. A discrete de Rham complex with enhanced smoothness. Calcolo, 43(4):287–306, 2006.
  37. L. Trefethen and D. Bau. Numerical Linear Algebra, volume 181. SIAM, 2022.
  38. F. White. Fluid Mechanics (7th Edition). McGraw-Hill, 2011.
  39. X. Ye and C. Hall. A discrete divergence-free basis for finite element methods. Numer. Algorithms, 16(3):365–380, 1997.

Summary

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

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