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.

The Athena Astrophysical MHD Code in Cylindrical Geometry

The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests

Analysis of Liquid Metal MHD Fluid Flow and Heat Transfer Using the KAT Code

The magnetohydrodynamic-particle-in-cell module in athena plus plus : implementation and code tests

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

A Proper Discretization of Hydrodynamic Equations in Cylindrical Coordinates for Astrophysical Simulations

A new MHD code with adaptive mesh refinement and parallelization for astrophysics

Conductivity of Metallic Si:B near the Metal-Insulator Transition: Comparison between Unstressed and Uniaxially Stressed Samples

A Multifluid Dust Module in Athena++: Algorithms and Numerical Tests

HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications

GR-Athena++: magnetohydrodynamical evolution with dynamical space-time

PARTICLE–GAS DYNAMICS WITH ATHENA: METHOD AND CONVERGENCE

The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity

A globally conservative finite element MHD code and its application to the study of compact torus formation, levitation and magnetic compression

Development of a gyrokinetic-MHD energetic particle simulation code Part I: MHD version

DISCO: A 3D MOVING-MESH MAGNETOHYDRODYNAMICS CODE DESIGNED FOR THE STUDY OF ASTROPHYSICAL DISKS

A Computational Code for Two-Dimensional Unsteady Magnetohydrodynamics by the Method of Characteristics

Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry

Pegasus: A New Hybrid-Kinetic Particle-in-cell Code for Astrophysical Plasma Dynamics.

Local Spherical Collapsing Box in Athena++: Numerical Implementation and Benchmark Tests