Study of the accuracy and applicability of the difference scheme for solving the diffusion-convection problem at large grid P'eclet numbers

Автор: Sukhinov Aleksandr Ivanovich, Kuznetsova Inna Yuryevna, Chistyakov Aleksandr Evgenyevich, Protsenko Elena Anatolyevna, Belova Yuliya Valeriyevna

Журнал: Вычислительная механика сплошных сред @journal-icmm

Статья в выпуске: 4 т.13, 2020 года.

Бесплатный доступ

The work is devoted to the study of a difference scheme for solving the diffusion-convection problem at large grid Péclet numbers. The suspension transport problem numerical solving is carried out using the improved Upwind Leapfrog difference scheme. Its difference operator is a linear combination of the operators of Upwind and Standard Leapfrog difference schemes, while the modified scheme is obtained from schemes with optimal weighting coefficients. At certain values of the weighting coefficients, this combination leads to mutual compensation of approximation errors, and the resulting scheme gets better properties than the original schemes. In addition, it includes a cell filling function that allows simulating problems in areas with complex geometry. Computational experiments were carried out to solve the suspension transport problem, which arises, for example, during the propagation of suspended matter plumes in an aquatic environment and changes in the bottom topography due to the deposition of suspended soil particles into the sediment during soil unloading into a reservoir (dumping). The results of modeling the suspension transport problem at various values of the grid Péclet number are presented. The algorithm implementation was carried out using the software and hardware architecture of parallel computing: on a central processing unit (Central Processing Unit - CPU) and on a graphics accelerator (Graphics Processing Unit - GPU). The solution to the applied problem has shown its efficiency on the CPU with small computational grids and, if it is necessary to decrease the space steps, then the GPU solution is preferable. It was found that, when using the modified Upwind Leapfrog scheme, an increase in the speed of the water flow does not lead to a loss of solution accuracy due to dissipative sources and is accompanied by an insignificant increase in computational labor costs.

Еще

Mathematical model, suspension transport, diffusion-convection problem, numerical simulation, upwind leapfrog difference scheme, grid péclet number, parallel computing

Короткий адрес: https://sciup.org/143172508

IDR: 143172508   |   DOI: 10.7242/1999-6691/2020.13.4.34

Статья научная