CC BY
38
2015
НАУЧНЫЙ ВЕСТНИК МГТУ ГА
№ 212
УДК 532.59; 532.527
РАСЧЕТ ВЕТРОВОЙ ОБСТАНОВКИ В ОКРЕСТНОСТИ КОРАБЛЯ
И.С. БОСНЯКОВ, А.А. КОРНЯКОВ, Г.Г. СУДАКОВ
В работе описывается задача расчета возмущенного поля скоростей, возникающего при движении корабля в приземном слое атмосферы, с учетом градиентного ветра и индуцированной турбулентности.
Ключевые слова: надводный корабль, поле скоростей, ветер, турбулентный след, расчет ветровой обстановки.
Введение
При создании системы посадки самолета или вертолета на движущийся корабль возникает задача расчета ветровой обстановки в окрестности корабля. Эта аэродинамическая задача состоит из трех подзадач:
1) генерация систематического градиентного ветра в приземном слое турбулентной атмосферы (модель приземного слоя атмосферы);
2) расчет стационарного (осредненного по времени) возмущенного поля скоростей воздушного потока, возникающего в зоне корабля при его движении в приземном слое атмосферы, которая решается с использованием стационарных уравнений Рейнольдса;
3) расчет индуцированной турбулентности, которая вычисляется с помощью специального алгоритма по известному осредненному полю скоростей и известным турбулентным характеристикам течения (уровню и масштабу турбулентности), которые определяются в рамках двухпа-раметрических моделей турбулентности, используемых для замыкания уравнений Рейнольдса.
Полное поле скоростей получается суперпозицией осредненного и турбулентного полей скоростей (так называемая синтетическая турбулентность). Ниже дано описание методов решения всех трех подзадач, а также пример расчета ветровой обстановки в окрестности конкретного корабля. Оказывается, что градиентный характер ветра в приземном слое атмосферы может оказывать заметное влияние на характеристики течения в окрестности корабля. Имеющиеся в литературе публикации ограничивались лишь случаем вертикального профиля ветра по высоте, а турбулентные характеристики получались с помощью прямого численного моделирования нестационарной задачи [1], что требует значительно больших вычислительных ресурсов, чем предлагаемый подход.
1. Модель Монина-Обухова приземного слоя атмосферы
Модель приземного слоя атмосферы характеризуется двумя параметрами: скоростью ветра на заданной высоте (обычно 10 м) и длиной Монина-Обухова Ьш, которая характеризует состояние атмосферы (устойчивая, нейтральная, неустойчивая). По величинам этих параметров можно восстановить все необходимые данные, характеризующие приземный слой атмосферы: скорость систематического ветра по высоте, распределение температуры по высоте, а также распределение по высоте параметров турбулентности (уровня q и масштаба L). Детали модели приземного слоя подробно изложены в работах [2; 3]. В качестве примера на рис. 1 показаны профили ветра в сравнении с экспериментом.
56
И.С. Босняков, А.А. Корняков, Г.Г. Судаков
Рис. 1. Профили скорости ветра: эксперимент и модель Монина-Обухова [3]
2. Расчет осредненных полей скоростей в окрестности корабля
Целью данной работы являлось получение данных о характеристиках вихревого следа в окрестности корабля. Для определения параметров ближнего следа решалась краевая задача для 3D RANS, в качестве модели замыкания использована модель турбулентности к-ш SST. Граничные условия ставились следующим образом: на поверхности корабля задавались условия прилипания; на подстилающей поверхности использовалось условие проскальзывания; на входе в расчетную область задавались компоненты скорости; на выходе и на боковых границах расчетной области задавалось давление. Решение проводилось с использованием компьютерного кода ANSYS CFX 12 (университетская версия). Для дискретизации по пространству использовалась схема второго порядка точности. Установление по времени осуществлялось с использованием неявной схемы второго порядка. Расчет был выполнен для воздуха при температуре 25°С.
Математическая модель корабля включала все достаточно крупные элементы конструкции. Моделировались только те части корабля, которые находятся выше ватерлинии. Многоблочная гексаэдральная расчетная сетка состояла из 55 млн. элементов (рис. 2).
Рис. 2. Расчетная сетка на поверхности корабля
Исследовалось изменение вихревой обстановки за кораблем в зависимости от профиля ветра. Рассмотрены варианты:
- с вертикальным профилем ветра (H) = const;
- с нейтральным профилем ветра Ц, (И) = 1,1258-1п(Н) + 7,7747.
В качестве исходных данных были использованы профиля приземного ветра, как описано в п. 1. Были приняты следующие параметры течения:
- скорость набегающего потока над палубой 25 - 30 м/с;
- скорость корабля была принята равной 10 м/с (навстречу ветру).
Основной эффект, который оказывает введение в постановку задачи градиентного профиля ветра, заключается во "всплытии" вихрей за кораблем (рис. 3).
б
Рис. 3. Отношение продольной скорости на линии глиссады к скорости набегающего потока: а - вертикальный профиль ветра; б - нейтральный профиль ветра
3. Расчет турбулентных порывов скоростей в окрестности корабля
Турбулентность в общем случае является анизотропной и неоднородной по пространству. Тем не менее, если рассматривать достаточно малый объем, то с определенной точностью можно считать, что в нем турбулентность однородна (модель локально-однородной турбулентности).
Для генерирования поля турбулентных порывов {у[} будем применять следующий подход [4]. Компоненты пульсаций скорости V- можно представить виде ряда Фурье
N
V! (г)=£ V; (кп ). ОС8(кпГ + фп ). (1)
п=1
Случайная величина V- задается так, чтобы корреляционный тензор (у- V^ соответствовал
известной матрице ковариации. Величина фп - является случайным числом на множестве [0,2^) с равномерным распределением вероятности.
Матрица ковариации е^ ^ имеет следующий вид
<И^> = Ь;Дк)Ак, (2)
где Ак - элемент объема фазового пространства; Ьу - корреляционный тензор [5]
Ьй(к ) = М 8а-I
4лк
Л" к2 ,
V к У
58
И.С. Босняков, А.А. Корняков, Г.Г. Судаков
Здесь E(k) - энергетический спектр турбулентной энергии. В данной работе использован кармановский энергетический спектр, который для единичных уровня и масштаба турбулентности можно представить в следующем виде [3]
Е(к) =
55 27п
(а к • ка
1+(°к • к^
17/
(4)
где а,к ~ 1,339. Чтобы случайная величина V имела заданную матрицу ковариации } достаточно взять случайное число ej с нормальным распределением вероятности, единичной дисперсией и нулевым математическим ожиданием и умножить на результат разложения матрицы ко-вариации по Холецкому [4]
( Г~ п Л
Н
к }=
0
С21/Ь11 д/
С31/Ь11 ('
с - Ь2
с22 1 21
0 0
С32 1
Ь21Ь31 )/Ь22
с — ь 2 — Ь 2
с33 Ь31 Ь32
(5)
После этого искомая величина находится следующим образом V[ = Ц • е . В каждой точке
физического пространства выделим достаточно малый объем так, чтобы в нем считать турбулентность однородной. В этом объеме турбулентность имеет масштаб LA(r) и уровень qA(r). С учетом этого формула (1) принимает следующий вид
N
^ (г )= £ V1(к П ^ я А СОз(к JLA • г + фп ) .
(6)
п=1
В качестве примера рассмотрим случай движения корабля, описанный в п. 2. На рис. 4 показан результат генерации поля турбулентных пульсаций в следе за кораблем в связанной (с кораблем) системе координат. Следует отметить достаточно высокий уровень пульсаций, амплитуда которых достигает 4 м/с при скорости ветра над палубой ~25 м/с (в связанной системе координат).
Рис. 4. Поле амплитуды (м/с) турбулентных пульсаций в следе за кораблем
Работа выполнена при финансовой поддержке министерства образования и науки в рамках договора № 700013728 от 21.11.2012 г. "Разработка моделирующего комплекса реалистичного восприятия оператором (летчиком) сложных режимов полета и оценки его психофизиологического состояния" по 218 постановлению правительства РФ.
ЛИТЕРАТУРА
1. Polsky S.A. A computational study of unsteady ship airwake. AIAA 2002-1022. 2002. P. 11.
2. Bobylev A.V., Vyshinsky V.V., Soudakov G.G., Yaroshevsky V.A. Aircraft Vortex Wake and Flight Safety Problems. Journal of Aircraft. Vol. 47. № 2. March-April 2010.
3. Вышинский В.В., Судаков Г.Г. Вихревой след самолета в турбулентной атмосфере // Труды ЦАГИ. 2005. Вып. 2667. С. 155.
4. Гарбарук А.В., Стрелец М.Х., Шур М.Л. Моделирование турбулентности в расчетах сложных течений: учеб. пособие. СПб.: Изд-во Политехн. ун-та, 2012. 88 с.
5. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Гидродинамика. М.: Наука, 1988. Т. 6. 736 с.
AIRCRAFT-CARRIER VORTEX WAKE CALCULATIONS
Bosnyakov I.S., Kornyakov A.A., Soudakov G.G.
Calculations of the perturbed velocity field and induced turbulence due to the motion of the ship in the atmospheric surface layer based on gradient wind are presented in this paper.
Keywords: surface ship, the velocity field, the wind, the turbulent wake, the calculation of wind conditions.
REFERENCES
1. Polsky S.A. A computational study of unsteady ship airwake. AIAA 2002-1022. 2002. P. 11.
2. Bobylev A.V., Vyshinsky V.V., Soudakov G.G., Yaroshevsky V.A. Aircraft Vortex Wake and Flight Safety Problems. Journal of Aircraft. Vol. 47. № 2. March-April 2010.
3. Vihshinskiyj V.V., Sudakov G.G. Vihrevoj sled samoleta v turbulentnoj atmosphere. Trudih CAGI. 2005. Vihp. 2667. P. 155. (In Russian).
4. Garbaruk A.V., Strelec M.Kh., Shur M.L. Modelirovanie turbulentnosti v raschetakh slozhnihkh techeniyj: uchebnoe posobie. SPb.: Izd-vo Politekhn. un-ta. 2012. 88 p. (In Russian).
5. Landau L.D., Lifshic E.M. Teoreticheskaya fizika. Gidrodinamika. M.: Nauka. 1988. T. 6. 736 p. (In Russian).
Сведения об авторах
Босняков Игорь Сергеевич, 1988 г.р., окончил МФТИ (2011), аспирант МФТИ, младший научный сотрудник ЦАГИ, автор 4 научных работ, область научных интересов - атмосферная турбулентность, спутные и вихревые следы за кораблями и самолетами, методы ВАД высокого порядка точности.
Корняков Антон Андреевич, 1984 г.р., окончил МАИ (2007), инженер ЦАГИ, автор 6 научных работ, область научных интересов - математические модели аэродинамики летательных аппаратов в возмущенном потоке, спутные и вихревые следы за кораблями и самолетами.
Судаков Георгий Григорьевич, 1947 г.р., окончил МФТИ (1971), доцент, доктор технических наук, начальник отдела ЦАГИ, ведущий научный сотрудник МФТИ, автор более 100 научных работ, область научных интересов - вычислительная аэрогидромеханика, аэродинамика самолета, струйно-вихревой след.
