Применение системы интегральных уравнений к решению плоских задач теории крыла Текст научной статьи по специальности «Физика»

Вычислительные технологии

Том 4, № 5, 1999

ПРИМЕНЕНИЕ СИСТЕМЫ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ К РЕШЕНИЮ ПЛОСКИХ ЗАДАЧ ТЕОРИИ КРЫЛА *

Д.Н. Горелов, Ю.С. Смолин Институт информационных технологий и прикладной математики СО РАН, Омск, Россия e-mail: gorelov@iitam.omsk.net.ru, smolin@iitam.omsk.net.ru

A system of two integral equations is suggested for the intensities of the vortex layers on the upper and lower sides of the airfoil making it possible to solve two-dimensional problems of stationary and non-stationary flow past the body allowing for the potential hydrodynamic interaction between the airfoil and other bodies and the flow boundaries.

The efficiency of this system of equations if applied instead of the input singular equations of the first and second kind has been estimated.

Введение

Плоские краевые задачи теории крыла сводятся обычно к решению сингулярных интегральных уравнений первого или второго рода. Эти уравнения имеют параметрическую особенность, связанную с толщиной профиля, что может привести к большой погрешности численного решения вблизи его кромок (см., например, [1, 3]).

В работе [2] для решения задач стационарного обтекания профиля безграничным потоком несжимаемой жидкости получена система интегральных уравнений (СИУ), не имеющих указанной параметрической особенности.

В настоящей работе предложенная в [2] система интегральных уравнений обобщена на случай обтекания профиля с учетом возможного гидродинамического взаимодействия профиля с другими телами и границами потока. Путем численного эксперимента оценена эффективность применения системы интегральных уравнений вместо исходных уравнений первого и второго рода.

1. Система интегральных уравнений

Рассмотрим задачу обтекания профиля нестационарным потенциальным потоком идеальной несжимаемой жидкости, в которой могут находиться и перемещаться другие тела, границы раздела сред и иные источники возмущений. При сделанных предположениях поле скоростей жидкости может быть описано комплексной скоростью v(z,t),z = x + iy,

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, грант №96-01-00093.

© Д.Н. Горелов, Ю.С. Смолин, 1999.

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

■ф.О = ^ + Дф,() + 2^/, С е ОД- (1)

Ь

Здесь — комплексная скорость жидкости в бесконечно удаленной точке, в — дуговая координата точки £(в,*), 7(в,*) — интенсивность вихревого слоя, моделирующего контур Ь(*), Д'(г,*) — комплексная скорость возмущенного движения жидкости, создаваемого другими телами, границами потока, а также вихревыми следами (линиями контактного разрыва), сходящими с профиля при изменении циркуляции скорости вокруг него с течением времени *

Пусть теперь точка г е Ь(*). Применяя в (1) к интегралу типа Коши формулы Племе-ля, получим

'(ОД = И (г,*) - 17 (в,*)е-*(в’4), г е ОД, (2)

где интеграл по Ь(*) в (1) понимается в смысле главного значения по Коши, и — комплексная скорость движения точек контура £(*), 0(в,*) — угол между касательной к Ь(*) и осью х. Умножая (2) на ехр{г0(в, *)} и разделяя в полученном выражении вещественную и мнимую части, придем к двум сингулярным интегральным уравнениям относительно функции 7(в,*):

Ме^ИОД - и (г,*)]} = 0, г е ОД, (3)

2 7(в,*) + Кв{егб(5’*) [V (г, *) — и (г,*)]} = 0, г е Ь(*). (4)

Уравнение (3) определяет собой условие непротекания жидкости через контур Ь(*) (равенство нулю нормальной составляющей относительной скорости) и является интегральным уравнением первого рода. Уравнение (4) налагает ограничения на касательную составляющую относительной скорости жидкости вдоль контура Ь(*) и является интегральным уравнением второго рода.

В предельном случае бесконечно тонкого профиля каждое из уравнений (3), (4) принимает одинаковые значения на верхней и нижней сторонах профиля. Эту параметрическую особенность, связанную с толщиной профиля, можно исключить, переходя от независимых уравнений (3), (4) к системе совместных интегральных уравнений специального вида. Следуя работе [2], запишем уравнения (3), (4) на верхней Ь^*) и нижней Ь2(*) сторонах

профиля. Беря разность соответствующих выражений, придем к следующей системе ин-

тегральных уравнений относительно интенсивностей вихревых слоев на верхней 71 (в,*) и нижней 72(в,*) сторонах профиля:

1ш{е^1’*)'ф 1,*) — е^’^'ф,*)} = 0, (5)

2[71(в1,«) — 72(82, *)] + — еЛ’<*’’,)ф2,ед = 0. (6)

Здесь точки г1,г2 должны соответствовать друг другу в том смысле, что в предельном случае бесконечно тонкого профиля |г2 — г11 ^ 0; в1,в2 — дуговые координаты точек г1,г2; 71,72 — значения 7 на Ь1,Ь2. Можно показать, что уравнения (5), (6) не вырождаются при уменьшении толщины профиля до нуля.

Следует отметить, что комплексная скорость Д'(х,*), входящая в выражение (1) для г;(х,*), в общем случае является неизвестной функцией, для определения которой требуются дополнительные соотношения. Ими могут быть, в частности, граничные условия на других телах, вихревых следах и границах раздела сред.

2. Решение интегральных уравнений методом панелей

Интегральные уравнения (3), (4) и систему уравнений (5), (6) будем решать методом панелей с линейным распределением интенсивности вихревых слоев. Выберем на контурах

г Г (!) (2) П ЛТ (!) (2) (!) (2)

Ьі,ь2 точки Хт , хт ,т = , полагая х0 = х0 = хА, х^ = = хв, где хА,хв —

комплексные координаты передней и задней кромок профиля. Следуя методу панелей,

(г) (г)

заменим контуры Ь1,Ь2 панелями КГ3-, соединяющими точки х3_1,х3 Є Ьг; і = ,

(г)

г = 1, 2. Обозначим через значения интенсивности вихревого слоя на стыках панелей

(в точках х3г)). Комплексная скорость г(х,і), индуцированная панелями, вычисляется по формулам

г(х, і) = + Дг(х, і) + ^ ^ (і) \ Р)Г> (х) + 3 (х)

г=1 3=0

где р№ (х) = о, д0г)

5(г)/

і(г)/

(7)

0,

Р}Г)(Х)

Є гвг^+\

2пі

1 +

X — X.

(г)

'3+1

х(г) — х(г) х3+1 х3

1п

А1 ( X — X

(г)

3+1

X —

(г)

х3)

«Г(Х) = —

_ гй ■

Є гйг^

2пі

1 +

X — X.

(г)

3-1

(г) Лг)

1п

А2 ( X — X

(г)

X

— X

3-1

X — X

(г)

3-1

і = 0,...,Ж — 1,

N.

(8)

Здесь А1 = —1, если х Є Кг3+1, иначе А1 = 1; А2 = —1, если х Є КГ3-, иначе А2 = 1.

Лр)

Контрольные точки х0к) на панелях Крк в общем случае будем размещать на расстоянии

у. Дкр) от начала панели:

х0к = 4-1(1 — Ук) + хк Ч, 0 < Ук < 1; Дк) =

ІР)

Лр~)

хк — хк-1

, к = 1,..., Ж, р = 1, 2. (9)

Подставляя (7), (8) в уравнения (3), (4) и (5), (6) и требуя выполнения этих уравнений в контрольных точках (9), получим соответственно следующие системы из 2Ж линейных алгебраических уравнений с вещественными коэффициентами относительно 2(Ж + 1) ве-

личин 70

(1)

(1) (2) >тУ ,т0 ,

2 N

(2)

■ тУ:

(х0кМ = —М с(х0*) к р = 1,2;

(10)

г=1 3=0

2 N

(1 — укЬк_1 + Укткр)] +^5]т3г)4г) (х0?) = —Ке {с (х0?) }, к= ^..., ^ р =1,2;

2 N

ЕЕтГ

г=1 3=0

г=1 3=0

Ь(г) ( х(1М Ь(г) Г х(2) °3 I х0к ) — °3 I х0к

—1т <1 с (х0кМ — с(х0Л \ , к = 1,... , N; (12)

(2)

(11)

1

і

2 N

+ 7-

(г)

3

г=1 3=0

а

1 2

(-) Л(1)

(1 — ук) (тк-1 — тГЛ) + ук (тГ — 7к‘

(2)

(1)

(2)

+

(х0?) — 3 (х0?)] = —Ке {с (х0?) — с (х0?)}, к = ^..., Ж. (13)

Здесь

(г)

а3

Др)

0

+ 3 (х0р?) = (—1)Р+Уйрк с(х0ї) = (— 1)р+1егйр‘ [ггте + Дгг (х0Рр))

Р(

(гМ ^(р) 0к

(г)

3

У(Р)

"0к

Єгй рк

1р) Лр)

-к

— х(Р)И /Д

4Р)-

К каждой системе (10), (11) и (12), (13) следует добавить дополнительные уравнения, связанные с решением конкретной задачи. В частности, при решении задачи стационарного обтекания профиля такими соотношениями являются

(1) (2) То — т0

0, 7м) + 7^ =

(14)

Первое из них определяет собой условие непрерывности вихревого слоя в передней кромке профиля, а второе следует из постулата Кутты — Жуковского.

3. Результаты численного эксперимента

Перейдем к сравнительной оценке эффективности применения интегральных уравнений

(3), (4) и системы уравнений (5), (6) к решению плоских задач теории крыла. В качестве тестовой выберем задачу стационарного обтекания профилей Жуковского, для которой известно точное решение [2]. Численное решение этой задачи с помощью интегрального уравнения первого рода приводит к системе алгебраических уравнений (10), (14), с помощью интегрального уравнения второго рода — к СЛАУ (11), (14), а с помощью системы интегральных уравнений (5), (6) — к СЛАУ (12) — (14).

Предварительные расчеты показали, что выбор контрольных точек в середине панелей Ч = 1/2, к = 1,... , N) дает хорошие результаты для тех систем алгебраических уравнений, которые соответствуют интегральному уравнению первого рода. Для решения же систем уравнений, соответствующих интегральному уравнению второго рода, контрольные точки следует выбирать не в середине, а на 3/4 длины каждой панели.

Численный эксперимент проводился для несимметричных профилей Жуковского с относительной толщиной 12% и 0,5%. Все СЛАУ для каждого профиля решались с одинаковым выбором панелей и контрольных точек. Число разбиений N варьировалось от

(г)

10 до 80. При N > 20 относительная погрешность расчета величин 7З практически не менялась.

На рис. 1 представлены результаты расчета распределения по профилю квадрата безразмерной скорости жидкости (и/и^)2 для угла атаки 10°, N = 40 по трем вариантам, соответствующим применению сингулярных уравнений первого рода (а), второго рода (б) и предлагаемой системы СИУ (в). Относительная погрешность расчета скорости жидкости вдоль верхней стороны профиля для тех же вариантов показана на рис. 2.

Приведенные результаты показывают, что точность решения задачи обтекания профиля с помощью интегральных уравнений первого и второго рода существенно зависит от толщины профиля. А именно, с уменьшением толщины профиля растет погрешность расчета. При этом для уравнения первого рода максимальная погрешность имеет место в

Рис. 1. Распределение квадрата безразмерной скорости вдоль контура профилей Жуковского. Сплошная линия — точное решение, точки — приближенное решение СИУ первого рода (а), СИУ второго рода (б) и системы СИУ (в).

5

И, %

1 1 ■Л 101 и, % 5 - П Л * б \ / \ / N

1 </> в 1 0 - I ^ 6 ч . N 1 \ | \ 1 л

;/6

с/Ь

Рис. 2. Относительная погрешность расчета распределения скорости жидкости на верхней стороне профиля для вариантов а, б, в.

окрестности задней кромки, где толщина профиля наименьшая, а для уравнения второго рода — вблизи передней кромки, где радиус кривизны наименьший. Что касается расчета с помощью предлагаемой системы интегральных уравнений, то он сохраняет высокую точность для всех толщин профиля.

Список литературы

[1] БЕлоцЕрковский С. М., Котовский В.Н., Ништ М. И., Федоров Р. М. Математическое моделирование плоскопараллельного отрывного обтекания тел. Наука, М., 1988.

[2] Горелов Д. Н. Об интегральных уравнениях задачи обтекания профиля. Изв. РАН, МЖГ, №2, 1992, 173-177.

[3] РявчЕнко В. П., Сарен В. Э. К расчету аэродинамических характеристик решеток профилей произвольной формы. Изв. АН СССР, МЖГ, №2, 1972, 105-112.

[4] YUKMELEN M. A., Erim M. Z. A general iterative method to design Karman —Trefftz and Joukowsky airfoils. Num. Meth. in Engineering, Int. J., 20, No. 7, 1984, 1349-1368.

Поступила в редакцию 17 ноября 1998 г., в переработанном виде 6 января 1999 г.