Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh order volume preserving integrators for three kinds of divergence-free vector fields via commutator
92
0
0.0
(
0
)
Authors
Huiyan Xue
Ask ChatGPT about the research
No Arabic abstract
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 prove that the commutators of elementary divergence-free vector fields (EDFVF) of those three kinds are still divergence-free vector fields of the same kind. Assuming then there is only diagonal part of divergence-free vector field of the monomial basis, for those three kinds of divergence-free vector fields, we construct high order volume-preserving inte- grators using the multi-commutators for EDFVFs. Moreover, we consider the ordering of the EDFVFs and their commutators to reduce the error of the schemes, showing by numerical tests that the strategy in [9] works very well.
rate research
Read More
A novel class of high-order linearly implicit energy-preserving exponential integrators are proposed for the nonlinear Schrodinger equation. We firstly done that the original equation is reformulated into a new form with a modified quadratic energy by the scalar auxiliary variable approach. The spatial derivatives of the system are then approximated with the standard Fourier pseudo-spectral method. Subsequently, we apply the extrapolation technique to the nonlinear term of the semi-discretized system and a linearized system is obtained. Based on the Lawson transformation, the linearized system is rewritten as an equivalent one and we further apply the symplectic Runge-Kutta method to the resulting system to gain a fully discrete scheme. We show that the proposed scheme can produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution and is extremely efficient in the sense that only linear equations with constant coefficients need to be solved at every time step. Numerical results are addressed to demonstrate the remarkable superiority of the proposed schemes in comparison with other high-order structure-preserving method.
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/2rfloor$.As an application, we derive high order energy-preserving methods for stochastic Poisson systems as well as further geometric numerical schemes for this wide class of Stratonovich SDEs.
The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and analyze a time integrator having the same property for all times. Furthermore, strong and weak convergence of the numerical scheme along with efficient multilevel Monte Carlo estimators are studied. Finally, extensive numerical experiments illustrate the performance of the proposed numerical scheme.
We perform a numerical analysis of a class of randomly perturbed {H}amiltonian systems and {P}oisson systems. For the considered additive noise perturbation of such systems, we show the long time behavior of the energy and quadratic Casimirs for the exact solution. We then propose and analyze a drift-preserving splitting scheme for such problems with the following properties: exact drift preservation of energy and quadratic Casimirs, mean-square order of convergence one, weak order of convergence two. These properties are illustrated with numerical experiments.
The goal of this paper is to develop energy-preserving variational integrators for time-dependent mechanical systems with forcing. We first present the Lagrange-dAlembert principle in the extended Lagrangian mechanics framework and derive the extended forced Euler-Lagrange equations in continuous-time. We then obtain the extended forced discrete Euler-Lagrange equations using the extended discrete mechanics framework and derive adaptive time step variational integrators for time-dependent Lagrangian systems with forcing. We consider two numerical examples to study the numerical performance of energy-preserving variational integrators. First, we consider the example of a nonlinear conservative system to illustrate the advantages of using adaptive time-stepping in variational integrators. We show a trade-off between energy-preserving performance and accurate discrete trajectories while choosing an initial time step. In addition, we demonstrate how the implicit equations become more ill-conditioned as the adaptive time step increases through a condition number analysis. As a second example, we numerically simulate a damped harmonic oscillator using the adaptive time step variational integrator framework. The adaptive time step increases monotonically for the dissipative system leading to unexpected energy behavior.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
