Symplectic schemes and symmetric schemes for nonlinear Schrödinger equation in the case of dark solitons motion

In this paper, symplectic schemes and symmetric schemes are presented to simulate Nonlinear Schrödinger Equation (NLSE) in case of dark soliton motion. Firstly, by Ablowitz–Ladik model (A–L model), the NLSE is discretized into a non-canonical Hamiltonian system. Then, different kinds of coordinate transformations can be used to standardize the non-canonical Hamiltonian system. Therefore, the symplectic schemes and symmetric schemes can be employed to simulate the solitons motion and test the preservation of the invariants of the A–L model and the conserved quantities approximations of the original NLSE. The numerical experiments show that symplectic schemes and symmetric schemes have similar simulation effect, and own significant superiority over non-symplectic and non-symmetric schemes in long-term tracking the motion of solitons, preserving the invariants and the approximations of conserved quantities. Moreover, it is obvious that coordinate transformations with more symmetry have a better simulation effect.

Wave Propagation in Rocks – Investigating the Effect of Rheology

Rocks exhibit beyond-Hookean, delayed and damped elastic, behaviour (creep, relaxation etc.). In many cases, the Poynting–Thomson–Zener (PTZ) rheological model proves to describe these phenomena successfully. A forecast of the PTZ model is that the dynamic elasticity coefficients are larger than the static (slow-limit) counterparts. This prediction has recently been confirmed on a large variety of rock types. Correspondingly, according to the model, the speed of wave propagation depends on frequency, the high-frequency limit being larger than the low-frequency limit. This frequency dependence can have a considerable influence on the evaluation of various wave-based measurement methods of rock mechanics. As experience shows, commercial finite element softwares are not able to properly describe wave propagation, even for the Hooke model and simple specimen geometries, the seminal numerical artefacts being instability, dissipation error and dispersion error, respectively. This has motivated research on developing reliable numerical methods, which amalgamate beneficial properties of symplectic schemes, their thermodynamically consistent generalization (including contact geometry), and spacetime aspects. The present work reports on new results obtained by such a numerical scheme, on wave propagation according to the PTZ model, in one space dimension. The simulation outcomes coincide nicely with the theoretically obtained phase velocity prediction.

The symbolic problems associated with Runge-Kutta methods and their solving in Sage

Runge-Kutta schemes play a very important role in solving ordinary differential equations numerically. At first we want to present the Sage routine for calculation of Butcher matrix, we call it an rk package. We tested our Sage routine in several numerical experiments with standard and symplectic schemes and verified our result by corporation with results of the calculations made by hand.Second, in Sage there are the excellent tools for investigation of algebraic sets, based on Gröbner basis technique. As we all known, the choice of parameters in Runge- Kutta scheme is free. By the help of these tools we study the algebraic properties of the manifolds in affine space, coordinates of whose are Butcher coefficients in Runge-Kutta scheme. Results are given both for explicit Runge-Kutta scheme and implicit Runge-Kutta scheme by using our rk package. Examples are carried out to justify our results. All calculation are executed in the computer algebra system Sage.

REBOUNDx: a library for adding conservative and dissipative forces to otherwise symplectic N-body integrations

ABSTRACT Symplectic methods, in particular the Wisdom–Holman map, have revolutionized our ability to model the long-term, conservative dynamics of planetary systems. However, many astrophysically important effects are dissipative. The consequences of incorporating such forces into otherwise symplectic schemes are not always clear. We show that moving to a general framework of non-commutative operators (dissipative or not) clarifies many of these questions, and that several important properties of symplectic schemes carry over to the general case. In particular, we show that explicit splitting schemes generically exploit symmetries in the applied external forces, which often strongly suppress integration errors. Furthermore, we demonstrate that so-called ‘symplectic correctors’ (which reduce energy errors by orders of magnitude at fixed computational cost) apply equally well to weakly dissipative systems and can thus be more generally thought of as ‘weak splitting correctors’. Finally, we show that previously advocated approaches of incorporating additional forces into symplectic methods work well for dissipative forces, but give qualitatively wrong answers for conservative but velocity-dependent forces like post-Newtonian corrections. We release REBOUNDx, an open-source C library for incorporating additional effects into REBOUNDN-body integrations, together with a convenient python wrapper. All effects are machine independent and we provide a binary format that interfaces with the SimulationArchive class in REBOUND to enable the sharing and reproducibility of results. Users can add effects from a list of pre-implemented astrophysical forces, or contribute new ones.

High-Order Symplectic Schemes for Stochastic Hamiltonian Systems

The construction of symplectic numerical schemes for stochastic Hamiltonian systems is studied. An approach based on generating functions method is proposed to generate the stochastic symplectic integration of any desired order. In general the proposed symplectic schemes are fully implicit, and they become computationally expensive for mean square orders greater than two. However, for stochastic Hamiltonian systems preserving Hamiltonian functions, the high-order symplectic methods have simpler forms than the explicit Taylor expansion schemes. A theoretical analysis of the convergence and numerical simulations are reported for several symplectic integrators. The numerical case studies confirm that the symplectic methods are efficient computational tools for long-term simulations.