Related papers: Explicit methods in extended phase space for insep…
We consider solutions of the $2\times 2$ matrix Hamiltonian of physical systems within the context of the asymptotic iteration method. Our technique is based on transformation of the associated Hamiltonian in the form of the first order…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
Recent well-posedness results have identified the Hardy space $L^2_+$ as the natural phase space for continuum Calogero-Moser models, both focusing and defocusing, on the line and on the torus. In this paper, we introduce a symplectic form…
The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the…
The present work provides a comprehensive study of symmetric-conjugate operator splitting methods in the context of linear parabolic problems and demonstrates their additional benefits compared to symmetric splitting methods. Relevant…
We present a new higher-order accurate finite difference explicit jump Immersed Interface Method (HEJIIM) for solving two-dimensional elliptic problems with singular source and discontinuous coefficients in the irregular region on a compact…
In this work we propose a method to perform optimization on manifolds. We assume to have an objective function $f$ defined on a manifold and think of it as the potential energy of a mechanical system. By adding a momentum-dependent kinetic…
An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as…
We develop an approach to affine symplectic invariant geometry of Lagrangian surfaces by the method of moving frames. The fundamental invariants of elliptic Lagrangian immersions in affine symplectic four-space are derived together with…
Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical…
It is common in classical mechanics to encounter systems whose Hamiltonian $H$ is the sum of an often exactly integrable Hamiltonian $H_0$ and a small perturbation $\epsilon H_1$ with $\epsilon\ll1$. Such near-integrability can be exploited…
In this article we introduce a low order implicit symplectic integrator designed to follow the Hamiltonian flow as close as possible. This integrator is obtained by the method of Liouvillian forms and does not require particular hypotheses…
In this paper, we propose two arbitrary order eXtended hybridizable Discontinuous Galerkin (X-HDG) methods for second order elliptic interface problems in two and three dimensions. The first X-HDG method applies to any piecewise $C^2$…
Solving quaternion kinematical differential equations is one of the most significant problems in the automation, navigation, aerospace and aeronautics literatures. Most existing approaches for this problem neither preserve the norm of…
A new immersed finite element (IFE) method is developed for second-order elliptic problems with discontinuous diffusion coefficient. The IFE space is constructed based on the rotated Q1 nonconforming finite elements with the integral-value…
In this introductory work I will present the Finite Difference method for hyperbolic equations, focusing on a method which has second order precision both in time and space (the so-called staggered leapfrog method) and applying it to the…
This paper introduces weighted finite difference methods for numerically solving dispersive evolution equations with solutions that are highly oscillatory in both space and time. We consider a semiclassically scaled cubic nonlinear…
In this paper, we consider a Bayesian inverse problem modeled by elliptic partial differential equations (PDEs). Specifically, we propose a data-driven and model-based approach to accelerate the Hamiltonian Monte Carlo (HMC) method in…
We leverage the connections between nonexpansive maps, monotone Lipschitz operators, and proximal mappings to obtain near-optimal (i.e., optimal up to poly-log factors in terms of iteration complexity) and parameter-free methods for solving…
It has been recently argued that near-integrable nonautonomous one-degree-of-freedom Hamiltonian systems are constrained by KAM theory even when the time-dependent (nonintegrable) part of the Hamiltonian is given in the form of a…