Related papers: Mimetic Metrics for the DGSEM
Continuum computational kinetic plasma models evolve the distribution function of a plasma species $f_s$ on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence,…
We present a hybrid continuous and discontinuous Galerkin spectral element approximation that leverages the advantages of each approach. The continuous Galerkin approximation is used on interior element faces where the equation properties…
Entropy stable schemes ensure that physically meaningful numerical solutions also satisfy a semi-discrete entropy inequality under appropriate boundary conditions. In this work, we describe a discretization of viscous terms in the…
In this paper we present a new high order semi-implicit DG scheme on two-dimensional staggered triangular meshes applied to different nonlinear systems of hyperbolic conservation laws such as advection-diffusion models, incompressible…
The immersed interface method (IIM) for models of fluid flow and fluid-structure interaction imposes jump conditions that capture stress discontinuities generated by forces that are concentrated along immersed boundaries. Most prior work…
Computational electrodynamics (CED), the numerical solution of Maxwell's equations, plays an incredibly important role in several problems in science and engineering. High accuracy solutions are desired, and the discontinuous Galerkin (DG)…
Some hyperbolic systems are known to include implicit preservation of differential constraints: these are for example the time conservation of the curl or the divergence of a vector that appear as an implicit constraint. In this article, we…
In this work, we develop a novel numerical scheme to solve the classical Keller--Segel (KS) model which simultaneously preserves its intrinsic mathematical structure and achieves optimal accuracy. The model is reformulated into a gradient…
Unstructured grid ocean models are advantageous for simulating the coastal ocean and river-estuary-plume systems. However, unstructured grid models tend to be diffusive and/or computationally expensive which limits their applicability to…
We propose a structure-preserving parametric finite element method (SP-PFEM) for discretizing the surface diffusion of a closed curve in two dimensions (2D) or surface in three dimensions (3D). Here the "structure-preserving" refers to…
In this work a non-conservative balance law formulation is considered that encompasses the rotating, compressible Euler equations for dry atmospheric flows. We develop a semi-discretely entropy stable discontinuous Galerkin method on…
Accurate representation of vertical turbulent fluxes is crucial for numerical ocean modelling, both in global and coastal applications. The state-of-the-art approach is to use two-equation turbulence closure models which introduces two…
Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not…
In this study, we propose a unified, general framework for the direct discontinuous Galerkin methods. In the new framework, the antiderivative of the nonlinear diffusion matrix is not needed. This allows a simple definition of the numerical…
This paper deals with the numerical resolution of kinetic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. We focus on the kinetic model considered in [18, 17] where alignment is…
A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in our previous computational work. This paper presents a numerical…
High-order discontinuous Galerkin spectral element methods (DGSEM) have received growing attention and development, especially in the regime of computational fluid dynamics in recent years. The inherent flexibility of the discontinuous…
Employing physically-consistent numerical methods is an important step towards attaining robust and accurate numerical simulations. When addressing compressible flows, in addition to preserving kinetic energy at a discrete level, as done in…
The Ericksen model for nematic liquid crystals couples a director field with a scalar degree of orientation variable, and allows the formation of various defects with finite energy. We propose a simple but novel finite element approximation…
In this work, we develop a localized numerical scheme with low regularity requirements for solving time-fractional integro-differential equations. First, a fully discrete numerical scheme is constructed. Specifically, for temporal…