Related papers: Construction and analysis of higher order Galerkin…
This paper is a summary of the theory of discrete embeddings introduced in [5]. A discrete embedding is an algebraic procedure associating a numerical scheme to a given ordinary differential equation. Lagrangian systems possess a…
An integro-differential equation, modeling dynamic fractional order viscoelasticity, with a Mittag-Leffler type convolution kernel is considered. A discontinuous Galerkin method, based on piecewise constant polynomials is formulated for…
In this paper, we develop a class of high order conservative semi-Lagrangian (SL) discontinuous Galerkin (DG) methods for solving multi-dimensional linear transport equations. The methods rely on a characteristic Galerkin weak formulation,…
A wide range of implicit time integration methods, including multi-step, implicit Runge-Kutta, and Galerkin finite-time element schemes, is evaluated in the context of chaotic dynamical systems. The schemes are applied to solve the Lorenz…
A structure preserving proper orthogonal decomposition reduce-order modeling approach has been developed in [Gong et al. 2017] for the Hamiltonian system, which uses the traditional framework of Galerkin projection-based model reduction but…
Gaussian process regression is increasingly applied for learning unknown dynamical systems. In particular, the implicit quantification of the uncertainty of the learned model makes it a promising approach for safety-critical applications.…
Variational integrators are derived for structure-preserving simulation of stochastic forced Hamiltonian systems. The derivation is based on a stochastic discrete Hamiltonian which approximates a type-II stochastic generating function for…
Optimal control problems for underactuated mechanical systems can be seen as a higher-order variational problem subject to higher-order constraints (that is, when the Lagrangian function and the constraints depend on higher-order…
We study the numerical approximation of a coupled hyperbolic-parabolic system by a family of discontinuous Galerkin space-time finite element methods. The model is rewritten as a first-order evolutionary problem that is treated by the…
In this paper we introduce a procedure, based on the method of equivariant moving frames, for formulating continuous Galerkin finite element schemes that preserve the Lie point symmetries of initial value problems for ordinary differential…
We present a high order scheme for approximating kinetic equations with stiff relaxation. The objective is to provide efficient methods for solving the underlying system of conservation laws. The construction is based on several…
We thoroughly investigate Discontinuous Galerkin (DG) discretizations as time integrators for second-order oscillatory systems, considering both second-order and first-order formulations of the original problem. Key contributions include…
We present a variational integrator based on the Lobatto quadrature for the time integration of dynamical systems issued from the least action principle. This numerical method uses a cubic interpolation of the states and the action is…
We construct quasi-Monte Carlo methods to approximate the expected values of linear functionals of Galerkin discretizations of parametric operator equations which depend on a possibly infinite sequence of parameters. Such problems arise in…
The application of discontinuous Galerkin (DG) schemes to hyperbolic systems of conservation laws requires a careful interplay between space discretization, carried out with local polynomials and numerical fluxes at inter-cells, and…
Projection-based nonlinear model order reduction methods can be used to reduce simulation times for the solution of many PDE-constrained problems. It has been observed in literature that such nonlinear reduced-order models (ROMs) based on…
In this paper, we introduce two types of variational integrators, one originating from the discrete Hamilton's principle while the other from Galerkin variational approach. It turns out that these variational integrators are equivalent to…
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…
This paper presents a method to construct variational integrators for time-dependent lagrangian systems. The resulting algorithms are symplectic, preserve the momentum map associated with a Lie group of symmetries and also describe the…
We propose and compare several projection methods applied to variational integrators for degenerate Lagrangian systems, whose Lagrangian is of the form $L = \vartheta(q) \cdot \dot{q} - H(q)$ and thus linear in velocities. While previous…