Related papers: Application of an Upwind Integration Method to Pla…
Magnetohydrodynamic (MHD) simulations are indispensable research infrastructure in astrophysics today. In order to satisfy the solenoidal constraint of the MHD equations on discretized grids, modern simulation codes often employ either…
In this work, we introduce three algorithmic improvements to reduce the cost and improve the scaling of orbital space variational Monte Carlo (VMC). First, we show that by appropriately screening the one- and two-electron integrals of the…
Recent developments in high peak-power table-top laser systems reaching highly relativistic light intensities have led to significant advances in laser-driven particle acceleration schemes (mainly the laser wakefield acceleration, LWFA)…
We investigate particle acceleration in an MHD-scale system of multiple current sheets by performing 2D and 3D MHD simulations combined with a test particle simulation. The system is unstable for the tearing-mode instability, and magnetic…
Through a direct comparison between numerical simulations in two and three dimensions, we investigate topological effects in reconnection. A simple estimate on increase in reconnection rate in three dimensions by a factor of $\sqrt{2}$,…
Many astrophysical simulations involve extreme dynamic range of timescales around 'special points' in the domain (e.g. black holes, stars, planets, disks, galaxies, shocks, mixing interfaces), where processes on small scales couple strongly…
We present a new numerical scheme which combines the Spectral Difference (SD) method up to arbitrary high order with \emph{a-posteriori} limiting using the classical MUSCL-Hancock scheme as fallback scheme. It delivers very accurate…
Equilibrium configurations of the internal magnetic field of a pulsar play a key role in modeling astrophysical phenomena, from glitches to gravitational wave emission. In this paper we present a numerical scheme for solving the…
We design a conservative finite difference scheme for ideal magnetohydrodynamic simulations that attains high-order accuracy, shock-capturing, and divergence-free condition of the magnetic field. The scheme interpolates pointwise physical…
We present a new Eulerian framework for the computation of turbulent compressible multiphase channel flows, specifically to assess turbulence modulation by dispersed particulate matter in dilute concentrations but with significant mass…
The numerical simulation of the 3D incompressible Euler equation is analyzed with respect to different integration methods. The numerical schemes we considered include spectral methods with different strategies for dealiasing and two…
We consider the effect of Hall-drift on the evolution of magnetic fields in neutorn star crusts using a finite difference code. Previously, most studies have focused on the situation where the ohmic term dominates the evolution. Here, we…
We present a new high-order accurate discretisation on unstructured meshes of quadrilateral elements. Our Face Upwinded Spectral Element (FUSE) method uses the same node distribution as a high-order continuous Galerkin (CG) method, but with…
The Reynolds-Averaged Navier-Stokes equations and the Large-Eddy Simulation equations can be coupled using a transition function to switch from a set of equations applied in some areas of a domain to the other set in the other part of the…
The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model [Park et. al., 1992] has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling…
A transient magneto-quasistatic vector potential formulation involving nonlinear material is spatially discretized using the finite element method of first and second polynomial order. By applying a generalized Schur complement the…
The Eulerian space-time correlation of strong Magnetohydrodynamic (MHD) turbulence in strongly magnetized plasmas is investigated by means of direct numerical simulations of Reduced MHD turbulence and phenomenological modeling. Two new…
We have extended the ANTARES code to simulate the coupling of pulsation with convection in Cepheid-like variables in an increasingly realistic way, in particular in multidimensions, 2D at this stage. Present days models of radially…
We present calculations of magnetic field evolution by the Hall effect and Ohmic decay in neutron star crusts. We first discuss when and for which neutron stars the Hall effect is important relative to Ohmic decay. We compute the evolution…
For the first time, the actuator line model (ALM) and the immersed boundary method (IBM) are integrated into the high-order gas-kinetic scheme (GKS) to simulate the wind turbine with the nacelle and tower. The high-order GKS is extended to…