Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,580

РАЗВИТИЕ ТУРБУЛЕНТНОЙ НЕИЗОТЕРМИЧЕСКОЙ СТРУИ В ОБЪЕМЕ АМБАРНОЙ НЕФТИ

Айсаев С.У. 1 Диярова Л.Д. 1 Жапбасбаев У.К. 2 Аринов Е. 3
1 Каспийский государственный университет технологий и инжиниринга
2 Казахстанско-Британский технический университет
3 Жезказганский государственный университет
В работе рассматривается математическая модель развития турбулентной неизотермической струи в объеме амбарной нефти. Предполагается, что жидкость струи имеет такие же физико-химические свойства, как и амбарная нефть, реологические свойства которой удовлетворяют модели ньютоновской жидкости. По мере развития струи в резервуаре ее скорость затухает. Тепло, вносимое горячей струей, подогревает окружающую высоковязкую среду и вызывает снижение ее вязкости, способствуя вовлечению амбарной нефти в движение. Движение струи считается турбулентным, а процесс нагрева амбарной нефти вышедшим на регулярный режим, т.е. количество вносимого в амбар тепла равно, отводимого через ограничивающие поверхности. Неизотермическое течение жидкости в объеме амбарной нефти описывается системой уравнений Навье−Стокса, осредненных по Рейнольдсу. Уравнение переноса тепла учитывает диссипацию механической энергии.
амбарная нефть
турбулентная неизотермическая струя
математическая модель
численное моделирование
1. Патент № 10116 на изобретение. Способ сбора амбарной нефти и устройство для его осуществления (авторы: Ершин Ш.А., Жапбасбаев У.К., Айсаев С.У., Утегалиев С.А, Хаиров Г.Б.), 9.10.2002 г.
2. Абрамович Г.Н. и др. Теория турбулентных струй. – М.: Наука, 1984. – 715 с.
3. Pope S.P. Turbulent Flows. Cambridge: Cambridge University Press, 2005. 771 p.
4. Chen H.C. and Patel V.C. Near-wall turbulence models for complex flows including separation. AIAA J. – 1988. –V.26. – pp. 641-648.
5. Menter F.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA J. –1993. – V.32. –pp. 1598-1604.
6. Craft T.J., Gerasimov A.V., Iacovides H., Kidger J.W. and Launder B.E. The Negetively Buoyant Turbulent wall Jet: Performance of Alternative Options in RANS Modelling. Proc. 4th Int. Symp. On Turbulence, Heat and Mass Transfer, Antaly, Turkey. –2003. –pp. 813-820.
7. He S., Xu Z. and Jackson J.D. An experimental investigation of buoyancy – opposed wall jet flow. Int. J. Heat and Fluid Flow. – 2002. V. 23. – pp. 487-496.
8. Chien, K.Y. Predictions of channel and boundary layer flows with low Reynolds-number two –equation model of turbulence. AIAA J. - 1982. –V.20, №1. –pp. 33-38.
9. Бэтчелор Дж. Введение в динамику жидкости. – М.: Мир, 1973. –757 с.
10. Госмен А.Д., Пан В.И., Ранчел А.К., Сполдинг Д.В., Вольфштейн М. Численные методы течений вязкой жидкости: Пер. с англ. – М.: Мир, 1972. – 320 с.
11. Роуч П. Вычислительная гидродинамика: Пер. с англ. – М.: Мир, 1980. 616 с.
12. Chung, T.J. Computational Fluid Dynamics. Cambridge: Cambridge University Press, 2002. – 520 р.
13. Kuznetsov, G.V., Sheremet, M.A. Numerical simulation of turbulent natural convection in a rectangular enclosure having finite thickness walls // Int. J. Heat and Mass Transfer. 2010. V.53. P. 163–177.
14. Gardon R., Akifrat J.C. The role of turbulence in determining the heat transfer characteristics of impinging jets. Int. J. Heat Mass Transfer. –1965. – V. 8. – pp. 1261-1272.

Все возрастающие требования к защите окружающей среды остро ставят перед предприятиями нефтяной промышленности Республики Казахстан задачу снижения до минимума числа аварий при добыче и транспортировке нефти. В нефтяных месторождениях Узень, Жетыбай, Каражанбас и др. Мангистауской области скопились нефтяные отходы, образовавшиеся в результате аварий, нарушении технологического режима подготовки нефти и технической эксплуатации нефтяных скважин. Аэрофотосъемка показывает наличие большого количества амбаров разлитой нефти, которые загрязняют территорию нефтяных месторождений и оказывают негативное воздействие на окружающую среду. Амбарная нефть является полезным углеводородным сырьем, хотя и претерпела серьезные структурные изменения. Верхний слой затвердел и не поддается излечению. Для сбора конденсированной среды предложен термомеханический способ разжижения амбарной нефти в комплексе с передвижной установкой [1].

Опытно-промышленные испытания установки, проведенные в амбарах месторождения Жалгыз-Тобе (декабрь, 2003 г.), Каражанбас (2004-2005 гг.), Узень (октябрь, 2006 г.) Казахстана показывают, что гидродинамика взаимодействия паровой сверхзвуковой струи с амбарной нефтью подобна развитию неизотермической струи жидкости в конденсированной среде (см. рис. 1). В этой связи ниже рассматривается задача распространение турбулентной неизотермической, круглой струи жидкости в объеме амбарной нефти.

aisa1.tif

Рис. 1. Фотография струи в амбарной нефти

Математическая модель процесса. Круглая струя жидкости с начальной температурой То, скоростью Uo и расходом Go развивается в высоковязкой среде (резервуар амбарной нефти). Скорость истечения – дозвуковая, температура струи высокая (То=373К) и позволяет подогревать высоковязкую среду. Для простаты анализа предполагается, что жидкость струи имеет такую же физико-химическую свойству, как и амбарная нефть, реологические свойства которой, удовлетворяет модель ньютоновской жидкости. Амбарная нефть имеет температуру (Тw=303К), чуть выше ее температуры застывания, что позволяет не учитывать теплоту плавления высоковязкой среды. По мере развития струи в резервуаре ее скорость затухает, тепло, вносимое горячей струей, подогревает высоковязкую среду и вызывает снижение ее вязкости, тем самым вовлекает амбарную нефть в движение. Теплофизические свойства амбарной нефти считаются функциями от температуры и находятся путем лабораторных исследований. Схема течения и область интегрирования приведены на рис. 2. Для простоты анализа считается, что амбарная нефть находится в емкости цилиндрической формы. Дно и боковые стенки емкости являются твердыми поверхностями, а верхняя часть – свободной, граничащей атмосферным воздухом. Круглая струя горячей жидкости подается по центру емкости и ее движения рассматривается в цилиндрической системе координат, ось OZ направляется по оси движение струи, а ось OR – по радиусу круглой струи. Вследствие осевой симметрии задачи относительно оси OZ рассматривается область с размерами asaev1.wmf; asaev2.wmf (см. рис. 2).

aisa2.tif

Рис. 2. Схема струи в емкости

Движение струи считается турбулентным и процесс нагрева амбарной нефти – стационарным. Неизотермическое движение жидкости в амбарной нефти описывается системой уравнении Навье-Стокса, осредненной по Рейнольдсу [2, 3], и, имеющей вид.

asaev3.wmf (1)

asaev4.wmf (2)

asaev5.wmf. (3)

Уравнение переноса тепла с учетом диссипации кинетической энергии движения можно записать в виде:

asaev6.wmf, (4)

где asaev7.wmf +asaev8.wmf – диссипация кинетической энергии движения в тепло.

В уравнениях (1) – (4): z, r – цилиндрические координаты; u, asaev9.wmf – компоненты вектора скорости asaev10.wmf; asaev11.wmf – давления, плотность, температура и теплоемкость жидкости; asaev12.wmf, asaev13.wmf – коэффициент динамической вязкости жидкости, asaev14.wmf – коэффициент турбулентной вихревой вязкости; asaev15.wmf, l – коэффициент теплопроводности жидкости, asaev16.wmf, asaev17.wmf – турбулентный аналог числа Прандтля.

Теплофизические свойства жидкости были найдены по результатам лабораторных исследований амбарной нефти в интервале температуры asaev18.wmf и обобщены в виде эмпирических зависимостей:

asaev19.wmf (кг/м3);

asaev20.wmf (кг/(м·с)); (5)

asaev21.wmf (Вт/(м·град);

Теплоемкость жидкости asaev22.wmf в интервале температуры asaev23.wmf мало меняется и считается постоянной asaev24.wmf=0,23 кдж/(кг·град).

Известно, что в инженерных задачах широко используется asaev25.wmf – модель турбулентности и позволяет получить достаточные надежные данные [3-8]. Например, в работе [6] даны результаты сравнения расчета пристенной турбулентной, плоской струи, вытекающей встречный поток в канале, полученные тремя моделями турбулентности: 1) LES-модели турбулентности; 2) asaev26.wmf – модели турбулентности для низких чисел Рейнольдса Re; 3) модели турбулентных напряжений с использованием граничных условий в виде закона стенки. Расчеты показывают, что данные полученные asaev27.wmf – модели турбулентности для низких чисел Рейнольдса Re, находятся в хорошем количественном согласии с данными LES-модели турбулентности и результатами экспериментальных исследований [7]. В некоторой степени рассматриваемая задача подобна развитию турбулентной струи встречном потоке, что позволяет выбрать для моделирования asaev28.wmf – модель турбулентности. Коэффициент турбулентной динамической вязкости asaev29.wmf находится на основе asaev30.wmf – модели турбулентности, предназначенной для низких чисел Рейнольдса Re, и выражается формулой [8]:

asaev31.wmf, (6)

где asaev32.wmf – скорость диссипации кинетической энергии турбулентности, asaev33.wmf – пристеночная функция, равная [8]: asaev34.wmf.

Дифференциальные уравнения кинетической энергии турбулентности и скорости ее диссипации имеют вид:

asaev35.wmf (7)

asaev36.wmf (8)

где asaev37.wmf – пристеночные функции, равные [8]:

asaev38.wmf, asaev39.wmf. (9)

Константы asaev40.wmf – модели равны asaev41.wmf.

Число Прандтля asaev42.wmf находится при начальной температуре жидкости asaev43.wmf и принимается, равным asaev44.wmf = 0,9. Турбулентный аналог числа Прандтля asaev45.wmf по данным [9,12] можно взять, равным asaev46.wmf= 0,9. Для удобства решения уравнение переноса тепла записывается относительно избыточной температуры asaev47.wmf, где asaev48.wmf – температура нефти в ложе амбара.

Система уравнений (1) – (9) приводится к безразмерным переменным. Координаты z, r делятся на радиус входного сечения струи; компоненты скорости u, v – на максимальную скорость струи на входе; давление р – на максимальное значение динамического напора струи; температура Т – на То; плотность, коэффициенты динамической вязкости, теплопроводности, теплоемкости – на значения этих величин при температуре То.

Граничные условия задачи имеют вид:

при z = 0:

asaev49.wmf; asaev50.wmf; asaev51.wmf=0; asaev52.wmf; asaev53.wmf; asaev54.wmf (10)

при z = 0:

asaev55.wmf; asaev56.wmf; – asaev57.wmf; (11)

при z = L:

asaev58.wmf; asaev59.wmf; – asaev60.wmf;

при r = 0:

asaev61.wmf; asaev62.wmf; (12)

при r = Rw: asaev63.wmf; asaev64.wmf; – asaev65.wmf.

где t – касательное напряжение; Bio – число Био.

Граничные условия на верхней границе (z = 0) для осевой компоненты скорости, кинетической энергии турбулентности и скорости ее диссипации, избыточной температуры во входном сечении струи (10) соответствуют развитому течению круглой турбулентной струи, а вне области втекания струи (11) ставятся условия на свободной поверхности (равенство нулю касательного напряжения [9]). На нижней твердой поверхности (z = L) ставятся условия прилипания и теплообмен с окружающей средой. На левой границе (r = 0) – условия симметричности течения, а на правой границе (r = Rw) – условия прилипания и теплообмен с окружающей средой.

Система уравнений (1) – (9) с граничыми условиями (10) – (12) решается численным методом в переменных функция тока и напряженность вихря [10 – 13]. Функция тока asaev66.wmf вводится, удовлетворяя уравнения неразрывности движения, посредством равенств:

asaev67.wmf, asaev68.wmf, (13)

а напряженность вихря asaev69.wmf – стандартным выражением:

asaev70.wmf. (14)

Подставляя (13) в (14), получим уравнение для функции тока:

asaev71.wmf. (15)

Из системы уравнения движения, исключая давления, стандартным образом можно получить уравнение для напряженности вихря [10-12]:

asaev72.wmf, (16)

где asaev73.wmf – источниковый член определяется также, как в [10].

Граничные условия для уравнения (15), (16) получаются из (10) – (12):

при z = 0: asaev74.wmf; asaev75.wmf asaev76.wmf, asaev77.wmf, asaev78.wmf; asaev79.wmf

при z = 0:

asaev80.wmf; asaev81.wmf; –asaev82.wmf; (17)

при z = L: asaev83.wmf; asaev84.wmf; –asaev85.wmf; asaev86.wmf

при r = 0: asaev87.wmf; asaev88.wmf;

при r = Rw:

asaev89.wmf; asaev90.wmf; –asaev91.wmf; asaev92.wmf. (18)

Как видно из (18), на твердых поверхностях граничные условия для напряженности вихря определяются по формуле Тома [10, 11].

Уравнения для напряженности вихря, переноса тепла, кинетической энергии турбулентности и скорости ее диссипации были аппроксимированы гибридной схемой [10-13]. Разностные аналоги уравнения рассчитывались методом верхней релаксации [11,12].

Для апробации математической модели и численного метода решена задача – расчет турбулентной, плоской струи воздуха набегающей на плоскую нагретую стенку. Расчетные данные по определению местного числа Стэнтона, выражающего теплообмен между стенкой и обтекаемым потоком, были сопоставлены с расчетами Сполдинга [10] и экспериментальными данными Гардона и Акифрата [14] (см. рис. 3).

aisa3.tif

Рис. 3. Верификация численного метода расчета

Результаты сравнения показывают, что при больших числах Рейнольдса (Re=11000, 22000) расчетные данные находятся в удовлетворительном качественном и количественном согласии с опытными данными [14]. При значении числа Рейнольдса (Re=5500) имеется количественное отклонение расчета с экспериментом, хотя качественное согласие зависимости местного числа Стэнтона довольно хорошее (см. рис. 3). Следует отметить также, что закономерности турбулентных струйных течений, как: 1) линейный закон развития границы струи; 2) подобия по числу Рейнольдса динамических и тепловых характеристик осредненного и пульсационного движения хорошо выполняются в расчете.

Обсуждение расчетных данных. Основные режимные параметры течения – число Рейнольдса Re (рассчитанное по средней скорости и диаметру на входном сечении струи); число Bio (выражающее теплообмен жидкости на дне и боковой стенке резервуара); число Biof (выражающее теплообмен жидкости на открытой поверхности); отношения hc/Rc, Rw/Rc, где hc – расстояние между соплом и дном резервуара, а Rc, Rw, – радуисы сопла и резервуара, соответственно; ламинарное Pr и турбулентное Prt числа Прандтля.

На рис. 4 – 6 приведены изолиний напряженности вихря asaev94.wmf, функций тока asaev95.wmf, вектора скорости asaev96.wmf, избыточной температуры asaev97.wmf, кинетической энергии турбулентности k и турбулентной вязкости asaev98.wmf при Re = 11000, Bio = 1, Biof = 5, Prt = 0,9, hc/Rc = 29 и Rw/Rc = 5,6. Теплообмен на свободной поверхности гораздо выше, чем на стенках резервуара, поэтому в расчетах Biof принимается гораздо больше, чем Bio.

Структура течения состоит из трех частей: 1) струйная часть течения, направленная на дно амбара; 2) область возвратных течений в вихре; 3) течение в пограничном слое боковой стенки резервуара, направленное на свободную поверхность резервуара.

По распределениям функций тока asaev99.wmf, вектора скорости asaev100.wmf, избыточной температуры asaev101.wmf можно судить о закономерности развития неизотермической турбулентной струи в амбарной нефти. В отличие от известного линейного закона изменения границы турбулентных струйных течений в затопленном пространстве [2], граница исследуемого струйного течения сильно искривлена и расширена (см. рис. 4,а). Такое поведение границы струи в амбарной нефти можно объяснить сильной зависимости вязкости нефти от температуры. Амбарная нефть находится при температуре окружающей среды (Тw = 303K) и постепенно вовлекается в струйное движение за счет силы трения. Горячая струя (То=373K) нагревает амбарную нефть путем конвективного перемешивания, молярной и молекулярной теплопроводности и диссипацией механической энергии в тепло. Величина избыточной температуры максимальна на оси и снижается к границам струи, а вязкость нефти, наоборот, принимает минимальное значение на оси струи и возрастает вне струйной части течения. Причем в начале расчетной области отношение максимального и минимального значения вязкости достигает по сечению струи asaev102.wmf. Высокая вязкость нефти и ее зависимость от температуры и является причиной резкого расширения границы струи. Кроме того, в систему уравнений входит турбулентная вязкость asaev103.wmf, характеризующая влияние турбулентных напряжений на структуру течения. Картина изолиний asaev104.wmf показывает распределения турбулентной вязкости в расчетной области (см. рис. 4,b). Нетрудно заметить, что существует поле максимальных значений asaev105.wmf, охватывающее, зону взаимодействия струей с возвратным течением и область торможения струи. Большие значения asaev106.wmf и тем самым турбулентных напряжений обусловлены генерацией энергии турбулентности в этих областях.

На процессы переноса существенное влияние оказывает структура течения в вихре, подобная тору, расположенная с центром в правой части расчетной области (см. рис. 5,a). Замкнутые изолиний функции тока asaev107.wmf характеризуют вовлечение в вихревое движение массу жидкости из амбарной нефти. В застойной зоне жидкость слабо вращается, перенося тепло горячей струи в холодную часть, и, подогревая амбарную нефть. Изолиний напряженности вихря asaev108.wmf заметны в струйной части течения и на боковой стенке, а в остальной части их величины незначительны (см. рис. 5,b).

Картину течения более детально раскрывают изолиний вектора скорости asaev109.wmf (см. рис. 6). В струйной части вектор скорости направлен вниз, вблизи ее границ наблюдается вовлечение жидкости в движение. Струя натекает на дно резервуара, что приводит к изменению направления вектора скорости на 90о и иллюстрирует обтекание твердой поверхности. В застойной зоне вихря видно возвратное течение жидкости. Вдоль боковой стенки вектор скорости направлен вертикально вверх и вблизи свободной поверхности изменяет свое направление на обратное и показывает взаимодействие струи с амбарной нефтью.

а b

aisa4.tiff

Рис. 4. Изолинии избыточной температуры (а) и турбулентной вязкости (b)

a b

aisa5.tif

Рис. 5. Изолинии функции тока (а) и завихренности (b)

aisa6.tif

Рис. 6. Картина вектора скорости

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

Аналогичный характер поведения наблюдается в распределениях изолиний скорости диссипации кинетической энергии турбулентности (см. рис. 7,b). Здесь также максимальное значение скорости диссипации достигается в зоне вовлечения струей амбарной нефти, а в остальной части их величины малы. Расчетные данные турбулентных характеристик течения показывают, что молярный перенос существенную роль играет в струйной части и вихревой зоне. В остальной части течения превалируют конвекция и молекулярный механизм переноса импульса и тепла. Коэффициент турбулентной теплопроводности asaev110.wmf, определяемый по asaev111.wmf модели турбулентности, выражающий молярный перенос тепла, изменяется подобным образом, как asaev112.wmf. Эти данные подтверждают подобия процессов переноса тепла и импульса в турбулентных потоках. В расчетах найдено влияние число Рейнольдса Re = 22000, Re = 44000 на процессы переноса в амбарной нефти без изменения других режимных параметров.

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

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

a b

aisa7.tiff

Рис. 7. Изолинии кинетической энергии турбулентности (а) и скорости диссипации кинетической энергии (b)

Заключение. В заключение можно сделать следующие выводы.

1. Изучено развитие турбулентной неизотермической струи в амбарной нефти с использованием RANS модели. Система уравнений RANS модели решается с известным численным методом.

2. Верификация численного метода решения показывает хорошее согласие расчетных данных с экспериментальными данными. Двухпараметрическая asaev113.wmf – модель турбулентности при низких числах Рейнольдса качественно хорошо раскрывает основные закономерности процессов переноса импульса и тепла при взаимодействии неизотермической струи с амбарной нефтью.

3. Показана, что структура течения состоит из трех частей: 1) струйная часть течения; 2) возвратное течение в вихре тороидальной формы; 3) течение в пограничном слое вязкой жидкости вблизи боковой стенки емкости.

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

5. Расчетные данные подтверждают подобия процессов переноса в турбулентном режиме потока и показывают правильность созданной математической модели течения, удовлетворяющей закон подобия турбулентности по числу Рейнольдса.

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


Библиографическая ссылка

Айсаев С.У., Диярова Л.Д., Жапбасбаев У.К., Аринов Е. РАЗВИТИЕ ТУРБУЛЕНТНОЙ НЕИЗОТЕРМИЧЕСКОЙ СТРУИ В ОБЪЕМЕ АМБАРНОЙ НЕФТИ // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 8-1. – С. 104-112;
URL: https://www.applied-research.ru/ru/article/view?id=7050 (дата обращения: 24.06.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074