Related papers: High-order energy-conserving Line Integral Methods…
Path integral-based simulation methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical…
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…
In this work, we develop a class of up to third-order energy-stable schemes for the Cahn--Hilliard equation. Building on Lawson's integrating factor Runge--Kutta method, which is widely used for stiff semilinear equations, we discuss its…
Most classical mechanical systems are based on dynamical variables whose values are real numbers. Energy conservation is then guaranteed if the dynamical equations are phrased in terms of a Hamiltonian function, which then leads to…
Magneto-hydrodynamics is one of the foremost models in plasma physics with applications in inertial confinement fusion, astrophysics and elsewhere. Advanced numerical methods are needed to get an insight into the complex physical phenomena.…
In this paper, we develop a class of high-order conservative methods for simulating non-equilibrium radiation diffusion problems. Numerically, this system poses significant challenges due to strong nonlinearity within the stiff source terms…
In this paper, we apply the Boole discrete line integral to solve the Lorentz force system which is written as a non-canonical Hamiltonian system. The method is exactly energy-conserving for polynomial Hamiltonians of degree $\nu \leq 4$.…
We present a quantum algorithm to achieve higher-order transformations of Hamiltonian dynamics. Namely, the algorithm takes as input a finite number of queries to a black-box seed Hamiltonian dynamics to simulate a desired Hamiltonian. Our…
We present high-order variational Lagrangian finite element methods for compressible fluids using a discrete energetic variational approach. Our spatial discretization is mass/momentum/energy conserving and entropy stable. Fully implicit…
In ordinary turbulence research it has been a long standing tradition to solve the equations in spectral space giving the best possible accuracy. This is indeed a natural choice for incompressible problems with periodic boundaries, but it…
Exponential time differencing methods is a power tool for high-performance numerical simulation of computationally challenging problems in condensed matter physics, fluid dynamics, chemical and biological physics, where mathematical models…
It is shown analytically that the energy-conserving implicit nonsymplectic scheme of Bacchini, Ripperda, Chen and Sironi provides a first-order accuracy to numerical solutions of a six-dimensional conservative Hamiltonian system. Because of…
We propose an open-boundary molecular dynamics method in which an atomistic system is in contact with an infinite particle reservoir at constant temperature, volume and chemical potential. In practice, following the Hamiltonian adaptive…
This paper discusses a novel fully implicit formulation for a 1D electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov-Poisson…
In this paper we construct high order finite volume schemes on networks of hyperbolic conservation laws with coupling conditions involving ODEs. We consider two generalized Riemann solvers at the junction, one of Toro-Castro type and a…
We investigate discretization strategies for a recently introduced class of energy-based models. The model class encompasses classical port-Hamiltonian systems, generalized gradient flows, and certain systems with algebraic constraints. Our…
We provide a framework for high-order discretizations of nonlinear scalar convection-diffusion equations that satisfy a discrete maximum principle. The resulting schemes can have arbitrarily high order accuracy in time and space, and can be…
Energy-preserving numerical methods for solving the Hodge wave equation is developed in this paper. Based on the de Rham complex, the Hodge wave equation can be formulated as a first-order system and mixed finite element methods using…
In this paper, a novel high-order, mass and energy-conserving scheme is proposed for the regularized logarithmic Schr\"{o}dinger equation(RLogSE). Based on the idea of the supplementary variable method (SVM), we firstly reformulate the…
In this paper, we introduce and discuss an exactly energy-conserving Particle-in-Cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates enables the study of plasmas in complex-shaped domains…