Related papers: Piecewise Parabolic Method on a Local Stencil for …

In this paper, a perfectly matched layer (PML) method is proposed to solve the time-domain electromagnetic scattering problems in 3D effectively. The PML problem is defined in a spherical layer and derived by using the Laplace transform and…

The numerical simulation of multiphase flows involving dispersed components with large scale disparities, such as the collisions between millimeter-sized bubbles and micron-sized mineral particles in flotation, poses a significant…

We develop a positivity-preserving (PP) PAMPA (Point-Average-Moment PolynomiAl-interpreted) scheme that enforces a discrete divergence-free (DDF) magnetic field for ideal MHD on Cartesian grids. Extending our 1D invariant-domain-preserving…

Magnetohydrodynamic (MHD) turbulence is a crucial component of the current paradigms of star formation, dynamo theory, particle transport, magnetic reconnection and evolution of structure in the interstellar medium (ISM) of galaxies.…

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent…

Multiscale problems are computationally costly to solve by direct simulation because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and analyzed new numerical methods…

In this work we introduce a novel semi-implicit structure-preserving finite-volume/finite-difference scheme for the viscous and resistive equations of magnetohydrodynamics (MHD) based on an appropriate 3-split of the governing PDE system,…

Turbulence simulations play a key role in advancing the general understanding of the physical properties turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence…

We present simulations of coherent structures in compressible flows near the transition to turbulence using the Dissipative Particle Dynamics (DPD) method. The structures we find are remarkably consistent with experimental observations and…

The inverse transfer of magnetic helicity is investigated through direct numerical simulations of large-scale-mechanically-driven turbulent flows in the isothermal ideal magnetohydrodynamics (MHD) framework. The mechanical forcing is either…

In this work we describe a $\delta f$ particle simulation method where the bulk density is periodically remapped on a coarse spline grid using a Forward-Backward Lagrangian (FBL) approach. This method is designed to handle plasma regimes…

In fluid dynamical simulations in astrophysics, large deformations are common and surface tracking is sometimes necessary. Smoothed Particle Hydrodynamics (SPH) method has been used in many of such simulations. Recently, however, it has…

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 Pseudo-Marginal (PM) algorithm is a popular Markov chain Monte Carlo (MCMC) method used to sample from a target distribution when its density is inaccessible, but can be estimated with a non-negative unbiased estimator. Its performance…

Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence…

Lattice Boltzmann Methods (LBM) stand out for their simplicity and computational efficiency while offering the possibility of simulating complex phenomena. While they are optimal for Cartesian meshes, adapted meshes have traditionally been…

We develop a provably energy stable discontinuous Galerkin spectral element method (DGSEM) approximation of the perfectly matched layer (PML) for the three and two space dimensional (3D and 2D) linear acoustic wave equations, in first order…

We develop a structure-preserving computational framework for acoustic wave scattering by moving objects, comprising a new PML-domain-embedding model and a compatible numerical approximation. The model couples a perfectly matched layer…

A modified method is presented to generate artificial magnetic turbulence that is used for test-particle simulations. Such turbulent fields are obtained from the superposition of a set of wave modes with random polarizations and random…

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy…