CC BY
86
Для цитирования: Хурамшина Р.А., Ямалов А.И. Влияние состава сжиженного природного газа и технологических характеристик хранилищ стартового оборудования на время наступления явления ролловера. Вестник Дагестанского государственного технического университета. Технические науки. 2019;46(1):66-78. DOI:10.21822/2073-6185-2019-46-1-66-78
For citation: Khuramshina R.A., Yamalov A.I. The impact of liquefied natural gas composition and technological characteristics of the starting storage of the equipment at the time of occurrence of the phenomenon of rollover. Herald of Daghestan State Technical University. Technical Sciences. 2019; 46(1):66-78.(In Russ.) DOI:10.21822/2073-6185-2019-46-1-66-78
ЭНЕРГЕТИЧЕСКОЕ, МЕТАЛЛУРГИЧЕСКОЕ И ХИМИЧЕСКОЕ МАШИНОСТРОЕНИЕ
УДК 629.7.036
DOI: 10.21822/2073-6185-2019-46-1-66-78
ВЛИЯНИЕ СОСТАВА СЖИЖЕННОГО ПРИРОДНОГО ГАЗА И ТЕХНОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК ХРАНИЛИЩ СТАРТОВОГО ОБОРУДОВАНИЯ НА ВРЕМЯ НАСТУПЛЕНИЯ
ЯВЛЕНИЯ РОЛЛОВЕРА
2 1
Хурамшина Р.А. , Ямалов А.И.
Уфимский государственный нефтяной технический университет,
450062, г.Уфа, ул.Космонавтов 1,Россия,
1 2
e-mail: [email protected], e-mail: [email protected],
Резюме. Цель. Способы предотвращения наступления ролловера в хранилище сжиженного природного газа. Метод.В статье рассматривается современные методы математического моделирования и программного обеспечения при решении задач в нелинейной постановке. Проведен нормативный анализ существующих методов транспортировки СПГ. Результат. С помощью математической программы ANSYS смоделирован процесс ролловера в хранилище сжиженного природного газа. В програмном комплексе ANSYS было произведено 3d моделирование процесса тепло и массообмена страфицированной жидкости в хранилище сжиженного природного газа. При моделировании явления ролловера время до его наступления составило 30 минут, что отличается от экспериментального значения на 1,7 % По результатам проведенных расчетов можно определить, что время до наступления ролловера зависит от концентрации компонентныхв слоев, от заполнения и геометрии резервуара. Эти факторы приводят к увеличению испарения сжиженного природного газа 0,04 %. Вывод. Метод численного эксперимента позволяет определить распределение давления, температуры, плотности, концентрации в исследуемом объекте, не прибегая к реальным экспериментам. Применение программного комплеса ANSYS вычислительной гидрогазодинамики при процессе тепло и массообмена страфицированной жидкости в хранилище сжиженного природного газа необходима при расчете повышенной точности.
Ключевые слова: математическая модель, сжиженный природный газ, программный комплекс ANSYS , газодинамика, процесс тепло и массообмена, страфицированная жидкость, хранилище сжиженного природного газа, ролловер, двухтопливные двигатели
Вестник Дагестанского государственного технического университета. Технические науки. Том 46, №1, 2019 Herald ofDaghestan State Technical University.Technical Sciences. Vol.46, No.1, 2019 _http://vestnik.dgtu.ru/ISSN (Print) 2073-6185 ISSN (On-line) 2542-095Х_
POWER, METALLURGICAL AND CHEMICAL MECHANICAL ENGINEERING
THE IMPACT OF LIQUEFIED NATURAL GAS COMPOSITION AND TECHNOLOGICAL CHARACTERISTICS OF THE STARTING STORAGE OF THE EQUIPMENT
AT THE TIME OF OCCURRENCE OF THE PHENOMENON OF ROLLOVER
• * 2 1
Regina A.Khuramshina , Azat I. Yamalov
Ufa State Petroleum Technological University,
1 Kosmonavtov Str., Ufa 450062, Russia,
1 2
e-mail: [email protected], e-mail: [email protected]
Abstract Objectives Ways to prevent the onset of rollover in the liquefied natural gas storage. Method. The article deals with modern methods of mathematical modeling and software when solving problems in a nonlinear formulation. A regulatory analysis of existing methods of LNG transportation has been carried out. Result. Using the mathematical program ANSYS, the rollover process in the liquefied natural gas storage was modeled. In the ANSYS software package, 3d modeling of the process of heat and mass transfer of the stratified liquid in the liquefied natural gas storage was performed. When simulating a rollover phenomenon, the time to its occurrence was 30 minutes, which differs from the experimental value by 1.7%. Based on the results of the calculations, it can be determined that the time to the onset of rollover depends on the concentration of the component layers, on the filling and geometry of the reservoir. These factors lead to an increase in evaporation of liquefied natural gas of 0.04%. Conclusion. The method of numerical experiment allows to determine the distribution of pressure, temperature, density, concentration in the object under study, without resorting to real experiments. The use of the ANSYS software package of computational fluid dynamics in the process of heat and mass transfer of the stratified liquid in the liquefied natural gas storage is necessary when calculating the increased accuracy.
Keywords: mathematical model, liquefied natural gas, ANSYS program complex, gas dynamics, heat and mass transfer process, stratified liquid, liquefied natural gas storage, rollover, dual-fuel engines
Введение Природный газ является сегодня и останется на обозримую перспективу жизненно важным компонентом в обеспечении глобальных потребностей из-за своих преимуществ перед другими видами ископаемого топлива и в силу постоянной растущей в нем потребности. Большая часть газа доставляется потребителям по трубопроводам в газообразной форме. В некоторых случаях транспортировка сжиженного природного газа в танкерах предпочтительнее, чем другие способы доставки. Использование СПГ потребует существенного обновления судов, работающих в Арктике. Но это в любом случае придется делать, учитывая значительный возраст арктического флота России. СПГ позволит снизить риск разливов судового остаточного топлива, уменьшить выбросы в воздух загрязняющих веществ. Активные действия в этой области позволят усилить позиции российских производителей на быстро развивающемся мировом рынке СПГ топлива, продолжить освоение самых современных технологий мирового судостроения.
Переход на использование СПГ в качестве бункерного топлива позволит решить еще одну важнейшую экологическую проблему - аварийных разливов судового топлива, особенно при наличии ледового покрова. При использовании СПГ в системах хранения на стартовом комплексе существует возможность протекания негативного явления, получившего название «ролловер», возникающего при дозаправке емкости хранилища новой дозой СПГ, физические параметры которой отличаются от имеющейся, с возникновением через некоторое время перемешивания новой и имеющейся в хранилище порций СПГ с резким увеличением давления в газовой подушке. Данный рост давления может привести к разрушению хранилища и разливу компонента, а при возникновении искры и к пожару на стартовом комплексе космодрома.
Время наступления процесса перемешивания СПГ с резким увеличением давления после дозаправки хранилища по проведенным оценкам [4, 6] и экспериментальным данным [7] может составлять от нескольких часов до десятков часов. Определение времени наступления роллове-ра является важным вопросом безопасности при хранении СПГ на стартовых комплексах. Из-за разных технологических систем хранения, в случае малого времени до протекания явления, необходимо применение специальных средств для исключения последствий ролловера, которые должны предусматриваться при проектировании данной системы.
На изменение кондиционности СПГ могут влиять: расслоение СПГ в емкостях хранилища в результате периодического пополнения их, в том числе при обратном сливе из баков потребителя продуктом, имеющим иную плотность и состав по сравнению с остатком СПГ в резервуарах; попадание примесей из газов подготовки к заправке емкостей хранилища, заправочных трубопроводов и баков потребителя; попадание примесей из газов наддува баков потребителя.
Одним из важных факторов опасности СПГ являются его низкая температура кипения, что не позволяет обеспечить его хранение без потерь, а также изменение состава (плотности) и температуры кипения СПГ при хранении за счет испарения более легкого компонента (метана). Данные факторы при пополнении хранилищ, имеющих остатки топлива, новой порцией СПГ, с плотностью и температурой, отличающейся от аналогичных параметров в остатке топлива, могут приводить к образованию стратифицированных макро слоев СПГ в хранилище [3,4]. При последующем протекании тепло- и массообменных процессов в слоях СПГ возможно возникновение режима интенсивного перемешивания (явление «ролловер») с практически мгновенным испарением больших масс СПГ и резким повышением давления в хранилище. Все это способствует и разрушению хранилища и интенсификации неблагоприятных технологических процессов [5]. Схематическая картина механизмов ролловера указанного явления кратко проиллюстрирована на рис. 1.
Рис. 1. Перемешивание слоев при закачке продукта в резервуар с процессами тепло- и
массообмена
Fig. 1. Mixing of layers during the injection of the product into the tank with the processes of heat and
mass transfer
Заметим, что физическая модель процесса отличается тем, что при загрузке резервуара свойства продукта, уже хранящегося в изотермическом резервуаре (ИР), отличаются от свойств СПГ, закачиваемого в резервуар. В ИР, оборудованных системой нижнего налива, СПГ, имеющий большую плотность, образует нижний слой. СПГ, имеющий меньшую плотность, формирует верхний слой, который располагается над поверхностью нижнего, не перемешиваясь с ним. Перегреваясь относительно температуры насыщения за счет внешнего теплопритока через днище и стенки, продукт в нижнем слое при этом не имеет возможности компенсировать теп-лоприток за счет испарения, т.к. поверхность теплообмена закрыта нижним слоем. Однако имеют место процессы тепло - и массообмена двух слоев между собой и верхнего слоя с парогазовой фазой. В конечном итоге, в определенный момент времени после загрузки и образования стратификации, за счет описанных выше процессов плотности слоев выравниваются с последующей резкой интенсификацией процессов испарения продукта (суть явления ролловер). Математическое моделирование явления ролловера целесообразно проводить на базе уравнений математической физики, характеризующих взаимосвязанные нелинейные процессы пере-
носа теплоты, массы и импульса. В частности, моделирование тепловой части задачи будем проводить с использованием уравнения энергии в энтальпийной форме.
Постановка задачи. Имеется сферический резервуар, объемом 100000 м . Резервуар наполовину заполнен сжиженным природным газом, который представляет собой многокомпонентное рабочее тело, находящееся в жидком агрегатном состоянии. После закачки в резервуар дополнительного объема СПГ он становится заполнен полностью. При этом в начальный момент времени можно выделить два слоя СПГ: нижний слой и верхний слой. Рабочее тело в этих слоях отличается компонентным составом, а также давлением и температурой. Параметры рабочего тела в верхнем и нижнем слое представлены в табл. 1.
Таблица 1. Параметры компонентов рабочего тела Table 1. Parameters of the working body components
Компонента (растворенного вещества) Component (solute) Количество молей вещества Moles of substance
Нижний слой Bottom layer Верхний слой Top layer
Метан Methane 0,6226 0,6362
Этан Ethane 0,2185 0,2416
Пропан+примеси Propane + impurities 0,1589 0,1222
Начальные свойства СПГ LNG Initial Properties
T, К 118,997 114,355
P, МПа 0,15 0,131
p, кг/м3 541,0316 536,9516
Требуется смоделировать явление ролловера, происходящее в резервуаре в течение 30
минут.
Методы исследования. Физические свойства всех компонент рабочего тела задавались в соответствии с данными работы [6]. Принималось, что все теплофизические свойства компо-нентпостоянны, что допустимо при условии небольших перепадов температуры внутри расчетной области. Физические свойства всех компонент приведены в табл. 2.
Таблица 2. Физические свойства компонентов рабочего тела Table 2. Physical properties of the components of the working fluid
Химическая компонента Chemical component Плотность, кг/м3 Density, kg / m3 Удельная теплоемкость, Дж/ktK Specific heat, j / kg'K Удельная теплопроводность, Вт/м/K Heat conductivity, W / m-K Вязкость, кг/ м/с Viscosity, kg / m-s Мольная масса, кг/моль Molar mass, kg / mol
CH4 415,2 5477 0,6 2,1 16,043
C2H6 561,0 2990 0,6 87,8 30,070
CзH8+примеси +impurities 493,1 2230 0,6 0,01 44,090
При определении осреднённых свойств гомогенного потока в контрольном объёме V, ограниченном замкнутой поверхностью А, система уравнений механики сплошной среды записывалась в интегральной декартовой форме в виде:
-л
— ^ШёУ + -С\М = (1)
^ V V
где вектора Ж, Г и G определяются следующим образом
W =
P pv 0
pu pvu + pi
< pv F = • Pv+pj G = • Ty.
pw pw + pk
PE. pE + pv Tjvj + q _
(2)
Вектор Н содержит источниковые члены.
где р, V, Е, р- плотность, скорость, полная энергия на единицу массы и давление среды в контрольном объеме; т - тензор вязких напряжений; q - плотность теплового потока.
Полная энергия Е связана с полной энтальпией Н следующим соотношением
Е = Н - р / р, где Н = и + 2 /2.
Для каждого из компонентов многокомпонентного рабочего тела решалось от-отдельное уравнение переноса массы в виде
(3)
(4)
д
—{pC, )+V-
дтк }
(5)
где С— локальная массовая концентрация /-го компонента;
g¡ - диффузионный поток /-го компонента;
а — скорость образования /-го компонента в химических реакциях.
В уравнении (5) члены слева направо учитывают: нестационарность процесса переноса массы вещества, перенос массы вещества за счет конвекции, перенос массы вещества диффузией, источник массы вещества, обусловленный наличием химических реакций. Последний член в правой части уравнения (5) принимался равным нулю.
Диффузионный поток массы g¡ вычислялся по закону Фика. Для моделирования диффузии массы выражение для диффузионного потока может быть записано в следующем виде (для турбулентных потоков)
V Т
УУ - А.,-, (6)
g, =
T ,
T
где - коэффициент ламинарной диффузии массы для /-го компонента;
турбулентная вязкость; турбулентное число Шмидта; Уг - массовая концентрация /-го компонента; Ат ^ - коэффициент термической диффузии (коэффициент Соре); Т- статическая температура.
В данной математической модели коэффициент В/т, который представляет собой коэффициент ламинарной диффузии массы для каждого /-го компонента, был заменен на один бинарный коэффициент Б.
Такой подход гарантирует то, что сумма диффузионных потоков будет равна нулю. Значение коэффициента D вычислялось в предположении, что число Льюиса равно единице.
А = —, (7)
Рср
где к- удельная теплопроводность смеси;
р - плотность смеси;
ср - удельная теплоёмкость смеси.
Турбулентное число Шмидта записывается следующим образом:
V
Sct —
pDt
(8)
где Б- коэффициент турбулентной диффузии массы;
р - плотность смеси.
В рамках данной работы рассматриваются течения с большими числами Рейнольдса. По этой причине с целью моделирования турбулентных эффектов вышеописанная математическая модель была дополнена гибридной моделью турбулентности ТгапвШопЗБТ, которая показывает достаточно достоверные результаты для большинства задач, связанных с естественной конвекцией.
Модель турбулентности ТгапвШопЗБТ представляет собой гибридную модель, которая была получена путем модификации модели турбулентности ББТк- т. Известно [6], что модель турбулентности ББТк- т основана на осреднении по Рейнольдсу и относится к классу моделей вихревой вязкости. Она решает два уравнения переноса скалярной величины: одно для турбулентной кинетической энергии к и одно для частоты турбулентных пульсаций т.
Модель турбулентности ТгапвШопББТ была получена путем дополнения модели ББТк-т, ещё двумя уравнениями: одно уравнение для перемещаемости и одно для критерия начала перехода в членах числа Рейнольдса по толщине потери импульса. Эмпирические коэффициенты модели турбулентности были подобраны таким образом, чтобы обеспечить достоверные результаты во всем диапазоне чисел Рейнольдса, особенно при моделировании обтекания тел потоками газа с низкой начальной турбулизацией.
Уравнение переноса для перемещаемости у имеет следующий вид
djpy) | d(pUjY)
—P-Е,+P,-£..„+
d
dt x 7 71 72 72 dx]
72
"72
Л
V
V + —
V °7 J
d7
dx,.
(9)
Переходные источниковые члены в уравнении (9) имеют следующий вид
Р = СА Е^рБ [уЕоте( р3; (10)
ЕУ1 = САРГ#, (11)
где ^-величина сдвиговых напряжений;
^«^-эмпирическая функция, отвечающая за длину переходной области течения:Са1= 2,0; Се1= 1,0- эмпирические константы.
Уравнение переноса для переходного числа Рейнольдса по толщине потери импульса имеет следующий вид
d(^pReet} d^pUjReet} Q
I — Pq+ +
dt
dx,.
dx,.
(v + Vt )
dRe
dx,.
et
(12)
3 3
Присутствующие в уравнении (22) эмпирические константы имеют следующие значения
са = 0.03, = 2.0. (13)
Выбор эмпирических функций для переменных Яеви Р^л и Яеос подробно описывается в работе [7]. Известно, для того чтобы решить дискретные аналоги уравнений Навье-Стокса и всех дополнительных моделей необходимо задать граничные условия, которые определят поведение всех уравнений на границах расчётной области.
В рамках данной работы сферическая расчётная область ограничивалась поверхностью, представляющей собой твердую стенку для которой задавались следующие параметры: для уравнений количества движения - условие прилипания (принималось, что на поверхности нормальная и касательная компоненты вектора скорости потока равны нулю (их = 0, иу = 0, и2 = 0)); для уравнения энергии - условие адиабатичности (данное условие является граничным условием второго рода; предполагается, что тепловые потоки в стенку равны нулю, теплообмен
между пространством резервуара и окружающей средой отсутствует; для уравнений переноса химических компонентов концентрации всех компонентов на поверхности определялись путем интерполяции из ячеек вверх по потоку.
Обсуждение результатов. Для решения поставленной задачи была построена неструктурированная расчетная сетка размерностью 830000 ячеек, которая показана на рис. 2. Для достоверного моделирования градиентов физических величин в пограничном слое на поверхно-
а б
а - общий вид; б - призмослой укрупнено
Рис. 2. Расчетная сетка (сечение по плоскости OXY) a - general view; b - prism layer enlarged Fig. 2. Calculated grid (section along the OXY plane) Выбор начального приближения может оказать значительное влияние на устойчивость расчета, а также на скорость его сходимости. Неправильный выбор может вызвать неустойчивость расчёта, особенно в начальный момент итерационного процесса. Это в свою очередь будет требовать проведения итерационного процесса с меньшими значениями критерия Куранта -Фридрихса - Леви. Что выразится, в конечном счете, в большем потребном числе итераций и, следовательно, в увеличении необходимых расчётных ресурсов. В рамках данной работы предполагалось, что заполнение резервуара происходит много быстрее, чем время протекания рол-ловера. В связи с этим в качестве начального приближения принималось такое состояние, в котором параметры рабочего тела внутри резервуара соответствуют параметрам в табл. 1. То есть имеется нижний слой и верхний слой. Концентрация всех компонентов в нижнем и верхнем слое в начальный момент времени показана на рис. 3.
а б в
а - концентрация CH4. б - концентрация C2H6; в - концентрация C3H8 + примеси Рис. 3 Мольные концентрации всех химических компонент в начальный момент времени
a - concentration of CH4; b — C2H6 concentration; в - concentration of C3H8 + impurity Fig. 3 Mole concentrations of all chemical components at the initial moment of time
В ходе расчетов были получены все нестационарные распределения компонент рабочего тела внутри резервуара, а также нестационарные поля всех физических величин
На рис. 4 - 6 приведено поле векторов скорости в плоскости сечения OXY в моменты времени. Приведены поля мольных концентраций всех компонент в плоскости сечения OXY в моменты времени
Рис. 4. Поле векторов скорости в момент времени 180 с Fig. 4.Field of velocity vectors at time instant 180 s
Рис.5. Поле векторов скорости в момент времени 720 с Fig.5. The field of velocity vectors at the moment of time 720 s
а б в
а - концентрация СН4; б - концентрация С2Н6;в - концентрация СзИ8 включая примеси Рис. 6. Мольные концентрации всех химических компонент a - concentration [CH] 4; b - concentration C_2 H_6; c - concentration C3H8 including impurities
Fig. 6. The molar concentrations of all chemical components
Результаты моделирования распределение скорости, температуры, плотности и мольных концентраций в верхнем и нижнем слое в зависимости от времени представлены на рис.7 - 12.
U, м/с
0,25 0,2 0,15 0,1 0,05 0
-0,05
\
A •A
ft V —■
¿У ч><
0 20 00 40 00 60 00 80 00 100
■Верхний слой ■Нижний слой
t, с
Рис.7. Распределение скорости в верхнем и нижнем слое в зависимости от времени
Рис.8. Распределение температуры в верхнем и нижнем слое в зависимости от времени Fig.8.Temperature distribution in the upper and lower layers depending on time
472
■a
I <2 471,5 X
H to u n В
471
Верхний слой Нижний слой
500
1000
1500
2000
t, с
0
Рис.9. Распределение плотности в верхнем и нижнем слое в зависимости от времени Fig.9. Density distribution in the upper and lower layer depending on time
Y, GH
0,638 0,636 0,634 0,632 0,63 0,628 0,626 0,624 0,622 0,62
500
1000
1500
2000
Верхний слой Нижний слой
t, с
0
Рис.10. Распределение мольной концентрации СЩв верхнем и нижнем слое в зависимо сти от времени
Рис.11. Распределение мольной концентрации С2Н6в верхнем и нижнем слое в зависимости от времени Fig.11. The distribution of the molar concentration of C2H6 in the upper and lower layers depending on time
G3H8
SOO
1OOO
1SOO
2OOO
Верхний слой
Нижний слой
t, с
O
Рис.12. Распределение мольной концентрации C3H в верхнем и нижнем слое в зависимости от времени Fig.12. The distribution of the molar concentration of C3H8 in the upper and lower layer depending on time
По результатам проведенных расчетов можно определить, что время до наступления ролловера зависит от состава сжиженного природного газа, от заполнения и геометрии резервуара. Проведено численное моделирование процессов перемешивания СПГ в резервуаре объемом 100000 м . Дискретизованные по методу конечного объема уравнения газодинамики решались на трехмерной неструктурированной расчетной сетке.
При численном моделировании явления ролловера время до его наступления составило 30 минут, что отличается от экспериментального значения на 1,3 %. Остальные расчетные характеристики по сравнению с экспериментальными данными имеют следующие погрешности: температура верхнего слоя - 1 %; температура нижнего слоя - 0,4 %; плотность верхнего слоя -0,7 %; плотность нижнего слоя - 0,7 %; разница плотностей слоев - 2,7 %. По результатам выполненных расчетов выявлено, что время до наступления явления зависит от состава СПГ, а так же степени заполнения, характеристик изоляции и геометрии хранилища.
Вывод. При использовании СПГ в системах хранения на стартовом комплексе существует возможность протекания негативного явления, получившего название «ролловер», возникающего при дозаправке емкости хранилища новой дозой СПГ, физические параметры которой отличаются от имеющейся, с возникновением через некоторое время перемешивания новой и имеющейся в хранилище порций СПГ с резким увеличением давления в газовой подушке [4, 25]. Можно сделать следующий вывод: метод численного эксперимента позволяет определить распределение давления, температуры, плотности, концентрации в исследуемом объекте, не прибегая к реальным экспериментам.
Меры по предотвращению ролловера: партии СПГ с различной плотностью, хранить отдельно; загрузка резервуара специальным оборудованием, таким, как насадки способствующие смешиванию; СПГ при загрузке (используется для береговых резервуаров); следует избегать продолжительных остановок подачи СПГ при загрузке резервуара; постоянно контролировать уровень испарения СПГ.
Как было отмечено выше, причина явления переворачивания слоев - недостаточное смешивание поступающей в резервуар партии продукта с уже находящейся в нем жидкостью. Однако стратификация с дальнейшим вскипанием может также произойти при выдерживании СПГ в течение длительного времени в резервуаре без рециркуляции или при содержании азота в составе СПГ более 1 % (такая пороговая концентрация установлена по данным международных исследований явления rollover).
Для предотвращения стратификации резервуар оснащается системой рециркуляции с использованием насосов, а налив продукта должен быть организован предпочтительно таким образом, чтобы более тяжелый продукт образовывал верхний слой и смешение фаз, происходило за счет естественной гравитации. Поэтому конструкция резервуара должна предусматривать системы верхнего и нижнего налива, а также средства контроля плотности и температуры на разных уровнях. Кроме того, как правило, предусматривается многоуровневая защита резервуара от превышения внутреннего избыточного давления газовой фазы, обеспечивающая сброс газа в атмосферу или на факел [8]. Чтобы предсказать параметры, при которых может произойти явление ролловера, предложено провести математическое моделирование процесса. Математическое моделирование явления ролловера предполагает использование уравнений математической физики, характеризующих процессы переноса теплоты, массы и импульса в рамках термодинамики необратимых процессов.
В програмном комплексе ANSYS было произведено 3d моделирование процесса тепло и массообмена страфицированной жидкости в хранилище сжиженного природного газа. При моделировании явления ролловера время до его наступления составило 30 минут, что отличается от экспериментального значения на 1,7 % По результатам проведенных расчетов можно определить, что время до наступления ролловера зависит от концентрации компонентныхв слоев, от заполнения и геометрии резервуара. Эти факторы приводят к увеличению испарения сжиженного природного газа 0,04 %.
Метод численного эксперимента позволяет определить распределение давления, температуры, плотности, концентрации в исследуемом объекте, не прибегая к реальным экспериментам. Применение программного комплеса ANSYS вычислительной гидрогазодинамики при процессе тепло и массообмена страфицированной жидкости в хранилище сжиженного природного газа необходимо при расчете повышенной точности.
Библиографический список:
1. Перспективы использования СПГ в Арктическом регионе РФ по оценкам WWF [Электронный ресурс]. -Режим доступа:http://proarctic.ru/04/07/2017/ecology/27146(дата обращения: 25.04.2018).
2. Власов А.В. Придет ли природный газ на смену мазуту? - Новые технологии. - Выпуск 4 (10). - 2010. - с. 67-70
3. Королев Н.С. К построению математической модели явления ролловер в хранилище СПГ // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон.журн. 2012. № 3.
4. Королев Н.С. Математическое моделирование явления ролловер в системах хранения сжиженного природного газа // Актуальные проблемы российской космонавтики: материалы XXXVII академических чтений по космонавтике. М.: КомиссияРАН, 2013. C. 383-384.
5. May E.F. Fluid Science for improved LNG Production and shipping. Presentation of the University Of Western Australia, 2010.
6. Варгафтик Н.Б., Справочник по теплофизическим свойствам газов и жидкостей // М. Наука. 1972.
7. Menter F. R., Langtry R. B., Likki S. R., Suzen Y. B., Huang P. G., Volker S.«A Correlation-Based Transition Model Using Local Variables: Part I — Model Formulation».(ASME-GT2004-53452),2004.pplications //AIAA Journal.August 1994. №32 (8).pp. 1598-1605.
8. Васильев, Г. Г. Особенности обеспечения безопасной эксплуатации крупногабаритных изотермических резервуаров для хранения сжиженного природного газа / Г. Г. Васильев, С. Г. Иванцова, А. И. Рахманин // Газовая промышленность. - Выпуск: 11 11. - Москва: РГУ нефти и газа им. И.М. Губкина, 2013 - С. 5761.
9. Рахимов В.О. Особенности термодинамических процессов при хранении сжиженного природного газа: автореф. дисс. на соиск. уч. степ. кандидата техн. наук: 25.00.19/ Рахимов Вадим Олегович -М., -2013, -24 с.
10. Шаммазов А.М. Производство, хранение и транспорт сжиженного природного газа/ А.М. Шаммазов, Р.К. Терегулов, Б.Н.Мастобаев, Г.Е. Коробков. -СПб.: Недра, 2007. -152 с.
11. Liquefied Natural Gas Worldwide [Электронный ресурс].- Режим доступа:http://www.energy.ca.gov/lng/intemational.html (дата обращения: 25.04.2018)
12. Liquefied Natural Gas Worldwide [Электронный ресурс].- Режим доступа: http://www.energy.ca.gov/lng/international.html (дата обращения: 25.04.2018).
13. De Wit, Johannes. Liquefied Gas Storage Tank Design, Construction Addressed/ Johannes de Wit// Oil & Gas Journal. -1987. - 6 July
14. Королев Н.С. Анализ возникновения явления «ролловер» в системах хранениясжиженного природного газа // Актуальные проблемы российской космонавтики: Материалы XXXVI академических чтений по космонавтике. М.:Комиссия РАН, 2012. C. 353-354.
15. MenterF. R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering AApplications //AIAA Journal.August 1994.№ 32(8).pp. 1598-1605.
16. Menter F. R., Langtry R. B., Likki S. R., Suzen Y. B., Huang P. G., Volker S.«A Correlation-Based Transition Model Using Local Variables: Part I — Model Formulation».(ASME-GT2004-53452),2004.pplications //AIAA Journal.August 1994. №32 (8).pp. 1598-1605.
17. Kharlamov S.N. Mathematical Modelling of Thermo- and Hydrodynamical Processes in Pipelines. Rome, Italy: Publ. House "Ionta", 2010. 263p.
18. Issledovanie voprosov proizvodstva i ispol'zovaniia SPG v kachestve al'ternativnogo topliva v razlichnykh ot-rasliakh narodnogo khoziaistva: otchet [Study on the production and use of LNG as an alternative fuel in various sectors of economy. Report]. Moscow, MEI, «Krion», 1994. 125 p. Non published
19. Kuzin A.I., Lozin S.N., Lekhov P.A., Semenov A.I., Mamin V.V. Issledovaniia GKNPTs im. M. V. Khrunicheva po obosnovaniiu mnogorazovoi trebuemoi razmernosti raketno-kosmicheskoi sistemy [Research of the Khrunichev State space scientific-production center on the substantiation of reusable required dimensionality of space-rocket system]. Aviakosmicheskaia tekhnika i tekhnologiia, 2010, no. 1, pp. 3-12.
20. Kuzin A.I., Rachuk V.S., Koroteev A.S., et al. Obosnovanie vybora komponentov raketnogo topliva dlia dvigatel'nykh ustanovok pervoi stupeni mnogorazovoi raketnokosmicheskoi sistemy [Justification of the choice of components of rocket fuel for propulsion of the first stage of reusable space-rocket systems]. Aviakosmicheskaia tekhnika i tekhnologiia , 2010, no. 1, pp. 19-55.
21. Korolev N.S. Matematicheskoe modelirovanie iavleniia rollover v sistemakh khraneniia szhizhennogo prirodnogo gaza [Mathematical modeling of phenomena of rollover in the systems of storage of liquefied natural gas]. Ak-tual'nye problemy rossiiskoi kosmonavtiki: materialy 37 akademicheskikh chtenii po kosmonavtike [Actual problems of Russian astronautics: Materials 37th academic readings on Astronautics]. Moscow, RAS Commission Publ., 2013, pp. 383-384.
22. Uznanski D., Château E., Gorieu O., Legrandais J.-P. Investigations into the aging extent of the 9% nickel 500 mlng storage tank of the nantes cryogenic testing station. Paris, France, Gaz de France, 2003.
23. Chugunkov V.V. Teploperedacha pri slozhnom teploobmene na poverkhnostiakh konstruktsii [Heat transfer in complex heat transfer on surfaces of structures]. Moscow, Bauman MSTU Publ., 2001. 28 p.
24. The European fuels conference 11th annual meeting // HYDROCARBON WORLD. Paris, France: World Refining Association, 2010. Vol. 4, iss. 2. C. 35-38.
25. Бармин И.В., Кунис И.Д. Сжиженный природный газ вчера, сегодня, завтра / Под ред. А.М. Архарова. М.: Изд-во МГТУ им. Н.Э. Баумана, 2009. 256 с.
References:
1. Perspektivy ispol'zovaniya SPG v Arkticheskom regione RF po otsenkam WWF [Elektron-nyy resurs]. - Rezhim dostupa:http://proarctic.ru/04/07/2017/ecology/27146(data obrashcheniya: 25.04.2018). [Prospects of LNG use in the Arctic region of the Russian Federation according to WWF [Electronic resource]. - Mode of access:http://pro arctic.ru/04/07/2017/ecology/27146(date accessed: 25.04.2018). (In Russ)]
2. Vlasov A.V. Pridet li prirodnyy gaz na smenu mazutu? - Novyyetekhnologii. - Vypusk 4 (10). - 2010. - s. 67-70 [Vlasov A. V. will the natural gas to replace fuel oil? - New technology. - Issue 4 (10). - 2010. - p. 67-70(In Russ)]
3. Korolev N.S. K postroyeniyu matematicheskoy modeli yavleniya rollover v khranilishche SPG // Nauka i obra-zovaniye. MGTU im. N.E. Baumana. Elektron.zhurn. 2012. № 3. [Queens N. S. To the construction of a mathematical model of the rollover phenomenon in LNG storage / / Science and education. MGTU im. N. Uh. Bauman. Electron.journal. 2012. No. 3. (In Russ)]
4. Korolev N.S. Matematicheskoye modelirovaniye yavleniya rollover v sistemakh khraneniya szhi-zhennogo prirodnogo gaza // Aktual'nyye problemy rossiyskoy kosmonavtiki: materialy XXXVII akademi-cheskikh chteniy po kosmonavtike. M.: KomissiyaRAN, 2013. C. 383-384. [Queens N. S. Mathematical modeling of the rollover phenomenon in liquefied natural gas storage systems // Actual problems of Russian cosmonautics: materials of XXXVII academic readings on cosmonautics. M.: Commission of RAS, 2013. 383-384. (In Russ)]
5. May E.F. Fluid Science for improved LNG Production and shipping. Presentation of the University Of Western Australia, 2010.
6. Vargaftik N.B., Spravochnik po teplofizicheskim svoystvam gazov i zhidkostey // M. Nauka. 1972. [Vargaftik N. B. Reference book on thermophysical properties of gases and liquids // M. Nauka. 1972. (In Russ)]
7. Menter F. R., Langtry R. B., Likki S. R., Suzen Y. B., Huang P. G., Volker S.«A Correlation-Based Transition Model Using Local Variables: Part I — Model Formulation».(ASME-GT2004-53452),2004.pplications //AIAA Journal.August 1994. №32 (8).pp. 1598-1605.
8. Vasil'yev, G. G. Osobennosti obespecheniya bezopasnoy ekspluatatsii krupnogabaritnykh izo-termicheskikh rezervuarov dlya khraneniya szhizhennogo prirodnogo gaza / G. G. Vasil'yev, S. G. Ivantsova, A. I. Rakhmanin // Gazovaya promyshlennost'. - Vypusk: 11 11. - Moskva: RGU nefti i gaza im. I.M. Gubkina, 2013 - S. 57-61. [Vasil'ev, G. G. Features ensure the safe operation of large insulated storage tanks for liquefied natural gas / G. G. Vasilyev, S. G. Ivantsov, A. I. Rachmaninov // Gas industry. - Issue: 11 11. - Moscow: Russian state University of oil and gas. I. M. Gubkina, 2013-P. 57-61. (In Russ)]
9. Rakhimov V.O. Osobennosti termodinamicheskikh protsessov pri khranenii szhizhennogo prirodnogo gaza: avtoref. diss. na soisk. uch. step. kandidata tekhn. nauk: 25.00.19/ Rakhimov Vadim Olegovich -M., -2013, -24 s. [Rakhimov V. O. features of thermodynamic processes during storage of liquefied natural gas: autoref. dis. on competition of a scientific degree. academic step. the candidate tehn. Sciences: 25.00.19 / Rakhimov Vadim Ole-govich-M., -2013, -24 p. (In Russ)]
10. Shammazov A.M. Proizvodstvo, khraneniye i transport szhizhennogo prirodnogo gaza/ A.M. Shammazov, R.K. Teregulov, B.N.Mastobayev, G.Ye. Korobkov. -SPb.: Nedra, 2007. -152 s. [Shammazov A. M. Production, storage and transport of liquefied natural gas/ by A. M. Shammazov, R. K. Teregulov, B. N. Mastobayev, G. E. Box. -SPb.: Nedra, 2007. -152 p.vvv(In Russ)]
11. Liquefied Natural Gas Worldwide [Электронный ресурс].- Режим доступа:http://www.energy.ca.gov/lng/intemational.html (дата обращения: 25.04.2018)
12. Liquefied Natural Gas Worldwide [Электронный ресурс].- Режим доступа: http://www.energy.ca.gov/lng/international.html (дата обращения: 25.04.2018).
13. De Wit, Johannes. Liquefied Gas Storage Tank Design, Construction Addressed/ Johannes de Wit// Oil & Gas Journal. -1987. - 6 July
14. Korolev N.S. Analiz vozniknoveniya yavleniya «rollover» v sistemakh khraneniyaszhizhennogo prirodnogo gaza // Aktual'nyye problemy rossiyskoy kosmonavtiki: Materialy XXXVI akademicheskikh chteniy po kosmonavtike. M.:Komissiya RAN, 2012. C. 353-354. [Korolev N. S. Analysis of the occurrence of the "rollover" phenomenon in liquefied natural gas storage systems // Actual problems of Russian cosmonautics: materials of XXXVI academic readings on cosmonautics. M.: Commission of RAS, 2012. C. 353-354(In Russ)]
15. MenterF. R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering AApplications //AIAA Journal.August 1994.№ 32(8).pp. 1598-1605.
16. Menter F. R., Langtry R. B., Likki S. R., Suzen Y. B., Huang P. G., Volker S.«A Correlation-Based Transition Model Using Local Variables: Part I — Model Formulation».(ASME-GT2004-53452),2004.pplications //AIAA Journal.August 1994. №32 (8).pp. 1598-1605.
17. Kharlamov S.N. Mathematical Modelling of Thermo- and Hydrodynamical Processes in Pipelines. Rome, Italy: Publ. House "Ionta", 2010. 263p.
18. Issledovanie voprosov proizvodstva i ispol'zovaniia SPG v kachestve al'ternativnogo topliva v razlichnykh ot-rasliakh narodnogo khoziaistva: otchet [Study on the production and use of LNG as an alternative fuel in various sectors of economy. Report]. Moscow, MEI, «Krion», 1994. 125 p. Non published (In Russ)]
19. Kuzin A.I., Lozin S.N., Lekhov P.A., Semenov A.I., Mamin V.V. Issledovaniia GKNPTs im. M. V. Khrunicheva po obosnovaniiu mnogorazovoi trebuemoi razmernosti raketno-kosmicheskoi sistemy [Research of the Khrunichev State space scientific-production center on the substantiation of reusable required dimensionality of space-rocket system]. Aviakosmicheskaia tekhnika i tekhnologiia , 2010, no. 1, pp. 3-12.
20. Kuzin A.I., Rachuk V.S., Koroteev A.S., et al. Obosnovanie vybora komponentov raketnogo topliva dlia dvigatel'nykh ustanovok pervoi stupeni mnogorazovoi raketnokosmicheskoi sistemy [Justification of the choice of components of rocket fuel for propulsion of the first stage of reusable space-rocket systems]. Aviakosmicheskaia tekhnika i tekhnologiia , 2010, no. 1, pp. 19-55. (In Russ)]
21. Korolev N.S. Matematicheskoe modelirovanie iavleniia rollover v sistemakh khraneniia szhizhennogo prirodnogo gaza [Mathematical modeling of phenomena of rollover in the systems of storage of liquefied natural gas]. Ak-tual'nye problemy rossiiskoi kosmonavtiki: materialy 37 akademicheskikh chtenii po kosmonavtike [Actual problems of Russian astronautics: Materials 37th academic readings on Astronautics]. Moscow, RAS Commission Publ., 2013, pp. 383-384. (In Russ)]
22. Uznanski D., Château E., Gorieu O., Legrandais J.-P. Investigations into the aging extent of the 9% nickel 500 mlng storage tank of the nantes cryogenic testing station. Paris, France, Gaz de France, 2003.
23. Chugunkov V.V. Teploperedacha pri slozhnom teploobmene na poverkhnostiakh konstruktsii [Heat transfer in complex heat transfer on surfaces of structures]. Moscow, Bauman Chugunkov V.V. Teploperedacha pri slozhnom teploobmene na poverkhnostiakh konstruktsii [Heat transfer in complex heat transfer on surfaces of structures]. Moscow, Bauman MSTU Publ., 2001. 28 p. (In Russ)]
24. The European fuels conference 11th annual meeting // HYDROCARBON WORLD. Paris, France: World Refining Association, 2010. Vol. 4, iss. 2. C. 35-38.
25. Barmin I.V., Kunis I.D. Szhizhennyy prirodnyy gaz vchera, segodnya, zavtra / Pod red. A.M. Arkharova. M.: Izd-vo MGTU im. N.E. Baumana, 2009. 256 s. [Barmin I. V., Kunis I. D. Liquefied natural gas yesterday, today, tomorrow / ed. by A. M. Arkharov. M.: Izd-vo MGTU im. N. Eh. Bauman, 2009. 256 p. (In Russ)]
Сведения об авторах.
Хурамшина Регина Азатовна - магистрант. Ямалов Азат Ильфатович - инженер второй категории Information about the authors. Regina A. Khuramshina - Master Student. Azat I.Yamalov- Second category Engineer Конфликт интересов
Авторы заявляют об отсутствии конфликта интересов. Поступила в редакцию 23.12.2018. Принята в печать 23. 01.2019.
Conflict of interest.
The authors declare no conflict of interest.
Received 23.12.2018.
Accepted for publication 23. 01.2019.
