Mathematical modeling of eutrophication processes in Azov Sea on supercomputers

Introduction. The paper is devoted to the research of effective numerical methods for solving eu-trophication problem of shallow waters taking into account water environment, spatially-nonuniform distribution of temperature and salinity, microturbulent diffusion, gravitational sedi-mentation, and spreading of biogenic pollution, oxygen, phyto- and zooplankton, etc. The simula-tion objects are shallow waters – the Azov Sea and Taganrog Bay. Materials and Methods. The mathematical model of eutrophication of shallow waters was devel-oped. Parallel implementation was performed for computationally laborious convection-diffusion problems, taking into account the architecture and parameters of supercomputers based on the de-composition methods of grid domains. We determined that the maximum acceleration was 43 times on 128 computational nodes. We developed two algorithms including the algorithm based on the k-means method for data distribution between processors in parallel implementation. Due to the using these algorithms, the efficiency of algorithm for solving the problem is increased on 15% compared to the algorithm of the standard partition of computational domain. Results. New mathematical models and software complex were developed for mathematical model-ing of eutrophication processes in shallow waters. The concentrations of pollutants and plankton calculated for different wind situations were taken into consideration, if the relative error did not exceed 30%. Due to expedition researches the primary verification of the model of ecosystem of the Azov Sea was performed. The problem of modeling and forecasting the state of water ecosystems of the Azov Sea in conditions of anthropogenic influence and comprehensive research of the unique water ob-ject was implemented. Because of the water object is shallow, it’s more affected by anthropogenic influence. The software complex, combining mathematical models and databases, was designed. Using this complex we researched conditions which are contributed by the eutrophication processes in shallow waters. Discussion and Conclusions. Due to the solving the water ecology problem we can forecast differ-ent scenario of changing the water quality in shallow waters, and to investigate the mechanisms of formation of zones with low oxygen content.

1. Samarsky, A.A., Nikolaev, E.S. Methods of solving grid equations. Moscow, Science, 1978, pp. 588.

2. Sukhinov, A.I., Chistyakov, A.E. Adaptive modified alternating triangular iterative method for solving grid equations with non-selfadjoint operator. Mathematical modeling, 2012, vol. 24, no. 1, pp. 3–20.

3. State Research Center "Planeta", available at: http://planet.iitp.ru/english/index_eng.htm.

4. Samarsky, A.A. Theory of difference schemes. Moscow, Nauka, 1989.

5. Konovalov, A.N. The method of steepest descent with adaptive alternately-triangular pre-amplification. Differential equations, 2004, vol. 40, no. 7, pp. 953.

6. Konovalov, A.N. The theory of alternating-triangular iterative method. Siberian mathematical journal, 2002, vol. 43, no. 3, pp. 552.

7. Sukhinov, A.I., Chistyakov, A.E., Shishenya, A.V. Error estimate of the solution of the diffusion equation on the basis of the schemes with weights. Mathematical modeling, 2013, vol. 25, no. 11, pp. 53–64.

8. Chetverushkin, B., Gasilov, V., Iakobovski, M., Polyakov, S., Kartasheva, E., Boldarev, A., Abalakin, I., Minkin, A. Unstructured mesh processing in parallel CFD project GIMM. Parallel Computational Fluid Dynamics, Amsterdam, Elsevier, 2005, pp. 501–508.

9. Petrov, I.B., Favorsky, A.V., Sannikov, A.V., Kvasov, I.E. Grid-characteristic method using high order interpolation on tetrahedral hierarchical meshes with a multiple time step. Mathematical modeling, 2013, vol. 25, no. 2, pp. 42–52.

10. Sukhinov, A.I., Chistyakov, A.E., Semenyakina, A.A., Nikitina, A.V. Parallel realization of the tasks of the transport of substances and recovery of the bottom surface on the basis of schemes of high order of accuracy. Computational methods and programming: new computing technology, 2015, vol. 16, no. 2, pp. 256–267.

11. Chistyakov, A.E., Hachunts, D.S., Nikitina, A.V., Protsenko, E.A., Kuznetsova, I.Y. Parallel Library of iterative methods of the SLAE solvers for problem of convection -diffusion-based decomposition in one spatial direction. Modern problems of science and education, 2015, no. 1–1, pp. 1786.

12. Sukhinov, A.I., Nikitina, A.V., Semenyakina, A.A., Protsenko, E.A. Complex programs and algorithms to calculate sediment transport and multi-component suspensions on a multiprocessor computer system. Engineering journal of Don. 2015, vol. 38, no. 4 (38), pp. 52.

13. Nikitina, A.V., Abramenko, Y.A., Chistyakov, A.E. Mathematical modeling of oil spill in shallow waters. Informatics, computer science and engineering education, 2015, № 3 (23), pp. 49–55.

14. Chistyakov, A.E., Nikitina, A.V., Ougolnitsky, G.A., Puchkin, V.M., Semenov, I.S., Sukhinov, A.I., Usov, A.B. A differential game model of preventing fish kills in shallow waterbodies. Game Theory and Applications, 2015, vol. 17, pp. 37–48.

15. Sukhinov, A.I., Nikitina, A.V., Semenyakina, A.A., Chistyakov, A.E., A set of models, explicit regularized schemes of high order of accuracy and programs for predictive modeling of consequences of emergency oil spill. Parallel computational technologies (PCT ' 2016), proceedings of the international scientific conference, 2016, pp. 308–319.

16. Nikitina, A.V., Semenyakina, A.A., Chistyakov, A.E. Parallel implementation of the tasks of diffusion-convection-based schemes of high order of accuracy. Vestnik of computer and information technology, 2016, no. 7 (145), pp. 3–8.

17. Sukhinov, A.I., Chistyakov, A.E., Semenyakina, A.A., Nikitina, A.V. Numerical modeling of an ecological condition of the Sea of Azov with application of schemes of the raised accuracy order on the multiprocessor computing system. Computer researches and modeling, 2016, vol. 8, no. 1, pp. 151–168.

18. Sukhinov, A.I., Nikitina, A.V., Semenyakina, A.A., Chistyakov, A.E. Complex of models, explicit regularized schemes of high-order of accuracy and applications for predictive modeling of aftermath of emergency oil spill. CEUR Workshop Proceedings. 10th Annual International Scientific Conference on Parallel Computing Technologies, PCT 2016, Arkhangelsk, Russian Federation, 2016, vol. 1576, pp. 308–319.

19. Sukhinov, A.I., Nikitina, A.V., Chistyakov, A.E., Semenov, I.S., Semenyakina, A.A., Khachunts, D.S. Mathematical modeling of eutrophication processes in shallow waters on multiprocessor computer system. CEUR Workshop Proceedings. 10th Annual International Scientific Conference on Parallel Computing Technologies, PCT 2016, Arkhangelsk; Russian Federation, 2016, vol. 1576, pp. 320-333.

20. Nikitina, A.V., Sukhinov, A.I., Ougolnitsky, G.A., Usov, A.B., Chistyakov, A.E., Puchkin, M.V., Semenov, I.S. Optimal management of sustainable development at the biological rehabilitation of the Azov Sea. Mathematical modeling, 2016, vol. 28, no. 7, pp. 96–106.

21. Nikitina, A.V., Sukhinov, A.I., Ougolnitsky, G.A., Usov, A.B., Chistyakov, A.E., Puchkin, M.V., Semenov, I.S. Optimal control of sustainable development in the biological rehabilitation of the Azov Sea. Mathematical Models and Computer Simulations, 2017, no. 9 (1), pp. 101–107.