В работе [1] была предложена реологическая модель жидкости, которая при определенных значениях скоростей сдвига демонстрирует проявление эффекта “отвердевания”. В [2] рассмотрена задача о течении такого рода жидкости в цилиндрическом канале. Там же было показано, что в зависимости от диапазона изменения перепада давления могут быть реализованы три схемы течения. Задача о конвективном теплопереносе для первой схемы течения, когда жидкость всюду внутри канала демонстрирует псевдопластическое поведение со степенным законом вязкости, была рассмотрена в [3]. В данной работе рассматривается подобная задача, но для второй схемы течения.
В [1] было показано, что вторая схема течения реализуется при значениях перепада давления 
на длине 
цилиндрического канала радиуса 
, которое удовлетворяет следующему условию
;
где 
- константа реологической модели суспензии, представляющая собой предельное значение касательного напряжения; 
, 
, 
- коэффициент консистенции, индекс течения и пороговое значение скорости сдвига, соответственно.
При этом область течения в канале разбивается на две зоны – псевдопластического течения в центральной части и дилатантного в окрестности стенки. Граница раздела между зонами течения представляет собой цилиндрическую поверхность радиуса 
.
Дальнейшее увеличение 
и, соответственно, выполнение условия
,
приводит к формированию в окрестности стенки канала третьей зоны течения, заполненной “отвердевшей” жидкостью.
Рассмотрим вторую схему течения в цилиндрическом канале жидкости демонстрирующей эффект “отвердевания”. Согласно реологической модели, представленной в [1], в данном случае вязкость жидкости на некотором интервале изменения скорости сдвига уменьшается, а на другом - по мере увеличения модуля скорости сдвига возрастает до некоторого максимального значения.
Решение задачи по нахождению распределения температур с учетом разбиения области течения на зоны предлагается искать в виде:
где T(k), (r,z) - распределения температур жидкости в зонах течения, представляющие собой неизвестные функции радиальной 
и продольной 
координат.
Здесь и далее верхний индекс 
в круглых скобках принимает значения 
или 
, соответственно, для первой и второй зон течения.
Решение поставленной задачи проводилось в безразмерном виде, с учетом соотношений, представленных в [3]. При этом гидродинамическая часть задачи решалась отдельно от тепловой. Выражения для распределения скоростей по сечению цилиндрического канала, а также средней скорости и объемного расхода, с учетом разбиения на зоны, представлены в работе [1]. В ходе решения этой части задачи было получено выражение для определения границы раздела зон течения
;
; 
; 
;
где 
, 
- геометрический критерий подобия и критерий подобия Лагранжа.
Здесь и далее верхним штрихом обозначены безразмерные величины.
Предполагалось, что конвективный теплоперенос вдоль оси канала допустимо определять по средней скорости 
потока. Диссипативную же составляющую принимали с учетом точного распределения скорости [2] в поперечном сечении канала. С учетом этих допущений уравнение конвективного теплопереноса [4] в цилиндрическом канале с учетом диссипации в безразмерной форме может быть представлено в виде
; (1)
;
; 
; 
;
; 
; 
; 
; 
;
; 
;
; 
; 
;
где 
, 
- параметры уравнения, определяемые через основные критерии подобия; 
- температура жидкости в канале, представляющая собой неизвестную функцию радиальной и продольной координат; 
, 
- касательное напряжение и пороговое значение скорости сдвига; 
, 
, 
, 
- некоторые характерные значения температуры среды в канале, а также скорости жидкости и динамической вязкости, принимаемые в качестве масштабных; 
, 
, 
- принимаемые постоянными плотность, теплоемкость, теплопроводность среды, соответственно; 
, 
, 
- критерии подобия Прандтля, Эккерта и Рейнольдса, соответственно.
С учетом реологической модели [1] выражения для касательных напряжений записываются в виде:
;
;
,
где 
- параметр реологической модели.
Принимая во внимание полученные выражения для касательных напряжений, диссипативная функция 
представляется в виде:
Учитывая выражение для распределения скорости в первой зоне [2] находим
;
.
Прямое определение диссипативной функции для второй зоны в аналитическом виде представляется затруднительным. В этой связи представим диссипативную функцию в виде интерполирующего полинома порядка 
,
где 
- коэффициенты полинома.
При решении задачи использовались следующие температурные граничные условия
(2)
(3)
(4)
(5)
; 
,
где 
- температура жидкости на входе в канал; 
- принимаемая постоянной температура стенки канала.
Решая (1) с учетом граничных условий (2) – (5), распределение температуры в канале может быть представлено в виде
; (6)
,
где 
- функция Бесселя первого рода нулевого порядка [5]; 
- корни характеристического уравнения 
; 
- коэффициенты разложения.
Для краткости записи в (6) приняты следующие обозначения
;
;
;
.
Учитывая условие ортогональности базисных функций, коэффициенты разложения 
в выражении (6) принимают вид
,
где 
- функция Бесселя первого рода первого порядка [5].
С полученными выражениями для распределения температуры в канале были проведены численные эксперименты по анализу влияния параметров математической модели на характеристики теплопереноса. В качестве базовых параметров модели были приняты следующие значения: 
; 
; 
; 
; 
; 
; 
; 
; 
; 
; 
. Заметим, что для данного набора параметров вторая схема течения реализуется при условии изменения критерия подобия Лагранжа в диапазоне от 
до 
.
Для примера, на рис.1 представлено распределение безразмерной температуры по радиальной координате в различных поперечных сечениях канала.
Рис. 1. Распределение безразмерной температуры в различных поперечных сечениях канала для 
0,05 (1); 0,1 (2); 0,2 (3); 0,4 (4); 0,6 (5); 0,8 (6); 1 (7).
Штриховой линией на этом и следующем рисунках отмечена граница раздела 
первой и второй зон течения. Из зависимостей, представленных на рис. 1, видно, что по мере течения жидкости по цилиндрическому каналу происходит ее прогрев за счет диссипации механической энергии. При этом наблюдается экстремум температуры, который по мере продвижения жидкости по каналу смещается из второй в первую зону течения.
На рис. 2 представлено влияние критерия подобия Эккерта на распределение температуры в выходном сечении канала.
Рис. 2. Влияние критерия подобия Эккерта на распределение температуры в выходном сечении канала, при 
(1); 
(2); 
(3); 
(4); 
(5).
Из представленных данных следует, что увеличение значений критерия подобия Эккерта приводит к более интенсивному диссипативному разогреву жидкости, который сопровождается появлением пика температуры в первой зоне течения.
Работа выполнена при поддержке гранта РФФИ, проект № 12-08-00629.
Список литературы
- 1.^Колодежнов В. Н. Математическое моделирование реологического поведения нелинейно-вязких жидкостей, которые демонстрируют проявление эффекта “отвердевания” // Вестник Воронежского государственного университета инженерных технологий. – 2012. – № 4. – С. 35-38.
- 2.^Колодежнов В. Н. Математическое моделирование течения в цилиндрическом канале жидкости, которая демонстрирует проявление эффекта “отвердевания” // Вестник Воронежского государственного технического университета. – 2013. – Т.9. № 2. – С. 118 – 122.
- 3.^Колодежнов В. Н., Капранчиков С. С., Веретенников А. С. Математическое моделирование диссипативного разогрева в цилиндрическом канале для жидкости, демонстрирующей проявление эффекта «отвердевания», при реализации первой схемы течения // Фундаментальные исследования. – 2013. №10. – с. 21 – 24.
- 4.^Фройштетер Г.Б., Данилевич С.Ю., Радионова Н.В. Течение и теплообмен неньютоновских жидкостей в трубах. – Киев.: Наукова думка, 1990. – 216 с.
- 5.^Янке Е., Эмде Ф., Лёш Ф. Специальные функции. – М.: Наука. 1977. – 344 С.
References
- 1.^Kolodezhnov V.N. Mathematical modeling of the rheological behavior of nonlinear-viscous liquids, which show a manifestation of the effect “solidification” // Proceedings of the Voronezh State University of Engineering Technologies. 2012. Vol. 4. p. 35 – 38.
- 2.^Kolodezhnov V.N. Mathematical modeling of fluid in a cylindrical channel, which shows a manifestation of the effect “solidification” // Proceedings of the Voronezh state technical University. 2013. Vol. 9. № 2. p. 118 – 122.
- 3.^Kolodezhnov V.N., Kapranchikov S.S., Veretennikov A.S. Mathematical modeling of dissipative heating in a cylindrical fluid passages that demonstrates the manifestation of the “solidification” AT the realization of the first scheme current // Technical Sciences, Vol .10, 2013. – p. 21-24.
- 4.^Frojshteter G.B., Danilevich S.Ju., Radionova N.V. Flow and heat transfer of non-Newtonian fluids in pipes. Kiev, Naukova dumka, 1990. 216 p.
- 5.^Janke E., Jemde F., Ljosh F. Special functions. Moscow, Nauka, 1977. 344 p.
Рецензенты:
Шашкин А. И., д.ф.-м.н., профессор, заведующий кафедрой «Математического и прикладного анализа», ГОУ ФБГОУ ВПО «Воронежский государственный университет», г. Воронеж;
Буховец А. Г., д.т.н., профессор кафедры «Прикладной математики и применение математических методов в экономике», ГОУ ФБГОУ ВПО «Воронежский государственный аграрный университет им. императора Петра I» г. Воронеж.
Библиографическая ссылка
Гаврилюк М.Н, Семериков И.С СИНТЕЗ И ИЗУЧЕНИЕ СВОЙСТВ ПОРТЛАНДЦЕМЕНТА ИЗ ИЗВЕСТНЯКА, ФЕЛЬЗИТА И МИНЕРАЛИЗАТОРА // Фундаментальные исследования. 2013. № 9-4. ;URL: https://fundamental-research.ru/ru/article/view?id=33337 (дата обращения: 30.06.2026).