Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
156 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 structure-preserving semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach and Alfvén numbers (2403.04517v1)

Published 7 Mar 2024 in math.NA and cs.NA

Abstract: We present a divergence-free semi-implicit finite volume scheme for the simulation of the ideal magnetohydrodynamics (MHD) equations which is stable for large time steps controlled by the local transport speed at all Mach and Alfv\'en numbers. An operator splitting technique allows to treat the convective terms explicitly while the hydrodynamic pressure and the magnetic field contributions are integrated implicitly, yielding two decoupled linear implicit systems. The linearity of the implicit part is achieved by means of a semi-implicit time linearization. This structure is favorable as second-order accuracy in time can be achieved relying on the class of semi-implicit IMplicit-EXplicit Runge-Kutta (IMEX-RK) methods. In space, implicit cell-centered finite difference operators are designed to discretely preserve the divergence-free property of the magnetic field on three-dimensional Cartesian meshes. The new scheme is also particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, since no explicit numerical dissipation is added to the implicit contribution and the time step is scale independent. Likewise, highly magnetized flows can benefit from the implicit treatment of the magnetic fluxes, hence improving the computational efficiency of the novel method. The convective terms undergo a shock-capturing second order finite volume discretization to guarantee the effectiveness of the proposed method even for high Mach number flows. The new scheme is benchmarked against a series of test cases for the ideal MHD equations addressing different acoustic and Alfv\'en Mach number regimes where the performance and the stability of the new scheme is assessed.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (62)
  1. A local time-stepping discontinuous galerkin algorithm for the mhd system. In ESAIM: proceedings, volume 28, pages 33–54, 2009.
  2. Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations. Applied Numerical Mathematics, 25:151–167, 1982.
  3. Linearly implicit all Mach number shock capturing schemes for the Euler equations. Journal of Computational Physics, 393:278–312, 2019.
  4. D. Balsara. Second-Order Accurate Schemes for Magnetohydrodynamics with Divergence-Free Reconstruction. The Astrophysical Journal Supplement Series, 151:149–184, 2004.
  5. D. S. Balsara. Multidimensional hlle riemann solver: application to euler and magnetohydrodynamic flows. Journal of Computational Physics, 229(6):1970–1993, 2010.
  6. D. S. Balsara and M. Dumbser. Divergence-free mhd on unstructured meshes using high order finite volume schemes based on multidimensional riemann solvers. Journal of Computational Physics, 299:687–715, 2015.
  7. A staggered mesh algorithm using high order godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. Journal of Computational Physics, 149(2):270–292, 1999.
  8. W. Barsukow. Truly multi-dimensional all-speed schemes for the euler equations on cartesian grids. Journal of Computational Physics, 435:110216, 2021.
  9. W. Barsukow. All-speed numerical methods for the euler equations via a sequential explicit time integration. Journal of Scientific Computing, 95(2):53, 2023.
  10. A well-balanced semi-implicit imex finite volume scheme for ideal magnetohydrodynamics at all mach numbers. Journal of Scientific Computing, 98(2):34, 2024.
  11. A low mach two-speed relaxation scheme for the compressible euler equations with gravity. Communications in Mathematical Sciences, 21, 2023.
  12. High Order Semi-implicit Schemes for Time Dependent Partial Differential Equations. Journal of Scientific Computing, 68:975–1001, 2016.
  13. S. Boscarino and L. Pareschi. On the asymptotic properties of IMEX Runge-Kutta schemes for hyperbolic balance laws. Journal of Computational and Applied Mathematics, 316:60–73, 2017.
  14. A unified IMEX Runge-Kutta approach for hyperbolic systems with multiscale relaxation. SIAM Journal of Numerical Analysis, 55(4):2085–2109, 2017.
  15. A high order semi-implicit IMEX WENO scheme for the all-Mach isentropic Euler system. Journal Of Computational Physics, 392:594–618, 2019.
  16. S. Boscarino and G. Russo. On a class of uniformly accurate IMEX Runge-Kutta schemes and applications to hyperbolic systems with relaxation. SIAM Journal of Scientific Computing, 31:1926–1945, 2009.
  17. Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes. Journal of Computational Physics, 486:112130, 2023.
  18. High-order ader-weno ale schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics. International Journal for Numerical Methods in Fluids, 76(10):737–778, 2014.
  19. High order cell-centered lagrangian-type finite volume schemes with time-accurate local time stepping on unstructured triangular meshes. Journal of Computational Physics, 291:120–150, 2015.
  20. W. Boscheri and L. Pareschi. High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers. Journal Of Computational Physics, 434:110206, 2021.
  21. W. Boscheri and M. Tavelli. High order semi-implicit schemes for viscous compressible flows in 3d. Applied Mathematics and Computation, 434:127457, 2022.
  22. A multiwave approximate riemann solver for ideal mhd based on relaxation. i: theoretical framework. Numerische Mathematik, 108:7–42, 2007.
  23. The effect of nonzero ∇⋅B⋅∇𝐵\nabla\cdot B∇ ⋅ italic_B on the numerical solution of the magnetohydrodynamic equations. Journal of Computational Physics, 35(3):426–430, 1980.
  24. L. Brugnano and V. Casulli. Iterative solution of piecewise linear systems. SIAM Journal of Scientific Computing, 30:463–472, 2007.
  25. High Order Asymptotic Preserving Finite Difference Weno Schemes with Constrained Transport for Mhd Equations in All Sonic Mach Numbers. Journal Of Computational Physics, 488:112240, 2023.
  26. An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations. Journal of Computational Physics, 231:5685–5704, 2012.
  27. Evolution of the orszag–tang vortex system in a compressible medium. i. initial average subsonic flow. Physics of Fluids B: Plasma Physics, 1(11):2153–2171, 1989.
  28. Hyperbolic divergence cleaning for the mhd equations. Journal of Computational Physics, 175(2):645–673, 2002.
  29. S. Dellacherie. Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number. Journal of Computational Physics, 229:978–1016, 2010.
  30. M. Dumbser and D. S. Balsara. A new efficient formulation of the hllem riemann solver for general conservative and non-conservative hyperbolic systems. Journal of Computational Physics, 304:275–319, 2016.
  31. A divergence-free semi-implicit finite volume scheme for ideal, viscous, and resistive magnetohydrodynamics. International Journal for Numerical Methods in Fluids, 89:16–42, 2019.
  32. M. Dumbser and V. Casulli. A conservative, weakly nonlinear semi-implicit finite volume scheme for the compressible Navier-Stokes equations with general equation of state. Applied Mathematics and Computation, 272:479–497, 2016.
  33. M. Dumbser and E. F. Toro. On universal osher-type schemes for general nonlinear hyperbolic conservation laws. Communications in Computational Physics, 10(3):635–671, 2011.
  34. S. A. E. G. Falle. Rarefaction Shocks, Shock Errors, and Low Order of Accuracy in ZEUS. The Astrophysical Journal, 577(2):L123–L126, October 2002.
  35. S. A. E. G. Falle and S. S. Komissarov. On the inadmissibility of non-evolutionary shocks. Journal of Plasma Physics, 65(1):29–58, 2001.
  36. F. Fambri. A novel structure preserving semi-implicit finite volume method for viscous and resistive magnetohydrodynamics. International Journal for Numerical Methods in Fluids, 93:3447–3489, 2021.
  37. A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume / finite element scheme for the incompressible mhd equations. Journal of Computational Physics, 493:112493, 2023.
  38. M. Fey and M. Torrilhon. A constrained transport upwind scheme for divergence-free advection. In Hyperbolic Problems: Theory, Numerics, Applications: Proceedings of the Ninth International Conference on Hyperbolic Problems held in CalTech, Pasadena, March 25–29, 2002, pages 529–538, 2003.
  39. An unsplit godunov method for ideal mhd via constrained transport. Journal of Computational Physics, 205(2):509–539, 2005.
  40. S. Godunov. Finite difference methods for the computation of discontinuous solutions of the equations of fluid dynamics. Mat. Sb., 47:271–306, 1959.
  41. H. Guillard and A. Murrone. On the behavior of upwind schemes in the low Mach number limit : II. Godunov type schemes. Computers & Fluids, 33:655–675, 2004.
  42. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev., 25:35–61, 1983.
  43. An unstaggered constrained transport method for the 3d ideal magnetohydrodynamic equations. Journal of Computational Physics, 230(10):3803–3829, 2011.
  44. A high-order unstaggered constrained-transport method for the three-dimensional ideal magnetohydrodynamic equations based on the method of lines. SIAM Journal on Scientific Computing, 35(2):A623–A651, 2013.
  45. J.M. Hyman and M. Shashkov. Natural discretizations for the divergence, gradient, and curl on logically rectangular grids. Computers and Mathematics with Applications, 33:81–104, 1997.
  46. R. Klein. Semi-implicit extension of a Godunov-type scheme based on low Mach number asymptotics I: One-dimensional flow. J. Comp. Phys., 121:213–237, 1995.
  47. A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors. Astronomy and Astrophysics, 668:A143, 2022.
  48. Mimetic finite difference method. Journal of Computational Physics, 257:1163–1227, 2014.
  49. An implicit-explicit solver for a two-fluid single-temperature model. Journal of Computational Physics, 498:112696, 2024.
  50. A discrete operator calculus for finite difference approximations. Computer Methods in Applied Mechanics and Engineering, 187:365–383, 2000.
  51. Divergence correction techniques for maxwell solvers based on a hyperbolic model. Journal of Computational Physics, 161(2):484–511, 2000.
  52. S. Osher and F. Solomon. A partially implicit method for large stiff systems of Ode’s with only few equations introducing small time-constants. SIAM Journal of Numerical Analysis, 13:645–663, 1976.
  53. L. Pareschi and G. Russo. Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation. Journal of Scientific Computing, 25:129–155, 2005.
  54. Evolution of the orszag–tang vortex system in a compressible medium. ii. supersonic flow. Physics of Fluids B: Plasma Physics, 3(1):29–44, 1991.
  55. Notes on the eigensystem of magnetohydrodynamics. SIAM Journal on Applied Mathematics, 56(1):57–67, 1996.
  56. C-M. Tang S. A. Orszag. Small-scale structure of two-dimensional magnetohydrodynamic turbulence. Journal of Fluid Mechanics, 90(1):129–143, 1979.
  57. Y. Saad and M. Schultz. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput., 7:856–869, 1986.
  58. E.F. Toro. Riemann Solvers and Numerical Methods for Fluid Dynamics: a Practical Introduction. Springer, 2009.
  59. Flux splitting schemes for the Euler equations. Computers & Fluids, 70:1–12, 2012.
  60. M. Torrilhon. Non-uniform convergence of finite volume schemes for riemann problems of ideal magnetohydrodynamics. Journal of Computational Physics, 192(1):73–94, 2003.
  61. Towards a new generation of multi-dimensional stellar evolution models: development of an implicit hydrodynamic code. Astronomy and Astrophysics, 531:A86, 2011.
  62. K. Yee. Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media. IEEE Transactions on antennas and propagation, 14(3):302–307, 1966.
Citations (1)
View on Semantic Scholar

Summary

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

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