A fast multi-dimensional magnetohydrodynamic formulation of the transition region adaptive conduction (TRAC) method
arXiv:2106.03989 · doi:10.1051/0004-6361/202140987
Abstract
We have demonstrated that the Transition Region Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, successfully removing the influence of numerical resolution on the coronal density response to impulsive heating. This is achieved by adjusting the parallel thermal conductivity, radiative loss, and heating rates to broaden the transition region (TR), below a global cutoff temperature, so that the steep gradients are spatially resolved even when using coarse numerical grids. Implementing the original 1D formulation of TRAC in multi-dimensional magnetohydrodynamic (MHD) models would require tracing a large number of magnetic field lines at every time step in order to prescribe a global cutoff temperature to each field line. In this paper, we present a highly efficient formulation of the TRAC method for use in multi-dimensional MHD simulations, which does not rely on field line tracing. In the TR, adaptive local cutoff temperatures are used instead of global cutoff temperatures to broaden any unresolved parts of the atmosphere. These local cutoff temperatures are calculated using only local grid cell quantities, enabling the MHD extension of TRAC to efficiently account for the magnetic field evolution, without tracing field lines. Consistent with analytical predictions, we show that this approach successfully preserves the properties of the original TRAC method. In particular, the total radiative losses and heating remain conserved under the MHD formulation. Results from 2D MHD simulations of impulsive heating in unsheared and sheared arcades of coronal loops are also presented. These simulations benchmark the MHD TRAC method against a series of 1D models and demonstrate the versatility and robustness of the method in multi-dimensional magnetic fields. We show, for the first time, that pressure differences, ...
arXiv comment18 pages, 11 figures, accepted for publication in A&A
References in corpus (26)
- On Solving the Coronal Heating Problem
- Highly Efficient Modeling of Dynamic Coronal Loops
- The Dynamic Formation of Prominence Condensations
- The influence of numerical resolution on coronal density in hydrodynamic models of impulsive heating
- Simulations of Prominence Formation in the Magnetized Solar Corona by Chromospheric Heating
- Enthalpy-based Thermal Evolution of Loops: II. Improvements to the Model
- On the occurrence of thermal non-equilibrium in coronal loops
- An MHD avalanche model of multi-threaded coronal loop
- Impact of Type II Spicules in the Corona: Simulations and Synthetic Observables
- The Role of Asymmetries in Thermal Non-Equilibrium
- 3D MHD modeling of twisted coronal loops
- A fully self-consistent model for solar flares
- The effects of numerical resolution, heating timescales and background heating on thermal non-equilibrium in coronal loops
- Plasma sloshing in pulse-heated solar and stellar coronal loops
- A Fast and Accurate Method to Capture the Solar Corona/Transition Region Enthalpy Exchange
- Phase mixing and wave heating in a complex coronal plasma
- Transition region adaptive conduction (TRAC) in multidimensional magnetohydrodynamic simulations
- A new approach for modelling chromospheric evaporation in response to enhanced coronal heating: I The method
- The Onset of 3D Magnetic Reconnection and Heating in the Solar Corona
- Modelling the solar transition region using an adaptive conduction method
- Self-consistent 3D radiative MHD simulations of coronal rain formation and evolution
- Identifying Observables that can Differentiate Between Impulsive and Footpoint Heating: Time Lags and Intensity Ratios
- The Role of Magnetic Helicity in Coronal Heating
- A new approach for modelling chromospheric evaporation in response to enhanced coronal heating: II Non-uniform heating
- Magnetohydrodynamic waves in braided magnetic fields
- Current systems of coronal loops in 3D MHD simulations