background image

 

 

 

 

W

W

Y

Y

K

K

Ł

Ł

A

A

D

D

 

 

1

1

4

4

 

 

R

R

U

U

C

C

H

H

 

 

T

T

U

U

R

R

B

B

U

U

L

L

E

E

N

N

T

T

N

N

Y

Y

 

 

 

 

C

C

.

.

D

D

.

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

“Gallery of Fluid Motion”-M. Samimy, K.S. Breuer

 

background image

R

R

Ó

Ó

W

W

N

N

A

A

N

N

I

I

E

E

 

 

R

R

E

E

Y

Y

N

N

O

O

L

L

D

D

S

S

A

A

 

 

Dowolna wielkość w opisie ruchu turbulentnego może być 
przedstawiona w postaci wolnozmiennej „średniej” i 
szybkozmiennej, niewielkiej oscylacji. Oscylacja (= pulsacja = 
fluktuacja) jest wie

lkością losową. 

 

 

 
 
 
 
 
Definicja wartości średniej 
 
 
 
 
 
 

srednia

oscylacja

f

f

f

 

czas 

średnia 

srednia

z oscylacji

f

0

 

 

 

t T

t T

1

f

f

d

2T

 

 

Czas uśredniania 2T powinien być 

większy od czasu charakteryzującego 

zmiany losowe i mniejszy od czasu, w 

którym istotnie zmienia się średnia 

background image

Pochodna 

f

 

jest funkcją miejsca i czasu.  

Pochodna średniej względem współrzędnej x ma postać:

 

 

 
 
 
 
 

Pochodna średniej względem położenia jest średnią pochodnej. 

 
Pochodną średniej względem czasu liczymy następująco: 

 

 

 

 

 

 

 

 

 

 

t T

t T

t T

t T

f

1

1

1

f

f ( ) d

f (t

T) f (t

T)

d

t

2T t

2T

2T

 

 

t T

t T

t T

t T

f

1

1

f ( )

f

f ( ) d

d

x

2T x

2T x

x

x

 

 

background image

Przedstawmy prędkość i ciśnienie w postaci sum: 

 

 

 

 

 
 
Podstawmy  powyższe  zależności  do  równań  Naviera  –  Stokesa 
dla cieczy: 
 
 
 

 

 

W rezultacie otrzymujemy równanie, które nosi nazwę 
 

R

ównania Reynoldsa dla średnich 

 
 
 
 
 

k

k

i

k

i

k

v

1

p

v v

v

t

x

x

 

 

 

i

k

i

k

k

i

k

i

k

v

1

p

v v

v

v

v

t

x

x

x

 

  

 

k

k

k

v

v

v

 

p

p

p

 

 

background image

T

T

E

E

N

N

S

S

O

O

R

R

 

 

R

R

E

E

Y

Y

N

N

O

O

L

L

D

D

S

S

A

A

 

 

 
Zauważmy, że równanie Reynoldsa różni się od równania Naviera- 

Stokesa dodatkowym członem 

i

k

i

( v v )

x

 

 . 

 
Oznaczmy: 
 
 
 
 

 

Tensor Reynoldsa

 

 

 
 
Średnie spełniają równanie ciągłości 

 

ii

k

ik

i

ik

ki

v v

v v

R

R

 

 

 

 

 

 

k

k

v

diw v

0

x

 

background image

 
Mamy układ opisujący średnie 

k

v

 i 

p

Trzeba określić składowe 

tensora Reynoldsa. 
 
 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

1

2

3

1

2

v

v

v , v , v ,

,

,

x

x

 

 

 

zależy od ruchu, a więc odpowiedni związek nie będzie 

wyrażał własności fizycznych płynu, lecz cechy ruchu!

 

Podanie powyższego  związku w jawnej formie nazywamy 

Hipotezą domknięcia

 

background image

Zapiszmy 

 

w sposób następujący: 

 
 
 
 
 

Tr(

)  

to ślad tensora Reynoldsa 

 

 

 
 

Gdzie  

 

nosi nazwę 

energii kinetycznej turbulencji 

 i jest to 

uśredniony kwadrat oscylacji turbulentnych. 
 
 

 

 
 

   

 

2

2

2

1

2

3

1 1

2

3

3

2

v

v

v

Tr( )

v v

v v

v v

2

2

2

 

 

 

 

 

 

ik

ik

ik

ik

1

1

R

Tr

R

Tr

3

3

 

Tr( )

2

 

 

ik

ik

ik

ik

2

2

R

R

3

3





 

 

background image

W równaniu Reynoldsa występują pochodne wielkości

i

k

( v v )

 

.

 

Zatem 
 
 
 
 
Wstawmy to wyrażenie do równania Reynoldsa. Dostaniemy 
wtedy:  
 
 
 
 
 
gdzie  
 
 
 
 
 

t

ik

i

ik

k

i

i

k

i

R

2

v v

T

x

x

3 x

x

 

 

 

t

k

ik

k

i

k

i

k

i

v

T

1

2

v

v

p

v

t

x

x

3

x



 

 

  

 

t

2

p

p

3



 

 

t

ik

ik

ik

2

T

R

3



 

ciśnienie turbulentne 

t

ensor naprężeń 

turbulentnych 

background image

Można pokazać, ze tensor naprężeń turbulentnych ma zerowy 
ślad. Taki ślad ma też tensor określony następująco: 
 
 
 
 
 
 
Hipotezą wykorzystującą powszechnie w opisach ruchu 
turbulentnego jest równanie: 
 
 
 
 
 

µ

turb

 

– nazywa się 

lepkością turbulentną

 

i zależy od rodzaju  

 

 

  

ruchu, miejsca i „zwykłej” lepkości .

 

 
 

i

k

sr

k

i

v

v

1

2

x

x

 

t

turb

sr

T

2

 

lub w składowych 

 

t

i

k

ik

turb

k

i

v

v

T

x

x

 

background image

Gdyby udało się określić 

µ

turb

, to równania dla średnich 

prędkości i ciśnienia byłyby takie: 
 
 
 
 
 
 
 
 
 
 
 
 

Jeśli znamy 

µ

turb

  

to  z powyższego układu równań możemy 

wyznaczyć 

k

turb

v

i

p

 

 
 

k

k

i

turb

k

turb

sr ik

i

k

i

v

1

v

v

p

v

t

x

x

x

 

  

 

i

i

v

0

x

 

równanie Reynoldsa 

równanie ciągłości dla średnich 

background image

Hipotezy określające 

µ

turb

  : 

 

1. 

 

 
 
 
 
2. 

Hipoteza „

κ - ε

” 

 
 

 

 

 

 

 

 

κ 

i

 

ε

 

wynikają z 

dwu dodatkowych równań różniczkowych cząstkowych

 

 

 

 

2

turb

v

l

n


 

v

n


 - 

 pochodna w kierunku normalnym  

 

    

dominującej składowej prędkości 

 

l

droga mieszania wyznaczana doświadczalnie

 

turb

turb

( , )

 

 

κ 

-  energia kinetyczna turbulencji

 

 

ε

 

– moc dyssypowana na ciepło skutkiem 

turbulencji