Перевод названия: Semi-Lagrangian method for two dimensional advection problem with discrete balance equation
Тип публикации: статья из журнала
Год издания: 2017
Ключевые слова: уравнение неразрывности, полуЛагранжевый метод, закон сохранения, балансовое соотношение, Semi-Lagrangian method, Balance equation, advection problem
Аннотация: Рассмотрен полулагранжевый метод численного решения уравнения неразрывности. Он основан на тождестве, связывающем двумерные интегралы решения на соседних слоях по времени и интегралы потоков на входе и выходе из области. Описана технология аппроксимации и вычисления интегралов, позволяющая найти численное решение с первым порядком Показать полностьюсходимости на гладких решениях. Обоснован дискретный закон сохранения для численного решения при переходе с одного слоя по времени на следующий слой без использования поправочных коэффициентов. 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.
Журнал: Вычислительные технологии
Выпуск журнала: Т. 22, № 5
Номера страниц: 27-38
ISSN журнала: 15607534
Место издания: Новосибирск
Издатель: Федеральное государственное бюджетное учреждение науки Институт вычислительных технологий Сибирского отделения Российской академии наук