Колебания тел с острыми кромками в несжимаемой маловязкой жидкости и некоторые задачи гидродинамики космических аппаратов тема диссертации и автореферата по ВАК РФ 01.02.05, доктор физико-математических наук Бужинский, Валерий Алексеевич

Оглавление диссертации доктор физико-математических наук Бужинский, Валерий Алексеевич

1 КОЛЕБАНИЯ ТЕЛА С ОСТРЫМИ КРОМКАМИ В

НЕСЖИМАЕМОЙ МАЛОВЯЗКОЙ ЖИДКОСТИ.

1.1 Сингулярные решения уравнения Лапласа.

1.2 Энергия вихреобразования за период колебаний пластинки в несжимаемой маловязкой жидкости.

1.3 Гидродинамический аналог формулы Ирвина.

1.4 Определение коэффициентов сопротивления и декрементов колебаний.

1.5 Экспериментальное определение универсальной постоянной.

1.6 Влияние вихревого течения в окрестности острых кромок на присоединенные к пластинкам массы жидкости.

1.7 Применение гидродинамического аналога формулы

Ирвина для получения аналитических решений.

1.8 Аналитические решения методом конформных отображений.

1.9 Аналитические решения методом интеграла типа Коши.

1.10 Некоторые аналитические решения пространственной задачи.

1.11 Общий случай колебаний в жидкости тела с двугранными кромками. щ 2 МЕТОД ГРАНИЧНЫХ ЭЛЕМЕНТОВ ДЛЯ ПЛОСКИХ

ЗАДАЧ О ПОТЕНЦИАЛЕ С НЕЗАМКНУТЫМИ

2.1 Сингулярное граничное интегральное уравнение.

2.2 Постоянные граничные элементы.

2.3 Квадратичные граничные элементы.

2.4 Вычисление интегралов.

2.5 Численные результаты и анализ сходимости.

2.6 Влияние границ на коэффициенты сопротивления и декременты колебаний.

2.6.1 Влияние стенки. Эффект зазора.

2.6.2 Влияние другой пластинки. Эффект затенения. 2.6.3 Влияние свободной поверхности жидкости.

2.6.4 Колебания заполненной жидкостью цилиндрической

• емкости с радиальными ребрами вокруг продольной оси.

2.7 Плоская стационарная задача гидродинамики.

2.7.1 Сингулярное граничное интегральное уравнение.

2.7.2 Метод дискретных вихрей в стационарных задачах гидродинамики.

2.7.3 Некоторые аналитические и численные решения.

3 МЕТОД ГРАНИЧНЫХ ЭЛЕМЕНТОВ ДЛЯ

ПРОСТРАНСТВЕННЫХ ЗАДАЧ О ПОТЕНЦИАЛЕ С

НЕЗАМКНУТЫМИ ГРАНИЧНЫМИ 133 ПОВЕРХНОСТЯМИ

3.1 Интегральные представления гармонической функции и щ производных гармонической функции.

3.2 Сингулярное граничное интегральное уравнение.

3.3 Определение предельных значений интегралов.

3.4 Аппроксимация граничных поверхностей.

3.4.1 Основные зависимости.

3.4.2 Плоский четырехугольный элемент поверхности.

3.4.3 Плоский треугольный элемент поверхности.

3.4.4 Квадратичный четырехугольный элемент поверхности.

3.4.5 Квадратичный треугольный элемент поверхности.

3.5 Приведение СГИУ к системе линейных алгебраических уравнений.

3.5.1 Постоянные граничные элементы.

3.5.2 Линейные граничные элементы.

3.6 Результаты расчетов для колебаний пластин прямоугольной формы в маловязкой жидкости. Сравнение с экспериментальными данными.

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

3.6.2 Некоторые результаты расчетов МГЭ.

3.6.3 Экспериментальное исследование сопротивления прямоугольных пластин. . .«•

3.7 Влияние воздушной среды при частотных испытаниях

3.8 О применении метода граничных элементов к задачам обтекания тел потоком несжимаемой жидкости.

3.8.1 Синтез метода граничных элементов с методом дискретных вихрей.

3.8.2 Сила давления жидкости на тело в потоке с поверхностями тангенциального разрыва.

4 ДВИЖЕНИЕ ТЕЛА С ПОЛОСТЬЮ, ЧАСТИЧНО ЗАПОЛНЕННОЙ ЖИДКОСТЬЮ, С

4.1 Уравнения возмущенного движения.

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

4.3 Интегральная формула для коэффициентов интенсивности скоростей.

4.4 Применение метода конечных элементов к решению краевых задач для резервуаров вращения с жидкостью.

4.4.1 Вариационная формулировка задач.

4.4.2 Линейные конечные элементы.

4.4.3 Квадратичные конечные элементы.

4.4.4 Задача на собственные значения для отсеков резервуаров вращения.

4.5 Определение гидродинамических параметров.

4.6 Характеристика программного обеспечения.

4.7 Анализ эффективности применения МКЭ для определения гидродинамических параметров.

4.7.1 Исследование влияния поперечных перегородок на инерционные гидродинамические параметры.

4.7.2 Исследование влияния поперечных перегородок на демпфирование колебаний жидкости.

4.8 Расчет гидродинамических параметров для топливных баков разгонного блока.

4.9 Определение требований к демпфированию колебаний жидкого топлива путем анализа устойчивости движения частотным методом.

4.9.1 Основные уравнения.

4.9.2 Частотный метод Найквиста.

4.9.3 Требования к демпфированию колебаний жидкости в баках РБ "Бриз-КМ" при полете в составе II ступени РН

4.9.4 Анализ структурных свойств III ступени

РН "Протон" с РБ ДМ.

Рекомендованный список диссертаций по специальности «Механика жидкости, газа и плазмы», 01.02.05 шифр ВАК

Введение диссертации (часть автореферата) на тему «Колебания тел с острыми кромками в несжимаемой маловязкой жидкости и некоторые задачи гидродинамики космических аппаратов»

Актуальность темы. Интерес к задаче о колебаниях в жидкости тела с острыми кромками резко возрос в начале второй половины XX века в связи с развитием ракетно-космической техники (РКТ). Жидкое топливо в баках ракет-носителей (РН) может составлять до 90% от их общей массы. Подвижность жидкости существенно влияет на движение РН, разгонных блоков (РБ) и космических аппаратов (КА) и может приводить к потере устойчивости их движения. Практически единственный способ подавления такой динамической неустойчивости состоит в резком увеличении демпфирования колебаний жидкого топлива. С этой целью в топливные баки устанавливают конструктивные элементы типа продольных и поперечных перегородок (пластин), частично перекрывающих продольное или поперечное сечение полости. При колебаниях жидкости на острых кромках этих конструктивных элементов периодически происходит образование и срыв вихрей, за счет чего и обеспечивается высокое демпфирование. Задачи о колебаниях жидкости в баках с демпфирующими перегородками и колебаний пластин, а в общем случае тел с острыми кромками, в несжимаемой маловязкой жидкости тесно связаны.

Выбрать демпфирующие перегородки с минимальным весом, которые обеспечивают устойчивость движения объектов РКТ, не имея общей теории, достаточно сложно: объем материальных затрат и времени на проведение испытаний чрезмерно велик.

Несколько позже возникла задача об учете влияния воздушной среды при наземных частотных испытаниях панелей солнечных батарей (ПСБ) КА. ПСБ КА представляют собой относительно легкие, но крупногабаритные конструкции. Влияние воздушной среды проявляется в появлении дополнительных инерционных и диссипативных сил, которые искажают определение динамических характеристик ПСБ, необходимых для настроек системы управления КА. Возникающие диссипативные силы обусловлены, главным образом, сложным вихревым движением воздуха около кромок ПСБ.

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

Определение сил, действующих на движущиеся в жидкости тела, является одной из центральных проблем гидромеханики. В общем случае решение уравнений Навье-Стокса dv 1 (vV)v = — VP + vAv dt p для определения этих сил при больших числах Рейнольдса наталкивается на трудности принципиального характера [27,51]. Поэтому развитие теоретической гидромеханики определялось, прежде всего, созданием приближенных моделей движения жидкостей, отражающих наиболее существенные черты рассматриваемых в них явлений, а также разработкой численных методов их исследования. Настоящая работа посвящена разработке полуэмпирической асимптотической теории вихревого сопротивления тел с острыми кромками при колебаниях в несжимаемой маловязкой жидкости и ее приложению к решению сформулированных выше задач гидродинамики космических аппаратов.

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

Одна из наиболее близких к теме исследования задач - это задача о колебаниях в несжимаемой жидкости твердого тела с гладкой граничной поверхностью. Она имеет общее асимптотическое решение [51], которое справедливо, если мала амплитуда колебаний тела и толщина пограничного слоя много меньше характерного размера тела (точнее, радиусов кривизны его поверхности). Первое из допущений необходимо, чтобы пренебречь конвективным ускорением и считать течение жидкости безотрывным, а второе - для пренебрежения кривизной течения жидкости в пограничном слое. Поэтому в первом приближении для любого малого элемента поверхности твердого тела можно использовать точное решение уравнений Навье-Стокса для колебаний плоскости под безграничным слоем жидкости, полученное Стоксом в 1851 году.

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

Уравнения движения тела с полостью, частично заполненной идеальной несжимаемой жидкостью, в начале второй половины XX века получили Моисеев [57], Нариманов [61], Охоцимский [62], а также зарубежные ученые [88]. Эквивалентные уравнения в несколько ином виде предложил Рабинович [65]. Определение коэффициентов полученных обыкновенных дифференциальных уравнений было сведено к решению двух типов краевых задач. Первая из них такая же, как для тела с полостью, полностью заполненной жидкостью. Фундаментальное исследование о движении тела с полостью, содержащей жидкость, было выполнено Жуковским

39] без связи с какими-либо практическими запросами в то время и опубликовано в 1885 году. Вторая краевая задача описывает колебания жидкости в полости неподвижного тела. Некоторые решения этой задачи были известны также еще в XIX веке [50].

Строгое решение задачи о движении тела с полостью, заполненной маловязкой жидкостью, получил Черноусько [86]. В частности, он разработал теорию малых колебаний жидкости, имеющей свободную поверхность, в резервуарах с гладкими стенками путем решения уравнений На-вье-Стокса в виде асимптотических разложений в ряд по степеням вязко

1/9 сти v . Однако демпфирование колебаний жидкости, обусловленное диссипацией энергии в пограничном слое на стенках баков объектов РКТ, оказывается слишком малым. Во многих случаях при таком демпфировании устойчивость движения РН, РБ и КА не обеспечивается. В окрестности острых кромок конструктивных элементов, устанавливаемых для резкого повышения демпфирования, допущения теории пограничного слоя не выполняются.

Течения жидкости около тел с острыми кромками изучаются в теории крыла самолета. Наличие острой кромки - принципиальная особенность крыла. В 1906 году Жуковский [40] установил вихревую природу подъемной силы, а в 1909 году, когда был сформулирован постулат Чаплыгина-Жуковского, стала ясной роль острой кромки. Книга Белоцерков-ского [6] положила начало интенсивному развитию вихревых моделей течений идеальной жидкости. Численные решения задач обтекания тел потоком идеальной жидкости получают методом дискретных вихрей (МДВ). Обоснование МДВ дано Лифановым [7,52]. Применение МДВ для решения задач о колебаниях тела с острыми кромками в жидкости и, в особенности, для решения задач о колебаниях жидкости в баках с демпфирующими перегородками связано с большими трудностями. О полученных этим методом некоторых частных решениях задач о колебаниях пластинок в жидкости будет сказано ниже.

Седов [74] получил асимптотическое решение плоской задачи о колебаниях пластинки, расположенной под малым углом атаки в стационарном на бесконечности потоке жидкости с большой скоростью. Это решение имеет некоторое практическое приложение для анализа процессов при вибрациях крыла самолета. Задача о колебаниях пластинки в покоящейся на бесконечности жидкости такого решения не имеет.

Одно из первых систематических экспериментальных исследований по малым колебаниям различных тел в жидкости, в частности пластин разной формы, выполнили Риман и Крепе [69]. Однако в их работе определялись только коэффициенты присоединенных масс жидкости.

Впервые глубокие экспериментальные исследования по сопротивлению длинной прямоугольной пластинки, расположенной перпендикулярно периодическому потоку жидкости, были проведены Келеганом и Карпен-тером [98] и опубликованы в 1958 году. Сопротивление тонкой пластинки при колебаниях в безграничной жидкости может зависеть только от двух критериев подобия: числа Струхаля или относительной амплитуды колебаний и числа Рейнольдса. Испытания показали, что коэффициент сопротивления практически не зависит от числа Рейнольдса при его изменении в очень широком диапазоне больших значений. Это указывает на то, что значение имеет энергия образующихся за период колебаний вихрей, а не эффекты, связанные с их диффузией и диссипацией в объеме жидкости.

Келеган и Карпентер представили результаты своих исследований в виде зависимости коэффициента сопротивления от безразмерного параметра vT/R , где v - амплитуда скорости, R - ширина пластинки, Т - период колебаний. Иногда в литературе этот параметр, обратный числу Струхаля, называется числом Келегана-Карпентера и обозначается Кс. В проведенных Келеганом и Карпентером испытаниях это число изменялось в диапазоне 2<v77i?<50 . Даже для нижней границы указанного диапазона это соответствует довольно большой амплитуде колебаний, составляющей около 1/3 ширины пластинки.

В том же 1958 году Майлс [100], основываясь на этих экспериментальных данных, аппроксимировал зависимость коэффициента сопротивления пластинки от числа Кс простейшей степенной функцией с показателем -1/2 и применил ее для оценки демпфирования колебаний жидкости в цилиндрическом баке с кольцевой перегородкой. В результате демпфирование оказалось пропорционально амплитуде колебаний жидкости в степени 1/2.

В СССР обширные экспериментальные исследования колебаний пластин в жидкости и колебаний жидкости в баках с демпфирующими перегородками были проведены Микишевым [54,55] и Бенедиктовым [29,30]. В этих исследованиях больше внимания уделялось малым амплитудам колебаний, представляющим практический интерес в связи с задачей о демпфировании колебаний жидкого топлива в баках РН и КА. Путем аппроксимации экспериментальных данных были получены различные частные эмпирические зависимости, рекомендованные для применения в конкретных случаях.

Первое основательное теоретическое исследование по сопротивлению пластинки, совершающей колебания в безграничной жидкости, принадлежит Грахаму [92] и опубликовано им в 1980 году в ведущем журнале по механике жидкости "Journal of Fluid Mechanics". Он рассмотрел плоскую задачу о гармонических колебаниях пластинки по нормали к ее плоскости при малых числах Келегана-Карпентера Кс. Из ряда других эту работу выделяет то, что в ней впервые был сделан правильный общий вывод о характере сопротивления. Применяя соображения подобия, Грахам показал, что при малых амплитудах колебаний показатель степени в зависимости сопротивления от числа Келегана-Карпентера Кс равен -1/3. Позже

Бернардинес, Грахам и Паркер [91] аналогично исследовали осесиммет-ричные задачи о колебаниях тонкого круглого диска и о периодическом течении жидкости через круглое отверстие.

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

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

Основная идея работы. Общее теоретическое решение проблемы совершенно другим путем было получено автором в 1987 году после знакомства с небольшой книгой Работнова [67] по линейной механике разрушения. В этой книге привлекло внимание то, что, основываясь на соображениях подобия, по значению коэффициента интенсивности напряжений (КИН) оценивался размер зоны пластичности в окрестности вершины трещины, т.е. методы подобия и размерностей соединялись с методом сингулярных решений. Трещина в упругом теле и пластина в идеальной несжимаемой жидкости имеют не только очевидное внешнее сходство. Задачи теории упругости, как и задачи о потенциальном движении жидкости, относятся к эллиптическим краевым задачам. С математической точки зрения решения этих задач имеют одинаковые особенности в окрестности острых кромок областей. Возникла идея перенесения методов теории трещин механики деформируемого твердого тела в область гидродинамики при рассмотрении колебаний пластин в несжимаемой маловязкой жидкости. Эта идея оказалась плодотворной, поэтому кратко остановимся на основных моментах развития и принципах линейной механики разрушения.

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

В 1957 году Ирвин [96] ввел в качестве меры сингулярности напряжений на острых кромках трещин коэффициенты интенсивности напряжений (КИН) и связал с ними силовой критерий прочности. Он показал, что в малой окрестности острой кромки любой трещины напряжения выражаются зависимостями вида a = K/yf2nr /(9), где / - некоторые известные функции полярного угла 9 , г - расстояние до острой кромки, К - коэффициент интенсивности напряжений (трем разным типам деформаций соответствуют коэффициенты Kj,Kjj,Kjjj). Ирвин получил также важные формулы, носящие его имя, которые устанавливают связь КИН с изменением энергии при росте трещин. С этого момента началось интенсивное развитие линейной механики трещин, которая к настоящему времени сформировалась в обширный самостоятельный раздел науки [68,85], имеющий важнейшие практические приложения.

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

Заимствуя идеи линейной механики разрушения, автор [12,13] предложил новую приближенную модель движения несжимаемой маловязкой жидкости, предназначенную для определения сил сопротивления, действующих на тела с острыми кромками при их чисто колебательном движении в жидкости либо при колебательном движении жидкости относительно тел с острыми кромками. В этой модели область течения жидкости разбивается на следующие подобласти:

• область основного потенциального движения жидкости;

• область малой "дальней" окрестности острой кромки, где течение жидкости описывается главными сингулярными членами решения задачи о потенциальном движении, в котором скорости жидкости на острой кромке обращаются в бесконечность;

• область ближней окрестности острой кромки и существенного вихревого движения жидкости, где происходит образование и периодическое изменение вихревых структур;

• область пограничного слоя.

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

Разработанная автором [12,13,17] асимптотическая теория вихревого сопротивления применима при больших числах Рейнольдса и малых числах Келегана-Карпентера. Она основана на введении количественной меры сингулярности скорости жидкости на острой кромке согласно зависимости v = К у / v 2лг / (0) , где Kv - коэффициент интенсивности скоростей (КИС), /(0) - некоторая известная функция полярного угла 0, отсчитываемого от одной стороны кромки к другой по дуге в жидкости, г - расстояние до острой кромки. С привлечением соображений подобия и размерностей показано, что с точностью до универсального постоянного множителя энергия вихреобразования за период колебаний зависит вполне определенным образом только от плотности жидкости, частоты колебаний и значений КИС на контуре острой кромки.

Основной результат теории может быть выражен зависимостью dE/de = Bp<a-2/2Kl'3(£) (0.1) где dEl di - производная энергии вихреобразования за период колебаний по длине острой кромки, т.е. энергия вихреобразования на единицу длины кромки, р - плотность жидкости, ю - частота колебаний, Kv (£) - значение

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

Для определения сил сопротивления, которые действует на пластину или панель, совершающую колебания в жидкости или воздушной среде, выполняются следующие операции [12,17]:

• решается задача о потенциальном движении жидкости;

• определяются коэффициенты при сингулярных членах потенциального решения, названные коэффициентами интенсивности скоростей (КИС), которые характеризуют течение жидкости в малых "дальних" окрестностях острых кромок;

• вычисляются интегралы по контурам острых кромок, в которых подынтегральные функции зависят от КИС.

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

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

Твердое тело может иметь не только острые кромки, но и кромки с отличным от нуля двугранным углом. Примером такого тела может быть, например, куб, двугранные углы всех кромок которого равны к/2. Если тело имеет кромку, двугранный угол которой меньше я, то в потенциальном движении скорость жидкости в вершине кромки обращается в бесконечность, а течение в некоторой ее окрестности при малых амплитудах колебаний описывается сингулярной функцией. Обобщение теории дано и на этот общий случай формы поверхности тела [13]. Однако в полученное конечное решение входит неизвестная функция В(а) , где а - величина двугранного угла. При а = 0, что соответствует острой кромке пластинки, эта функция принимает значение указанной выше универсальной постоянной. Функция В(а) должна быть определена в диапазоне 0 < а < п путем специально поставленных экспериментов. Дальнейшее развитие теории возможно только после накопления новых экспериментальных фактов.

Собственно течение жидкости в малых ближних окрестностях острых кромок, имеющее весьма сложную вихревую структуру, в описанной модели остается неизвестным. Впрочем, такое положение, когда определяются только интегральные характеристики явления, довольно типично даже для фундаментальных научных теорий. Так как используемые эмпирические соображения имеют не частный, а достаточно общий характер, и многочисленные экспериментальные данные подтверждают теорию, то она не хуже и не лучше некоторых других моделей гидромеханики, ценность и полезность которых не вызывает сомнений. Так в теории несущей поверхности [6,103] существование вполне определенной вихревой структуры в потоке жидкости принимается из физических соображений, основанных на наблюдениях, а не получено из уравнений Навье-Стокса. Определенной платой за такой подход является большая размытость границ области применимости.

Коэффициенты интенсивности скоростей (КИС), через которые определяется энергия вихреобразования за период колебаний, являются количественной мерой сингулярности скорости жидкости на острых кромках. Требуется не только решить краевые задачи для уравнения Лапласа, имеющие сингулярные особенности на некоторых линиях, но и достаточно точно вычислить КИС, характеризующие количественно эти особенности. Обычные методы для этого не подходят, поэтому значительная часть работы посвящена методам вычисления КИС. Используются как известные методы, заимствованные из линейной механики разрушения, так и развиваются новые, основанные на полученных автором [16,17] новых сингулярных граничных интегральных уравнениях. Применяются различные аналитические методы и общие численные методы конечных элементов и граничных элементов.

Краткое содержание и предварительная сводка основных результатов. Диссертация состоит из четырех глав (разделов). В первой главе излагается новая асимптотическая теория вихревого сопротивления тел с острыми кромками при малых колебаниях в несжимаемой маловязкой жидкости [12,13]. Рассматриваются различные аналитические методы определения вихревого сопротивления на основе этой теории [17]. Приводятся многочисленные примеры решенных задач, некоторые из которых имеют важное практическое значение. Исследуется также влияние вихревого движения жидкости в окрестности острых кромок на присоединенные массы жидкости [14].