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

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ, Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 48211. ISSN 1994-0408

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

# 08, август 2014

DOI: 10.7463/0814.0725648 Клюквин А. Д.1,а

УДК 532.5 : 519.6

Проведено исследование влияния зависимости теплофизических воздуха от температуры на точность решения задач теплообмена в турбулентном потоке при использовании различных методов осреднения уравнений движения в форме Навье-Стокса. Получена форма дополнительных членов, возникающих в осредненных уравнения Навье-Стокса в случае учета пульсаций теплофизических свойств воздуха. На примере затопленной подогретой струи проведена оценка погрешностей при определении интенсивности процессов переноса, возникающих при пренебрежении зависимостью теплофизических свойств воздуха от температуры осредненного потока. Проведено сравнение эффективности использования различных методов осреднения уравнений Навье-Стокса при учете переменности теплофизических свойств воздуха.

Введение

Одной из насущных задач современной теплотехники является численное моделирование процессов теплообмена в турбулентных течениях. При решении данного класса задач необходимо учитывать, что в исследуемых течениях возникает поток энергии, непрерывно передаваемой от осредненного потока к мелкомасштабным пульсациям, в которых, в основном, происходит диссипация энергии (ее переход в теплоту) [1]. Таким образом, для наиболее точного описания процессов теплообмена в турбулентном потоке необходима подробная информация о пространственном положении энергосодержащих структур турбулентного потока.

Однако эта информация принципиально не может быть получена с помощью подходов RANS (Reynolds-averaged Navier-Stocks equations - уравнения Навье-Стокса, осред-ненные по Рейнольдсу) или FANS (Favre-averaged Navier-Stocks equations - уравнения На-вье-Стокса, осредненные по Фавру) [1], в рамках которых происходит моделирование

всех, в том числе и самых крупных масштабов турбулентности. Основное ограничение данного подхода заключается в том, что крупные вихри имеют анизотропную структуру и сильно зависят от геометрических параметров течения [2].

К примеру, в работе [3] показано, что RANS дают неудовлетворительный результат в возвратно-рециркуляционной зоне исследуемого течения.

Наиболее полную информацию о турбулентном потоке можно получить, применяя прямое численное моделирование (Direct numerical simulation, DNS), однако при существующем сегодня уровне развития вычислительной техники этот подход может быть реализован только в условиях простой геометрии и низких чисел Рейнольдса.

Учитывая, с одной стороны, важность детального описания турбулентного движения жидкости при анализе процессов теплообмена, а, с другой, - возрастающую при детализации вычислительную сложность, возникает необходимость в подходах, более универсальных и строгих, чем RANS/FANS, и менее дорогостоящих в вычислительном плане, чем DNS.

Одним из таких подходов является моделирование крупных вихрей (Large eddy simulation, LES) - метод, основанный на наложении на осредненную картину движения крупных и мелких вихрей. Крупные вихри, находящиеся под прямым воздействием граничных условий и несущие в себе максимум рейнольдсовых напряжений, рассчитываются. Мелкомасштабная турбулентность считается изотропной и имеющей универсальные характеристики, а потому более легко поддающейся моделированию. Взаимодействия между крупными и мелкими вихрями аппроксимируются только по крупным вихрям с использованием подсеточных моделей. Иначе говоря, принимается гипотеза о статистической независимости крупных и мелких вихрей.

В то время как DNS отображает весь диапазон размеров вихрей, в LES считаются наиболее важными крупные вихри.

Решение, полученное с помощью LES, содержит более подробную информацию по сравнению с решением на основе уравнений Рейнольдса, и, следовательно, позволяет более точно описать процесс теплообмена с окружающей средой.

Однако, несмотря на меньшую вычислительную сложность LES в сравнении с DNS, применение этого метода также ограничивается относительно низкими числами Рей-нольдса. Так, в работе [4] показано, что для расчета работы турбины с помощью LES потребовалось в 1750 раз больше машинного времени, чем с помощью RANS.

Постановка задачи

В отличие от DNS, при решении задач вычислительной гидродинамики с помощью LES или RANS/FANS теплоемкость, теплопроводность и вязкость воздуха обычно принимаются постоянными в рамках исследуемой системы [5]. Однако опытные данные показывают, что эти характеристики зависят от температуры и давления воздуха (зависимость от давления обычно оказывается очень слабой [6]) в системе и могут существенно колебаться в случае больших градиентов температуры.

В данной работе рассматривается влияние переменности теплофизических параметров воздуха на результат фильтрации и осреднения уравнений Навье-Стокса. Причем это влияние разбивается на две составляющие:

появление дополнительных моделируемых слагаемых в осредненных и фильтрованных уравнениях Навье-Стокса, вызванных пульсациями параметров потока;

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

В качестве исходной принимается система уравнений Навь-Стокса в индексной форме записи (под уу, уу, , у п понимается суммирование по повторяющемуся индеку,к,т,п=1,2,3). 1. Уравнение неразрывности

др + д

£ р )=0

дг дхЛг ]> (1)

где р - плотность газа, кг/м ; г - время, с; ху - декартова координата, м; V у - проекция

вектора скорости газа, м/с; 2. Уравнение движения

а р),) = +дту

Г а,. ~ ~ Л

где Ти = Ц

дх,. дх 3 у дх,

V 7 1 к

fi - проекция единичного вектора объемных сил; Н/м , р - статическое давление газа,

Па; Ту - составляющая тензора вязких напряжений, Па; Ц- молекулярная вязкость, Па*с;

дре)^^)_ д{ру]) дд^+ ду)

дг дху дху дху дху & (ъ)

_ , дТ

где Яу = Я

дху &

е - удельная внутренняя энергия газа, Дж/кг; Я у - проекция теплопроводной компоненты теплового потока, Дж/мъ; Sу - составляющая тензора скоростей деформации, с-1; Я - коэффициент теплопроводности, Вт/(м*К); Т - температура газа, К.

Осреднение уравнений Навье-Стокса

В рамках данной работы наибольший интерес представляют течения воздуха с существенными градиентами температуры, где воздух ведет себя как сжимаемый газ. При работе с течениями такого вида может использоваться как осреднение по Рейнольдсу, так и осреднение по Фавру. Однако при учете сжимаемости воздуха последнее позволяет получать уравнения в более компактной форме, чем при использовании первого и с меньшим количеством дополнительных неизвестных [1].

В рамках данного раздела будем использовать следующие обозначения: щ, р -функции, осредненные по Фавру, щ", р" - величины пульсаций функций, осредненных по Фавру, щ , р - функции, осредненные по Рейнольдсу, щ&, р - величины пульсаций функций, осредненных по Рейнольдсу, Щ - величина, осредненная по всем потоку.

Проведя математические выкладки (аналогично работам, [1] и [6]), получим осредненные по Фавру уравнения Навье-Стокса:

Р ) = о;

а дхл^ >> <4)

—(p W^M = f.-dp+A

Tij-pvivi

^ [Р+P) v -dj Qj Р),

dx■ dx 3 dxh iJ

V J k У

i dT qJ Щ

Qj =-cp p(vjT+VjT) ,

Dj=\\{vjvkvk-~vJvkvk),

Qj - проекция конвективного подсеточного теплового потока, Дж/м3; cр - удельная теплоемкость воздуха при постоянном давлении, Дж/(кг*К).

Получение дополнительных членов в FANS

Уравнения (4) - (6) были получены без учета изменения теплофизических свойств воздуха при изменении его температуры. Рассмотрим слагаемые, зависящие от этих свойств воздуха, и определим дополнительные члены, возникающие при учете зависимости теплопроводности, теплоемкости и вязкости воздуха от его температуры.

Проекция теплопроводной компоненты теплового потока описывается уравнением Фурье

Ч ~=~лщ

Представим величины, входящие в правую часть данного уравнения как сумму средней величины и пульсационной составляющих по Рейнольдсу

= —(л+л)

Т+Т&)

дх,^л дх дх■ дх дху у у у у

Учитывая, что (р&у = 0, выделим из выражения, полученного выше, пульсацион-ную ( ) и осредненную ( Ч у ) по Рейнольдсу составляющие Чу :

лд^+Л&дТ

дХ: дХ:

лё+лдг ^

дхдху у У

Средние и пульсационные величины по Рейнольдсу и по Фавру связаны следующим соотношением [1]:

С его помощью получим среднюю по Фавру величину проекции теплопроводной

компоненты теплового потока

-г дТ „дТ& рЛ дТ ЛрдТ& дх- дх- р дх■ р дх7

у у ^ у Г- у у

Определим величину дополнительных членов, возникающих при учете пульсаций теплопроводности воздуха:

Чу - Чу

,дТ& + рЛ дТ

дхр дх

где Ч

_ - проекция теплопроводной компоненты теплового потока без учета пульсаций

теплопроводности.

Г — Гу у

дху дх■ 3 дхк V у

дх дх: 3 дх,

V у & К У

где Tjj

_ - составляющая тензора вязких напряжений без учета пульсаций молекулярной

вязкости.

с - - - л

где Qj _ - конвективная компонента теплового потока без учета пульсаций теплоемкости при постоянном давлении.

Фильтрация уравнений Навье-Стокса

Как было сказано во введении, одним из подходов к расчету параметров турбулентных течений является LES, при использовании которого происходит разделение вихрей на крупные и мелкие. Математическим эквивалентом такого разделения является операция фильтрации. Ее свойства (общие, а также в приложении к уравнениям Навье-Стокса) подробно рассмотрены в работе [1].

Применим оператор фильтрации к уравнениям Навье-Стокса (в рамках данного раздела и далее под у будем понимать величину, фильтрованную со взвешиванием по плотности, под у - величину, фильтрованную без взвешивания по плотности, под у - величину, осредненную по всем потоку, а под у& - пульсации соответствующих величин). Выполним математические преобразования аналогично работе [5], получим фильтрованные уравнения Навье-Стокса.

др+ д (— dt дх

А(ри£№=jf+А

дх,У&> дх, J& дх, дх,

vv- — vv ■

где Ту - величина подсеточных напряжений;

д(ре) | д(рвУу) = д~ру] _д_ъ +<Т)_A(Q )

гН дх. дх. дх ■ дх. дх Л 3 3&

р( Rij+l+sj )! &

Чу =—Л

Ц = V/V/ -\\>,У,,

= & + V Ч , 0 = —^Р(^Т+) ■

Получение дополнительных членов при фильтрации уравнений НавьеСтокса

Аналогично осреднению по Фавру получим дополнительные члены, возникающие при фильтровании (при выводе моменты пульсаций второго порядка и выше считались пренебрежимо малыми):

Чу— Чу

г \\ , дТ

о 8кк$у

Выражения (7) - (9) и (13) - (15) представляют собой дополнительные члены, которые возникают в случае учета пульсаций теплофизических свойств газа при осреднении и фильтрации уравнений Навье-Стокса, соответственно.

Проблема влияния данных членов на точность получаемых расчетов в литературе освещена весьма слабо, поэтому количественная оценка их важности затруднительна. Однако структура полученных выражений позволяет утверждать, что их учет будет наиболее оправдан в случае расчета интенсивно охлаждаемых потоков газа, где пульсации теплоемкости, теплопроводности и вязкости могут достигать довольно больших значений.

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

В предыдущем разделе проводился анализ влияния дополнительных членов, возникающих при осреднении и фильтрации уравнений Навье-Стокса вследствие учета пульсаций теплофизических свойств воздуха. В данном разделе, будет исследована зависимость теплофизических параметров воздуха от температуры осредненного потока.

Проведем оценку изменения теплофизических параметров воздуха на примере истечения неизотермической дозвуковой струи из круглого сопла, расчет которого приведен в [5]. Для этого, пользуясь данными [7], разложим теплоемкость, теплопроводность и вязкость воздуха в ряд Тейлора в районе 400 К (данная температура принимается за среднюю по всему потоку). В результате получим

Л(Т& )-105 = 3,365 + 7,096-10-3Т& + 3,884-10-6Т&2-1,735Т3 +..., ср {Т )-10-3 = 1,014+1,250-10-4Т - 3,187-10-7Т2 + 4,731Т3 +..., М(Т& )-105 = 4,846 + 4,155-10-3Т - 3,643-10-6Т2 - 3,575Т3 +...,

где Т& - отклонение температуры от среднего значения.

Определим в процентном соотношении величину изменения теплофизических свойств газов в случае колебаний температуры на 160 К (характерная величина изменения осредненной температуры).

с & 2& &

ср. = 3,7%, 4 = 31,3%, 4 = 32,0% с р Л 4

Из уравнений (7) - (9) и (13) - (15) следует, что максимальная относительная погрешность, вызванная переменностью теплофизических свойств воздуха, имеет величины:

* \\М=М

= 4 = 32,0%, 4

а -а\\

ю = ср_ = 3,7%,

= Ат = 31,3% Я

где Т

составляющая тензора вязких напряжений, конвективная и

теплопроводная компоненты теплового потока при постоянных вязкости, теплоемкости и давлении, соответственно.

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

решения

Начнем с анализа вязких напряжений. Как видно из выражения (16), принятие гипотезы о постоянстве вязкости воздуха ведет к недооценке величины вязких напряжений вплоть до 32 %. Данный эффект наиболее сильно проявляется в непосредственной близости от границы ядра подогретой струи, так как при удалении от него происходит активное перемешивание масс газа различной температуры, и она довольно быстро выравнивается. Однако, как показано в [8], в данной области силы молекулярной вязкости малы в сравнении с силами инерции. Следовательно, учет зависимости их величины от температуры воздуха не приведет к значительному повышению результатов расчета.

Из формулы (17) видно, что учет зависимости теплоемкости от температуры газа приводит к уточнению величины конвективной компоненты теплового потока лишь на 3,7%, что находится в рамках допустимой погрешности. То есть в данном случае гипотеза о постоянстве теплоемкости воздуха при постоянном давлении оказывается оправданной.

Далее рассмотрим переменность коэффициента теплопроводности. Как можно видеть из (18), в случае принятия его постоянным по всему течению, модуль теплопроводной составляющей теплового потока может быть недооценен на 31,3%. С одной стороны, даже такие большие колебания интенсивности процессов теплопроводности не должны оказывать заметного влияния на процессы теплообмена в газе при больших числах Рей-нольдса (в этом случае наиболее важную роль играют процессы конвективного теплообмена).

Но в то же время, одним из специфических свойств рассматриваемого случая является доминирующая роль крупных турбулентных структур в зоне смешивания подогретой струи газа с окружающей средой. Таким образом, большие массы газа оказываются связанными крупными вихрями и более слабо подвергаются процессам теплообмена за счет уменьшения конвективной составляющей теплового потока. Следовательно, повышается важность процессов теплопередачи, и возникает возможность существенного уточнения решения задач теплообмена с помощью учета температурной зависимости теплопроводности воздуха. Однако проведение такого уточнения практически невозможно с помощью RANS/FANS, так эти методы плохо описывают зоны газа с крупномасштабными турбулентными структурами.

То есть при попытках уточнения расчетов процессов теплообмена в турбулентном потоке путем учета реальных теплофизических свойств воздуха перед нами встает вопрос о переходе от RANS/FANS к более сложным методам, что вызвано необходимостью учета поведения крупномасштабных турбулентных структур потока.

Однако использование более сложных методов расчета влечет за собой увеличение требуемых вычислительных мощностей, что не всегда является возможным и целесообразным как с технической, так и с экономической точки зрения.

Заключение

В рамках данной работы был проведен качественный анализ влияния температурной зависимости теплофизических свойств воздуха на форму осредненных и фильтрованных уравнений Навье-Стокса и точность получаемого численного решения.

В результате была определена форма добавочных членов в исследуемых уравнениях, обусловленных пульсациями вязкости, теплопроводности и теплоемкости воздуха.

Показано, что учет распределения теплофизических свойств по осредненному потоку может оказать существенное влияние на результат решения задач теплообмена в турбулентном течении. При этом, согласно проведенным расчетам, погрешности, вызванные пренебрежением этим влиянием, могут достигать

Список литературы

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. ISSN 1994-0408

Analysis of Influence of the Thermal Dependence of Air

Thermophysical Properties on the Accuracy of Simulation of Heat

Transfer in a Turbulent Flow in Case of Applying Different

Methods of Averaging Navier-Stokes Equations

# 08, August 2014

DOI: 10.7463/0814.0725648

A.D. Kliukvin1a

aaklyukvin@yandex.ru

There is theoretically investigated the influence of thermal dependence of air thermophysical properties on accuracy of heat transfer problems solution in a turbulent flow when using different methods of averaging the Navier-Stokes equations.

There is analyzed the practicability of using particular method of averaging the Navier-Stokes equations when it&s necessary to clarify the solution of heat transfer problem taking into account the variability of air thermophysical properties.

It&s shown that Reynolds and Favre averaging (the most common methods of averaging the Navier-Stokes equations) are not effective in this case because these methods inaccurately describe behavior of large scale turbulent structures which strongly depends on geometry of particular flow. Thus it&s necessary to use more universal methods of turbulent flow simulation which are not based on averaging of all turbulent scales.

In the article it&s shown that instead of Reynold and Favre averaging it&s possible to use large eddy simulation whereby turbulent structures are divided into small-scale and large-scale ones with subsequent modelling of small-scale ones only. But this approach leads to the necessarity of increasing the computational power by 2-3 orders.

For different methods of averaging the form of additional terms of averaged Navier-Stokes equations in case of accounting pulsation of thermophysical properties of the air is obtained.

On the example of a submerged heated air jet the errors (which occur when neglecting the thermal dependence of air thermophysical properties on averaged flow temperature) in determination of convectional and conductive components of heat flux and viscous stresses are evaluated. It&s shown that the greatest increase of solution accuracy can be obtained in case of the flows with high temperature gradients.

Finally using infinite Teylor series it&s found that underestimation of convective and conductive components of heat flux and viscous stresses can reach 3,7%, 31,3% and 32%, respectively.

References