Related papers: Conservative Integrators for Many-body Problems
Parameterization (closure) schemes in numerical weather and climate prediction models account for the effects of physical processes that cannot be resolved explicitly by these models. Methods for finding physical parameterization schemes…
We propose an efficient way of solving optimal control problems for rigid-body systems on the basis of inverse dynamics and the multiple-shooting method. We treat all variables, including the state, acceleration, and control input torques,…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
We extend the self-consistent Green's functions formalism to take into account three-body interactions. We analyze the perturbative expansion in terms of Feynman diagrams and define effective one- and two-body interactions, which allows for…
We consider the solution of the fully kinetic (including electrons) Vlasov-Amp\`ere system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian…
This paper introduces a new symbolic-numeric strategy for finding semidiscretizations of a given PDE that preserve multiple local conservation laws. We prove that for one spatial dimension, various one-step time integrators from the…
We propose a new approach for solving systems of conservation laws that admit a variational formulation of the time-discretized form, and encompasses the p-system or the system of elastodynamics. The approach consists of using constrained…
Multi-body mechanical systems have rich internal dynamics, whose solutions can be exploited as efficient control targets. Yet, solutions non-trivially depend on system parameters, obscuring feasible properties for use as target…
We design a variational asymptotic preserving scheme for the Vlasov-Poisson-Fokker-Planck system with the high field scaling, which describes the Brownian motion of a large system of particles in a surrounding bath. Our scheme builds on an…
In this paper we consider the spatial semi-discretization of conservative PDEs. Such finite dimensional approximations of infinite dimensional dynamical systems can be described as flows in suitable matrix spaces, which in turn leads to the…
We develop a general framework for designing conservative numerical methods based on summation by parts operators and split forms in space, combined with relaxation Runge-Kutta methods in time. We apply this framework to create new classes…
Interacting spin systems are of fundamental relevance in different areas of physics, as well as in quantum information science, and biology. These spin models represent the simplest, yet not fully understood, manifestation of quantum…
It is well known, thanks to Lax-Wendroff theorem, that the local conservation of a numerical scheme for a conservative hyperbolic system is a simple and systematic way to guarantee that, if stable, a scheme will provide a sequence of…
We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping,…
We are interested in the numerical approximation of discontinuous solutions in non conservative hyperbolic systems. An extension to second-order of a new strategy based on in-cell discontinuous reconstructions to deal with this challenging…
Periodically driven nonequilibrium many-body systems are interesting because they have a quasi-energy spectra, which can be tailored by controlling the external driving fields. We derive the general spectral representation of retarded Green…
A new method is proposed to numerically integrate a dynamical system on a manifold such that the trajectory stably remains on the manifold and preserves first integrals of the system. The idea is that given an initial point in the manifold…
Slow-fast dynamical systems, i.e., singularly or non-singularly perturbed dynamical systems possess slow invariant manifolds on which trajectories evolve slowly. Since the last century various methods have been developed for approximating…
The exact equations of motion for microscopic density of classical many-body system with account of inter-particle retarded interactions are derived. It is shown that interactions retardation leads to irreversible behaviour of many-body…
We propose a locally conservative enriched Galerkin scheme that preserves the physical bounds for an elliptic problem. To this end, we use a substantial over-penalization of the discrete solution's jumps to obtain optimal convergence. To…