Related papers: The 3D MHD code GOEMHD3 for large-Reynolds-number …
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…
This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl. \textbf{182}, 310…
Large Eddy Simulations (LES) of gyrokinetic plasma turbulence are investigated as interesting candidates to decrease the computational cost. A dynamic procedure is implemented in the GENE code, allowing for dynamic optimization of the free…
In many plasmas, in particular in space science, protons govern much of the essential physics. Minority ions, suprathermal tails, and electrons at times account for additional details. However, electron effects usually appear on much…
We present results of large-scale three-dimensional simulations of weakly magnetized supersonic turbulence at grid resolutions up to 1024^3 cells. Our numerical experiments are carried out with the Piecewise Parabolic Method on a Local…
Numerical simulations of magneto-convection have greatly expanded our understanding of stellar interiors and stellar magnetism. Recently, fully compressible hydrodynamical simulations of full-star models have demonstrated the feasibility of…
(Abridged) In the implicit large eddy simulation (ILES) paradigm, the dissipative nature of high-resolution shock-capturing schemes is exploited to provide an implicit model of turbulence. Recent 3D simulations suggest that turbulence might…
The core collapse events observed in a stellarator experiment are studied by a three-dimensional nonlinear MHD simulations. In the low magnetic shear configuration like the Wendelstein 7-X, the rotational transform profile is very sensitive…
We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms…
Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects…
Fully realizing the potential of multigrid solvers often requires custom algorithms for a given application model, discretizations and even regimes of interest, despite considerable effort from the applied math community to develop fully…
Three-dimensional (3D) radiative magnetohydrodynamics (MHD) simulations are the only way to model stellar atmospheres without any ad hoc parameterisations. Several 3D radiative MHD codes have achieved good quantitative agreement with…
Traditional 1D system thermal hydraulic analysis has been widely applied in nuclear industry for licensing purposes due to its numerical efficiency. However, such codes are inherently deficient for modeling of multiscale multidimensional…
AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous…
We have performed high resolution 3D simulations with adaptive mesh refinement, following the ISM evolution in a star forming galaxy both on small (<1 pc) and large (>10 kpc) scales, enabling us to track structures in cooling shock…
Giant Molecular Clouds (GMCs) are dominated by supersonic turbulence, creating a complex network of shocks and filaments that regulate star formation. While the global inefficiency of star formation is well-observed, predicting exactly…
We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code…
Simulating plasmas in the Hall-MagnetoHydroDynamics (Hall-MHD) regime represents a valuable {approach for the investigation of} complex non-linear dynamics developing in astrophysical {frameworks} and {fusion machines}. Taking into account…
This article presents the development and validation of a hybrid multi-sphere discrete element framework - Rigid3D, for the simulation of granular systems with arbitrarily shaped particles in 3D space. In this DEM framework, a non-spherical…
Numerical methods to improve the treatment of magnetic fields in smoothed field magnetohydrodynamics (SPMHD) are developed and tested. Chapter 2 is a review of SPMHD. In Chapter 3, a mixed hyperbolic/parabolic scheme is developed which…