Računanje turbulentnih flukseva količine kretanja i toplote primenom teorije Monin-Obuhov | seminarski diplomski
Ovo je pregled DELA TEKSTA rada na temu "Računanje turbulentnih flukseva količine kretanja i toplote primenom teorije Monin-Obuhov". Rad ima 6 strana. Ovde je prikazano oko 500 reči izdvojenih iz rada.
Napomena: Rad koji dobjate na e-mail ne izgleda ovako, ovo je samo DEO TEKSTA izvučen iz rada, da bi se video stil pisanja. Radovi koje dobijate na e-mail su uređeni (formatirani) po svim standardima. U tekstu ispod su namerno izostavljeni pojedini segmenti.
Uputstvo o načinu preuzimanja rada možete pročitati OVDE.
RAČUNANJE TURBULENTNIH FLUKSEVA KOLIČINE KRETANJA I TOPLOTE PRIMENOM TEORIJE MONIN-OBUHOV
1 1
UVOD
U ovom radu biće prikazana jedna od mogućih parametrizacija turbulentnih procesa u površinskom sloju atmosfere, koja se oslanja na teoriju sličnosti M onin-Obuhov (1). Naime od interesa nam je da primenom ove teorije odredimo turbulentne flukseve količine kretanja i toplote. Dobra strana ove teorije je da se na prilično jednostavan način dolazi do turbulentnih karakteristika prizemnog sloja atmosfere, odnosno fluksevi se odredjuju samo na osnovu razlike tempeature i vetra, na površini i nekoj zadatoj visini. I pored svoje jednostavnosti, odredjivanje fluskeva na ovaj način pokazalo se prilično pouzdanim u većini slučajeva. Nedostatak teorije je da vrednosti dobijenih flukseva mogu biti potpuno netačni u odredjenim ekstremnim situacijama kao sto je slobodna konvekcija (2).
MATERIJAL I METODE
∂Θ θ = − * Φ h (ζ ) ∂z kz
(2)
gde je u * brzina trenja, θ* = wθ / u * razmer temperature, odnosno wθ turbulentni fluks toplote, k=0.4 Fon-Karmanova konstanta, ζ= z L
L je tzv. dužina M onin-Obuhova data izrazom u *3 L= kgβ wθ (3)
konstante koje se pojavljuju u (3) imaju vrednosti g=9.81, β =1/300. Funkcije Φm i Φ h su empiriski odredjene, merenjima na odgovarajućim lokacijama i to za slučaj nestabilna stratifikacija Φ m (ζ ) = (1 + bζ )-1/3 ; stabilna stratifikacija Φ m (ζ ) = (1 − 5ζ ) ; Φ h (ζ ) = (1 − 5ζ ) (4.2) Φ h (ζ ) = (1 + dζ ) -1/3 (4.1)
gde je Prt=0.9 Prantlov turbulentni broj a konstante b i d imaju vrednosti b=3A, A=5.5, d=3B, B=6.2. Da bi smo uz pomoć j-na (1) i (2) izračunali flukseve u površinskom sloju prvo ih moramo integraliti i to od z 0 , visine hrapavosti, do z, visine na kojoj jos uvek važi teorija M onin-Obuhov. Konkretno u programu je korišćena vrednost z=10m. Vrednost visine trenja biće konstantna u slučaju j-ne (1) i to z 0m =0.1m, dok u slučaju j-ne (2) vrenost z 0 je data izrazom z 0h 3Prt 0.067 2 .25 Pr 3 gβ wθ 1/ 4 = , z 0r = ( )( ) ( ) ku * 1 k 3Bk αt + αt z 0r 1
α t = 2.15 10 -6 , Pr = 0.609 , Prantlov broj, a ostale konstante uzimaju predhodne vrednosti.
Dakle integraljenjem (1) za slučaj nestabilne stratifikacije dobijamo
k z 3 x2 + x + 1 2x + 1 2x 0 + 1 U( z) = ln − ln 2 + 3[arctg( ) − arctg( )] u* z 0m 2 x 0 + x 0 + 1 3 3 Ako odavde izrazimo fluks količine kretanja u* = kU(z ) / {ln z 3 x2 + x + 1 2x + 1 2x 0 + 1 − ln 2 + 3[ arctg( ) − arctg( )]} (5) z0 m 2 x 0 + x 0 + 1 3 3
gde je x = (1 + dζ )1 / 3 , x 0 = (1 + dζ 0 )1 / 3 , ζ 0 = z0 h / L . slično dobijamo i za fluks toplote wθ = − {k[Θ( z) − Θ( z 0 )]u* } / {Pr t[ln z 3 y2 + y+ 1 2y + 1 2y + 1 − ln 2 + 3[arctg ( ) − arctg( 0 )]]} (6) z 0 h 2 y0 + y 0 + 1 3 3
...
CEO RAD MOŽETE PREUZETI NA SAJTU: WWW.MATURSKIRADOVI.NET