Related papers: MultiWave: A computational lab for adaptive numeri…
We study a recent timestep adaptation technique for hyperbolic conservation laws. The key tool is a space-time splitting of adjoint error representations for target functionals due to S\"uli and Hartmann. It provides an efficient choice of…
We derive formulas for the Coulomb matrix within the full-potential linearized augmented-plane-wave (FLAPW) method. The Coulomb matrix is a central ingredient in implementations of many-body perturbation theory, such as the Hartree-Fock and…
Discontinuous Galerkin finite element schemes exhibit attractive features for accurate large-scale wave-propagation simulations on modern parallel architectures. For many applications, these schemes must be coupled with non-reflective…
This is the fourth installment in our series on implementing the discontinuous Galerkin (DG) method as an open source MATLAB /GNU Octave toolbox. Similarly to its predecessors, this part presents new features for application developers…
Weight-adjusted inner products are easily invertible approximations to weighted $L^2$ inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time-domain method for wave propagation…
It is well-known that reliable and efficient domain truncation is crucial to accurate numerical solution of most wave propagation problems. The perfectly matched layer (PML) is a method which, when stable, can provide a domain truncation…
We develop an approach for simulating acousto-elastic wave phenomena, including scattering from fluid-solid boundaries, where the solid is allowed to be anisotropic, with the Discontinuous Galerkin method. We use a coupled first-order…
Next generation computer architectures will include order of magnitude more intra-node parallelism; however, many application programmers have a difficult time keeping their codes current with the state-of-the-art machines. In this context,…
The $M_1$ minimum entropy moment system is a system of hyperbolic balance laws that approximates the radiation transport equation, and has many desirable properties. Among them are symmetric hyperbolicity, entropy decay, moment…
In this article, we present Bembel, the C++ library featuring higher order isogeometric Galerkin boundary element methods for Laplace, Helmholtz, and Maxwell problems. Bembel is compatible with geometries from the Octave NURBS package and…
Stiff hyperbolic balance laws exhibit large spectral gaps, especially if the relaxation term significantly varies in space. Using examples from rarefied gases and the general form of the underlying balance law model, we perform a detailed…
We first investigate the structure of the systems derived from the gPC based stochastic Galerkin method for the nonlinear hyperbolic systems with random inputs. This method adopts a generalized Polynomial Chaos (gPC) approximations in the…
The development and application of the Discontinuous Galerkin (DG) method have attracted great attention in computational fluid dynamics (CFD) com- munity in the past decades. The underlying reason for such an intensive investigation is due…
We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra space.…
For hyperbolic conservation laws, traditional methods and physics-informed neural networks (PINNs) often encounter difficulties in capturing sharp discontinuities and maintaining temporal consistency. To address these challenges, we…
In this paper, we propose a class of adaptive multiresolution (also called adaptive sparse grid) ultra-weak discontinuous Galerkin (UWDG) methods for solving some nonlinear dispersive wave equations including the Korteweg-de Vries (KdV)…
Merlin++ is a C++ charged-particle tracking library developed for the simulation and analysis of complex beam dynamics within high energy particle accelerators. Accurate simulation and analysis of particle dynamics is an essential part of…
A coupling approach is presented to combine a wave-based method to the standard finite element method. This coupling methodology is presented here for the Helmholtz equation but it can be applied to a wide range of wave propagation…
We present the design and implementation of a \texttt{C++} class for reliability analysis of multi-state systems using an algebraic approach based on monomial ideals. The class is implemented within the open-source \verb|CoCoALib| library…
The two-fluid plasma model has a wide range of timescales which must all be numerically resolved regardless of the timescale on which plasma dynamics occurs. The answer to solving numerically stiff systems is generally to utilize…