Vestnik КRAUNC. Fiz.-Mat. nauki. 2024. vol. 49. no. 4. P. 135-156. ISSN 2079-6641

INFORMATION AND COMPUTING TECHNOLOGIES
https://doi.org/10.26117/2079-6641-2024-49-4-135-156
Research Article
in Russian
MSC 85-04; 49N45; 26A33

Contents of this issue

Read Russian Version

Some Aspects of the Implementation of the PRPHMM 1.0 Software Package for Refining the Parameters of Hereditary Mathematical Models of Radon Transfer in a Storage Chamber

D. A. Tverdyi¹^{\ast}, E. O. Makarov²

¹Institute of Cosmophysical Research and Radio Wave Propagation FEB RAS, 684034, Paratunka village, Mirnaya str., 7, Russia
²Kamchatka Branch of the Federal Research Center «Unified Geophysical Service of the Russian Academy of Sciences», 683023, Petropavlovsk-Kamchatsky, Piipa Boulevard st., 9, Russia

Abstract. Mathematical models of some dynamic processes can be significantly enhanced by using
derivatives and integrals of non-integer order in them, taking into account effects that cannot be described by ordinary derivatives. For example, by using fractional Gerasimov- Caputo derivatives of constant and variable order, it is possible to take into account the memory effect in the process model, and the order of the derivative will be related to the intensity of the process. In particular, the authors have previously developed an hereditary α-model of the volumetric activity of radon, where the parameter α is related to the permeability of the medium. However, the question arises about determination of optimal values of both α and other parameters of the model. To solve the problem, it is possible to solve the inverse problem, a common type of problem in many scientific fields, where it is necessary to determine the values of model parameters from observed data, but it is impossible to make direct measurements of these parameters. The need for such an approach often arises when working with geological data. The article describes the software implementation of the PRPHMM 1.0 software package which can clarifying optimal values of hereditary mathematical models based on the Gerasimov- Caputo derivative. The Levenberg-Marquardt unconditional Newtonian optimisation algorithm is adapted and implemented in MATLAB language. Subroutines for reading, processing and visualisation of experimental and model data are implemented. A test case solving the inverse problem for the hereditary α-model for the parameters α and λ0-air exchange coefficient on the basis of experimental radon monitoring data is presented. It is shown that PRPHMM 1.0 allows for the clarify of parameter values close to the optimum values for the hereditary mathematical models.

Key words: мathematical modeling, fractional derivatives, Gerasimov-Caputo, memory effect, nonlocality, nonlinear equations, inverse problems, unconditional optimization, Levenberg-Marquardt algorithm, MATLAB.

Received: 02.11.2024; Revised: 10.11.2024; Accepted: 24.11.2024; First online: 28.11.2024

For citation. Tverdyi D. A., Makarov E. O. Some aspects of the implementation of the PRPHMM 1.0 software package for refining the parameters of hereditary mathematical models of radon transfer in a storage chamber. Vestnik KRAUNC. Fiz.-mat. nauki. 2024, 49: 4, 135-156. EDN: FMWIIQ . https://doi.org/10.26117/2079-6641-2024-49-4-135-156.

Funding. The research was funded by a grant from the Russian Science Foundation, project number 22-11-00064, which can be found at https://rscf.ru/project/22-11-00064/.

Competing interests. The authors declare that there are no conflicts of interest regarding authorship and publication.

Contribution and Responsibility. All authors contributed to this article. Authors are solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by all authors.

^{\ast}Correspondence: E-mail: tverdyi@ikir.ru

The content is published under the terms of the Creative Commons Attribution 4.0 International License

© Tverdyi D. A., Makarov E. O., 2024

© Institute of Cosmophysical Research and Radio Wave Propagation, 2024 (original layout, design, compilation)

References

  1. Nakhushev A. M. Fractional calculus and its application. Moscow: Fizmatlit, 2003, 272 pp., (In Russian)
  2. Uchaikin V. V. Fractional Derivatives for Physicists and Engineers. Vol. I. Background and Theory. Berlin/Heidelberg, Springer, 2013, 373 pp. DOI: 10.1007/978-3-642-33911-0.
  3. Kilbas A. A., Srivastava H. M., Trujillo J. J. Theory and Applications of Fractional Differential Equations, 1st ed. Amsterdam, Elsevier, 2006, 523 pp.
  4. Samko S. G., Kilbas A. A., Marichev O. I. Integraly i proizvodnye drobnogo poryadka i nekotorye ih prilozheniya [Fractional integrals and derivatives and some of their applications]. Science and tech: Minsk, 1987, 688 pp. (In Russian)
  5. Tverdyi D. A., Parovik R. I. Application of the Fractional Riccati Equation for Mathematical Modeling of Dynamic Processes with Saturation and Memory Effect, Fractal and Fractional, 2022, vol. 6, no. 3, pp. 163. DOI: 10.3390/fractalfract6030163.
  6. Tverdyi D. A., Parovik R. I. Fractional Riccati equation to model the dynamics of COVID-19 coronavirus infection, Journal of Physics: Conference Series, 2021, vol. 2094, pp. 032042. DOI: 10.1088/1742-6596/2094/3/032042.
  7. Tverdyi D. A., Makarov E. O., Parovik R. I. Hereditary Mathematical Model of the Dynamics of Radon Accumulation in the Accumulation Chamber, Mathematics, 2023, vol. 11, no. 4, pp. 850. DOI: 10.3390/math11040850.
  8. Volterra V. Sur les ´equations int´egro-diff´erentielles et leurs applications, Acta Mathematica, 1912, vol. 35, no. 1, pp. 295–356. DOI: 10.1007/BF02418820.
  9. Parovik R. I. Tverdyi D. A. Some Aspects of Numerical Analysis for a Model Nonlinear Fractional Variable Order Equation, Mathematical and Computational Applications, 2021, vol. 26, no. 3, pp. 55. DOI: 10.3390/mca26030055.
  10. Tverdyi D. A., Parovik R. I. Investigation of Finite-Difference Schemes for the Numerical Solution of a Fractional Nonlinear Equation, Fractal and Fractional, 2022, vol. 6, no. 1, pp. 23. DOI: 10.3390/fractalfract6010023.
  11. Tverdyi D. A., Makarov E. O., Parovik R. I. Identification of parameters of the mathematical α-model of radon transport in the accumulation chamber based on data from the Karymshina site in Kamchatka, Bulletin KRASEC. Physical and Mathematical Sciences, 2024, vol. 48, no. 3, pp. 95–119. DOI: 10.26117/2079-6641-2024-48-3-95-119.(In Russian)
  12. Tarantola A. Inverse problem theory: methods for data fitting and model parameter estimation, Amsterdam and New York: Elsevier Science Pub. Co., 1987, 613 pp., isbn: 0444427651.
  13. Lailly P. The seismic inverse problem as a sequence of before stack migrations, Conference on Inverse Scattering, Theory and application, 1983, pp. 206–220.
  14. Firstov P.P., Makarov E. O. Dynamics of subsurface radon in Kamchatka and strong earthquakes. Petropavlovsk-Kamchatsky, Vitus Bering Kamchatka State University, 2018, 148 pp., isbn: 978-5-7968-0691-3 (In Russian)
  15. Firstov P.P., Rudakov V.P. Results from observations of subsurface radon in 1997-2000 at the Petropavlovsk-Kamchatskii geodynamic site. Journal of Volcanology and Seismology, 2003, no. 1, pp. 26–41 (In Russian)
  16. Utkin V. I., Yurkov A. K. Radon as a tracer of tectonic movements, Russian Geology and Geophysics, 2010, vol. 51, no. 2, pp. 220–227. DOI: 10.1016/j.rgg.2009.12.022
  17. Biryulin S. V., Kozlova I. A., Yurkov A. K. Investigation of informative value of volume radon activity in soil during both the stress build up and tectonic earthquakes in the South Kuril region, Bulletin of KRASEC. Earth Sciences, 2019, vol. 4, no. 44, pp. 73–83. DOI: 10.31431/1816-5524-2019-4-44-73-83.
  18. Gerasimov A. N. Generalization of linear deformation laws and their application to internal friction problems, Applied Mathematics and Mechanics, 1948, vol. 12, pp. 529–539.
  19. Caputo M. Linear models of dissipation whose Q is almost frequency independent – II, Geophysical Journal International, 1967, vol. 13, no. 5, pp. 529–539. DOI: 10.1111/j.1365-246X.1967.tb02303.x.
  20. Dennis J. E., Robert Jr., Schnabel B. Numerical methods for unconstrained optimization and nonlinear equations. Philadelphia, SIAM, 1996, 394 pp., isbn: 9781611971200
  21. Gill P. E., Murray W., Wright M. H. Practical Optimization. Philadelphia, SIAM, 2019, 421 pp.
  22. Levenberg K. A method for the solution of certain non-linear problems in least squares, Quarterly of applied mathematics, 1944, vol. 2, no. 2, pp. 164–168. DOI: 10.1090/qam/10666.
  23. Marquardt D. W. An algorithm for least-squares estimation of nonlinear parameters, Journal of the society for Industrial and Applied Mathematics, 1963, vol. 11, no. 2, pp. 431–441. DOI: 10.1137/0111030.
  24. Tverdyi D. A., Parovik R. I. The optimization problem for determining the functional dependence of the variable order of the fractional derivative of the Gerasimov-Caputo type, Bulletin KRASEC. Physical and Mathematical Sciences, 2024, vol. 47, no. 2, pp. 35–57. DOI: 10.26117/2079-6641-2024-47-2-35-57.(In Russian)
  25. Ford W. Numerical linear algebra with applications: Using MATLAB, 1st edition. Massachusetts, Academic Press, 2014, 628 pp., isbn: 978-0123944351. DOI: 10.1016/C2011-0-07533-6

Information about the authors

Tverdyi Dmitrii Alexsandrovich – PhD (Phys. & Math.), Researcher, Electromagnetic Radiation Laboratory, Institute of Cosmophysical Research and Radio Wave Propagation, FEB RAS, Paratunka village, Russia, ORCID 0000-0001-6983-5258.

Makarov Evgeny Olegovich – PhD (Phys. & Math.), Senior Researcher, Acoustic and Radon Monitoring Laboratory, Kamchatka Branch of the Federal Research Centre “Unified Geophysical Service of the Russian Academy of Sciences Petropavlovsk-Kamchatsky, Russia, ORCID 0000-0002-0462-3657.

Read articleDownload article