Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry
arXiv:2203.06272 · doi:10.1051/0004-6361/202141449
Abstract
Ambipolar diffusion is a process occurring in partially ionised astrophysical systems that imparts a complicated mathematical and physical nature to Ohm's law. The numerical codes that solve the magnetohydrodynamic (MHD) equations have to be able to deal with the singularities that are naturally created in the system by the ambipolar diffusion term. The global aim is to calculate a set of theoretical self-similar solutions to the nonlinear diffusion equation with cylindrical symmetry that can be used as tests for MHD codes which include the ambipolar diffusion term. First, following the general methods developed in the applied mathematics literature, we obtained the theoretical solutions as eigenfunctions of a nonlinear ordinary differential equation. Phase-plane techniques were used to integrate through the singularities at the locations of the nulls, which correspond to infinitely sharp current sheets. In the second half of the paper, we consider the use of these solutions as tests for MHD codes. To that end, we used the Bifrost code, thereby testing the capabilities of these solutions as tests as well as (inversely) the accuracy of Bifrost's recently developed ambipolar diffusion module. The obtained solutions are shown to constitute a demanding, but nonetheless viable, test for MHD codes that incorporate ambipolar diffusion. The Bifrost code is able to reproduce the theoretical solutions with sufficient accuracy up to very advanced diffusive times. Using the code, we also explored the asymptotic properties of our theoretical solutions in time when initially perturbed with either small or finite perturbations. The functions obtained in this paper are relevant as physical solutions and also as tests for general MHD codes. They provide a more stringent and general test than the simple Zeldovich-Kompaneets-Barenblatt-Pattle solution.
arXiv comment15 pages, 10 figures and accompanying animations. Accepted for publication in Astronomy & Astrophysics, Section Astrophysical Processes
References in corpus (37)
- Theory of Star Formation
- The stellar atmosphere simulation code Bifrost
- Global simulations of protoplanetary disks with ohmic resistivity and ambipolar diffusion
- Numerical simulation of the three-dimensional structure and dynamics of the non-magnetic solar chromosphere
- A publicly available simulation of an enhanced network region of the Sun
- 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
- Partially Ionized Plasmas in Astrophysics
- Non-equilibrium hydrogen ionization in 2D simulations of the solar atmosphere
- 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
- Emergence of a flux tube through a partially ionised solar atmosphere
- Three-dimensional simulations of solar magneto-convection including effects of partial ionization
- Ionized Plasma and Neutral Gas Coupling in the Sun's Chromosphere and Earth's Ionosphere/Thermosphere
- Heating mechanisms in the low solar atmosphere through magnetic reconnection in current sheets
- Formation and Collapse of Nonaxisymmetric Protostellar Cores in Planar Magnetic Molecular Clouds
- Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere
- Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES
- On the minimum temperature of the quiet solar chromosphere
- A three-dimensional numerical method for modelling weakly ionized plasmas
- Magnetohydrodynamics of the Weakly Ionized Solar Photosphere
- Ion-neutral interactions and non-equilibrium ionization in the solar chromosphere
- Fast Reconnection in a Two-Stage Process
- Three-dimensional Simulation of Magnetized Cloud Fragmentation Induced by Nonlinear Flows and Ambipolar Diffusion
- Ambipolar Drift Heating in Turbulent Molecular Clouds
- An Explicit Scheme for Incorporating Ambipolar Diffusion in a Magnetohydrodynamics Code
- Numerical simulations of quiet Sun magnetic fields seeded by Biermann battery
- MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
- Effect of ion-neutral collisions in simulations of emerging active regions
- Nonequilibrium ionization and ambipolar diffusion in solar magnetic flux emergence processes
- A Magnetic Reconnection model for Hot Explosions in the Cool Atmosphere of the Sun
- Ambipolar diffusion in the Bifrost code
- A Simflowny-based high-performance 3D code for the generalized induction equation
- Joint action of Hall and ambipolar effects in 3D magneto-convection simulations of the quiet Sun. I. Dissipation and generation of waves
- The challenges of modelling microphysics: ambipolar diffusion, chemistry, and cosmic rays in MHD shocks
- Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry
Cited by in corpus (3)
- Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry
- Implementation of thermal conduction energy transfer models in the Bifrost Solar atmosphere MHD code
- On Thermal Conduction in the Solar Atmosphere: An Analytical Solution for Nonlinear Diffusivity without Compact Support