
- •В. М. Шаповалов математическое моделирование процессов переноса
- •Волгоградский государственный технический университет
- •В. М. Шаповалов
- •Введение……………………………………………………………………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.6.3. Совместное течение двух несмешивающихся жидкостей в плоском канале (стратифицированное течение)
Математический анализ восполняет краткость нашей жизни и несовершенство наших чувств. Он идет одной и той же дорогой в изучении всех явлений; он объясняет их одним языком, как бы для того, чтобы подчеркнуть единство и простоту устройства вселенной.
Фурье Жане Батист (1768-1830гг).
Экспериментально обнаружено, что при течении двух несмешивающихся жидкостей в канале ограниченного сечения маловязкая жидкость обволакивает высоковязкую, отделяет ее от стенок как показано на рис. 1.9. При этом затраты энергии на транспортировку значительно снижаются. Этот эффект используется при транспортировке нефти, путем добавления воды - как маловязкой жидкости. Кроме того, существует технологический процесс производства многослойных полимерных пленок из расплавов полимеров посредством плоскощелевой формующей головки. Таким образом, рассматриваемая задача имеет широкое приложение.
Схема
течения и система координат показаны
на рис. 1.10. Две несмешивающиеся жидкости
различной вязкости и расходом движутся
п
од
действием избыточного давления вдоль
плоского канала. Ось х соответствует
направлению течения. Начало декартовой
системы координат поместим на поверхности
соприкосновения двух жидкостей. Режим
течения - ламинарный. Жидкости ньютоновские.
Течение изотермическое, установившееся.
Гидродинамическое влияние боковых
стенок не учитываем. Требуется найти
профили скоростей, расходы жидкостей
и соотношение толщины слоев.
Течение описывается системой уравнений движения в декартовой системе координат
;
;
;
Здесь индексом 1 отмечены параметры верхнего слоя, а индексом 2 – нижнего.
Граничные условия задачи. Принимаем на поверхностях канала условия прилипания:
y=1, 1=0; y=-2, 2=0. (1.4)
На границе раздела жидкостей имеет место равенство скоростей
y=0, 1=2 (1.5)
и равенство (или непрерывность) касательных напряжений
y=0,
.
(1.6)
Проинтегрируем уравнения движения
,
.
Из условия равенства касательных напряжений (1.6) следует
с1=с2.
Выполнив повторное интегрирование уравнений, имеем
,
.
Используем условие равенства скоростей на стыке жидкостей (1.5). Откуда имеем 1=c3; 2=c4. Следовательно с3=с4, При этом для распределения скоростей имеем
,
Из полученных выражений видно, что профили скорости описываются параболами.
Постоянные с1 и с3 найдем из условий прилипания жидкости к стенкам канала (1.4)
,
.
Вычитая из первого равенства второе, находим постоянную с1
.
После преобразований можем записать следующие выражения для постоянных:
,
.
Профили скоростей в каждом из слоев описываются следующими выражениями:
,
.
Общий объемный расход Q равен сумме расходов жидкостей
Q=Q1+Q2,
где
-расход верхнего слоя жидкости;
,
b
– ширина щели.
Выполнив интегрирование, найдем
,
.
Из полученных выражений видно, что профили скоростей в слоях жидкостей зависят от соотношения их толщины и вязкости.
Аналогично может быть решена задача течения любого конечного числа слоев разнородных жидкостей.
Пример.
Посредством плоскощелевой головки формуется бикомпонентная полимерная пленка. Давление в головке 20 Мпа, протяженность губок фильеры 30 мм, толщины слоев различных полимеров 0,5 мм и 0,5 мм. Ширина формуемого полотна 1 м. Вязкость первого полимера 102 Па.с, второго – 103 Па.с. Найти расходы полимеров.
Решение.
Используем расчетные формулы
,
.
Подставив численные значения, найдем
,
.
Задачи.
-
Как изменятся расходы компонентов в условиях рассмотренного примера, если давление в головке увеличить в 2 раза?
-
Какому виду течения будет соответствовать решение задачи, если в расчетных формулах положить δ1=-δ2=δ, μ1= μ2= μ?