ciśnienie, prawo Pascala i wypórlepkośćliczba Reynoldsarównanie Bernoulliegosiła nośna, opór i współczynniki aerodynamicznewarstwa przyścienna i turbulencja
Fizyka Dział 10
Płyny: od ciśnienia do siły nośnej
Co utrzymuje w powietrzu latawiec, a w wodzie stalowy statek? Dlaczego w przewężeniu rury ciśnienie spada, mgła opada tak wolno, a piłka golfowa ma dołki? Ten dział zaczyna od płynu w spoczynku, przechodzi przez przepływ idealny i równanie Bernoulliego do lepkości, liczby Reynoldsa i oporu, a kończy na sile nośnej skrzydła i równowadze latawca. Na tym fundamencie stoi dynamika lotu w automatyce.
Czas
ok. 4 h teorii + ok. 5 h zadań
Zadania
25 zadań w treści + 8 feynmanowskich, z odpowiedziami
Źródło
Feynmana wykłady z fizyki, t. II, rozdz. 40–41; aerodynamika skrzydła: materiały uzupełniające
0 Prolog: co trzyma latawiec?
Latawiec o powierzchni 0,8 m² i masie 250 g wisi w wietrze wiejącym z prędkością 6 m/s (umiarkowany wiatr, 4 w skali Beauforta). Ciężar latawca to 2,5 N. Powietrze ma gęstość prawie tysiąc razy mniejszą od wody, a mimo to linka ciągnie dłoń z siłą prawie 9 N, a latawiec stoi wysoko nad głową.
Skala sił wynika z jednej wielkości: ciśnienia dynamicznegoq=21ρv2. Przy 6 m/s to 21⋅1,2⋅36≈22 Pa, a pomnożone przez powierzchnię daje ok. 18 N, siedem razy więcej niż ciężar. Wszystkie siły aerodynamiczne zapisuje się jako q razy powierzchnia razy współczynnik zależny od kształtu i ustawienia. Ten dział tłumaczy, skąd biorą się te współczynniki, a przykład 10.6 i symulacja 4 ustawią latawiec w równowadze.
Feynman poświęca płynom dwa rozdziały tomu II: „przepływ suchej wody” (ciecz bez lepkości, idealna dla matematyków) i „przepływ mokrej wody” (prawdziwa ciecz z lepkością). Pokazuje, że teoria suchej wody przewiduje zerowy opór i zerową siłę nośną dla wielu kształtów, co jest jawnie fałszywe, i że kluczem jest cienka warstwa płynu przy powierzchni ciała.FLP II-40-4, 41-4
Pytanie na start. Dlaczego nie da się wypić napoju słomką o długości 11 m, nawet z najlepszą pompą? Odpowiedź w części 1 i w zadaniu feynmanowskim 2.
1 Płyn w spoczynku: ciśnienie i wypór
Płyn (ciecz lub gaz) w spoczynku nie przenosi sił stycznych: gdyby je przenosił, zacząłby płynąć. Na każdy kawałek powierzchni działa więc tylko siła prostopadła, proporcjonalna do pola. Stosunek siły do pola to ciśnieniep, wielkość skalarna: w danym punkcie jest takie samo we wszystkich kierunkach.FLP II-40-1
Weźmy poziomą warstwę cieczy o grubości dh i polu A. Ciśnienie od dołu musi podtrzymać jej ciężar: Adp=ρgAdh. Dla cieczy o stałej gęstości:
Prawo · ciśnienie hydrostatyczne i wypór
p=p0+ρgh,Fwyporu=ρpłynuVzanurzoneg.
h to głębokość pod powierzchnią, na której panuje ciśnienie p0. Ciśnienie zależy tylko od głębokości, nie od kształtu naczynia. Zmiana ciśnienia w jednym miejscu przenosi się na całą ciecz (prawo Pascala). Wypór (prawo Archimedesa) wynika z różnicy ciśnień: na dolną powierzchnię ciała ciecz naciska mocniej niż na górną, a wypadkowa równa się ciężarowi wypartej cieczy.
Liczby warte zapamiętania: 1 atm = 101 325 Pa to ciśnienie słupa wody o wysokości 10,3 m albo słupa rtęci 760 mm. Na każde 10 m głębokości w wodzie przybywa ok. 1 atm; na dnie Rowu Mariańskiego (11 km) panuje ok. 1100 atm. W gazach ciśnienie maleje z wysokością wykładniczo, bo gęstość zależy od ciśnienia (dział 8 „Ciepło i teoria kinetyczna”).
Ciało pływa, gdy jego średnia gęstość jest mniejsza od gęstości cieczy; zanurza się wtedy w takiej części, żeby wyparta ciecz ważyła tyle co ono. Stalowy statek pływa, bo jego kadłub otacza ogromną objętość powietrza: średnia gęstość kadłuba z zawartością jest mniejsza od gęstości wody.
Zadania · część 1
5 zadań · ok. 30 min
Jakie ciśnienie (bezwzględne i nadciśnienie) panuje na głębokości 30 m w wodzie słodkiej? Z jaką siłą woda naciska na tamę o szerokości 100 m i wysokości wody 30 m?
Wskazówka
Ciśnienie rośnie liniowo z głębokością, więc średnie nadciśnienie na tamie to połowa wartości przy dnie.
Odpowiedź
Nadciśnienie ρgh≈294 kPa, bezwzględne ok. 395 kPa (3,9 atm). Siła: 21ρgh2⋅b≈441 MN (licząc tylko nadciśnienie, bo atmosfera naciska z obu stron).
Jaka część objętości góry lodowej (ρ=917kg/m3) wystaje nad wodę morską (1025kg/m3)?
Odpowiedź
Zanurzona część to 917/1025≈89,5%, wystaje ok. 10,5%.
Podnośnik hydrauliczny ma tłoki o polach 2 cm² i 200 cm². Jaką siłę trzeba przyłożyć do małego tłoka, by unieść samochód o masie 1200 kg? O ile musi przesunąć się mały tłok, by samochód podniósł się o 10 cm?
Odpowiedź
Ciśnienie jest to samo: F=11770N⋅2/200≈118 N. Objętość cieczy się nie zmienia: mały tłok przesuwa się o 10 m (w praktyce w wielu ruchach pompki). Praca jest ta sama po obu stronach.
Kostka drewna o krawędzi 10 cm i gęstości 600 kg/m³ pływa w wodzie. Jak głęboko jest zanurzona? Jaką siłą trzeba ją wcisnąć, by zanurzyła się całkowicie?
Odpowiedź
Zanurzenie 0,6⋅10=6 cm. Pełne zanurzenie: wypór 1000⋅10−3⋅9,81=9,81 N minus ciężar 5,89 N, czyli 3,9 N.
Dlaczego słomką nie da się wypić napoju z wysokości większej niż ok. 10 m? Co by się zmieniło na Evereście?
Odpowiedź
Ssanie tylko zmniejsza ciśnienie nad napojem w słomce; do góry pcha go ciśnienie atmosferyczne. Nawet przy próżni w słomce podniesie go najwyżej na p0/ρg≈10,3 m. Na Evereście (p≈34 kPa) tylko na ok. 3,5 m.
2 Przepływ idealny: ciągłość i Bernoulli
Ruch płynu opisujemy polem prędkości v(r,t): w każdym punkcie podajemy prędkość płynu, który akurat tam jest. Linie styczne do tego pola to linie prądu; w przepływie ustalonym (niezależnym od czasu) cząstki płynu poruszają się wzdłuż nich.FLP II-40-2, 40-3
2.1 Równanie ciągłości
Masa się nie tworzy ani nie znika, więc to, co wypływa z małej objętości, ubywa z jej wnętrza: ∂ρ/∂t+∇⋅(ρv)=0 (dywergencja, lekcja 3.5 Strumień, dywergencja i rotacja). Ciecze, a także gazy przy prędkościach znacznie mniejszych od prędkości dźwięku (do ok. 0,3 Macha gęstość zmienia się o kilka procent), są praktycznie nieściśliwe: ∇⋅v=0. W rurze lub strudze o przekroju A oznacza to, że strumień objętości Q=Av jest wszędzie taki sam: tam, gdzie rura się zwęża, płyn płynie szybciej.
Dla przepływu ustalonego wzdłuż linii prądu z tego równania (albo prościej, z bilansu energii dla odcinka strugi: praca sił ciśnienia na końcach równa się zmianie energii kinetycznej i potencjalnej) wynika:FLP II-40-3
Prawo · twierdzenie Bernoulliego
p+21ρv2+ρgz=constwzdłuz˙ linii prądu.
Warunki: płyn nieściśliwy, bez lepkości, przepływ ustalony. Jeśli przepływ jest bezwirowy, stała jest ta sama w całym przepływie. 21ρv2 to ciśnienie dynamiczne: o tyle ciśnienie w rurce zwróconej otworem pod prąd, gdzie płyn zostaje zatrzymany, przewyższa ciśnienie statyczne.
Zastosowania: woda wypływa z otworu na głębokości h pod powierzchnią z prędkością 2gh, jak ciało spadające z tej wysokości (prawo Torricellego). W przewężeniu rury ciśnienie spada (zwężka Venturiego). Rurka Pitota na samolocie mierzy różnicę między ciśnieniem całkowitym (otwór zwrócony do przodu) i statycznym (otwór boczny), równą 21ρv2; stąd prędkościomierz (przykład 10.2).
Symulacja 1 Zwężka Venturiego z dwiema rurkami manometrycznymi. W szerokiej części woda płynie z prędkością v1, w przewężeniu szybciej. Wysokość słupa w rurce pokazuje ciśnienie (nadciśnienie względem atmosfery). Zmieniaj prędkość i stopień przewężenia.
2.3 Cyrkulacja i wiry
Cyrkulacja to całka prędkości wzdłuż zamkniętej krzywej, Γ=∮v⋅dl (lekcja 3.4 Całki krzywoliniowe, pola potencjalne i różniczki zupełne), a wirowość to rotacja pola prędkości, Ω=∇×v. W płynie idealnym cyrkulacja wokół krzywej poruszającej się razem z płynem się nie zmienia (twierdzenie Kelvina), a linie wirowe płyną razem z płynem: dlatego pierścień dymu tak długo zachowuje kształt. Cyrkulacja wróci w części 5 jako klucz do siły nośnej.FLP II-40-4, 40-5
Zadania · część 2
5 zadań · ok. 30 min
Woda płynie rurą o średnicy 4 cm z prędkością 1 m/s, a potem przewężeniem o średnicy 2 cm. Jaka jest prędkość w przewężeniu i różnica ciśnień?
Odpowiedź
Pole maleje 4 razy: v2=4 m/s. Δp=21⋅1000⋅(16−1)=7,5 kPa.
Z jaką prędkością wypływa woda z otworu 1 m pod powierzchnią w beczce? Na jaką odległość poziomą doleci, jeśli otwór jest 1 m nad ziemią?
Odpowiedź
v=2gh≈4,4 m/s; czas spadania 2⋅1/9,81≈0,45 s, zasięg ok. 2 m.
Wiatr o prędkości 30 m/s (wichura) wieje nad płaskim dachem o powierzchni 100 m². Oszacuj, z jaką siłą może go podrywać, przyjmując różnicę ciśnień rzędu 21ρv2.
Odpowiedź
21⋅1,225⋅900≈550 Pa, razy 100 m²: ok. 55 kN, ciężar ok. 5,6 tony. Dlatego wichury zrywają dachy, a nie wgniatają ich do środka.
Pokaż, że strumień wody z kranu zwęża się ku dołowi. Jeśli na wylocie ma prędkość 0,5 m/s i średnicę 1 cm, jaką średnicę ma 10 cm niżej?
Odpowiedź
Woda przyspiesza w spadku: v=v02+2gh≈1,49 m/s. Av stałe: d=d0v0/v≈0,58 cm.
Wyprowadź twierdzenie Bernoulliego z zasady zachowania energii dla odcinka strugi między przekrojami 1 i 2, w czasie, w którym przez każdy przekrój przepływa objętość ΔV.
Odpowiedź
Praca sił ciśnienia: p1ΔV−p2ΔV. Zmiana energii: masa ρΔV „przeniesiona” z 1 do 2: ρΔV[21(v22−v12)+g(z2−z1)]. Przyrównanie daje p+21ρv2+ρgz stałe.
3 Lepkość
Teoria płynu idealnego ma poważną wadę, znaną jako paradoks d'Alemberta: ciało poruszające się ze stałą prędkością w takim płynie nie doznaje żadnego oporu. Prawdziwe płyny mają lepkość: sąsiednie warstwy poruszające się z różnymi prędkościami ciągną się nawzajem. Między dwiema płytami, z których górna przesuwa się z prędkością u nad dolną w odległości d, potrzebna jest siła na jednostkę pola ηu/d.FLP II-41-1
Prawo · lepkość newtonowska
τ=ηdydu,ν=ρη.
τ to naprężenie styczne (siła na jednostkę pola), η lepkość dynamiczna (Pa·s), ν lepkość kinematyczna (m²/s). Woda w 20 °C: η=1,0 mPa·s, ν=1,0⋅10−6m2/s. Powietrze: η=1,8⋅10−5 Pa·s, ale ν=1,5⋅10−5m2/s, piętnaście razy więcej niż woda. Miód: kilka Pa·s.
Lepkość gazu wynika z ruchu cząsteczek (dział 8 „Ciepło i teoria kinetyczna”): cząsteczki przeskakujące między warstwami przenoszą pęd. Dlatego lepkość gazów rośnie z temperaturą, a cieczy maleje. Dopisanie siły lepkości do równania Eulera daje równanie Naviera–Stokesa:FLP II-41-2
ρ(∂t∂v+(v⋅∇)v)=−∇p+η∇2v+ρg,∇⋅v=0.
To równanie opisuje prawie wszystkie przepływy, od kropli po huragany, ale rozwiązuje się je dokładnie tylko w nielicznych przypadkach. Dwa najważniejsze:
Przepływ Poiseuille'a w rurze o promieniu R i długości L: profil prędkości jest paraboliczny, v(r)=4ηLΔp(R2−r2), a strumień objętości Q=8ηLπR4Δp. Czwarta potęga promienia sprawia, że zwężenie naczynia krwionośnego o 20% zmniejsza przepływ ponad dwukrotnie (przykład 10.3).
Opór Stokesa na małą kulę o promieniu a poruszającą się powoli: F=6πηav (dział 8, ruchy Browna). Kropelka mgły o średnicy 10 μm opada w powietrzu z prędkością ok. 3 mm/s, dlatego mgła wisi godzinami (przykład 10.4).
Zadania · część 3
5 zadań · ok. 30 min
Płyta o polu 0,5 m² sunie z prędkością 0,2 m/s po warstwie oleju (η=0,1 Pa·s) grubości 1 mm. Jaka siła jest potrzebna?
Odpowiedź
F=ηAu/d=0,1⋅0,5⋅0,2/0,001=10 N.
O ile zmieni się przepływ przez rurę, gdy podwoimy jej średnicę przy tym samym spadku ciśnienia? A gdy podwoimy długość?
Odpowiedź
Q∝R4: 16 razy więcej. Podwojenie długości: dwa razy mniej.
Kulka stalowa o średnicy 2 mm (ρ=7800kg/m3) opada w glicerynie (ρ=1260kg/m3, η=1,4 Pa·s). Oblicz prędkość graniczną ze wzoru Stokesa i sprawdź, czy wzór był uprawniony (Re < 1).
Odpowiedź
v=9η2a2(ρk−ρc)g=12,62⋅10−6⋅6540⋅9,81≈1,0 cm/s. Re=1260⋅0,010⋅0,002/1,4≈0,018 ✓. Tak mierzy się lepkość wiskozymetrem kulkowym.
Wyprowadź z profilu parabolicznego, że średnia prędkość w przepływie Poiseuille'a to połowa prędkości na osi.
Dlaczego olej silnikowy gęstnieje na mrozie, a lepkość powietrza w upale rośnie?
Odpowiedź
W cieczy cząsteczki są stale w kontakcie i przeskok warstwy wymaga pokonania bariery energetycznej; w wyższej temperaturze częściej się to udaje (czynnik Boltzmanna, dział 8), więc lepkość maleje. W gazie pęd przenoszą swobodne cząsteczki, a w wyższej temperaturze są szybsze i przenoszą go więcej: lepkość rośnie (w najprostszym modelu jak T, w rzeczywistości nieco szybciej).
4 Liczba Reynoldsa, warstwa przyścienna i opór
Porównajmy w równaniu Naviera–Stokesa wyraz bezwładności ρ(v⋅∇)v∼ρv2/L z wyrazem lepkości η∇2v∼ηv/L2. Ich stosunek to jedyna bezwymiarowa liczba, od której zależy charakter przepływu wokół ciała o danym kształcie (lekcja 0.5 Jednostki, analiza wymiarowa, szacowanie i niepewność pomiaru):FLP II-41-3
Definicja · liczba Reynoldsa
Re=ηρvL=νvL.
Dwa przepływy wokół geometrycznie podobnych ciał przy tej samej liczbie Reynoldsa (i małej liczbie Macha) są podobne: różnią się tylko skalą. Na tym opiera się badanie modeli w tunelach aerodynamicznych. Bakteria pływa przy Re∼10−4, człowiek idący w powietrzu przy ∼105, samochód przy ∼107.
Feynman opisuje, jak zmienia się przepływ wokół walca, gdy rośnie Re.FLP II-41-4 Przy Re<1 linie prądu opływają walec symetrycznie (przepływ pełzający, rządzi lepkość). Około Re≈10 za walcem pojawiają się dwa stojące wiry. Powyżej ok. 40–50 wiry odrywają się na przemian z jednej i drugiej strony, tworząc ścieżkę wirową Kármána. (Poza Feynmanem: odrywają się z częstotliwością ok. 0,2v/d, co słychać jako „śpiew” drutów na wietrze.) Przy Re∼103–105 ślad za walcem jest turbulentny, a w okolicy Re≈3⋅105 opór gwałtownie spada (kryzys oporu, symulacja 2).
4.1 Warstwa przyścienna
Nawet przy ogromnych Re lepkość nie znika z gry. Płyn tuż przy powierzchni ciała przylega do niej (ma prędkość zerową), a w cienkiej warstwie przyściennej prędkość rośnie do wartości przepływu zewnętrznego. Pojęcie warstwy przyściennej wprowadził Prandtl (1904); dla warstwy laminarnej jej grubość to ok. 5x/Rex (rozwiązanie Blasiusa, 1908), gdzie x to odległość od krawędzi natarcia: na skrzydle o cięciwie 2 m przy 50 m/s wychodzi kilka milimetrów, a warstwa turbulentna, jaka rzeczywiście tam występuje, ma kilka centymetrów. Poza tą warstwą płyn zachowuje się jak idealny, ale to warstwa decyduje o oporze: gdy przy tylnej części ciała ciśnienie rośnie, spowolniony płyn w warstwie nie daje rady płynąć „pod górę ciśnienia” i warstwa się odrywa, zostawiając za ciałem szeroki, wirujący ślad o niskim ciśnieniu.FLP II-41-4, 41-5
4.2 Opór
Siłę oporu przy dużych Re zapisuje się wzorem z analizy wymiarowej:
FD=21CDρv2A,
gdzie A to pole przekroju poprzecznego, a współczynnik oporu CD zależy od kształtu i (słabo) od Re. Typowe wartości: płaska płyta prostopadła do przepływu ok. 1,2; walec 1,0–1,2; kula 0,4–0,5 (przed kryzysem oporu); samochód osobowy 0,25–0,35; kształt opływowy (kropla) ok. 0,04. Kula w przepływie pełzającym ma CD=24/Re, co jest wzorem Stokesa zapisanym w tej samej postaci.
Symulacja 2 Współczynnik oporu kuli w funkcji liczby Reynoldsa (dopasowanie do pomiarów, skala logarytmiczna) i prawo Stokesa. Wybierz płyn, średnicę kuli i prędkość: punkt pokazuje, w jakim reżimie jesteś i jaka jest siła oporu. Spróbuj: kropelka mgły w powietrzu, kulka w miodzie, piłka do golfa przy 70 m/s.
Zadania · część 4
5 zadań · ok. 35 min
Oblicz Re dla: samochodu (4 m, 30 m/s, powietrze), pływaka (2 m, 1,5 m/s, woda), bakterii (2 μm, 30 μm/s, woda).
Odpowiedź
8⋅106; 3⋅106; 6⋅10−5. Bakteria żyje w świecie, w którym bezwładność nie gra roli: przestaje płynąć natychmiast, gdy przestaje poruszać wicią.
Samochód ma CD=0,3 i pole przekroju 2,2 m². Jaką moc pochłania opór powietrza przy 108 km/h, a jaką przy 144 km/h?
Odpowiedź
P=21CDρAv3: przy 30 m/s ok. 11 kW, przy 40 m/s ok. 26 kW. Moc rośnie jak v3: o jedną trzecią szybciej to 2,4 raza więcej mocy.
Skoczek spadochronowy o masie 80 kg w pozycji „na brzuchu” ma CDA≈0,45 m². Jaką osiąga prędkość graniczną? Ile razy musi wzrosnąć CDA po otwarciu spadochronu, by spadał z prędkością 5 m/s?
Odpowiedź
v=2mg/(ρCDA)≈53 m/s (ok. 190 km/h). Dla 5 m/s: CDA≈51 m², ponad 100 razy więcej.
Linka latawca o średnicy 1,5 mm brzęczy w wietrze 6 m/s. Oszacuj częstotliwość dźwięku. Jaka jest liczba Reynoldsa?
Odpowiedź
f≈0,2v/d=0,2⋅6/0,0015=800 Hz. Re=6⋅0,0015/1,5⋅10−5=600: zakres ścieżki wirowej Kármána ✓.
Oszacuj grubość warstwy przyściennej na końcu skrzydła o cięciwie 2 m przy 50 m/s ze wzoru δ≈5x/Rex. Dlaczego wynik, choć mały, ma tak duże znaczenie?
Odpowiedź
Re=50⋅2/1,5⋅10−5≈6,7⋅106, δ≈10/2580≈3,9 mm (w praktyce warstwa jest już turbulentna i grubsza, kilka centymetrów). W tej warstwie gradient prędkości jest ogromny, więc to w niej działają siły lepkości, a jej oderwanie decyduje o oporze i o przeciągnięciu skrzydła.
5 Siła nośna: cyrkulacja, skrzydło i latawiec
Ta część wykracza poza wykłady Feynmana, który na przykładzie walca pokazuje tylko, że bez cyrkulacji nie ma siły nośnej, a z cyrkulacją się ona pojawia. Skrzydło ustawione pod niewielkim kątem do przepływu odchyla strumień powietrza w dół; z trzeciej zasady dynamiki powietrze pcha skrzydło w górę. Matematycznie: przepływ wokół skrzydła to opływ bez cyrkulacji plus wir z cyrkulacją Γ, przez co nad skrzydłem powietrze płynie szybciej, a pod nim wolniej, i z twierdzenia Bernoulliego ciśnienie na górze jest mniejsze.FLP II-40-4
Prawo · twierdzenie Kutty–Żukowskiego i współczynnik siły nośnej
L′=ρvΓ,L=21CLρv2S,CL≈2πα(cienki profil, mały kąt α w radianach).
L′ to siła nośna na jednostkę rozpiętości, S powierzchnia skrzydła. Cyrkulacja ustala się tak, by powietrze schodziło gładko z ostrej krawędzi spływu (warunek Kutty); przy starcie skrzydło zostawia za sobą „wir startowy” o przeciwnej cyrkulacji. Dla skrzydła o skończonym wydłużeniu AR=b2/S nachylenie jest mniejsze, ok. 2πAR/(AR+2), a powyżej kąta krytycznego (zwykle 12–16°) przepływ odrywa się od górnej powierzchni i siła nośna spada: przeciągnięcie.
Skończone skrzydło płaci za siłę nośną dodatkowym oporem indukowanym: na końcach skrzydła powietrze przepływa spod spodu na górę, tworząc wiry brzegowe, które odchylają strumień w dół i pochylają siłę aerodynamiczną do tyłu. Współczynnik oporu to w przybliżeniu
CD=CD0+πeARCL2,
gdzie CD0 to opór przy zerowej sile nośnej, a e≈0,7–0,9 to współczynnik Oswalda. Stosunek L/D (doskonałość) mówi, ile metrów samolot przeleci w locie ślizgowym na każdy metr wysokości: Cessna ok. 9, szybowiec 40–60, albatros ok. 20.
Symulacja 3 Skrzydło lekkiego samolotu (powierzchnia 16,2 m², masa samolotu 1100 kg, model uproszczony): współczynniki siły nośnej i oporu w funkcji kąta natarcia, z przeciągnięciem przy 14°. Zmieniaj kąt natarcia i prędkość i znajdź lot poziomy (siła nośna równa ciężarowi). Zmieniaj wydłużenie skrzydła i obserwuj doskonałość.
5.1 Latawiec
Latawiec to skrzydło przy dużym kącie natarcia, zwykle 20–45°, uwiązane linką. Siła aerodynamiczna na płaską powierzchnię działa w przybliżeniu prostopadle do niej; jej składowa prostopadła do wiatru to siła nośna L, a składowa wzdłuż wiatru to opór D. W równowadze siła aerodynamiczna, ciężar i naciąg linki sumują się do zera. Kąt natarcia ustala uzda (miejsce przyczepienia linki), a ogon stabilizuje latawiec, przesuwając środek oporu do tyłu.
Prosty model płaskiej płyty przy dużych kątach natarcia: siła normalna FN=qSCN z CN≈1,2sinα (przy 90∘ płyta ma CD≈1,2). Wtedy L=FNcosα, D=FNsinα. To przybliżenie rzędu kilkudziesięciu procent, ale dobrze pokazuje kompromis: umiarkowanie mały kąt natarcia (ok. 15°) daje wysoko stojący latawiec, duży daje silny naciąg, a zbyt mały – za małą siłę nośną (symulacja 4, przykład 10.6).
Symulacja 4 Latawiec 0,8 m² o masie 250 g w równowadze sił, w modelu płaskiej płyty. Strzałki: siła nośna (prostopadle do wiatru), opór (z wiatrem), ciężar i naciąg linki (wzdłuż linki, której ciężar pomijamy). Zmieniaj prędkość wiatru i kąt natarcia ustawiony uzdą.
Zadania · część 5
5 zadań · ok. 35 min
Samolot o masie 1100 kg i powierzchni skrzydeł 16,2 m² leci poziomo z prędkością 55 m/s na poziomie morza. Jaki jest współczynnik siły nośnej?
Odpowiedź
CL=2mg/(ρv2S)=2⋅10791/(1,225⋅3025⋅16,2)≈0,36.
Ten sam samolot ma CL,max=1,6 (z klapami). Z jaką najmniejszą prędkością może lecieć poziomo?
Odpowiedź
vmin=2mg/(ρSCL,max)≈26 m/s (ok. 94 km/h): prędkość przeciągnięcia.
Jaka cyrkulacja jest potrzebna wokół skrzydła o rozpiętości 11 m, by przy 55 m/s utrzymać 1100 kg (rozkład siły nośnej przyjmij jednorodny)?
Odpowiedź
L′=mg/b≈981 N/m, Γ=L′/(ρv)≈14,6m2/s.
Jak zmienia się opór indukowany, gdy samolot leci dwa razy szybciej przy tym samym ciężarze? A opór pasożytniczy CD0? Dlaczego istnieje prędkość najmniejszego oporu?
Odpowiedź
CL∝1/v2, więc Di=qSCL2/(πeAR)∝v2⋅v−4=v−2: cztery razy mniejszy. Opór pasożytniczy ∝v2: cztery razy większy. Suma ma minimum, gdy oba są równe.
Latawiec z prologu (0,8 m², 0,25 kg, CN=1,2sinα) przy kącie natarcia 45°. Jaki jest najmniejszy wiatr, przy którym się utrzyma?
Odpowiedź
Siła nośna L=qS⋅1,2sinαcosα jest największa przy 45∘: L=0,6qS. Warunek L=mg: q=2,45/(0,6⋅0,8)≈5,1 Pa, v=2q/ρ≈2,9 m/s.
6 Metoda i przykłady rozwiązane
Metoda dla zadań z płynów
Spoczynek czy ruch? W spoczynku: p=p0+ρgh i wypór. W ruchu: równanie ciągłości Av=const i Bernoulli wzdłuż linii prądu.
Policz Re, zanim wybierzesz wzór na opór: Re<1 Stokes, Re≫1000 wzór kwadratowy z CD.
Siły aerodynamiczne: zawsze q=21ρv2 razy powierzchnia razy współczynnik; sprawdź, którą powierzchnię (przekrój czy skrzydło) wiąże podany współczynnik.
Równowaga: rysuj wszystkie siły (nośna prostopadle do przepływu, opór wzdłuż, ciężar, naciąg lub ciąg) i rozkładaj na składowe.
Sprawdź założenia: nieściśliwość (Mach < 0,3), brak lepkości poza warstwą przyścienną, przepływ ustalony.
Sprawdzenie: jednostki, rząd wielkości (1 atm = 10 m wody, q przy 30 m/s to ok. 550 Pa).
Przykład 10.1 Dlaczego pływa stalowy statek
Barka ma prostopadłościenny kadłub 50 m × 10 m. Ile ładunku przyjmie, zanurzając się z 1 m do 3 m w wodzie morskiej? Jaka jest średnia gęstość barki z ładunkiem przy zanurzeniu 3 m, jeśli kadłub ma wysokość 5 m?
Wyparta objętość przy zanurzeniu 3 m: 50⋅10⋅3=1500 m³, masa wody 1500⋅1025≈1,54⋅106 kg. Przy 1 m: 5,1⋅105 kg.
Ładunek to różnica: ok. 1000 ton.
Średnia gęstość całości (objętość kadłuba 2500 m³): 1,54⋅106/2500≈615kg/m3, mniej niż woda.
Wynik: barka przyjmie ok. 1000 t ładunku; pływa, bo jej średnia gęstość z powietrzem w środku jest mniejsza od gęstości wody.
Sprawdzenie: każdy centymetr zanurzenia to 50⋅10⋅0,01⋅1025≈5 t ładunku; 200 cm to ok. 1000 t ✓.
Przykład 10.2 Prędkościomierz samolotu
Rurka Pitota pokazuje różnicę ciśnienia całkowitego i statycznego 600 Pa. Jaka jest prędkość samolotu na poziomie morza, a jaka na wysokości 3000 m (gęstość 0,909 kg/m³)?
Bernoulli między punktem spiętrzenia i przepływem swobodnym: Δp=21ρv2.
Poziom morza: v=2⋅600/1,225≈31,3 m/s.
3000 m: v=2⋅600/0,909≈36,3 m/s.
Wynik: ten sam odczyt to 31 m/s przy ziemi i 36 m/s na 3000 m. Prędkościomierz jest wyskalowany dla gęstości przy ziemi i pokazuje prędkość przyrządową (IAS); prędkość rzeczywista (TAS) jest na wysokości większa.
Sprawdzenie: dla pilota to zaleta: siła nośna zależy od 21ρv2, czyli od wskazania prędkościomierza, więc prędkość przeciągnięcia „na wskaźniku” jest taka sama na każdej wysokości ✓.
Przykład 10.3 Zwężone naczynie krwionośne
Blaszka miażdżycowa zmniejsza promień tętnicy o 20% na pewnym odcinku. Ile razy wzrasta opór hydrauliczny tego odcinka? O ile musiałaby wzrosnąć różnica ciśnień, by utrzymać ten sam przepływ, gdyby cały opór tętnicy pochodził z tego odcinka?
Poiseuille: Q∝R4Δp, opór ∝1/R4.
1/0,84≈2,44: opór rośnie prawie dwa i pół raza.
Przy tej samej różnicy ciśnień przepływ spada do 0,84≈41%; żeby go utrzymać, różnica ciśnień musiałaby wzrosnąć 2,44 raza.
Wynik: 20% mniejszy promień to 2,4 raza większy opór. Organizm kompensuje to częściowo rozszerzaniem innych naczyń i podnoszeniem ciśnienia.
Sprawdzenie: prawo Poiseuille'a dotyczy przepływu laminarnego w sztywnej rurze; krew nie jest cieczą newtonowską, a tętnice są sprężyste, więc to oszacowanie, ale kierunek i siła efektu są właściwe ✓.
Przykład 10.4 Mgła i deszcz
Oblicz prędkość opadania kropelki mgły o średnicy 10 μm i kropli deszczu o średnicy 2 mm w powietrzu. Sprawdź, który wzór na opór stosuje się w każdym przypadku.
Deszcz: przy Re rzędu setek wzór kwadratowy z CD≈0,5: mg=21CDρav2πa2, czyli v=3CDρa8aρwg≈6,5 m/s; Re≈900 ✓.
Wynik: mgła opada z prędkością 3 mm/s (ok. 10 m na godzinę), kropla deszczu 6,5 m/s. Promień różni się 200 razy, prędkość ponad 2000 razy.
Sprawdzenie: zmierzona prędkość graniczna kropli o średnicy 2 mm to ok. 6,5 m/s ✓. Duże krople (powyżej ok. 5 mm) spłaszczają się i rozpadają, więc deszcz rzadko spada szybciej niż 9 m/s.
Przykład 10.5 Lekki samolot w locie poziomym
Samolot: 1100 kg, skrzydło 16,2 m², wydłużenie 7,3, e=0,8, CD0=0,03. Lot poziomy z prędkością 55 m/s na poziomie morza. Oblicz kąt natarcia, opór, doskonałość i moc potrzebną do lotu.
CL=0,36 (zadanie 5.1). Nachylenie krzywej nośnej 2π⋅7,3/9,3≈4,93 /rad, więc α≈0,073 rad ≈4,2∘ (od kąta zerowej siły nośnej).
CD=0,03+0,362/(π⋅0,8⋅7,3)≈0,03+0,007=0,037.
L/D≈9,7; opór mg/(L/D)≈1,1 kN; moc Dv≈61 kW.
Wynik: ok. 4° kąta natarcia, opór 1,1 kN, doskonałość prawie 10 i ok. 61 kW mocy użytecznej, czyli ok. 75 kW na wale silnika przy sprawności śmigła 0,8.
Sprawdzenie: lekkie samoloty tej klasy mają silniki o mocy 110–130 kW i przy ok. 60% mocy przelatują z prędkością ok. 55–60 m/s ✓.
Przykład 10.6 Latawiec w równowadze
Latawiec z prologu (0,8 m², 0,25 kg) przy wietrze 6 m/s i kącie natarcia 30°, w modelu CN=1,2sinα. Oblicz siłę nośną, opór, naciąg linki i kąt, pod jakim linka wznosi się nad ziemię (linkę traktujemy jako prostą i nieważką).
q=21⋅1,225⋅36≈22,1 Pa; FN=qS⋅1,2sin30∘≈10,6 N.
L=FNcos30∘≈9,2 N, D=FNsin30∘≈5,3 N, ciężar 2,45 N.
Linka musi zrównoważyć wypadkową: składowa pozioma 5,3 N, pionowa 9,2−2,45=6,7 N. Naciąg ≈8,6 N, kąt arctg(6,7/5,3)≈52∘.
Wynik: naciąg ok. 8,6 N, linka pod kątem ok. 52° do ziemi. Przy kącie natarcia 15° naciąg spada do ok. 3 N, a linka staje bardziej stromo (ok. 64°).
Sprawdzenie: przy wietrze 3 m/s siła nośna przy 30° (2,3 N) nie równoważy ciężaru, a model daje najmniejszy wiatr ok. 2,9 m/s przy 45° (zadanie 5.5) ✓.
7 Błędne intuicje, których trzeba się oduczyć
„Powietrze nad skrzydłem musi dogonić powietrze pod spodem.”
Nie musi i tego nie robi: dociera do krawędzi spływu wcześniej. Siła nośna wynika z kąta natarcia i cyrkulacji.
„Szybciej płynący płyn zawsze ma niższe ciśnienie.”
Tylko wzdłuż jednej linii prądu (lub w przepływie bezwirowym) i bez lepkości. Struga z wentylatora ma ciśnienie atmosferyczne.
„Ciężkie przedmioty toną, lekkie pływają.”
Decyduje średnia gęstość, nie ciężar. Lotniskowiec o masie 100 000 t pływa, gwóźdź tonie.
„Ciśnienie w wodzie naciska tylko w dół.”
Ciśnienie działa jednakowo we wszystkich kierunkach; dlatego woda wypływa z boku butelki i naciska na dno łodzi od spodu.
„Gęsta ciecz jest lepka.”
Rtęć jest 13,6 razy gęstsza od wody, a ma podobną lepkość (1,5 mPa·s). Olej jest lżejszy od wody i sto razy bardziej lepki.
„Im szybciej, tym zawsze większy opór.”
W kryzysie oporu kuli (Re ok. 3⋅105) siła oporu w pewnym zakresie maleje przy rosnącej prędkości.
„Skrzydło musi być wypukłe od góry.”
Płaska płyta też daje siłę nośną, a samoloty akrobacyjne mają symetryczne profile i latają do góry nogami, zwiększając kąt natarcia.
8 Zadania feynmanowskie
Te zadania wymagają więcej myślenia niż liczenia.
Zadania feynmanowskie
8 zadań · ok. 2 godziny
Myszy i ludzie. Przy oporze kwadratowym prędkość graniczna to 2mg/(ρCDA). Pokaż, że dla zwierząt podobnych kształtem rośnie jak L (gdzie L to rozmiar). Jeśli człowiek (1,8 m) spada z prędkością ok. 55 m/s, to z jaką prędkością spada mysz (3 cm)? Haldane pisał („On Being the Right Size”, 1926), że mysz zrzucona do szybu kopalni wstaje i odchodzi. Dlaczego?
Rozwiązanie
m∝L3, A∝L2, więc v∝L. Mysz: 550,03/1,8≈7 m/s (w rzeczywistości jeszcze mniej, bo futro i rozłożone łapki zwiększają CDA).
Energia na jednostkę masy przy uderzeniu ∝v2∝L, a wytrzymałość kości i tkanek na jednostkę masy rośnie przy małych rozmiarach (przekrój do masy ∝1/L). Mysz uderza 60 razy łagodniej na kilogram niż człowiek, a jest dużo bardziej odporna.
Słomka i pompa. Studniarze i górnicy w XVII w. nie mogli wypompować wody z głębokości większej niż ok. 10 m pompami ssącymi. Wyjaśnij to. Jak rozwiązano problem?
Rozwiązanie
Pompa ssąca tylko obniża ciśnienie nad wodą; wodę podnosi ciśnienie atmosferyczne, a ono zrównoważy słup wody o wysokości najwyżej p0/ρg≈10,3 m. Ten problem doprowadził Torricellego do barometru (1643).
Rozwiązanie: pompy tłoczące umieszczone na dole (pchające wodę do góry) albo kilka pomp ssących jedna nad drugą, każda podnosząca wodę o kilka metrów.
Podkręcona piłka. Piłka obracająca się w locie skręca (efekt Magnusa). Wyjaśnij to za pomocą cyrkulacji. W którą stronę skręci piłka nożna kopnięta tak, że patrząc z góry obraca się zgodnie z ruchem wskazówek zegara?
Rozwiązanie
Obracająca się powierzchnia wlecze za sobą powietrze w warstwie przyściennej i przesuwa punkty oderwania: przepływ wokół piłki dostaje cyrkulację. Z twierdzenia Kutty–Żukowskiego pojawia się siła prostopadła do prędkości. Prościej: ślad za piłką jest odchylony w jedną stronę, więc piłka jest pchana w drugą.
Obrót zgodny z ruchem wskazówek zegara (patrząc z góry) przy locie „do przodu”: lewa strona piłki porusza się do przodu (pod wiatr), prawa do tyłu (z wiatrem). Powietrze szybciej opływa stronę, która porusza się z wiatrem, czyli prawą; siła działa w prawo.
Start jumbo jeta. Boeing 747 o masie 400 t i powierzchni skrzydeł 511 m² odrywa się od ziemi przy ok. 80 m/s. Jaki CL jest potrzebny? Dlaczego przy starcie wysuwa się klapy, a w przelocie je chowa?
Rozwiązanie
CL=2mg/(ρv2S)≈2,0. Gładkie skrzydło osiąga najwyżej ok. 1,4–1,6, zanim nastąpi przeciągnięcie.
Klapy i sloty zwiększają krzywiznę i powierzchnię skrzydła oraz opóźniają oderwanie: CL,max rośnie do 2,5–3. W przelocie (250 m/s na 11 km, gdzie ρ≈0,36) potrzebne CL to ok. 0,5–0,7, a klapy dawałyby tylko dodatkowy opór.
Tunel wodny. Chcemy zbadać opływ samochodu (4 m, 30 m/s w powietrzu) na modelu w skali 1:10. Jaką prędkość powietrza w tunelu trzeba by ustawić, by zachować Re? Dlaczego to niemożliwe i jak pomaga tunel wodny?
Rozwiązanie
Re=8⋅106; model 0,4 m w powietrzu wymagałby 300 m/s, blisko prędkości dźwięku: powietrze przestaje być nieściśliwe i przepływ jest zupełnie inny.
W wodzie (ν 15 razy mniejsze): v=8⋅106⋅10−6/0,4=20 m/s, osiągalne. W praktyce często godzi się na mniejsze Re i poprawia wyniki, bo CD powyżej 106 zmienia się powoli.
Śpiewające druty. Przewody trakcji elektrycznej i liny masztów „śpiewają” na wietrze. Przy jakiej częstotliwości śpiewa drut o średnicy 2 mm przy wietrze 10 m/s? Co się stanie, jeśli ta częstotliwość pokryje się z częstotliwością własną liny?
Rozwiązanie
f≈0,2v/d=1000 Hz (liczba Strouhala ok. 0,2 w szerokim zakresie Re).
Odrywające się na przemian wiry działają na linę okresową siłą poprzeczną. Przy zgodności częstotliwości pojawia się rezonans (dział 6 „Drgania, rezonans i układy liniowe”); wiry potrafią się wtedy „zablokować” na częstotliwości drgań. Dlatego na wysokich kominach montuje się spiralne żebra, które psują regularność wirów.
Latawiec w ciszy. Jest bezwietrznie. Dlaczego latawiec da się podnieść, biegnąc? Z jaką prędkością trzeba biec dla latawca z przykładu 10.6? Co się dzieje z latawcem, gdy przestajesz biec?
Rozwiązanie
Siły aerodynamiczne zależą od prędkości względem powietrza. Biegnąc z prędkością u, tworzysz „wiatr” u względem latawca.
Według modelu potrzeba ok. 2,9 m/s przy najlepszym kącie natarcia (zadanie 5.5): szybki marsz lub trucht.
Po zatrzymaniu prędkość względna maleje, latawiec ślizga się w dół, a jego opór przez chwilę jeszcze napina linkę. Bez wiatru spadnie.
Kropla, która nie spada. Chmura składa się z kropelek o średnicy ok. 10–20 μm. Prądy wznoszące w chmurze kłębiastej mają prędkość kilku metrów na sekundę. Dlaczego chmury nie opadają i co musi się stać, żeby spadł deszcz?
Rozwiązanie
Prędkość opadania kropli 10–20 μm to 3–12 mm/s (Stokes, v∝a2). Najsłabszy prąd wznoszący ją przewyższa.
Deszcz wymaga wzrostu kropli milion razy w objętości (do ok. 1–2 mm), przez zderzenia i zlewanie się kropel lub przez wzrost kryształków lodu kosztem kropel wody. Dopiero takie krople spadają z prędkością metrów na sekundę i pokonują prąd wznoszący.
9 Sprawdź się
Odpowiedz w myślach, zanim rozwiniesz odpowiedź.
1 Od czego zależy ciśnienie na dnie naczynia z wodą?
Tylko od głębokości (i ciśnienia nad powierzchnią): p=p0+ρgh; nie od kształtu naczynia.
2 Ile wynosi siła wyporu?
Ciężar wypartej cieczy: ρpłynuVzanurzoneg.
3 Jak zmienia się prędkość wody, gdy rura zwęża się do połowy pola przekroju?
Rośnie dwukrotnie (równanie ciągłości).
4 Co mierzy rurka Pitota?
Ciśnienie dynamiczne 21ρv2, czyli różnicę ciśnienia całkowitego i statycznego.
5 Co mówi liczba Reynoldsa?
Stosunek sił bezwładności do sił lepkości; przepływy o tej samej Re wokół podobnych ciał są podobne.
6 Jak opór kropelki mgły zależy od prędkości, a jak opór samochodu?
Przepływ odrywa się od górnej powierzchni, siła nośna spada, a opór rośnie.
8 Dlaczego szybowce mają długie, wąskie skrzydła?
Duże wydłużenie zmniejsza opór indukowany CL2/(πeAR), więc rośnie doskonałość L/D.
10 Karta wzorów i co dalej
Spoczynek i przepływ idealny
p=p0+ρgh; Fw=ρpłVg; 1 atm ≈ 10,3 m wody
ciągłość: ∇⋅v=0, Av=const
Bernoulli: p+21ρv2+ρgz=const
Torricelli: v=2gh; Pitot: Δp=21ρv2
cyrkulacja Γ=∮v⋅dl; Kutta–Żukowski L′=ρvΓ
Lepkość, opór, siła nośna
τ=ηdu/dy; ν=η/ρ; Re=vL/ν
Poiseuille: Q=πR4Δp/(8ηL); Stokes: F=6πηav
opór: FD=21CDρv2A; warstwa: δ≈5x/Rex
nośna: L=21CLρv2S, CL≈2παAR+2AR
CD=CD0+CL2/(πeAR); Strouhal: f≈0,2v/d
Do przeczytania u Feynmana
Obowiązkowo: t. II, rozdz. 40 (przepływ „suchej wody”) i 41 (przepływ „mokrej wody”).
Dla chętnych: aerodynamika skrzydła wykracza poza Feynmana; przystępne wprowadzenie dają materiały NASA Glenn Research Center („Beginner's Guide to Aeronautics”) i podręczniki aerodynamiki, np. J. D. Anderson, Fundamentals of Aerodynamics.
Następny dział: Fale i dźwięk. Ściśliwość powietrza, którą tu pomijaliśmy, pozwala mu przenosić fale ciśnienia. Zobaczymy, jak z równań ruchu ośrodka wynika równanie falowe, skąd bierze się prędkość dźwięku i dlaczego struny i piszczałki brzmią harmonicznie. A siły aerodynamiczne z tego działu wrócą w automatyce, w równaniach ruchu samolotu.