Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MATHEMATICAL MODELING OF DISSIPATIVE HEATING IN A CYLINDRICAL CHANNEL FOR FLUID EFFECT OF THE «SOLIDIFICATION» AT THE REALIZATION OF THE SECOND SCHEME CURRENT

Kolodezhnov V.N. 1 Kapranchikov S.S. 2 Veretennikov A.S. 1
1 Federal State Budget Educational Institution of Higher Professional Education «Voronezh State University of Engineering Technologies»
2 Federal State Military Institution of Higher Professional Education «Military Educational and Scientific Center of the Air Force
Consider the second flow pattern in a cylindrical channel for fluid rheological model which takes into account the manifestation of shear thinning fluid (STF) or the effect of «solidification». In accordance with the rheological model for such scheme currents the flow region is divided into two zones with pseudoplastic and dilatant fluid behavior . For this scheme we consider the problem of steady flow of convective heat transference, taking into account the dissipation of mechanical energy. For this there were placed at the channelthe wall temperature boundary condition of the first kind. Convective heat transference along the axis of the channel is determined by the average flow velocity. Dissipative component taken into account with the exact velocity distribution in the channel cross section An expression for the temperature distribution in the channel. Numerical experiments and analysis of the impact of model parameters on the characteristics of the process of dissipative heating. The obtained results can be used for modeling the flow of non-Newtonian fluids of such elements in the flow of the process equipment.
dissipation
non-Newtonian fluid
temperature
the effect of «solidification
shear thickening fluid
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. no. 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.

В работе [1] была предложена реологическая модель жидкости, которая при определенных значениях скоростей сдвига демонстрирует проявление эффекта «отвердевания». В [2] рассмотрена задача о течении такого рода жидкости в цилиндрическом канале. Там же было показано, что в зависимости от диапазона изменения перепада давления могут быть реализованы три схемы течения. Задача о конвективном теплопереносе для первой схемы течения, когда жидкость всюду внутри канала демонстрирует псевдопластическое поведение со степенным законом вязкости, была рассмотрена в [3]. В данной работе рассматривается подобная задача, но для второй схемы течения.

В [1] было показано, что вторая схема течения реализуется при значениях перепада давления Δp на длине L цилиндрического канала радиуса R, которое удовлетворяет следующему условию:

kolod01.wmf,

где τ2 – константа реологической модели суспензии, представляющая собой предельное значение касательного напряжения; kolod02.wmf, kolod03.wmf, kolod04.wmf – коэффициент консистенции, индекс течения и пороговое значение скорости сдвига соответственно.

При этом область течения в канале разбивается на две зоны – псевдопластического течения в центральной части и дилатантного в окрестности стенки. Граница раздела между зонами течения представляет собой цилиндрическую поверхность радиуса kolod05.wmf.

Дальнейшее увеличение Δp и, соответственно, выполнение условия

kolod06.wmf

приводит к формированию в окрестности стенки канала третьей зоны течения, заполненной «отвердевшей» жидкостью.

Рассмотрим вторую схему течения в цилиндрическом канале жидкости демонстрирующей эффект «отвердевания». Согласно реологической модели, представленной в [1], в данном случае вязкость жидкости на некотором интервале изменения скорости сдвига уменьшается, а на другом – по мере увеличения модуля скорости сдвига возрастает до некоторого максимального значения.

Решение задачи по нахождению распределения температур с учетом разбиения области течения на зоны предлагается искать в виде

kolod07.wmf,

где T(k)(r, z) – распределения температур жидкости в зонах течения, представляющие собой неизвестные функции радиальной r и продольной z координат.

Здесь и далее верхний индекс k в круглых скобках принимает значения k = 1 или k = 2 соответственно для первой и второй зон течения.

Решение поставленной задачи проводилось в безразмерном виде, с учетом соотношений, представленных в [3]. При этом гидродинамическая часть задачи решалась отдельно от тепловой. Выражения для распределения скоростей по сечению цилиндрического канала, а также средней скорости и объемного расхода, с учетом разбиения на зоны, представлены в работе [1]. В ходе решения этой части задачи было получено выражение для определения границы раздела зон течения

kolod08.wmf

kolod09.wmf kolod10.wmf

kolod11.wmf

где G, La – геометрический критерий подобия и критерий подобия Лагранжа.

Здесь и далее верхним штрихом обозначены безразмерные величины.

Предполагалось, что конвективный теплоперенос вдоль оси канала допустимо определять по средней скорости kolod12.wmf потока. Диссипативную же составляющую принимали с учетом точного распределения скорости [2] в поперечном сечении канала. С учетом этих допущений уравнение конвективного теплопереноса [4] в цилиндрическом канале с учетом диссипации в безразмерной форме может быть представлено в виде

kolod13.wmf (1)

kolod14.wmf

kolod15.wmf kolod16.wmf

kolod17.wmf kolod18.wmf

kolod19.wmf kolod20.wmf kolod21.wmf kolod22.wmf

kolod23.wmf

kolod24.wmf

kolod25.wmf

kolod26.wmf kolod27.wmf

где S1, S2 – параметры уравнения, определяемые через основные критерии подобия; T = T(r, z) – температура жидкости в канале, представляющая собой неизвестную функцию радиальной и продольной координат; τrz, kolod28.wmf – касательное напряжение и пороговое значение скорости сдвига; T**, T*, uS, μS – некоторые характерные значения температуры среды в канале, а также скорости жидкости и динамической вязкости, принимаемые в качестве масштабных; ρ, c, λ – принимаемые постоянными плотность, теплоемкость, теплопроводность среды соответственно; Pr, Ec, Re – критерии подобия Прандтля, Эккерта и Рейнольдса соответственно.

С учетом реологической модели [1] выражения для касательных напряжений записываются в виде

kolod29.wmf

kolod30.wmf

kolod31.wmf

kolod32.wmf

где n2 – параметр реологической модели.

Принимая во внимание полученные выражения для касательных напряжений, диссипативная функция W′(r′) представляется в виде

kolod33.wmf

Учитывая выражение для распределения скорости в первой зоне [2], находим

kolod34.wmf

kolod35.wmf

Прямое определение диссипативной функции для второй зоны в аналитическом виде представляется затруднительным. В этой связи представим диссипативную функцию в виде интерполирующего полинома порядка I:

kolod36.wmf

где wi – коэффициенты полинома.

При решении задачи использовались следующие температурные граничные условия:

kolod37.wmf (2)

kolod38.wmf (3)

kolod39.wmf (4)

kolod40.wmf (5)

kolod41.wmf kolod42.wmf

где T0 – температура жидкости на входе в канал; Tw – принимаемая постоянной температура стенки канала.

Решая (1) с учетом граничных условий (2)–(5), распределение температуры в канале может быть представлено в виде

kolod43.wmf (6)

kolod45.wmf

где J0 – функция Бесселя первого рода нулевого порядка [5]; εj – корни характеристического уравнения J0(ε) = 0; Cj – коэффициенты разложения.

Для краткости записи в (6) приняты следующие обозначения:

kolod46.wmf

kolod47.wmf

kolod48.wmf

kolod49.wmf

kolod50.wmf

Учитывая условие ортогональности базисных функций, коэффициенты разложения Cj в выражении (6) принимают вид

kolod51.wmf

где J1 – функция Бесселя первого рода первого порядка [5].

С полученными выражениями для распределения температуры в канале были проведены численные эксперименты по анализу влияния параметров математической модели на характеристики теплопереноса. В качестве базовых параметров модели были приняты следующие значения: n1 = 0,7; kolod52.wmf; B = 0,752; n2 = 0,233; G = 0,133; Pr = 2,476·105; Ec = 1,731·10–6; La = 50; Re = 2,07·10–3; kolod53.wmf; kolod54.wmf. Заметим, что для данного набора параметров вторая схема течения реализуется при условии изменения критерия подобия Лагранжа в диапазоне от La1 = 34,247 до La2 = 60.

Для примера на рис. 1 представлено распределение безразмерной температуры по радиальной координате в различных поперечных сечениях канала.

pic_1.wmf

Рис. 1. Распределение безразмерной температуры в различных поперечных сечениях канала для z′ = 0,05 (1); 0,1 (2); 0,2 (3); 0,4 (4); 0,6 (5); 0,8 (6); 1 (7)

Штриховой линией на этом и следующем рисунках отмечена граница раздела kolod55.wmf первой и второй зон течения. Из зависимостей, представленных на рис. 1, видно, что по мере течения жидкости по цилиндрическому каналу происходит ее прогрев за счет диссипации механической энергии. При этом наблюдается экстремум температуры, который по мере продвижения жидкости по каналу смещается из второй в первую зону течения.

На рис. 2 представлено влияние критерия подобия Эккерта на распределение температуры в выходном сечении канала.

pic_2.wmf

Рис. 2. Влияние критерия подобия Эккерта на распределение температуры в выходном сечении канала при Ec = 1·10–6 (1); 1,4·10–6 (2); 1,731·10–6 (3); 2·10–6 (4); 2,4·10–6 (5)

Из представленных данных следует, что увеличение значений критерия подобия Эккерта приводит к более интенсивному диссипативному разогреву жидкости, который сопровождается появлением пика температуры в первой зоне течения.

Работа выполнена при поддержке гранта РФФИ, проект № 12-08-00629.

Рецензенты:

Шашкин А.И., д.ф.-м.н., профессор, заведующий кафедрой «Математический и прикладной анализ», ГОУ ФБГОУ ВПО «Воронежский государственный университет», г. Воронеж;

Буховец А.Г., д.т.н., профессор кафедры «Прикладная математика и применение математических методов в экономике», ГОУ ФБГОУ ВПО «Воронежский государственный аграрный университет им. императора Петра I», г. Воронеж.

Работа поступила в редакцию 10.12.2013.