Related papers: Structure-preserving algorithms with uniform error…
We consider a class of relaxation problems mixing slow and fast variations which can describe population dynamics models or hyperbolic systems, with varying stiffness (from non-stiff to strongly dissipative), and develop a multi-scale…
Designing the topology of three-dimensional structures is a challenging problem due to its memory and time consumption. In this paper, we present a robust and efficient algorithm for solving large-scale 3D topology optimization problems.…
The rich non-linear dynamics of the coupled oscillators (under second harmonic injection) can be leveraged to solve computationally hard problems in combinatorial optimization such as finding the ground state of the Ising Hamiltonian. While…
A non-equilibrium open-dissipative neural network, such as a coherent Ising machine based on mutually coupled optical parametric oscillators, has been proposed and demonstrated as a novel computing machine for hard combinatorial…
In this paper an asymptotic expansion of the global error on the stepsize for partitioned linear multistep methods is proved. This provides a tool to analyse the behaviour of these integrators with respect to error growth with time and…
We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the…
We propose a high-order finite element method for linear fourth-order elliptic problems that is both nodally bound-preserving and mass-conservative, based on a variational inequality formulation. The method admits an equivalent strictly…
Recently, there has been an increasing interest in modelling and computation of physical systems with neural networks. Hamiltonian systems are an elegant and compact formalism in classical mechanics, where the dynamics is fully determined…
Discovering a suitable coordinate transformation for nonlinear systems enables the construction of simpler models, facilitating prediction, control, and optimization for complex nonlinear systems. To that end, Koopman operator theory offers…
Hamilton's equations of motion form a fundamental framework in various branches of physics, including astronomy, quantum mechanics, particle physics, and climate science. Classical numerical solvers are typically employed to compute the…
We survey the numerical stability of some fast algorithms for solving systems of linear equations and linear least squares problems with a low displacement-rank structure. For example, the matrices involved may be Toeplitz or Hankel. We…
We consider a dynamical system, possibly infinite dimensional or non-autonomous, with fast and slow time scales which is oscillatory with high frequencies in the fast directions. We first derive and justify the limit system of the slow…
This paper is concerned with the construction of high order schemes on irregular grids for balance laws, including a discussion of an a-posteriori error indicator based on the numerical entropy production. We also impose well-balancing on…
We introduce and analyze a family of heterogeneous multiscale methods for the numerical integration of highly oscillatory systems of delay differential equations with constant delays. The methodology suggested provides algorithms of…
In this work we demonstrate that SVD-based model reduction techniques known for ordinary differential equations, such as the proper orthogonal decomposition, can be extended to stochastic differential equations in order to reduce the…
High order energy-preserving methods for Hamiltonian systems are presented. For this aim, an energy-preserving condition of continuous stage Runge--Kutta methods is proved. Order conditions are simplified and parallelizable conditions are…
We propose a perturbation algorithm for Hamiltonian systems on a Lie algebra $\mathbb{V}$, so that it can be applied to non-canonical Hamiltonian systems. Given a Hamiltonian system that preserves a subalgebra $\mathbb{B}$ of $\mathbb{V}$,…
In this work, we propose a structure-preserving discretisation for the recently studied Cahn-Hilliard-Biot system using conforming finite elements in space and problem-adapted explicit-implicit Euler time integration. We prove that the…
We study numerically classical 1-dimensional Hamiltonian lattices involving inter-particle long range interactions that decay with distance like 1/r^alpha, for alpha>=0. We demonstrate that although such systems are generally characterized…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…