XVII школа-семинар по численным методам механики вязкой жидкости
На основе метода расщепления по физическим процессам (Марчук Г.И., 1988) предлагается численный метод решения начально-краевых задач для уравнений свободной конвекции в замкнутых областях с твердыми границами.
Рассматривается два варианта метода. В одном случае расщепление по физическим процессам осуществляется для уравнения переноса вихря, а в другом - для уравнения сохранения импульса. Таким образом, этап конвекции формулируется либо в переменных "вихрь - функция тока", либо в физических переменных (компоненты скорости). Этап диффузии в обоих случаях формулируется в переменных "вихрь - функция тока".
Численная реализация этапа конвекции осуществляется с помощью двуциклической редукции на основе схем Кранка-Николсона, что обеспечивает энергетическую нейтральность искомых полей. Для расчета уравнений на этапе диффузии используется модификация так называемого "двухполевого" метода (Тарунин Е.Л., 1990), позволяющая разделить задачи вычисления вихря и функции тока и, при этом, точным образом (в разностном смысле) выполнить условия прилипания на твердых стенках (Воеводин А.Ф., Юшкова Т.В., 1999). Этап диффузии реализуется на основе схем повышенного порядка точности.
Результаты расчетов по предложенным методикам сравниваются между собой и с результатами других авторов. Сравнение демонстрирует эффективность применения метода для расчета конвективных течений вязкой жидкости.
Примечание. Тезисы докладов публикуются в авторской редакции
Ваши комментарии Обратная связь |
[Головная страница] [Конференции] |
© 1996-2000, Институт вычислительных технологий СО РАН, Новосибирск
© 1996-2000, Сибирское отделение Российской академии наук, Новосибирск