Моделирование криволинейных поверхностей в задачах газовой динамики

Представлен способ задания начальных условий на границе тела произвольной формы на прямоугольной сетке. В качестве рассматриваемого тела выбран сферический объем сжатого газа, образовавшийся в результате взрыва над поверхностью Земли. Так как ячейки расчетной сетки прямоугольные, а контур криволинейный, то для задания условий на границе используются дробные ячейки. Давление и плотность внутри сферы известны и распределены равномерно по всему объему. Параметры на границе тела предлагается рассчитывать пропорционально объему, который занимает тело в каждой ячейке, через которую проходит контур. Такой объем может быть найден интегрированием по области, отсекаемой кривой от прямоугольной ячейки сетки. Тестирование алгоритма проводилось на численном решении задачи о разлете шара в чистом газе методом крупных частиц. Граница шара является контактным разрывом, поэтому для демонстрации работы метода приведены графики положения изолиний плотности в процессе расширении сферы. Результаты расчетов показали, что описанный механизм обеспечивает сохранение сферической границы в процессе счета: отклонение от значений, удовлетворяющих уравнению окружности, составило менее 1 %.