Была проведена доработка методики для расчета нестационарного переноса в многомерных областях с заданными в начальный момент разрывами или тонкими погранслоями ([1]) для применения к расчетам течений вязкой несжимаемой вязкости. По полученной схеме было проведено моделирование различных течений вязкой жидкости при числах Рейнольдса, больших подсеточного: начальные участки истечения струй и развития вихревых структур с наличием внутренних погранслоев. Показано, что обычные сеточные методы с расщеплением по пространственным направлениям дают при этом существенные ошибки в решениях за счет повышенной схемной вязкости, превышающей моделируемую как по направлению переноса, так и в боковом направлении, а также не обладают инвариантностью относительно поворота регулярной сетки из-за зависимости схемной вязкости от направления скорости переноса к узловым линиям сетки.
В полной мере эти недостатки могут быть устранены схемой средней точности с алгоритмом уменьшения схемной вязкости. Для одного течения с числом Re = 1000 показано хорошее совпадение расчета по обычному методу расщепления при числе узлов 1280 в каждом пространственном направлении, когда схемная вязкость мала по сравнению с моделируемой, с расчетом по алгоритму уменьшения схемной вязкости ([2]) при числе узлов 40 и 80. Экономичность по машинным ресурсам у последнего расчета по сравнению с первым составляет более трех порядков. Во многих же публикациях расчеты проводятся с числом узлов 30 или 60 в каждом пространственном направлении, и для задачи с разрывами первого рода свойственно ожидать больших искажений в численном решениии.
1. Паничкин А.В. Численный расчет задач нестационарного конвективно- диффузионного переноса на регулярных сетках при наличии погранслоев //Вычисл. технологии, т. 8, совместный выпуск, часть 3, 2003. -- С.15-20.
2. Паничкин А.В. Моделирование течений вязкой жидкости по улучшенным схемам с алгоритмом уменьшения схемной вязкости // Вычисл. технологии, 6, спец. выпуск, часть 2, 2001. -- С.485-495.
Примечание. Тезисы докладов публикуются в авторской редакции
Ваши комментарии Обратная связь |
[Головная страница] [Конференции] |
© 1996-2000, Институт вычислительных технологий СО РАН, Новосибирск
© 1996-2000, Сибирское отделение Российской академии наук, Новосибирск