Related papers: Stabilization in relation to wavenumber in HDG met…
Multi-scale wave propagation problems are computationally costly to solve by traditional techniques because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and…
High-order Discontinuous Galerkin (DG) methods offer excellent accuracy for turbulent flow simulations, especially when implemented on GPU-oriented architectures that favor very high polynomial orders. On modern GPUs, high-order polynomial…
Wave propagation problems for heterogeneous media are known to have many applications in physics and engineering. Recently, there has been an increasing interest in stochastic effects due to the uncertainty, which may arise from impurities…
We propose and analyze a hybridized discontinuous Galerkin (HDG) method for the spherically symmetric Einstein--scalar system in Bondi gauge. After rewriting the model as a local first-order PDE--ODE system by introducing suitable scaled…
This paper introduces a numerical scheme for time harmonic Maxwell's equations by using weak Galerkin (WG) finite element methods. The WG finite element method is based on two operators: discrete weak curl and discrete weak gradient, with…
We develop a high order accurate numerical method for solving the elastic wave equation in second-order form. We hybridize the computationally efficient Cartesian grid formulation of finite differences with geometrically flexible…
In this paper, the authors devise a new discretization scheme for div-curl systems defined in connected domains with heterogeneous media by using the weak Galerkin finite element method. Two types of boundary value problems are considered…
The $hp$ local discontinuous Galerkin (LDG) method proposed by Castillo et al. [Math. Comp.,~71 (238): 455-478, 2002] has been shown to be an efficient approach for solving convection-diffusion equations. However, theoretical analysis…
This paper introduces a high order numerical framework for efficient and robust simulation of compressible flows. To address the inefficiencies of standard hybridized discontinuous Galerkin (HDG) methods in large scale settings, we develop…
The Sagdeev pseudo-potential approach has been employed extensively in theoretical studies to determine large-amplitude (fully) nonlinear solutions in a variety of multi-species plasmas. Although these solutions are repeatedly considered as…
This paper presents an $hp$ a posteriori error analysis for the 2D Helmholtz equation that is robust in the polynomial degree $p$ and the wave number $k$. For the discretization, we consider a discontinuous Galerkin formulation that is…
We study an embedded Trefftz discontinuous Galerkin method for the Helmholtz equation. The method starts from a polynomial DG space and enforces the Trefftz property through local constraints, avoiding an explicit construction of Trefftz…
Solving time-harmonic wave propagation problems in the frequency domain within heterogeneous media poses significant mathematical and computational challenges, particularly in the high-frequency regime. Among the available numerical…
The propagation of electromagnetic waves in general media is modeled by the time-dependent Maxwell's partial differential equations (PDEs), coupled with constitutive laws that describe the response of the media. In this work, we focus on…
We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a…
We study the unique solvability of the discretized Helmholtz problem with Robin boundary conditions using a conforming Galerkin $hp$-finite element method. Well-posedness of the discrete equations is typically investigated by applying a…
An error analysis of a mixed discontinuous Galerkin (DG) method with Brezzi numerical flux for the time-harmonic Maxwell equations with minimal smoothness requirements is presented. The key difficulty in the error analysis for the DG method…
The accurate numerical simulation of turbulent incompressible flows is a challenging topic in computational fluid dynamics. For discretisation methods to be robust in the under-resolved regime, mass conservation as well as energy stability…
An hyperelastic biphasic model is presented. For slow-draining problems (permeability less than 1\times10-2 mm4 N-1 s-1), numerical instabilities in the form of non-physical oscillations in the pressure field are observed in 3D problems…
Collocation boundary element methods for integral equations are easier to implement than Galerkin methods because the elements of the discretization matrix are given by lower-dimensional integrals. For that same reason, the matrix assembly…