Ambipolar diffusion in the Bifrost code
arXiv:2004.11927 · doi:10.1051/0004-6361/202037809
Abstract
Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key in many different astrophysical systems. However, understanding its effects is challenging due to basic uncertainties concerning relevant microphysical aspects and the strong constraints it imposes on the numerical modeling. Our aim is to introduce a numerical tool that allows us to address complex problems involving ambipolar diffusion in which, additionally, departures from ionization equilibrium are important or high resolution is needed. The primary application of this tool is for solar atmosphere calculations, but the methods and results presented here may also have a potential impact on other astrophysical systems. We have developed a new module for the stellar atmosphere Bifrost code that improves its computational capabilities of the ambipolar diffusion term in the Generalized Ohm's Law. This module includes, among other things, collision terms adequate to processes in the coolest regions in the solar chromosphere. As a key feature of the module, we have implemented the Super Time-Stepping (STS) technique, that allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code. We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium but also the main electron donors like sodium, silicon and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution.
Accepted in A&A, 10 pages, 7 figures
References in corpus (48)
- The stellar atmosphere simulation code Bifrost
- Global simulations of protoplanetary disks with ohmic resistivity and ambipolar diffusion
- The Balbus-Hawley instability in weakly ionised discs
- Dynamic Hydrogen Ionization
- Heating of the magnetized solar chromosphere by partial ionization effects
- Magnetically Regulated Star Formation in 3D: The Case of Taurus Molecular Cloud Complex
- On the generation of solar spicules and Alfvénic waves
- Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks
- Nonideal MHD Effects and Magnetic Braking Catastrophe in Protostellar Disk Formation
- Partially Ionized Plasmas in Astrophysics
- Non-equilibrium hydrogen ionization in 2D simulations of the solar atmosphere
- Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods
- Magnetic Braking and Protostellar Disk Formation: Ambipolar Diffusion
- Protostellar Disk Formation Enabled by Removal of Small Dust Grains
- Fast magnetic reconnection in the solar chromosphere mediated by the plasmoid instability
- 2D Radiative MHD Simulations of the Importance of Partial Ionization in the Chromosphere
- Collisional dissipation of Alfvén waves in a partially ionised solar chromosphere
- Multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma
- Collisions, magnetization, and transport coefficients in the lower solar atmosphere
- Ambipolar Drift in a Turbulent Medium
- Onset of Fast Reconnection in Hall Magnetohydrodynamics Mediated by the Plasmoid Instability
- Magnetohydrodynamic Waves in a Partially Ionized Filament Thread
- Energy flux of Alfven waves in weakly ionized plasma
- Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere
- Formation and Collapse of Nonaxisymmetric Protostellar Cores in Planar Magnetic Molecular Clouds
- Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects
- Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES
- Heating of the partially ionized solar chromosphere by waves in magnetic structures
- Si IV Resonance Line Emission During Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations
- On the minimum temperature of the quiet solar chromosphere
- A three-dimensional numerical method for modelling weakly ionized plasmas
- Damping of Alfven waves in solar partially ionized plasmas: effect of neutral helium in multi-fluid approach
- Magnetohydrodynamics of the Weakly Ionized Solar Photosphere
- Ion-neutral interactions and non-equilibrium ionization in the solar chromosphere
- Non-equilibrium helium ionization in an MHD simulation of the solar atmosphere
- Three-dimensional Simulation of Magnetized Cloud Fragmentation Induced by Nonlinear Flows and Ambipolar Diffusion
- Ambipolar Drift Heating in Turbulent Molecular Clouds
- Effect of coronal temperature on the scale of solar chromospheric jets
- An Explicit Scheme for Incorporating Ambipolar Diffusion in a Magnetohydrodynamics Code
- MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
- The Role of Partial Ionization Effects in the Chromosphere
- Nonequilibrium ionization and ambipolar diffusion in solar magnetic flux emergence processes
- On the Misalignment between Chromospheric Features and the Magnetic Field on the Sun
- On Molecular Hydrogen Formation and the Magnetohydrostatic Equilibrium of Sunspots
- The Formation Height of Millimeter-wavelength Emission in the Solar Chromosphere
- On the importance of the nonequilibrium ionization of Si IV and O IV and the line-of-sight in solar surges
- The challenges of modelling microphysics: ambipolar diffusion, chemistry, and cosmic rays in MHD shocks
- Non-Equilibrium Spectrum Formation Affecting Solar Irradiance
Cited by in corpus (18)
- MPI-AMRVAC 3.0: updates to an open-source simulation framework
- Probing the physics of the solar atmosphere with the Multi-slit Solar Explorer (MUSE): I. Coronal Heating
- Transition region adaptive conduction (TRAC) in multidimensional magnetohydrodynamic simulations
- Numerical simulations and observations of Mg II in the solar chromosphere
- ALMA and IRIS Observations of the Solar Chromosphere I: an On-Disk Type II Spicule
- Evidence of multithermal nature of spicular downflows. Impact on solar atmospheric heating
- Efficient numerical treatment of ambipolar and Hall drift as hyperbolic system
- Shocks and instabilities in the partially ionised solar atmosphere
- Multi-fluid Simulation of Solar Chromospheric Turbulence and Heating Due to the Thermal Farley-Buneman Instability
- Chromosphere of the quiet sun -- I. Shock and current-sheet dynamics and heating
- Time-Adaptive PIROCK Method with Error Control for Multi-Fluid and Single-Fluid MHD Systems
- Designing wavelength sampling for Fabry-Pérot observations. Information-based spectral sampling
- Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry
- Fast magneto-acoustic waves in the solar chromosphere: Comparison of single-fluid and two-fluid approximations
- Systematic 2.5 D resistive MHD simulations with ambipolar diffusion and Hall effect for fast magnetic reconnection
- A new view of the solar interface region from the Interface Region Imaging Spectrograph (IRIS)
- Radiation-magnetohydrodynamics with MPI-AMRVAC using flux-limited diffusion
- Effects of ambipolar diffusion on waves in the solar chromosphere