Метод расчета нестационарного обтекания произвольного профиля Текст научной статьи по специальности «Физика»

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ОБРАБОТКА ДАННЫХ

УДК 51-7

DOI: 10.18101/2304-5728-2019-4-65-72

МЕТОД РАСЧЕТА НЕСТАЦИОНАРНОГО ОБТЕКАНИЯ ПРОИЗВОЛЬНОГО ПРОФИЛЯ

© Пурбо Ламажапович Абидуев

кандидат физико-математических наук, доцент,

Бурятская государственная сельскохозяйственная академия

имени В. Р. Филиппова

Россия, 670024, г. Улан-Удэ, ул. Пушкина, 8

E-mail: apl087@yandex.ru

© Дармаев Тумэн Гомбоцыренович

кандидат физико-математических наук, доцент, Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: dtg@bsu.ru

© Дамбаев Жаргал Гомбоевич

доктор технических наук, профессор,

Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: g.dambaev@rambler.ru

В данной работе рассматривается нестационарное безотрывное обтекание плоским потенциальным потоком несжимаемой жидкости произвольного профиля, колеблющегося как твердое тело с малой амплитудой по некоторому гармоническому закону. Для устранения особенности ядра интегрального уравнения Фредгольма второго рода относительно величины нестационарной скорости в точке задней кромки профиля используется постулат Кутта — Жуковского. Гидродинамический смысл проведенного устранения особенности ядра состоит в том, что решение ищется в классе функций, обеспечивающих непрерывность давления всюду в потоке, в том числе и в острой задней кромке. Критическая линия тока стационарного потока, необходимая для вычисления несобственных интегралов, определяется как решение задачи Коши. Заменяя интегралы в уравнениях конечными суммами, получаем системы из N алгебраических уравнений относительно искомых величин. Ключевые слова: произвольный профиль; нестационарное обтекание; идеальная несжимаемая жидкость; интегральное уравнение Фредгольма второго рода.

Для цитирования:

Абидуев П. Л., Дармаев Т. Г., Дамбаев Ж. Г. Метод расчета нестационарного обтекания произвольного профиля // Вестник Бурятского государственного университета. Математика, информатика. 2019. № 4. С. 65-72.

Введение

В данной работе рассматривается нестационарное безотрывное обтекание плоским потенциальным потоком несжимаемой жидкости профиля, колеблющегося как твердое тело с малой амплитудой по некоторому гармоническому закону. По данной проблеме опубликовано большое число специальных монографий и журнальных статей, но практически все имеющиеся результаты получены для слабоизогнутых и тонких профилей [1-10]. В данной же работе предлагается метод расчета нестационарного обтекания произвольного или телесного профиля, приведенного ранее авторами в работе [11] к интегральным уравнениям Фредгольма второго рода относительно величины стационарной скорости у0(£ ) и амплитудного значения нестационарной части относительной скорости у(я).

1 Устранение особенности ядра

Введем прямоугольную систему координат ОХУ, так, что ось ОХ совпадает с геометрической хордой профиля, в среднем положении, причем х=0 соответствует носику, а координата х=1 — хвостику профиля. Обозначим через Ц — контур профиля в момент времени 1;, через Ь — контур профиля в среднем положении и, соответственно, Сг и С — вихревые следы. Предполагаем, что в бесконечном удалении перед профилем жидкость имеет постоянную скорость У¥ .

Рассмотрим полученное ранее авторами [11] интегральное уравнение Фредгольма второго рода относительно величины нестационарной скорости у(я):

1 у(5) - —V(а)Я^а^а (у,(я),я) (1)

2 2—1

¥

Я( я,а) = К (я, а) - | е~ И"К (я, и )с1и ,

К (я,а) = - Яе,

0

ia(а)

2 (я) -£(а)

(х(я) - £(а))8та(я) - (у(я) - г](а))оо8а(я)

(х(я)-С(а))2 + (у (я)-ц(а))

1 "

^(у0(я),я) = Яе, Ф(я) +—q2Д х |е(я,и)с1и ,

2— п

Ф(я) = е

,а( 5 )

( ^ . г,1^\„а(а)

щг ааа __!_ г у^g(») -«(а) „„.

2—1 2 -д 2— Ц (2 - д)

у0(я) - у ^(я) 0 А = Яе, |g(а)еа(а)аа = | ^^х + gydy)

где 1 — мнимая единица, ] — мнимая единица, связанная с временными процессами, а(я) — угол наклона к оси ОХ касательной к контуру Ь,

аДя, t) = а(я) + Р(я)в3'м — угол между касательной к контуру Ц и осью ОХ, д = — число Струхаля, ю — круговая частота колебаний,

g = gx - igy — комплексная форма гармонических колебаний профиля.

Интегрирование по контуру ведется в положительном направлении (против часовой стрелки). Все линейные размеры нормированы к длине хорды профиля Ь, а величины, имеющие размерность скорости, — к скорости набегающего потока У¥ .

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

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

Интегральное уравнение (1) с данным ядром справедливо для всех точек контура Ь за исключением задней кромки В, в которой 8=о=о и ядро имеет полюс первого порядка. Для однозначного определения решения уравнения (1) используется известный постулат Кутта — Жуковского, который формулируется так: скачок касательной составляющей скорости жидкости в задней кромке равен интенсивности сходящего с нее вихря:

где ^(0) — скорость стационарного потока в точке В, у(+о) и у(-о) — предельные значения относительной скорости жидкости при подходе к точке В снизу и сверху, соответственно.

,2

(2)

Условие (2) используется для устранения особенности ядра уравнения

(1) в задней кромке. Из (1) следует, что:

у(±0) = — Гу(а)R(±0, ст)dст + 2F(±0).

- L

Тогда уравнение (1) с учетом (2) и (3) запишется в виде:

(3)

v(s) - 2-Гу(а) I R(s, а) - 2 R(+0, а) + R(-0, а) +

\\

^(0)

dа =

1

(

2 Л А

(4)

= F (s) — F (+0) + F (-0) -- А1 21 2^,(0)

/

Ядро полученного уравнения (4) не имеет особенностей на контуре Ь, и его решение совпадает с решением (1).

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

2 Метод расчета

Интегральное уравнение (4) разделим на действительную и мнимую по j части. Для этого запишем его в операторном виде:

(А + ]В№) = а + ]Ь, = ^(з) + р2^), (5)

где

I Г|*(*„) - *(+0'а) + К(-"'а) -

-ту » 1 1

А = Е--|

-

Г 8т qu I К(s, и)

К (+0, и) + К (-0, и)

Л

dа,

(6)

в=~— Ц !о 008 qu (* и) -

- L N

а = 2Яе, 1 (

К (+0, и) + К (-0, и)

—

( е'а( з) 2—

| У0(а)

g (з) - g (а) (г -V)2

dа

2у0(0)у q 2 А м

-- Г cos(qu) К (з, и^и

dа, (7)

а(+0) + а(-0) -

<2 А1 Л 2у0(0) ^

(8)

Ь = д Яе

РI ^^^ - я!5) ^

г - д

/0

^I¿ш(ди)К(5,и)ёи - Ь(+0) + Ь(-0).

р а 2

(9)

Критическая линия тока стационарного потока, необходимая для вычисления несобственных интегралов в формулах (6)-(9), определяется как решение следующей задачи Коши:

ёх

— = ^ х(х у\

(10)

— = V) У(х У> х(0) = хв = 1, у(0) = ув = 0.

Здесь (хв, ув) — координата выходной кромки профиля в среднем положении. Функции У0х(х,у), У0у (х,у) определяются по формуле:

= I у0(°)аа

0 х 0 у

2я1{г (г) -£(ст)

(11)

(12)

где г(г)= (х(г)-£(ст)) + /(у(г)г — длина дуги контура, отсчитываемая вдоль линии вихревого следа от задней кромки, 0 — геометрический угол атаки, образуемый вектором скорости потока на

бесконечности V« и осью ОХ.

Разделяя уравнение (5) на действительную и мнимую части, получим систему двух уравнений:

Г АуДя) - Ву2(я ) = а, I ВУ1(5) + Ау2(5) = Ь. Из первого уравнения находим:

у1(5) = А_1Ву2(5) + А'1 а. (13)

Подставляя (13) во второе уравнение (12) и решая полученное уравнение относительно ^(я), имеем:

У2(5) = (ВА-'В + А)(Ь - ВА-'а). (14)

Откуда, подставляя в (13),

v1 (5) = А-1 В(ВА-в + А)(Ь - ВА-'а) + А'1 а. (15)

Далее заменой интегралов в уравнениях (6)-(9), (11) конечными суммами методом, предложенным в работах [12-14], получим системы из N алгебраических уравнений относительно искомых величин, которые можно решить с применением ЭВМ.

Заключение

Интегральное уравнение Фредгольма второго рода относительно величины нестационарной скорости v(s) справедливо для всех точек контура L за исключением задней кромки В, в которой 8=о=о и ядро имеет полюс первого порядка. Для устранения особенности ядра уравнения используется постулат Кутта — Жуковского. Гидродинамический смысл проведенного устранения особенности ядра состоит в том, что решение ищется в классе функций, обеспечивающих непрерывность давления всюду в потоке, в том числе и в острой задней кромке. Критическая линия тока стационарного потока, необходимая для вычисления несобственных интегралов, определяется как решение задачи Коши (10). Заменяя интегралы в уравнениях (6)-(9), (11) конечными суммами, получаем системы из N алгебраических уравнений относительно искомых величин, которые предполагается в дальнейшем решить с применением ЭВМ.

Литература

1. Prandtl L. Uber die Enstehung von Wirbeln in einer idealen Flüssigkeit // Vortrage zur Hydro-und Aerodynamik. Berlin, 1924. P. 18-33.

2. Чаплыгин С. А. О влиянии плоскопараллельного потока воздуха на движущее в нем цилиндрическое крыло // Труды ЦАГИ. 1926. Вып. 19. 74 с.

3. Келдыш М. В, Лаврентьев М. А. К теории колеблющегося крыла // Технические заметки ЦАГИ. 1935. № 45. С. 48-52.

4. Седов Л. И. К теории неустановившихся движении внутри жидкости // Труды ЦАГИ. 1935. Вып. 229. 40 с.

5. Седов Л. И. Теория плоских движений идеальной жидкости. М.: ГТТИ, 1939. 144 с.

6. Поляков Н. Ф. Теория нестационарных движений несущей поверхности. Л., 1960. 84 с.

7. Горелов Д. Н. Теория крыла в нестационарном потоке. Новосибирск, 1975. 152 с.

8. Горелов Д. Н., Куляев Р. Л. Нелинейная задача о нестационарном обтекании профиля несжимаемой жидкостью // Известия АН СССР. МЖГ. 1971. № 6. С. 38-47.

9. Головкин В. А. Нелинейная задача о неустановившемся обтекании произвольного профиля со свободно деформирующимся вихревым следом // Ученые записки ЦАГИ. 1972. Т. 3, № 3. С. 1-11.

10. Павловец Г.А. Методы расчета обтекания сечения крыла идеальным несжимаемым потоком // Труды ЦАГИ. 1971. Вып. 1344. 72 с.

11. Дармаев Т. Г., Дамбаев Ж. Г. О нестационарном обтекании произвольного профиля // Вестник Бурятского государственного университета. Математика, информатика. 2013. № 2. С. 60-69.

12. Дармаев Т. Г., Цыдыпов Б. Д. Метод расчета стационарного обтекания произвольного профиля // Вестник Бурятского государственного университета. Математика, информатика. 2016. № 3. С. 51-56.

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

14. Рябченко В. П. Нестационарные аэродинамические характеристики решеток произвольных профилей, вибрирующих в потенциальном потоке несжимаемой жидкости // Известия АН СССР. МЖГ. 1974. № 1. С. 15-20.

A METHOD FOR CALCULATING NON-STATIONARY FLOW OF AN ARBITRARY PROFILE

Purbo L. Abiduev

Cand. Sci. (Phys. and Math.), A/Prof., Fillipov Buryat State Academy of Agriculture 8 Pushkina St., Ulan-Ude 670024, Russia E-mail: apl087@yandex.ru

Tumen G. Darmaev Cand. Sci. (Phys. and Math.), A/Prof., Dorzhi Banzarov Buryat State University 24a Smolina St., Ulan-Ude 670000, Russia E-mail: dtg@bsu.ru

Zhargal G. Dambaev

Dr. Sci. (Engineering), Prof.,

Dorzhi Banzarov Buryat State University

24a Smolina St., Ulan-Ude 670000, Russia

E-mail: g.dambaev@rambler.ru

The article deals with non-stationary continuous flow of an arbitrary profile fluctuating as a solid with small amplitude under some harmonic law by potential flow of incompressible liquid. We have used the Kutta-Zhukovsky postulate to eliminate the features of a kernel of Fredholm integral equation of the second kind with respect to the value of non-stationary speed at the point of profile trailing edge. The hydrodynamic meaning of the carried-out elimination of the features of a kernel lies in the fact that the decision is sought in the class of functions that ensure the continuity of pressure everywhere in the flow, including the sharp trailing edge. The critical streamline of a stationary flow necessary for calculating improper integrals is defined as a solution to the Cauchy problem. Replacing the integrals in the equations with finite sums, we obtain systems of N algebraic equations with respect to the sought quantities.

Keywords: arbitrary profile; unsteady flow; ideal incompressible liquid; Fredholm integral equations of second kind.

References

1. Prandtl L. Uber die Enstehung von Wirbeln in einer idealen Flussigkeit. Vortrage zur Hydro-undAerodynamik. Berlin, 1924. Pp. 18-33.

2. Chaplygin S. A. O vliyanii ploskoparallelnogo potoka vozdukha na dviz-hushchee v nem tsilindricheskoe krylo [About the Influence of Plane-Parallel Air Flow on the Cylindrical Wing Moving in It]. Trudy TsAGI. 1926. Vol. 19. 74 p.

3. Keldysh M. V, Lavrentyev M. A. K teorii koleblyushchegosya kryla TOn the Theory of Oscillating Wing]. Tekhnicheskie zametki TsAGI. 1935. No. 45. Pp. 48-52.

4. Sedov L. I. K teorii neustanovivshikhsya dvizhenii vnutri zhidkosti [On the Theory of Unsteady Motions in Liquidl. Trudy TsAGI. 1935. Vol. 229. 40 p.

5. Sedov L. I. Teoriya ploskikh dvizhenii idealnoi zhidkosti [Theory of Planar Motions of an Ideal Fluidl. GTTI, 1939. 144 p.

6. Polyakov N. Ph. Teoriya nestatsionarnykh dvizhenii nesuschei poverkhnosti [Theory of Non-Stationary Motions of a Carrying Surface]. Leningrad, 1960. 84 p.

7. Gorelov D. N. Teoriya kryla v nestatsionarnom potoke [Wing Theory for Unsteady Flow]. Novosibirsk, 1975. 152 p.

8. Gorelov D. N., Kulyaev R. L. Nelineinaya zadacha o nestatsionarnom obte-kanii profilya neszhimaemoi zhidkostiyu [A Nonlinear Problem of Unsteady Flow of a Profile with an Ideal Liquid]. Izvestiya AN SSSR. Mekhanika zhidkosti i gaza. 1971. No. 6. Pp. 38-47.

9. Golovkin V. A. Nelinaya zadacha o neustanovivshemsya obtekanii proizvol-nogo profilya so svobodno deformiruyuschimsya vikhrevym sledom [A Nonlinear Problem of Unsteady Flow of An Arbitrary Profile with a Freely Deformed Vortex Sheet]. TsAGI Science Journal. 1972. Vol. 3. No. 3. Pp. 1-11.

10. Pavlovets G. A. Metody rascheta obtekanya secheniya kryla idealnym nesz-himaemym potokom [Methods for Calculating the Flow of Wing Section with an Ideal Incompressible Flow]. Trudy TsAGI. 1971. Vol. 1344. 72 p.

11. Darmaev T. G., Dambaev Zh. O nestatsionarnom obtekanii proizvolnogo profilya [About the Unsteady Flow of an Arbitrary Profile]. Vestnik Buryatskogo gosu-darstvennogo universiteta. Matematika, informatika. 2013. No. 2. Pp. 60-69.

12. Darmaev T. G., Tsydypov B. D. Metod rascheta statsionarnogo obtekaniya proizvolnogo profilya [A Method for Calculating the Stationary Flow of an Arbitrary Profile]. Vestnik Buryatskogo gosudarstvennogo universiteta. Matematika, informatika. 2016. Vol. 3. Pp. 51-56.

13. Ryabchenko V. P., Saryan V. E. K raschetu aerodinamicheskikh kharakteristik reshetok profilei proizvolnoi formy [On Calculation of Aerodynamic Characteristics of the Lattices of Arbitrarily-Shaped Profiles]. Izvestiya AN SSSR. Mekhanika zhidkosti i gaza. 1972. No. 2. Pp.105-112.

14. Ryabchenko V. P. Nestatsionarnye aerodynamicheskie kharakteristiki reshe-tok proizvolnykh profilei, vibriruyushikh v potentsialnom potoke neszhimaemoi zhid-kosti [Non-Stationary Aerodynamic Characteristics of Lattices of Arbitrary Profiles Vibrating in a Potential Flow of Incompressible Liquid]. Izvestiya AN SSSR. Mekhanika zhidkosti i gaza. 1974. No. 1. Pp. 15-20.