Научная статья на тему 'КРЕЙСЕРСКАЯ СКОРОСТЬ ЦИЛИНДРИЧЕСКОГО КРЫЛА ПРИ МАЛЫХ ПОСТУПАТЕЛЬНО-ВРАЩАТЕЛЬНЫХ КОЛЕБАНИЯХ'

КРЕЙСЕРСКАЯ СКОРОСТЬ ЦИЛИНДРИЧЕСКОГО КРЫЛА ПРИ МАЛЫХ ПОСТУПАТЕЛЬНО-ВРАЩАТЕЛЬНЫХ КОЛЕБАНИЯХ Текст научной статьи по специальности «Физика»

CC BY
33
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАШУЩЕЕ КРЫЛО / ГАРМОНИЧЕСКИЕ КОЛЕБАНИЯ / КРЕЙСЕРСКАЯ СКОРОСТЬ / УРАВНЕНИЕ НАВЬЕ - СТОКСА / АСИМПТОТИЧЕСКИЙ АНАЛИЗ

Аннотация научной статьи по физике, автор научной работы — Егоров Андрей Геннадьевич, Нуриев Артем Наилевич

Работа посвящена исследованию пропульсивного движения машущего крыла кругового поперечного сечения. Решена задача о гармонических поступательно-вращательных колебаниях крыла с произвольным сдвигом фаз в вязкой несжимаемой жидкости, движение которой описывается нестационарным уравнением Навье - Стокса. С применением метода последовательных асимптотических разложений для случая малых амплитуд колебаний построено аналитическое решение задачи в первых двух членах. Показано, что нелинейное взаимодействие временных гармоник поступательных и вращательных колебаний создает вторичные потоки, которые заставляют крыло двигаться в направлении, перпендикулярном к оси поступательных колебаний. Для крейсерского движения, когда средняя гидродинамическая сила, действующая на крыло, равна нулю, получены зависимости средней скорости от безразмерных параметров колебаний.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Егоров Андрей Геннадьевич, Нуриев Артем Наилевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

CRUISING SPEED OF A CYLINDRICAL WING PERFORMING SMALL TRANSLATIONAL-ROTATIONAL OSCILLATIONS

This work considers the propulsive motion of a flapping wing of a circular cross section. The problem of harmonic translational-rotational oscillations of the wing with an arbitrary phase shift in a viscous incompressible fluid, the motion of which is described by the non-stationary Navier-Stokes equation, is handled. An analytical solution of the problem is obtained in the first two terms by using the method of successive asymptotic expansions for the case of small oscillation amplitudes. It is shown that the nonlinear interaction of time harmonics of translational and rotational oscillations creates secondary flows that make the wing to move in the direction perpendicular to the axis of translational oscillations. For the cruising motion regime, when the average hydrodynamic force acting on the wing is equal to zero, the dependences of the average speed on the parameters of dimensionless oscillation are found.

Текст научной работы на тему «КРЕЙСЕРСКАЯ СКОРОСТЬ ЦИЛИНДРИЧЕСКОГО КРЫЛА ПРИ МАЛЫХ ПОСТУПАТЕЛЬНО-ВРАЩАТЕЛЬНЫХ КОЛЕБАНИЯХ»

УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

2022, Т. 164, кн. 2-3 С. 170-180

ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)

ОРИГИНАЛЬНАЯ СТАТЬЯ

УДК 534.121

doi: 10.26907/2541-7746.2022.2-3.170-180

КРЕЙСЕРСКАЯ СКОРОСТЬ ЦИЛИНДРИЧЕСКОГО КРЫЛА ПРИ МАЛЫХ ПОСТУПАТЕЛЬНО-ВРАЩАТЕЛЬНЫХ

КОЛЕБАНИЯХ

А.Г. Егоров, А.Н. Нуриев

Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия

Работа посвящена исследованию пропульсивного движения машущего крыла кругового поперечного сечения. Решена задача о гармонических поступательно-вращательных колебаниях крыла с произвольным сдвигом фаз в вязкой несжимаемой жидкости, движение которой описывается нестационарным уравнением Навье-Стокса. С применением метода последовательных асимптотических разложений для случая малых амплитуд колебаний построено аналитическое решение задачи в первых двух членах. Показано, что нелинейное взаимодействие временных гармоник поступательных и вращательных колебаний создает вторичные потоки, которые заставляют крыло двигаться в направлении, перпендикулярном к оси поступательных колебаний. Для крейсерского движения, когда средняя гидродинамическая сила, действующая на крыло, равна нулю, получены зависимости средней скорости от безразмерных параметров колебаний.

Ключевые слова: машущее крыло, гармонические колебания, крейсерская скорость, уравнение Навье-Стокса, асимптотический анализ

Полет птиц, при котором крыло совершает сложное поступательно-вращательное периодическое движение, издавна привлекал внимание людей. Математическая теория машущего полета начала развиваться примерно сто лет назад. В 1924 г. Прандтль [1] сформулировал проблему нестационарного движения крыла, которая в линейном приближении была решена в работе Бирнбаума [2]. Подход Прандт-ля - Бирнбаума, основанный на введении поверхностей разрыва скоростей внутри идеальной жидкости, нашел свое дальнейшее развитие в работах [3-6] и др. До настоящего времени основой большинства аналитических исследований машущего крыла является модель идеальной жидкости, дополняемая специальными предположениями о характере отрыва течения [1, 2, 4, 6], о структуре вихревого следа [7, 8] и т. п. В настоящем исследовании развивается иной подход, основанный на методе асимптотических разложений решения полной системы уравнений Навье-Стокса по малому параметру [9-13]. В качестве малого параметра выбирается отношение амплитуды А поперечных колебаний к радиусу а цилиндра: 7 = А/а.

Асимптотическое исследование задачи о вращательно-поступательных колебаниях круглого цилиндра впервые было проведено в [14]. В работе главным образом рассматривались структурные особенности вторичных течений в условиях нулевой средней скорости движения цилиндра. Для случая высоких частот авторами было

Аннотация

Введение

отмечено появление ненулевой средней силы, действующей на цилиндр. Этот факт, однако, остался практически незамеченным.

В действительности, как показано в [15], асимптотическое решение задачи открывает широкие возможности в изучении пропульсивного движения цилиндра в режиме машущего крыла. В отличие от [14], в указанной работе аналитическое решение строится для случая движения тела с крейсерской скоростью в условиях нулевой средней силы, действующей за период колебаний. Крейсерская скорость при этом определяется из решения задачи.

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

1. Постановка задачи

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

1 дф дф 1 (д (гу9) дуг

г г дв1 9 дг1 г \ дг дв

запишем определяющие уравнения в виде

л , дт 7 д(ф,т) „ . .„.

г > 1: т = Аф, — -1 =в-1 Ат, (1)

д1 г д(г, в) д(ф, т) дф дт дф дт

д(г,в) дг дв дв дг'

где Д - оператор Лапласа, в и Y - безразмерные частота и амплитуда поперечных колебаний. Эти параметры определяются как

в - a2^ - A

v ' Y а

Далее параметр y будет считаться малым.

Уравнения движения жидкости дополним следующими граничными условиями:

г -1: ф -0, дф - f cost (2)

дг

г ^ ж : ф ~—г sin в cos (t + ф) + ифг cos (в + а). (3)

На границе цилиндра (г - 1) задаются условия прилипания (2). Параметр f -это третий безразмерный параметр задачи. Он определяет отношение амплитуд вращательного (©) и поступательного (A) колебаний

„ а©

f - A.

Рис. 1. Схема синфазных поперечно-крутильных колебаний цилиндра. Пунктирной линией показано крайнее положение цилиндра. А, © - амплитуды поперечных и крутильных колебаний, Н = А/© - расстояние между центром цилиндра О и центром скоростей О\ в выделенном положении

Правая часть условия на бесконечности (3) содержит два слагаемых. Первое определяет однородный осцилляционный поток в направлении оси Ох колебаний цилиндра, - сдвиг фаз между крутильными и поперечными колебаниями. Второе слагаемое в (3) задает стационарный поток, направленный под углом а к оси Оу, он определяет среднюю скорость движения цилиндра в жидкости. Далее рассмотрим случай движения цилиндра с крейсерской скоростью. В условиях малости амплитуд колебания порядок крейсерской скорости, очевидно, не может быть выше чем ), поэтому далее примем, что п*ж = . Величины и а будем счи-

тать неизвестными и определять их в ходе решения задачи из условия равенства нулю средних за период сил, действующих на цилиндр:

(Ш) = {¡у (г)) =0. (4)

Угловыми скобками здесь и далее обозначается среднее за период колебаний, через 1х, ¡у - нормированные на 2пра2иш (р - плотность жидкости) компоненты сил, представимые в терминах завихренности как [15]:

2п

¡у - Их = 2 ^ + ^ + {дт -_ егв(5)

о Г_1

Первое слагаемое в выражении для сил есть сила Крылова - Фруда, появляющаяся как результат перехода в подвижную систему координат. Слагаемые, пропорциональные т и Вт/Вт в (5), представляют вклад в силы и момент касательных напряжений и давления соответственно.

2. Асимптотическое решение

2.1. Разложение по малому параметру. Не делая пока никаких допущений относительно частоты колебаний в, предположим малость амплитуд поперечных колебаний 7 ^ 1 и малость углов поворота © ~ 7 (то есть £ ~ 1). Будем искать решение поставленной задачи в виде ряда по степеням 7, ограничиваясь первыми двумя членами

ф = ф(0) + ГФ(1) + ■■■ , т = т(0) + + ■■■ . (6)

В главном члене получаем стоксово приближение исходной задачи

д_ at

-w(0) — в-1Ды(0) -0, w(0) - Дф(0)

дф(0) (7) г - 1: ф(0) -0, —-— - f cos t, г ^ж : ф(0) — — cos (t + ф) ■ г sin в. дг

Стоксово решение задачи описывает чисто периодическое движение жидкости. Для следующего члена имеем из (1)

ф,<1> — в-1Д»<1> - С, w(n - Дф<1>, G - IМ^;

т г д(г,в) (g)

г - 1: ф(1) - д^.— - о, г ^ж : ф(1) — Си^г cos (в + а). дг

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

2.2. Стоксово приближение. В стоксовом приближении линейная задача (7) расщепляется на две. Первая из них - задача обтекания неподвижного цилиндра колебательным потоком жидкости, вторая - задача вращательных колебаний цилиндра в неподвижной жидкости. Таким образом, общий поток представляется в виде суммы двух потоков

- sinвКе(Ф1(г)ег(4+у)) + £Ке(Ф2(г)вгг:), w(0) - sin в Яе^ффф^^) + Сяе^ффф*).

Решение соответствующих одномерных задач дает

*, - —г + i + ^ + .-2w„ W - -^фф!-, SM- Щфф,

г г K0 (z) zK0 (z)

Ф2 - Sg + z-2W2, W2 - —zK-^T, Sg(z)- Kt), z - у/ф.

K (z) zK1 (z)

Здесь Kj - модифицированные функции Бесселя. Как видно, вне стоксова пограничного слоя толщины порядка в-1/2 завихренность экспоненциально быстро затухает и течение становится потенциальным. Дополнительные вязкие члены Sxf-1 и Sg, исчезающие в эффективно невязкой жидкости (при в - ж), корректируют вне пограничного слоя бесциркуляционный поток идеальной жидкости.

В практически важном случае больших в поведение Sx , Sg дается асимптотиками

113 3

Sg(z) - 1 — + 8^ — 8^ + °(z-5), z ^ ж : 2 1 1 1

Sx(z) - z + — ^ + ^ + °(z-5).

2.3. Задача о вторичных течениях. Обратимся к нахождению членов первого порядка в асимптотическом представлении (6). Нас будут интересовать лишь их стационарные компоненты, описывающие вторичные течения. Такие же компоненты следует удерживать в правой части G уравнения (8).

2.3.1. Вычисление правой части. Выделим стационарные компоненты в О. Используя полученные выше в п. 2.2. выражения для ф(0),ш(0), получим

rG =1 sin 2в

Re (Viei(i+^) Re (Viei(t+v)) - Re (^iei(t+v)) Re (Vi + £ cos в Re (Ф2eit) Re (Wiei(t+v)) - Re Weit) Re (^iei(t+v))

+

Штрихи здесь и далее обозначают дифференцирование по т. Возьмем среднее от произведений в правой части. Обозначая чертой сверху комплексное сопряжение, получим

^Ие (Ф2еЙ) ■Ие (ж^^) ^ = 1 ^ (^'2еи + Ф2е_") + ^ =

= ^'^е^2^ +Ф2^1е-4(24+^) 1 е-11* + Ф2=

= 4 (О + 0 + 1е-1* + = 2ие (Ф2Ж ^^) .

Проводя аналогичные выкладки для остальных произведений, найдем

(G) = Gx sin 2в + £Gg cos в, Gx(r; в) = 4rRe (Vi^i - ^iWi) , Ge(r; в, y) = 1 Re (e-i^ (Wi*2 - )) . Используя полученные в п. 2.2. формулы для ^i, Ф2, Wi, W2 и обозначая

Ф, (r; в) = Kk(,f^, 5 (в ) = Sx + 1 = ,

kK)l> Kk (уДву x^ ЫУФУ

найдем Gx, Gg . Для Gg имеем

Gg(r; в, У) = T^Re (e-inWiV - W!2 (-r + Sr-i - W-i

2r V V гв 2 \ гв

= -^Re ( e-ivWU 2Wi + (r - -

2r у 2 у гв i 1 V r

С учетом

W2 = ¿вФь Wi = -¿вФ^- - 1)

получим

Gg(r; в, y) = Re (гвe-iíp ( 1 |Ф1 |2 (- - l) + ^i (1 - -r-2)^

Аналогичные выкладки для Gx дают

в f — Фп —

Gx(r; в) = в Im + 2—Ф2ФП - -Ф2

2.3.2. Решение задачи. Представляя решение в виде

= %X(r)sin2e + (r)cos в, (9)

получим для нахождения Фх и % следующие задачи

1 d d k2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

¿2% = -13GX, L1% = -£Ge, Lk = --r----2; (10)

r dr dr r2

r = 1: Фх = %X = % = % = 0, r ^^ : % =0 , % = (11)

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

Í. = -£ jojrz) (x2 -1)3 *

~ Br4 {

Фв = --— Ge(rx)(x4 - 4x2 lnx - 1)dx. 16 /

(12)

Добавляя к Фx, 'в общие решения однородных уравнений, запишем

Bx „л 1А

(13)

'x - + С*[ 1 - Г2

5=Фв -Б+с [г - Г)+»к1пг-2+

Граничные условия при г = 1 удовлетворяются за счет выбора произвольных констант Бхд, Сх,е

в,= 5х,(1), С = -М+М>, Сх = - «да^УИ- (14)

Функция 5х при г ^ ж выходит на константу Сх, что не противоречит граничным условиям на бесконечности. Функция 5е имеет классическую логарифмическую особенность на бесконечности. Определить константу Б и удовлетворить условию для функции тока при г ^ ж для произвольной скорости стационарного потока без специального учета конвективных членов невозможно (парадокс Стокса). Лишь при специальном выборе = Се (в) и а = 0 условие на бесконечности выполняется, а константа Б обращается в нуль. Оказывается, что указанный выбор в точности соответствует принятому нами требованию равенства нулю средних за период сил, действующих на тело со стороны жидкости (см. [15]).

Отсюда следует, что крейсерская скорость, с которой будет перемещаться в неподвижной жидкости свободный цилиндр, направлена вдоль оси Оу и равна "то = Се.

3. Анализ результатов

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

uто = Ux,c(P) cos ф - Ux,s(P) sin ф.

1

0

-1

100 102 104 Р

Рис. 2. Зависимость крейсерской скорости от безразмерной частоты в при ф = 0 (их = = их,с) и ф = —п/2 (иж = их,а). Штриховыми линиями показаны асимптотики различных порядков

Величины иж,с и подсчитываются по формулам (14) с заменой в них Фд на

Фдс и Фд}В соответственно. Экспоненциально убывающие функции Фдс и Фд}В задаются интегралом (12) с заменой в нем ядра Од на вещественнозначные функции

Од,с и Од,в

Од,в + гОдс = в |Ф1|2 (5 - 1) + 2 Ф1 (1 - Бт-2) . т 2

Зависимости ижс(в), их,е(в) показаны сплошными линиями на рис. 2. Первая из них отвечает случаю р = 0 синфазных колебаний. Она была получена ранее в [15]. Вторая - случаю р = —п/2.

Как было установлено в [15], для синфазных колебаний (когда иж = ижс) имеются два вида среднего стационарного движения цилиндра: прямое (ито > 0) и обратное (ито < 0). Первый из них реализуется при в > вд = 13.5355 и характеризуется монотонным ростом скорости ито с ростом в. Максимальная размерная

33

скорость такого движения (см рис. 2), равная 4 ^Си = 4 ©и, вполне соизмерима с максимальной скоростью и колебаний цилиндра. Обратное движение происходит при в < вд . Его безразмерная скорость достигает наименьшего значения иж = = —0.0534 при в = 2.18. Эта скорость более чем на порядок меньше наибольшей безразмерной скорости иж = 0.75 прямого движения.

Для случая иж = (когда фазы вращательных и поступательных колеба-

ний смещены на р = —п/2) полученные результаты показывают наличие лишь одного вида среднего стационарного движения цилиндра - обратного (иж < 0). Безразмерная скорость монотонно убывает с ростом в и достигает наименьшего 55

значения = — 4 ^Си = — 4 ©и. С изменением фазы колебаний на противоположную (р = п/2) крейсерская скорость, очевидно, меняет свой знак.

Асимптотическое поведение крейсерской скорости в практически важном случае в ^ ж можно найти с использованием известного асимптотического разложения модифицированных функций Бесселя при больших значениях аргумента

и леммы Ватсона:

Сд,с = 4 — 7 /2 в-1/2 + ^ в-1 — 8§ в-3/2 + О (в-2),

2 (15)

Сд а = — 5 + / в-1/2 + 187 в-1 — ^ в-3/2 + О (в-2) . д'3 4 4^ 8 ^ 32 л/2 ^ ^

Для сравнения они, наряду с точным результатом, показаны штриховыми линиями на рис. 2.

На практике для вычисления иж,с(в), иж,3(в) удобно пользоваться аппроксимациями Паде. Для иж,с имеем

(д) в — 3.7 /в

иж, с(в) =

1.33в + 12.7 л/в + 41.7'

Представленная аппроксимация точна в главных двух членах асимптотики в ^ ж и верно определяет нуль их,с(в). При в > 0.5 ее относительная погрешность не превышает 5%. Аналогичная аппроксимация

и ^ = - в +1.4 /в и™'3(Р) 0.8в + 3.5 /в + 3.6

также точна в главных двух членах асимптотики в ^ ж>. При в > 10-3 ее относительная погрешность не превышает 3%.

Заключение

В работе проведено исследование задачи о направленном движении цилиндрического крыла круглого сечения, совершающего поступательно-вращательные колебания в вязкой несжимаемой жидкости. Аналитически найдены первые два члена асимптотического разложения решения задачи по малым амплитудам колебаний (7 ^ 1): первый член описывает первичное нестационарное стоксово течение, второй - установившееся течение, формируемое в результате нелинейного взаимодействия временных гармоник. Показано, что взаимодействие временных гармоник поступательных и вращательных колебаний приводит к формированию ненулевой средней скорости движения. Полученное решение позволяет описать режим крейсерского движения, когда средняя сила, действующая на тело в период колебаний, равна нулю. Результаты исследования показывают, что максимальная по модулю безразмерная крейсерская скорость составляет 3/4© (три четверти максимального угла поворота) при синфазных поступательно-вращательных колебаниях и 5/4© для случая, когда фазы вращательных и поступательных колебаний смещены на п/2.

Благодарности. Работа выполнена при финансовой поддержке Российского научного фонда в рамках научного проекта № 22-79-10033.

Литература

1. Prandtl L. Über die Entstehung von Wirbeln in der idealen Flüssigkeit, mit Anwendung auf die Tragflügeltheorie und andere Aufgaben // v. Kärmän T., Levi-Civita T. (Hrsg.) Vortrage aus dem Gebiete der Hydro- und Aerodynamik (Innsbruck 1922). - Berlin, Heidelberg: Springer, 1924. - S. 18-33. - doi: 10.1007/978-3-662-00280-3_2.

2. Birnbaum W. Der Schlagflugelpropeller und die kleinen Schwingungen elastisch befestigter Tragflügel // Z. Flugtech. Motorluftschiffahrt. - 1924. - Bd. 15, H. 11-12. - S. 128-134.

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

4. Theodorsen T. General Theory of Aerodynamic Instability and the Mechanism of Flutter: NACA Rep. No 496. - Washington: Supt. Doc., 1935. - 23 p.

5. Некрасов А.И. Теория крыла в нестационарном потоке. - М.: Изд-во АН СССР, 1947. - 258 с.

6. Седов Л.И. Плоские задачи гидродинамики и аэродинамики. - М.: Наука, 1966. -448 с.

7. Голубев В.В. Лекции по теории крыла. - М.-Л.: Гос. изд-во техн.-теорет. лит., 1949. -480 с.

8. Голубев В.В. О некоторых вопросах теории машущего крыла // Учен. зап. Моск. гос. ун-та. - 1951. - № 152. - С. 3-12.

9. Schlichting H. Berechnung ebener periodischer Grenzschichtstromungen // Phys. Z. -1932. - Bd. 33. - S. 327-335.

10. Holtsmark J., Johnsen I., Sikkeland T., Skavlem S. Boundary layer flow near a cylindrical obstacle in an oscillating, incompressible fluid // J. Acoust. Soc. Am. - 1954. - V. 26, No 1. - P. 26-39. - doi: 10.1121/1.1907285.

11. Riley N. The steady streaming induced by a vibrating cylinder // J. Fluid Mech. - 1975. -V. 68, No 4. - P. 801-812. - doi: 10.1017/S0022112075001243.

12. Wang Ch.-Y. On high-frequency oscillatory viscous flows //J. Fluid Mech. - 1968. -V. 32, No 1. - P. 55-68. - doi: 10.1017/S0022112068000583.

13. Nuriev A.N., Egorov A.G. Asymptotic investigation of hydrodynamic forces acting on an oscillating cylinder at finite streaming Reynolds numbers // Lobachevskii J. Math. -2019. - V. 40, No 6. - P. 794-801. - doi: 10.1134/S1995080219060180.

14. Riley N., Watson E.J. Eccentric oscillations of a circular cylinder in a viscous fluid // Mathematika. - 1993. - V. 40, No 2. - P. 187-202. - doi: 10.1112/S0025579300006975.

15. Nuriev A.N., Egorov A.G. Asymptotic theory of a flapping wing of a circular cross-section // J. Fluid Mech. - 2022. - V. 941. - Art. A23, P. 1-18. - doi: 10.1017/jfm.2022.287.

Поступила в редакцию 29.07.2022

Егоров Андрей Геннадьевич, доктор физико-математических наук, ведущий научный сотрудник НИЛ «Интеллектуальные биомиметические и природосообразные системы» Института математики и механики им. Н.И. Лобачевского Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: aegorov0@gmail.com

Нуриев Артем Наилевич, кандидат физико-математических наук, доцент кафедры аэрогидромеханики

Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: nuriev an@mail.ru

ISSN 2541-7746 (Print)

ISSN 2500-2198 (Online)

UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI

(Proceedings of Kazan University. Physics and Mathematics Series)

2022, vol. 164, no. 2-3, pp. 170-180

ORIGINAL ARTICLE

doi: 10.26907/2541-7746.2022.2-3.170-180

Cruising Speed of a Cylindrical Wing Performing Small Translational-Rotational Oscillations

A.G. Egorov* , A.N. Nuriev**

Kazan Federal University, Kazan, 420008 Russia E-mail: * aegorov0@gmail.com, ** nuriev_an@mail.ru

Received July 29, 2022 Abstract

This work considers the propulsive motion of a flapping wing of a circular cross section. The problem of harmonic translational-rotational oscillations of the wing with an arbitrary phase shift in a viscous incompressible fluid, the motion of which is described by the non-stationary Navier-Stokes equation, is handled. An analytical solution of the problem is obtained in the first two terms by using the method of successive asymptotic expansions for the case of small oscillation amplitudes. It is shown that the nonlinear interaction of time harmonics of translational and rotational oscillations creates secondary flows that make the wing to move in the direction perpendicular to the axis of translational oscillations. For the cruising motion regime, when the average hydrodynamic force acting on the wing is equal to zero, the dependences of the average speed on the parameters of dimensionless oscillation are found.

Keywords: flapping wing, harmonic oscillations, cruising speed, Navier-Stokes equation, asymptotic analysis

Acknowledgments. This study was supported by the Russian Science Foundation (project no. 22-79-10033).

Figure Captions

Fig. 1. Scheme of in-phase lateral-torsional oscillations of the cylinder. The extreme position of the cylinder is shown by the dotted line. A, © are amplitudes of the lateral and torsional oscillations, H = A/© is the distance between the center of the cylinder O and the center of speeds Oi in the given position.

Fig. 2. Dependence of the cruising speed on the dimensionless frequency ß at p = 0 (= ) and p = —n/2 (= ). Asymptotics of different orders are indicated by the dashed lines.

References

1. Prandtl L. Uber die Entstehung von Wirbeln in der idealen Flussigkeit, mit Anwendung auf die Tragflugeltheorie und andere Aufgaben. In: v. Karman T., Levi-Civita T. (Hrsg.) Vorträge aus dem Gebiete der Hydro- und Aerodynamik (Innsbruck 1922). Berlin, Heidelberg, Springer, 1924, S. 18-33. doi: 10.1007/978-3-662-00280-3_2. (In German)

2. Birnbaum W. Der Schlagflügelpropeller und die kleinen Schwingungen elastisch befestigter Tragflügel. Z. Flugtech. Motorluftschiffahrt, 1924, Bd. 15, H. 11-12, S. 128-134. (In German)

3. Keldysh M.V., Lavrent'ev M.A. On the theory of the oscillating wing. Tekh,. Zametki, 1935, no. 45: Collection of articles by the General Theoretical Group at Central Aerohy-drodynamic Institute, pp. 48-52. (In Russian)

4. Theodorsen T. General Theory of Aerodynamic Instability and the Mechanism of Flutter. NACA Rep. no. 496. Washington, Supt. Doc., 1935. 23 p.

5. Nekrasov A.I. Teoriya kryla v nestatsionarnom potoke [Wing Theory in Unsteady Flow]. Moscow, Izd. Akad. Nauk SSSR, 1947. 258 p. (In Russian)

6. Sedov L.I. Ploskie zadachi gidrodinamiki i aerodinamiki [Flat Problems of Hydrodynamics and Aerodynamics]. Moscow, Nauka, 1966. 448 p. (In Russian)

7. Golubev V.V. Lektsii po teorii kryla [Lectures on the Theory of Wing]. Moscow, Leningrad, Gos. Izd. Tekh.-Teor. Lit., 1949. 480 p. (In Russian)

8. Golubev V.V. On some problems of the theory of a flapping wing. Uch. Zap. Mosk. Gos. Univ., 1951, no. 152, pp. 3-12. (In Russian)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

9. Schlichting H. Berechnung ebener periodischer Grenzschichtstromungen. Phys. Z., 1932, Bd. 33, S. 327-335. (In German)

10. Holtsmark J., Johnsen I., Sikkeland T., Skavlem S. Boundary layer flow near a cylindrical obstacle in an oscillating, incompressible fluid. J. Acoust. Soc. Am., 1954, vol. 26, no. 1, pp. 26-39. doi: 10.1121/1.1907285.

11. Riley N. The steady streaming induced by a vibrating cylinder. J. Fluid Mech., 1975, vol. 68, no. 4, pp. 801-812. doi: 10.1017/S0022112075001243.

12. Wang Ch.-Y. On high-frequency oscillatory viscous flows. J. Fluid Mech., 1968, vol. 32, no. 1, pp. 55-68. doi: 10.1017/S0022112068000583.

13. Nuriev A.N., Egorov A.G. Asymptotic investigation of hydrodynamic forces acting on an oscillating cylinder at finite streaming Reynolds numbers. Lobachevskii J. Math., 2019, vol. 40, no. 6, pp. 794-801. doi: 10.1134/S1995080219060180.

14. Riley N., Watson E.J. Eccentric oscillations of a circular cylinder in a viscous fluid. Mathematika, 1993, vol. 40, no. 2, pp. 187-202. doi: 10.1112/S0025579300006975.

15. Nuriev A.N., Egorov A.G. Asymptotic theory of a flapping wing of a circular cross-section. J. Fluid Mech., 2022, vol. 941, art. A23, pp. 1-28. doi: 10.1017/jfm.2022.287.

Для цитирования: Егоров А.Г., Нуриев А.Н. Крейсерская скорость цилиндри-/ ческого крыла при малых поступательно-вращательных колебаниях // Учен. зап. \ Казан. ун-та. Сер. Физ.-матем. науки. - 2022. - Т. 164, кн. 2-3. - С. 170-180. - doi: 10.26907/2541-7746.2022.2-3.170-180.

For citation: Egorov A.G., Nuriev A.N. Cruising speed of a cylindrical wing perfor-/ ming small translational-rotational oscillations. Uchenye Zapiski Kazanskogo Universi-\ teta. Seriya Fiziko-Matematicheskie Nauki, 2022, vol. 164, no. 2-3, pp. 170-180. doi: 10.26907/2541-7746.2022.2-3.170-180. (In Russian)

i Надоели баннеры? Вы всегда можете отключить рекламу.