Related papers: Application of an Upwind Integration Method to Pla…
The constrained transport (CT) method reflects the state of the art numerical technique for preserving the divergence-free condition of magnetic field to machine accuracy in multi-dimensional MHD simulations performed with Godunov-type, or…
This paper introduces multidimensional algorithms for simulating multiphase flows, leveraging the wave structure of the Euler equations in characteristic space and the physical properties of variables in physical space. The algorithm…
In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed particle magnetohydrodynamics (Spmhd) simulations and study the…
We present 2D MHD numerical simulations of tearing-unstable current sheets coupled to a population of non-thermal test-particles, in order to address the problem of numerical convergence with respect to grid resolution, numerical method and…
High-order accurate summation-by-parts (SBP) finite difference (FD) methods constitute efficient numerical methods for simulating large-scale hyperbolic wave propagation problems. Traditional SBP FD operators that approximate first-order…
The nonlinear magnetic induction equation with Hall effect can be used to model magnetic fields, e.g. in astrophysical plasma environments. In order to give reliable results, numerical simulations should be carried out using effective and…
A general method for constructing high order upwind schemes for multidimensional magnetohydrodynamics (MHD), having as a main built-in condition the divergence-free constraint $\divb=0$ for the magnetic field vector $\bb$, is proposed. The…
A numerical analysis of the effect of artificial viscosity is undertaken in order to understand the effect of numerical diffusion on numerical boundary feedback control. The analysis is undertaken on the linear hyperbolic systems…
We investigate magnetic reconnection in systems simultaneously containing asymmetric (anti-parallel) magnetic fields, asymmetric plasma densities and temperatures, and arbitrary in-plane bulk flow of plasma in the upstream regions. Such…
We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation…
Magnetised turbulence is ubiquitous in astrophysical systems, where it notoriously spans a broad range of spatial scales. Phenomenological theories of MHD turbulence describe the self-similar dynamics of turbulent fluctuations in the…
Numerical schemes for the solution of the Euler equations have recently been developed, which involve the discretisation of the internal energy equation, with corrective terms to ensure the correct capture of shocks, and, more generally,…
The particle-in-cell approach has proven effective at modeling neutron star and black hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where…
Recently, several authors have suggested the use of first order methods, such as fast dual ascent and the alternating direction method of multipliers, for embedded model predictive control. The main reason is that they can be implemented…
We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements…
The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere…
The numerical simulation of three-dimensional charged-particle dynamics (CPD) under strong magnetic field is a basic and challenging algorithmic task in plasma physics. In this paper, we introduce a new methodology to design two-scale…
In this paper, we present a fully local second-order upwind scheme, applicable on generic meshes. This is done by hybridisation, which is achieved by introducing unknowns on each edge of the mesh. By doing so, fluxes only depend on values…
Large eddy simulations with three-dimensional finite-volume Conditional Moment Closure (CMC) model are performed for a hydrogen / nitrogen lifted flame with detailed chemical meachanism. The emphasis is laid on the influences of mesh…
Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic ray (CR) modified plane shock evolution using the two-fluid model for diffusive particle acceleration. The numerical shocks evolve to…