Praca doktorska - Instytut Geofizyki

Praca doktorska - Instytut Geofizyki Praca doktorska - Instytut Geofizyki

igf.fuw.edu.pl
from igf.fuw.edu.pl More from this publisher
23.01.2014 Views

Prędkość powietrza na wlocie 70 m/s 0.25 0.2 Ciśnienie odsysania stosowane w pomiarach [62 hPa] (d) ( 1[ m] przedział ) -1 0.15 0.1 Ciśnienie odsysania po optymalizacji [30 hPa] 0.05 0 5 10 15 20 25 Średnica kropli [ m] Prędkość powietrza na wlocie 90 m/s 0.25 0.2 Ciśnienie odsysania stosowane w pomiarach [117 hPa] (d) ( 1[ m] przedział ) -1 0.15 0.1 Ciśnienie odsysania po optymalizacji [50 hPa] 0.05 0 5 10 15 20 25 Średnica kropli [ m] Rys. 5.10 Funkcja gęstości prawdopodobieństwa Φ(2r) policzona dla różnych rozmiarów kropel, dwóch różnych prędkości 70 i 90 m/s i dla wartości ciśnień użytych w eksperymentach (62 i 117 hPa) oraz dla optymalnych (30 i 50 hPa) znalezionych na drodze symulacji. Wyniki symulacji uzyskane dotychczas są poprawne tylko dla dwóch wymiarów i nie uwzględniają zależności od czasu. Funkcja (r) zależy tylko od promienia. Dla uogólnienia wniosków na trzeci wymiar oraz wyznaczenia częstotliwości zachodzenia kolizji konieczna jest znajomość postaci widma kropel chmurowych. Dalsze obliczenia będą prowadzone w oparciu o widmo kropel przedstawione na rysunku 5.1. Pierwszym krokiem w opisie ilościowym zjawiska jest znalezienie całkowitego czasu T, po którym 10 000 kropel chmurowych przejdzie przez prostokąt o polu wyznaczonym przez długość drutu i szerokość wlotu kropel, dla którego wyznaczono funkcję (r) czyli 0,2 x 0,5 cm 2 . Czas ten dany jest wzorem 56 1 1 1 1 T U l C 10 000 0.2 [s], (5.4) gdzie U [cm/s] jest prędkością napływu powietrza (prędkość samolotu), l = 0,5 cm długość drutu oporowego, C to koncentracja kropel na jednostkę objętości, a to 0,2 cm szerokość wlotu kropel. Kolejny etap to wyprowadzenie wzoru, który będzie określał, ile spośród tych kropel, które przeleciały przez prostokąt w czasie T, trafi w czujnik temperatury. Dane o koncentracji kropel podane w (Warner, 1969 z Pruppachera i Kletta 1978, § 2) są dyskretne. Każda wartość odpowiada przedziałowi, w którym promień zawiera się pomiędzy r i a r i dr , gdzie dr jest

szerokością przedziału i wynosi 1 ,25 m . Wewnątrz każdego przedziału ( ri , ri dr) koncentracja kropel jest stała. Interesujące wartości wielkości promienia, dla których zarówno koncentracja, jak i funkcja (r) są większe od zera to r1 2, 5 m ; r2 3, 75 m ; r3 5 m ; r4 6, 25 m . Koncentracja kropel w poszczególnych przedziałach ( ri , ri dr) wynosi: 3 C 1 r1 , r1 dr 100cm , 3 C 2 r2 , r2 dr 40cm , 3 C 3 r3 , r3 dr 40cm , 3 C 4 r4 , r4 dr 100cm . Dla każdego z przedziałów promieni liczba kolizji n(r i , r i + dr) została wyznaczona na podstawie funkcji przedstawionej na wykresie 5.10 n r i , r dr * * d 2r , (5.5) i i i gdzie n(r i , r i + dr) to liczba kropel o średnim promieniu r i trafiających w drut o średnicy d = 2,5 m, a i ( ri , ri dr) to średnia wartość prawdopodobieństwa w przedziale ( r i , r i dr) . Ostatecznie liczba kropel o promieniu z przedziału ( r i , r i dr) trafiających w drut na jednostkę czasu dana jest wzorem n ri , ri dr 5 1 N i U l Ci ( d 2ri ) i 2 10 [s ]. (5.6) T Całkowita liczba kropel trafiająca w drut na jednostkę czasu wynosi 5 N 2 10 U l C ( d 2r ). (5.7) c i i i 5 Stała 2 10 wynika z przyjętej definicji funkcji (r) . Można też zdefiniować funkcję gęstości prawdopodobieństwa w taki sposób, aby uniknąć tej stałej, jednak przyjęta w tej formie definicja funkcja (r) jest łatwo interpretowalna. N c Dla 110[s 1 i i odpowiadającej prędkości U = 70 m/s i wartości odsysania 62 hPa otrzymujemy ]. Dla tej samej prędkości i zoptymalizowanego ciśnienia 30hPa częstotliwość zderzeń jest c 1 mniejsza prawie o czynnik 2 i wynosi N 60[s ]. Dla większych prędkości U 90 m/s i ciśnienia odsy- 1 sania użytego w eksperymentach 117 hPa otrzymujemy N c 130[s ]. Redukcja ciśnienia do minimalnej wartości 50 hPa przy prędkości U= 90 m/s ogranicza liczbę kolizji do wartości nieco mniejszej niż w poprzednim przypadku. N 80[s ] 1 . c Objętości wody, która osiada na drucie po zderzeniu z kroplą, nie można oszacować na podstawie liczby zderzeń i znajomości rozmiarów kropel. Krople mają większe średnice niż drut i podczas zderzeń część wody ulega rozproszeniu. Poza tym, ilość wody osiadającej na drucie zależy ściśle od parametru zderzenia. Nawet przy centralnym zderzeniu małe krople o średnicy 5 μm nie osiadają w całości. Możliwa jest również taka sytuacja, w której krople są przecinane na pół, a następnie odrywają się w całości i tworzą się z nich dwie nowe, mniejsze krople. Wówczas istotnym parametrem jest czas w którym kropla styka się z drutem. Czas ten jest zależy od prędkości samolotu oraz od parametru zderzenia. Zgrubne oszacowania dają wynik w granicach ułamków mikro, sekund czyli poniżej czasu inercji termicznej. Dodatkowo podczas zderzeń istotne jest działanie napięcia powierzchniowego oraz prędkości, dlatego w celu dokładniejszego przeanalizowania problemu konieczne jest przeprowadzenie modelowania w mikro skali, które pozwoli na uwzględnienie chwilowego kształtu kropel. W przeprowadzonych symulacjach nie uwzględniono wirowości przepływu. W pobliżu warstwy przyściennej krople mogą doznawać rotacji, która również finalnie może zmieniać obraz zderzenia kropli z drutem. 57

szerokością przedziału i wynosi 1 ,25 m . Wewnątrz każdego przedziału ( ri , ri<br />

dr)<br />

koncentracja kropel jest<br />

stała. Interesujące wartości wielkości promienia, dla których zarówno koncentracja, jak i funkcja (r)<br />

są<br />

większe od zera to r1 2, 5 m ; r2 3, 75 m ; r3 5 m ; r4 6, 25 m . Koncentracja kropel w poszczególnych<br />

przedziałach ( ri , ri<br />

dr)<br />

wynosi:<br />

3<br />

C<br />

1<br />

r1<br />

, r1<br />

dr 100cm<br />

,<br />

3<br />

C<br />

2<br />

r2<br />

, r2<br />

dr 40cm<br />

,<br />

3<br />

C<br />

3<br />

r3<br />

, r3<br />

dr 40cm<br />

,<br />

3<br />

C<br />

4<br />

r4<br />

, r4<br />

dr 100cm<br />

.<br />

Dla każdego z przedziałów promieni liczba kolizji n(r i , r i + dr) została wyznaczona na podstawie<br />

funkcji przedstawionej na wykresie 5.10<br />

n r i<br />

, r dr * * d 2r<br />

, (5.5)<br />

i<br />

i<br />

i<br />

gdzie n(r i , r i + dr) to liczba kropel o średnim promieniu r<br />

i<br />

trafiających w drut o średnicy d = 2,5 m, a<br />

i<br />

( ri<br />

, ri<br />

dr)<br />

to średnia wartość prawdopodobieństwa w przedziale ( r i<br />

, r i<br />

dr)<br />

.<br />

Ostatecznie liczba kropel o promieniu z przedziału ( r i<br />

, r i<br />

dr)<br />

trafiających w drut na jednostkę czasu<br />

dana jest wzorem<br />

n ri<br />

, ri<br />

dr<br />

5 1<br />

N<br />

i<br />

U l Ci<br />

( d 2ri<br />

)<br />

i<br />

2 10 [s ].<br />

(5.6)<br />

T<br />

Całkowita liczba kropel trafiająca w drut na jednostkę czasu wynosi<br />

5<br />

N 2 10 U l C ( d 2r<br />

).<br />

(5.7)<br />

c<br />

i<br />

i<br />

i<br />

5<br />

Stała 2 10 wynika z przyjętej definicji funkcji (r)<br />

. Można też zdefiniować funkcję gęstości prawdopodobieństwa<br />

w taki sposób, aby uniknąć tej stałej, jednak przyjęta w tej formie definicja funkcja (r)<br />

jest<br />

łatwo interpretowalna.<br />

N<br />

c<br />

Dla<br />

110[s<br />

1<br />

i<br />

i<br />

odpowiadającej prędkości U = 70 m/s i wartości odsysania 62 hPa otrzymujemy<br />

]. Dla tej samej prędkości i zoptymalizowanego ciśnienia 30hPa częstotliwość zderzeń jest<br />

c<br />

1<br />

mniejsza prawie o czynnik 2 i wynosi N 60[s ]. Dla większych prędkości U 90 m/s i ciśnienia odsy-<br />

1<br />

sania użytego w eksperymentach 117 hPa otrzymujemy N<br />

c<br />

130[s ]. Redukcja ciśnienia do minimalnej<br />

wartości 50 hPa przy prędkości U= 90 m/s ogranicza liczbę kolizji do wartości nieco mniejszej niż w poprzednim<br />

przypadku. N 80[s ]<br />

1<br />

.<br />

c<br />

Objętości wody, która osiada na drucie po zderzeniu z kroplą, nie można oszacować na podstawie<br />

liczby zderzeń i znajomości rozmiarów kropel. Krople mają większe średnice niż drut i podczas zderzeń<br />

część wody ulega rozproszeniu. Poza tym, ilość wody osiadającej na drucie zależy ściśle od parametru zderzenia.<br />

Nawet przy centralnym zderzeniu małe krople o średnicy 5 μm nie osiadają w całości. Możliwa jest<br />

również taka sytuacja, w której krople są przecinane na pół, a następnie odrywają się w całości i tworzą się z<br />

nich dwie nowe, mniejsze krople. Wówczas istotnym parametrem jest czas w którym kropla styka się z drutem.<br />

Czas ten jest zależy od prędkości samolotu oraz od parametru zderzenia. Zgrubne oszacowania dają<br />

wynik w granicach ułamków mikro, sekund czyli poniżej czasu inercji termicznej. Dodatkowo podczas zderzeń<br />

istotne jest działanie napięcia powierzchniowego oraz prędkości, dlatego w celu dokładniejszego przeanalizowania<br />

problemu konieczne jest przeprowadzenie modelowania w mikro skali, które pozwoli na<br />

uwzględnienie chwilowego kształtu kropel. W przeprowadzonych symulacjach nie uwzględniono wirowości<br />

przepływu. W pobliżu warstwy przyściennej krople mogą doznawać rotacji, która również finalnie może<br />

zmieniać obraz zderzenia kropli z drutem.<br />

57

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!