Related papers: Construction and analysis of higher order Galerkin…
This overview is devoted to splitting methods, a class of numerical integrators intended for differential equations that can be subdivided into different problems easier to solve than the original system. Closely connected with this class…
We propose a block finite difference, error inhibiting scheme that is fourth-order accurate for short to moderate times and has a six-order convergence rate for long times. This scheme outperforms the standard fourth-order Finite Difference…
We study higher-order space-time variational discretisations for modeling complex processes in porous media that include fluid and structure interactions which are of fundamental importance in many engineering fields with applications in…
This work proposes a suite of numerical techniques to facilitate the design of structure-preserving integrators for nonlinear dynamics. The celebrated LaBudde-Greenspan integrator and various energy-momentum schemes adopt a difference…
A weak Galerkin (WG) finite element method without stabilizers was introduced in [J. Comput. Appl. Math., 371 (2020). arXiv:1906.06634] on polytopal mesh. Then it was improved in [arXiv:2008.13631] with order one superconvergence. The goal…
In this paper, we construct an efficient numerical scheme for full-potential electronic structure calculations of periodic systems. In this scheme, the computational domain is decomposed into a set of atomic spheres and an interstitial…
In this paper, we develop a high order structure-preserving local discontinuous Galerkin (DG) scheme for the compressible self-gravitating Euler equations, which pose great challenges due to the presence of time-dependent gravitational…
This paper introduces an extension of the Morley element for approximating solutions to biharmonic equations. Traditionally limited to piecewise quadratic polynomials on triangular elements, the extension leverages weak Galerkin finite…
This work focuses on the conservation of quantities such as Hamiltonians, mass, and momentum when solution fields of partial differential equations are approximated with nonlinear parametrizations such as deep networks. The proposed…
The discontinuous Galerkin time-stepping method has many advantageous properties for solving parabolic equations. However, it requires the solution of a large nonsymmetric system at each time-step. This work develops a fully robust and…
We develop 3rd order maximum-principle-satisfying direct discontinuous Galerkin methods [8, 9, 19, 21] for convection diffusion equations on unstructured triangular mesh. We carefully calculate the normal derivative numerical flux across…
The high-order accurate continuous Galerkin finite element method offers attractive computational efficiency for computational fluid dynamics. A challenge is however spurious oscillations which result for convection dominated flows over…
Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved…
This paper focuses on the numerical approximation of the solutions of non-local conservation laws in one space dimension. These equations are motivated by two distinct applications, namely a traffic flow model in which the mean velocity…
We develop efficient numerical integration methods for computing an integral whose integrand is a product of a smooth function and the Gaussian function with a small standard deviation. Traditional numerical integration methods applied to…
We consider the systematic numerical approximation of Biot's quasistatic model for the consolidation of a poroelastic medium. Various discretization schemes have been analysed for this problem and inf-sup stable finite elements have been…
The problem of optimal recovering high-order mixed derivatives of bivariate functions with finite smoothness is studied. Based on the truncation method, an algorithm for numerical differentiation is constructed, which is order-optimal both…
We introduce and analyze a class of Galerkin-collocation discretization schemes in time for the wave equation. Its conceptual basis is the establishment of a direct connection between the Galerkin method for the time discretization and the…
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…
Multibody dynamics simulators are an important tool in many fields, including learning and control for robotics. However, many existing dynamics simulators suffer from inaccuracies when dealing with constrained mechanical systems due to…