Related papers: Pseudo-spectral solver versus grid-based solver: A…
Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales.…
This paper presents a new code for performing multidimensional radiation hydrodynamic (RHD) simulations on parallel computers involving anisotropic radiation fields and nonequilibrium effects. The radiation evolution modules described here…
Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these…
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy…
We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
High-order Discontinuous Galerkin (DG) methods offer excellent accuracy for turbulent flow simulations, especially when implemented on GPU-oriented architectures that favor very high polynomial orders. On modern GPUs, high-order polynomial…
The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an…
The dynamics of energetic particle (EP) species, born from fusion reactions or plasma heating schemes, are critical for predicting the behavior of magnetic confinement fusion experiments and future fusion reactors. Because energetic…
The calculation of a three-dimensional underwater acoustic field has always been a key problem in computational ocean acoustics. Traditionally, this solution is usually obtained by directly solving the acoustic Helmholtz equation using a…
A method for modelling non-Newtonian fluids (dilatants and pseudoplastics) by a power law under the Godunov-Peshkov-Romenski model is presented, along with a new numerical scheme for solving this system. The scheme is also modified to solve…
Since numerical computing with MATLAB offers a wide variety of advantages, such as easier developing and debugging of computational codes rather than lower-level languages, the popularity of this tool is significantly increased in the past…
Magnetohydrodynamic turbulence is a ubiquitous phenomenon in solar physics, plasma physics and astrophysics and governs many properties of the flows of well-conductive fluids. Recently, conflicting spectral slopes for the inertial range of…
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,…
This paper extends the author's previous two-dimensional work with Ou and LeVeque to high-resolution finite volume modeling of systems of fluids and poroelastic media in three dimensions, using logically rectangular mapped grids. A method…
We present three-dimensional solutions of the magnetohydrostatic equations in the co-rotating frame of reference outside a magnetized rigidly rotating cylinder. We make no symmetry assumption for the magnetic field, but to be able to make…
This paper proposes a novel GPU-based active nematic flow solver based on the smoothed particle hydrodynamics (SPH) method. Nematohydrodynamics equations are discretized using the SPH algorithm, and the periodic domain is enforced using the…
A new numerical code, called SFUMATO, for solving self-gravitational magnetohydrodynamics (MHD) problems using adaptive mesh refinement (AMR) is presented. A block-structured grid is adopted as the grid of the AMR hierarchy. The total…
Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon…