Powstawanie siły nośnej
Powstawanie siły nośnej – opisy głównych modeli fizycznych i niektórych metod obliczeń siły nośnej.
Siła nośna to składowa siły aerodynamicznej powstająca przy ruchu otoczonego płynem ciała (profilu) względem tego płynu, prostopadła do kierunku ruchu. Powstaje w wyniku korzystnego rozkładu ciśnień na profilu, a ten jest wynikiem tworzącego się wokół ciała wiru cyrkulacyjnego – dodatkowego pola prędkości przyspieszającego strugi na górnym obrysie profilu, a spowalniającego te na dolnym, co w odniesieniu do prawa Bernoulliego oznacza, że ciśnienie na górnym obrysie będzie mniejsze niż na dolnym – ciało jest „zasysane w górę”.
Wzór fizyczny
Siłę nośną określa wzór:
gdzie:
- – wytworzona siła nośna,
- – kąt natarcia, czyli kąt między pochyleniem profilu w stosunku do kierunku napływających strug powietrza,
- – współczynnik siły nośnej, dla profili skrzydeł samolotów wyliczalny teoretycznie, ale także wyznaczany empirycznie, zależny jest głównie od kąta natarcia ale także od kształtu profilu. Dla niewielkich kątów natarcia, poniżej 0,2 radiana, nachylenie charakterystyki jest stałe, bliskie 6 – teoretycznie
- – gęstość płynu (powietrze na poziomie morza 1,225 kg/m³),
- – powierzchnia ciała,
- – prędkość ciała względem płynu.
Wstępne wnioski z zasad dynamiki Newtona
Zgodnie z trzecią zasadą dynamiki Newtona, jeżeli płyn wywiera na ciało siłę nośną Pz (będącą wypadkową sił od ciśnień) w kierunku prostopadłym do prędkości płynu niezakłóconego, to ciało wywiera taką samą reakcję na płyn i powoduje zmianę pędu płynu w kierunku
gdzie:
- – masa płynu, której pęd uległ zmianie,
- – przyrost prędkości masy w kierunku
Nie znając ani nie można obliczyć z tego wzoru siły nośnej. Można jednak oszacować, jakiego rzędu jest wielkość obszaru płynu, który bierze udział w jej powstawaniu. W czasie przez obszar ten przepływa masa płynu m:
gdzie:
Załóżmy, że rozpatrywany jest ruch samolotu i przyjmijmy dopuszczalne oszacowanie, że za skrzydłem występuje odchylenie strumienia opływu o kąt rzędu kąta natarcia Mamy wtedy:
Niech to będzie duży samolot, lecący dość wolno:
- waga = = 200 000 kG,
- rozpiętość skrzydeł = 50 m,
- kąt natarcia = 0,1 rad,
- prędkość = 100 m/s.
Otrzymamy wtedy:
- = 1600 m².
Czyli dla rozpiętości skrzydeł 50 m obszar, którego wielkość oceniamy, będzie miał ponad 30 m szerokości w kierunku prostopadłym do rozpiętości skrzydeł.
Reasumując: Proces powstawania siły nośnej nie jest ograniczony do zjawisk na styku profilu z płynem, lecz zachodzi w obszarze, który ma wymiary rzędu wymiarów ciała, na które ta siła działa. Powyższe rozważania mają tylko poglądowy charakter i nie pozwalają określić rozkładu prędkości i ciśnień w tym obszarze.
Ogólne sformułowanie zadania
![](http://upload.wikimedia.org/wikipedia/commons/thumb/0/03/Mechanical_Potato.svg/250px-Mechanical_Potato.svg.png)
Aby określić oddziaływania płynu na izolowany profil należy rozpatrywać nieograniczony obszar płynu będącego w spoczynku w nieskończenie wielkiej odległości od ciała.
W układzie odniesienia związanym z nieruchomym płynem przyjmuje się następujące warunki brzegowe:
- składowe prędkości
gdy
- składowe tych prędkości prostopadłe do powierzchni ciała muszą być na tej powierzchni zerowe.
W wyniku zastosowania matematycznych narzędzi teorii pola i równań dynamiki płynów oblicza się pole ciśnień i pole prędkości
Oznaczając powierzchnię ciała przez a przez – normalną zewnętrzną, otrzymamy całkowitą reakcję ze wzoru:
przy czym oznacza tu ciśnienie w punktach powierzchni natomiast oznacza elementarny wektor siły stycznej do tej powierzchni, o kierunku i zwrocie wektora prędkości płynu w tym miejscu.
Składowa prostopadła do kierunku ruchu ciała, nazywana jest siła nośną.
Metody aerodynamiki klasycznej, zasada Kutty-Żukowskiego
Analityczne rozwiązania zadań teorii pola i dynamiki płynów są w większości przypadków niemożliwe, wymagane są zaawansowane metody numeryczne.
Dlatego już kilkadziesiąt lat temu zagadnieniu temu poświęciła głównie uwagę aerodynamika klasyczna – dział aerodynamiki, w którym jako podstawę wszystkich rozważań przyjęto model cieczy doskonałej (to znaczy płynu nieściśliwego i nielepkiego). W połączeniu z pewnymi dodatkowymi założeniami, wynikającymi z obserwacji rozpatrywanych zjawisk, umożliwiło to sformułowanie zasadniczych praw dotyczących sił wywieranych przez płyn na poruszające się względem niego ciała, w tym siły nośnej. Mimo niedoskonałości przyjętego modelu, doświadczenie potwierdza w większości przypadków uzyskane wyniki.
Dla cieczy doskonałej której ruch jednostajny rozpoczął się od stanu spoczynku, według zasady Thomsona ruch jest ruchem potencjalnym. W przyjętym jak uprzednio układzie odniesienia istnieje potencjał prędkości a składowe prędkości są:
Warunki brzegowe są jak uprzednio, z tym, że
Równanie ciągłości redukuje się do wzoru:
co prowadzi do:
- czyli równania Laplace’a.
Zadanie sprowadza się zatem do znalezienia funkcji spełniającej równania Laplace’a i wymienione już warunki brzegowe, jest więc zagadnieniem Neumana.
W wyniku jego rozwiązania otrzymujemy pole prędkości Wykorzystując całkę Cauchy-Lagrange’a, oznaczając przez ciśnienie w przepływie niezakłóconym, – prędkość w dowolnym punkcie pola, – potencjał jednostkowych sił masowych, a – gęstość, obliczymy pole ciśnień:
Wartości ciśnień wykorzystane w znanym już wzorze:
w którym pominięto już całkę reprezentującą siły tarcia pozwalają obliczyć wartość siły oddziaływania cieczy na ciało.
Ruch w układzie związanym z cieczą jest z natury ruchem nieustalonym, mimo stałej prędkości ciała względem tego układu. W celu dalszego uproszczenia pod względem formalnym wprowadza się zwykle układ ze związany z ciałem. Ruch jest wtedy ustalony, jeżeli ciało porusza się ruchem prostoliniowym ze stałą prędkością. Potencjał nie zależy od czasu, a warunki na pochodne mają teraz postać:
Układy i są równoległe, a ciało porusza się z prędkością:
względem układu nieruchomego.
Przepływ płaski
Wiele ważnych zagadnień aerodynamiki można traktować w przybliżeniu jako zagadnienia dwuwymiarowe, płaskie. Jest to szczególnie ważne i korzystne z punktu widzenia metod obliczeniowych, ponieważ można tu stosować potężny aparat teorii funkcji zmiennej zespolonej.
Przepływ płaski wokół płata nośnego ma podstawowe znaczenie w teorii profilu lotniczego. Ponieważ wszystkie parametry płaskiego ustalonego ruchu cieczy doskonałej zależą tylko od dwóch zmiennych niezależnych, można wprowadzić pewne szczególne funkcje tego ruchu: potencjał prądu i potencjał zespolony.
Jeśli ruch cieczy zachodzi w płaszczyźnie można napisać różniczkowe równanie każdej linii prądu w postaci:
czyli:
Lewa strona stanowi tu różniczkę zupełną pewnej funkcji pod warunkiem:
który jest spełniony w każdym płaskim bezźródłowym przepływie cieczy doskonałej.
Istnieje zatem funkcja Ψ, która nazywana będzie potencjałem prądu. Mamy dalej:
Udowadnia się, że jeżeli ruch jest bezwirowy, czyli rotacja jest zero:
to:
czyli pole płaskie bezwirowe jest potencjalne. Istnieje potencjał prędkości oznaczany i dalej:
Składowe te wyrażają się jednak również przez pochodne funkcji prądu. Otrzymuje się zatem:
Są to znane z teorii funkcji zespolonej warunki Cauchy-Riemanna.
Wynika z nich, że potencjał prędkości jest częścią rzeczywistą, a potencjał prądu – częścią urojoną pewnej funkcji holomorficznej. Funkcję tę nazywamy potencjałem zespolonym i oznaczamy zwykle w(z); zatem:
przy czym:
gdzie oznacza jednostkę urojoną.
Każda funkcja holomorficzna określa pewien płaski, ustalony przepływ bezwirowy cieczy doskonałej; równanie rodziny linii tego przepływu otrzymuje się porównując część urojoną funkcji do stałej, stanowiącej parametr.
Linie prądu są ortogonalne do linii stałego potencjału prędkości. Dla każdego przepływu, określonego potencjałem zespolonym, można stworzyć przepływ sprzężony, którego linie prądu są liniami stałego potencjału prędkości zadanego przepływu, a linie stałego potencjału prędkości – liniami prądu zadanego przepływu.
Przepływy, określone prostą pod względem analitycznym funkcją, nazywane są przepływami elementarnymi.
Przepływy określone potencjałem zespolonym, podlegają superpozycji:
Można nakładać na siebie dowolną liczbę przepływów, aby uzyskać pożądany przepływ wynikowy.
Dowolną linię prądu możemy traktować jako ściankę materialną, gdyż spełnia warunek styczności prędkości do ścianki.
Jeżeli w danym przepływie płaskim cieczy doskonałej pojawi się zamknięta linia prądu, można traktować przepływ na zewnątrz niej jako opływ profilu określonego tą linią.
Obliczając reakcję wywieraną na profil przez ciecz doskonałą, gdy znany jest potencjał zespolony w(z) określający opływ profilu, otrzymuje się dla siły nośnej wzór:
będący matematyczną formą prawa Kutty-Żukowskiego.
Prawo to mówi, że reakcja wywierana przez ciecz doskonałą na profil izolowany w ruchu ustalonym jest proporcjonalna do sumy cyrkulacji wszystkich wirów, objętych konturem profilu, do gęstości cieczy i do modułu prędkości cieczy w nieskończoności.
Odwzorowanie konforemne
![](http://upload.wikimedia.org/wikipedia/commons/thumb/e/e9/Odwzorowanie_op%C5%82ywu_okr%C4%99gu_na_op%C5%82yw_badanego_profilu_teoretycznego.svg/350px-Odwzorowanie_op%C5%82ywu_okr%C4%99gu_na_op%C5%82yw_badanego_profilu_teoretycznego.svg.png)
Znane jest twierdzenie, że można zawsze odwzorować konforemnie dowolną zamkniętą, nieprzecinającą się linię wraz z jej zewnętrzem, na inną zamkniętą linię i jej zewnętrze.
Niech będzie znany jest opływ pewnego profilu, np. kołowego, będący w nieskończoności przepływem jednorodnym. Płaszczyznę tego znanego opływu nazwijmy płaszczyzną pomocniczą zeta, a znany potencjał oznaczmy w(zeta).
Jeżeli w innej płaszczyźnie, zwanej płaszczyzną przepływu, jest pewien zadany profil, to na mocy przytoczonego stwierdzenia istnieje funkcja holomorficzna która odwzorowuje konforemnie opływ profilu kołowego na opływ profilu zadanego.
Taką samą nazwę, funkcji odwzorowującej, nadaje się też funkcji odwrotnej.
Jeśli zatem znany jest potencjał i funkcja odwzorowująca, to znany jest tym samym opływ danego profilu.
Opływ profilu kołowego o promieniu może być traktowany jako superpozycja trzech przepływów elementarnych: przepływu jednorodnego, wiru płaskiego i dipola.
Jedną z funkcji odwzorowujących jest funkcja:
zwana w aerodynamice funkcją Żukowskiego. Funkcja ta przekształca przepływ wokół profilu kołowego w przepływ wokół profilu Żukowskiego – pierwszego historycznie profilu teoretycznego.
Wartość cyrkulacji w potencjale profilu kołowego nie jest określona żadnym fizycznym uwarunkowaniem.
Warunek taki wynika dopiero z tak zwanej hipotezy Żukowskiego, czyli „hipotezy spływu na ostrzu”, która oznacza matematycznie warunek ograniczenia prędkości opływu na profilu.
Na tej podstawie można wyznaczyć jedyną wartość cyrkulacji, która spełnia ten warunek, pozostając przy tym na gruncie cieczy doskonałej. Warunek ten odnosi się do wszystkich profilów, nie tylko profilu Żukowskiego. Jeżeli profil nie ma ostrza, co się zawsze dzieje w praktyce, odpowiada mu punkt o największej krzywiźnie.
Przedstawiona metoda pozwala obliczyć rozkład prędkości i ciśnień na całym profilu teoretycznym.
Jeżeli zaś wrócimy do samej siły nośnej, to wzór na siłę nośną profilu Żukowskiego jest:
Pamiętając o wzorze z początku artykułu:
Uzyskujemy:
co bardzo dobrze zgadza się z rzeczywistością i stanowi jeden z najważniejszych wyników teorii Żukowskiego.
Tworzenie się opływu cyrkulacyjnego i ustanawianie warunków spływu na ostrzu może być obserwowane w czasie doświadczenia. W przepływie, w którym została w pewien sposób zapewniona wizualizacja, zwiększa się stopniowo prędkość. Widać, że w pobliżu ostrza profilu, ustawionego pod pewnym kątem natarcia, tworzy się wir, który w pewnej chwili odrywa się od profilu i odpływa niesiony prądem cieczy. Opływ od tego momentu oznacza się spływem na ostrzu.
Taki wir nosi nazwę wiru początkowego. Każda zmiana kąta natarcia lub prędkości powoduje oderwanie się wiru o odpowiednim kierunku, zmniejszającego lub zwiększającego cyrkulację.
Z punktu widzenia potrzeb praktyki profil Żukowskiego ma szereg wad i nie był stosowany na szerszą skalę, powstało potem wiele innych profili teoretycznych, eliminujących wady – na przykład moment przy zerowej sile nośnej. Ponieważ jednak był historycznie pierwszy profil, to przyczynił się do rozwoju olbrzymiej liczby profili teoretyczno-eksperymentalnych.
Wyniki aerodynamiki klasycznej, wraz z bogatym materiałem doświadczalnym, są aktualne do dziś w zakresie małych prędkości lotu. Mogą być też stosowane przy użyciu odpowiednich transformacji (np. transformacja Prandtla-Glauerta) także w zakresie dużych, ale poddźwiękowych prędkości lotu.
Pewne przejściowe trudności sprawiała tu tak zwana „bariera dźwięku”.
Równanie siły nośnej
![](http://upload.wikimedia.org/wikipedia/commons/thumb/8/80/Szkic_do_wyprowadzenia_rownania_sily_nosnej.svg/350px-Szkic_do_wyprowadzenia_rownania_sily_nosnej.svg.png)
Istnieje mniej formalnie poprawna od metody Kutty-Żukowskiego, lecz łatwiejsza droga prowadząca do określenia współczynnika siły nośnej, stosująca bezpośrednio prawa Newtona, podobnie jak w pierwszej części tego artykułu, ale przy zastosowaniu pojęcia cyrkulacji, czyli okrężnej składowej ruchu płynu wokół profilu. Metoda ta nie pozwala na uzyskanie innych korzyści – poza potwierdzeniem teoretycznej wartości współczynnika siły nośnej.
Wykorzystuje się tu hipotezę Kutty-Żukowskiego o spływie na ostrzu profilu, co pociąga za sobą konieczność istnienia okrężnego ruchu cieczy wokół takiego profilu. Warunek równowagi sił bezwładności i siły od gradientu ciśnienia pozwala określić prędkość ruchu okrężnego w funkcji odległości od środka krzywizny – prędkość ta jest odwrotnie proporcjonalna do odległości od tego punktu.
Prędkość obwodowa na krawędzi spływu i przed krawędzią natarcia, w odległości od środka obrotu jest:
a w punkcie o współrzędnych x,z jest zatem:
W płaszczyźnie zmiana pędu w kierunku osi skierowana w dół, wynosi:
gdzie przez oznaczono rozpiętość płata.
Po zmianie zmiennych przy zastosowaniu zależności:
otrzymuje się dla małych kątów natarcia:
przy czym wynik nie zależy od wyboru wartości współrzędnej
W płaszczyźnie będącej odpowiednikiem ale znajdującej się przed krawędzią natarcia mamy te same zależności i zmianę pędu o tym samym module, lecz skierowaną do góry osi
Sumaryczna zmiana pędu w jednostce czasu, równa działającej w kierunku dodatnim osi sile, czyli sile nośnej, jest więc:
Otrzymano dla dokładnie taką samą zależność jaka wynika z teorii Żukowskiego.
Uwagi końcowe
W zastosowaniu do lotnictwa przepływ płaski i wynikająca z niego teoria siły nośnej Kutty-Żukowskiego może być bezpośrednio stosowana tylko dla skrzydeł o nieskończonym wydłużeniu (to znaczy stosunku rozpiętości do długości profilu).
Ta teoria jest także przydatna do oszacowania siły nośnej skrzydła o dostatecznie dużym wydłużeniu i stałym profilu, bez zmiennego wzdłuż rozpiętości kąta ustawienia.
Należy zwrócić uwagę na to, że powstawanie siły nośnej nie jest tu związane z powstawaniem oporu, który – zgodnie z tą teorią – jest zerowy.
Gdy skrzydło ma skończone wydłużenie, przepływ staje się trójwymiarowy i występuje opór indukowany, będący funkcją siły nośnej.
Właściwości skrzydła o skończonym wydłużeniu określa aerodynamika skrzydła.
Bibliografia
- wykłady prof. J. Bukowskiego, prof. W. Fiszdona, prof. W. Prosnaka, prof. J. Rościszewskego, wygłoszone na MEiL Politechniki Warszawskiej.
- Batchelor G.K., An Introduction to Fluid Dynamics, Cambridge University Press, 1967.
Linki zewnętrzne
- http://www.regenpress.com/Gamma.htm