Обтекание сферы потоком идеальной несжимаемой жидкости

Министерство образования  и науки РФ

Федеральное Государственное  Бюджетное Образовательное Учреждение Высшего Профессионального Образования

Тульский Государственный  Университет

 

 

 

 

 

 

 

 

 

 

 

Математическое моделирование

Лабораторная работа №2

Обтекание сферы потоком  идеальной несжимаемой жидкости.

 

 

 

 

 

 

 

 

 

 

 

 

 

Выполнила : студентка группы  530292               Горлова Т.И.

Принял: доцент кафедры пмии                                                  Скобельцын С.А.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тула 2013 г.

 

1. Цель работы.

Приобретение навыков математического  моделирования процесса обтекания  сферы потоком несжимаемой идеальной  жидкости и аналитического решения  этой задачи.

 

2. Задание.

1. Построить модель процесса обтекания сферы потоком жидкости.
2. Провести подробное аналитическое решение задачи.
3. Определить компоненты вектора скорости частиц жидкости.
4. Найти распределение скоростей частиц жидкости на поверхности сферы.
5. С помощью интеграла Бернулли определить давление на поверхности сферы.
6. Рассчитать зависимость давления на поверхности сферы от полярного угла q.
7. Рассчитать зависимость величин V_r и V_q от угла q при r = 1,2a; 5a; 10a.
8. Рассчитать зависимости V_r и V_q от r в диапазоне изменения r от a до 10a при любых значениях q.

 

3. Выполнение.

1. Математическая  модель процесса обтекания сферы  потоком жидкости.

Движение  несжимаемой идеальной жидкости описывается системой уравнений  гидродинамики, состоящей из уравнения  Эйлера

 

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

 

Используя формулу векторного анализа

 

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

 

следовательно,

 

где – константа, имеющая одно и то же значение во всех точках движущейся жидкости. Это равенство называется интегралом Бернулли.

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

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

 

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

Из уравнения неразрывности  следует, что потенциал скоростей  для несжимаемой жидкости должен удовлетворять всюду вне сферы уравнению Лапласа

 

Граничное условие на поверхности  сферы состоит в равенстве  нулю нормальной составляющей скорости :

 

Кроме того, потенциал  в бесконечности должен удовлетворять условию

 

Потенциал представим в виде суммы двух слагаемых

 

Где потенциал скоростей однородного потока.

  потенциал скоростей возмущенного потока.

Тогда получаем математическую модель

 

 

 

Таким образом, задача состоит в нахождении решения уравнения,  удовлетворяющего граничному условию и условию на бесконечности.

 

2. Аналитическое решение задачи.

В сферической системе координат

 

Так как рассматриваем осесимметричную  задачу, то функция  не зависит от координаты . Поэтому уравнение принимает вид

 

Будем решать уравнение методом  разделения переменных. Для этого  полагаем

 

 

Так как левая часть  последнего равенства не зависит  от , а правая от , то можем записать

 

где постоянная разделения, которую считаем неотрицательной.

Положим , где – целое число. Тогда для первого уравнения будем иметь

 
 

Положим , получим:

 

Подставим в исходное уравнение, получим:

 
 

 

Окончательно получаем решение:

 

где постоянные.

Для выполнения условия , необходимо положить

 

Решим второе уравнение:

 
 

Положим , тогда . Отсюда следует:

 
 

Подставим полученные выражения:

 

Получили  дифференциальное уравнение Лежандра.

Окончательно  получаем решение:

 

 

где постоянные, полином Лежандра порядка .

Таким образом,

 

Будем искать потенциал  в виде ряда

 

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

 

 

Учитывая, что , находим

 

Таким образом,

 

Окончательно получаем решение задачи

 

3. Определение компонент вектора скорости частиц жидкости

Определив потенциал скоростей  , можем найти компоненты вектора скорости частиц жидкости по формулам

 

Окончательно получаем:

 

 

4. Нахождение распределения скоростей частиц жидкости на поверхности сферы

На поверхности сферы, т. е. при r = a, компоненты скорости равны:

 

Тогда модуль скорости частиц жидкости на поверхности сферы равен:

 

5. Определение давления на поверхности сферы.

Для определения давления воспользуемся интегралом Бернулли:

 

где .

На поверхности сферы  будем иметь

 

Отсюда видно, что распределение  давления относительно плоскости XOY, перпендикулярной к направлению потока, симметрично. Тогда ясно, что давления, приложенные к поверхности сферы, взаимно уравновешиваются. Это, так называемый, парадокс Даламбера. Он объясняется тем, что в действительности безотрывное безвихревое обтекание сферы не имеет места. С поверхности сферы срываются вихри, которые видоизменяют как картину течения, так и распределение давления на поверхности сферы.

 

Пункты 6,7 и 8 представлены на графиках: