О ТОЧНЫХ РЕШЕНИЯХ κ-ε - МОДЕЛИ ТУРБУЛЕНТНОСТИ

УДК 514.7

О ТОЧНЫХ РЕШЕНИЯХ к —£-МОДЕЛИ ТУРБУЛЕНТНОСТИ

Н.Г. ХОРЬКОВА

В статье метод поиска инвариантных решений применяется к системе уравнений, описывающих к — £ -модель турбулентности. Найдены точные решения этой системы. Методика построения точных решений, использованная в данной работе, применима и к другим моделям турбулентности.

Ключевые слова: системы дифференциальных уравнений в частных производных, локальные симметрии, инвариантные решения, к — £ -модель турбулентности.

1. ВВЕДЕНИЕ

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

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

В данной работе мы строим многопараметрические семейства инвариантных решений диссипативной двухпараметрической модели турбулентности (к — £ -модель).

2. ОСНОВНЫЕ ПОНЯТИЯ ТЕОРИИ СИММЕТРИЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ

Пусть Е = {¥ = 0} — система дифференциальных уравнений в частных производных (далее "дифференциальное уравнение" или просто "уравнение") в векторном расслоении п: Еп+т ^ Мп. В рамках алгебро-геометрической теории любое дифференциальное уравнение рассматривается как подмногообразие пространства джетов к -го порядка Jk (п) расслоения п, где к — максимальный порядок уравнений, входящих в систему, п — число независимых переменных, а т — число неизвестных функций (зависимых переменных). На многообразии бесконечных джетов J " (п) имеется п -мерное интегрируемое распределение С (распределение Картана), задаваемое операторами полных производных Д,..., Д, которое допускает ограничение С на бесконечно

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

В алгебро-геометрической теории высшая инфинитезимальная симметрия отождествляется с эволюционным дифференцированием

э,=Е Д—

которое однозначно определяется своей производящей функцией ( = ((,..., (т), ( е С- (J- (п)). В формуле (1), определяющей эволюционные дифференцирования в координатной форме, использованы следующие обозначения: х1У...,хп, и1,.,ит, и— (и^ = и1) — канонические координаты в пространстве бесконечных джетов J~ (п), Д. = + ^ . и— -Ц-

оператор полной производной по переменной хг, Да — композиция операторов полных производных, соответствующая мультииндексу —, черта обозначает ограничение оператора на бесконечно продолженное уравнение Е ~ (подробнее см., например, [1], [6]).

Для вычисления алгебры симметрий уравнения Е : ¥ = 0 надо решить систему

Те

где 1¥ — оператор универсальной линеаризации, отвечающий функции ¥ = (¥р...,¥г), ¥ е С - (Jk (п)), те(() = Э((¥).

Множество всех высших инфинитезимальных симметрий данного уравнения является алгеброй Ли относительно скобки Якоби {(,\} = Э((\) - ЭДф) . Компоненты скобки Якоби

вычисляются по формуле:

(¥у = ДО \ - ДАГ)(^

—, а

V ди а ди— J

Решения уравнения Е = {¥ = 0}, инвариантные относительно симметрии (, определяются следующей системой дифференциальных уравнений:

Г ¥1 = 0,., ¥г = 0,

I ( = 0.

Решение уравнения Е = {¥ = 0}, инвариантное относительно симметрий ( и \, будет

также инвариантно относительно скобки {(\}. Поэтому следует искать решения уравнения,

инвариантные относительно некоторой подалгебры g алгебры симметрий. Если (, ... ( —

образующие этой подалгебры, отыскание инвариантных решений сводится к решению переопределенной системы

Г ¥1 = 0,., ¥г = 0, 1(1 = 0—( = 0

3. к -£ -МОДЕЛЬ ТУРБУЛЕНТНОСТИ

Для расчета турбулентных течений с использованием к — £ -модели турбулентности система осредненных по Рейнольдсу уравнений Навье-Стокса замыкается двумя уравнениями для турбулентной кинетической энергии к и скорости диссипации турбулентной энергии £. В настоящее время используются различные варианты к — £ -уравнений. В данной работе рассматривается следующая система, которая описывает движение несжимаемой жидкости или газа с учетом турбулентности при больших числах Рейнольдса ([4], [3]):

д-

дх

д(Р-1) , д(р-1-1) _ др

_ 0;

дг

-+

дхг

= -дГ+д7 Ту-рп>п11' 1_1'2'3;

дк+- д^ __д

дг ' дхг дхг

У + -

С„ к

г IV

д- 1

&к £ у-,

дк I

-> - — -. -;

дхг 1 дхг

д£ + - д£ дг ' дхг дхг

Г (

у +

с^ к

2 Л

г IV

сг£ £

д£ У дхг

д- .

с — - , -

£ к г 1 дхг

где

Т _РУ

(д_ , д_

дх. + дх,г

V 1 г У

- -. - . _

г 1

с к

2

+

д-

_ Л

дх. дх гг

V 1 1 J

- ,

3 г

у,р,с1и,с£,с£ ,ак,ае-константы.

Это система на шесть неизвестных функций - _ - -2 _ - 2, - _ - 3, -4 _ р, -5 _ к,

- _ £, зависящих от четырех переменных х1

х3, х4 _ г. При описании классических

симметрий используются канонические координаты в пространстве джетов [1]. 4. ИНВАРИАНТНЫЕ РЕШЕНИЯ

Алгебра классических симметрий к -£ -модели турбулентности (3) - (6) была вычислена в [3]. В данной работе мы ищем решения, инвариантные относительно подалгебры, порожденной симметриями Х1(/), Х2(§), Х3(И), где

X (/) _ (/-1 - /, /-2, /-, М + рх, , /£), Х2 (§) _ (§-2 > §-2 - §> §-\, §Р2 + Р§х2 > §к2 , §£2 X

Х3 (И) _ (И-, к-1, к-1 - И., Нр3 + рИх3, кк3, к£3),

х

2

/, §, И

— некоторые функции переменной г.

Рассматриваемые симметрии при подходящем выборе функций I, g, Н включают в себя трансляции по переменным хг и галилееву симметрию. Симметрии Х1(/), Х2(g), Хз(Н) попарно коммутируют:

X (I), Х2 (g)} = {X2 (g), Хз (Н) } = Х1 (I), Хз( Н) } = 0.

Покажем на простом примере, как можно с помощью симметрий свести исходную систему к системе уравнений с меньшим числом независимых переменных. Для отыскания решений системы (3) - (6), инвариантных, например, относительно симметрии Х1(1), сначала

решим систему уравнений Х1 (I) = 0 :

д-1 - I

I% -1 = 0 ^ и1 = 1x1 + А($, х2, Х3);

д—2

I^ = 0 ^ и2 = ВЦ,Х2,Хз);

• дц"3

1= 0 ^ и' = С^их2,Хз); Iдр = 0 ^ Р = -р1 х2 +ДЦ,х2,хз);

Iддк = 0 ^ к = Е«,х2,хз);

I= 0 ^ £ = ¥Ц,х2,хз).

Затем полученные функции и1, и2, и3, р, к, £ подставим в систему (з) - (6), в результате чего данная система преобразуется в систему уравнений на шесть функций трех независимых переменных t, х2, хз.

Легко видеть, что функции и1, и2, и3, р , к, £, удовлетворяющие системе уравнений

Х1С/) = Х2 (g) = Хз( Н) = 0,

имеют следующий вид:

/ g Н

и1 = х + а

I g Н

г - \

Р

Р = ~ 4х2 + ^х2 +Нхз2 + р

2 [I g Н )

к = к Ц ),£ = £

После подстановки (7) - (8) в уравнения (з) - (4) получаем, что

а = о^/"1, в = а.2g-1, у = аъНт\ М = А,

где А, аг, — произвольные константы, и

u

1 fxCj I 2 I ^2 3 I

f

g

/

— du

Вычислим сумму —u uj -g^- для функций u, определяемых формулой (ii):

"^du;

f ck2 '

—u. u.

1 1 dx.

_ О M

i, j=1

к2 з (Л =i Idx 1 J

du. + du.

dx, dx ,

V j 1 J

—

3 1

du,

J

dx.

-2 к ± *=к1 l,

3 7=i dx. £

где

L = 2c,.

f Л2 (gЛ2 (h^

f

+

\J J

g_

g J

+

V h J

Теперь система уравнений (5) - (6) с учетом (9), (12) принимает вид:

f т к2 к = L--£;

£

£ = c£LA — c£ —

Рассмотрим сначала случай L = 0, то есть f = g = h = 0 и u = const (в этом случае симметрии Xi(f), X2(g), X3(h) являются трансляциями):

£ к

£

:£2 к

2

Исключим £ и получим уравнение для к:

к к" к = c£ —г

£2 к2

или

откуда получаем

к к -

— = cr (ln к) — cr (ln к) = 0, к £ к £

к к

const -= Cit + C 2 (ce, * 1).

к2 i — c

Итак,

1 с,

1 - с . C 1 - с

k = (C1t + C2) £, £ = -k = —^(C1t + C2) £

1 - с„

где С1, С2 — произвольные постоянные.

Рассмотрим теперь общий случай и преобразуем теперь систему (1з) - (14) к одному уравнению. Умножив уравнения системы на £ и к соответственно, вычислим производную

k \ £k - k£ (1 - cJLk2 + (с -1)£> Г k4 2

= (1 - с£LI — I + (С£2 -1).

V £) £ £ " V £.

Введем новую функцию w = k. Тогда

£

w = -l

l (t) = (с£ -1) L > 0, a = с£-1 > 0.

Подставив выражение k = w£ в систему (13) - (14), убеждаемся, что она сводится к одному уравнению

£ =

из которого по решению w уравнения (15) находим функцию £.

1 v

Заметим, что уравнение (15) с помощью подстановки w =— сводится к линейному

l v

уравнению на функцию v :

lv - iv = al2v.

Таким образом, поиск инвариантных решений для рассматриваемой подалгебры проводится по следующей схеме. Сначала выбираем функции f, g, h, удовлетворяющие

условию (10), и по формулам (11), (8) находим и' и p. Затем вычисляем функцию l(t) и,

1 v

подставив ее в уравнение (17), находим его решение v. Затем вычисляем w = —, подставляем

l v

w в (16) и находим функцию £, а затем функцию k по формуле k = £w.

Например, для функций (выбраны с учетом (10)) f = b1eM, g = b2e^, h = b3e-(Ä+M)t, где X, ¡1, bi — произвольные постоянные, получаем:

и1 = Xx + — e~Xt, и2 = ¡x2 + ^ e~1t, и3 = -(X + ¡)Xj + ^ e(Ä+1)t, p = p(t).

b1 b2 b3

Функция l = const, l = 0 , и общее решение уравнения (i7) имеет следующий вид

V = Qe^ + C4 e~4alt

(C3, C4 — произвольные постоянные). Подставляем функцию

w

a C3el4alt - C4

l C,el4alt + C4

в уравнение (i6) и после интегрирования получаем выражения для £:

. _Í£L_

Ce^ + CAe~4alt'^

* = C5'

C e^ait L c 1 cei-1

a fQe^ + C4e-^lt ЛСЛ-1

1 íC.e^ - C4e~^lt

c~-1

и для к:

к _ С5

( С5 — произвольная постоянная). 5. ЗАКЛЮЧЕНИЕ

В настоящий момент имеется большое количество моделей для расчёта турбулентных течений. Для проверки пригодности той или иной модели для решения конкретной задачи иногда требуется проведение дорогостоящих экспериментов. Наличие инвариантных решений может помочь сделать выбор адекватной модели [6]. Методика построения точных решений, использованная в данной работе, применима и к другим моделям турбулентности.

ЛИТЕРАТУРА

1. Бочаров А.В., Вербовецкий A.M., Виноградов А.М., Дужин С.В., Красильщик И.С., Торхов Ю.Н., Самохии А.В., Хорькова Н.Г., Четвериков В.Н. Симметрии и законы сохранения уравнений математической физики. 2-е изд. М.: Изд-во Факториал-Пресс, 2005.

2. Khor'kova N.G., Verbovetsky A.M On symmetry subalgebras and conservation laws for k-e turbulence model and the Navier-Stokes equation. Amer. Math. Soc. Transl. Series 2. i995. V. i67. P. 6i-90.

3. Jones W.P., Launder B.E. The prediction of laminarization with a two-equation model of turbulence. Internat. J. Heat Mass Transfer i5 (i972), no. 2, p. 30i-3i4.

4. Kollmann W. (ed.) Prediction method for turbulent flows, Hemisphere, Washington, i980, 3i2 p.

5. Oberlack M. Symmetries and Invariant Solutions of Turbulent Flows and their Implications for Turbulence Modelling. — Theories of Turbulence. International Centre for Mechanical Sciences. Vol. 442. 2002. P. 30i-366.

6. Вербовецкий A.M., Хорькова Н.Г., Четвериков В.Н. Симметрии дифференциальных уравнений: учеб. пособие.- М.: Изд-во МГТУ им. Н. Э. Баумана, 2002.

ON EXACT SOLUTIONS OF k — £ TURBULENCE MODEL

Khorkova N.G.

A method of constructing invariant solutions is applied to k — £ turbulence model. Exact solutions for k — £ turbulence model are obtained. The method for construction of exact solutions used in this paper may be applied to other models of turbulence.

Keywords: nonlinear differential equation, local infinitesimal symmetries, invariant solution, k — £ turbulence model.

REFERENCES

1. Bocharov A.V., Chetverikov V.N., Duzhin S.V., Khor'kova N.G., Krasil'shchik I.S., Samokhin A.V., Torkhov Yu.N., Verbovetsky A.M., Vinogradov A.M. Symmetries and Conservation Laws for Differential Equation of Mathemetical Phisics // Translations of Mathematical Monographs, vol. 182 - Providence, RI: AMS, 1999. 333 p.

2. Khor'kova N.G., Verbovetsky A.M. On symmetry subalgebras and conservation laws for k-e turbulence model and the Navier-Stokes equation. Amer. Math. Soc. Transl. Series 2, 1995, v. 167, pp. 61-90.

3. Jones W.P., Launder B.E. The prediction of laminarization with a two-equation model of turbulence. Internat. J. Heat Mass Transfer, 15 (1972), no. 2, pp. 301-314.

4. Kollmann W. (ed.) Prediction method for turbulent flows, Hemisphere, Washington, 1980, 312 p.

5. Oberlack M. Symmetries and Invariant Solutions of Turbulent Flows and their Implications for Turbulence Modelling. Theories of Turbulence. International Centre for Mechanical Sciences. — Volume 442, 2002, pp. 301-366.

6. Verbovetsky A.M., Khor'kova N.G., Chetverikov V.N. Simmetrii differenzialnych uravneniy. (Symmetries of Differential Equations). Uchebnoe posobiye. Moscow, MGTU, 2002, 36 p.

Сведения об авторе

Хорькова Нина Григорьевна, окончила МГУ им. М.В. Ломоносова (1983), кандидат физико-математических наук, доцент кафедры прикладной математики МГТУ им. Н. Э. Баумана, автор 17 научных работ, область научных интересов - алгебро-геометрическая теория дифференциальных уравнений.