Научная статья на тему 'МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛООБМЕННОГО АППАРАТА С УЧЕТОМ СИЛЬНОЙ ЗАВИСИМОСТИ ВЯЗКОСТИ НЕФТИ ОТ ТЕМПЕРАТУРЫ'

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛООБМЕННОГО АППАРАТА С УЧЕТОМ СИЛЬНОЙ ЗАВИСИМОСТИ ВЯЗКОСТИ НЕФТИ ОТ ТЕМПЕРАТУРЫ Текст научной статьи по специальности «Физика»

CC BY
78
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭНЕРГЕТИКА / ТЕПЛООБМЕН / ВЯЗКОСТЬ / ГИДРОДИНАМИКА / НЕФТЕПРОДУКТЫ / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ / ЛАМИНАРНО-ТУРБУЛЕНТНЫЙ ПЕРЕХОД

Аннотация научной статьи по физике, автор научной работы — Курманова Д. Е., Джайчибеков Н. Ж.

Введение. Подогрев нефти и нефтепродуктов широко применяется для уменьшения энергопотерь при транспортировке. Течение в межтрубном пространстве теплообменника имеет сложный характер и зависит от многих факторов. Использование тонких трубок в теплообменных аппаратах геликоидного типа приводит к необходимости учета перехода режима течения от ламинарного к турбулентному. Традиционно используемые в численных расчетах полуэмпирические модели турбулентности не учитывают ламинарно-турбулентный переход. Разработан подход к определению эффективной длины теплообменного аппарата и температуры холодного теплоносителя на его выходе в случае сильной зависимости вязкости нефти от температуры с учетом возможности ламинарно-турбулентного перехода. В качестве нагреваемого теплоносителя рассмотрена нефть, а нагревающего компонента - вода. Метод. Новизна разработанного подхода заключается в применении модели турбулентности, учитывающей ламинарно-турбулентный переход, к расчету теплообменных аппаратов геликоидного типа. Для численного моделирования применены осредненные по Рейнольдсу уравнения Навье-Стокса, замкнутые при помощи модели турбулентности γ-Reθ t , учитывающей ламинарно-турбулентный переход. Основные результаты. Выполнено сравнение результатов численных расчетов с данными, полученными на основе метода среднелогарифмической разницы температур при постоянной и переменной вязкостях. В случае переменной вязкости нефти обнаружен переход ламинарного режима течения в турбулентный, который оказывает существенное влияние на эффективную длину теплообменного аппарата. Обсуждение. Результаты численных расчетов могут быть полезны при проектировании теплообменных аппаратов геликоидного типа.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Курманова Д. Е., Джайчибеков Н. Ж.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MODELING AND SIMULATION OF HEAT EXCHANGER WITH STRONG DEPENDENCE OF OIL VISCOSITY ON TEMPERATURE

Heating of oil and oil products is widely used to reduce energy losses during transportation. The flow in the annular space of the heat exchanger is complex and depends on many factors. The use of thin tubes in helicoid-type heat exchangers makes it necessary to take into account the transition of the flow regime from laminar to turbulent. The semi-empirical turbulence models traditionally used in numerical calculations do not take into account the laminar-turbulent transition. An approach is developed to determine the effective length of the heat exchanger and the temperature of the cold coolant at its outlet in the case of a strong dependence of oil viscosity on temperature, taking into account the possibility of a laminar-turbulent transition. Oil is considered as a heated coolant, and water is considered as a heating component. The novelty of the developed approach lies in the application of the turbulence model, which takes into account the laminar-turbulent transition, to the calculation of helicoid-type heat exchangers. For numerical simulation, the Reynolds-averaged Navier-Stokes equations are used which are closed using γ-Reθ t turbulence model that takes into account the laminar-turbulent transition. The results of numerical calculations are compared with the data obtained on the basis of the log-mean temperature difference method at constant and variable viscosity. In the case of variable oil viscosity, a transition from the laminar flow regime to the turbulent one is manifested which has a significant effect on the effective length of the heat exchanger. The results of numerical calculations can be useful in designing helicoid-type heat exchangers.

Текст научной работы на тему «МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛООБМЕННОГО АППАРАТА С УЧЕТОМ СИЛЬНОЙ ЗАВИСИМОСТИ ВЯЗКОСТИ НЕФТИ ОТ ТЕМПЕРАТУРЫ»

НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ июль-август 2023 Том 23 № 4 http://ntv.ifmo.ru/

I/ITMO SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ. МЕХАНИКИ И ОПТИКИ

July-August 2023 Vol. 23 No 4 http://ntv.ifmo.ru/en/

ISSN 2226-1494 (print) ISSN 2500-0373 (online)

doi: 10.17586/2226-1494-2023-23-4-820-827 УДК 536.24

Математическое моделирование теплообменного аппарата с учетом сильной зависимости вязкости нефти от температуры

Динара Есентаевна Курманова1, Нурболат Жумабекович Джайчибеков2Н

!>2 Евразийский национальный университет им. Л.Н. Гумилева, Астана, 010008, Республика Казахстан

1 [email protected], https://orcid.org/0009-0009-9787-7426

2 [email protected]и, https://orcid.org/0000-0002-3053-8288

Аннотация

Введение. Подогрев нефти и нефтепродуктов широко применяется для уменьшения энергопотерь при транспортировке. Течение в межтрубном пространстве теплообменника имеет сложный характер и зависит от многих факторов. Использование тонких трубок в теплообменных аппаратах геликоидного типа приводит к необходимости учета перехода режима течения от ламинарного к турбулентному. Традиционно используемые в численных расчетах полуэмпирические модели турбулентности не учитывают ламинарно-турбулентный переход. Разработан подход к определению эффективной длины теплообменного аппарата и температуры холодного теплоносителя на его выходе в случае сильной зависимости вязкости нефти от температуры с учетом возможности ламинарно-турбулентного перехода. В качестве нагреваемого теплоносителя рассмотрена нефть, а нагревающего компонента — вода. Метод. Новизна разработанного подхода заключается в применении модели турбулентности, учитывающей ламинарно-турбулентный переход, к расчету теплообменных аппаратов геликоидного типа. Для численного моделирования применены осредненные по Рейнольдсу уравнения Навье-Стокса, замкнутые при помощи модели турбулентности у-Яеи, учитывающей ламинарно-турбулентный переход. Основные результаты. Выполнено сравнение результатов численных расчетов с данными, полученными на основе метода среднелогарифмической разницы температур при постоянной и переменной вязкостях. В случае переменной вязкости нефти обнаружен переход ламинарного режима течения в турбулентный, который оказывает существенное влияние на эффективную длину теплообменного аппарата. Обсуждение. Результаты численных расчетов могут быть полезны при проектировании теплообменных аппаратов геликоидного типа. Ключевые слова

энергетика, теплообмен, вязкость, гидродинамика, нефтепродукты, численное моделирование, ламинарно-турбулентный переход

Ссылка для цитирования: Курманова Д.Е., Джайчибеков Н.Ж. Математическое моделирование теплообменного аппарата с учетом сильной зависимости вязкости нефти от температуры // Научно-технический вестник информационных технологий, механики и оптики. 2023. Т. 23, № 4. С. 820-827. doi: 10.17586/2226-1494-2023-23-4-820-827

Modeling and simulation of heat exchanger with strong dependence of oil viscosity on temperature Dinara E. Kurmanova1, Nurbolat Zh. Jaichibekov2®

L.N. Gumilyov Eurasian National University, Astana, 010008, Kazakhstan

1 [email protected], https://orcid.org/0009-0009-9787-7426

2 [email protected]®, https://orcid.org/0000-0002-3053-8288

Abstract

Heating of oil and oil products is widely used to reduce energy losses during transportation. The flow in the annular space of the heat exchanger is complex and depends on many factors. The use of thin tubes in helicoid-type heat exchangers makes it necessary to take into account the transition of the flow regime from laminar to turbulent. The semi-empirical turbulence models traditionally used in numerical calculations do not take into account the laminar-turbulent transition.

© KypMaHOBa ftE., ^^aHHHSeKOB H.®., 2023

An approach is developed to determine the effective length of the heat exchanger and the temperature of the cold coolant at its outlet in the case of a strong dependence of oil viscosity on temperature, taking into account the possibility of a laminar-turbulent transition. Oil is considered as a heated coolant, and water is considered as a heating component. The novelty of the developed approach lies in the application of the turbulence model, which takes into account the laminarturbulent transition, to the calculation of helicoid-type heat exchangers. For numerical simulation, the Reynolds-averaged Navier-Stokes equations are used which are closed using y-Re0i turbulence model that takes into account the laminarturbulent transition. The results of numerical calculations are compared with the data obtained on the basis of the log-mean temperature difference method at constant and variable viscosity. In the case of variable oil viscosity, a transition from the laminar flow regime to the turbulent one is manifested which has a significant effect on the effective length of the heat exchanger. The results of numerical calculations can be useful in designing helicoid-type heat exchangers. Keywords

energy, heat transfer, viscosity, hydrodynamics, oil products, simulation, laminar and turbulent transition

For citation: Kurmanova D.E., Jaichibekov N.Zh. Modeling and simulation of heat exchanger with strong dependence of oil viscosity on temperature. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2023, vol. 23, no. 4, pp. 820-827 (in Russian). doi: 10.17586/2226-1494-2023-23-4-820-827

Введение

Изменение качественного состояния сырьевой базы приводит к освоению и вовлечению в эксплуатацию месторождений нефти с высоким содержанием парафинов, смол, асфальтенов. Разработка таких месторождений требует применения нетрадиционных методов добычи нефти и ее подготовки к транспортировке. Светлые нефтепродукты (бензин, керосин) легко транспортируются по трубопроводам в любое время года и операции с ними не вызывают особых затруднений. Операции с темными нефтепродуктами (мазут, смазочные масла) и сырой нефтью вызывают значительные трудности, связанные тем, что темные нефтепродукты при понижении температуры воздуха становятся более вязкими, и их транспортировка без подогрева невозможна. Для трубопроводной транспортировки нефти и нефтепродуктов используется подход, основанный на регулировании реологических свойств нефти, например, при помощи нагрева нефти с ее последующей транспортировкой по трубопроводу с повышенной теплоизоляцией (горячая перекачка нефти). В некоторых случаях увеличение вязкости нефти при понижении температуры приводит к недопустимым напряжениям на стенках трубы и остановке транспортировки.

Процессы теплообмена осуществляются в теплоо-бменных аппаратах (heat exchanger) различных типов и конструкций [1-5]. Для подогрева применяются различные теплоносители, например, горячая вода или водяной пар. Потребление энергии является одним из важных факторов, который оказывает существенное влияние на конструкцию теплообменника [6]. В нефтегазовой отрасли применение находят кожухотрубные теплообменники, которые обеспечивают хорошие характеристики производительности в широком диапазоне рабочих условий, высокую надежность и низкую стоимость.

Состав нефти (в частности, содержание асфальте-нов, смол, парафинов) оказывает существенное влияние на зависимость вязкости от температуры [7, 8]. Эмпирические формулы, описывающие изменение кинематической вязкости в зависимости от температуры, имеют вид различных функций (экспоненциальные, полиноминальные, степенные и т. д.), которые характеризуются наличием коэффициентов, зависящих от

свойств жидкости. Постоянные коэффициенты определяются на основе значений измеренных кинематических вязкостей в экспериментальных точках. Для расчетов термодинамических параметров нефти, газовых конденсатов и их фракций используется обобщенное уравнение состояния Ли-Кесслера [9]. В работе [10] проведены исследования зависимости кинематической вязкости нефти и смесей нефтей от температуры, а также проанализированы существующие формулы для расчета кинематической вязкости нефти в магистральных трубопроводах.

Течение жидкости в межтрубном пространстве теплообменника имеет сложный характер и зависит от многих факторов. Численное моделирование теплообмена в теплообменных устройствах различных конструкций проведено в работах [11, 12]. Результаты численных расчетов применяются для поиска оптимальных способов интенсификации процессов теплообмена [13-15]. Полученные результаты указывают на уменьшение влияния вязкости перекачиваемой нефти на гидравлическую характеристику трубопровода при перекачке на развитых турбулентных режимах.

Сравнение точности различных моделей турбулентности, используемых для замыкания уравнений Рейнольдса, выполнено в [16-19]. В основном в расчетах используются k-e и k-ю модели турбулентности, а также Shear Stress Transport (SST) модель k-ю. Полученные результаты сравнения моделей позволяют добиться удовлетворительного совпадения результатов расчетов с данными промышленных экспериментов. В результате возможно решить задачи управления процессами, повысить эффективность производства и определить оптимальные режимные параметры технологических процессов. Вместе с тем имеющиеся расчеты с использованием двухпараметрических моделей турбулентности не учитывают ламинарно-тур-булентный переход, что сказывается на определении эффективной длины теплообменного аппарата [20].

В классических теплообменных аппаратах пучок труб для одного теплоносителя помещается внутрь кожуха, по которому движется другой теплоноситель. В конструкции геликоидных теплообменников применяются профилированные трубки и ребра винтового профиля, с помощью которых улучшаются условия теплообмена. Трубки в таких аппаратах имеют малый

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

В настоящей работе разработана математическая модель теплообменного аппарата, учитывающая лами-нарно-турбулентный переход. Выбрана схема теплообменника «труба в трубе» с тонкой и гладкой внутренней трубкой. Приведен метод расчета теплообменника прямоточного типа, в котором во внутреннем трубопроводе рабочей жидкостью является нефть (холодный теплоноситель), а во внешней трубе — вода (горячий теплоноситель). Расчеты проведены для модельной конструкции теплообменного аппарата на базе двух подходов, основанных на методе среднелогарифми-ческой разницы температур (Log-Mean Temperature Difference, LMTD) при постоянной и переменной вязкости и на основе средств вычислительной гидродинамики (Computational Fluid Dynamics, CFD). Осредненные по Рейнольдсу уравнения Навье-Стокса замкнуты при помощи модели турбулентности, учитывающей лами-нарно-турбулентный переход. Выполнено сравнение данных, полученных в рамках двух подходов.

Математическая модель

В расчетах использованы методы: среднелогариф-мической разницы температур при постоянной и переменной вязкости; вычислительной гидродинамики.

Метод среднелогарифмической разницы температур. Для оценки тепловых потоков от горячего теплоносителя к холодному применяется модель теплоносителя с постоянной или переменной по длине вязкостью, основанная на использовании среднелога-рифмической разницы температур. В случае сильной зависимости вязкости от температуры теплообменник разбивается на элементарные участки по длине. На каждом участке делается предположение о малом изменении вязкости.

Число Нуссельта зависит от режима течения (ламинарный, переходный, турбулентный) и режима теплообмена (нагрев или охлаждение). Локальные числа Нуссельта в ламинарном и турбулентном режимах течения рассчитываются при помощи соотношений, приведенных в работах [21-23]. Для кольцевого канала в турбулентном режиме течения используется соотношение, как и для течения в трубе, но со своим эквивалентным гидравлическим диаметром. Переход между различными критериальными соотношениями, соответствующими ламинарному, переходному и турбулентному режимам течения, реализуется в зависимости от числа Рейнольдса, которое находится по среднемас-совой скорости (при Re ~ 2 103 и Re ~ 104).

Метод вычислительной гидродинамики. Результаты расчетов при помощи метода среднелогарифмиче-ской разницы температур сравниваются с данными, полученными методами вычислительной гидродинамики. Течение считается стационарным и осесимметричным,

нефть — ньютоновской жидкостью с постоянной плотностью и переменными теплофизическими свойствами. Расчеты проводятся при помощи численного решения осредненных по Рейнольдсу уравнений Навье-Стокса (Reynolds-Averaged Navier-Stokes, RANS) для вязкой несжимаемой жидкости, замкнутой благодаря модели турбулентности, учитывающей ламинарно-турбулент-ный переход [24].

Модель турбулентности SST k-ю разработана для эффективного сочетания надежной и точной модели k-ю в пристеночной области и модели k-e в свободном потоке [25, 26]. Для переключения между моделями применяется специальная функция, которая принимает единичное значение в пристеночной области (стандартная модель k-ю) и нулевое значение вдали от стенки (модель k-e).

Модель, учитывающая ламинарно-турбулентный переход (Local-Correlation Transition Model, y-Re0i transition model), основана на сочетании уравнений SST k-ю модели турбулентности с двумя дополнительными уравнениями переноса для параметра перемежаемости у и критического числа Рейнольдса Re0i, построенного по толщине потери импульса [27, 28]. Для упрощения модели уравнение для Reef не рассматривается, а в уравнении для параметра перемежаемости делается предположение о малости конвективных слагаемых [29]. Такой подход приводит к алгебраическим соотношениям для нахождения параметра перемежаемости.

Модель турбулентности y-Reei, используемая в расчетах, была разработана для решения задач внешней аэродинамики. В дальнейшем модель успешно применялась для моделирования широкого круга течений. В частности, в работе [30] проведено сравнение различных моделей турбулентности (k-e, k-ю, k-ю SST, v2-f, y-Reef) применительно к течениям в каналах и показаны преимущества модели y-Reei перед другими моделями турбулентности.

В научных работах применяются различные зависимости вязкости от температуры. В нефтяной отрасли при расчете кинематической вязкости v, зависящей от температуры T, используется формула Вальтера [10]

lg[lg(v + 0,8)] = a + blgT,

(1)

где а и Ь — эмпирические коэффициенты, определяемые для исследуемой жидкости экспериментальным путем. Коэффициенты а и Ь в формуле (1) находятся

из соотношении

b =

a = lg[lg(v i + 0,8)] - blgT, lg[lg(v1 + 0,8)]-lg[lg(v2 + 0,8)] Ig^-lg^ ,

где V! и V2 — значения кинематическои вязкости жидкости при температурах и Т2.

Численное моделирование

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

Уравнения решены методом конечных разностей на интервале x Е [0, L], где L — длина теплообменника, что позволяет применить разработанный подход для расчета прямоточных и противопоточных теплообменников. Для стабилизации итерационного процесса при линеаризации системы применен метод нижней релаксации. Длина промежутка интегрирования L заранее неизвестна, которая определена методом Ньютона.

Выполнена дискретизация осредненных по Рей-нольдсу уравнений Навье-Стокса и уравнений модели турбулентности с помощью методов конечных объемов на неструктурированных сетках и коррекции давления Semi-Implicit Method for Pressure Linked Equations (SIMPLE) [31]. Дискретизация невязких потоков осуществлена с использованием монотонной противо-поточной схемы для законов сохранения (Monotonic Upstream Schemes for Conservation Laws, MUSCL), а вязких потоков — центрированная схема второго порядка точности. Схема MUSCL позволила повысить порядок аппроксимации по пространственным переменным без потери монотонности решения, соответствует условию Total Variation Diminishing (TVD) и представляет собой комбинацию центрированных конечных разностей второго порядка и диссипативного члена, для переключения между которыми служит ограничитель потока, построенный на основе характеристических переменных. Для решения системы разностных уравнений использован геометрический многосеточный метод [31].

Для расчетов ядра потока выбрана сетка, состоящая из 19 461 ячеек, из которых 500 х 24 ячеек размещены в области, заполненной нефтью, 500 х 5 ячеек — в области из стали, а 500 х 13 — в области, заполненной водой. Сгущение ячеек сетки произведено около стенок трубы таким образом, чтобы y+ <2, где y+ - безразмерная пристеночная координата. Для этого в пристеночном слое размещено 15 узлов сетки, расстояние между которыми изменяется по закону геометрической прогрессии.

Результаты и обсуждение

В качестве примера рассмотрим нефть Узеньского нефтегазового месторождения, расположенного в Мангистауской области Казахстана. Залежи нефти находятся на глубине 0,9-2,4 км. Нефть имеет следующие параметры: плотность составляет 844-874 кг/м3; вязкость — 3,4-8,15 мПас; содержание серы — 0,16-2 %, парафинов 16-22 % и смол — 8-20 %.

Выполним сравнение результатов расчетов по формуле Вольтера с экспериментальными значениями динамической вязкости д = pv, где р — плотность (рис. 1) для нефти Узеньского месторождения. Температура изменяется от 283 до 383 К.

Для расчета эффективной длины теплообменника и температуры горячего теплоносителя на выходе заданы следующие параметры: Tci = 303 К и Tco = 328 К (входная и выходная температуры холодного теплоносителя); Thi = 423 К (входная температура горячего теплоносителя); vc = 4 мс и vh = 4 мс (скорости холодного и горячего теплоносителя); Gc = 0,3814 кг/с и Gh = 0,6386 кг/с

ц, Па-с

6

4

2

0

283 303 323 343 363 Г, К

Рис. 1. Зависимость динамической вязкости нефти Узеньского месторождения от температуры.

Экспериментальные данные [20] (треугольные значки) и расчеты по формуле Вальтера (1) (сплошная линия) Fig. 1. The dependence of the dynamic viscosity of oil from the

Uzen field vs. temperature. Triangular symbols correspond to experimental data [20], solid line corresponds to calculations by the Walther formula (1)

(расходы горячего и холодного теплоносителя); геометрические характеристики (толщина стенок трубок 5W = 1 мм, внутренний и внешний диаметры трубок d = 12 мм и dc = 14 мм, внутренний диаметр кожуха D = 20 мм); физические свойства материала трубок.

В результате решения уравнения теплового баланса методом конечных разностей получено, что длина теплообменника составляет L = 5,28 м (для постоянной вязкости) и L = 4,26 м (для переменной вязкости). При этом температура горячего теплоносителя на выходе достигла 416 К.

Выполним сравнение полученных результатов с данными численного моделирования. На рис. 2, a показаны распределения среднемассовой температуры нефти (холодного теплоносителя) (Tm) вдоль длины теплообменника, полученные при помощи метода конечных разностей и на основе численного моделирования. Среднемассовая температура нефти возрастает вдоль длины за счет нагревания от источника тепла (горячего теплоносителя). Результаты аналитических и численных расчетов достаточно хорошо согласуются между собой.

Распределения среднемассовой температуры воды (горячего теплоносителя) вдоль длины, полученные на основе аналитических и численных расчетов, показаны на рис. 2, b. Среднемассовая температура воды убывает вдоль длины трубы за счет передачи тепла от нее к холодному теплоносителю. Следует отметить хорошее согласование результатов расчетов, полученных на основе различных подходов.

На рис. 2 видны характерные изменения кривизны линий на расстоянии порядка 2,5 м от входного сечения, где заметны резкие изменения градиентов температур. Изменения происходят на расстоянии, на котором ламинарный режим течения переходит в турбулентный. На рис. 3 приведены графики изменения числа Рейнольдса вдоль длины теплоносителей, полученные при аналитическом расчете с моделью переменной

д

к

Рис. 2. Распределения среднемассовых температур нефти (а) и воды (b) вдоль длины теплообменника. Результаты, полученные на основе теоретического подхода (сплошная линия) и численных расчетов (пунктирная линия с треугольными

значками)

Fig. 2. Distribution of average mass temperatures of oil (a) and water (b) along the length of the heat exchanger. The results obtained on the basis of a theoretical approach (solid line) and numerical calculations (dotted line with triangular icons)

вязкости. Для воды число Рейнольдса практически не меняется при постоянном значении вязкости и незначительном изменении температуры. Обратная картина наблюдается для холодного теплоносителя (нефти). Число Рейнольдса существенно возрастает, что связано с резким уменьшением вязкости нефти при ее практически неизменной плотности. Отметим, что наблюдается переходный участок на расстоянии порядка от 1,85 до 3,8 м, где происходит переход ламинарного течения в турбулентный.

На рис. 4 приведен график изменения числа Рейнольдса нефти (холодного теплоносителя) вдоль

длины теплообменного аппарата при переменных значениях вязкости нефти (линия 1). Видно, что значения числа Рейнольдса резко возрастают вдоль продольной координаты, что объясняется уменьшением вязкости нефти за счет ее нагревания. На рисунке приведен также график поведения числа Рейнольдса нефти при средней температуре теплоносителя (линия 2). При постоянной вязкости число Рейнольдса холодного теплоносителя постоянно и не может правильно охарактеризовать теплообмен в данном теплообменном аппарате.

Re 106Г

102

10° I-'-■-'-'-'-■-'-'-

0 12 3 4

L, м

Рис. 3. Распределения чисел Рейнольдса по длине, полученные при аналитическом расчете для воды (линия 1)

и нефти (линия 2) Fig. 3. Length distributions of Reynolds numbers obtained in an analytical calculation. Line 1 corresponds to water and line 2 to oil

Re-104

L, m

Рис. 4. Распределения чисел Рейнольдса нефти вдоль длины теплоносителя при численном расчете при переменной вязкости (линия 1) и теоретическом расчете (линия 2) Fig. 4. Distributions of oil Reynolds numbers along length. Line 1 corresponds to the numerical calculation at variable viscosity, and line 2 corresponds to the theoretical calculation

Заключение

Снижение вязкости нефти при нагреве является одним из способов повышения энергоэффективности процесса перекачки высоковязкой нефти при добыче и транспортировке. Численное моделирование позволяет сократить сроки и трудоемкость испытаний новых конструкций теплообменных аппаратов в нефтегазовой отрасли.

Проведено численное моделирование кожухотруб-чатого теплообменника средствами вычислительной гидродинамики. Получены распределения скорости и

температуры в теплообменном аппарате. Выполнено сравнение результатов численных расчетов с данными, полученными на основе метода среднелогариф-мической разницы температур в случае постоянной и переменной вязкости. При переменной вязкости нефти осуществлен переход от ламинарного режима в турбулентный. Результаты эксперимента показали, что модель с постоянной вязкостью приводит к заниженной примерно на 20 % длине теплообменного аппарата по сравнению с расчетами, учитывающими зависимость вязкости нефти от температуры.

Литература

1. Nakaso K., Mitani H., Fukai J. Convection heat transfer in a shell-and-tube heat exchanger using sheet fins for effective utilization of energy // International Journal of Heat and Mass Transfer. 2015. V. 82. P. 581-587. https://doi.org/10.1016/j.ijheatmasstransfer.2014.11.033

2. Wang Y., Gu X., Jin Z., Wang K. Characteristics of heat transfer for tube banks in crossflow and its relation with that in shell-and-tube heat exchangers // International Journal of Heat and Mass Transfer. 2016. V. 93. P. 584-594. https://doi.org/10.1016/j. ijheatmasstransfer.2015.10.018

3. Araavind S., Athreya A.S. CFD analysis of shell and tube heat exchanger for pre-heating of biodiesel // International Journal of Advance Research in Science and Engineering. 2017. V. 6. N 7. P. 687-693.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4. Abda A.A., Kareema M.Q., Najib S.Z. Performance analysis of shell and tube heat exchanger: parametric study // Case Studies in Thermal Engineering. 2018. V. 12. P. 563-568. https://doi.org/10.1016/j. csite.2018.07.009

5. Лычаков В.Д., Егоров М.Ю., Щеглов А.А., Сивоволов А.С., Матяш А.С., Балунов Б.Ф. Анализ теплоотдачи сборок оребрен-ных теплообменных элементов // Теплоэнергетика. 2022. № 3. С. 63-71.

6. Darbandi M., Abdollahpour M.-S., Hasanpour-Matkolaei M. A new developed semi-full-scale approach to facilitate the CFD simulation of shell and tube heat exchangers // Chemical Engineering Science. 2021. V. 245. P. 116836. https://doi.org/10.1016/j.ces.2021.116836

7. Bizhan K.G., Mohammad R.D., Hossein P. Prediction of kinematic viscosity of petroleum fractions using artificial neural networks // Iranian Journal of Oil and Gas Science and Technology. 2014. V. 3. N 2. P. 51-65.

8. Boda M.A., Bhasagi P.N., Sawade A.S., Andodgi R.A. Analysis of kinematic viscosity for liquids by varying temperature // International Journal of Innovative Research in Science, Engineering and Technology. 2015. V. 4. N 4. P. 1951-1954. https://doi.org/10.15680/ IJIRSET.2015.0404020

9. Lee B.I., Kesler M.G. A generalized thermodynamic correlation based on three-parameter corresponding states // AIChE Journal. 1975. V. 21. N 3. P. 510-527. https://doi.org/10.1002/aic.690210313

10. Аралов О.В., Буянов И.В., Саванин А.С., Иорданский Е.И. Исследование методов расчета кинематической вязкости нефти в магистральном нефтепроводе // Наука и технологии трубопроводного транспорта нефти и нефтепродуктов. 2017. Т. 7. № 5. С. 97-105.

11. Yogesh S.S., Selvaraj A.S., Ravi D.K., Rajagopal T.K.R. Heat transfer and pressure drop characteristics of inclined elliptical fin tube heat exchanger of varying ellipticity ratio using CFD code // International Journal of Heat and Mass Transfer. 2018. V. 119. P. 26-39. https://doi. org/10.1016/j.ijheatmasstransfer. 2017.11.094

12. Chen K., Mohammed H.I., Mahdi J.M., Rahbari A., Cairns A., Talebizadehsardari P. Effects of non-uniform fin arrangement and size on the thermal response of a vertical latent heat triple-tube heat exchanger // Journal of Energy Storage. 2022. V. 45. P. 103723. https://doi.org/10.1016/j.est.2021.103723

13. Osley W.G., Droegemueller P., Ellerby P. CFD investigation of heat transfer and flow patterns in tube side laminar flow and the potential for enhancement // Chemical Engineering Transactions. 2013. V. 35. P. 997-1002. https://doi.org/10.3303/CET1335166

References

1. Nakaso K., Mitani H., Fukai J. Convection heat transfer in a shell-and-tube heat exchanger using sheet fins for effective utilization of energy. International Journal of Heat and Mass Transfer, 2015, vol.82, pp. 58 1-587. https ://doi. org/1 0. 1 0 1 6/j. ijheatmasstransfer.2014.11.033

2. Wang Y., Gu X., Jin Z., Wang K. Characteristics of heat transfer for tube banks in crossflow and its relation with that in shell-and-tube heat exchangers. International Journal of Heat and Mass Transfer, 2016, vol. 93, pp. 584-594. https://doi.org/10.1016/j. ijheatmasstransfer.2015.10.018

3. Araavind S., Athreya A.S. CFD analysis of shell and tube heat exchanger for pre-heating of biodiesel. International Journal of Advance Research in Science and Engineering, 2017, vol. 6, no. 7, pp. 687-693.

4. Abda A.A., Kareema M.Q., Najib S.Z. Performance analysis of shell and tube heat exchanger: parametric study. Case Studies in Thermal Engineering, 2018, vol. 12, pp. 563-568. https://doi.org/10.1016/j. csite.2018.07.009

5. Lychakov V.D., Scheglov A.A., Sivovolov A.S., Matyash A.S., Balunov B.F., Egorov M.Y. An analysis of heat transfer in assemblies of finned heat-exchange elements. Thermal Engineering, 2022, vol. 69, no. 3, pp. 202-209. https://doi.org/10.1134/ S0040601522030089

6. Darbandi M., Abdollahpour M.-S., Hasanpour-Matkolaei M. A new developed semi-full-scale approach to facilitate the CFD simulation of shell and tube heat exchangers. Chemical Engineering Science, 2021, vol. 245, pp. 116836. https://doi.org/10.1016/j.ces.2021.116836

7. Bizhan K.G., Mohammad R.D., Hossein P. Prediction of kinematic viscosity of petroleum fractions using artificial neural networks. Iranian Journal of Oil and Gas Science and Technology, 2014, vol. 3, no. 2, pp. 51-65.

8. Boda M.A., Bhasagi P.N., Sawade A.S., Andodgi R.A. Analysis of kinematic viscosity for liquids by varying temperature. International Journal of Innovative Research in Science, Engineering and Technology, 2015, vol. 4, no. 4, pp. 1951-1954. https://doi. org/10.15680/IJIRSET. 2015.0404020

9. Lee B.I., Kesler M.G. A generalized thermodynamic correlation based on three-parameter corresponding states. AIChE Journal, 1975, vol. 21, no. 3, pp. 510-527. https://doi.org/10.1002/aic.690210313

10. Aralov O.V., Buyanov I.V., Savanin A.S., Iordansky E.I. Research of methods for oil kinematic viscosity calculation in the oil-trunk pipeline. Science & Technologies: Oil and Oil Products Pipeline Transportation, 2017, vol. 7, no. 5, pp. 97-105. (in Russian)

11. Yogesh S.S., Selvaraj A.S., Ravi D.K., Rajagopal T.K.R. Heat transfer and pressure drop characteristics of inclined elliptical fin tube heat exchanger of varying ellipticity ratio using CFD code. International Journal of Heat and Mass Transfer, 2018, vol. 119, pp. 26-39. https:// doi.org/10.1016/j.ijheatmasstransfer.2017.11.094

12. Chen K., Mohammed H.I., Mahdi J.M., Rahbari A., Cairns A., Talebizadehsardari P. Effects of non-uniform fin arrangement and size on the thermal response of a vertical latent heat triple-tube heat exchanger. Journal of Energy Storage, 2022, vol. 45, pp. 103723. https://doi.org/10.1016Zj.est.2021.103723

13. Osley W.G., Droegemueller P., Ellerby P. CFD investigation of heat transfer and flow patterns in tube side laminar flow and the potential

14. Karar O., Emani S., Gounder S.M., Myo Thant M.M., Mukhtar H., Sharifpur M., Sadeghzadeh M. Experimental and numerical investigation on convective heat transfer in actively heated bundle- 14. pipe // Engineering Applications of Computational Fluid Mechanics.

2021. V. 15. N 1. P. 848-864. https://doi.org/10.1080/19942060.202 1.1920466

15. Rana S., Zunaid M., Kumar R. CFD approach for the enhancement of thermal energy storage in phase change material charged heat exchanger // Case Studies in Thermal Engineering. 2022. V. 33. 15. P. 101921. https://doi.org/10.1016/j.csite.2022.101921

16. Allouche Y., Varga S., Bouden C., Oliveira A.C. Validation of a CFD model for the simulation of heat transfer in a tubes-in-tank PCM storage unit // Renewable Energy. 2016. V. 89. P. 371-379. https:// 16. doi.org/10.1016/j.renene.2015.12.038

17. Balaji D., Prakash L.S.S. CFD analysis of a pressure drop in a staggered tube bundle for a turbulent cross flow // International Advanced Research Journal in Science, Engineering and Technology. 17. 2016. V. 3. N 2. P. 35-40. https://doi.org/10.17148/ IARJSET. 2016.3209

18. Czarnota T., Wagner C. Turbulent convection and thermal radiation in a cuboidal Rayleigh-Benard cell with conductive plates // International Journal of Heat and Fluid Flow. 2016. V 57. P. 150-172. 18. https://doi.org/10.1016/j.ijheatfluidflow.2015.10.006

19. Mohanan A.K., Prasad B.V., Vengadesan S. Flow and heat transfer characteristics of a cross-flow heat exchanger with elliptical tubes //

Heat Transfer Engineering. 2021. V. 42. N 21. P. 1846-1860. https:// 19. doi.org/10.1080/01457632.2020.1826742

20. Тугунов П.И., Новоселов В.Ф., Коршак А.А., Шаммазов А.М. Типовые расчеты при проектировании и эксплуатации нефтебаз

и нефтепроводов. М.: Дизайн Полиграф Сервис, 2002. 234 с. 20.

21. Кириллов П.Л., Юрьев Ю.С., Бобков В.П. Справочник по тепло-гидравлическим расчетам (ядерные реакторы, теплообменники, парогенераторы). М.: Энергоатомиздат, 1990. 360 с.

22. Петухов Б.С., Генин Л.Г., Ковалев С.А. Теплообмен в ядерных 21. энергетических установках. М.: Атомиздат, 1974. 408 с.

23. Кутателадзе С.С. Теплопередача и гидродинамическое сопротивление. М.: Энергоатомиздат, 1990. 367 с.

24. Menter F.R. Two-equation eddy-viscosity turbulence models for 22. engineering applications // AIAA Journal. 1994. V. 32. N 8. P. 15981605. https://doi.org/10.2514Z3.12149 23.

25. Menter F., Esch T., Kubacki S. Transition modelling based on local variables // Engineering Turbulence Modelling and Experiments. 24. Elsevier, 2002. P. 555-564. https://doi.org/10.1016/b978-008044114-6/50053-3

26. Langtry R.B., Menter F.R. Correlation-based transition modeling for 25. unstructured parallelized computational fluid dynamics codes // AIAA Journal. 2009. V. 47. N 12. P. 2894-2906. https://doi. org/10.2514/1.42362

27. Coder J.G., Maughmer M.D. Computational fluid dynamics 26. compatible transition modeling using an amplification factor transport equation // AIAA Journal. 2014. V. 52. N 11. P. 2506-2512. https:// doi.org/10.2514/1j052905

28. Menter F.R., Smirnov P.E., Liu T., Avancha R. A one-equation local 27. correlation-based transition model // Flow, Turbulence and Combustion. 2015. V. 95. N 4. P. 583-619. https://doi.org/10.1007/ s10494-015-9622-4

29. Gorji S., Seddighi M., Ariyaratne C., Vardy A.E., O'Donoghue T., 28. Pokrajac D., He S. A comparative study of turbulence models in a transient channel flow // Computers and Fluids. 2014. V. 83. P. 111—

123. https://doi.org/10.1016/j.compfluid.2013.10.037

30. Volkov K. Numerical analysis of Navier-Stokes equations on 29. unstructured meshes // Handbook on Navier-Stokes Equations: Theory and Analysis. Nova Science, 2016. P. 365-442.

31. Volkov K. Multigrid and preconditioning techniques in CFD applications // CFD Techniques and Thermo-Mechanics Applications. 30. Springer, 2018. P. 83-149. https://doi.org/10.1007/978-3-319-70945-

1 6

for enhancement. Chemical Engineering Transactions, 2013, vol. 35, pp. 997-1002. https://doi.org/10.3303/CET1335166 Karar O., Emani S., Gounder S.M., Myo Thant M.M., Mukhtar H., Sharifpur M., Sadeghzadeh M. Experimental and numerical investigation on convective heat transfer in actively heated bundle-pipe. Engineering Applications of Computational Fluid Mechanics, 2021, vol. 15, no. 1, pp. 848-864. https://doi.org/10.1080/19942060 .2021.1920466

Rana S., Zunaid M., Kumar R. CFD approach for the enhancement of thermal energy storage in phase change material charged heat exchanger. Case Studies in Thermal Engineering, 2022, vol. 33, pp. 101921. https://doi.org/10.1016/j.csite.2022.101921 Allouche Y., Varga S., Bouden C., Oliveira A.C. Validation of a CFD model for the simulation of heat transfer in a tubes-in-tank PCM storage unit. Renewable Energy, 2016, vol. 89, pp. 371-379. https:// doi.org/10.1016/j.renene.2015.12.038

Balaji D., Prakash L.S.S. CFD analysis of a pressure drop in a staggered tube bundle for a turbulent cross flow. International Advanced Research Journal in Science, Engineering and Technology, 2016, vol. 3, no. 2, pp. 35-40. https://doi.org/10.17148/ IARJSET.2016.3209

Czarnota T., Wagner C. Turbulent convection and thermal radiation in a cuboidal Rayleigh-Benard cell with conductive plates.

International Journal ofHeat and Fluid Flow, 2016, vol. 57, pp. 150172. https://doi.org/10.1016/j.ijheatfluidflow.2015.10.006 Mohanan A.K., Prasad B.V., Vengadesan S. Flow and heat transfer characteristics of a cross-flow heat exchanger with elliptical tubes. Heat Transfer Engineering, 2021, vol. 42, no. 21, pp. 1846-1860. https://doi.org/10.1080/01457632.2020.1826742 Tugunov P.I., Novoselov V.F., Korshak A.A., Shammazov A.M. Standard Calculations for the Design and Operation Of Petrolium Storage Depots and Oil Pipelines. Moscow, Dizajn Poligraf Servis Publ., 2002, 234 p. (in Russian)

Kirillov P.L., Iurev Iu.S., Bobkov V.P. Handbook of Thermal-Hydraulic Calculations (Nuclear-Power Reactors, Heat Exchangers, Steam Generating Units). Moscow, Jenergoatomizdat, 1990, 360 p. (in Russian)

Petukhov B.S., Genin L.G., Kovalev S.A. Heat transfer in nuclear power plants. Moscow, Atomizdat Publ., 1974, 408 p. (in Russian) Kutateladze S.S. Heat Transfer and Hydrodynamical Resistance. Moscow, Jenergoatomizdat Publ., 1990, 367 p. (in Russian) Menter F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal, 1994, vol. 32, no. 8, pp. 1598-1605. https://doi.org/10.2514Z3.12149 Menter F., Esch T., Kubacki S. Transition modelling based on local variables. Engineering Turbulence Modelling and Experiments, Elsevier, 2002, pp. 555-564. https://doi.org/10.1016/b978-008044114-6/50053-3

Langtry R.B., Menter F.R. Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes. AIAA Journal, 2009, vol. 47, no. 12, pp. 2894-2906. https://doi. org/10.2514/1.42362

Coder J.G., Maughmer M.D. Computational fluid dynamics compatible transition modeling using an amplification factor transport equation. AIAA Journal, 2014, vol. 52, no. 11, pp. 2506-2512. https:// doi.org/10.2514/1j052905

Menter F.R., Smirnov P.E., Liu T., Avancha R. A one-equation local correlation-based transition model. Flow, Turbulence and Combustion, 2015, vol. 95, no. 4, pp. 583-619. https://doi. org/10.1007/s10494-015-9622-4

Gorji S., Seddighi M., Ariyaratne C., Vardy A.E., O'Donoghue T., Pokrajac D., He S. A comparative study of turbulence models in a transient channel flow. Computers and Fluids, 2014, vol. 83, pp. 111123. https://doi.org/10.1016/j.compfluid.2013.10.037 Volkov K. Numerical analysis of Navier-Stokes equations on unstructured meshes. Handbook on Navier—Stokes Equations: Theory and Analysis. Nova Science, 2016, pp. 365-442. Volkov K. Multigrid and preconditioning techniques in CFD applications. CFD Techniques and Thermo-Mechanics Applications, Springer, 2018, pp. 83-149. https://doi.org/10.1007/978-3-319-70945-1 6

Авторы

Курманова Динара Есентаевна — докторант, Евразийский национальный университет им. Л.Н. Гумилева, Астана, 010008, Республика Казахстан, https://orcid.org/0009-0009-9787-7426, dikonya89_29@ mail.ru

Джайчибеков Нурболат Жумабекович — доктор физико-математических наук, профессор, профессор, Евразийский национальный университет им. Л.Н. Гумилева, Астана, 010008, Республика Казахстан, sc 57195809348, https://orcid.org/0000-0002-3053-8288, [email protected]

Authors

Dinara E. Kurmanova — Doctoral Student, L.N. Gumilyov Eurasian National University, Astana, 010008, Kazakhstan, https://orcid.org/0009-0009-9787-7426, [email protected]

Nurbolat Zh. Jaichibekov — D.Sc. (Physics & Mathematics), Full Professor, L.N. Gumilyov Eurasian National University, Astana, 010008, Kazakhstan, sc 57195809348, https://orcid.org/0000-0002-3053-8288, [email protected]

Статья поступила в редакцию 10.03.2023 Одобрена после рецензирования 11.05.2023 Принята к печати 16.07.2023

Received 10.03.2023

Approved after reviewing 11.05.2023

Accepted 16.07.2023

Работа доступна по лицензии Creative Commons «Attribution-NonCommercial»

i Надоели баннеры? Вы всегда можете отключить рекламу.