Related papers: General local energy-preserving integrators for so…
In this paper we design discrete port-Hamiltonian systems systematically in two different ways, by applying discrete gradient methods and splitting methods respectively. The discrete port-Hamiltonian systems we get satisfy a discrete notion…
We present a new parallel numerical method for solving the non-stationary Schr\"odinger equation with linear nonlocal condition and time-dependent potential which does not commute with the stationary part of the Hamiltonian. The given…
In this paper, a time-domain discontinuous Galerkin (TDdG) finite element method for the full system of Maxwell's equations in optics and photonics is investigated, including a complete proof of a semi-discrete error estimate. The new…
A new approach for the construction of high order A-stable explicit integrators for ordinary differential equations (ODEs) is theoretically studied. Basically, the integrators are obtained by splitting, at each time step, the solution of…
We propose a high order adaptive-rank implicit integrators for stiff time-dependent PDEs, leveraging extended Krylov subspaces to efficiently and adaptively populate low-rank solution bases. This allows for the accurate representation of…
This work presents a novel formulation and numerical strategy for the simulation of geometrically nonlinear structures. First, a non-canonical Hamiltonian (Poisson) formulation is introduced by including the dynamics of the stress tensor.…
We present two novel classes of fully discrete energy-preserving algorithms for the sine-Gordon equation subject to Neumann boundary conditions. The cosine pseudo-spectral method is first used to develop structure-preserving spatial…
We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit, K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original…
We design a consistent Galerkin scheme for the approximation of the vectorial modified Korteweg-de Vries equation. We demonstrate that the scheme conserves energy up to machine precision. In this sense the method is consistent with the…
A family of arbitrarily high-order fully discrete space-time finite element methods are proposed for the nonlinear Schr\"odinger equation based on the scalar auxiliary variable formulation, which consists of a Gauss collocation temporal…
Structure-preserving particle methods have recently been proposed for a class of nonlinear continuity equations, including aggregation-diffusion equation in [J. Carrillo, K. Craig, F. Patacchini, Calc. Var., 58 (2019), pp. 53] and the…
Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian…
Numerical schemes for the solution of the Euler equations have recently been developed, which involve the discretisation of the internal energy equation, with corrective terms to ensure the correct capture of shocks, and, more generally,…
We investigate a two-state conformational conversion system and introduce a novel structure-preserving numerical scheme that couples a local discontinuous Galerkin space discretization with the backward Euler time-integration method. The…
We study the recently-proposed hyperbolic approximation of the Korteweg-de Vries equation (KdV). We show that this approximation, which we call KdVH, possesses a rich variety of solutions, including solitary wave solutions that approximate…
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are…
This paper presents high-order Runge-Kutta (RK) discontinuous Galerkin methods for the Euler-Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to…
In this work, we introduce a structure-preserving local discontinuous Galerkin (LDG) method \cite{cockburn1998local} for solving the non-local non-linear Fokker-Planck-Landau (FPL) equations. We rephrase the structure-preserving strategy of…
Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general…
We study relaxation-based approaches for conserving mass and energy in the numerical solution of Schr\"odinger-Poisson (SP) type systems. Relaxation-based methods offer a general approach that can be applied as post-time step processing to…