Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D

Abstract

Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.

arXiv comment11 pages, 8 figures

References in corpus (21)

• The stellar atmosphere simulation code Bifrost
• MPI-AMRVAC 2.0 for Solar and Astrophysical Applications
• Extension of the MURaM radiative MHD code for coronal simulations
• Preserving Monotonicity in Anisotropic Diffusion
• Intermittent heating in the solar corona employing a 3D MHD model
• Nonlinear Evolution of the Magnetothermal Instability in Two Dimensions
• Magneto-acoustic waves in sunspots: first results from a new 3D nonlinear magnetohydrodynamic code
• Three-dimensional simulations of solar magneto-convection including effects of partial ionization
• Diagnosing the time-dependence of active region core heating from the emission measure: I. Low-frequency nanoflares
• Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: I. The Numerical Model
• Observationally driven 3D MHD model of the solar corona above an active region
• A model for the formation of the active region corona driven by magnetic flux emergence
• Nanoflare statistics in an active region 3D MHD coronal model
• MHS sunspot model from deep sub-photospheric to chromospheric layers
• Numerical Simulation of Solar Microflares in a Canopy-Type Magnetic Configuration
• MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
• A Fast Semi-implicit Method for Anisotropic Diffusion
• Thermal conductivity of micro/nano-porous polymers: Prediction models and applications
• Generalized Fluid Models of the Braginskii Type
• MAGNUS: A new resistive MHD code with heat flow terms
• Thermal conduction effects on formation of chromospheric solar tadpole-like jets

Cited by in corpus (5)

• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Two-fluid numerical model of chromospheric heating and plasma outflows in a quiet-Sun
• Time-Adaptive PIROCK Method with Error Control for Multi-Fluid and Single-Fluid MHD Systems
• Implementation of thermal conduction energy transfer models in the Bifrost Solar atmosphere MHD code
• Future Prospects for Partially Ionised Solar Plasmas: the Prominence Case