相关论文: An energetically balanced, quasi-Newton integrator…
A new horizontally explicit/vertically implicit (HEVI) time splitting scheme for atmospheric modelling is introduced, for which the horizontal divergence terms are applied within the implicit vertical substep. The new HEVI scheme is…
Atmospheric systems incorporating thermal dynamics must be stable with respect to both energy and entropy. While energy conservation can be enforced via the preservation of the skew-symmetric structure of the Hamiltonian form of the…
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…
We introduce energy-preserving integrators for nonholonomic mechanical systems. We will see that the nonholonomic dynamics is completely determined by a triple $({\mathcal D}^*, \Pi, \mathcal{H})$, where ${\mathcal D}^*$ is the dual of the…
A very simple explicit integrator for the rotational motion of rigid linear molecules is presented which can preserve the rigidity of the molecules without requiring any constraint force. The integrator is time-reversible and symplectic,…
A model of the three-dimensional rotating compressible Euler equations on the cubed sphere is presented. The model uses a mixed mimetic spectral element discretization which allows for the exact exchanges of kinetic, internal and potential…
We introduce a second-order numerical scheme for compressible atmospheric motions at small to planetary scales. The collocated finite volume method treats the advection of mass, momentum, and mass-weighted potential temperature in…
Dynamic simulation of elastic bodies is a longstanding task in engineering and computer graphics. In graphics, numerical integrators like implicit Euler and BDF2 are preferred due to their stability at large time steps, but they tend to…
Global non-hydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the…
Energy-correction method is proposed as an addition to mainstream integrators for equations of motion of systems of classical spins. This solves the problem of non-conservation of energy in long computations and makes mainstream integrators…
We present a collection of well-conditioned integral equation methods for the solution of electrostatic, acoustic or electromagnetic scattering problems involving anisotropic, inhomogeneous media. In the electromagnetic case, our approach…
This work proposes a suite of numerical techniques to facilitate the design of structure-preserving integrators for nonlinear dynamics. The celebrated LaBudde-Greenspan integrator and various energy-momentum schemes adopt a difference…
The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. Here, a novel, implicit, particle-in-cell based…
We develop a variational integrator for the shallow-water equations on a rotating sphere. The variational integrator is built around a discretization of the continuous Euler-Poincar\'{e} reduction framework for Eulerian hydrodynamics. We…
Coherent or exact equations of motion for a post-Newtonian Lagrangian formalism are the Euler-Lagrange equations without any terms truncated. They naturally conserve energy {and} angular momentum. Doubling the phase-space variables of…
Non-hydrostatic atmospheric models often use semi-implicit temporal discretisations in order to negate the time step limitation of explicitly resolving the fast acoustic and gravity waves. Solving the resulting system to machine precision…
By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for…
We introduce novel high order well-balanced finite volume methods for the full compressible Euler system with gravity source term. They require no a priori knowledge of the hydrostatic solution which is to be well-balanced and are not…
In this paper, we formulate and analyse exponential integrations when applied to nonlinear Schr\"{o}dinger equations in a normal or highly oscillatory regime. A kind of exponential integrators with energy preservation, optimal convergence…
Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved…