
- •В. М. Шаповалов математическое моделирование процессов переноса
- •Волгоградский государственный технический университет
- •В. М. Шаповалов
- •Введение……………………………………………………………………6
- •Глава 1. Моделирование отдельных задач гидродинамики…….… 12
- •1.6.2. Течение вязкой жидкости в плоской щели………………………31
- •Глава 2. Отдельные задачи теории теплопроводности………….……79
- •2.7. Теория пленочной конденсации Нуссельта……..…………………106
- •Приложение 6. Фрагмент выпускной бакалаврской работы Шахминой е.В. На тему «Построение математической модели течения материала в резиносмесителе «Бенбери»»….………..……….154
- •Введение
- •Глава 1. Моделирование отдельных задач гидродинамики
- •1.1. Реологические уравнения
- •Истечение жидкости из отверстия в дне сосуда
- •1.3. Истечение при переменном уровне
- •1.4. Истечение высоковязкой жидкости из сосуда
- •1.5. Течение жидкости по вертикальной поверхности
- •Проинтегрируем уравнение движения
- •1.6. Установившееся движение вязкой несжимаемой жидкости по цилиндрическим и призматическим трубам
- •1.6.1. Течение в трубе эллиптического сечения
- •Лука Пачоли
- •1.6.2. Течение вязкой жидкости в плоской щели
- •1.6.3. Совместное течение двух несмешивающихся жидкостей в плоском канале (стратифицированное течение)
- •1.6.4. Аксиальное течение вязкой жидкости в кольцевом канале
- •1.6.5. Напорное течение вязкой жидкости в прямоугольном канале
- •1.6.6. Течение степенной жидкости в плоской щели
- •1.6.7. Течение среды Бингама в круглой трубе
- •1.7. Течение в кольцевом зазоре при поступательном движении внутреннего цилиндра
- •1.8. Сдвиговое течение вязкой жидкости в клинообразном зазоре
- •1.9. Течение вязкой жидкости в коаксиальном зазоре при вращении одного из цилиндров
- •1.10. Фильтрация через неподвижные пористые слои. Закон Дарси
- •1.10.1. Фильтрация через плоскую пористую стенку
- •1.10.2. Фильтрация через пористую цилиндрическую стенку
- •Глава 2. Отдельные задачи теории теплопроводности
- •2.1. Диссипативный саморазогрев жидкости в условиях простого сдвига
- •2.2. Диссипативный саморазогрев при напорном течении
- •2.3. Теплопроводность охлаждающего ребра
- •2.4. Нестационарная теплопроводность пластины
- •2.5. Нестационарная теплопроводность неограниченного цилиндра
- •Нестационарная теплопроводность шара при граничных условиях 3-го рода
- •2.7. Теория пленочной конденсации Нуссельта
- •2.9. Нестационарная массопроводность плоской стенки
- •Литература
- •Приложения
- •Соответственно, на валу второго ротора
- •Потребляемая мощность первым ротором . Вторым
- •Аналогично
- •Литература
- •400131 Волгоград, просп. Им. В. И. Ленина, 28
- •400131 Волгоград, ул. Советская, 35.
1.4. Истечение высоковязкой жидкости из сосуда
Число лежит в основе вселенной.
Пифагор.
Сила рождается из ограничений, а свобода ее губит.
Леонардо да Винчи
Полученные выше расчетные выражения правомерны только в случае течения маловязких жидкостей в турбулентном режиме (т.е. при больших числах Рейнольдса). Формулой Торричелли можно пользоваться лишь при больших числах Рейнольдса. В случае значительной вязкости жидкости коэффициент расхода непостоянен и зависит от критерия Рейнольдса.
Выше было показано, что объемный расход жидкости при истечении из открытого сосуда определяется формулой
.
При
этом для турбулентного режима мы полагали
постоянное значение коэффициента
расхода =const
0,63. В широком интервале чисел Рейнольдса
экспериментально установлена следующая
зависимость (Re),
которая представлена на рис. 1.5. Для
Re50,
что соответствует ламинарному течению
высоковязкой жидкости, зависимость
носит линейный характер и аппроксимируется
формулой =Re/25,2.
Соответственно, область непосредственного
использования формулы Торричелли (с
=const,
горизонтальный участок линии рис. 1.5)
ограничена значениями Re>50.
Рассмотрим истечение высоковязкой жидкости конечного объема из сосуда с отверстием в дне. Условие равенства элементарных объемов (см. предшествующую задачу) записывается
,
где
,
- число Рейнольдса для условий слива.
Отметим, что в качестве характерной
скорости в числе Рейнольдса используется
скорость слива по Торричелли.
Таким образом, для скорости истечения в ламинарном режиме можно записать
.
Соответственно, дифференциальное уравнение сохранения объема жидкости, примет вид
.
Разделим переменные и проинтегрируем, с учетом начального условия t=0, h=h1
.
В результате интегрирования находим время слива
.
Согласно полученному выражению время слива пропорционально вязкости жидкости (). Для сокращения времени слива достаточно, например, жидкость подогреть. Необходимо отметить, что при полном сливе h2=0, . Полученное решение имеет ограничение, поскольку при малых уровнях жидкости происходит прорыв мениска и картина течения существенно меняется.
Если над поверхностью жидкости избыточное давление (сосуд закрытый PPатм) то время слива сокращается, и расчетная формула приобретает вид
.
Пределы применяемости полученных формул
определяются неравенством
.
Пример.
Цилиндрический сосуд диаметром D=0,1 м и высотой h=0,1 м заполнен глицерином (вязкость μ=1 Па.с, плотность ρ=1260 кг/м3). За какое время вытечет половина объема жидкости через отверстие в дне сосуда, если диаметр отверстия d=5 мм?
Решение.
Используем формулу
.
Подставив численные значения, получим
.
Задачи.
-
Как повлияет на продолжительность слива, в условиях предыдущей задачи, подогрев глицерина?
-
Определить кинематическую вязкость жидкости (ν=μ/ρ), если время половинного опорожнения сосуда диаметром 2 см. с отверстием диаметром 2 мм составило 300 секунд.