Related papers: Explicit methods in extended phase space for insep…
Dissipation can be represented in Hamiltonian mechanics in an extended phase space as a symplectic process. The method uses an auxiliary variable which represents the excitation of unresolved dynamics and a Hamiltonian for the interaction…
As part of our development of a computer code to perform 3D `constrained evolution' of Einstein's equations in 3+1 form, we discuss issues regarding the efficient solution of elliptic equations on domains containing holes (i.e., excised…
We propose an explicit iterative leap-frog discontinuous Galerkin method for time-domain Maxwell's equations in anisotropic materials and derive its convergence properties. The a priori error estimates are illustrated by numerical means in…
A practical method to solve cut-off Coulomb problems of two-cluster systems in the momentum space is given. When a sharply cut-off Coulomb force with a cut-off radius $\rho$ is introduced at the level of constituent particles, two-cluster…
We consider the numerical solution of partial differential equations with coefficients that are strongly heterogeneous in space. We provide an overview of higher-order localized orthogonal decomposition (LOD) methods for the elliptic…
Many experimental techniques aim at determining the Hamiltonian of a given system. The Hamiltonian describes the system's evolution in the absence of dissipation, and is often central to control or interpret an experiment. Here, we…
By introducing an integration factor to the differential one-form of contact dynamics, equations of motion are derived variationally, and contact Poisson bracket and contact Lagrangian are formulated. Discrete symplectic integrator, named…
The recent development of spectral method has been praised for its high-order convergence in simulating complex physical problems. The combination of embedded boundary method and spectral method becomes a mainstream way to tackle…
Using the framework of metriplectic systems on $\R^n$ we will describe a constructive geometric method to add a dissipation term to a Hamilton-Poisson system such that any solution starting in a neighborhood of a nonlinear stable…
We consider the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients (methods with real coefficients of order greater than two necessarily have negative…
The dilation method is a practical way to experimentally simulate non-Hermitian, especially $\cal PT$-symmetric quantum systems. However, the time-dependent dilation problem cannot be explicitly solved in general. In this paper, we present…
We prove in this work that the well-known lasso problem can be solved exactly without homotopy using novel differential inclusions techniques. Specifically, we show that a selection principle from the theory of differential inclusions…
This article introduces the Modified Parameterized Leapfrog Hamiltonian Monte Carlo (MPL-HMC) method, a novel extension of HMC addressing key limitations through tunable integration parameters $\alpha(\delta t)$ and $\beta(\delta t)$,…
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 present an example of an integrable Hamiltonian system with scalar potential in the three-dimensional Euclidean space whose integrals of motion are quadratic polynomials in the momenta, yet its Hamilton-Jacobi / Schrodinger equation…
This paper presents a hybrid numerical method to solve efficiently a class of highly anisotropic elliptic problems. The anisotropy is aligned with one coordinate-axis and its strength is described by a parameter $\eps \in (0,1]$, which can…
We study the finite element approximation of linear second-order elliptic partial differential equations in nondivergence form with highly heterogeneous diffusion and drift coefficients. A generalized Cordes condition is imposed to…
An adequate characterization of the dynamics of Hamiltonian systems at physically relevant scales has been largely lacking. Here we investigate this fundamental problem and we show that the finite-scale Hamiltonian dynamics is governed by…
We investigate the connection between the linear harmonic oscillator equation and some classes of second order nonlinear ordinary differential equations of Li\'enard and generalized Li\'enard type, which physically describe important…
In this paper, we develop a new type of accelerated algorithms to solve some classes of maximally monotone equations as well as monotone inclusions. Instead of using Nesterov's accelerating approach, our methods rely on a so-called…