Related papers: High-order symplectic Lie group methods on $SO(n)$…
This paper addresses the problem of efficiently computing higher-order variational integrators in simulation and trajectory optimization of mechanical systems as those often found in robotic applications. We develop $O(n)$ algorithms to…
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…
Discrete variational methods show excellent performance in numerical simulations of mechanical systems. In this paper, we adapt discrete variational integrators for the case of mechanical systems with double-bracket dissipation. In…
Computer simulations in QCD are based on the discretization of the theory on a Euclidean lattice. To compute the mean value of an observable, usually the Hybrid Monte Carlo method is applied. Here equations of motion, derived from an…
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…
Isospectral flows appear in a variety of applications, e.g. the Toda lattice in solid state physics or in discrete models for two-dimensional hydrodynamics, with the isospectral property often corresponding to mathematically or physically…
A new procedure for the construction of higher-dimensional Lie-Hamilton systems is proposed. This method is based on techniques belonging to the representation theory of Lie algebras and their realization by vector fields. The notion of…
In this article, a brief description of Discrete Mechanics and Variational Integrators which preserve the symplectic structure of the flow will be provided and a Newton-Raphson algorithm that can be used to solve implicit equations on the…
Since they were introduced in the 1990s, Lie group integrators have become a method of choice in many application areas. These include multibody dynamics, shape analysis, data science, image registration and biophysical simulations. Two…
To provide tools, especially L-operators, for use in studies of rational Yang-Baxter algebras and quantum integrable models when the Lie algebras so(N) (b_n, d_n) or sp(2n) (c_n) are the invariance algebras of their R matrices, this paper…
Exponential integrators based on contour integral representations lead to powerful numerical solvers for a variety of ODEs, PDEs, and other time-evolution equations. They are embarrassingly parallelizable and lead to global-in-time…
In this paper, we describe a geometric setting for higher-order lagrangian problems on Lie groups. Using left-trivialization of the higher-order tangent bundle of a Lie group and an adaptation of the classical Skinner-Rusk formalism, we…
We develop Lie-Poisson integrators for general Hamiltonian systems on $\mathbf{R}^{3}$ equipped with the rigid body bracket. The method uses symplectic realisation of $\mathbf{R}^{3}$ on $T^{*}\mathbf{R}^{2}$ and application of symplectic…
On the basis of the previous work by Tang \& Zhang (Appl. Math. Comput. 323, 2018, p. 204--219), in this paper we present a more effective way to construct high-order symplectic integrators for solving second order Hamiltonian equations.…
In this paper, new Levin methods are presented for calculating oscillatory integrals with algebraic and/or logarithmic singularities. To avoid singularity, the technique of singularity separation is applied and then the singular ODE…
In this paper we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for…
We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These…
This paper proposes a strong second-order two-step explicit/implicit technique with spectral orthogonal basis Galerkin finite element method for solving a two-dimensional Gray-Scott model subject to appropriate initial and boundary…
A hybridized discontinuous Galerkin method is proposed for solving 2D fractional convection-diffusion equations containing derivatives of fractional order in space on a finite domain. The Riemann-Liouville derivative is used for the spatial…
We propose a method of constructing completely integrable systems based on reduction of bihamiltonian structures. More precisely, we give an easily checkable necessary and sufficient conditions for the micro-kroneckerity of the reduction…