Related papers: WHAM: A WENO-based general relativistic numerical …
Many numerical schemes for hyperbolic systems require a piecewise polynomial reconstruction of the cell averaged values, and to simulate perturbed steady states accurately we require a so called 'well balanced' reconstruction scheme. For…
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…
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…
In this paper, we develop high-order, conservative, non-splitting Eulerian-Lagrangian (EL) Runge-Kutta (RK) finite volume (FV) weighted essentially non-oscillatory (WENO) schemes for convection-diffusion equations. The proposed…
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…
A non-hydrostatic depth-averaged model for dry granular flows is proposed, taking into account vertical acceleration. A variable friction coefficient based on the $\mu(I)$ rheology is considered. The model is obtained from an asymptotic…
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to…
In subsurface multiphase flow simulations, poor nonlinear solver performance is a significant runtime sink. The system of fully implicit mass balance equations is highly nonlinear and often difficult to solve for the nonlinear solver,…
A self-consistent nonlinear hydrodynamic theory is presented of the propagation of a long and thin relativistic electron beam, for a typical plasma wake field acceleration configuration in an unmagnetized and overdense plasma. The random…
For the simulation of compressible flow with a broadband of length scales and discontinuities, the WENO schemes using incremental stencil sizes other than uniform ones are promising for more robustness and less numerical dissipation.…
General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole shadows from the…
In this paper, a third order gas kinetic scheme is developed on the three dimensional hybrid unstructured meshes for the numerical simulation of compressible inviscid and viscous flows. A third-order WENO reconstruction is developed on the…
Achieving accurate numerical results of hydrodynamic loads based on the potential-flow theory is very challenging for structures with sharp edges, due to the singular behavior of the local-flow velocities. In this paper, we introduce the…
Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was…
Steady plasma flows have been studied almost exclusively in systems with continuous symmetry or in open domains. In the absence of continuous symmetry, the lack of a conserved quantity makes the study of flows intrinsically challenging. In…
Non-linear wave-driven processes in plasmas are normally described by either a monochromatic pump wave that couples to other monochromatic waves, or as a random phase wave coupling to other random phase waves. An alternative approach…
In the study of plasma, particularly in applications involving strong laser-plasma interactions, the propagation of a strong electromagnetic wave induces relativistic velocities in the electron flow. Given such conditions, the wave…
At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be…
We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity…
Fluid discontinuities, such as shock fronts and vortex sheets, can reflect waves and become unstable to corrugation. Analytical calculations of these phenomena are tractable in the simplest cases only, while their numerical simulations are…