Курсовая работа: Продольное и поперечное обтекание тел вращения
получим (3)
Будем искать частное решение этого уравнения в виде произведения двух функций от переменных l и m в отдельности
j = L( l ) M( m ) ; (4)
тогда в уравнении (2) переменные разделятся и из равенства
в силу независимости l и m будет следовать, что каждая из частей равенства должна быть постоянной. Полагая эту постоянную равной n (n+1) , где n – целое положительное число, получим для определения L(l) и М(m) два обыкновенных линейных уравнения второго порядка лежандрова типа
(5)
Этим уравнениям удовлетворяют[3] два класса независимых решений:
1) функции Лежандра 1-го рода – полиномы Лежандра Pn (х) , определяемые равенствами
2)
P0 (x) = 1, Р1 (х) = х, P2 (x) = 0.5 (Зх2 -1), P3 (x) = 0.5 (5x3 -3x),…
и рекуррентным соотношением для вычисления последующих полиномов
(n + 1) Pn +1 (х) = (2n + 1) хРn (х) – nРn-1 (х);
2) функции Лежандра 2-го рода Qn (х) , определяемые равенствами
и рекуррентным соотношением
(n + 1) Qn+1 (х) = (2n + 1) xQn (х) – nQn-1 (х),
совпадающим с предыдущим соотношением для полиномов Лежандра.
Представим решение уравнения (3) как сумму двух потенциалов: 1) потенциала j¥ однородного потока, набегающего на тело со скоростью U ¥ ; этот потенциал по первой из формул (1) будет равенj ¥ = U ¥ x = U ¥ c l m . и 2) потенциала j' скоростей возмущений, который выразим суммой частных решений (4).
Функция Pn (х) , как полином n -й степени, обращается в бесконечность при бесконечно возрастающем аргументе, функция же Qn (х) при этом стремится к нулю, но зато логарифмически бесконечна при х = ± 1 . В случае внешнего обтекания тела координата l = ch x может достигать бесконечных значений, а координата m ограничена. Примем во внимание, что потенциал скоростей возмущенного движения (т.е. обтекания за вычетом однородного потока со скоростью, равной скорости на бесконечности) должен стремиться к нулю при удалении от поверхности тела, причем .
Из приведенных соображений следует, что искомые частные решения должны иметь вид произведений Qn ( l ) Pn ( m ) (n = 1, 2,…);
подчеркнем, отсчет n при суммировании начинается с единицы, а не с нуля. Это подтверждается наличием следующих очевидных асимптотических равенств, справедливых при больших значениях l, а, следовательно, согласно (1), и R = = Ö х2 + r2 , имеющего тот же порядок, что иl:
Таким образом, будем иметь правильный порядок убывания j' на бесконечности, если положим
, (6)
где An - постоянные коэффициенты, зависящие от формы поверхности тела.
Складывая потенциалы j¥ и j', получим искомый потенциал скоростей продольного обтекания тела вращения со скоростью на бесконечности, равной U ¥ ,
(7)
Для определения коэффициентов An найдем выражение функции тока y. По формуле (2) будем иметь