Accuracy of the 3rd order Runge-Kutta method for a cubic dynamical systems

2018 ◽  
Vol 2018 (3) ◽  
pp. 29-36
Author(s):  
O.S. Akhmedov
Author(s):  
Yu Ying ◽  
Mikhail D. Malykh

We implement several explicit Runge-Kutta schemes that preserve quadratic invariants of autonomous dynamical systems in Sage. In this paper, we want to present our package ex.sage and the results of our numerical experiments. In the package, the functions rrk_solve, idt_solve and project_1 are constructed for the case when only one given quadratic invariant will be exactly preserved. The function phi_solve_1 allows us to preserve two specified quadratic invariants simultaneously. To solve the equations with respect to parameters determined by the conservation law we use the elimination technique based on Grbner basis implemented in Sage. An elliptic oscillator is used as a test example of the presented package. This dynamical system has two quadratic invariants. Numerical results of the comparing of standard explicit Runge-Kutta method RK(4,4) with rrk_solve are presented. In addition, for the functions rrk_solve and idt_solve, that preserve only one given invariant, we investigated the change of the second quadratic invariant of the elliptic oscillator. In conclusion, the drawbacks of using these schemes are discussed.


2017 ◽  
Vol 12 (5) ◽  
Author(s):  
Guang-Da Hu

In this paper, explicit Runge–Kutta methods are investigated for numerical solutions of nonlinear dynamical systems with conserved quantities. The concept, ε-preserving is introduced to describe the conserved quantities being approximately retained. Then, a modified version of explicit Runge–Kutta methods based on the optimization technique is presented. With respect to the computational effort, the modified Runge–Kutta method is superior to implicit numerical methods in the literature. The order of the modified Runge–Kutta method is the same as the standard Runge–Kutta method, but it is superior in preserving the conserved quantities to the standard one. Numerical experiments are provided to illustrate the effectiveness of the modified Runge–Kutta method.


Author(s):  
Jialin Tian ◽  
Jie Wang ◽  
Yi Zhou ◽  
Lin Yang ◽  
Changyue Fan ◽  
...  

Abstract Aiming at the current development of drilling technology and the deepening of oil and gas exploration, we focus on better studying the nonlinear dynamic characteristics of the drill string under complex working conditions and knowing the real movement of the drill string during drilling. This paper firstly combines the actual situation of the well to establish the dynamic model of the horizontal drill string, and analyzes the dynamic characteristics, giving the expression of the force of each part of the model. Secondly, it introduces the piecewise constant method (simply known as PT method), and gives the solution equation. Then according to the basic parameters, the axial vibration displacement and vibration velocity at the test points are solved by the PT method and the Runge–Kutta method, respectively, and the phase diagram, the Poincare map, and the spectrogram are obtained. The results obtained by the two methods are compared and analyzed. Finally, the relevant experimental tests are carried out. It shows that the results of the dynamic model of the horizontal drill string are basically consistent with the results obtained by the actual test, which verifies the validity of the dynamic model and the correctness of the calculated results. When solving the drill string nonlinear dynamics, the results of the PT method is closer to the theoretical solution than that of the Runge–Kutta method with the same order and time step. And the PT method is better than the Runge–Kutta method with the same order in smoothness and continuity in solving the drill string nonlinear dynamics.


Sign in / Sign up

Export Citation Format

Share Document