Файл: В.Д. Богатырев Моделирование течения жидкости при действии силы тяжести.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.05.2024
Просмотров: 27
Скачиваний: 0
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
Государственное образовательное учреждение высшего профессионального образования
«КУЗБАССКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Кафедра теоретической и геотехнической механики
МОДЕЛИРОВАНИЕ ТЕЧЕНИЯ ЖИДКОСТИ ПРИ ДЕЙСТВИИ СИЛЫ ТЯЖЕСТИ
Методические указания по выполнению лабораторной работы по курсу «Подземная гидрогазодинамика» для студентов специальности 070600 «Физические процессы горного производства»
Составители В.Д. Богатырев Е.В. Гурский
Утверждены на заседании кафедры Протокол № 9 от 17.03.03
Рекомендованы к печати учебнометодической комиссией по специальности 070600 Протокол № 9 от 17.03.03
Электронная копия находится в библиотеке главного корпуса ГУ КузГТУ
Кемерово 2003
1
ВВЕДЕНИЕ
Физические процессы, изучаемые горной наукой, отличаются большой сложностью и описываются системами дифференциальных уравнений, а также начальными и граничными условиями с большим числом переменных. Найти решения таких задач довольно трудно. Поэтому на практике большое значение приобретает применение методов физического моделирования и опытных исследований, с помощью которых можно получить численные значения искомых переменных, а затем подобрать уравнения, описывающие результаты опытов.
Для изучения влияния на процесс какой-либо одной величины остальные следует сохранять неизменными, однако из-за большого числа переменных добиться этого часто очень сложно. Кроме того, при этом нужно быть уверенным, что результаты, получаемые с помощью какойлибо установки (модели), можно перенести и на другие аналогичные процессы (образец). Эти проблемы позволяет решить теория подобия.
С помощью теории подобия ряд размерных физических величин можно объединить в безразмерные комплексы. При этом число этих комплексов меньше числа исходных величин. Эти комплексы можно рассматривать как новые переменные. В результате введения в уравнения, описывающие исследуемый процесс, безразмерных комплексов число величин под знаком искомой функции значительно сокращается, что упрощает изучение физического процесса. Таким образом, теория
подобия и, в частности, данная лабораторная работа являются теоретической базой эксперимента и практически необходимы при теоретических исследованиях.
2
1. ЦЕЛЬ РАБОТЫ
Дать представление о правильном планировании эксперимента или лабораторной работы при рассмотрении стационарного течения несжимаемой вязкой жидкости для следующих случаев: при наличии или отсутствии влияния силы тяжести, вязкостном или инерционном течениях.
2.ТЕОРЕТИЧЕСКИЕ ПОЛОЖЕНИЯ
2.1.Условные обозначения, принятые в работе
V – величина вектора скорости;
Vx, Vy, Vz – проекции вектора v на координатные оси x, y, z (причем ось x направлена вертикально);
P – давление в некоторой точке потока; ρ – плотность жидкости;
ν = µ/ρ – коэффициент кинематической вязкости жидкости; µ – коэффициент вязкости жидкости;
g – ускорение свободного падения;– оператор Лапласа;
"м" – индекс, относящийся к параметрам физической модели; "н" – индекс, относящийся к параметрам натурного образца. Величины ρ, ν, g рассматриваются как постоянные во всем поле.
2.2. Течение при наличии влияния силы тяжести
Уравнения движения:
Vx |
∂Vx |
|
|
∂Vy |
|
|
|
∂Vz |
|
|
|
|
|
1 ∂P |
|
|
2 |
|||||||||||||||
|
|
|
+Vy |
|
|
|
|
+Vz |
|
|
|
= |
|
|
|
∂x + g +ν |
Vx ; |
|||||||||||||||
∂x |
∂y |
∂z |
|
|
ρ |
|||||||||||||||||||||||||||
Vx |
∂Vy |
|
∂Vy |
|
∂Vy |
|
1 ∂P |
+ν 2Vy ; |
|
|
||||||||||||||||||||||
|
+Vy |
|
+Vz ∂z |
= |
|
|
|
∂y |
|
|
|
|||||||||||||||||||||
∂x |
∂y |
|
ρ |
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
∂Vz |
|
|
|
|
∂Vz |
|
|
|
|
∂Vz |
|
|
|
1 ∂P |
|
2 |
|
|||||||||||
|
Vx |
|
+Vy |
|
|
+Vz |
|
|
= |
|
|
|
∂z |
+ν |
Vz . |
|||||||||||||||||
|
∂x |
∂y |
∂z |
|
ρ |
|
||||||||||||||||||||||||||
Уравнение неразрывности: |
|
∂Vy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂V |
x |
+ |
|
+ |
∂V |
z |
= 0. |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
∂y |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
(1)
(2)
3
Уравнения движения (1) и неразрывности (2) могут быть записаны и в векторной форме:
уравнение движения
|
|
|
1 |
|
|
|
|
|
|
|
(V gradV ) = − |
gradP + g +ν 2V ; |
(3) |
||||||||
ρ |
|
|||||||||
уравнение неразрывности |
|
|
|
|
|
|
|
|||
|
|
divV |
|
= 0. |
(4) |
|||||
|
|
|
Краевые условия в данной задаче могут содержать только скорость, которая задается на границе в виде уравнения, выражающего скорость как функцию координат точек границы. Если рассматривается течение жидкости внутри какого-либо аппарата, то границами потока следует считать его входное сечение, где имеют место управляемые краевые условия и поверхности твердых стенок, где существуют неуправляемые граничные условия. В случае неподвижных стенок V=0. Поскольку это условие автоматически соблюдается и в модели и в натурном образце, то задание краевых условий сводится лишь к заданию скорости на входном сечении изучаемого объекта. Приведем рассматриваемые уравнения и граничные условия к безразмерному виду. Для этого выберем масштабы приведения, в качестве которых целесообразно принимать величины, входящие в граничные условия. Для линейных величин выберем какой-либо характерный размер l0, например гидравлический радиус для скорости V0 . После масштабных преобразований
( x = l0X; y = l0Y; z = l0Z; v = V0V; ν =ν0N; p = p0P; g = g0G; ρ = ρ0R )
и подстановки новых значений в уравнения (1) и (2) получим соотношения между величинами, при которых уравнения (1) и (2), а также граничные условия будут приведены к безразмерному виду. Будем иметь
V 2 |
|
P |
|
ν |
V |
|
|||
|
0 |
= |
0 |
= g0 = |
|
|
0 0 |
. |
(5) |
l0 |
ρ0l0 |
|
|
|
|||||
|
|
|
|
|
l02 |
|
Из условий (5) можно получить три уравнения связи, которые содержат шесть масштабов. Следовательно, любые три из них, имеющие независимые размерности, можно выбрать, например l0, ρ0, g0, а остальные три определить из уравнений связи (5):
|
|
|
|
|
|
|
V0 = |
|
; |
P = ρgl0 ; ν0 = gl03 . |
(6) |
||
gl0 |
Используя соотношения (6), можно записать выражения для безразмерных величин, например безразмерный коэффициент вязкости
N =ν gl03 , |
(7) |
4
а ускорение свободного падения G = 1.
Из выражений (5) и (7) следует, что осуществление требования
Nм = Nн при условии, что жидкость в модели и образце одна и та же, возможно только тогда, когда размеры образца и модели одинаковы,
т.е. l0м = l0н, только в этом случае будут соблюдаться равенства по скоростям граничных условий в модели и образце.
2.3. Течение при отсутствии влияния силы тяжести
Этот процесс характеризуется отсутствием свободных поверхностей, т.е. имеет место напорное течение. В данном случае влиянием сил тяжести можно пренебречь по сравнению с влиянием остальных сил, действующих на элементарную частицу жидкости. На основании анализа дифференциальных уравнений (1) и (2) и граничных условий при отсутствии в уравнениях параметра g получим систему уравнений связи
V 2 |
|
P |
|
ν |
V |
|
|||
|
0 |
= |
0 |
= |
|
|
0 0 |
. |
(8) |
l0 |
ρ0l0 |
|
|
|
|||||
|
|
|
|
|
l02 |
|
Используя уравнения (8), можно показать, что для достижения подобия нужно на границах модели обеспечить распределение скорости V, тождественно одинаковое с тем, которое имеется в образце. Обычно результаты опытов при рассмотрении данного вида течения жидкости обрабатывают для установления эмпирических законов движения, переходя к усредненным по сечению значениям скоростей и давлений, и заменяют их значения критериями подобия Eu и Re для усредненных величин, т.е. получают зависимости вида
Eu=ψ( Re).
2.4. Вязкостный режим течения
Этот процесс характеризуется исчезающе малым влиянием инерционных сил по сравнению с влиянием сил трения и давления. В данном случае на основании анализа уравнений:
движения –
− gradP + µ 2V = 0 ;
неразрывности –
divV = 0
и граничных условий –
Vг = ϕ(lг);
5
получаем одно уравнение связи
P0 = µ0V0 , l0 l02
используя которое можно установить подобие двух потоков.
Для осуществления подобия необходимо на границах модели обеспечить распределение скорости V, тождественно одинаковое с тем, которое имеется в образце. Отсюда следует, что подобие может быть осуществлено независимо от размеров модели.
2.5. Инерционный режим течения
В этом случае силы трения настолько малы, что их влиянием по сравнению с влиянием сил инерции и давления на режим течения можно пренебречь. Уравнения значительно упрощаются:
движения –
VgradV = ρ1 gradP ;
неразрывности –
divV = 0
и граничных условий –
Vг = ϕ(lг).
После приведения уравнений к безразмерному виду получаем одно уравнение связи
V02 = P0 .
l 0
Следовательно, условия подобия будут те же, что и при вязкостном течении. Условие подобия может быть достигнуто независимо от размеров физической модели, которая, разумеется, должна быть геометрически подобна образцу.
Таким образом, при физическом моделировании гидродинамических процессов, механизм которых для нас недостаточно ясен, необходимо соблюдать динамическое подобие, как показано в примере.
6
3. ПРИМЕР. Моделирование глубины центральной воронки в сосуде с мешалкой
Найти глубину центральной воронки при стационарном движении жидкости в большом резервуаре без отражательных перегородок (рис.1), заполненном маслом, как функцию скорости мешалки. Сделать это предполагается путем проведения модельных опытов в геометрически подобном сосуде меньших размеров. Поэтому определите условия, при которых следует осуществлять модельные испытания, чтобы обеспечить правильность предсказания результатов.
N1
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
ν =ν1 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ν = ν1 / 8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Н1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Н |
|
Т |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
D1 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
||
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
Рис.1
Решение. Структура потоков в аппарате настолько сложна, что не допускает выполнения точных расчетов. Следовательно, нужно использовать методы теории подобия и анализа размерностей. Форма центральной воронки для любого сосуда с мешалкой будет одна и та же при одинаковых безразмерных уравнениях и безразмерных граничных условиях, описывающих течение. Дифференциальные уравнения (движения и неразрывности) мы уже рассмотрели. Граничные условия следующие: