相关论文: A sixth order averaged vector field method
New families of composition methods with processing of order 4 and 6 are presented and analyzed. They are specifically designed to be used for the numerical integration of differential equations whose vector field is separated into three or…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
We discuss the efficient implementation of Hamiltonian BVMs (HBVMs), a recently introduced class of energy preserving methods for canonical Hamiltonian systems, via their blended formulation. We also discuss the case of separable problems,…
Recently, a whole class of evergy-preserving integrators has been derived for the numerical solution of Hamiltonian problems. In the mainstream of this research, we have defined a new family of symplectic integrators depending on a real…
As is well known, energy is generally deemed as one of the most important physical invariants in many conservative problems and hence it is of remarkable interest to consider numerical methods which are able to preserve it. In this paper,…
We consider structure-preserving methods for conservative systems, which rigorously replicate the conservation property yielding better numerical solutions. There, corresponding to the skew-symmetry of the differential operator, that of…
In this paper we are concerned with the analysis of a class of geometric integrators, at first devised in [14, 18], which can be regarded as an energy-conserving variant of Gauss collocation methods. With these latter they share the…
In this paper, we focus on the construction of high order volume preserving in- tegrators for divergence-free vector fields: the monomial basis, the exponential basis and tensor product of the monomial and the exponential basis. We first…
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…
The active flux (AF) method is a compact high-order finite volume method that simultaneously evolves cell averages and point values at cell interfaces. Within the method of lines framework, the existing Jacobian splitting-based point value…
We introduce an algorithmic approach based on generating tree method for enumerating the inversion sequences with various pattern-avoidance restrictions. For a given set of patterns, we propose an algorithm that outputs either an accurate…
In this article, we construct a representation formula for stochastic B-series evaluated in a B-series. This formula is used to give for the first time the order conditions of implicit Taylor methods in terms of rooted trees. Finally, as an…
In this paper, a class of arbitrarily high-order linear momentum-preserving and energy-preserving schemes are proposed, respectively, for solving the regularized long-wave equation. For the momentum-preserving scheme, the key idea is based…
The tree code for the approximate evaluation of gravitational forces is extended and substantially accelerated by including mutual cell-cell interactions. These are computed by a Taylor series in Cartesian coordinates and in a completely…
We show that applying any deterministic B-series method of order $p_d$ with a random step size to single integrand SDEs gives a numerical method converging in the mean-square and weak sense with order $\lfloor p_d/2\rfloor$.As an…
A rank-adaptive integrator for the approximate solution of high-order tensor differential equations by tree tensor networks is proposed and analyzed. In a recursion from the leaves to the root, the integrator updates bases and then evolves…
A novel hybrid spectral difference/embedded finite volume method is introduced in order to apply a discontinuous high-order method for large scale engineering applications involving discontinuities in the flows with complex geometries. In…
The effective Hamiltonian of $m\alpha^6$ and $m\alpha^6(m/M)$ order corrections for hydrogen molecular ions has been derived in [ Z.-X. Zhong, \emph{et al.}, Phys. Rev. A {\bf98}, 032502(2018).], in this work we express the energy…
In this paper, a centred universal high-order finite volume method for solving hyperbolic balance laws is presented. The scheme belongs to the family of ADER methods where the Generalized Riemann Problems (GRP) is a building block. The…
In this paper, high order well-balanced finite difference weighted essentially non-oscillatory methods to solve general systems of balance laws are presented. Two different families are introduced: while the methods in the first one…