Опыт использования FlowVision на ТАНТК им. Г.М. Бериева для моделирования различных гидродинамических задач. - презентация

1 Опыт использования FlowVision на ТАНТК им. Г.М. Бериева для моделирования различных гидродинамических задач

2 Далее представлены результаты экспериментов: 1.Моделирование плавания объёмного тела. 2.Установление уровня жидкости в баке с переборкой. 3.Обтекание плоскокилеватого днища. 4.Сброс жидкости из бака в набегающий воздушный поток.

3 Моделирование плавания объёмного тела. Используется математическая модель со свободной поверхностью (Free Surface). Включены уравнения Навье-Стокса (velocity) и уравнение переноса функции заполнения VOF (метод отслеживания свободной поверхности), уравнение турбулентности не включено, т.е. вихреобразование не учитывалось в данном эксперименте. В модели рассчитывается только жидкая фаза. Учитывается гравитация, g=9,81 кг/м2. Плотность жидкости составляет 1000 кг/м3, температура 273 К, атмосферное давление Па.

4 Моделирование плавания объёмного тела. Геометрия задачи – расчётная область и тело.

5 Моделирование плавания объёмного тела. Объектом исследования выбран параллелепипед. Размеры: 0,2х0,05х0,1 метра (XYZ). Объём тела 0,001м3 (1 литр). Масса тела составляет 0,5 кг. Массовый центр тяжести тела расположен ниже геометрического центра тяжести на 0,01 м (см. рис.3). Массовый центр тяжести совпадает с центром вращения тела. Моменты инерции тела (относительно массового ц.т.): Ixx= , Iyy= , Izz= ,кг м2. Размеры бассейна (расчетной зоны) (размеры): 1,8х0,4х0,9 метров (XYZ). Уровень воды в бассейне 0,2 метра (от дна). Исследуемое тело расположено по центру бассейна. Начальное положение тела относительно границы раздела вода- воздух: геометрический центр тяжести тела находится на границе раздела сред, углы относительно осей – 0 градусов.

6 Геометрический центр тяжести Массовый центр тяжести Изолинии функции VOF (уровень свободной поверхности) Моделирование плавания объёмного тела. Геометрия модели – центр тяжести и свободная поверхность.

7 Граничные условия: на боковых стенках задано условие выход (Full Outlet), на поверхности тела – условие Стенка (Wall) («непротекание»), на остальных поверхностях – условие симметрия (Symmetry). Расчётная сетка: Составила элементов, из них расчётных , по поверхности тела задано измельчение сетки (адаптация с уровнем 2), разбиение начальной расчётной сетки по XYZ: 16х20х21. Моделирование плавания объёмного тела.

8 Расчётная сетка с измельчением по границе исследуемого тела.

9

10 Расчёт производился на ПК с процессором Pentium 1,5 ГГц, объём оперативной памяти RAM 1 Гб. Было просчитано итераций, расчётное время составило 32,953 сек., Время счёта составило 98 часов 42 минуты. Метод расчёта – неявный, алгоритм «удаление мелких частиц» был отключён, CFL=1, ограничение по максимальному шагу итерации 1 секунда. Результаты расчёта анализировались в постпроцессоре FlowVision. На основе анализа данных файла, созданного по супер-группе тела и слою Characteristics HydrostaticPressure, были построены графики нагрузок и перемещений плавающего тела. Моделирование плавания объёмного тела.

11 Результаты расчёта Далее представлены: График зависимости сил всплывания (выталкивающей силы) от времени. График зависимости сил всплывания от времени (увеличен масштаб по оси сил). График зависимости вертикальных перемещений по времени (перемещения массового центра тяжести). График изменения угла дифферента (продольные колебания тела) от времени. График изменения моментов вокруг основных осей (относительно массового центра тяжести) от времени. Учитывая, что геометрический центр тяжести тела был расположен выше уровня воды 10 мм (0,01 м) тело стало совершать вертикальные затухающие колебания. В результате отклонение от нулевого уровня составило 0, 52 мм. Моделирование плавания объёмного тела.

12 Результаты расчёта Конечные усреднённые (среднее арифметическое на текущий шаг) значения величин составили: Выталкивающая сила (усредненная) Н ; Положение центра тяжести тела по вертикали относительно равновесного положения (границей раздела сред) y= - 0,00052 м (отсчёт идёт от уровня раздел сред); Угол дифферента : градусов. Теоретическое значение выталкивающей силы Равновесное положение тело (в системе координат тела): т.е. соответствует координате 0 в абсолютной системе координат (относительно расчётной области). Моделирование плавания объёмного тела.

13 График изменения от времени вертикальных перемещений тела в абсолютной системе координат (0 – граница раздела сред)

14 Моделирование плавания объёмного тела.

15 График изменения от времени вертикальных перемещений тела в абсолютной системе координат (увеличен масштаб) Моделирование плавания объёмного тела.

16 График изменения от времени угла дифферента (точка вращения - ц.т. тела)

17 Моделирование плавания объёмного тела. График изменения момента mZ тела от времени

18 Моделирование плавания объёмного тела. График изменения момента mY тела от времени

19 Моделирование плавания объёмного тела. График изменения момента mX тела от времени

20 На рисунке приведена взволнованная поверхность жидкости – это модификация описываемого эксперимента, когда тело совершает значительные колебания. Моделирование плавания объёмного тела.

21 Выводы: Как видно из графиков результатов, отклонения от теоретических значений невелики. На момент прекращения эксперимента тело продолжает колебаться – причиной этому может служить близость стенок бассейна. Полученные значения величин выталкивающей силы, угла дифферента, глубины погружения тела удовлетворяют аналитическим значениям для условий данного эксперимента. Таким образом, на примере данного численного эксперимента подтверждается возможность использования вычислительного комплекса FlowVision для гидростатического расчёта. А в ходе проведения описываемого численного моделирования определена и отработана методика постановки таких экспериментов.. Моделирование плавания объёмного тела.

22 Моделирование процесса установления уровня жидкости в баке с переборкой. Постановка эксперимента: Бак моделировался внутри расчётной области как подвижный фильтр, в начальном положении бак повёрнут на угол 5 о. Объём жидкости внутри бака задавался как фильтр значения VOF. Модель: Расчётная модель - Сводная поверхность, включая уравнение Навье- Стокса, k- модель турбулентности, уравнение переноса функции заполнения VOF. Размеры бака 8х5х5 м (XYZ), объём бака 200 м 3, объём жидкости 98 м 3. Число расчётных ячеек шт. Моделирование процесса установления уровня жидкости в баке с переборкой.

23

24

25

26 Моделирование обтекания плоскокилеватого днища. Для отработки методики проведения численного гидродинамического эксперимента, в среде комплекса FlowVision моделировалось обтекание плоскокилеватого днища со скоростью 8 м/с без учёта набегающего потока воздуха при фиксированных угле дифферента и глубине погружения модели. Воспроизводилось глиссирование на неполной ширине. Характерные размеры построенной модели: длина 2200мм, килеватость 30 o, ширина транца 600 мм, высота 250 Размеры расчётной области 7,2 х 1,6 х 2,4 м, глубина канала составляет 0,8 м. Масса днища в эксперименте составляла 21 кг, модель имела свободу вертикальных угловых и линейных перемещений. Моделирование обтекания плоскокилеватого днища.

27 Для расчёта использовалась математическая модель Free Surface (Свободная поверхность), с уравнениями Навье-Стокса, турбуленции и уравнение переноса функции заполнения VOF. Плоскокилеватое днище моделировалось с помощью фильтра «подвижное тело» с заданной массой, моментами инерции и начальным положением (угол дифферента, глубина погружения). Расчётная сетка начального уровня имеет размерность:151х26х33, сетка адаптировалась до 2-го уровня в зонах треугольника глиссирования и следа, а также по днищу и транцу. В результате расчётная сетка составила ячеек, из них расчётных (calculation cells) Расчёт выполнялся на ПК с процессором Pentium 1500 МГц с RAM=1Гб и Pentium 2800 МГц с RAM=2Гб (последовательно), операционная система WindowsXP SP1. Далее представлен ряд результатов данного моделирования.

28 Моделирование обтекания плоскокилеватого днища. Расчётная сетка эксперимента

29 Моделирование обтекания плоскокилеватого днища. Распределение давлений в продольном сечении, вокруг днища

30 Моделирование обтекания плоскокилеватого днища. Распределение давлений по днищу

31 Моделирование обтекания плоскокилеватого днища. Смоченная поверхность днища (глиссирование на «треугольнике»)

32 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

33 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

34 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

35 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

36 Моделирование обтекание плоскокилеватого днища.

37 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

38 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

39 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF), векторы скорости.

40 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

41 Моделирование обтекания плоскокилеватого днища. Визуализация свободной поверхности (функция VOF)

42 Сброс жидкости из бака в набегающий воздушный поток. Расчётная область представляет собой параллелепипед с цилиндром на верхней грани параллелепипеда. Цилиндром имитируется бак, параллелепипед - свободная область под баком, верхняя грань бака имеет граничное условие непротекания. Размеры: цилиндр (высота - 5м, диаметр – 3 м, V = 35,34 м 3 ), параллелепипед ( 28х20х20 м). Объём жидкости задавался с помощью фильтра VOF, т.е. задавался первоначальный цилиндрический объём жидкости в баке-цилиндре, геометрические параметры этого водяного столба – диаметр 3м, высота 3м (V=21,206 м 3 ) Далее, под действием силы тяжести (вектор силы тяжести направлен против оси Y, g=9.81 м/с2) жидкость выливается из цилиндра в набегающий поток, средняя скорость набегающего воздуха 120 м/с. Моделирование сброса жидкости из бака в набегающий воздушный поток.

43 В расчёте использовалась математическая модель k-e турбулентность, причём жидкость несжимаемая, вязкая, эта модель позволяет моделировать два вещества (в данном эксперименте - вода и воздух). В данном эксперименты используются («включены») уравнения концентрации, скорости, турбулентности. Рассчитываются два вещества – воздух (Concentration index = 0), вода (Concentration index = 1), физические параметры веществ взяты из БазыДанных программы. Расчётная сетка составила, после наложения фильтров адаптации по пути распространения жидкости, ячейки, из них рассчитываемых ячеек шт., было просчитана 3750 итераций, время задачи (время эксперимента) 0, секунды; время счёта на ПК Pentium МГц, ОП 1 Гб составило 21 час. Жидкость имела физические свойства воды (вязкость, плотность), не учитывались сжимаемость и поверхностное натяжение. Движение воздуха моделировалось заданием скорости на входе (граничное условие inlet), 70 м/с. Моделирование сброса жидкости из бака в набегающий воздушный поток.

44 Расчётная сетка

45 Моделирование сброса жидкости из бака в набегающий воздушный поток. Изолинии свободной поверхности (VOF) и линии тока

46 Моделирование сброса жидкости из бака в набегающий воздушный поток. Изолинии свободной поверхности (VOF)