Related papers: Geometric numerical integrators for Hunter-Saxton-…
This paper is dedicated to provide theta function representations of algebro-geometric solutions and related crucial quantities for the two-component Hunter-Saxton (HS2) hierarchy through studying an algebro-geometric initial value problem.…
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…
This paper describes a novel numerical model aiming at solving moving-boundary problems such as free-surface flows or fluid-structure interaction. This model uses a moving-grid technique to solve the Navier--Stokes equations expressed in…
Many applications, such as optimization, uncertainty quantification and inverse problems, require repeatedly performing simulations of large-dimensional physical systems for different choices of parameters. This can be prohibitively…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
As is known that various dynamical systems including all Hamiltonian systems preserve volume in phase space. This qualitative geometrical property of the analytical solution should be respected in the sense of Geometric Integration. This…
Many important physical systems can be described as the evolution of a Hamiltonian system, which has the important property of being conservative, that is, energy is conserved throughout the evolution. Physics Informed Neural Networks and…
We have researched the condition for symplectic discretization to preserve local boundedness for the space of 2-dimensional Hamiltonian dynamical systems in this paper.
One-dimensional scattering mediated by non-Hermitian Hamiltonians is studied. A schematic set of models is used which simulate two point interactions at a variable strength and distance. The feasibility of the exact construction of the…
While symplectic integration methods based on operator splitting are well established in many branches of science, high order methods for Hamiltonian systems that split in more than two parts have not been studied in great detail. Here, we…
We present and analyze a framework for designing symplectic neural networks (SympNets) based on geometric integrators for Hamiltonian differential equations. The SympNets are universal approximators in the space of Hamiltonian…
The numerical analysis of higher-order mixed finite-element discretizations for saddle-point problems, such as the Stokes equations, has been well-studied in recent years. While the theory and practice of such discretizations is now…
The success of symplectic integrators for Hamiltonian ODEs has led to a decades-long program of research seeking analogously structure-preserving numerical methods for Hamiltonian PDEs. In this paper, we construct a large class of such…
In this paper, for a variety of nonholonomic (reducible) Hamiltonian systems, we first give to various distributional Hamiltonian systems, by analyzing carefully the dynamics and structures of the nonholonomic Hamiltonian systems. Secondly,…
We leverage the proximal Galerkin algorithm (Keith and Surowiec, Foundations of Computational Mathematics, 2024, DOI: 10.1007/s10208-024-09681-8), a recently introduced mesh-independent algorithm, to obtain a high-order finite element…
This paper presents an integral formulation for Helmholtz problems with mixed boundary conditions. Unlike most integral equation techniques for mixed boundary value problems, the proposed method uses a global boundary charge density. As a…
We study geodesic equations for a family of right-invariant Riemannian metrics on the group of diffeomorphisms of a compact manifold. The metrics descend to Fisher's information metric on the space of smooth probability densities. The right…
We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In…
We present a geometric multilevel optimization approach that smoothly incorporates box constraints. Given a box constrained optimization problem, we consider a hierarchy of models with varying discretization levels. Finer models are…
This paper studies diagonal implicit symplectic extended Runge--Kutta--Nystr\"{o}m (ERKN) methods for solving the oscillatory Hamiltonian system $H(q,p)=\dfrac{1}{2}p^{T}p+\dfrac{1}{2}q^{T}Mq+U(q)$. Based on symplectic conditions and order…