Related papers: Phurbas: An Adaptive, Lagrangian, Meshless, Magnet…
We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field ($\nabla \cdot \mathbf{B} = 0$) on adaptively refined, conformally…
Micromagnetics depends on high-fidelity numerical methods for magnetization dynamics. This work proposes a third-order temporal accuracy scheme for the Landau-Lifshitz-Gilbert equation, addressing accuracy-efficiency trade-offs in existing…
This paper presents an overview and introduction to Smoothed Particle Hydrodynamics and Magnetohydrodynamics in theory and in practice. Firstly, we give a basic grounding in the fundamentals of SPH, showing how the equations of motion and…
Smoothed particle hydrodynamics (SPH) has been extensively studied in computer graphics to animate fluids with versatile effects. However, SPH still suffers from two numerical difficulties: the particle deficiency problem, which will…
We present the implementation of a three-dimensional, second order accurate Godunov-type algorithm for magneto-hydrodynamic (MHD), in the adaptive-mesh-refinement (AMR) cosmological code {\tt CHARM}. The algorithm is based on the full…
We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves…
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…
We present several implementations of the Metropolis method, an adaptive Monte Carlo algorithm, which allow for the calculation of multi-dimensional integrals over arbitrary on-shell four-momentum phase space. The Metropolis technique…
Filtering in spatially-extended dynamical systems is a challenging problem with significant practical applications such as numerical weather prediction. Particle filters allow asymptotically consistent inference but require infeasibly large…
A new fully discrete linearized $H^1$-conforming Lagrange finite element method is proposed for solving the two-dimensional magneto-hydrodynamics equations based on a magnetic potential formulation. The proposed method yields numerical…
The Athena MHD code has been extended to integrates the motion of particles coupled with the gas via aerodynamic drag, in order to study the dynamics of gas and solids in protoplanetary disks and the formation of planetesimals. Our…
Invariant finite-difference schemes for the one-dimensional shallow water equations in the presence of a magnetic field for various bottom topographies are constructed. Based on the results of the group classification recently carried out…
Lagrangian particle tracking is essential for characterizing turbulent flows, but inferring particle acceleration from inherently noisy position data remains a significant challenge. Fluid particles in turbulence experience extreme,…
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…
A numerical method based on smoothed particle hydrodynamics with adaptive spatial resolution (SPH-ASR) was developed for simulating free surface flows. This method can reduce the computational demands while maintaining the numerical…
Methods to solve the relativistic hydrodynamic equations are a key computational kernel in a large number of astrophysics simulations and are crucial to understanding the electromagnetic signals that originate from the merger of…
Hydrodynamical simulations are the most accurate way to model structure formation in the universe, but they often involve a large number of astrophysical parameters modeling subgrid physics, in addition to cosmological parameters. This…
Relativistic strongly magnetized plasmas are produced in laboratories thanks to state-of-the-art laser technology but can naturally be found around compact objects such as neutron stars and black holes. Detailed studies of the behaviour of…
The local approach to computing electrostatic interactions proposed by Maggs and adapted by Pasichnyk for molecular dynamics simulations is extended to situations where the dielectric background medium is inhomogeneous. We furthermore…
Mesh-free Lagrangian methods are widely used for simulating fluids, solids, and their complex interactions due to their ability to handle large deformations and topological changes. These physics simulators, however, require substantial…