Расчет поля скоростей в окрестности корабля при его движении, наличии градиентного ветра и качки

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

• •    задача генерации систематического градиентного ветра в приземном слое турбулентной атмосферы (модель приземного слоя атмосферы);

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

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

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

2.    Модель Монина—Обухова приземного слоя атмосферы 3.    Расчет осреднённых полей скоростей в окрестности корабля 4.    Расчет осредненных полей скоростей в окрестности корабля

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

Модель приземного слоя атмосферы характеризуется двумя параметрами: скоростью ветра на заданной высоте (обычно 10 м) и длиной Монина-Обухова Lm, которая характеризует состояние атмосферы (устойчивая, нейтральная, неустойчивая). По величинам этих параметров можно восстановить все необходимые данные, характеризующие приземный слой атмосферы: скорость систематического ветра по высоте, распределение температуры по высоте, а также распределение по высоте параметров турбулентности (уровня q и масштаба L). Детали модели приземного слоя подробно изложены в работах [2], [3].

В качестве примера на рис. 1 показаны профили ветра в сравнении с экспериментом.

Рис. 1. Профили скорости ветра: эксперимент и модель Монина-Обухова [3]

Целью данной работы являлось получение данных о характеристиках вихревого следа в окрестности корабля. Для определения параметров ближнего следа решалась краевая задача для 3D RANS, в качестве модели замыкания использована модель турбулентности к — ш SST. Граничные условия ставились следующим образом: на поверхности корабля задавались условия прилипания; на подстилающей поверхности использовалось условие проскальзывания; на входе в расчетную область задавались компоненты скорости; на выходе и на боковых границах расчетной области задавалось давление. Решение проводилось с использованием компьютерного кода ANSYS CFX 12 (университетская версия). Для дискретизации по пространству использовалась схема второго порядка точности. Установление по времени осуществлялось с использованием неявной схемы второго порядка. Расчёт был выполнен для воздуха при температуре 25 °C.

Математическая модель корабля включала все достаточно крупные элементы конструкции (рис. 2). Моделировались только те части корабля, которые находятся выше ватерлинии. Многоблочная гексаэдральная расчетная сетка состояла из 55 млн элементов (рис. 3).

Рис. 2. Расчётная модель корабля

Исследовалось изменение вихревой обстановки за кораблем в зависимости от профиля ветра.

Рис. 3. Расчётная сетка, па. поверхности корабля

Рассмотрены варианты:

• -    с вертикальным профилем ветра: /inf (Н ) = const,

• -    с нейтральным профилем ветра: /inf (Н ) = 1.1258 • ln(H) + 7.7747.

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

• -    скорость набегающего потока над палубой 25-30 м/с,

• -    скорость корабля была принята равной 10 м/с,

• -    корабль движется навстречу ветру.

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

Основной эффект, который оказывает введение в постановку задачи градиентного профиля ветра, заключается во «всплытии» вихрей за. кораблем (рис. 4).

Рис. 4. Отношение продольной скорости на. линии глиссады к скорости набегающего потока: вверху - вертикальный профиль ветра, внизу - нейтральный профиль

В ходе работы исследовалось влияние килевой качки корабля на форму и характеристики спутного следа за ним, а также влияние самого следа на возможную область посадочного коридора. В качестве исследуемого режима была принята килевая качка корабля с дифферентом на нос в диапазоне от -1.5° до +1.5°, с полным периодом равномерных колебаний Т = 12 с с полным периодом р;звиомерпых колебаний ([фц (Н ) = const = 50 мм), с углами кажущегося ветра 3 = 0°, +7°, -7°, +14°, где отрицательньіе значения углов 3 соответствуют набегающему потоку со стороны угловой палубы корабля. Результаты получены посредством нестационарного решения 3D RANS-уравнений с использованием компьютерного кода ANSYS CFX 12 (университетская версия) на расчетной области длиной 500 м и шириной 600 м за кораблем. Так как прямой расчет нестационарного течения в окрестности корабля при наличии качки требует динамического перестроения сетки, что существенно увеличивает время счета и снижает точность вычислений, в данной работе принята упрощенная постановка задачи. Предполагается, что ввиду малых величин дифферента корабля поле характеристики течения над кораблем близко к стационарным значениям, рассматривалась только область за кормой корабля. При этом на входную границу расчетной области интерполировались значения компонент скорости и давления, полученные в стационарных расчетах компоновки корабля (см. выше), в районе его кормового среза, для соответствующих углов /3- Килевая качка корабля моделировалась смещением интерполированных значений в вертикальной плоскости на входной области по закону: Z — (0.5 • L • sin(g min ) • cos(2?r • t/Т )), г де L - длинна корабля по продольной оси; a_min -угол максимального отрицательного дифферента на корму; 0.5 • L • sin(amin) задано в расчетной модели как const = —3.5 м и соответствует нижнему положению линии кормового среза; Т - период колебаний; Z - исходная координата по высоте для интерполируемых данных при а = 0°. Модель турбулентности и остальные параметры задачи соответствуют принятым для стационарных расчетов. Расчетные исследования показывают, что течение за кораблем с учетом качки заметно изменяется. Наиболее выраженными эти изменения становятся на расстояниях >  250 метров от кормового среза корабля. Предположительно это вызвано сближением и взаимодействием сходящих с элементов конструкции корабля вихрей у подстилающей поверхности. Особенно заметно изменение характера течения при углах 3, отличных от нуля (рис. 5).

Рис. 5. Изоповерхности завихренности за кормовым срезом (вид сверху) при 3 = —7°: вверху -стационарный расчет, внизу - расчет с моделированием килевой качки, t = 2 • Т = 24 с

Как изменяется влияние вихрей спутного следа на характеристики течения на предполагаемой линии глиссады заходящего на посадку летателвного аппарата, с учетом фактора качки, можно видеть на рис. 6. При расчетах с углами /3 ~ — 7° вихри в следе за кораблем расположены параллельно области посадочного коридора и способны оказывать воздействие на большом удалении от корабля (рис. 7). Формирование и эволюция вихрей во времени показана на рис. 8.

Рис. 6. Отношение продольной скорости на. линии глиссады к скорости набегающего потока, при 3 = —7°: вверху - стационарный расчет, внизу - расчет с моделированием килевой качки, t = 2 • Т = 24 с

Рис. 7. Положение вихрей в спутном следе за. кораблем относительно предполагаемой линии глиссады (вид сверху) при t = 24 с

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

Рис. 8. Формирование и эволюция вихрей во времени

6.    Расчет турбулентных порывов скоростей в окрестности корабля

Турбулентность является анизотропной и неоднородной по пространству. Тем не менее если рассматривать достаточно малый объём, то с определенной точностью можно считать, что в нём турбулентность однородна (модель локально-однородной турбулентности).

Для генерирования поля турбулентных порывов {и'’ } будем применять следующий подход [5]. Компоненты пульсаций скорости и- можно представить в виде ряда Фурье:

N й-Д) = 52в-(к„) • cos(knr + yn).                              (1)

П=1

Случайная величина й- задаётся так, чтобы корреляционный тензор (й-Uj) соответствовал известной матрице ковариации. Величина у,, является случайным числом на множестве [0, 2тт) с равномерным распределением вероятности.

Матрица ковариации c-j = (й'Uj) имеет следующий вид:

^cij = ^(й1^) = btj ⁽k) M k,                                     (2)

где M k - элемент объема фазового пространства, a b-j - корреляционный тензор [4]:

6'j(k) = ^ 0« - ^).                     |3>

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

Е (к) = — v ’    27т

(«/г ^L • ^к/2_Т )

Г                                 "I 17 /б

1 + (а/L • к/2^) J где а^ ~ 1.339. Чтобы случайная величина У’ имела заданную матрицу ковариации {cij} достаточно взять случайное число ej с нормальным распределением вероятности, единичной дисперсией и нулевым математическим ожиданием и умножить на результат разложения матрицы ковариации по Холецкому [5]:

H = {Иц } =

V0H 021/Һ11 . О31/И11

V 022 ^{- h} 21 ^(о 32 ^{- И} 21 ^И 31 ) ^/И 22

V ° 33 ^{- И} 21 ^-

После этого искомая величина находится следующим образом: У- = h^j • ej. В каждой точке физического пространства выделим достаточно малый объём так, чтобы в нём считать турбулентность однородной. В этом объёме турбулентность имеет масштаб Ьд(г) и уровень 7д(г). С учетом этого формула (1) принимает следующий вид:

N v‘(r) = 52 У (k_n) • q_A cos(k_n/LA -г + ip_n ).                        (6)

n =1

В качестве примера рассмотрим случай движения корабля, описанный в п. 2. На рис. 9 показан результат генерации поля турбулентных пульсаций в следе за кораблем в связанной (с кораблем) системе координат. Следует отметить достаточно высокий уровень пульсаций, амплитуда которых достигает 4 м/с при скорости ветра над палубой 25 м/с (в связанной системе координат).

Рис. 9. Поле амплитуды (м/с) турбулентных пульсаций в следе за. кораблем

Работа, выполнена, при финансовой поддержке Министерства, образования и науки РФ в рамках договора. № 700013728 от 21.11.2012 «Разработка, моделирующего комплекса, реалистичного восприятия оператором (летчиком) сложных режимов полета и оценки его психофизиологического состояния» по 218 постановлению правительства. РФ.