Related papers: Numerical simulations with a first order BSSN form…
The paper proposes a new, conservative fully-discrete scheme for the numerical solution of the regularised shallow water Boussinesq system of equations in the cases of periodic and reflective boundary conditions. The particular system is…
High order methods based on diagonal-norm summation by parts operators can be shown to satisfy a discrete conservation or dissipation of entropy for nonlinear systems of hyperbolic PDEs. These methods can also be interpreted as nodal…
In this paper we propose a novel numerical approach for the Boltzmann equation with uncertainties. The method combines the efficiency of classical direct simulation Monte Carlo (DSMC) schemes in the phase space together with the accuracy of…
In this paper, discontinuous Galerkin finite element methods are applied to one dimensional Rosenau equation. Theoretical results including consistency, a priori bounds and optimal error estimates are established for both semidiscrete and…
This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations (PDEs). In the paper we present a general framework for…
Einstein's system of equations in the ADM decomposition involves two subsystems of equations: evolution equations and constraint equations. For numerical relativity, one typically solves the constraint equations only on the initial time…
A numerical study of the Einstein field equations, based on the 3+1 foliation of the spacetime, is presented. A pseudo-spectral technique has been employed for simulations in vacuum, within two different formalisms, namely the…
We consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry…
In this work, we study of the algebraic-hyperbolic formulation of the Einstein constraint equations for numerically constructing initial data sets for inhomogeneous cosmological space-times with $\mathbb{T}^3$ topology. We implement a…
We propose energy-conserving discontinuous Galerkin (DG) methods for symmetric linear hyperbolic systems on general unstructured meshes. Optimal a priori error estimates of order $k+1$ are obtained for the semi-discrete scheme in one…
This paper designs a high-order positivity-preserving discontinuous Galerkin (DG) scheme for a linear hyperbolic equation. The scheme relies on augmenting the standard polynomial DG spaces with additional basis functions. The purpose of…
We implement a high-order numerical scheme for the entropy-based moment closure, the so-called M$_N$ model, for linear kinetic equations in slab geometry. A discontinuous Galerkin (DG) scheme in space along with a strong-stability…
For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are…
We prove convergence of a class of space-time discontinuous Galerkin schemes for scalar hyperbolic conservation laws. Convergence to the unique entropy solution is shown for all orders of polynomial approximation, provided strictly monotone…
We experiment with several new modifications for the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein field equations and demonstrate how these modifications affect the stability of numerical black hole evolution…
We present a new, monolithic first--order (both in time and space) BSSNOK formulation of the coupled Einstein--Euler equations. The entire system of hyperbolic PDEs is solved in a completely unified manner via one single numerical scheme…
Design of modern nanostructured semiconductor devices often calls for simulation tools capable of modeling arbitrarily-shaped multiscale geometries. In this work, to this end, a discontinuous Galerkin (DG) method-based framework is…
This paper is concerned with developing accurate and efficient numerical methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in multiple spatial dimensions. It presents a general framework…
We introduce a space-time Trefftz discontinuous Galerkin method for the first-order transient acoustic wave equations in arbitrary space dimensions, extending the one dimensional scheme of Kretzschmar et al. (2016, IMA J. Numer. Anal., 36,…
Discontinuous Galerkin methods are developed for solving the Vlasov-Maxwell system, methods that are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in…