Abstract:A two-dimensional Eulerian finite difference code is developed for solving the non-ideal magnetohydrodynamic (MHD) equations including the effects of self-consistent magnetic field, thermal conduction, resistivity, gravity, and radiation transfer, which when combined with specified pulse current models and plasma equations of state, can be used as a unified lightning return stroke solver. The differential equations are written in the covariant form in the cylindrical geometry and kept in the conservative form which enables some high-accuracy shock capturing schemes to be equipped in the lightning channel configuration naturally. In this code, the 5-order weighted essentially non-oscillatory scheme combined with Lax-Friedrichs flux splitting method is introduced for computing the convection terms of the MHD equations. The 3-order total variation diminishing Runge-Kutta integral operator is also equipped to keep the time-space accuracy of consistency. The numerical algorithms for non-ideal terms, e.g., artificial viscosity, resistivity, and thermal conduction, are introduced in the code via operator splitting method. This code assumes the radiation is in local thermodynamic equilibrium with plasma components and the flux limited diffusion algorithm with grey opacities is implemented for computing the radiation transfer. The transport coefficients and equation of state in this code are obtained from detailed particle population distribution calculation, which makes the numerical model is self-consistent. This code is systematically validated via the Sedov blast solutions and then used for lightning return stroke simulations with the peak current being 20 kA, 30 kA, and 40 kA, respectively. The results show that this numerical model consistent with observations and previous numerical results. The population distribution evolution and energy conservation problems are also discussed.

Numerical code SELFAS-3 and electrodynamic aggregation of magnetized nanodust

Dynamic percolation of electric conductivity and a trend towards fractal skeletal structuring in a random ensemble of magnetized nanodust

Multiscale Modeling of 3-D Electromagnetic Fields With Magnetization Dynamics

Numerical study on the aggregation of magnetic particles with power law interactions

A self-consistent spin-diffusion model for micromagnetics

Simulating the Self-Assembly and Hysteresis Loops of Ferromagnetic Nanoparticles with Sticking of Ligands

Computational General Relativistic Force-Free Electrodynamics: II. Characterization of Numerical Diffusivity

Nonlinear magnetohydrodynamics simulation using high-order finite elements

Monte Carlo Simulations of Interacting Magnetic Nanoparticles

Numerical investigation of magnetic-field induced self-assembly of nonmagnetic particles in magnetic fluids

Large-scale micromagnetics simulations with dipolar interaction using all-to-all communications

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

A Unified Radiative Magnetohydrodynamics Code for Lightning-Like Discharge Simulations

LCODE: Quasistatic code for simulating long-term evolution of three-dimensional plasma wakefields

Micromagnetic simulations of thermally activated magnetization reversal of nanoscale magnets

Diffusive transport of a 2-D magnetized dusty plasma cloud

Quasi-2D Monte Carlo simulations of the regime change in the aggregates of magnetic cubic particles on a material surface

A Machine-Learning Augmented Mixed Method for Simulating Α Particle Transport in Inertial Confinement Fusion

Force Fields for Molecular Dynamics Simulations of Charged Dust Particles with Finite Size in Complex Plasmas

Particle-in-cell Simulation Method for Macroscopic Degenerate Plasmas.

Numerical simulation of cold emission in coaxial diode with magnetic isolation