[Текст] : статья / Александр Владимирович Вяткин, Елена Владимировна Кучунова, Владимир Викторович Шайдуров // Вычислительные технологии. - 2017. -
Т. 22,
№ 5. - С. 27-38
. - ISSN 1560-7534
Перевод заглавия: Semi-Lagrangian method for two dimensional advection problem with discrete balance equation
Аннотация: Рассмотрен
полулагранжевый метод численного решения уравнения неразрывности. Он основан на тождестве, связывающем двумерные интегралы решения на соседних слоях по времени и интегралы потоков на входе и выходе из области. Описана технология аппроксимации и вычисления интегралов, позволяющая найти численное решение с первым порядком сходимости на гладких решениях. Обоснован дискретный закон сохранения для численного решения при переходе с одного слоя по времени на следующий слой без использования поправочных коэффициентов.
In this paper, the semi-Lagrangian method (the generalized method of characteristics) is considered for the numerical solution of the two-dimensional (in space) continuity equation for the gas or fluid motion with given velocities. The geometric pattern of the method is represented by a fluid flow tube with an arbitrary two-dimensional rectangular cell at the upper time level constructed by passing flow lines (characteristics of the equation) backward in time to the intersection with the previous time level or the inflow surface. For the obtained three-dimensional (in the time-space coordinates) body, the Gauss -Ostrogradsky theorem is applied. Due to the absence of the gas flow across the side surfaces defined by trajectories, the relation between two-dimensional integrals of a solution at neighbouring time levels (in some cases on the inflow surface) is obtained. The paper primarily deals with the detailed approximation of these integrals resulting in explicit grid equations where the values of an approximate solution at the upper time level are expressed in terms of the values at the lower level and on the inflow surface. A fluid flow tube for gas issuing through a rectangular cell on the outflow surface is considered in a similar way. After the corresponding approximation, the quantity of issuing gas is expressed in terms of its integrals at the previous time level. The proposed methods of the approximation of integrals provide the fulfillment of the mass conservation law at the discrete level for the motion of gas between two time levels taking into account gas inflow and outflow. Numerical results confirming the fulfillment of the discrete conservation law are presented. Due to the dynamical adaptation of the grid pattern, the Courant -Friedrichs -Lewy condition for the time step size is removed.
РИНЦ Держатели документа: Институт вычислительного моделирования СО РАН
Сибирский федеральный университет
Доп.точки доступа: Вяткин, Александр Владимирович; Vyatkin Alexander Vladimirovich; Кучунова, Елена Владимировна; Kuchunova Elena Vladimirovna; Шайдуров, Владимир Викторович; Shaydurov Vladimir Viktorovich
