Уравнение теплопроводности решалось по неявной конечно-разностной схеме с использованием схемы с опережением для аппроксимации конвективного члена по методике, рассмотренной в работе [19]. Шаги Ах и Аг увеличивались по геометрической прогрессии от Ах = 1 км у оси до Ах = 70 км на правой границе области и от Аг = 1 км у поверхности до Аг = 6 км у нижней границы области. Шаг по времени выбирался из условия обеспечения устойчивости решения. Точность решения проверялась сравнением с аналитическими решениями (остывание однородного полупространства и теплообмен двух блоков с отличающимися температурами), а также сравнением последовательных решений, полученных для различных А/, Ах, Аг.
Скачать страницу
[Выходные данные]