Курсовая работа: Продольное и поперечное обтекание тел вращения

Определив коэффициенты A_n и C_n , найду выражения потенциалов и компонент скоростей для продольного и поперечного обтеканий, после чего уже нетрудно разыскать и распределение скоростей и давлений по поверхности заданного тела вращения или вне его при любом угле. Отмечу, что при всех вычислениях на поверхности удлиненного тела и вблизи нее можно пользоваться для Q_n и dQ_n /d l приближенными выражениями (18). Само собой разумеется, что при удалении от поверхности обтекаемого тела l возрастает и формулы (18) становятся все менее и менее точными.

4. Применение метода особенностей для расчета продольного и поперечного обтеканий тел вращения

Изложенный в предыдущих параграфах (§ 1 и § 2) метод исследования продольного и поперечного обтеканий тел вращения, основанный на непосредственном решении уравнения Лапласа в эллиптических координатах, не является единственным методом решения этой задачи. Первоначально формы обтекаемых тел вращения для дирижаблей определялись наложением однородного, параллельного некоторой оси потока на поток от системы источников (стоков), распределенных вдоль той же оси. Для этой цели применялись вначале дискретные особенности потока – системы источников (стоков) или диполей, а впоследствии – непрерывные их распределения.

Предположим для определенности, что на отрезке (– с, + с) оси Ох задано непрерывное распределение источников (стоков) интенсивности q(х). Тогда потенциал j' возмущенного движения, созданного этой системой особенностей, будет равен (знак минус введем в определение интенсивности q )

(24)

Если задаться видом функции q(x') , то, вычисляя интеграл (24), получим потенциал скоростей, а дифференцирование по r и x позволит вычислить и проекции скорости V_r и V_x . Наоборот, задаваясь формой обтекаемого тела, можно, переходя от потенциала скоростей возмущенного движения к полному потенциалу продольного обтекания тела однородным потоком с заданной скоростью на бесконечности и написав условие непроницаемости поверхности тела, получить интегральное уравнение, в котором q(х') будет неизвестной функцией. Заменяя потенциал скоростей на функцию тока, Карман[6] разработал метод приближенного интегрирования соответствующего интегрального уравнения, основанный на замене интеграла конечной суммой. Однако метод Кармана не был достаточно общим и, кроме того, требовал решения в каждом отдельном случае системы большого числа линейных алгебраических уравнений, что делало его на практике слишком трудоемким.

Аналогично, пользуясь выражением потенциала диполя: – m cos q / (4 p r²⁾ , можно составить и потенциал j¢ поперечного обтекания тела вращения, складывая потенциал однородного натекания с заданной скоростью на бесконечности с потенциалом скоростей возмущенного движения жидкости от непрерывно распределенных по отрезку – с < х < с диполей интенсивности m(х')

(25)

Здесь также можно задаваться распределением интенсивности m (х') или, наоборот, определять эту интенсивность из интегрального уравнения, представляющего условие непроницаемости заданной поверхности тела по отношению к потоку, складывающемуся из возмущенного и однородного на бесконечности.

Не останавливаясь на изложении этих в настоящее время уже малоупотребительных методов, укажем лишь на простую их связь с методами, изложенными в предыдущих параграфах. Покажем, что при заданной форме поверхностей обтекаемых тел вращения неизвестные функции q(х') и m(х') могут быть выражены через ранее введенные коэффициенты A_n и С_n .

Разобьем ось симметрии тела вращения Ох на две области: одну, определяемую интервалом – с £ х £ с , заполненным особенностями, и вторую, представляющую остальную часть оси Ох , где | x | > c. В эллиптических координатах l, m отрезок, на котором расположены особенности, можно представить, согласно формуле r = c Ö l ² – 1 Ö 1 – m ² , так:

l = 1, – 1 £ m £ 1,

а остальную часть оси Ох , как

m = ±1, 1 < l < ¥ .

Тогда, сравнивая между собой вне отрезка (– с < х' < с) выражения потенциалов возмущений (24) и (25) с соответственными выражениями тех же потенциалов и приняв во внимание, что Р_n (1) = 1 , получим следующие два равенства:

(26)

(27)

которые при заданных коэффициентах A_n и С_n можно рассматривать как два интегральных уравнения для определения неизвестных функций q и m.

Интегральное уравнение (26) может быть решено, если искать решение в виде ряда , – 1 £ m' £ 1. Подставляя это разложение в (26), получим

Замечая, что по известной формуле теории функций Лежандра[7]

перепишем предыдущее интегральное уравнение в виде

откуда будет сразу следовать искомое решение

a_n = 2 p cU _¥ A_n , (28)

Для разыскания второй неизвестной функции m(х') продифференцируем раз по l и другой раз по m' известное разложение[8]

тогда получим