Численно моделировалась задача о дифракции ударной волны на свободной сфере в ударной трубе. Использовалась модель вязкого газа в предположении ламинарности течения.
В рамках модели конечного объема использовалась различная аппроксимация для вязких и невязких потоков. Для аппроксимации невязких потоков использовалась TVD – схема Хартена с незначительной модификацией [1]. Для аппроксимации операторов второго порядка в тензоре вязкости использовался конечно - объемный подход, изложенный первоначально в [2]. Для повышения порядка аппроксимации использовался оператор вычисления вязких потоков, аналогичный предложенному в [3]. Существенно упрощается в рамках конечного объема процедура задания необходимых краевых условий, сводящаяся для вязкого обтекания к заданию на теле давления и экстраполяции нормальной проекции тензора вязкости и градиента температуры. Сравнение результатов расчетов данными экспериментов [4,5] показывают применимость использования численного алгоритма для расчета исследуемого класса течений.
1. Мартюшов C.H.. Вычисл. технологии. 1997. Т.2, № 6, С. 53-60.
2. Vinokur M. . J. of Comp. Phys. 1989.V. 81 pp.1-52.
3. Волков К.Н. Выч. методы и программ.-2004. Т.6, № 1. с.47-64.
4. Tanno H., Itoh K., Saito T., Abe A., Tokayama K. Shock Waves n.13, pp. 191-200, 2003.
5. Igra O., Tokayama K., Proc. Ro. Soc. London a.- 1993-V.442.-p.231.
Файл тезисов: | MartyushViscousMIT2013.doc |