Related papers: Metric-tensor flexible-cell algorithm for isotherm…
Isothermal-isobaric molecular dynamics simulations are used to examine the microscopic structure and some properties of water-methanol liquid mixture. The TIP4P/2005 and SPC/E water models are combined with the united atom TraPPE and the…
The Metropolis implementation of the Monte Carlo algorithm has been developed to study the equilibrium thermodynamics of many-body systems. Choosing small trial moves, the trajectories obtained applying this algorithm agree with those…
This paper presents a mixed-integer quadratic programming formulation of an existing data-driven approach to computational elasticity. This formulation is suitable for application of a standard mixed-integer programming solver, which finds…
Full relativistic simulations in three dimensions invariably develop runaway modes that grow exponentially and are accompanied by violations of the Hamiltonian and momentum constraints. Recently, we introduced a numerical method…
There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we…
We introduce a new and efficient numerical method for multicriterion optimal control and single criterion optimal control under integral constraints. The approach is based on extending the state space to include information on a "budget"…
A dynamic iteration scheme for linear differential-algebraic port-Hamil\-tonian systems based on Lions-Mercier-type operator splitting methods is developed. The dynamic iteration is monotone in the sense that the error is decreasing and no…
A general framework for performing event-driven simulations of systems with semi-flexible or rigid bodies interacting under impulsive torques and forces is outlined. Two different approaches are presented. In the first, the dynamics and…
We develop a linear method for solving the nonlinear differential equations of a lumped-parameter thermal model of a spacecraft moving in a closed orbit. Our method, based on perturbation theory, is compared with heuristic linearizations of…
In a series of recent papers we have shown how the dynamical behavior of certain classical systems can be analyzed using operators evolving according to Heisenberg-like equations of motions. In particular, we have shown that raising and…
Hamiltonian Monte Carlo (HMC) is widely used for sampling from high dimensional target distributions with densities known up to proportionality. While HMC exhibits favorable scaling properties in high dimensions, it struggles with strongly…
The Boris algorithm, a closely related variational integrator and a newly proposed filtered variational integrator are studied when they are used to numerically integrate the equations of motion of a charged particle in a non-uniform strong…
Numerical methods for developing port-Hamiltonian representations of general linear time-invariant systems are studied. The approach extends previous port-Hamiltonian characterizations to include the general non-minimal case and the case…
In this article we study intrusive uncertainty quantification schemes for systems of conservation laws with uncertainty. Standard intrusive methods lead to oscillatory solutions which sometimes even cause the loss of hyperbolicity. We…
We show that some modern geometric methods of Hamiltonian dynamics can be directly applied to the nonholonomic Heisenberg type systems. As an example we present characteristic Killing tensors, compatible Poisson brackets, Lax matrices and…
Modified Hamiltonians are used in the field of geometric numerical integration to show that symplectic schemes for Hamiltonian systems are accurate over long times. For nonlinear systems the series defining the modified Hamiltonian usually…
The self-organized Monte Carlo simulations of 2D Ising ferromagnet on the square lattice are performed. The essence of devised simulation method is the artificial dynamics consisting of the single-spin-flip algorithm of Metropolis…
We implement and investigate the numerical properties of a new family of integrators connecting both variants of the symplectic Euler schemes, and including an alternative to the classical symplectic mid-point scheme, with some additional…
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
An analytical method to compute thermodynamic properties of a given Hamiltonian system is proposed. This method combines ideas of both dynamical systems and ensemble approaches to thermodynamics, providing de facto a possible alternative to…