Вестник Московского государственного областного университета. Серия Естественные науки 2018 2

• Название: Вестник Московского государственного областного университета. Серия Естественные науки 2018 2
• Анкор: SRHSE
• Дата: 22.08.2020
• Размер: 0.93 Mb.
• Имя файла: chislennoe-modelirovanie-vozniknoveniya-konvektivn-h-techeniy-pr.doc
• Тип: Документы; #135918
• страница: 2 из 3

136




ISSN 2072-8360 Вестник Московского государственного областного университета. Серия: Естественные науки 2018 / № 2

где: t – время, u – скорость текучей среды, ρ – плотность текучей среды, p – дав-ление текучей среды, S_i – внешние массовые силы, действующие на единичную массу текущей среды, в нашем случае S _i S_i^gravity , действие гравитации S _i^gravity g_i , g_i– составляющая гравитационного ускорения в координатном направлении x_i.
Для ньютоновских сред тензор вязких сдвиговых напряжений определяется как:

u

i

u j

2

u

k

2

(5)

ij

ij

k ij

,

x

x

3

x

3

j

i

k

где: μ = μ_l + μ_t , μ_l – коэффициент динамической вязкости, μ_t – коэффициент турбулентной вязкости, δ_ij – дельта функция Кронекера (δ_ij = 1 при i = j; δ_ij = 0 при i≠ j), k – кинетическая энергия турбулентности. В соответствии с k – ε модельютурбулентности, μ определяется через величины кинетической энергии турбу-лентности k и диссипации этой энергии ε:

f

c k 2

,

(6)

где f

1 exp 0, 025 Ry 2

1

20,5

, Ry

k y

, RT

k 2

, y – расстояние от поверхно-

l

RT

l

сти стенки, c

0,09 .

Кинетическая энергия турбулентности k и диссипация этой энергии ε опреде-ляются из следующих уравнений:

k

u k k

l

t

k

Sk

,

(7)

t

x

x

x

k

k

k

k

u k

l

t

k

S ,

(8)

t

x

x

x

k

k

k

где S k

R

u

u

2

ij

i

l PB

, S

c 1

f i ijR

i

i c B PB

c 2 f

,

xi

2

k

x

k

j

u

u j

2

u

2

, PB

g

i

1

ijRt

i

ij

l

k ij

,

x

x

x

3

x

3

j

B

i

l

i

3

σ

= 0,9, c

=1приP

0 и c

=0при P 0,

0, 05

B

B

B

f11

,

B

B

f

f 21 exp RT2,c11, 44,c21,92,

1,3,

k

1.

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

		l			t		T		k = 1, 2, 3	(9)
qk						cp		,
		Pr
				c			xk