Related papers: $\nabla\cdot{B}=0$ versus the Universe
Understanding the origin and evolution of magnetic fields on cosmological scales opens up a window into the physics of the early Universe. Numerical simulations of such fields require a careful treatment to faithfully solve the equations of…
Smoothed particle hydrodynamics (SPH) offers distinct advantages for modeling many engineering problems, yet achieving high-order consistency in its conservative formulation remains to be addressed. While zero- and higher-order…
We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled…
Precise estimation of cosmological parameters from the cosmic microwave background (CMB) remains a central goal of modern cosmology and a key test of inflationary physics. However, this task is fundamentally limited by strong foreground…
We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because…
A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations…
We consider the problem of high-dimensional filtering of state-space models (SSMs) at discrete times. This problem is particularly challenging as analytical solutions are typically not available and many numerical approximation methods can…
We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is…
We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing…
We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms…
We introduce the concept of compressed convolution, a technique to convolve a given data set with a large number of non-orthogonal kernels. In typical applications our technique drastically reduces the effective number of computations. The…
In this paper, we present an advanced analysis of near optimal algorithms that use limited space to solve the frequency estimation, heavy hitters, frequent items, and top-k approximation in the bounded deletion model. We define the family…
We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. Local, third-order, least-squares, polynomial…
Finding fast yet accurate numerical solutions to the Helmholtz equation remains a challenging task. The pollution error (i.e. the discrepancy between the numerical and analytical wave number k) requires the mesh resolution to be kept fine…
A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained…
Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure…
We establish the geometric ergodicity of the preconditioned Hamiltonian Monte Carlo (HMC) algorithm defined on an infinite-dimensional Hilbert space, as developed in [Beskos et al., Stochastic Process. Appl., 2011]. This algorithm can be…
Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct…
In astrophysical magnetohydrodynamics (MHD) and electrodynamics simulations, numerically enforcing the divB=0 constraint on the magnetic field has been difficult. We observe that for point-based discretization, as used in finite-difference…
We perform a validation of estimates of the growth rate of structure, described by the parameter combination $f\sigma_8$, in modified gravity cosmologies. We consider an analysis pipeline based on the redshift-space distortion modelling of…