авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

Смешанные задачи удара твердых тел, плавающих на поверхности несжимаемой жидкости

На правах рукописи

НОРКИН Михаил Викторович

СМЕШАННЫЕ ЗАДАЧИ УДАРА ТВЕРДЫХ ТЕЛ,

ПЛАВАЮЩИХ НА ПОВЕРХНОСТИ

НЕСЖИМАЕМОЙ ЖИДКОСТИ

01.02.05 механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

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

Казань

2010

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

Научный консультант: доктор физико-математических наук, профессор Жуков Михаил Юрьевич

Официальные оппоненты: доктор физико-математических наук, про фессор Маклаков Дмитрий Владимирович доктор физико-математических наук, профессор Житников Владимир Павлович доктор физико-математических наук, профессор Наседкин Андрей Викторович

Ведущая организация: Институт механики МГУ

Защита состоится 27 мая 2010 г. в 14 часов 30 минут в аудитории мех. 2 на заседании диссертационного совета Д 212.081.11 при Казанском государственном университете по адресу: 420008, г. Казань, ул. Кремлев ская, 18.

C диссертацией можно ознакомиться в научной библиотеке им. Н. И. Ло бачевского Казанского государственного университета.

Автореферат разослан апреля 2010 г.

Ученый секретарь диссертационного совета кандидат физ.-мат. наук, доцент А. А. Саченков

Общая характеристика работы

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

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

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

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

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

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

Методы исследования. Главное внимание при исследовании линей ных задачах гидродинамического удара в ограниченных областях уделя ется развитию эффективного асимптотического метода, основанного на предположении о том, что стенки бассейна удалены от плавающего тела на большие расстояния. Такой подход позволяет свести решение исходной задачи для области сложной геометрической конфигурации к последова тельному решению задач в областях, имеющих более простые формы гра ниц. Это дает хорошую возможность использовать для решения поставлен ных задач методы разделения переменных в специальных ортогональных криволинейных координатах (тороидальных, биполярных, вырожденных биполярных и др.), а также технику парных интегральных уравнений, свя занных со специальными функциями (Бесселя, Лежандра). Для решения смешанных задач гидродинамического удара широко применяются методы граничных интегральных уравнений. В некоторых случаях, на их основе, разработаны специальные асимптотические методы, позволяющие пред ставить основные характеристики удара в простой аналитической форме, удобной для проведения инженерных расчетов.

Большое внимание в диссертации уделено развитию аналитических и численных методов решения задач об ударе с отрывом. Одним из наи более эффективных и, в тоже время, перспективных методов их реше ния является метод нелинейных граничных интегральных уравнений типа Гаммерштейна. Такой подход позволяет одновременно определить неиз вестную заранее зону отрыва частиц жидкости и течение жидкости после удара. Также применяются аналитические методы, основанные на исполь зовании техники конформных отображений в сочетании с математическим аппаратом парных интегральных уравнений. Для определения неизвест ной заранее области отрыва частиц жидкости, наряду с традиционными подходами, применяется вариационный принцип Огазо. Кроме этого ис пользуются специальные асимптотические методы.

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

В диссертационной работе получены следующие новые результаты:

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

2. Решен ряд конкретных смешанных задач в областях сложной формы.

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

4. Разработаны эффективные аналитические методы решения нелиней ных задач гидродинамического удара в ограниченных областях (задач об ударе с отрывом).

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

6. Дано обобщение задачи об отрывном ударе эллиптического цилиндра на случай неоднородной несжимаемой жидкости.

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

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

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

Апробация работы. Основные результаты диссертации докладыва лись на международных конференциях Современные проблемы механики сплошной среды (Ростов-на-Дону, 1995, 1997–2001, 2006–2009);

Нелиней ные задачи теории гидродинамической устойчивости и турбулентность (Москва, 2008);

на семинаре Неустойчивость и турбулентность (Москва, Институт механики МГУ, 2009);

на совместном заседании отделения меха ники НИИММ им. Н. Г. Чеботарева КГУ и кафедры аэрогидромеханики КГУ (Казань, 2009);

на семинаре кафедры математического моделирова ния ЮФУ (Ростов-на-Дону, 2009), неоднократно на семинаре Математи ческие вопросы гидродинамики (кафедра вычислительной математики и математической физики ЮФУ, Ростов-на-Дону, 1993–2009);

на заседании Ростовского математического общества (2008).

Результаты диссертации также докладывались на двух зарубежных на учных семинарах: School of Mathematics, University of East Anglia, Norwich, UK (2009);

Faculty of Mathematics, Physics and Computer Science. Branden burg University of Technology, Cottbus (BTU), Germany (2009).

В диссертацию вошли результаты, полученные автором как руководи телем проекта РФФИ Взаимодействие твердых и упругих тел с жидко стью (01-01-00105, 2001–2003) и как исполнителем проектов: РФФИ (00 15–96188), Президентской программы поддержки ведущих научных школ Российской Федерации (НШ.1768.2003.1, НШ.5747.2006.1), аналитической ведомственной целевой программы Развитие научного потенциала выс шей школы (2.1.1/554;

2.1.1/6095). Исследования по удару с отрывом были также поддержаны грантом The Royal Society grant, UK (Project JP080479) Free-surface separation from a body which starts to move suddenly.

Публикации. По результатам диссертации автором опубликовано работ в ведущих рецензируемых научных журналах и изданиях, входящих в перечень ВАК и одна научная монография, которые приводятся в конце автореферата.

Структура и обьем диссертации. Диссертация состоит из введения, семи глав, приложения и списка литературы. Полный объем диссертации составляет 288 страниц. Список литературы содержит 338 наименований.

СОДЕРЖАНИЕ РАБОТЫ

Введение Во введении обоснована актуальность тематики, сделан обзор литера туры по данной теме, описана структура работы, сформулированы цели диссертации, приведены сведения об апробации работы.

Глава I В этой главе дается постановка задачи о безотрывном, вертикальном ударе плавающего тела и излагается эффективный асимптотический ме тод, направленный на решение пространcтвенных задач гидродинамиче ского удара в областях сложной геометрической конфигурации. В каче стве примера рассматривается задача о вертикальном ударе круглого аб солютно твердого диска, плавающего на поверхности идеальной несжима мой жидкости, ограниченной произвольной поверхностью вращения. Так же приводятся некоторые другие результаты, имеющие общий характер.

В п. 1.1 этой главы, в классической постановке (Л. И. Седов, 1934 г.), рассматривается задача об ударе твердого тела, плавающего на поверхно сти идеальной несжимаемой жидкости, наполняющей ограниченный бас сейн. Удар предполагается таким, что в результате него не происходит отрыва частиц жидкости от смоченной поверхности тела (безотрывный удар). Отметим, что определение для данного тела системы импульсив ных сил и соответствующих им точек приложения, при которых возможен безотрывный удар (определение условий безотрывности удара), является одним из наиболее сложных вопросов теории гидродинамического удара.

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

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

P = 0, r D, (1.1) S = V n, r S1, (1.2) x n = 0, r S2, (1.3) S y = 0, r S3, (1.4) n S = 0, (1.5) z pt = 0, (1.6) Фиг. 1. V n = V0 + R, n = v0 nz + x (ynz zny ) + y (znx xnz ). (1.7) Здесь D область, занятая жидкостью (граница области D предполагает ся кусочно-гладкой);

S1, S2 и S3 соответственно смоченная поверхность твердого тела, свободная поверхность жидкости и неподвижная твердая граница бассейна;

V n нормальная компонента скоростей точек границы тела;

V0 = (0, 0, v0 ) и = (x, y, 0) векторы поступательной и угло вой скорости, приобретенные телом в результате удара;

R = (x, y, z) радиус-вектор точек границы;

n = (nx, ny, nz ) единичный вектор внеш ней нормали к поверхности тела;

pt импульсивное давление;

= const плотность жидкости. Условие (1.5) добавляется только в случае неограни ченной жидкости.

Кроме условий (1.1)–(1.7) нужно еще поставить условие регулярности вблизи линий раздела свободной границы с твердыми поверхностями, где граница области D теряет гладкость. Если особые линии и точки (ребра и конические точки) имеются на границах S1 и S3, то в них также ста вится условие регулярности. Достаточным условием является требование конечности кинетической энергии:

( )2 dD. (1.8) D Отметим, что ограничение (1.8), являясь естественным с физической точки зрения, вполне замыкает математическую постановку задачи.

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

Приведем эти уравнения в векторной форме:

G= n ds + P, M = R n ds + R0 P. (1.9) S1 S Здесь G и M импульс и момент импульса твердого тела;

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

P = (Px, Py, Pz ) импульс внешней ударной силы, приложенный к твердому телу в точке R0 = (x0, y0, z0 ).

Важно отметить, что решение задачи о безотрывном ударе плавающе го тела применимо только в том случае, когда импульсивное давление pt неотрицательно во всей области, занятой жидкостью (формула (1.6)), что в силу классического принципа максимума и минимума для гармониче ской функции равносильно неотрицательности функции pt на смоченной поверхности тела S1. Проверить выполнение последнего условия до ре шения поставленной задачи удается только в некоторых исключительных случаях. Поэтому обычно предположение о безотрывности удара делается априори и проверяется только после решения задачи (1.1)–(1.8), то есть только после нахождения потенциала скоростей. Если условие неотри цательности функции pt (неположительности функции ) выполнено на всей смоченной поверхности тела, то данная математическая модель кор ректна и полностью соответствует задаче о безотрывном ударе плавающего тела. Однако, предположение о безотрывности удара часто приводит к ре шениям с отрицательными импульсивными давлениями в жидкости, что недопустимо физически. В этом случае математическая модель (1.1)–(1.8) оказывается некорректной и строится другая модель, отвечающая задаче об ударе с отрывом. Таким образом, появление отрицательных импуль сивных давлений в жидкости объясняется наличием зон отрыва частиц жидкости от смоченной поверхности плавающего тела.

Дальнейшие рассуждения будем проводить при следующих упрощаю щих предположениях (смотреть п. 1.2). Считаем, что плавающее тело од нородно и имеет три взаимно перпендикулярные плоскости симметрии xz, yz, xy. В этом случае центр тяжести тела совпадает с началом координат, а его главные оси инерции с осями координат. Предположим также, что вся область, занятая жидкостью, обладает двумя взаимно перпендикуляр ными плоскостями симметрии xz и yz. Покомпонентная запись уравнений (1.9) с учетом сделанных предположений симметрии приводит к следую щим соотношениям:

(mT + m) v0 = Pz, (Ix + Jx ) x = y0 Pz z0 Py, (Iy + Jy ) y = z0 Px x0 Pz, Px = µ2 y, Py = µ1 x, x0 Py = y0 Px, (1.10) 1 nz ds, µ1 = 2 ny ds, µ2 = 3 nx ds, m = (1.11) S1 S1 S 2 (ynz zny ) ds, 3 (znx xnz ) ds, Jx = Jy = S1 S где mT, Ix, Iy соответственно масса тела и его моменты инерции отно сительно осей x и y (в дальнейшем считаем их заданными величинами);

коэффициенты присоединенных масс;

через 1, 2, m, Jx, Jy, µ1, µ 3 обозначены единичные потенциалы, соответствующие случаям v0 = 1, x = 0, y = 0, (центральный удар), v0 = 0, x = 1, y = 0 и v0 = 0, x = 0, y = 1.

Таким образом, в линейной постановке, дело сводится к решению сме шанной краевой задачи (1.1)–(1.8), проверке условия безотрывности удара ( 0 на S1 ) и определению коэффициентов присоединенных масс по формулам (1.11).

Подчеркнем, что коэффициенты присоединенных масс m, Jx, Jy, µ1, µ являются главными неизвестными в задаче о безотрывном и вертикальном ударе плавающего тела.

В п. 1.3 излагается специальный асимптотический метод, направленный на решение пространственных задач гидродинамического удара в областях сложной геометрической конфигурации. В его основе лежит предположе ние о том, что стенки бассейна удалены от плавающего тела на большие расстояния. В дальнейшем считаем, что неподвижная твердая граница бас сейна S3 получена в результате гомотетии с центром в начале координат 0 и коэффициентом h некоторой фиксированной поверхности S3 : S3 = hS (x = hx0, y = hy 0, z = hz 0 ). Собственно говоря, речь идет о построении асимптотики решения рассматриваемой смешанной задачи при h.

Идея этого метода своими корнями уходит к трудам классиков, кото рые при исследовании задач о движении в жидкости двух подобных тел, например, двух шаров или двух круглых цилиндров применяли метод по следовательных приближений Стокса. С другой стороны в вычислитель ной математике хорошо известен альтернирующий метод Шварца. В обоих случаях решение исходной задачи для области сложной геометрической конфигурации сводится к последовательному решению задач в областях, имеющих более простые формы границ. Таким образом, поочередно рас сматриваются две краевые задачи в областях D и G случай h = и задача в ограниченном бассейне при отсутствии плавающего тела. При этом каждый раз ликвидируются невязки, возникающие на неподвижной границе S3 и смоченной поверхности тела S1. После разложения получен ных приближений в ряды по степеням h1 и удержания необходимого ко личества членов, приходим к асимптотике для больших значений h. Отме тим, что такой подход позволяет построить асимптотику потенциала ско рости в любой фиксированной (не зависящей от h) окрестности смоченной поверхности тела и на ее основе определить асимптотики основных ха рактеристик удара коэффициентов присоединенных масс. В частности, асимптотика присоединенной массы плавающего тела, являющаяся основ ной характеристикой задачи о центральном ударе, имеет вид ( = 1):

1 (m + V )2 3 (m + V )(c2 1 + c3 2 ) m = m + h+ h+ 2 1 (m + V )3 h + O(h7 ), h.

+ (1.12) Здесь m присоединенная масса твердого тела, плавающего на поверх ности неограниченной жидкости (случай h = );

V объем погруженной части тела.

Важно отметить, что постоянные m, c2, c3, входящие в асимптоти ческую формулу (1.12), зависят только от геометрии плавающего тела, а другие коэффициенты 1, 1, 2 только от формы границы бассейна.

Первая группа характерных величин находится на основании решения 1 задачи об ударе твердого тела, плавающего на поверхности неограничен ной жидкости (случай h =, c1 = m + V ):

x2 z 31 z 3 nz ds 2xznx + x2 nz z 2 nz 1 ds, c2 = S1 S y 2 z 31 z 3 nz ds 2yzny + y 2 nz z 2 nz 1 ds.

c3 = S1 S Для определения второй группы характерных постоянных требуется ре шить краевые задачи в ограниченном бассейне при отсутствии плавающего тела:

fk fk = 0, r G0, = 0, = Qk, k = 1,..., 3, (1.13) fk n S z= z(4x2 y 2 z 2 ) z(4y 2 x2 z 2 ) z Q1 =, Q2 =, Q3 =, n R3 R7 R n n где G0 фиксированная область, ограниченная поверхностью S3 и плос- костью z = 0, R = x2 + y 2 + z 2. После решения краевых задач (1.13) постоянные 1, 1, 2 находятся по формулам:

1 = 35 + 2, 2 = 36 + 3, 3 f1 3 f f1 f2 f 1 =, 2 = 2, 3 = 2, 5 =, 6 =, z x z y z z z где частные производные вычисляются в точке M0 (0, 0, 0).

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

методы граничных интегральных урав нений (методы ГИУ);

вариационные методы.

Аналогичные асимптотики получены для коэффициентов J1, J2, µ1, µ2.

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

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

в противном случае возникает отрыв.

Впервые условие безотрывности удара круглого диска, плавающего на поверхности слоя жидкости конечной или бесконечной глубины, было най дено в статье И. И. Воровича и В. И. Юдовича (ПММ, 1957 г.). Для жидко сти малой глубины этот вопрос подробно изучен в статьях М. И. Чебакова (ПММ, 1974 г.) и Д. Б. Рохлина (Журнал ВМ и МФ, 1998 г.). В насто ящем параграфе дается обобщение результатов статьи И. И. Воровича и В. И. Юдовича на случай произвольной осесимметричной области.

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

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

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

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

В п. 2.1, в качестве первых примеров применения полученных в преды дущей главе формул, рассмотрены задачи удара в областях со сфериче скими границами. Изучены случаи, когда тело плавает на поверхности жидкости, наполняющей полушар и шар, полупогружен в жидкость, на полняющую ограниченный бассейн. Отмечено, что области, имеющие сфе рические границы, выделяются среди других областей, представляя собой некоторый вырожденный случай по отношению к величине погрешности остаточного члена. Именно для этих границ асимптотические формулы принимают наиболее простой вид.

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

В п. 2.3, в качестве нетривиального примера, рассмотрена задача об ударе твердого тела, плавающего на поверхности идеальной несжимае мой жидкости, ограниченной веретенообразной поверхностью вращения.

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

Фиг. 2.1. Твердое тело, плавающее на поверхности жидкости, наполняющей бассейн, ограниченный веретенообразной поверхностью вращения В п. 2.4 исследуется задача о центральном ударе вырожденного тора твердого тела, полученного вращением окружности вокруг своей касатель ной, и наполовину погруженного в идеальную несжимаемую жидкость бесконечной глубины. Для ее решения применяются методы разделения переменных для уравнений Лапласа и функции тока в вырожденных би полярных координатах. Для основных характеристик удара потенциала скоростей на смоченной поверхности вырожденного тора и присоединен ной массы, получены явные выражения.

В п. 2.5 изучена задача о вертикальном ударе кольца, плавающего на по верхности идеальной несжимаемой жидкости, занимающей все нижнее по лупространство. С помощью тороидальных координат данная задача сво дится к решению парных интегральных уравнений, связанных с функци ями Лежандра, а затем к интегральным уравнениям Фредгольма второго рода с гладкими ядрами. Найдены численные значения относительных ко эффициентов присоединенных масс, выведено условие безотрывности уда ра. Отдельно рассмотрены логарифмически-степенные асимптотики тон кого кольца и степенные асимптотики кольца с малым центральным от верстием. В осесимметричном случае асимптотики присоединенной массы кольца перекрываются во внутреннем диапазоне и дают исчерпывающее решение задачи с погрешностью менее, чем 0,5%.

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

Глава III В этой главе рассмотрена задача о вертикальном ударе твердого тора эллиптического поперечного сечения, плавающего на поверхности идеаль ной несжимаемой и неограниченной жидкости. Для решения данной зада чи применяется метод граничных интегральных уравнений (метод ГИУ).

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

В п 3.1 дается постановка задачи и кратко излагается метод граничных интегральных уравнений. Данная задача имеет два существенных пара метра: = a/b, = c/a, где a и c соответственно горизонтальная и вертикальная полуоси эллипса поперечного сечения тора, b расстояние от оси вращения z до центра этого эллипса, 0 1, 0.

В п. 3.2 для решения задачи о вертикальном ударе твердого тора кру гового поперечного сечения применяется метод граничных интегральных уравнений. С помощью параметризации поверхности тора в тороидальных координатах данная задача сводится к решению интегральных уравнений Фредгольма второго рода на конечном промежутке. Причем ядра получен ных интегральных уравнений выражаются через полные эллиптические интегралы I и II рода и имеют логарифмические особенности на диагонали.

В случае тонкого тора построены и строго обоснованы логарифмически степенные асимптотики решений интегральных уравнений. На их основе определены асимптотики коэффициентов присоединенных масс ( 0).

В частности, асимптотика присоединенной массы тора имеет вид:

1 3 ln 2 7 1 m = a3 2 1 + ln 1 + + 3 ln + 4 4 16 21 ln 2 9(ln 2) 3 ln 2 43 3 ln 1 3 + o(3 ).

+ (3.1) 16 256 64 32 Важно отметить, что асимптотика присоединенной массы тора позволя ет проводить численные расчеты во всем диапазоне изменения характер ного параметра задачи (то есть для любых торов, 0 1) с погрешно стью менее одного процента. Такой вывод был сделан на основе сравнения асимптотической формулы (3.1) с численными результатами, приведенны ми в статье T. Miloh, G. Waisman, D. Weils (J. Eng. Math. 1978), а также с присоединенной массой вырожденного тора.

В п. 3.3 асимптотики тонкого тора обобщены на случай тора эллиптиче ского поперечного сечения. Для основных характеристик удара найдены и строго обоснованы явные асимптотические формулы ( 0):

(1 + )2 (1 + )2 32 1 m = a + ln + ln 16 ln(1 + ) + 16 16 + O 2 ln 1, 3(1 + )2 5 ln 1 + O(2 ), J = a 2 (2 2 + 5) 4 µ2 = a + O().

В предельном случае, когда 0, приходим к асимптотикам при соединенной массы и присоединенного момента инерции тонкого кольца, плавающего на поверхности неограниченной жидкости.

В п. 3.4 рассматривается тор произвольного поперечного сечения.

Глава IV В этой главе дается постановка задачи об ударе с отрывом твердого те ла, плавающего на поверхности идеальной несжимаемой жидкости, напол няющей ограниченный бассейн. Особенность данной задачи заключается в том, что область контакта тела с жидкостью (равно как и зона отрыва) заранее неизвестна и подлежит определению вместе с течением жидкости после удара. Вследствие этого поставленная задача является нелинейной и относится к классу задач со свободными границами. Разработаны специ альные аналитические методы решения нелинейных задач гидродинами ческого удара в ограниченных областях. Дано решение ряда плоских сме шанных задач об отрывном ударе пластины, плавающей на поверхности идеальной и несжимаемой жидкости. Рассмотрены случаи, когда область, занятая жидкостью, представляет собой полуплоскость, полосу (отрывной удар пластины, плавающей на поверхности жидкости конечной глубины), усеченную круговую луночку и произвольную ограниченную область.

В п. 4.1 настоящей главы, в классической постановке (Л. И. Седов и дру гие), рассматривается плоская задача об ударе твердого тела, плавающего на поверхности идеальной несжимаемой жидкости, наполняющей ограни ченный бассейн. Предполагается, что в результате удара происходит от рыв частиц жидкости от поверхности тела (удар с отрывом). Вследствие этого на смоченной поверхности плавающего тела образуется зона отры ва, которая в общем случае представляет собой многосвязное множество.

Важно отметить, что зона отрыва (равно как и область контакта тела с жидкостью) заранее неизвестна и подлежит определению вместе с тече нием жидкости после удара, то есть вместе с потенциалом скоростей.

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

Пусть до удара тело и жидкость покоились. Тогда движение жидкости после удара будет потенциальным, причем потенциал скоростей, приоб ретенных частицами жидкости в результате удара, определяется решени ем смешанной краевой задачи теории потенциала с неизвестными априори областями контакта (Л. И. Седов, 1965 г.;

В. И. Юдович, 2005 г.):

= 0 r D, (4.1) P d d vr r j © = Vn, 0 r S11, (4.2) n S2 E x t c3 Vn, =0 r S12, (4.3) tc 1 n t c = 0 r S2, (4.4) y S c = 0 r S3, (4.5) n 0 r, (4.6) Vn = vx nx + vy ny + (ynx xny ). (4.7) Здесь D область, занятая жидкостью (граница области D предпола гается кусочно-гладкой);

S1 = S11 S12 смоченная поверхность тела, причем S11 часть границы на которой не происходит отрыва частиц жидкости, а S12 зона отрыва;

S2 свободная поверхность жидкости;

S3 неподвижная твердая граница бассейна;

n = (nx, ny ) вектор нор мали к поверхности S1, направленный внутрь области D;

Vn проекция на эту нормаль скоростей точек границы тела;

vx, vy, поступательные вдоль осей x и y и угловая скорости, приобретенные телом в результате удара;

r = (x, y). Декартовы координаты x, y введены таким образом, что ось x расположена вдоль линии свободной поверхности жидкости, ось y направлена вертикально вниз вглубь жидкости, начало координат совпа дает с некоторой точкой тела. В рассматриваемых далее частных случаях пластинки и эллипса начало координат выбирается в центре этих тел. Гра ничное условие (4.6) добавляется только в случае неограниченной жидко сти. Импульсивное давление pt связано с потенциалом скоростей соот ношением pt =, где плотность жидкости. После решения задачи (4.1)–(4.6) скорости жидких частиц сразу после удара находятся по фор муле V = grad.

Объясним физический смысл неравенств в (4.2), (4.3). Первое условие означает, что импульсивное давление pt должно быть направлено в сторо ну жидкости, должно сжимать, но не растягивать жидкость. Второе нера венство требует, чтобы в каждой точке смоченной границы твердого тела жидкая частица не входила внутрь твердого тела, хотя ей разрешается отрываться от этой границы. В тех точках границы S1, в которых 0, выполняется равенство = Vn, а в тех точках, где Vn, справедли n n во равенство = 0. Вместе с тем, в некоторых точках границы S1 могут одновременно выполняться равенства = 0 и = Vn. Это, например, n происходит в точках раздела зоны отрыва S12 и зоны безотрывного удара S11.

В статье В. И. Юдовича (Владик. матем. журнал, 2005 г.) доказано, что поставленная задача с односторонними неравенствами на границе имеет единственное решение при дополнительном условии конечности кинетиче ской энергии течения жидкости после удара.

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

В п. 4.2 при помощи метода парных интегральных уравнений, связан ных с тригонометрическими функциями, строится точное решение зада чи об отрывном ударе пластины, плавающей на поверхности идеальной несжимаемой и неограниченной жидкости. Для точки отрыва пластины найдено явное выражение (2a ширина пластины;

vy = v0, vx = 0):

4v0 a c = +.

3 Установлена связь между точкой приложения внешнего ударного им пульса x0 и точкой отрыва пластины c (массой и моментом инерции пла стины пренебрегаем). Получены явные формулы для определения полного ударного импульса и его момента относительно начала координат, подей ствовавшие на пластину со стороны жидкости во время удара.

В п. 4.3 рассматриваются различные способы определения точки отры ва c. Приводится формулировка вариационного принципа Огазо (см., на пример, Лионс Ж. Л. Оптимальное управление системами, описываемыми уравнениями с частными производными, 1972 г.) и дается его применение к решению задачи о вертикальном ударе пластины. Суть данного вариа ционного принципа заключается в том, чтобы вначале построить решение 2 (x, y) соответствующей линейной смешанной краевой задачи в области D с фиксированным разбиением границы S1 на области задания краевых условий типа Дирихле–Неймана (2 и 1 предполагаемые зоны отрыва и контакта, S1 = 1 2 ). После этого, для определения функции (x, y) в произвольной точке области D, необходимо взять точную нижнюю грань функции 2 (x, y) по всем таким разбиениям:

(x, y) = inf 2 (x, y) для всех (x, y) D.

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

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

P -a ca c r E x0 x T H y c c Решение данной задачи состоит из двух основных этапов. Вначале, при помощи метода парных интегральных уравнений, связанных с тригоно метрическими функциями, строится точное решение линейной смешанной краевой задачи в полосе, когда на одной из ее границ имеется отрезок, разделяющий краевые условия первого и второго рода. После этого точка отрыва пластины находится с помощью вариационного принципа Огазо, из которого автоматически следует выполнение граничных условий в ви де неравенств. Отметим, что решение соответствующей линейной задачи о центральном ударе пластинки в случае жидкости конечной глубины, впер вые, было получено М. В. Келдышем, 1935 г.

В результате для определения точки отрыва пластины выведено транс цендентное уравнение:

v0 a c H + f (b) = 0, (4.8) 2 b t2 dt ch b/2 a+c f (b) =, b= ·, K(th b/2) 2 H 2 ch b 2 ch t где K(k) полный эллиптический интеграл первого рода;

H глубина жидкости;

2a ширина пластины;

v0 и поступательная и угловая скорости, приобретенные пластиной в результате удара;

c точка отрыва пластины.

Уравнение (4.8) может быть эффективно решено асимптотически при больших значениях H. Для этого необходимо разложить функцию f (b) в ряд по степеням H 1 и искать решение уравнения (4.8) в виде ряда по таким же степеням. После применения метода неопределенных коэф фициентов, получим для точки отрыва c асимптотическую формулу:

2 3 (a+c )5 H 4 +O(H 6 ), c = c + (a+c ) H H, (4.9) 576 где c точка отрыва в случае неограниченной жидкости.

Для более полного исследования задачи также был найден главный член асимптотики точки отрыва пластины в случае жидкости малой глу бины (H 0):

3v0 a c= +. (4.10) 2 Подчеркнем, что трансцендентное уравнение (4.8) вместе с асимптоти ками для больших и малых H позволяет провести исчерпывающий каче ственный и количественный анализ задачи об отрывном ударе пластины в случае жидкости конечной глубины. Отметим также, что данное урав нение может быть достаточно эффективно решено методом секущих.

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

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

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

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

1 a + c [1 h2 ph3 + qh4 ] + O(h5 ), c = c + h, 3 где h безразмерный параметр, характеризующий удаление стенок бас сейна произвольной формы от плавающей пластины (коэффициент гомо тетии некоторой фиксированной кривой S3 на кривую S3 );

c точка отрыва пластины в случае неограниченной жидкости;

коэффициенты, p, q зависят исключительно от формы границы бассейна.

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

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

Метод нелинейных граничных интегральных уравнений типа Гаммер штейна был предложен Б. А. Галановым (ПММ, 1985 г.) для решения ста тических контактных задач теории упругости с неизвестными заранее об ластями контакта;

контактные задачи для решения которых применялся такой подход приводятся в монографиях В. М. Александрова, Д. А. По жарского (1988 г.) и В. М. Александрова и М. И. Чебакова (2004 г.).

В п. 5.1 настоящей главы дается постановка задачи об отрывном ударе эллиптического цилиндра и излагается способ ее сведения к нелинейному граничному интегральному уравнению типа Гаммерштейна. Далее более подробно остановимся на связи данной задачи с нелинейным интеграль ным уравнением.

Наряду с исходной нелинейной задачей (4.1)–(4.8) рассмотрим линей ную смешанную краевую задачу в области D:

W W W = 0, r D, = f, W = 0, = 0. (5.1) n n S1 S2 S Пусть K: L2 (S1 ) L2 (S1 ) линейный оператор, ставящий в соответ ствие нормальной производной функции W на S1 решение задачи (5.1) на S1 :

K(f ) = W. (5.2) С учетом соотношений (4.2), (4.3) и (5.2) нормальная производная по тенциала скоростей на границе S1 будет определяться в результате ре шения системы неравенств:

u= Vn, n K(u) K(Vn ), u = 0, r S11, (5.3) K(u) = K(Vn ), u 0, r S12. (5.4) Естественно предположить, что существует область S0 = {M : K(Vn ) 0}, K(Vn ) 0 при M S0 (S0 подобласть S1 ). Отметим, что условие K(Vn ) 0 на всей границе S1 соответствует задаче о безотрывном уда ре плавающего тела.

Далее, следуя работе Б. А. Галанова (ПММ, 1985 г.) введем нелинейные операторы v, v + :

v = v (M ) = inf{v(M ), 0}, v + = v + (M ) = sup{v(M ), 0}, и рассмотрим относительно неизвестной функции v нелинейное оператор ное уравнение:

T v = µv K(v + ) K(Vn ), T v = 0, (5.5) где параметр µ может принимать произвольные положительные значения.

Справедливо следующее утверждение.

Если v решение уравнения (5.5), то (u = v +, S12 = {(x, y) S1 :

v (x, y) 0}) решение системы (5.3)–(5.4), причем S12 = при S0 =, обратно, если (u, S12 ) решение системы (5.3)–(5.4), то функция v = µ1 K(Vn ) + u + µ1 K(u), M S решение уравнения (5.5). Область отрыва S12 может быть многосвязной.

Доказательство этого утверждения, а также существование единствен ного решения уравнения (5.5) в пространстве L2 (S1 ) при любом µ проводится в полной аналогии со статьей Б. А. Галанова (ПММ, 1985 г.).

Также объясняется произвольность параметра µ. Если v1 и v2 решения уравнения (5.5), соответствующие значениям µ = µ1 и µ = µ2, µ1 = µ2, + + то v1 = v2.

Таким образом, для решения системы (5.3)–(5.4) достаточно найти ре шение v нелинейного интегрального уравнения (5.4). При этом неизвест ная априори зона отрыва S12 определяется как множество всех точек M S1, для которых выполняется неравенство v 0.

После решения нелинейного интегрального уравнения (5.5) потенциал скоростей на границе S1 определяется по формуле:

= K(v + + Vn ). (5.6) Интересно обратить внимание на то, что в случае безотрывного удара, функция v 0 на S1 и, следовательно, v + = 0. В этом случае формула (5.6) дает решение задачи о безотрывном ударе плавающего тела.

В п. 5.2 описывается построение оператора K для случая жидкости ко нечной глубины. С помощью методов теории потенциала типа простого слоя нахождение оператора K сводится к решению интегрального уравне ния Фредгольма второго рода с гладким ядром и гладкой правой частью и вычислению на основе его решения интеграла, имеющего логарифмиче скую особенность.

В п. 5.3 рассматривается численная реализация метода. Для решения нелинейного интегрального уравнения (5.5) применяется метод М. А. Крас носельского, суть которого заключается в построении последовательных приближений по формулам:

vn+1 = vn (Q v0 ) T vn, n = 0, 1, 2, 3,..., где Q дифференцируемый оператор, хорошо аппроксимирующий опера тор T по равномерной метрике. Отметим, что производная Фреше опера тора T существует не на всех элементах пространства L2 (S1 ) (или C(S1 )), а только на некоторых специальных множествах.

Изучаются итерационные процессы, соответствующие различным на чальным данным. Важно отметить, что при специальном выборе началь ного приближения (v0 (x) c, c 0) решение задачи сводится к явному итерационному процессу:

vn+1 = vn T vn, n = 0, 1, 2, 3,...

µ В п. 5.4 проводится сопоставление численных результатов, полученных при помощи метода нелинейных интегральных уравнений, с известными точными решениями.

В п. 5.5 дается численный анализ задачи. Исследуется влияние дна, а также кинематических параметров и геометрических размеров на обра зующуюся на поверхности эллиптического цилиндра зону отрыва частиц жидкости. В частности, отметим случай, когда в результате удара эллипс начинает двигаться в горизонтальном направлении и вращаться вокруг своей оси (горизонтальный удар с вращением). Полагаем = a/h = 2/3, = b/a = 2, vy = 0, фиг. 5.1. Здесь основное внимание уделяется изучению влияния безразмерной угловой скорости 1 = a/vx на зону отрыва S в случае жидкости конечной глубины. Далее индекс 1 опускаем.

При горизонтальном ударе эллипса ( = 0) на смоченной границе S образуется только одна зона отрыва. По мере увеличения угловой скорости эта зона уменьшается и, начиная с некоторого места, появляется вторая зона отрыва, которая вначале очень быстро расширяется (точки отрыва C и C2 движутся в разные стороны). Затем точка отрыва C2 меняет направ ление. Дальнейшее увеличение угловой скорости вращения приводит к тому, что зона отрыва, расположенная в задней части эллипса исчезает и остается одна зона отрыва, находящаяся строго под водой. При больших ордината точки отрыва C1 становится сколь угодно маленькой, а точка отрыва C2 приближается к точке отрыва, соответствующей крутильному удару эллипса (vx = 0, vy = 0).

Аналогичная ситуация имеет место и в случае, когда эллипс больше вытянут в горизонтальном направлении ( = 0,5, = 2/3, vy = 0, фиг. 5.2).

Интересно обратить внимание на то, что образование двух зон отрыва происходит в случаях a b и a b при различных направлениях вращения эллипса. Отметим также образование двухсвязной зоны отрыва, а также зоны отрыва, расположенной строго в подводной части плавающего тела.

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

vx vx x x c3 c c c h h y y vx v = c1 x x c2 c h h y y Фиг. 5.1. Горизонтальный удар с вращением, b a, vx vx x x c c3 c c h h y y v =0 vx x x c c c2 h h y y Фиг. 5.2. Горизонтальный удар с вращением, b a, В п. 5.6 рассматривается другая постановка задачи, когда заданными следует считать импульс внешней ударной силы P = (Px, Py ) и точку его приложения R0 = (x0, y0 ) (последние величины должны быть связа ны уравнением эллипса). В этом случае скорости vx, vy,, приобретенные эллипсом в результате удара, находятся в результате решения системы нелинейных уравнений уравнений изменения импульса и момента им пульса плавающего тела при ударе. Для решения этой системы применяет ся метод простой итерации, на каждом шаге которого необходимо решить нелинейное интегральное уравнение типа Гаммерштейна.

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

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

Глава VI В этой главе рассматривается задача об ударе с отрывом эллиптиче ского цилиндра, плавающего на поверхности неоднородной несжимаемой жидкости, наполняющей ограниченный бассейн. Подробно исследуется случай экспоненциально стратифицированной жидкости. Данная задача сводится к одному нелинейному операторному уравнению, для решения которого применяется модифицированный метод Ньютона–Канторовича.

Такой подход позволяет одновременно определить неизвестную заранее зо ну отрыва частиц жидкости и течение жидкости после удара. Изучено влияние неоднородности жидкости, а также кинематических параметров и геометрических размеров на зону отрыва частиц жидкости от поверхно сти цилиндра.

В п. 6.1 этой главы рассматривается постановка задачи об отрывном ударе эллиптического цилиндра, наполовину погруженного в неоднород ную, несжимаемую жидкость, наполняющую ограниченный бассейн (В. И. Юдович, 2004 г.):

div = 0, (6.1) · = Vn, 0, r S11, (6.2) n · Vn, = 0, r S12, (6.3) n = 0, r S2, (6.4) = 0, r S3, (6.5) n 0, r, (6.6) Vn = vx nx + vy ny + (ynx xny ), (6.7) 1 ( )2 dD. (6.8) 2 D Здесь = (x, y) известная функция распределение плотности после удара, совпадающее с распределением плотности до удара, им пульсивное давление;

D область, занятая жидкостью;

S1 = S11 S погруженная в жидкость половина эллипса, причем S11 часть границы на которой не происходит отрыва частиц жидкости, а S12 зона отрыва;

S2 свободная поверхность жидкости;

S3 неподвижная твердая граница бассейна;

vx, vy, поступательные вдоль осей x и y и угловая скоро сти, приобретенные эллипсом в результате удара;

n = (nx, ny ) вектор нормали к границе тела, направленный внутрь области D. Условие (6.6) добавляется в случае неограниченной жидкости. Ограничение (6.8) явля ется достаточным и физически оправданным требованием регулярности решения данной задачи вблизи особых точек границы области D. Теоре ма существования и единственности решения задачи (6.1)–(6.8) доказана В. И. Юдовичем (2005 г.).

После решения задачи (6.1)–(6.8) скорости жидких частиц сразу после удара находятся из уравнения: u =, где u приращение скорости в результате удара. Если до удара жидкость покоилась, то функция u есть скорость приобретенная частицами жидкости после удара.

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

Важно подчеркнуть, что область контакта тела с жидкостью (равно как и зона отрыва) заранее неизвестна и подлежит определению вместе с импульсивным давлением. На основании этого поставленная задача является нелинейной и относится к классу задач со свободными граница ми.

Поставленная задача сводится к одному нелинейному операторному уравнению типа Гаммерштейна, для решения которого применяется мо дифицированный метод Ньютона–Канторовича. Такой подход позволяет одновременно определить движение жидкости после удара и неизвестную заранее зону отрыва частиц жидкости.

В п. 6.2 описывается построение оператора K для экспоненциально стратифицированной жидкости конечной глубины. Дело сводится к ре шению интегрального уравнения Фредгольма второго рода с логарифми ческой особенностью в ядре и вычислению интеграла, имеющего также логарифмическую особенность.

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

В п. 6.4 проводится численный анализ задачи. Исследуется влияние неоднородности жидкости на образующуюся на поверхности эллиптиче ского цилиндра зону отрыва частиц жидкости.

Глава VII В этой главе рассматривается существенно пространственная смешан ная задача об отрывном ударе круглого диска, плавающего на поверхности идеальной несжимаемой и неограниченной жидкости. Область контакта тела с жидкостью (равно как и зона отрыва) заранее неизвестна и зави сит от соотношения между поступательной и угловой скоростями, приоб ретенными диском в результате удара. Вследствие этого данная задача является нелинейной и относится к классу задач со свободными грани цами. Для решения поставленной задачи применяется метод нелинейных граничных интегральных уравнений типа Гаммерштейна, развитый в ра ботах Б. А. Галанова (ПММ, 1985;

ДАН, 1987). Такой подход позволяет одновременно находить течение жидкости после удара и неизвестную за ранее зону отрыва частиц жидкости.

В п. 7.1 дается постановка задачи об ударе с отрывом круглого дис ка и описывается ее сведение к нелинейному граничному интегральному уравнению типа Гаммерштейна.

В п. 7.2 рассматривается построение оператора K соответствующей ли нейной смешанной краевой задачи для случая неограниченной жидкости.

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

В заключении диссертации перечислены основные результаты работы.

В приложении собраны наиболее часто используемые в диссертации формулы: разрывные интегралы Вебера и Мелера;

интегральные пред ставления функций Бесселя и Лежандра;

разложения функций Лежандра в гипергеометрические ряды и различные другие соотношения.

Основные публикации автора по теме диссертации 1. Норкин М. В. Удар вырожденного тора о жидкость бесконечной глу бины // Изв. РАН. МЖГ. 1995. № 5. С. 161–165.

2. Норкин М. В. Удар по твердому телу веретенообразной формы, по груженному в жидкость бесконечной глубины // ПМТФ. 1996. Т. 37.

№ 1. С. 36–41.

3. Норкин М. В. О начале отрыва при гидродинамическом ударе по пла вающему телу // Изв. РАН. МЖГ. 1996. № 6. С. 99–104.

4. Норкин М. В. Удар тонкого тора о поверхность идеальной жидкости бесконечной глубины // Журнал вычисл. матем. и матем. физики.

1997. Т. 37. № 10. С. 1263–1268.

5. Норкин М. В. Вертикальный удар по твердому телу, плавающему на поверхности слоя идеальной несжимаемой жидкости конечной глуби ны // Изв. РАН. МЖГ. 1999. № 1. С. 74–81.

6. Норкин М. В. Вертикальный удар тонкого тора эллиптического попе речного сечения, плавающего на поверхности идеальной и несжимае мой жидкости // Изв. РАН. МЖГ. 2000. № 2. С. 144–152.

7. Норкин М. В. Об учете влияния стенок бассейна произвольной формы при безотрывном ударе плавающего тела // ПМТФ. 2001. Т. 42. № 1.

С. 77–81.

8. Норкин М. В. Вертикальный удар твердого тела, плавающего на по верхности идеальной несжимаемой жидкости в ограниченном бассей не произвольной формы // Изв. РАН. МЖГ. 2002. № 3. С. 114–122.

9. Норкин М. В. Отрывной удар пластины, плавающей на поверхности идеальной несжимаемой жидкости в ограниченном бассейне // ПМТФ.

2003. Т. 44, № 4. С. 43–50.

10. Норкин М. В. Методы решения нелинейных задач гидродинамическо го удара в ограниченных областях // Изв. РАН. МЖГ. 2005. № 4.

С. 135–147.

11. Норкин М. В. Смешанные задачи гидродинамического удара. Ростов на-Дону: Изд. ЦВВР, 2007. 136 с.

12. Норкин М. В. Отрывной удар эллиптического цилиндра, плавающего на поверхности идеальной несжимаемой жидкости конечной глубины // Изв. РАН. МЖГ. 2008. № 1. С. 120–132.

13. Норкин М. В. Отрывной удар круглого диска, плавающего на поверх ности идеальной несжимаемой жидкости бесконечной глубины // ПМТФ. 2009. Т. 50, № 4. С. 76–86.

14. Норкин М. В. Удар с отрывом эллиптического цилиндра, плавающего на поверхности несжимаемой, экспоненциально-стратифицированной жидкости // Изв. Вузов. Сев.-Кавк. регион. Естественные науки. Спец выпуск, посвященный 75-летию В. И. Юдовича. Ростов-на-Дону: Изд во СКНЦ ВШ, 2009. С. 168–173.



 




 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.