CC BY
7
УДК 519.64
doi: 10.18101/2304-5728-2016-3-51-56
© Т. Г. Дармаев, Б. Д. Цыдыпов
Метод расчета стационарного обтекания произвольного профиля
В данной работе для численного решения интегрального уравнения Фредгольма второго рода относительно величины стационарной скорости предложен метод замены интегрального уравнения системой алгебраических уравнений, основанный на замене профиля вписанным N-угольником с длиной стороны порядка величины Nи внутренним углом, близким к к. На основе предложенного метода созданы алгоритм решения и программа на ФОРТРАН. Проведены тестовые численные расчеты для эллипсов разных толщин и серии профилей «В» ЦАГИ.
Ключевые слова: произвольный профиль, стационарное обтекание, интегральное уравнение Фредгольма второго рода.
О Т. G. Darmaev, В. D. Tsydypov Method of calculation of a stationary flow of any profile
In this work for the numerical solution of the integrated equation of Fredholm of the second sort concerning the size of stationary speed the method of replacement of the integrated equation with system of the algebraic equations based on replacement of a profile with the entered N-square with a length of the party of an order of size and the internal corner close to к is offered. On the basis of the offered method the algorithm of the decision and the program on the FORTRAN are created. Test numerical calculations for ellipses of different thickness and a series of profiles "B" of TsAGI are carried out.
Keywords: any profile, stationary flow, integrated equation of Fredholm of the second sort.
Введение
Задача нестационарного обтекания произвольного профиля, колеблющегося в потоке идеальной несжимаемой жидкости в работе [1] приведена к интегральным уравнениям Фредгольма второго рода относительно величины стационарной скорости v0 (5) и амплитудного значения нестационарной части относительной скорости v(s). Решение полученных интегральных уравнений позволяет свести к квадратурам определение нестационарных характеристик потока и гидродинамических реакций на профиле. Для решения этих уравнений можно использовать метод итераций, примененный в работах [2-3]. Этим методом может быть получен результат с наперед заданной точностью для широкого класса
профилей. Но, как указывается в работе [4], для некоторых несимметричных профилей с большой величиной прогиба средней линии, а также для профилей с точкой возврата (профиль Жуковского с очень тонкой хвостовой частью, например) итерационный процесс расходится. Итерационный процесс также расходится в случае расчета тонких симметричных профилей (тонкая пластина) с профилированной головной и хвостовой частями.
В данной работе для вычисления стационарной скорости у0 (5) применяется метод замены интегрального уравнения системой алгебраических уравнений, предложенный в работах [5-6].
1. Численный метод
Метод основан на замене профиля вписанным Ы-угольником с длиной стороны порядка величины Nи внутренним углом, близким к к. Контур профиля Ь разбивается на N участков 1к длины с!1к. (к = 1.2.....Л'). Обозначим I)' и <$1 соответственно хорду, стягивающую дужку 1к и длину хорды ¡1, далее я°,(т0 - дуговые координаты точек на ломаной
= иС отсчитываемые в положительном направлении обхода контура Разбиение контура выбирается так, чтобы выходная кромка являлась общей дужек 1хи1ы. Целесообразность такого разбиения оправдывается характером течения в окрестности выходной кромки, который используется далее для устранения особенностей ядер данных уравнений. Заменяя интегралы в уравнениях конечными суммами, получим системы из N алгебраических уравнений относительно искомых величин, которые решаются с применением ЭВМ.
Рассмотрим интегральное уравнение Фредгольма второго рода ([1]) относительно величины стационарной скорости :
1 1 г
TvoOO -— v0(cr)K(s,cr)dcr =cos (a(s)-e) (1)
2 2ж i
(x(s) -¿f(cr))sina(s) - (y(s) - t](<j))cosa(s)
K(s,<j) = - Re -
где / - мнимая единица, а (5) - угол наклона касательной к контуру к оси ОХ, 0 - геометрический угол атаки, образуемый вектором скорости потока на бесконечности V, и осью ОХ. Интегрирование по контуру ведется в положительном направлении (против часовой стрелки). Все линейные размеры нормированы к длине хорды профиля Ь, а величины, имеющие размерность скорости, к скорости набегающего потока Vt . Имеют место следующие равенства:
lim Afr.cr) = Пт (Х(*}" ^ - (y(s)-*l(v))cosa(s) =
{x{s)-¿;{o))2+{y{s)-r,{o))2 2r(s)
где /'(.V) - кривизна контура Ь в точке э,
|соз(а(5)-0)Л = О. (3)
ь
Отсюда следует, что для гладких контуров ядро интегрального уравнения (1) может быть доопределено до непрерывного всюду на контуре Ь по обеим переменным 5 и а.
Пусть 5=0 фиксированная выходная кромка профиля. Из общей теории уравнений Фредгольма второго рода, равенств (2) - (3) и теорем 1 и 2 [5] заключаем, что интегральное уравнение (1) при дополнительном условии
= 0 (4)
имеет единственное решение в классе непрерывно дифференцируемых функций. При этом решение удовлетворяет следующему уравнению:
|уо " Гуо " >Аа = совСоф) - в) - сое(сф0) - в). (5)
2 2п\
Решение этого уравнения совпадает с решением уравнения (1) ([5]). Для искомого решения в любой окрестности точки
существуют точки ^ и sдля которых выполняется равенство:
При достаточно большом числе разбиений N на дужки (малых сИ1 и сП,:) можно считать приближенно /, = /, и /, . = /, ..
Тогда уравнение (5) можно записать в следующей эквивалентной форме:
"0<1)
2 2л
= соз(а(5)-0)-
соз(а(51) -в) + соз(а (%)-$)
ёа =
(6)
Заменяя далее интегралы в (6) конечной суммой, получим систему из N алгебраических уравнений относительно величин У0т, (т = 1.2..... А'):
^ "!>„>„„ =2со5(а(5)-0)-[со5(а(51)-0) + со5(а(^)-0)],
к=1 2
(8)
Я То я
ак) ^тк+1 + ^тк+1 С05(аш ак) аГС^ дус!^ ^ ^
V ^тк+1 ^ тк )
2 Х1к+1тк П
Т^Хтк =х(и0т)-%(о°к), ¥тк=у{э°т)-г1{<у1), - координаты
концов дужек 1к, - расчетные точки, в которых выполняется
уравнение (8), ак - угол, образуемый хордой I)' с положительным направлением оси ОХ.
Многозначная функция агс^, необходимая для вычисления элементов атк нормируется так, чтобы обеспечить ее непрерывное изменение по переменным э, о всюду, кроме 8=0. Нормировка функции
(<т°), хя=х(5°я), ут=у(*°т)
Л &-хт
зависит от знаков числителя и знаменателя, и от взаимного расположения
1 \п — У
точек ак и 5 . Обозначим =—аг^ —^-^<0.5 .
Непосредственной проверкой нетрудно убедиться, что требуемым условиям удовлетворяет функция:
[1, если (4к-хт)(г]к-ут)>0,
[-1, иначе
(1, если %к~хт> 0, где Л
[-1, иначе
Г1, если расположена на спинке профиля,
<Зт=<
[-1, если на корытце.
2. Расчеты
На основе предложенного метода созданы алгоритм решения и программа на ФОРТРАН. Проведены численные расчеты для эллипсов разных толщин и серии профилей «В» ЦАГИ со следующими геометрическими характеристиками:
у = + 7х2 - 6х3) ± х°'87(1- х2)0'56 ^ - уравнение верхнего и ниж-
него контуров (знак + относится к верхнему контуру), х и у - координаты точек профиля. На рисунке в нижней части изображены расчетные точки верхнего контура, а в верхней части - распределение стационарной скорости для профиля «В» ЦАГИ при с = 0.2.
Заключение
В данной работе для численного решения интегрального уравнения Фредгольма второго рода относительно величины стационарной скорости у0 (5) предложен метод замены интегрального уравнения системой алгебраических уравнений, основанный на замене профиля вписанным 14-угольником с длиной стороны порядка величины N~1 и внутренним углом, близким к к. На основе предложенного метода созданы алгоритм решения и программа на ФОРТРАН. Проведены тестовые численные расчеты для эллипсов разных толщин и серии профилей «В» ЦАГИ.
Литература
1. Дармаев Т. Г., Дамбаев Ж. Г. О нестационарном обтекании произвольного профиля // Вестник Бурятского государственного университета. Математика, информатика. — 2013. — № 2. — С.60 - 69.
2. Головкин В. А. Нелинейная задача о неустановившемся обтекании произвольного профиля со свободно деформирующимся вихревым следом // Учёные записки ЦАГИ. — 1972. — Т.З. — №3.
3. Шумский Г. М. Нелинейная задача о движении системы профилей вблизи волнистой стенки // Известия СО АН СССР. Сер. техн. наук. — 1982. — №3. — Вып. 1. — С. 24 - 27.
4. Воробьев Н. Ф., Шашкина Г. Н. Метод расчета обтекания профилей несжимаемым потоком. — Исследование обтекания тел численными методами. — Новосибирск ИТПМ СО АН СССР. —1976. — С. 66 - 76.
5. Рябченко В. П., Сарян В. Э. К расчёту аэродинамических характеристик решёток профилей произвольной формы // Известия АН СССР. МЖГ. — 1972. — №2. — С. 105 - 112.
6. Рябченко В. П. Нестационарные аэродинамические характеристики
решёток произвольных профилей, вибрирующих в потенциальном потоке несжимаемой жидкости // Известия АН СССР. МЖГ. — 1974. — № 1. — С. 15-20.
References
1. Darmaev Т. G., Dambaev Zh. G. О nestacionarnom obtekanii proiz-vol'nogo profilja // Vestnik Buijatskogo gosudarstvennogo universiteta. Mate-matika, informatika. — 2013. — № 2. — S.60 - 69.
2. Golovkin V. A. Nelinejnaja zadacha о neustanovivshemsja obtekanii pro-izvol'nogo profilja so svobodno deformirujushhimsja vihrevym sledom // Uch-jonye zapiski CAGI. — 1972. — T.3. — №3.
3. Shumskij G. M. Nelinejnaja zadacha о dvizhenii sistemy profilej vblizi volnistoj stenki // Izvestija SO AN SSSR Ser. tehn. nauk. — 1982. — №3. — Vyp. 1, —S.24-27.
4. Vorob'ev N. F., Shashkina G. N. Metod rascheta obtekanija profilej neszhimaemym potokom. — Issledovanie obtekanija tel chislennymi me-todami. — Novosibirsk ITPM SO AN SSSR. —1976. — S. 66 - 76.
5. Rjabchenko V. P., Sarjan V. Je. К raschjotu ajerodinamicheskih harakter-istik reshjotok profilej proizvol'noj formy // Izvestija AN SSSR. MZhG. — 1972. — №2, —S. 105- 112.
6. Rjabchenko V. P. Nestacionarnye ajerodinamicheskie harakteristiki reshjotok proizvol'nyh profilej, vibrirujushhih v potencial'nom potoke neszhi-maemoj zhidkosti // Izvestija AN SSSR. MZhG. — 1974. — № 1. — S. 15 -20.
Дармаев Тумэн Гомбоцыренович, кандидат физико-математических наук, заведующий лабораторией вычислительных и геоинформационных технологий Бурятского государственного университета, e-mail: [email protected].
Цыдыпов Балдандоржо Дашиевич, доктор технических наук, ведущий научный сотрудник Института физического материаловедения СО РАН, e-mail: [email protected].
Darmaev Tumen Gombotsyrenovich, PhD, Head of the laboratory of Buryat State University.
Tsydypov Baldandorzho Dashievich, DSc, Leading researcher, Institute of Physical Materials Science of SB RAS.
