Разработка программного комплекса для численного моделирования и визуального представления вязких теплопроводных течений в каналах с фиксированными граничными условиями

Тип работы:
Реферат
Предмет:
Механика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Секция «Моделирование систем»
Разработка программного комплекса для численного моделирования и визуального представления вязких теплопроводных течений в каналах с фиксированными
граничными условиями
Александр Бондарев Институт прикладной математики им. М. В. Келдыша РАН
Москва, Россия bond@keldysh. ru
Аннотация
В работе описаны основные свойства программного комплекса для решения задач численного моделирования и визуального представления течений вязких теплопроводных сред в каналах с фиксированными граничными условиями. Рассматриваются различные типы задания граничных и начальных условий. Приведены результаты тестовых расчетов.
Ключевые слова: вязкие теплопроводные течения, программный комплекс
1. Введение
Данная работа посвящена начальной фазе общего проекта, направленного на решение актуальной проблемы — управления свойствами нестационарных вязких теплопроводных течений в каналах и трубках с гибкими стенками с помощью комплексного решения прямых, сопряженных и обратных задач методами численного моделирования и научной визуализации. Разрабатываемые подходы, методы и алгоритмы позволят получать принципиально новые результаты, используемые для решения технологических инженерных задач в широком диапазоне применения.
Управление свойствами нестационарных вязких теплопроводных течений в каналах и трубках с гибкими стенками является одной из актуальных современных задач. Одной из важных целей подобного управления является выявление закономерностей возникновения различных аномалий в течениях, что имеет важное
прикладное значение в ряде задач современных индустриальных технологий и в ряде задач прикладной физиологии. Управление свойствами течения достигается с помощью комплексного решения прямых, сопряженных и обратных задач методами численного моделирования и научной визуализации. Об актуальности работы свидетельствует большой интерес, который выражен многочисленными публикациями по данной проблематике в таких журналах как Computer& amp-Mathematics with Application, Computer Methods in Appl. Mech. and Eng. и др. Также об актуальности тематики говорят обсуждения на конференциях NHTC ASME и многих других форумах, посвященных тепломассообмену и решению обратных задач.
Данная работа относится к начальной стадии проекта — разработке программных средств для решения прямых задач моделирования течения вязких теплопроводных сред в каналах с фиксированными граничными условиями.
Вязкие теплопроводные течения являются объектом изучения на протяжении многих десятилетий. Основные принципы и положения, посвященные этим течениям, достаточно полно изложены в [1−3]. Основная задача, решаемая на данной стадии проекта, — разработка программных средств моделирования и визуализации течений в канале, позволяющих получать результаты количественно и качественно согласующимися с известными решениями и закономерностями, изложенными в вышеупомянутых работах.
2. Математическая модель
В общем случае математическая задача о движении сплошной среды с учетом эффектов вязкости, теплопроводности и сжимаемости сводится к решению полной системы нестационарных уравнений Навье-Стокса. Для двумерного случая эта система записывается в декартовой (_/ = 0) или цилиндрической (_/ = 1) системе координат следующим образом:
Э/ Э/
— + и- + Эt Эх
v —
1 •• Ц
R-Ja2 —
Re р у
Э/
Э • Э/ Э , —
— ax ц^- ±-a2 Ц^- & gt- + F
Эу Re р I Эх Эх Эу
a2 ц Э/1 + F 2 Эу J
(2. 1)
где / = (р, u, v, T) — вектор искомых газодинамических функций.
. 4 g
ах = (0, —, 1,-P-)
3 Pr
. 4 g
а2 = (0,1, — ,~P~)
3 Pr
F = (Fp FuFVFT)
1
Рр _- р
Гди + дг^ дх ду
-
У
р _ 1 дР+
р дх
д дг 2 д ду.
ду дх 3 дх ду
2
3 Эх
У
т ду
У дх
р
р _ дР + Гэ_ ди — 2 д дм-. 2 у /
у р ду I дх ду 3 ду дх 3 у
дт+2 т
ду у
1 -а-| ^ р
рт _-(у — 1) т м? У (г — 1) т
дм ду .у — + - + ] -
дх ду у
+
+
р
ди 2 2
дх
+
ду
ду
+ 3
2
у
у
+
ду ди
— + -
дх ду
2
3
ди ду у -±+3-
дх ду у
Здесь х, у — декартовы (3 _ 0) или цилиндрические (3 _ 1) координаты, * - время, и, у — составляющие компоненты вектора скорости, р -плотность, Р-давление, Т -температура, у — соотношение удельных теплоемкостей, ц — коэффициент динамической вязкости, 1 — коэффициент теплопроводности. (При записи уравнения энергии считается, что 1 _ ц)•
Система (2. 1) записана в безразмерной форме, переход к которой осуществляется путем отнесения размерных величин к характерным геометрическим размерам и величинам в набегающем потоке. Также по характерным размерам и параметрам набегающего потока определяются характерные числа задачи, такие как число Маха М?, число Рейнольдса Яе?, число Прандтля Рг.
Система (2. 1) замыкается соотношениями для коэффициентов вязкости и теплопроводности 1 _ ц _ ц (Т) и уравнением состояния в безразмерной форме
_ рт
Р_
ум?
При решении задач в качестве определяющих газодинамических
функций выбираются Р, и, у, Т. 3. Метод решения
Задача решается в предположении, что для каждой из искомых газодинамических функций Р, и, у, Т соответствующее уравнение системы (2. 1) является определяющим. Каждое из уравнений (2. 1) представляется в общем виде
2
2
2
Э/. «Э/. а У
Э
--+ а--+ 6- = с-
Эt Эх Эу Эх
* У
Эх
Э
+ е Э& quot- ЭУ
'- У
я-

(3. 1),
где а, Ь, с, й, в^ - соответствующие коэффициенты. В Н включаются члены со смешанными производными, а также члены с производными по тем функциям, для которых данное конкретное уравнение из (2. 1) не является определяющим.
Для решения системы уравнений (2. 1) используется метод переменных направлений. Для примера приведем запись уравнений, рассматриваемых для направления вдоль оси X:
Э/ Э/ Э
— + а — - с — Эt Эх Эх
& lt- У
х
= А — А = -ЬУ + е — Эу Эу
?
У
Эу

(3. 2)
Для решения подобных уравнений по каждому координатному направлению в программном комплексе реализованы общие алгоритмы применения неявных конечно-разностных схем, разрешаемых путем скалярной прогонки. Таким образом, программный комплекс позволяет обеспечить численную реализацию большинства неявных конечно-разностных схем первого и второго порядка аппроксимации по времени и пространству, записанных на соответствующем сеточном шаблоне, согласно [4,5]. Также программный комплекс позволяет обеспечить численную реализацию многих гибридных конечно-разностных схем, описанных в [6]. Гибридные конечно-разностные схемы являются согласно [6] в общем случае линейными комбинациями разностных схем, обладающих разными вычислительными свойствами, и применяются с тем, чтобы максимально удовлетворить потребности вычислителя.
При проведении тестовых расчетов применялась гибридная неявная конечно-разностная схема, подробно представленная в [7,8,9]. Данная конечно-разностная схема на каждом направлении обладает вторым порядком аппроксимации по времени и пространству с добавочным членом, выполняющим функции искусственной вязкости. Управление искусственной вязкостью осуществляется с помощью весовых коэффициентов. Применяемая схема является безусловно устойчивой. Учет нелинейности коэффициентов организован с помощью экстраполяции нелинейных коэффициентов по предыдущим временным слоям.
4. Граничные условия
Для общей постановки задачи используется несколько разных типов граничных условий. В программном комплексе эти граничные условия реализуются по модульному принципу. В блоке граничных условий GRUS программные модули расположены по принципу: каждому типу граничных условий соответствует свой отдельный программный модуль. Это позволяет пользователю в управляющем блоке UPR легко присваивать различные типы граничных условий геометрическим границам в соответствии с конкретными задачами.
В данный период в блоке граничных условий реализованы следующие типы граничных условий:
а) для входной границы канала:
— условие задания характерного набегающего потока
Up = U"(y, t), = 0, Рр = Р_, Тр = Т_(y, t)
Эти условия задаются с помощью известных аналитических и экспериментальных зависимостей-
— условия в набегающем потоке, полученные из других расчетов
Up = U (y, t), Vrp = V (y, t), Рр = P (y, t), Тр = T (y, t)
Эти условия нужны для того, чтобы полученные в предыдущих расчетах профили газодинамических величин задать в качестве граничных условий на входной границе.
б) для выходной границы канала:
— условие типа линейной экстраполяции (или «мягкие» граничные условия)
Э2 f '-
^ = 0, f = (Р, u, v, Т) —
— условия «сноса» по течению
f = 0, f = (Р, u, v, Т) —
Эп
— условия входного потока на выходной границе
игр = U^(y, t), Vrp = 0, Рр = Р_, Тр = Т_(y, t).
в) для оси симметрии канала — условие симметрии
Эй = v = Эр=ЭТ = о*
Эп гр Эп Эп
г) для стенки канала:
— условие прилипания для стенки с постоянной температурой
& quot-гр = ^ = Щ- = 0 '- Тгр =сп —
— условие прилипания для теплоизолированной стенки
ЭР ЭТ
и* = ^ = эТ = эТ=0'-
— условие прилипания для стенки с теплоотводом или подводом тепла-
ЭР 0 ЭТ ЭТ (ч
& quot-гр = Угр =-г = 0, -г =-г (х, t).
Эп Эп Эп
Дальнейшие типа граничных условий могут добавляться в блок граничных условий в виде отдельных программных модулей по мере разработки. В настоящее время вопросы подвижности границ, а также вопросы моделирования граничных условий эластичности и частичной проницаемости границ находятся в стадии разработки.
5. Начальные условия
Блок задания начальных условий построен по модульному принципу аналогично боку задания граничных условий. В данный период реализованы следующие типы задания начальных условий:
— условия покоя
и (х, у^) = V (. х, у^) = 0, Р (. х, у^) = Р?, Т (. х, у^) = Т?(у^)
Это т.н. «плохие» начальные условия. Они нужны для того, чтобы контролировать и тестировать процессы эволюции и установления течения.
— условия «сноса»
и (х, у^) =, V (х, у^) = 0, Р (х, у^) = Р., Т (х, у^) =)
Эти начальные условия представляют собой «снос» граничных условий на входной границе вниз по потоку на все поле течения.
— условия импорта поля течения из других расчетов или экспериментов
и = и (х, у^), V = V (х, у^), Р = Р (х, у^), Т = Т (х, у^)
6. Средства анализа и визуального представления решения
В программном комплексе предусмотрен блок анализа решения ANALYSER. Его основная роль предназначена для решения обратных и сопряженных задач. На данном этапе он используется для контроля качества решения. Данный блок контролирует локальные экстремумы, анализируя поля первых и вторых производных определяющих газодинамических величин в расчетной области.
В настоящее время для рабочего контроля решения также используются инструменты визуализации, позволяющие получать распределения величин P, u, v, T как вдоль выбираемых пользователем линий, так и во всей расчетной области. Эти инструменты рабочего визуального контроля реализованы с помощью системы ГРАФОР [10]. В работе использовалась современная версия [11,12] реализованная на базе Windows на основе системы Microsoft Fortran Power Station 4.0 и Microsoft Developer Studio. Учитывая, что текущая версия описываемого программного комплекса также реализована под Microsoft Fortran Power Station 4.0, это довольно удобно.
Вопрос выбора платформы и системы визуализации для проекта в целом находится в настоящее время в стадии обсуждения и разработки.
7. Проведение тестовых расчетов
Для проверки работоспособности программного комплекса был проведен ряд стандартных тестовых расчетов. В качестве тестов рассматривались варианты течений в широком диапазоне чисел M? и Re?. Проверялись все вышеописанные типы задания граничных условий и начальных условий в счетной области. В первую очередь контролировалась физичность результатов и их количественное и качественное соответствие известным решениям из [1,2,3]. Практически все тесты показали хорошую согласованность с известными аналитическими и экспериментальными зависимостями [1,2,3].
Приведем пример результатов для случая ламинарного течения при M? = 3 и Re?= 1000. На рис. 1 представлены линии постоянных значений числа Маха для установившегося течения.
Рис. 1.
На рис. 2 представлены линии постоянных значений числа Маха для установившегося течения.
Рис. 2.
На рисунках 3, 4, 5 представлены распределения по направлению У величин и, Р, Т соответственно. Это распределения газодинамических величин на выходе из канала.
Рис. 3.
со)
о (X) у /
о щ*
о СМ
сз о N
& lt-=1.
ИШ
1051
по-.
1151
1201.
1251.
1301

Рис. 4
Рис. 5
На рисунке 6 представлено распределение давления вдоль стенки канала.
Рис. 6
Пример этих данных показывает, что результаты расчетов хорошо согласуются с известными аналитическими и эмпирическими зависимостями из [1,2,3]. Полученные в тестовых расчетах данные также хорошо согласуются с результатами расчетов, приведенных в [5,15].
Следует заметить, что тестовые расчеты показали необходимость при применении гибридных неявных конечно-разностных схем для оптимизации вычислительных свойств и точного определения границ весовых коэффициентов в рассматриваемом классе задач применять численно-визуальный методологический подход, подробно описанный в [13,14].
Применение подобного подхода позволит резко повысить эффективность применения программного комплекса к конкретным классам задач.
8. Заключение
Разработан первый вариант прототипа программного комплекса для решения прямых задач моделирования течения вязких теплопроводных сред в каналах с фиксированными граничными условиями. Данный прототип является начальной фазой общего проекта, направленного на решение актуальной проблемы -управления свойствами нестационарных вязких теплопроводных течений в каналах и трубках с гибкими стенками с помощью комплексного решения прямых, сопряженных и обратных задач методами численного моделирования и научной визуализации. Проведено тестирование и численное обоснование на известных решениях. Тестирование показало достаточно хорошее согласование с известными качественными и количественными параметрами подобных течений. В перспективе рассматривается разработка алгоритмов подвижности границ, а также граничных условий эластичности и частичной проницаемости границ.
Список литературы
[1] Шлихтинг Г. Теория пограничного слоя. М., Наука, 1974, 711 с.
[2] Абрамович Г. Н. Прикладная газовая динамика. М., Наука, 1976, 888 с.
[3] Бондарев Е. Н., Дубасов В. Т., Рыжов Ю. А., Свирщевский С. Б., Семенчиков М. В. Аэрогидромеханика. М., Машиностроение, 1993, 608 с.
[4] Самарский А. А., Попов Ю. П. Разностные методы решения задач газовой динамики. М., Наука, 1980, 352 с.
[5] Пасконов В. М., Полежаев В. И., Чудов Л. А. Численное моделирование процессов тепло- и массообмена. М., Наука, 1984, 288 с.
[6] Куликовский А. Г., Погорелов Н. В., Семенов А. Ю. Математические вопросы численного решения гиперболических систем уравнений. М., Физматлит, 2001, 608с.
[7] Бондарев А. Е. Численное решение уравнения Бюргерса в области высоких градиентов. Препринт ИПМ им. М. В. Келдыша РАН, М., N12, 1990, 13с.
[8] Бондарев А. Е. Разработка метода численного исследования отрывных течений вязкого газа. Отчет ИПМ им. М. В. Келдыша РАН, М., N170−90, 1990, 16с.
[9] Бондарев А. Е. Оптимизация вычислительных свойств гибридных разностных схем с помощью методов визуального представления для ламинарных и турбулентных течений в дальнем следе. Препринт ИПМ им. М. В. Келдыша РАН, М., N21, 2007, 21с.
[10] Баяковский Ю. М., Галактионов В. А., Михайлова Т. Н. «ГРАФОР. Графическое расширение Фортрана» М., Наука, 1985.
[11] Базаров С. Б., Баяковский Ю. М. «Комплекс графических программ ГРАФОР в среде WINDOWS» Препринт ИПМ им. М. В. Келдыша РАН, № 30, 2000, 10 с.
[12] Базаров С. Б., Баяковский Ю. М., Сейдалиева Ф. Ф., Скачков А. Ю. «Адаптация комплекса графических программ ГРАФОР в операционных системах WINDOWS и LINUX» Препринт ИПМ им. М. В. Келдыша РАН, № 27, 2002, 16 с.
[13] Бондарев А. Е. Применение методов научной визуализации для оптимизации вычислительных свойств конечно-разностных схем. Препринт ИПМ им. М. В. Келдыша РАН, № 79, 2006, 17 с.
[14] Бондарев А. Е. Применение методов научной визуализации для оптимизации вычислительных свойств конечно-разностных схем. Препринт ИПМ им. М. В. Келдыша РАН, № 79, 2006, 17 с.
[15] Елизарова Т. Г. Квазигазодинамические уравнения и методы расчета вязких течений. М., Научный мир, 2007, 352 с.

ПоказатьСвернуть
Заполнить форму текущей работой