Time-Adaptive PIROCK Method with Error Control for Multi-Fluid and Single-Fluid MHD Systems

Abstract

The solar atmosphere is a complex environment with diverse species and varying ionization states, especially in the chromosphere, where significant ionization variations occur. This region transitions from highly collisional to weakly collisional states, leading to complex plasma state transitions influenced by magnetic strengths and collisional properties. These processes introduce numerical stiffness in multi-fluid models, imposing severe timestep restrictions on standard time integration methods. New numerical methods are essential to address these computational challenges, effectively managing the diverse timescales in multi-fluid and multi-physics models. The widely used time operator splitting technique offers a straightforward approach but requires careful timestep management to avoid stability issues and errors. Despite some studies on splitting errors, their impact on solar and stellar astrophysics is often overlooked. We focus on a Multi-Fluid Multi-Species (MFMS) model, which presents significant challenges for time integration. We propose a second-order Partitioned Implicit-Explicit Runge-Kutta (PIROCK) method that combines efficient explicit and implicit integration techniques with variable time-stepping and error control. Compared to a standard third-order explicit method and a first-order Lie splitting approach, the PIROCK method shows robust advantages in accuracy, stability, and computational efficiency. Our results reveal PIROCK's capability to solve multi-fluid problems with unprecedented efficiency. Preliminary results on chemical fractionation represent a significant step toward understanding the First-Ionization-Potential (FIP) effect in the solar atmosphere.

References in corpus (21)

• PLUTO: a Numerical Code for Computational Astrophysics
• The stellar atmosphere simulation code Bifrost
• Extension of the MURaM radiative MHD code for coronal simulations
• The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
• Partially Ionized Plasmas in Astrophysics
• Multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma
• The Multi-slit Approach to Coronal Spectroscopy with the Multi-slit Solar Explorer (MUSE)
• The cool surge following flux emergence in a radiation-MHD experiment
• Surges and Si IV bursts in the solar atmosphere. Understanding IRIS and SST observations through RMHD experiments
• Chromospheric Extension of the MURaM Code
• MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
• Nonequilibrium ionization and ambipolar diffusion in solar magnetic flux emergence processes
• Ambipolar diffusion in the Bifrost code
• Spatio-temporal analysis of chromospheric heating in a plage region
• Analysis of operator splitting in the non-asymptotic regime for nonlinear reaction-diffusion equations. Application to the dynamics of premixed flames
• Adaptive time splitting method for multi-scale evolutionary partial differential equations
• Multi-Parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR
• Multi-Fluid Simulations of Upper Chromospheric Magnetic Reconnection with Helium-Hydrogen mixture
• Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D
• The Yarkovsky effect on the long-term evolution of binary asteroids
• Assessing the capability of a model-based stellar XUV estimation