Исследование вихревых систем на элементах, обтекаемых несжимаемой жидкостью

Соловьев, Сергей Юрьевич

Исследование вихревых систем на элементах, обтекаемых несжимаемой жидкостью тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Соловьев, Сергей Юрьевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Санкт-Петербург МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

Диссертация по механике на тему «Исследование вихревых систем на элементах, обтекаемых несжимаемой жидкостью»

Автореферат диссертации на тему "Исследование вихревых систем на элементах, обтекаемых несжимаемой жидкостью"

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

¿ЙУ

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

005058277

Соловьев Сергей Юрьевич

ИССЛЕДОВАНИЕ ВИХРЕВЫХ СИСТЕМ НАЭЛЕМЕНГАХ, ОКГЕКАЕМЫХ НЕСЖИМАЕМОЙЖИДКОСГЬЮ

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

АВТОРЕФЕРАТ

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

Санкт-Петербург ] р ¡МД!^ 2013

2013

005058277

Работа выполнена на кафедре гидроаэромеханики математико-механического факультета Санкт-Петербургского государственного университета.

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

профессор Матвеев Сергей Константинович

Официальные оппоненты: доктор технических наук, профессор

Гурьев Юрий Владимирович

(ВУНЦ ВМФ «Военно-морская академия им.

Н.Г. Кузнецова»), заведующий кафедрой «Механика

и гидромеханика»

кандидат технических наук Векслер Валерий Яковлевич (ОАО «Санкт-Петербургское морское бюро машиностроения «Малахит»), ведущий конструктор

Ведущая организация: ОАО «Центральное конструкторское бюро

морской техники «Рубин»

Защита состоится 2013 г. в [Ь часов на заседании

совета Д 212.232.30 по защите диссертаций на соискание ученой степени кандидата наук, на соискание ученой степени доктора наук при Санкт-Петербургском государственном университете по адресу: 198504, Санкт-Петербург, Петродворец, Университетский пр., 28, математико-механический факультет, ауд.

С диссертацией можно ознакомиться в Научной библиотеке им. М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская наб., 7/9.

Автореферат разослан " 2 У" а и2013 г.

405.

Ученый секретарь диссертационного совета

Кустова Е.В.

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

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

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

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

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

скорости в плоскости движителя, обусловленной этими вихревыми

структурами.

Научная новизна работы заключается в следующем:

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

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

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

4. разработаны рекомендации, позволяющие снизить неоднородность потока, натекающего на движитель.

Результаты, выносимые на защиту сводятся к следующему:

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

2. На основании опытных данных дана количественная оценка диссипации этих вихрей при их распространении вниз по потоку, а также установлено их количественное влияние на окружную неоднородность поля скорости в плоскости расположения движителя.

3. Найдены координаты продольных осей подковообразных вихрей от стабилизаторов. Из экспериментальных данных следует, что они вносят основной вклад в неоднородность поля скорости в плоскости расположения движителя.

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

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

6. Расчетными методами найдена форма обтекателя в виде прикорневой наделки к ВЧ, уменьшающего вихреобразование. Наиболее эффективной является форма обтекателя в плане в виде параболы. Максимальная толщина продольных сечений обтекателя должна быть меньше, чем максимальная толщина ВЧ. При этом отношение длины обтекателя к высоте х/у =3.

Достоверность полученных результатов диссертационной работы обеспечивается совместным использованием методов физического эксперимента и численного моделирования. Все экспериментальные исследования проводились на поверенном оборудовании ФГУП «Крыловский государственный научный центр» с действующими сертификатами о метрологической аттестации. Результаты численного моделирования показали удовлетворительное соответствие известным экспериментальным данным.

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

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

Апробация работы. Основные результаты работы докладывались на российских и международных конференциях:

1. X международная конференция по гидродинамики 1СНБ (Санкт-Петербург, 2012);

2. XII всероссийская школа-конференция молодых ученых «Актуальные вопросы теплофизики и физической гидрогазодинамики» (Новосибирск, 2012);

3. XXX и XXIX отраслевая научно-техническая конференция ОАО «Концерн «Морское подводное оружие - Гидроприбор» (Санкт-Петербург, 2011,2010);

4. X и IX научно-техническая конференция «Взгляд в будущее», ЦКБ морской техники «РУБИН» (Санкт-Петербург, 2012,2011);

5. XI международная конференция «Оптические методы исследования потоков» (Москва, 2011);

6. Международная конференция ICMAR «Методы аэрофизических исследований» (Новосибирск, 2010);

7. XXII международный семинар «Струйные, отрывные и нестационарные течения» (Санкт-Петербург, 2010).

Публикации по теме диссертации. Основные результаты работы изложены в 17 научных публикациях, 4 из них опубликованы в изданиях, рекомендованных ВАК РФ. Диссертант является соавтором патента на изобретение.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы из 81 наименования. Работа изложена на 134 страницах текста, включая 84 рисунка.

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

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

В нашей стране обстоятельные работы по исследованию полей скорости в плоскости движителя подводного аппарата выполнены В.А. Тюшкевичем, А.И. Короткиным, В.М. Котловичем в ЦНИИ им. акад. А.Н. Крылова, а также A.B. Жариновым в ЦАГИ им. Н.Е. Жуковского.

Среди работ, выполненных с использованием современных пакетов программ для численного моделирования, следует выделить работы Bull Р., Ballio F., Devenport W., Fureby С., Praisner Т., Chao-Ho Sung, Pankajakshan. R. В области численного моделирования обтекания подводного аппарата общая картина сводится к следующему: многим исследователям удалось достаточно точно рассчитать коэффициенты давления и трения, а также профили средней скорости на поверхности подводного аппарата, однако при расчете полей скорости в кормовой оконечности имеются существенные рассогласования с экспериментальными данными.

Среди российских исследователей по численному моделированию в судостроительной промышленности выделяются работы A.C. Воробьева, М.П. Лобачева, H.A. Овчинникова, A.B. Пустошного, И.В. Ткаченко, И.А. Чичерена, которые посвящены как методическим аспектам (правильному выбору моделей турбулентности), так и расчетам прикладных задач гидромеханики.

Во второй главе приведено описание экспериментальных установок и методов, использованных при физическом и численном моделировании.

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

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

Третья глава посвящена экспериментальному исследованию возникновения и развития вихревых систем вокруг ВЧ, расположенных на модели подводного аппарата. Для этого в аэродинамической трубе проводилось измерение полей полного вектора скорости вокруг модели конструктивно подобной современному подводному аппарату (рисунок 1). Число Рейнольдса Re = 6-106, за характерный размер принято расстояние (L0) от передней кромки ограждения рубки (ОР) до плоскости расположения движителя. Для нахождения оси вихря и оценки его интенсивности в каждом сечении на основе экспериментальных данных строились поля

завихренности со = rot U (рисунок 2). Положение вдоль корпуса оси одной из ветвей подковообразного вихря и концевого вихря от ОР, а также их интенсивности изучены вплоть до плоскости движителя. График диссипации подковообразного и концевого вихрей в зависимости от продольной координаты приведен на рисунке 3. Из приведенного сравнения следует, что диссипация концевого вихря происходит на порядок быстрее, чем подковообразного. В отличие от концевого, подковообразный вихрь не исчезает, а доходит до плоскости расположения движителя и вносит определенный вклад в неоднородность потока, хотя в плоскости движителя он слабее более чем на порядок по сравнению с первым измеренным сечением.

Рисунок 1 - Схема модели и положение сечений измерений около ее борта: а) вид сбоку; б) вид с кормы.

Экспериментальные данные по полям скорости, измеренные в нескольких сечениях за ОР, показывают, что к сечению 2/Ь0 = 0,68 скорость за ней восстановилась до значения и = 0,99, поэтому возмущения от вязкого следа за ОР не влияют на поле скорости в плоскости движителя.

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

1.0

0.8

0.6

* N N

З 0.4

0.2

0.0

0.0 0.2 0.4 0.6 0.8 1.0

Рисунок 3 - Зависимость изменения относительной интенсивности (ш*2 -интенсивность в первом измеренном сечении) концевого и подковообразного вихрей от продольной координаты. На рисунке цифрами обозначено: 1 - начало ОР, 2 - конец ОР, 3 - плоскость движителя

Для определения влияния ОР на неоднородность потока в плоскости движителя проведено измерение поля скорости за моделью схематизированного подводного аппарата в сечение с продольной координатой 2/Ь0 = 1 при наличии и отсутствии ОР (рисунок 4). На рисунке 5 приведены экспериментальные данные измерений на относительных радиусах Я = 0,6; 0,9.

Г 1 2

0° 0°

Рисунок 4 - Система координат и схемы моделей с ОР и без ОР, испытанных в аэродинамической трубе. Вид с кормы

• Диссипация кон йеною вихря —•— Диссипации иодконообра шок» пихря |

к >

\ ■ 1

к \ \

.........ж;

Т......"ж: #

1 О

—

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

где итах(1) , ит(П(к) , - максимальное и минимальное безразмерное значение скорости на данном относительном радиусе, ити1(Н.), - среднее значение скорости на данном относительном радиусе.

Разница между значением коэффициента осредненной по диску неоднородности поля скорости при наличии и отсутствии ОР составляет АХ= 15%. Ось продольной ветви подковообразного подпорного вихря от ОР проходит вблизи концевого радиуса движителя. Поэтому наибольшее значение прибавки = 28% в неоднородность поля скорости приходится

на радиус И = 0,9.

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

радиусах: а) И = 0,6; б) К = 0,9. Цифрами на рисунке обозначено: 1-модель с

ОР; 2-моделъ без ОР

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

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

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

е° а) ©° б)

Рисунок 6 — Значения скорости в плоскости движителя на относительном радиусе а) Я = 0,7; б) Я = 1. На рисунке цифрами обозначено положение горизонтальных стабилизаторов: 1 — выше линии вала; 2 — по линии вала;

3 — ниже линии вала

0.04

0.03

0.01

Рисунок 7 -а) коэффициенты неоднородности потока на различных относительных радиусах; б) коэффициенты а„ряда Фурье при разложении скорости на относительном радиусе /? = 1,0. На рисунках цифрами обозначено положение горизонтальных стабилизаторов: 1 - выше линии вала; 2 -по линии вала; 3 - ниже линии вала

Четвертая глава посвящена количественной оценке влияния основных параметров ВЧ на интенсивность подковообразных вихрей и их влияние на поле скорости в дальнем следе.

Для решения этой задачи проведена серия экспериментов в гидродинамической трубе и численное моделирование обтекания ВЧ, расположенных на плоской пластине. Расчеты выполнялись с использованием метода конечных объемов для решения трехмерных нестационарных уравнений Навье-Стокса (URANS) с привлечением модели турбулентности SST (модель переноса касательных напряжений). Расчеты проводились в ФГУП «Крыловский государственный научный центр» при помощи лицензионного прикладного вычислительного пакета Star ССМ+. Построение расчетной блочно-структурированной сетки осуществлялось в этом же пакете. Во всех расчетах число Рейнольдса, посчитанное по длине выступающей части Re > 105.

В работе получены расчетные данные об интенсивности подковообразного вихря при его распространении вниз по потоку для некоторых вариантов передней кромки ВЧ и обтекателей в виде прикорневой наделки. Из приведенных на рисунке 8 зависимостей следует, что применение эллиптических обводов передней кромки ВЧ позволяет несколько снизить интенсивность вторичных течений перед ней. Использование обтекателей с соотношениями длины к высоте х/у = 2 и х/у = 3 дает существенное (на 26% и 54% соответственно) уменьшение интенсивности подковообразного вихря в месте его зарождения. Также значительно быстрее диссипируют продольные ветви подковообразных вихрей, на расстоянии 22В (В — максимальная толщина выступающей части) для обтекателя с соотношением х/у =2 и 16В для обтекателя с соотношением х/у = 3. Полученные данные полезны проектантам подводных аппаратов. На ранних стадиях проектанту достаточно знать лишь два основных параметра: максимальную толщину ВЧ и расстояние от передней кромки до плоскости движителя, чтобы оценить интенсивность подковообразного вихря, который может попасть в область расположения движителя. Одновременно с этим легко провести анализ того, насколько необходимо изменить переднюю кромку или сместить ВЧ для уменьшения интенсивности вторичных структур, а, следовательно, и неоднородности по окружной координате поля скорости в плоскости расположения движителя.

-а/Ь=1

---а/Ь=2

—■— а/Ь=3

X/В б)

Рисунок 8 - Изменение интенсивности подковообразного вихря при распространении вниз по потоку: а) для ВЧс эллиптической формой передней кромки; б) для ВЧ с обтекателем в виде прикорневой наделки

В заключении кратко сформулированы основные результаты

диссертационной работы, которые сводятся к следующему:

2. На основании опытных данных дана количественная оценка диссипации этих вихрей при их распространении вниз по потоку, а также установлено их влияние на окружную неоднородность поля скорости в плоскости расположения движителя.

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

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

6. На основе результатов численного моделирования оптимизирована форма обтекателя в виде прикорневой наделки к ВЧ, уменьшающего вихреобразование. Наиболее эффективной является форма обтекателя в плане в виде параболы. Максимальная толщина продольных сечений обтекателя должна быть меньше, чем максимальная толщина ВЧ. При этом отношение длины обтекателя к высоте х/у ~3.

Основные результаты диссертации опубликованы в работах:

1. Соловьев С.Ю. Экспериментальное исследование вихревых структур, влияющих на поле скорости в плоскости движителя подводного аппарата // Сборник тр. «Гидродинамика и гидроакустика

перспективных движителей», выпуск 73 (357) ФГУП «ЦНИИ им. акад. А.Н. Крылова» - СПб, 2012. С. 105-113. (Из перечня ВАК)

2. Воробьев А.В., Гузеев А.С., Соловьев С.Ю. Распространение вихревых систем, возникающих на корпусе судна и его элементах И Морской вестник №4(40), СПб, 2012. (Из перечня ВАК)

3. Гузеев А.С., Соловьев С.Ю. Исследование структуры потока вблизи крыла, расположенного в пограничном слое пластины И Морской вестник. 2010. №1. С. 102-103. (Из перечня ВАК)

4. Гузеев А.С., Короткий А.И., Соловьев С.Ю. Об одном из методов определения числа Струхаля при обтекании контуров с фиксированными точками отрыва потока // Морской вестник. 2012. №2. С. 97-100. (Из перечня ВАК)

5. A.S. Guzeev, A.I. Korotkin, S.U. Soloviev. Wind flowing around buildings and bridges // 10-th International Conference on Hydrodynamics, 2012.

6. Соловьев С.Ю. Экспериментальное исследование вторичных течений, вызывающих неоднородность поля скорости в плоскости движителя подводного аппарата // Взгляд в будущее: тр. X молодежной науч.-техн. конф., ОАО «РУБИН» - СПб, 2012. С. 60-75.

7. Соловьев С.Ю. Исследование вихревого следа за крыловым профилем // Тр. XII Всероссийской школы-конференции молодых ученых «Актуальные вопросы теплофизики и физической гидрогазодинамики», Новосибирск, 2012. С. 45-47.

8. Соловьев С.Ю. Исследование структуры подковообразного вихря около выступающей части подводного // Тр. XXX отраслевой науч.-техн. конф. ОАО «Концерн «Морское подводное оружие - Гидроприбор» -СПб, 2011.С. 98-100.

9. Соловьев С.Ю. Влияние числа Рейнольдса и толщины препятствия на положение оси подпорного вихря // Взгляд в будущее: тр. IX молодежной науч.-техн. конф., ОАО «РУБИН» - СПб, 2011. С. 87-101.

10. Guzeev A.S., Soloviev S.U., Vorobyev A.S. A tip vorteces spreading downstream the wing attached to the flat surface // Optical Methods of Flow Investigation: 11-th International Conference - Moscow, 2011. С. 112-114.

11. Solovyev S.U., Vorobyev A.S., Gouzeyev A.S. Downstream vortices speading near the axisymmetric body // International conference on the methods of

aerophysical research «ICMAR» / Edit by acad. Fomin V.M. Novosibirsk, 2010. Page. 238-239.

12. Гузеев A.C., Соловьев С.Ю. Взаимодействие круглой струи и сносящего потока // Струйные, отрывные и нестационарные течения: тр. XXII международного семинара - СПб, 2010. С. 211-212.

13. Гузеев A.C., Соловьев С.Ю. Экспериментальное исследование обтекания выступающих частей судов // Тр. XXIX отраслевой науч.-техн. конф. ОАО «Концерн «Морское подводное оружие - Гидроприбор» - СПб, 2010. С. 87-90.

14. Соловьев С.Ю. Экспериментальное исследование обтекания выступающих частей судов И Моделирование явлений в технических и гуманитарных науках: тр. пятой студенческой конференции СПб ГМТУ -СПб,2010, с.44-47.

15. Соловьев С.Ю. Экспериментальное исследование обтекания выступающих частей судов // Молодежь. Техника. Космос: тр. второй общероссийской молодежной науч.-техн. конф. СПб, 2010. с. 136-139.

16. Соловьев С.Ю. О причинах образования вихрей около подводных аппаратов// IV Уткинские чтения: тр. науч.-техн. конф. СПб, 2009. Т. 3, с. 98-102.

17. Соловьев С.Ю. Экспериментальное исследование обтекания рулевой колонки потоком жидкости // IV Уткинские чтения: тр. науч.-техн. конф. СПб, 2009. Т. 2, с. 101-103.

В совместных публикациях автору принадлежат результаты экспериментов, проведенных в аэродинамической, гидродинамической трубах, а также результаты расчетов. A.C. Гузеев участвовал в "постановке и проведение исследований в гидродинамической трубе. A.C. Воробьеву принадлежит численное моделирование концевых вихрей на различных крыловых профилях. А.И. Короткин выступал в роли научного консультанта по отдельным вопросам.

Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Соловьев, Сергей Юрьевич, Санкт-Петербург

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

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

04201356514

Соловьев Сергей Юрьевич

ИССЛЕДОВАНИЕ ВИХРЕВЫХ СИСТЕМ НА ЭЛЕМЕНТАХ, ОБТЕКАЕМЫХ НЕСЖИМАЕМОЙ ЖИДКОСТЬЮ

Специальность 01.02.05 - «Механика жидкости, газа и плазмы» диссертации на соискание ученой степени кандидата физико-математических наук

Научный руководитель д.ф.-м.н., проф. С.К. Матвеев

Санкт-Петербург 2013

СОДЕРЖАНИЕ

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ..............................................................................4

ВВЕДЕНИЕ..............................................................................................................6

ГЛАВА 1. СОВРЕМЕННОЕ СОСТОЯНИЕ ВОПРОСА...................................10

ГЛАВА 2. ОПИСАНИЕ МЕТОДОВ ПРОВЕДЕНИЯ ИССЛЕДОВАНИЙ.....19

2.1 Описание экспериментальных установок и методик, использованных при физическом моделировании.......................................................................19

2.1.1 Аэродинамическая труба и методы проведения измерений.............19

2.1.2 Гидродинамическая труба и методы проведения визуализации .....24

2.2 Описание математических средств, использованных при численном моделировании....................................................................................................29

2.2.1 Описание математической модели.......................................................29

2.3 Обоснование методов проведения исследований.....................................38

ГЛАВА 3. ВОЗНИКНОВЕНИЕ И РАЗВИТИЕ ВИХРЕВЫХ СИСТЕМ ВБЛИЗИ ВЫСТУПАЮЩИХ ЧАСТЕЙ, РАСПОЛОЖЕННЫХ НА ТЕЛЕ ВРАЩЕНИЯ...........................................................................................................41

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

3.2Исследование подковообразных вихрей, возникающих около ограждения рубки...............................................................................................44

3.3 Исследование концевых вихрей и вязкого вихревого следа от ограждения рубки...............................................................................................60

3.4 Исследование вторичных течений от палубной надстройки...................75

3.5 Исследование подпорных вихрей, возникающих около стабилизаторов...................................................................................................76

3.6 Исследование взаимодействия продольных ветвей подковообразных вихрей, возникающих около стабилизаторов..................................................85

ГЛАВА 4. ИССЛЕДОВАНИЕ СПОСОБОВ УМЕНЬШЕНИЯ ВИХРЕОБРАЗОВАНИЯ ВБЛИЗИ ВЫСТУПАЮЩИХ ЧАСТЕЙ, РАСПОЛОЖЕННЫХ НА ТЕЛЕ ВРАЩЕНИЯ................................................103

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

4.2 Экспериментальное исследование влияния параметров выступающей части и набегающего потока на формирование подковообразного вихря. 105

4.3 Численное исследование способов уменьшение интенсивности образования вихрей перед препятствием, расположенном в пограничном слое пластины ..................................................................................................111

4.3.1 Численные аспекты задачи ................................................................111

4.3.2 Влияние формы передней кромки ВЧ на интенсивность подковообразного вихря..............................................................................118

4.3.3 Уменьшение интенсивности подковообразного вихря перед ВЧ с помощью применения обтекателя..............................................................120

ЗАКЛЮЧЕНИЕ....................................................................................................125

ЛИТЕРАТУРА.....................................................................................................127

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ

X, У, Z - оси декартовой системы координаты;

и — скорость потока в данной точке;

ио — скорость невозмущенного потока;

и = и/и0 ~ относительное значение скорости;

их - проекция вектора скорости на ось X;

иу - проекция вектора скорости на ось У;

и2 - проекция вектора скорости на ось Ъ\

а _ угол скоса потока в вертикальной плоскости, угол атаки;

угол скоса потока в горизонтальной плоскости, угол скольжения; Яе - Число Рейнольдса;

Вп — полуширина палубы модели;

^ расстояние от передней кромки ограждения рубки до

0 плоскости расположения движителя на модели;

О - полярная координата в плоскости расположения движителя;

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

И = г/г0 — относительная величина радиальной координаты; п - номер гармоники;

. расстояние от передней кромки модели до оси

подковообразного вихря; ^ коэффициент неоднородности поля скорости на

к соответствующем радиусе;

5 - толщина пограничного слоя;

„ максимальная толщина препятствия, установленного на

В —

гладкой пластине; А — размер ячейки расчетной сетки;

безразмерное расстояние от расчетной точки до ближайшей твердой стенки

Аббревиатуры и сокращения: ВЧ~ выступающая часть; ОР - ограждение рубки; ПС - пограничный слой;

PIV (Particle Image Velocimetry) - измерение скорости по фотографиям трассеров, запущенных в поток (экспериментальная методика); LDV {Laser Doppler Velocimetry) - лазерная доплеровская анемометрия (экспериментальная методика);

RANS (Reynolds Average Navier Stokes) - осредненные по Рейнольдсу уравнения Навье-Стокса

LES (Large Eddy Simulation) - моделирование крупных вихрей (расчетная методика);

DES (Detached Eddy Simulation) - моделирование отсоединенных вихрей (расчетная методика);

SST (Shear Stress Transport - перенос сдвиговых напряжений) - SST-версия модели турбулентности Ментера

Примечание: остальные обозначения вводятся и поясняются в тексте диссертации.

ВВЕДЕНИЕ

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

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

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

Основные цели настоящей диссертационной работы можно сформулировать следующим образом:

• исследование траектории распространения, интенсивности и диссипации вихревых структур вокруг подводного аппарата: подковообразного и концевого вихря от ограждения рубки (ОР), подпорных вихрей от стабилизаторов, вязких вихревых следов от ОР и стабилизаторов, вихрей от палубной надстройки;

• оценка вклада каждой из этих вихревых структур в неоднородность поля скорости в плоскости расположения движителя;

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

Во введении обосновывается актуальность выбранной темы

исследования, ее научное и практическое значение, сформулированы

основные задачи работы.

Первая глава посвящена современному состоянию вопроса, приведен

обзор ранее выполненных экспериментальных и численных исследований.

Во второй главе приведено описание экспериментальных установок и

методов, использованных при физическом и численном моделировании. В

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

аэродинамической трубе вокруг модели измерялись поля полного вектора

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

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

корпуса подводного аппарата. Также во второй главе приведено описание

использованной в расчете модели турбулентности Ментера (БЗТ) и методов

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

Уделено внимание обоснованию методики проведения исследований.

Третья глава посвящена экспериментальному исследованию

возникновения и развития вихревых систем вокруг выступающих частей,

расположенных на теле вращения. Найдены положение вихрей вдоль

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

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

скорости в плоскости расположения движителя.

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

ГЛАВА I. СОВРЕМЕННОЕ СОСТОЯНИЕ ВОПРОСА

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

Неоднородность по окружной координате поля скорости в плоскости движителя обусловлено, в основном, вязкостными эффектами. Они связаны либо с особенностями развития пограничного слоя на корпусе и вязкими гидродинамическими следами от потерь энергии потока в слоях за выступающими частями, либо с формированием свободных вихревых систем при обтекании элементов архитектуры подводного аппарата. Такие вихри образуются в результате проявления несущих свойств корпуса и выступающих частей, особенно крыльев под углом атаки, а также при локальных пространственных отрывах. В последнем случае одними из наиболее важных являются пары вихрей взаимно противоположных по знаку, представляющие собой свободные ветви подковообразных вихрей (рисунок 1.2), возникающих из-за отрыва пограничного слоя корпуса под воздействием сильного положительного градиента давления перед выступающими частями, расположенными на корпусе. Если смотреть на сформированный подпорный вихрь сверху, то он для всех препятствий имеет подковообразную форму, охватывает препятствие и распространяется в его следе на значительные расстояния.

Параметры течения в зоне такого отрыва и интенсивность

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

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

Рисунок 1.1 - Причины окружной неоднородности поля скорости в плоскости движителя подводного аппарата

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

\ / \ /

\ /

V _

г "о ~

Рисунок 1.3 — Схема распределения продольной скорости II на расстоянии И от стенки в исходном ПС (сплошная линия) и при влиянии на него вязкого вихревого следа от препятствия (штриховая линия)

У / н Ч

1 V' 4 У

11 ч ч / N ✓ Ч _ х \ ✓ ч ✓ V __ /

1 1

-ъ о

Рисунок 1.4 — Схема распределения продольной скорости и на расстоянии И от стенки в исходном ПС (сплошная линия) и при влиянии на него пары продольных вихрей (штриховая линия)

распределения скорости, который лишь в принципе связан с перепадом экстремальных значений, поскольку периодические силы и звук вращения, генерируемые n-лопастным гребным винтом, обуславливаются не всеми гармоническими составляющими рассматриваемого поля, а лишь номерами а-п и a-n ± 1, где а = 1,2...

В последний период, в связи с успехами, достигнутыми в деле снижения отмеченных выше причин шумности, внимание проектантов было привлечено к проблеме широкополосного излучения движителей, соответствующая часть спектра которого определяется уровнем турбулентных пульсаций. Прогресс в этой области теории корабля был достигнут в результате целенаправленных усилий многих специалистов ЦНИИ им. акад. А.Н. Крылова. Среди них следует отметить Ю.С. Базилевского, М.В. Бегака, Л.Д. Волкова, A.C. Горшкова, J1.H. Дубенского, В.М. Железного, О.П. Орлова, A.B. Пономарева, А.Ф. Пустошного, Б.Е. Тощева, С.П. Чекалова, этот перечень является далеко не полным.

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

Измерением полей скорости в кормовой оконечности крупномасштабных моделей подводных аппаратов в опытовом бассейне

ЦНИИ им. акад. А.Н. Крылова занимался В.М. Котлович.

В ЦАГИ подобными работами занимался A.B. Жаринов, особого внимания заслуживает его совместная работа со специалистами ЦНИИ им. акад. А.Н. Крылова. В ней было проведено три цикла испытаний по измерению поля скорости в плоскости дв