Related papers: FELIX-1.0: A finite element solver for the time de…
An operator-splitting finite element scheme for the time-dependent, high-dimensional radiative transfer equation is presented in this paper. The streamline upwind Petrov-Galerkin finite element method and discontinuous Galerkin finite…
We introduce a max-plus analogue of the Petrov-Galerkin finite element method, to solve finite horizon deterministic optimal control problems. The method relies on a max-plus variational formulation, and exploits the properties of…
In this paper we study the numerical method and the convergence for solving the time-dependent Maxwell-Schr\"{o}dinger equations under the Lorentz gauge. An alternating Crank-Nicolson finite element method for solving the problem is…
This paper introduces JAX-FEM, an open-source differentiable finite element method (FEM) library. Constructed on top of Google JAX, a rising machine learning library focusing on high-performance numerical computing, JAX-FEM is implemented…
This article establishes a discrete maximum principle (DMP) for the approximate solution of convection-diffusion-reaction problems obtained from the weak Galerkin finite element method on nonuniform rectangular partitions. The DMP analysis…
A numerical method is proposed to solve the full-Eulerian time-dependent Vlasov-Poisson system in high dimension. The algorithm relies on the construction of a tensor decomposition of the solution whose rank is adapted at each time step.…
The third paper in our series on open source MATLAB / GNU Octave implementation of the discontinuous Galerkin (DG) method(s) focuses on a hybridized formulation. The main aim of this ongoing work is to develop rapid prototyping techniques…
This paper introduces a new weak Galerkin (WG) finite element method for second order elliptic equations on polytopal meshes. This method, called WG-FEM, is designed by using a discrete weak gradient operator applied to discontinuous…
The Galerkin difference (GD) basis is a set of continuous, piecewise polynomials defined using a finite difference like grid of degrees of freedom. The one dimensional GD basis functions are naturally extended to multiple dimensions using…
Discretizing a solution in the Fourier domain rather than the time domain presents a significant advantage in solving transport problems that vary smoothly and periodically in time, such as cardiorespiratory flows. The finite element…
We propose and analyze a time-stepping discontinuous Petrov-Galerkin method combined with the continuous conforming finite element method in space for the numerical solution of time-fractional subdiffusion problems. We prove the existence,…
The Discontinuous Galerkin time-domain method is well suited for adaptive algorithms to solve the time-domain Maxwell's equations and depends on robust and economically computable drivers. Adaptive algorithms utilize local indicators to…
This paper proposes an efficient ADER (Arbitrary DERivatives in space and time) discontinuous Galerkin (DG) scheme to directly solve the Hamilton-Jacobi equation. Unlike multi-stage Runge-Kutta methods used in the Runge-Kutta DG (RKDG)…
This work presents GALAEXI as a novel, energy-efficient flow solver for the simulation of compressible flows on unstructured meshes leveraging the parallel computing power of modern Graphics Processing Units (GPUs). GALAEXI implements the…
We study a space-time finite element approach for the nonhomogeneous wave equation using a continuous time Galerkin method. We present fully implicit examples in 1+1, 2+1, and 3+1 dimensions using linear quadrilateral, hexahedral, and…
We design an arbitrary-order free energy satisfying discontinuous Galerkin (DG) method for solving time-dependent Poisson-Nernst-Planck systems. Both the semi-discrete and fully discrete DG methods are shown to satisfy the corresponding…
We present a new line-based discontinuous Galerkin (DG) discretization scheme for first- and second-order systems of partial differential equations. The scheme is based on fully unstructured meshes of quadrilateral or hexahedral elements,…
The Vlasov-Poisson equations describe the evolution of a collisionless plasma, represented through a probability density function (PDF) that self-interacts via an electrostatic force. One of the main difficulties in numerically solving this…
This study presents the Fourier-Gegenbauer Integral-Galerkin (FGIG) method, a novel and efficient numerical framework for solving the one-dimensional advection-diffusion equation with periodic boundary conditions. The FGIG method uniquely…
The discontinuous Galerkin (DG) method is an established method for computing approximate solutions of partial differential equations in many applications. Unlike continuous finite elements, in DG methods, numerical fluxes are used to…