Математическое моделирование течений вязкой жидкости


Математическая модель


осредненного турбулентного движения

Пусть имеем систему уравнений пограничного слоя:

Математическая модель
     (2.1)

Так как первый член в правой части первого уравнения системы (2.1) записан как

Математическая модель
, а не
Математическая модель
, то надо оставить  и второе уравнение
Математическая модель
, чтобы сохранилась корректность системы уравнений пограничного слоя.

Для описания турбулентного движения Рейнольдс предложил следующий  прием. Регистрируя во времени скорости и давления в данной точке  потока, можно их представить как:

Математическая модель
,

где

Математическая модель
 - действительно существующие в потоке мгновенные (актуальные) проекции скорости и давления;
Математическая модель
 - осредненные во времени их значения;
Математическая модель
 - пульсации проекций скорости и давления.

Под осредненным значением параметра понимается обычное интегральное  среднее по времени t за промежуток T , называемый, периодом осреднения:

Математическая модель
Математическая модель
Математическая модель
.

В турбулентном движении добавляется пульсационная составляющая скорости (рис.10), в результате чего наблюдается вихревое движение, при котором сопротивление значительно возрастает. Таким образом, турбулентное течение обладает бóльшим сопротивлением по сравнению с ламинарным движением.

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

Математическая модель

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

Математическая модель
.

(Если ограничиться несжимаемой однородной жидкостью, то r=const и, следовательно,  

Математическая модель
.

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

В настоящее время турбулентное движение принято характеризовать осредненным по времени значением величин. В уравнениях сохранения массы, количества движения и энергии в потоке вязкой жидкости истинные (мгновенные) величины заменяются осредненными во времени их значениями следующим образом. Истинные величины в данной точке турбулентного потока раскладываются  на осредненные и пульсационные их значения, что соответствует физическому представлению турбулентного движения. Тогда уравнения неразрывности, движения и энергии для осредненного турбулентного движения несжимаемой жидкости в общем случае получаются из исходных уравнений после замены в них истинных значений переменных осредненными их значениями и пульсациями с последующим осреднением этих параметров по времени. При введении в действие новых переменных добавляется три неизвестных:

Математическая модель
, и задача переходит в разряд неопределенных. Для устранения неопределенности и применяется усреднение по времени.

Рассмотрим решение задачи. Возьмем, например, уравнение:
Математическая модель
.

Проведя операцию осреднения, его можно записать следующим образом:

Математическая модель
 или
Математическая модель
. (здесь
Математическая модель
, т.к. второе осреднение по условию не меняет результата). Так как левая часть уравнения равна
Математическая модель
, то
Математическая модель
. По аналогии
Математическая модель
;
Математическая модель
. Следовательно, среднее значение пульсационных составляющих равно нулю. (Но надо учесть, что
Математическая модель
;
Математическая модель
 и т.д.). Применяя вышесказанное к исходной системе уравнений (2.1), можно после определенных преобразований получить уравнения турбулентного пограничного слоя в следующем виде:

Математическая модель
 (2.2)

Здесь а)
Математическая модель
 б)
Математическая модель
в)
Математическая модель


г)
Математическая модель
, где
Математическая модель
, д)
Математическая модель
,

 где
Математическая модель
.

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

Для вывода уравнений турбулентного пограничного слоя надо осреднить исходные уравнения погранслоя, несколько преобразовав первое уравнение - уравнение движения (аналогично случаю ламинарного пограничного слоя).



Для этого уравнение неразрывности умножим на ux

Математическая модель


 и добавим его в левую часть первого уравнения системы (2.1)

Математическая модель
.

В результате преобразований (как и в случае ламинарного погранслоя - уравнение (1.21)) первое уравнение системы (2.1) получим в виде:

Математическая модель
.

Здесь:
Математическая модель
;
Математическая модель
.

Проведем над обеими частями этого равенства операцию осреднения:

Математическая модель
  (2.3)

Математическая модель
 (для первого члена используется правило осреднения
Математическая модель
). Так как
Математическая модель
, то
Математическая модель
;

Математическая модель
.

Так как
Математическая модель
, то
Математическая модель
. Аналогично

Математическая модель
 и
Математическая модель


Подставляя значения
Математическая модель
и
Математическая модель
в уравнение (2.3), получим:

Математическая модель
 (2.4)

Учитывая уравнение неразрывности в осредненном виде:

Математическая модель
,  (2.5)

можно уравнение движения (2.4) записать так:

Математическая модель
  (2.6)

С этой целью левая часть уравнения (2.4) преобразовывается с учетом уравнения неразрывности следующим образом:

Математическая модель
.

Уравнения (2.5) и (2.6) входят в систему дифференциальных уравнений Рейнольдса осредненного турбулентного движения несжимаемой вязкой жидкости, которую можно окончательно представить в виде:

Математическая модель
 (2.7)

Эта система имеет одинаковый вид как для основного течения жидкости, так и для течения жидкости в погранслое.

Сопоставим первое уравнение системы (2.7) с уравнением движения вязкой жидкости в напряжениях, которое выглядит следующим образом:

Математическая модель
.

В случае одномерного стационарного движения и отсутствия массовых сил это уравнение имеет вид:

Математическая модель
.                         (2.8)

Сравнивая уравнение Рейнольдса с уравнением движения в напряжениях, можно представить себе правую часть уравнения Рейнольдса как результат подстановки в уравнение в напряжениях  вместо величин pxx и pxy суммы вязких напряжений, определяемых обобщенным законом Ньютона, и дополнительных турбулентных напряжений p'xx и p'xy, возникших за счет наличия в потоке пульсаций, т.е.:

Математическая модель
Математическая модель
.

В нашем случае

а)
Математическая модель
(т.к. в уравнении Прандтля пограничного слоя пропадает член
Математическая модель
при стремлении Re¥®¥).

Математическая модель
. Тогда
Математическая модель
.

б)
Математическая модель
 (т.к. в уравнении Прандтля пограничного слоя пропадает член
Математическая модель
 при стремлении Re¥®¥).
Математическая модель
. Тогда
Математическая модель
.

Таким образом получаем полную тождественность уравнений движения в напряжениях (2.8) и Рейнольдса (2.7).



В общем случае трехмерного движения эти дополнительные турбулентные напряжения p'xx, p'xy и т.д. образуют, так же как и вязкие напряжения, симметричный тензор второго ранга:

Математическая модель
 (2.9)

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

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

 Назовем тензором полного (суммарного) напряжения тензор P, равный

 
Математическая модель
,  (2.10)

и имеющий компоненты:

Математическая модель
 (2.11)

Не только вид уравнений движения, но и вид уравнений импульсов (интегральное соотношение Кармана) в турбулентном пограничном слое остается таким же, как и для ламинарного пограничного слоя:

Математическая модель
,  (2.12)

только значения d, d*, d** и tw (напряжение трения на твердой стенке) будут иными:

а) толщина вытеснения масс в пограничном турбулентном слое

Математическая модель
;  (2.13)

б) толщина потери импульса в турбулентном погранслое

Математическая модель
; (2.14)

в) напряжение трения на твердой стенке
Математическая модель
.

Граничные условия будут следующими :

а) на стенке:
Математическая модель


б) на внешней границе турбулентного погранслоя:

 
Математическая модель
,

Необходимо учесть, что уравнение Эйлера
Математическая модель
в случае обтекания плоской пластины преобразуется к виду: u'¥=0 (т.к. в этом случае u¥=ux,¥ постоянна вдоль оси Х и тогдаu'¥=u'x,¥=0 - нет изменения скорости вдоль пластины) и
Математическая модель
. Для плоской пластины уравнение импульсов имеет вид:

Математическая модель
 (2.16)


Содержание раздела