M.N. Afanaseva*, E.B. Kuznetsov**
Moscow Aviation Institute (National Research University), Moscow, 125993 Russia
E-mail: *mary.mai.8@yandex.ru, **kuznetsov@mai.ru
Received April 24, 2019
DOI: 10.26907/2541-7746.2019.2.181-190
For citation: Afanaseva M.N., Kuznetsov E.B. The numerical solution of the nonlinear boundary value problem with singularity for the system of delay integrodifferential-algebraic equations. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2019, vol. 161, no. 2, pp. 181–190. doi: 10.26907/2541-7746.2019.2.181-190.
(In Russian)
Abstract
The numerical method for solving the nonlinear boundary value problem for a delay system of integrodifierential-algebraic equations was discussed. The occurrence of a delay argument in the system characterizes the behavior of the studied parameters not only at the current, but also at the previous period of time. Of particular interest are the problems characterized by the existence of singular limit points. It is very di–cult to solve these problems using the classical methods. A numerical solution of the boundary value problem was constructed by the shooting method. The values of the \shooting" parameter were found using a combination of the method of continuation with respect to the best parameter, the method of continuation with respect to the parameter in the Lahaye form, and the Newton method. At each iteration of the shooting method, the initial problem was solved. The computation of the initial problem influences the flnding of the required solution and the continuation of the iterative process of the shooting method. The initial problem was rearranged based on the best parameter - the length of the curve of the solution set, and flnite-difierence representation of derivatives. The resulting problem was solved by the Newton method. The values of the functions at the deviation point, where the values are not deflned by conditions of the problem, were calculated with the help of the Lagrange polynomial with three points. To flnd the value of the integral components of the problem, the trapezoid method was used. The results of the numerical study conflrm the efficiency of the proposed algorithm for solving the studied problem. The obtained numerical solution of the nonlinear boundary value problem with delay has the equation that loses its meaning in singular limit points. Thus, using of the continuation with respect to the best parameter while solving the problem allows to flnd all possible values of the parameter of the shooting method and to solve the problem.
Keywords: boundary value problem, numerical solution, difierential equations with delay, shooting method, method of continuation with respect to best parameter, singular limit points
Acknowledgments. The study was supported by the Russian Science Foundation (project no. 18-19-00474).
References
1. Na T.Y. Computational Methods in Engineering Boundary Value Problems. New York, Acad. Press, 1980. 309 p.
2. Ivanova E.P. Continuous dependence of solutions of boundary-value problems for differential-difference equations on shifts of the argument. CMFD, 2016, vol. 59, pp. 74–96. (In Russian)
3. Kamenskii G.A., Skubachevskii A.L. Lineinye kraevye zadachi dlya differentsial'noraznostnykh uravnenii [Linear Boundary Value Problems for Differential-Difference Equations]. Moscow, Izd. MAI, 1992. 190 p. (In Russian)
4. Kuznetsov E.B., Mikryukov V.N. Numerical integration of systems of delay differentialalgebraic equations. Comput. Math. Math. Phys., 2007, vol. 47, no. 1, pp. 80–92. doi: 10.1134/S0965542507010095.
5. Krasnikov S.D., Kuznetsov E.B. On the parametrization of numerical solutions to boundary value problems for nonlinear differential equations. Comput. Math. Math. Phys., 2005, vol. 45, no. 12, pp. 2066–2076.
6. Budkina E.M., Kuznetsov E.B. Modeling the technological process for manufacturing of aircraft structural components based on the best parametrization of the boundary value problem for nonlinear differential-algebraic equations. Vestn. MAI, 2016, vol. 23, no. 1, pp. 189–196. (In Russian)
7. Dmitriev S.S., Kuznetsov E.B. Numerical solution to systems of delay integrodifferential algebraic equations. Comput. Math. Math. Phys., 2008, vol. 48, no. 3, pp. 406–419. doi: 10.1134/S096554250803007X.
8. Shalashilin V.I., Kuznetsov E.B. Parametric Continuation and Optimal Parametrization in Applied Mathematics and Mechanics. Springer, 2003. viii, 228 p. doi: 10.1007/978-94- 017-2537-8.
9. Bakhvalov N.S., Zhidkov N.P., Kobel'kov G.M. Chislennye metody [Numerical Methods]. Moscow, Nauka, 1987. 600 p. (In Russian)
10. Samoilenko A.M., Ronto N.I. Chislenno-analiticheskie metody issledovaniya reshenii kraevykh zadach [Numerical and Analytic Methods in Studying the Solutions of Boundary Value Problems]. Kiev, Naukova Dumka, 1986. 222 p. (In Russian)
11. Davidenko D.F. On a new method of numerical solution of systems of nonlinear equations. Dokl. Akad. Nauk SSSR, 1953, vol. 88, no. 4, pp. 601–602. (In Russian)
12. Lahaye M.E. Une metode de resolution d'une categorie d'equations transcendentes. Compter Rendus hebdomataires des seances de L'Academie des sciences. 1934, vol. 198, no. 21, pp. 1840–1842. (In French)
13. Kuznetsov E.B. Optimal parametrization in numerical construction of curve. J. Franklin Inst., 2007, vol. 344, pp. 658–671.
The content is available under the license Creative Commons Attribution 4.0 License.