Related papers: Prolongation-Collocation Variational Integrators
We present the symplectic algorithm in the Lagrangian formalism for the Hamiltonian systems by virtue of the noncommutative differential calculus with respect to the discrete time and the Euler--Lagrange cohomological concepts. We also show…
It is well known that symplectic integrators lose their near energy preservation properties when variable step sizes are used. The most common approach to combine adaptive step sizes and symplectic integrators involves the Poincar\'e…
We have applied a collocation approach to obtain the numerical solution to the stationary Schr\"odinger equation for systems of coupled oscillators. The dependence of the discretized Hamiltonian on scale and angle parameters is exploited to…
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…
In recent years, many exceptional orthogonal polynomials (EOP) were introduced and used to construct new families of 1D exactly solvable quantum potentials, some of which are shape invariant. In this paper, we construct from Hermite and…
We study multiple orthogonal polynomials exploiting their explicit determinantal representation in terms of moments. Our reasoning follows that applied to solve the Hermite-Pad\'{e} approximation and interpolation problems. We study also…
This article proposes a new numerical algorithm for second order elliptic equations in non-divergence form. The new method is based on a discrete weak Hessian operator locally constructed by following the weak Galerkin strategy. The…
We develop a variational integrator for the shallow-water equations on a rotating sphere. The variational integrator is built around a discretization of the continuous Euler-Poincar\'{e} reduction framework for Eulerian hydrodynamics. We…
We propose a novel algorithmic method for constructing invariant variational schemes of systems of ordinary differential equations that are the Euler-Lagrange equations of a variational principle. The method is based on the invariantization…
We propose a third-order numerical integrator based on the Neumann series and the Filon quadrature, designed mainly for highly oscillatory partial differential equations. The method can be applied to equations that exhibit small or moderate…
Variational integrators are well-suited for simulation of mechanical systems because they preserve mechanical quantities about a system such as momentum, or its change if external forcing is involved, and holonomic constraints. While they…
In this paper, we derive a variational integrator for certain highly oscillatory problems in mechanics. To do this, we take a new approach to the splitting of fast and slow potential forces: rather than splitting these forces at the level…
In recent years, two important techniques for geometric numerical discretization have been developed. In computational electromagnetics, spatial discretization has been improved by the use of mixed finite elements and discrete differential…
We address the problem of constructing numerical integrators for nonholonomic Lagrangian systems that enjoy appropriate discrete versions of the geometric properties of the continuous flow, including the preservation of energy. Building on…
The equations for the critical points of the action functional defined by a Lagrangian depending on higher-order derivatives of admissible curves on a Lie algebroid are found. The relation with Euler-Poincar\'e and Lagrange Poincar\'e type…
In this work, we develop a fully implicit Hybrid High-Order algorithm for the Cahn-Hilliard problem in mixed form. The space discretization hinges on local reconstruction operators from hybrid polynomial unknowns at elements and faces. The…
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…
A discrete version of Lagrangian reduction is developed in the context of discrete time Lagrangian systems on $G\times G$, where $G$ is a Lie group. We consider the case when the Lagrange function is invariant with respect to the action of…
The paper suggests a Hamiltonian formulation for delay ordinary differential equations (DODEs). Such equations are related to DODEs with a Lagrangian formulation via a delay analog of the Legendre transformation. The Hamiltonian delay…
We present geometric numerical integrators for contact flows that stem from a discretization of Herglotz' variational principle. First we show that the resulting discrete map is a contact transformation and that any contact map can be…