Related papers: A robust second-order low-rank BUG integrator base…
A new weak Galerkin (WG) finite element method for solving the second-order elliptic problems on polygonal meshes by using polynomials of boundary continuity is introduced and analyzed. The WG method is utilizing weak functions and their…
We establish rigorous \emph{a posteriori} error bounds for a space-time finite element method of arbitrary order discretising linear wave problems in second order formulation. The method combines standard finite elements in space and…
This paper establishes and analyzes a second-order accurate numerical scheme for the nonlinear partial integrodifferential equation with a weakly singular kernel. In the time direction, we apply the Crank-Nicolson method for the time…
This paper presents a new column-and-constraint generation method for two-stage robust mixed-integer programs with finite uncertainty sets. Our method combines and extends speed-up techniques used in previous column-and-constraint…
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…
Quantifying uncertainties in hyperbolic equations is a source of several challenges. First, the solution forms shocks leading to oscillatory behaviour in the numerical approximation of the solution. Second, the number of unknowns required…
This paper proposes an implicit family of sub-step integration algorithms grounded in the explicit singly diagonally implicit Runge-Kutta (ESDIRK) method. The proposed methods achieve third-order consistency per sub-step and thus the…
Embedded, or immersed, approaches have the goal of reducing to the minimum the computational costs associated with the generation of body-fitted meshes by only employing fixed, possibly Cartesian, meshes over which complex boundaries can…
Structured reduced-order modeling is a central component in the computer-aided design of control systems in which cheap-to-evaluate low-dimensional models with physically meaningful internal structures are computed. In this work, we develop…
We derive and analyze numerical methods for underdamped (kinetic) Langevin dynamics in a domain with elastic reflection at the boundary. First-order approximations are based on an Euler-type scheme incorporating collision-handling at the…
We present a new method for developing time step controllers based on a technique from the field of machine learning. This method is applicable to stable time integrators that have an embedded scheme, i.e., that have local error estimation…
Accurate and efficient simulation of fluid-structure interaction (FSI) problems remains a central challenge in computational physics. High-order discontinuous Galerkin (DG) methods offer low numerical errors and excellent scalability on…
The numerical method of dynamical low-rank approximation (DLRA) has recently been applied to various kinetic equations showing a significant reduction of the computational effort. In this paper, we apply this concept to the linear…
We combine the newly-constructed Galerkin difference basis with the energy-based discontinuous Galerkin method for wave equations in second order form. The approximation properties of the resulting method are excellent and the allowable…
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…
We introduce generalized Galerkin variational integrators, which are a natural generalization of discrete variational mechanics, whereby the discrete action, as opposed to the discrete Lagrangian, is the fundamental object. This is achieved…
The design of numerical integrators for solving stochastic dynamics with high weak order relies on tedious calculations and is subject to a high number of order conditions. The original approaches from the literature consider strong…
Following on our previous work [S. Delong and B. E. Griffith and E. Vanden-Eijnden and A. Donev, Phys. Rev. E, 87(3):033302, 2013], we develop temporal integrators for solving Langevin stochastic differential equations that arise in…
In this paper, we develop an adaptive Generalized Multiscale Discontinuous Galerkin Method (GMs-DGM) for a class of high-contrast flow problems, and derive a-priori and a-posteriori error estimates for the method. Based on the a-posteriori…
In this paper, we propose an adaptive approach, based on mesh refinement or parametric enrichment with polynomial degree adaption, for numerical solution of convection dominated equations with random input data. A parametric system emerged…