Related papers: Structure and structure-preserving algorithms for …
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…
In this paper, based on the weak form of the Hamiltonian formulation of the regularized long-wave equation and a novel approach of transforming the original Hamiltonian energy into a quadratic functional, a fully implicit and three…
We introduce a new class of parametricization structure-preserving partitioned Runge-Kutta ($\alpha$-PRK) methods for Hamiltonian systems with holonomic constraints. When the scalar parameter $\alpha=0$, the methods are reduced to the usual…
In this article we apply a discrete action principle for the Vlasov--Maxwell equations in a structure-preserving particle-field discretization framework. In this framework the finite-dimensional electromagnetic potentials and fields are…
The hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell…
We discuss a version of Hamiltonian (2+1)-dimensional dynamics, in which one allows nonvanishing Poisson brackets also between the coordinates, and between the momenta. The resulting equations of motion are not any more derivable from a…
We consider a class of multiscale parabolic problems with diffusion coefficients oscillating in space at a possibly small scale $\varepsilon$. Numerical homogenization methods are popular for such problems, because they capture efficiently…
We develop structure-preserving numerical methods for the compressible Euler equations, employing potential temperature as a prognostic variable. We construct three numerical fluxes designed to ensure the conservation of entropy and total…
We propose new algorithms for numerical integration of the equations of motion for classical spin systems with fixed spatial site positions. The algorithms are derived on the basis of a mid-point scheme in conjunction with the multiple time…
We propose explicit symplectic integrators of molecular dynamics (MD) algorithms for rigid-body molecules in the canonical and isothermal-isobaric ensembles. We also present a symplectic algorithm in the constant normal pressure and lateral…
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…
In this paper, we develop Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations by applying conforming finite element methods in space and splitting methods in time. For the spatial discretisation, the criteria for choosing…
Numerical schemes that preserve the structure of the kinetic equations can provide stable simulation results over a long time. An electromagnetic particle-in-cell solver for the Vlasov-Maxwell equations that preserves at the discrete level…
This paper proposes a novel numerical integrator for modeling multispecies Coulomb collisions in kinetic plasmas. The proposed scheme provides an energy-, momentum-, and positivity-preserving particle discretization of the nonlinear Landau…
Recent results will be presented on steady states of magnetically confined plasmas obtained by conventional and Hamiltonian methods. The presentation consists of three parts. The first one concerns the derivation of a generalized…
An inclusive framework for joined Hamiltonian and dissipative dynamical systems, which preserve energy and produce entropy, is given. The dissipative dynamics of the framework is based on the metriplectic 4-bracket, a quantity like the…
Stability conditions of magnetized plasma flows are obtained by exploiting the Hamiltonian structure of the magnetohydrodynamics (MHD) equations and, in particular, by using three kinds of energy principles. First, the Lagrangian variable…
Port-Hamiltonian (PH) systems provide a framework for modeling, analysis and control of complex dynamical systems, where the complexity might result from multi-physical couplings, non-trivial domains and diverse nonlinearities. A major…
We present a geometric Particle-in-Cell (PIC) algorithm on two-dimensional (2D) unstructured meshes for studying electrostatic perturbations in magnetized plasmas. In this method, ions are treated as fully kinetic particles, and electrons…
We develop a comprehensive Hamiltonian formulation for plasma dynamics that unifies collisionless gyrokinetic and collisional processes. Our framework rigorously describes the evolution of free energy and entropy during the transition from…