Praca doktorska - Instytut Geofizyki
Praca doktorska - Instytut Geofizyki Praca doktorska - Instytut Geofizyki
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
- Page 6 and 7: 1 Wstęp 1.1 Przedstawienie problem
- Page 8 and 9: Realizacja przedstawionych celów b
- Page 10 and 11: 2 Ultraszybki termometr chmurowy 2.
- Page 12 and 13: obserwowane zmiany oporu są spowod
- Page 14 and 15: Temperatura [K] du samolotu można
- Page 16 and 17: Fluktuacje temperatury [K] Fluktuac
- Page 18 and 19: 3 Modelowanie przepływu 3.1 Wprowa
- Page 20 and 21: gdzie K=t n+1 -t n , oznacza całko
- Page 22 and 23: 3.4.3 Konstrukcja siatki obliczenio
- Page 24 and 25: 4 Wyniki dwuwymiarowej symulacji dr
- Page 26 and 27: Ciśnienie [Pa] Rys. 4.3 Rozkład p
- Page 28 and 29: Ciśnienie [Pa] Ciśnienie [Pa] Pr
- Page 30 and 31: e ij 1 2 u x i j u x i j (4.4) 5 1
- Page 32 and 33: le, 2002, § 5). Dla czasu t 0,14μ
- Page 34 and 35: Fluktuacja temperatury [K] Fluktuac
- Page 36 and 37: Widmowa gęstość mocy [Pa 2 /Hz]
- Page 38 and 39: Częstość [kHz] Liczba Reynoldsa
- Page 40 and 41: Średnie ciśnienie [Pa] Średnie c
- Page 42 and 43: cona technicznemu udoskonaleniu ult
- Page 44 and 45: ozkładu sił działających na pow
- Page 46 and 47: - Niepożądane działanie odsysani
- Page 48 and 49: Prędkość powietrza na wlocie 90
- Page 50 and 51: Prawa statystyczne będące narzęd
- Page 52 and 53: Prędkość powietrza na wlocie 90
- Page 54 and 55: Prędkość na wlocie 70 [m/s] Prę
- Page 58 and 59: 5.7 Wnioski Odsysanie stabilizuje p
- Page 60 and 61: w przekroju ciśnienie statyczne je
- Page 62 and 63: wiry Rys. 6.4 Rozkład ciśnienia
- Page 64 and 65: Rys. 6.8 Linie prądu — widok z g
- Page 66 and 67: 7 Inne kształty osłon Przeprowadz
- Page 68 and 69: Ciśnienie [Pa] Funkcja prądu Rys.
- Page 70 and 71: 7.2 Osłona składająca się z trz
- Page 72 and 73: Rys. 7.8 Obrazy przedstawiają czas
- Page 74 and 75: danych producenta (H&N) przeznaczon
- Page 76 and 77: 8 Wizualizacja przepływu płynów
- Page 78 and 79: 9 Optyczne metody wizualizacji wir
- Page 80 and 81: gdzie K [m 3 /kg] jest stałą Glad
- Page 82 and 83: 1 x , y dx HT x [ x, y ] , (9.14) x
- Page 84 and 85: ną. Światło, które nie podlega
- Page 86 and 87: jedności. Zatem nieliniowa, przez
- Page 88 and 89: Ogniskowa soczewki [m] Rys. 9.6 Prz
- Page 90 and 91: Do wykonania symulacji użyto pokaz
- Page 92 and 93: Natężenie d) Filtr Zernikego 1 0.
- Page 94 and 95: Płaszczyzna wejściowa Soczewka Fi
- Page 96 and 97: Contrast a b Rys. 9.13 Rozkład nat
- Page 98 and 99: 10 Końcowe podsumowanie i wnioski
- Page 100 and 101: mania powietrza wywołanego fluktua
- Page 102 and 103: 25. Grabowski W., 1983: Measurement
- Page 104: 71. SCIFEAT pakiet procedur graficz
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