Related papers: Semi-Lagrangian 4d, 5d, and 6d kinetic plasma simu…
The advent of 3D Gaussian Splatting has revolutionized graphics rendering by delivering high visual quality and fast rendering speeds. However, training large-scale scenes at high quality remains challenging due to the substantial memory…
3D Gaussian Splatting (3DGS) combines classic image-based rendering, pointbased graphics, and modern differentiable techniques, and offers an interesting alternative to traditional physically-based rendering. 3DGS-family models are far from…
The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, kinetic Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy…
We present a new high-order accurate Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate material dynamics (for e.g., gasses, fluids, and solids) with up to fourth-order accuracy on cubic meshes. The variables, such as…
High-performance computing is used for diverse simulations, some of which parallelize over the Message Passing Interface (MPI) with ease, whilst others may have challenges related to uniform balancing of computational load and communication…
The geometric multigrid method (GMG) is one of the most efficient solving techniques for discrete algebraic systems arising from elliptic partial differential equations. GMG utilizes a hierarchy of grids or discretizations and reduces the…
In this paper we discuss a new and very efficient implementation of high order accurate ADER discontinuous Galerkin (ADER-DG) finite element schemes on modern massively parallel supercomputers. The numerical methods apply to a very broad…
Simulating fully kinetic, two-species plasmas is computationally challenging due to the stiff multiscale dynamics of electrons and ions. While enforcing a quasi-neutral time evolution mitigates this stiffness, it requires an electric…
A tridiagonal matrix algorithm (TDMA), Pipelined-TDMA, is developed for multi-GPU systems to resolve the scalability bottlenecks caused by the sequential structure of conventional divide-and-conquer TDMA. The proposed method pipelines…
Numerical solution of reaction-diffusion equations in three dimensions is one of the most challenging applied mathematical problems. Since these simulations are very time consuming, any ideas and strategies aiming at the reduction of CPU…
This article presents new algorithms for massively parallel granular dynamics simulations on distributed memory architectures using a domain partitioning approach. Collisions are modelled with hard contacts in order to hide their…
3D Gaussian Splatting (3DGS) has garnered significant attention due to its superior scene representation fidelity and real-time rendering performance, especially for dynamic 3D scene reconstruction (\textit{i.e.}, 4D reconstruction).…
A new parallel algorithm utilizing partitioned global address space (PGAS) programming model to achieve high scalability is reported for particle tracking in direct numerical simulations of turbulent flow. The work is motivated by the…
In this paper, a physics-oriented stochastic kinetic scheme will be developed that includes random inputs from both flow and electromagnetic fields via a hybridization of stochastic Galerkin and collocation methods. Based on the BGK-type…
This paper presents a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for solving the Poisson equation in electrostatic plasma simulations using particle-in-cell (PIC) schemes. This approach…
We present a MATLAB-based framework for two- and three-dimensional fast Fourier transforms on multiple GPUs for large-scale numerical simulations using the pseudo-spectral Fourier method. The software implements two complementary multi-GPU…
In this article we introduce a novel coupled algorithm for massively parallel direct numerical simulations of electrophoresis in microfluidic flows. This multiphysics algorithm employs an Eulerian description of fluid and ions, combined…
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical…
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…
In this work, a discontinuous Galerkin scheme is employed to perform LES simulations of supersonic jet flows. A total of four simulations are performed with different meshes and order of accuracy. The number of degrees of freedom from the…