Related papers: Novel Kinetic 3D MHD Algorithm for High Performanc…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
We discuss the problem of ultra high energy particles propagation in astrophysical backgrounds. We present two different computational schemes based on both kinetic and Monte Carlo approaches. The kinetic approach is an analytical…
Numerical studies of shock waves in large scale systems via kinetic simulations with millions of particles are too computationally demanding to be processed in serial. In this work we focus on optimizing the parallel performance of a…
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…
Many applications of computational fluid dynamics require multiple simulations of a flow under different input conditions. In this paper, a numerical algorithm is developed to efficiently determine a set of such simulations in which the…
Using methods of kinetic theory and liquid state theory we propose a description of the non-equilibrium behavior of molecular fluids which takes into account their microscopic structure and thermodynamic properties. The present work…
We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU…
A new scheme that tightly couples kinetic turbulence codes across a spatial interface is introduced. This scheme evolves from considerations of competing strategies and down-selection. It is found that the use of a composite kinetic…
We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites in 2D and 3D large scale relativisticmagnetohydrodynamic numerical simulations. This has been implemented in the PLUTO…
A general, fast, and effective approach is developed for numerical calculation of kinetic plasma dispersion relations. The plasma dispersion function is approximated by $J$-pole expansion. Subsequently, the dispersion relation is…
Magnetohydrodynamic (MHD) simulations based on the ideal MHD equations have become a powerful tool for modeling phenomena in a wide range of applications including laboratory, astrophysical, and space plasmas. In general, high-resolution…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
High-intensity laser plasma interactions create complex computational problems because they involve both fluid and kinetic regimes, which need models that maintain physical precision while keeping computational speed. The research…
Cosmological field-level inference requires differentiable forward models that solve the challenging dynamics of gas and dark matter under hydrodynamics and gravity. We propose a hybrid approach where gravitational forces are computed using…
The purpose of this paper is to study the properties of kinetic models for traffic flow described by a Boltzmann-type approach and based on a continuous space of microscopic velocities. In our models, the particular structure of the…
The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. Here, a novel, implicit, particle-in-cell based…
As compute power increases with time, more involved and larger simulations become possible. However, it gets increasingly difficult to efficiently use the provided computational resources. Especially in particle-based simulations with a…
Lattice Boltzmann method (LBM) is particularly well-suited for implementation on quantum circuits owing to its simple algebraic operations and natural parallelism. However, most quantum LBMs fix $\tau$ = 1 to avoid nonlinear collision,…
A novel parallel algorithm for matrix multiplication is presented. The hyper-systolic algorithm makes use of a one-dimensional processor abstraction. The procedure can be implemented on all types of parallel systems. It can handle…
We propose and investigate a new multi-level Monte Carlo scheme for numerical solutions of the kinetic Boltzmann equation for neutral species in edge plasmas. In particular, this method explicitly exploits a key structural property of…