- •В. М. Шаповалов математическое моделирование процессов переноса
- •Волгоградский государственный технический университет
- •В. М. Шаповалов
- •Введение……………………………………………………………………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. Моделирование отдельных задач гидродинамики
1.1. Реологические уравнения
Искусство – это я,
наука – это мы.
В.Гюго
Для описания течений материалов необходимо знать зависимость их текучих свойств от приложенных напряжений. Указанная зависимость позволяет рассчитывать потребляемую мощность (или затраты энергии) при движении жидкости.
Наиболее распространен и давно применяется ньютоновский закон течения для так называемых ньютоновских или вязких жидкостей
,
где - касательное напряжение, - вязкость, - скорость сдвига. Часто скорость сдвига обозначают так: =.
Свойства большинства жидкостей отличаются от ньютоновского закона. Реологические свойства изучают с помощью вискозиметров.
При этом обычно организуют течение простого сдвига, сущность которого поясняется на рис. 1.1.
Реологические свойства жидкостей представляют графически с помощью «кривой течения» - зависимости касательного напряжения от скорости сдвига (см. рис. 1.2). Для описания аномальных свойств предложена и широко используется в инженерных приложениях «степенная жидкость» или жидкость Оствальда- де Виля, имеющая следующее уравнение
,
где n- индекс течения жидкости, - реологическая константа.
В зависимости от численного значения индекса течения, различают: n1 -псевдопластики (большинство полимеров, резиновые смеси и т.д.). С увеличением скорости деформации эффективная вязкость таких жидкостей снижается, n=1- ньютоновские жидкости (низкомолекулярные жидкости: вода, керосин, глицерин, спирт, но не растворы и расплавы полимеров), их вязкость постоянна в изотермических условиях течения с любой скоростью деформации, n1- дилатантные жидкости. У этих жидкостей с увеличением скорости деформации вязкость повышается. Встречаются редко. В качестве примера дилатантной среды – мокрый песок.
Кроме того, достаточно распространены вязкопластические тела Бингама, которые иногда называются средой Шведова - Бингама. К ним относятся высоконаполненные системы, например, зубная паста, краска, сметана, наполненные полимеры, резина и т.д. Реологическое уравнение для таких систем имеет вид
0, 0 ,
0, =0,
где 0-начальное напряжение сдвига, - пластическая вязкость. Течение таких систем начинается () когда напряжение сдвига превышает некоторое значение 0 При меньших напряжениях течение отсутствует (=).
Если предел текучести равен нулю (0=0), то модель Шведова-Бингама переходит в ньютоновскую.
Можно степенной закон записать в «ньютоновской» форме
,
где эф= -эффективная вязкость степенной жидкости, справедливая только в узком диапазоне скоростей деформаций. Такая запись используется для упрощенного решения задач течения аномальных жидкостей в перерабатывающем оборудовании (валках, шнеках, каландрах и др.)..
В табл. 1.1 представлены теплофизические и реологические свойства наиболее распространенных полимеров, необходимые для технологических расчетов. В скобках указаны температуры, при которых измерены параметры.
В зависимости от вязкости материала применяется различное перерабатывающее оборудование. Так, роторные смесители – для материалов с вязкостью 7-1000 Па.с, шнековые – для 2-1000 Па.с, валковые – для 10-10000 Па.с.(имеется ввиду эффективная вязкость).
В условиях переработки большинство полимеров ведет себя как аномально-вязкие (псевдопластические) жидкости. Наиболее простой зависимостью, связывающей напряжения сдвига τ с градиентом скорости деформации γ, является степенной закон или закон Освальда-де Виля
τ =μоγn,
где μо- коэффициент консистентности, Па.сn, n – индекс течения (зависят от типа полимера и температуры).
Таблица 1.1
Ориентировочные технологические свойства полимеров
Показатель |
Резина |
ПЭНП |
ПЭВП |
ПП |
ПС |
ПКА |
n |
0,21 (60оС) |
0,4 (190оС) |
0,47(230оС) |
0,27 (200оС) |
0,4 (260оС) |
0,8 (288оС) |
o,Паcn |
5x104 |
1,4x104 (190оС) |
6,6x103 (230оС) |
1,9x104 |
4000 (260оС) |
300 (288оС) |
,кг/м3 |
1200 |
920 |
951 |
900 |
1100 |
970 (288оС) |
ах107, м2/с |
1 |
1,2 |
1,5 |
1 |
1,1 |
1(220оС) |
, Вт/мК |
0,14 |
0,3 |
0,4 |
0,24 |
0,15 |
0,24 (240oC) |
С, кДж/кгК |
2 |
2,8-3 |
2,5-3 |
3,2 (170оС) |
1,5 |
2,6 (240оС) |
Е, кДж/моль |
34-70 (40оС-140оС) |
50 |
50 |
43 |
94 |
44 |
Температура (оС): в зоне дозирования |
|
150-180 |
170-190 |
170-250 |
190-240 |
300 |
В головке |
|
160-210 |
170-210 |
190-260 |
220-260 |
280 |
Средний градиент скорости при экструзии (с-1) |
|
50-240 |
50-230 |
55-110 |
50-200 |
50-110 |
При инженерных расчетах широко используют эффективную вязкость μэф, которая определяется так:
μэф=μoγсn-1,
где γс – среднее значение градиента скорости сдвига (см. табл.1.1). При этом степенной закон превращается в ньютоновский
τ=μэфγ.
Это открывает возможность использовать ньютоновское приближение для расчетов течения материала в перерабатывающем оборудовании. Следует отметить, что вопрос о среднем значении градиента скорости сдвига γс достаточно сложен и μэф сильно отличается от параметра μо.
Если имеет место неизотермическое течение, то необходимо учитывать зависимость текучих свойств жидкости от температуры. Широко используется закон Аррениуса
,
где Е – энергия активации вязкого течения (см. табл. 1.1), R=8,3 Дж/моль – универсальная газовая постоянная, μо- вязкость (или показатель консистенции) при температуре То. Используются абсолютные температуры в градусах Кельвина.
Пример.
Найти эффективную вязкость резины при ее переработке на экструдере, если реологические константы степенной модели μо=105 Пасn, n=0,25. Средний градиент скорости γ=100 с-1.
Решение.
Используем формулу
μэф=μоγn=105*1000,25=3,16*105 Па.с.
Задачи.
-
При исследовании на вискозиметре вязкопластической жидкости получены следующие данные: γ1=50 с-1, τ1=103 Па, γ2=100 с-1, τ2=1,5*103 Па. Найти начальное напряжение сдвига.
-
Чему равно касательное напряжение при течении степенной жидкости если: μо=104 Пасn, n=0,5, γ=200 с-1?