Файл: Численное моделирование ламинарного потока в канале.docx

Санкт-Петербургский политехнический университет Петра ВеликогоИнститут машиностроения, материалов и транспортаВысшая школа машиностроенияОтчётДисциплина: «Связанные задачи гидро и аэроупругости»Тема: «Численное моделирование ламинарного потока в канале»Студенты гр. 3331501/70601 Осман А. А.Преподаватель Гатаулин Я. А. « »________2022г.Санкт-Петербург2022

Оглавление

Цель работы

Провести численный расчет во Fluent стационарного ламинарного течения вязкой несжимаемой жидкости в трубки кольцеобразного сечения для разных значений числа Рейнольдса: 100, 200, 400.

Задача

1) Построить геометрию и расчетную сетку;2) Выполнить расчет потока для заданных чисел Рейнольдса;3) Сравнить рассчитанный коэффициент сопротивления с теоретическими значениями.

Исходные данные

П араметры геометрии трубы:

• 
  Поперечное сечение: кольцо;
  
• 
  Внешний диаметр
  
• 
  Внутренний диаметр
  
  
  Рисунок 1 – сечение трубы
  
• 
  Длина канала:
  

Параметры жидкости:

• 
  Плотность
  
• 
  Коэффициент динамической вязкости
  
• 
  Числа Рейнольдса: .
  

Граничные условия:

1. 
   Скорость на входе: м/с;
   
2. 
   Давление на выходе: ;
   
3. 
   Условие прилипания на стенках: ;
   

---

Рассчитываемые параметры

В рассматриваемой задаче течение определяется числом Рейнольдса, определяемым по формуле:где плотность жидкости, скорость потока, м/с; диаметр трубы, м; коэффициент динамической вязкости жидкости, .Параметры течения, соответствующие заданным числам Рейнольдса представлены в таблице 1.Таблица 1 – Параметры потока

Re	100	200	400
, Па с	1	1	1
, м	1	1	1
d, м	0.5	0.5	0.5
, кг/м3
V, м/с	1	2	4

Математическая модель

В качестве математической модели выбрана модель несжимаемой жидкости, течение – ламинарное, задача – стационарная.Уравнение несжимаемости:Уравнения Навье-Стокса:

---

Создание геометрии

Для создания геометрии трубы использовалась программа DesignModeler. Геометрия была получена с получена с помощью эскиза и операции вытягивания.Рисунок 2 – геометрия трубы для расчёта

Создание сетки

Расчётная сетка создана в программе ANSYS Meshing. Для построения сетки использован метод протяжки (Sweep), на торце задано разбиение линий, как показано на рисунке 3 (а), также использовалась опция Inflation для создания сгущения сетки к стенкам трубы (рисунок 3 (б)).Рисунок 3 (а, б)По длине трубы в направлении оси Z задано количество разбиений и сгущение сетки к началу трубы (рисунок 4).Рисунок 4Количество элементов в расчётной сетке составило 440440 шт.Входу в канал, выходу и стенкам даны названия.

Расчёт в Fluent

В программе Fluent настроена точность, заданы параметры течения, материала, а также граничные условия.Рисунок 5 – граничные условия.Проведён расчёт для 100 итераций, сходимость решения для Re = 400 представлена на рисунке 6.Рисунок 6 – сходимость решения.

Анализ результатов решения

Векторное поле скоростей.

На рисунках ниже представлены векторные поля скоростей исследуемых потоков в центральном сечении (Х=0).Рисунок 7 - векторное поле скоростей для Re = 100Рисунок 8 - векторное поле скоростей для Re = 200

---

Рисунок 9 - векторное поле скоростей для Re = 400Максимальные скорости потока для разных чисел Рейнольдса представлены в таблице 2.Таблица 2 – максимальные расчётные скорости потока

Re	100	200	400
Vmax, м/с	1.501	2.997	5.893

Профили скорости

На рисунках ниже представлены профили скоростей для разных чисел Рейнольдса для половины сечения тубы.Рисунок 10 - профили скоростей для Re = 100 Рисунок 11 - профили скоростей для Re = 200Рисунок 11 - профили скоростей для Re = 400Из графиков видно, что на вход задаётся поток скорости равный расчётному значению, следуя далее по каналу он приобретает М-образную форму, а затем приобретает вид параболы Пуазейля.

Графики скоростей, давления и касательных напряжений

На рисунках ниже представлен графики распределения максимальных скоростей, давления и касательных напряжений по всей длине трубы вдоль проекции серединной цилиндрической поверхности трубы, изображённой на рисунке 12 (далее упоминаемой как средняя линия течения).Рисунок 12 – линия сбора исследуемых показателей расчётаРисунок 13 – распределение скоростей потоков вдоль средней линии течения.Рисунок 14 – распределение давления вдоль средней линии течения.Рисунок 15 - распределение касательных напряжений на внутренней стенке трубы вдоль средней линии течения.

---

Рисунок 15 - распределение касательных напряжений на внутренней стенке трубы вдоль средней линии течения.Из полученных графиков скорости можно определить длину начального участка, как расстояние от входа до места, где скорость потока составляет 98% от максимальной. Для Re=100 длина начального участка состаляет примерно 0.78 м. Для Re=200 – примерно 1.2 м. Для Re=400 течение при длине трубы 2 м не успело стать установившимся. Удлинение трубы для более подробного изучения потока с Re=400 осложнено ограничениями студентческой версии ANSYS по количеству элементов.Из графиков распределения касательных напряжений на внешней и внутренней стенках трубы можно заметить, что на входе напряжения на внешней стенке трубы значительно меньше напряжений на внутренней стенке, однако на установившемся участке значения напряжений становятся близкими. Данная ассиметрия на входе наиболее ярко заметна на поле скоростей:Рисунок 16 – поле скоростей для Re=100

Сравнение аналитического решения с полученным в программе ANSYS Fluent

Аналитические решение для скорости в установившемся течении , а также для средней скорости вычисляются по формуле:где — разница давлений (рассчитывается для участка с установившимся течением); — длина исследуемого участка с установившемся течением (для длины начального участка взято значение рассчитанное в ANSYS за неимение точной расчётной формулы); Па/м. Па/м. Па/м.Тогда, учитывая параболическую форму профиля скорости на установившемся течении, аналитические значения скорости для потоков: Re=100:

---