Related papers: A Tutorial on Variational Integrators
In this paper we consider the numerical approximation of systems of Boussinesq-type to model surface wave propagation. Some theoretical properties of these systems (multi-symplectic and Hamiltonian formulations, well-posedness and existence…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
Particle-based simulations of the Vlasov equation typically require a large number of particles, which leads to a high-dimensional system of ordinary differential equations. Solving such systems is computationally very expensive, especially…
Models involving hybrid systems are versatile in their application but difficult to optimize efficiently due to their combinatorial nature. This work presents a method to cope with hybrid optimal control problems which, in contrast to…
A Lagrangian formulation is constructed for particle interpretations of quantum mechanics, a well-known example of such an interpretation being the Bohm model. The advantages of such a description are that the equations for particle motion,…
It is well known that the number of particles should be scaled up to enable industrial scale simulation. The calculations are more computationally intensive when the motion of the surrounding fluid is considered. Besides the advances in…
This paper develops a structure-preserving numerical integration scheme for a class of higher-order mechanical systems. The dynamics of these systems are governed by invariant variational principles defined on higher-order tangent bundles…
In this paper, we present a new variational integrator for problems in Lagrangian mechanics. Using techniques from Galerkin variational integrators, we construct a scheme for numerical integration that converges geometrically, and is…
We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The…
We report on a lattice based algorithm, completely vectorized for molecular dynamics simulations. Its algorithmic complexity is of the order $O(N)$, where $N$ is the number of particles. The algorithm works very effectively when the…
For the calculation of the partition function $\mathcal{Z}$ of small, isolated and interacting many body systems an improvement with respect to previous formulations is presented. By including anharmonicities and employing a variational…
The paper intends to lay out the first steps towards constructing a unified framework to understand the symplectic and spectral theory of finite dimensional integrable Hamiltonian systems. While it is difficult to know what the best…
We develop a geometric framework for the numerical integration of mechanical systems evolving on manifolds. After briefly reviewing classical numerical methods and highlighting their limitations and shortcomings in non-flat (non-Euclidean)…
In this paper, we propose a variational Lagrangian scheme for a modified phase-field model, which can compute the equilibrium states for the original Allen-Cahn type model. Our discretization is based on a prescribed energy-dissipation law…
A nonlinear MPC framework is presented that is suitable for dynamical systems with sampling times in the (sub)millisecond range and that allows for an efficient implementation on embedded hardware. The algorithm is based on an augmented…
The Hamiltonian formalism is extremely elegant and convenient to mechanics problems. However, its application to the classical field theories is a difficult task. In fact, you can set one to one correspondence between the Lagrangian and…
Particle-in-cell methods couple mesh-based methods for the solution of continuum mechanics problems, with the ability to advect and evolve particles. They have a long history and many applications in scientific computing. However, they have…
Two specialized algorithms for the numerical integration of the equations of motion of a Brownian walker obeying detailed balance are introduced. The algorithms become symplectic in the appropriate limits, and reproduce the equilibrium…
A simple mathematical procedure is introduced which allows redefining in an exact way divergent integrals and limits that appear in the basic equations of classical electrodynamics with point charges. In this way all divergences are at once…
We report on a lattice based algorithm, completely vectorized for molecular dynamics simulations. Its algorithmic complexity is of the order O(N), where $N$ is the number of particles. The algorithm works very effectively when the particles…