Related papers: Sapphire++: A particle transport code combining a …
The discontinuous Galerkin finite element method (DG-FEM) is successfully applied to treat a broad variety of transport problems numerically. In this work, we use the full capacity of the DG-FEM to solve the radiative transfer equation in…
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics…
In highly diffusion regimes when the mean free path $\varepsilon$ tends to zero, the radiative transfer equation has an asymptotic behavior which is governed by a diffusion equation and the corresponding boundary condition. Generally, a…
We present an approach for simulating the collisional evolution of spherical isotropic stellar systems based on the one-dimensional Fokker-Planck equation. A novel aspect is that we use the phase volume as the argument of the distribution…
This paper develops a fully discrete Fourier spectral Galerkin (FSG) method for the fractional Zakharov--Kuznetsov (fZK) equation posed on a two-dimensional periodic domain. The equation generalizes the classical ZK model by replacing the…
High-dimensional transport equations frequently occur in science and engineering. Computing their numerical solution, however, is challenging due to its high dimensionality. In this work we develop an algorithm to efficiently solve the…
We discuss the development, analysis, implementation, and numerical assessment of a spectral method for the numerical simulation of the three-dimensional Vlasov-Maxwell equations. The method is based on a spectral expansion of the velocity…
In this work, we extend the solid harmonics derivation, which was used by Ackroyd et al to derive the steady-state SP$_N$ equations, to transient problems. The derivation expands the angular flux in ordinary surface harmonics but uses…
The high-order accurate continuous Galerkin finite element method offers attractive computational efficiency for computational fluid dynamics. A challenge is however spurious oscillations which result for convection dominated flows over…
The Fokker-Planck equation describing the transport of energetic particles interacting with turbulence is difficult to solve analytically. Numerical solutions are of course possible but they are not always useful for applications. In the…
We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a…
We propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite volume (FV) method. The strategy relies on a simple operator splitting…
This paper deals with uncertainty propagation of general stochastic hybrid systems (GSHS) where the continuous state space is a compact Lie group. A computational framework is proposed to solve the Fokker-Planck (FP) equation that describes…
This paper is concerned with structure-preserving numerical approximations for a class of nonlinear nonlocal Fokker-Planck equations, which admit a gradient flow structure and find application in diverse contexts. The solutions,…
Discontinuous Galerkin (DG) methods are widely adopted to discretize the radiation transport equation (RTE) with diffusive scalings. One of the most important advantages of the DG methods for RTE is their asymptotic preserving (AP)…
We consider a class of time dependent second order partial differential equations governed by a decaying entropy. The solution usually corresponds to a density distribution, hence positivity (non-negativity) is expected. This class of…
The numerical solution of time-dependent radiative transfer problems is challenging, both, due to the high dimension as well as the anisotropic structure of the underlying integro-partial differential equation. In this paper we propose a…
The semi-Lagrangian discontinuous Galerkin method, coupled with a splitting approach in time, has recently been introduced for the Vlasov--Poisson equation. Since these methods are conservative, local in space, and able to limit numerical…
This spreading of prion proteins is at the basis of brain neurodegeneration. This paper deals with the numerical modelling of the misfolding process of $\alpha$-synuclein in Parkinson's disease. We introduce and analyze a discontinuous…
In this work, we propose a new Galerkin-Petrov method for the numerical solution of the classical spatially homogeneous Boltzmann equation. This method is based on an approximation of the distribution function by associated Laguerre…