Научная статья на тему 'Pseudospectral Methods for Nonlinear Pendulum Equations'

Pseudospectral Methods for Nonlinear Pendulum Equations Текст научной статьи по специальности «Физика»

CC BY
6
2
Поделиться
Ключевые слова
chaotic pendulum / Chebyshev / pseudospectral methods / differentiation matrices / collocation method / nonlinear equations.

Аннотация научной статьи по физике, автор научной работы — Le Anh Nhat

The article searched on mathematics and numerical solutions for the nonlinear pendulum (Chaotic pendulum). The numerical solution that was used for our research suitably the pseudospectral methods. With these equations, we studied and calculated on the interval [−1, 1], with boundary conditions already known. We used the software Mathematica 10.4 to calculate the results of the problems.

Псевдоспектральные методы нелинейных уравнений маятника

В статье исследованы математические и численные решения для нелинейного маятника (хаотический маятник). Решение, которое использовалось для наших исследований, соответствовало псевдоспектральным методам. Вычисления проводили на интервале [–1,1] с уже известными граничными условиями. Для расчета использовалось программное обеспечение Mathematica 10.4.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Текст научной работы на тему «Pseudospectral Methods for Nonlinear Pendulum Equations»

Journal of Siberian Federal University. Mathematics & Physics 2019, 12(1), 79—84

YAK 517.9

Pseudospectral Methods for Nonlinear Pendulum Equations

Le Anh Nhat*

RUDN University 6 Miklukho-Maklaya str., Moscow, 117198

*

Russia

Tan Trao University Tuyen Quang 22227

Vietnam

Received 29.05.2018, received in revised form 10.08.2018, accepted 20.10.2018

The article searched on mathematics and numerical solutions for the nonlinear pendulum (Chaotic pendulum). The numerical solution that was used for our research suitably the pseudospectral methods. With these equations, we studied and calculated on the interval [—1,1], with boundary conditions already known. We used the software Mathematica 10.4 to calculate the results of the problems.

Keywords: chaotic pendulum, Chebyshev, pseudospectral methods, differentiation matrices, collocation

method, nonlinear equations.

DOI: 10.17516/1997-1397-2019-12-1-79-84.

The differential equation of a single pendulum is

here v is the time, d is value of the angular displacement, the number value v2 = g/l, the clarify the value g is the acceleration due to gravity and l is the length of the pendulum.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

The fact that the pendulum is placed in a friction environment, the drag force on the pendulum is calculated accordance with the Stokes law and is propertied to the instantaneous velocity of the pendulum because of the environment has value friction the equation of motion (1) becomes

and here by a is the coefficient of friction of the environment and show that is a > 0.

To resist the friction, we need to add some external force, to simplified we choose the external periodic force and so the equation of motion becomes the equation:

d2

—r e(t) = -v2 sin e(t),

(1)

d2 d

—2 e(t) = -v2 sin e(t) - a—0(t),

(2)

(3)

here 3 is amplitude and w is frequency. When a = 0 and then (3) becomes

—2 e(t) = -v2 sin e(t) + ß cos(wt),

(4)

* leanhnhat@mail.ru © Siberian Federal University. All rights reserved

all the above equations are nonlinear second order differential equations. The nonlinear equation (3) has no analytic solution. However, if 0 ^ 0 then sin(0) « 0, so the nonlinear differential equations above can be transformed into linear differential equations. This has got a lot of solutions, so we do not study.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

With these nonlinear equations, we will be study and compute by the pseudospectral method.

1. Pseudospectral method (PSM)

Let p(x) is a polynomial of degree n, and we know that it is values at the points p(x0), p(x1), ..., p(xn), then the first and second derivatives p(x) at the same points are expressed in matrix form:

ip'(xo)\ (p(xo)\ (p"(x0)\ (p(x0y

D

V (xn )/

D

V '(xn)/

(5)

p(xn)

\p(xn))

where D = {ditj} is the so-called differentiation matrix [4].

In case when the Chebyshev-Gauss-Lobatto points are chose as the collocation points yk = cos(kn/n) [5],

2n2 + 1

Di

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

6

(-1)i+j

2cj sin[n(i + j)/(2n)]sin[n(i - j)/(2n)]

cos(nj/n) 2sin(nj/n) 2n2 + 1

6

i = j =0 i = j

0 < i = j < n i=j=n

(6)

here ck

1

if k if k

0, n

1, 2,

, n — 1.

The application of differential algebra in ordinary differential equations can also extend to nonlinear differential equations, so we transformed the matrix D into matrices [3]:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

E(1) = {dij}, 1 < i, j < n - 1, eo = {di,o}, en = {di,n}, 0 < i <n

(7)

for a first-order differential element, the form u'(x ) = Ed) u(xi) + beo + aen here a and b are the two-point boundary-value on the range [-1; 1] of problem.

For a second-order differential element, we use D2 = {d^j} and define the matrices:

E(2) = d2)}, 1 < i, j < n - 1, e0 = K5L4 = K(2n}, 0<i<n

(8)

has the form u"(xj) = E(2)u(xj) + be02) + aen\

When we applied Pseudospectral method to solve nonlinear differential equations, nonlinear equations may or may not have a unique solution.

If we have identified a solution, we will proceed with an iterative procedure. Therefore, it is important to determine the iterative equation. The iterative procedure is simple, we assume u(0) = const, then find u(1), u(2), ..., stop it until the error e = |u(fc) — u(k-1)| < e0.

c

2

2. Results

Equation (1), we shall consider the two-point boundary-value of problem on the range [—1,1], 0( — 1) = a, 0(1) = b. We can find and transfer to form of method:

E(2)0(tj) = —v2 sin 0(tj) — be{2) — aeff, j = 1,n — 1. (9)

Equation (9) are iterative equation. We give the error e < 10~12. In the Fig. 1 and Tab. 1 are the result are calculated based on the program by pseudospectral method, and then the solid line show the result calculated by Mathematica 10.4.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1.0 -0.5 0.5 1.0

Fig. 1. Graphic of equation with v = 1.7, n = 64 in the case a = b = 0.9

Table 1. Competition the numerical results and error with Mathematica's calculations of the equation (1) with v = 1.7, n = 64 in the case a = b = 0.9

k y(k) PSM Mathematica 10.4 Error

1 0.998795 0.903167 0.903167 1.53307x10" -7

10 0.881921 1.19371 1.19371 1.33523x10" -7

20 0.55557 1.80554 1.80554 7.05996x10" 8

30 0.0980171 2.17294 2.17294 4.01061x10" 8

40 -0.382683 2.00793 2.00793 2.20135x10" 8

50 -0.77301 1.43174 1.43174 1.19427x10" 8

60 -0.980785 0.950125 0.950125 1.65847x10" -9

Remasks: by this case, when 0 < v < 1.8, the results are correct and convergence. When v < 1.8, means that the length of the string is at least g/v2. When v > 1.8, equation (1) non-compliance boundary conditions 0(—1) = a, 0(1) = b.

Equation (2), we will be consider on the range [—1,1] and the boundary conditions 0( —1) = a, 0(1) = b. We can find iterative equation and transfer to form of method:

(E(2) + a E(1))0 = f — b(e0 + aeo) — a(e^ + aen), j = 1,n — 1 (10)

here 0 and f denote the vectors with the elements {tj} and {—v2 sin 0(tj)}.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

In this problem, we give the error e < 10~8. In the Fig. 2, points of value are the calculated by result of the program by pseudospectral method, and then the solid line show the result

1.0 -0.5 0.5 1.0

Fig. 2. Graphic of equation (2) with a = 10, v = 1.7, n = 96 in the case a = b = —0.3

Table 2. Competition the numerical results and error with Mathematica's calculations of the equation (2) with a =10, v = 1.7, n = 96 in the case a = b = —0.3

k y(k) PSM Mathematica 10.4 Error

5 0.986643 -0.301648 -0.301648 6.69452x10- -10

20 0.793353 -0.326501 -0.326501 3.08531x10- -9

35 0.412707 -0.381289 -0.381289 5.35524x10- -y

50 -0.0654031 -0.462363 -0.462363 6.22867x10- -9

65 -0.528068 -0.551825 -0.551825 6.90365x10- -9

80 -0.866025 -0.533131 -0.533131 1.19349x10- -7

95 -0.999465 -0.301757 -0.301757 7.68573x10- -7

calculated by Mathematica. In the Tab. 2 are the competition the numerical results and error with Mathematica's calculations with a =10, v = 1.7, n = 96 in the case a = b = —0.3.

Remasks: The fixed point, O(t) = dO/dt = 0 is linearly stable, i.e., small perturbations from this point will decay in time. The fixed point O(t) = n, dO/dt = 0 is linearly unstable, which means that small perturbations from this point will grow exponentially in time [2]. When 0 < v < 2.7, the results are correct and convergence; when v ^ 2.7 method for big error and not convergence. When v fixed, a increase, then graph are left-leaning.

Equation (3), we will be consider on the range [—1,1] and in the case of boundary conditions O( — 1) = a, O(1) = b. Equation (3) is non-Hamiltonian and it does not have an analytical solution. This is a nonlinear equation and has three dynamic variables. We found iterative equation and transferred to form of the method:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+ a = f — b(e0 + aeo) — a(e^ + aen), j = 1,n — 1 (11)

here O and f denote the vectors with the elements {tj} and {ft coswtj — v2 sinO(tj)}.

With v, a, ft,w > 0, we give the error e ^ 10~12. In the Fig. 3, points of value are calculated the result of the program by the PSM, and then the solid line show the result calculated by Mathematica. Numerical results are given in Tab. 3.

Remasks: When 0 < v< 1.7 the results are correct and convergence; when ft cos(wt) descended then O(t) ascended; when a ascended then O(t) descended; the complex motion one would expect when the three forces are comparable; from here we can orient the pendulum control.

0.10 -

C.OB

\ o.oe

1 I 1 I 0.02 .....

-1.0 -0.5 0.5 1.0

Fig. 3. Graphic of equation (3) with a = 0.5, v = 1.1, w = 2n, ¡3 = 1.3, n = 128 in the case a = b = 0

Table 3. Competition the numerical results and error with Mathematica's calculations of the equation (3) with a = 0.5, v = 1.1, w = 2n, 3 =1.3, n = 128 in the case a = b = 0

k y(k) PSM Mathematica 10.4 Error

1 0.999699 0.0000121827 0.0000121743 8.4617x10- 10

25 0.817585 0.027557 0.027557 1.50768x10 -8

50 0.33689 0.0854559 0.0854558 2.98413x10 -8

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

75 -0.266713 0.0722918 0.0722917 3.26543x10 -8

100 -0.77301 0.0487747 0.0487747 4.7381 x10- 9

125 -0.99729 0.000283445 0.000283453 7.51629x10 -9

Equation (4), we will be consider on the range [—1,1] and in the case of boundary conditions 0( — 1) = a, 0(1) = b. We found iterative equation and transferred to form of method:

E(2)0 = f — be{2) — ae(2\ j = 1,n — 1 (12)

here 0 and f denote the vectors with the elements {tj} and {3coswtj — v2 sin0(tj)} with v,3,w > 0, we give the error e < 10~12. In the Fig. 4, the value are the calculated base on the program by the PSM, and then the solid line show the result calculated by Mathematica. The calculation results are shown in Tab. 4.

Remasks: When 0 < v < 1.9 the results is show correct and convergence; when 3 cos(wt) ascended then 0(t) ascended; when w ascended then 0(t) descended.

Fig. 4. Graphic of equation (4) with v = 1.0, ß = 20, w = 20, n = 128 in the case a = b = 0

Table 4. Competition the numerical results and error with Mathematica's calculations of the equation (4) with v = 1.0, ¡3 = 20, w = 20, n = 128 in the case a = b = 0

k y(k) PSM Mathematica 10.4 Error

1 0.999699 -0.000265697 -0.000265549 1.47924x10" ■7

30 0.740951 0.0595138 0.059514 1.47715x10" ■7

60 0.0980171 0.0566888 0.0566889 1.26391x10" ■7

90 -0.595699 -0.0085085 -0.00850845 5.6271 x10"8

127 -0.999699 -0.000265697 -0.000265696 9.03764x10" 10

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3. Summary

We have monitored four equations of the nonlinear pendulum with the range [—1,1] and the know boundary conditions by the pseudospectral methods. The method has issued good results. And then we can survey the variation of the function through our parameters and design of force acting on the pendulum to control the operation the pendulum.

The publication was prepared with the support of the "RUDN University Program 5-100".

References

[1] A.Belendez, E.Arribas, Approximate solutions for the nonlinear pendulum equation using a rational harmonic representation, Article in Computers & Mathematics with Applications, 9(2012), no. 64, 1602-1611.

[2] G.Moshe, The Chaotic Pendulum, World Scientific Publishing Co. Pte. Ltd, 2010.

[3] G.L.Baker, J.A.Blackburn, The Pendulum — a case study in physics, New York, Oxford University Press Inc., 2005.

[4] J.C.Mason, D.C.Handscomb, Chebyshev Polynomials, CRC Press LLC, 2003.

[5] W.S.Don, A.Solomonoff, Accuracy and speed in computing the Chebyshev collocation de-vivative, SIAM Juarnal of Scientific Computing, 16(1991), no. 6, 1253-1268.

Псевдоспектральные методы нелинейных уравнений маятника

Ле Ань Ньат

Российский университет дружбы народов Миклухо-Маклая, 6, Москва, 117198

Россия

В статье исследованы математические и численные решения для нелинейного маятника (хаотический маятник). Решение, которое использовалось для наших исследований, соответствовало псевдоспектральным методам. Вычисления проводили на интервале [-1, 1] с уже известными граничными условиями. Для расчета использовалось программное обеспечение Mathematica 10.4.

Ключевые слова: хаотический маятник, Чебышев, псевдоспектральные методы, дифференцирующие матрицы, метод коллокации, нелинейные уравнения.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.