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

УДК 622.691

ОСОБЕННОСТИ ПОСТАНОВКИ ЗАДАЧИ РАСЧЕТА ПРОЦЕССОВ ДЛЯ

УНИВЕРСАЛЬНОГО ГАЗОВОГО БАЛЛОНА КАК НЕКЛАССИЧЕСКОЙ ДВУХФАЗНОЙ СИСТЕМЫ И ОЦЕНКА ЕГО ТЕМПЕРАТУРНОГО СОСТОЯНИЯ

© 2006 Т. В. Ашихмина, А. И. Довгялло Самарский государственный аэрокосмический университет

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

Особенности предлагаемого по патенту [1] баллона заключаются в его универсальности по заправке. Баллон, содержащий внутренний теплоизолированный сосуд и систему экранов в газовой полости (рис. 1), может заправляться как заливкой (во внутренний сосуд) криогенной жидкости, так и компри-мированием (заполняется весь объем баллона). Такое устройство баллона и ожидаемые особенности эксплуатации требуют изучения процессов, происходящих при его заправке, хранении продукта в безрасходном состоянии и отборе продукта.

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

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

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

Рис. 1. Конструктивная схема универсального газового баллона:

1 - вентиль; 2 - клапан предохранительный;

3 - баллон; 4 - экраны; 5 - теплоизоляция;

6 - термос; 7 - тарелка прижимная; 8 - фильтр; 9 - вставка амортизационная; 10 - трубка газозаборная; 11-сильфон; 12 - змеевик;

13 - стакан; 14 - заправочный штуцер

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

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

Жидкость начинает испаряться, и испарившаяся доля криопродукта перераспределяется по слоям газовой полости. Давление в начальный момент может даже снизиться, что будет способствовать дальнейшему испарению жидкости. Однако через какое-то время давление может возрасти и кипение прекратится. Особенности начального этапа и последующего поведения криопродукта будут обусловлены соотношением объемов жидкой и газовой фаз. Так, если удельный объем двухфазной смеси (при изохорном процессе) больше V 1д, то жидкость будет непрерывно испаряться вплоть до полного перехода криопродукта в состояние сухого насыщенного пара.

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

Тогда для такой термодинамической системы (всего баллона) уравнение Клапей-

dP г

рона-Клаузиуса в виде dT = Х

может быть неприменимым и аналитическое решение задачи будет некорректным.

Наиболее близкой из рассмотренных подобных систем будет система при неодинаковых давлениях фаз [2], в которой на каждую из двух сосуществующих фаз действуют

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

Тогда при общем изменении давления в баллоне и воздействии его на жидкую фазу

по уравнению Д. Г. Пойтинга

д£2

Э£\

Р

VI

и2

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

Вместе с этим относительно представленных выводов необходимо сделать следующие замечания. Во-первых, в отличие от системы с сосуществующими фазами и инертным газом по [2] в рассматриваемом случае в газовой полости баллона и термосе находится одно и то же рабочее тело, и эти разделенные по объемам составляющие не отвечают закону Дальтона. Во-вторых, уравнение Д. Г. Пойтинга выведено для условия Т=const в сосуществующих фазах, а это условие в рассматриваемом случае сохраняется только в объеме термосной емкости.

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

В настоящее время такой численный метод расчета реализован с применением вво -да теплофизических характеристик. Предва-

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

Полный объем универсального газового баллона находится из условия заполнения его газом при максимально допустимых давлении и температуре:

т • К • Т

V, ____ с тах

А _ Ртах

По значению массы заправки тс определяется объем внутреннего сосуда (термоса):

т

Р а

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

Тогда объем газовой полости определится как добавочный объем к объему термоса:

т • К •Т т

V _ V- — V _____________________________2____________тах_______________

У Р1 ї V А ' Г)

Р

Ра

_т

Я

К •От Ртах

тах г а л .

\

Ра

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

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

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

Оо + О, _ До • Р  О^Оа- +

О

и

Протекание процессов в баллоне и изменение параметров будут зависеть от режима работы баллона.

Основными такими режимами будут: режим заправки (заполнение криопродуктом); безрасходное хранение; работа в расход -ном режиме.

При этом все три режима могут иметь особенности по начальному состоянию баллона . Например, заправку криопродуктом можно осуществлять в «теплый» баллон (тем -

где 1Т, 1С - теплопроводность материала трубки и сильфона; ¥Т, ¥С - площади сечения стенки и сильфона; О ^ , Оа - текущие значения температуры стенки баллона и криопродукта (температура стенки термоса принимается равной температуре криопродукта

Оо = Оа ).

Тепло, подводимое к боковой поверхности термосной емкости, обусловлено теплопроводностью газа и тепловыми мостами (для предварительного этапа расчета последними можно пренебречь - это доказывается ниже) и определяется следующим выражением:

— Т — О

О _ Д • Р а а

\£ёс ёс уд

1

Здесь Дёс - эквивалентная теплопроводность системы экранов, теплоизоляция термоса и др.; Руд - среднее значение эффективной площади прохождения теплового потока; 5ёс - толщина газовых прослоек (радиальный зазор между стенкой баллона и стенкой термоса).

Тогда общие теплопритоки определятся суммой:

Q _ QT + Qc + Qeg _ /

_(oa -o„ )•

It/ ^c/ deg,

/AT /Ac / Aec

lT lc

до критического значения ( ГёЗ _ 0). Для ме-

тана С _ 511

112

ёАа _ ^ГГ ёАа ґ

------, r _ 255 ---------- (для

ёй ёй

предварительных оценок такое осреднение допустимо).

Из предыдущего выражения эквивалентная теплопроводность изоляции и газовых прослоек определится следующим образом:

A _______r-mcА________

ec (дл - da )• Fye • At'

Несложно показать, что комплекс

РТ + Рс

I / I / составляет менее 1 % от вели-

*Т/ с/

Температура баллона ОА понижается с 300 К до 200 К, т. е. среднее значение температуры баллона ОА = 250 К, а температура

криопродукта Оа повышается от 112 К до 190 К, т. е. среднее значение температуры криопродукта Оа = 150 К.

Вычислим значение теплопроводности изоляции:

F,.

Ув

чины § / . Таким образом, основные теп-

лопритоки будут осуществляться от боковых стенок баллона:

Q _ K'°'A S 0a • Fye 'At _ r-m,

eg

A _ 255-12-0,05-103

ec _ (250 -150)-1-3600-24

_ 0,018

Ад і -Е

Если время полной газификации криопродукта уменьшится до 2,4 часа, то теплопроводность изоляции, соответственно, изменит-

(7) ся до значения Aeg @ 0,2 -

Предположим, что время полной газификации криопродукта составляет 1 сутки, т. е. 24 часа, и за это время термические деформации и напряжения находятся в допустимых пределах.

Ставится вопрос: какого уровня должно быть значение коэффициента теплопроводности , обеспечивающее заданное вре-

мя газификации. Энергетический баланс будет следующим:

Здесь г - среднее значение теплоты испарения криопродукта в диапазоне давлений

Ад

I -Ё

Для метана критические значения температуры и давления: Т = 190 К, Р =

^ ^ ^ кр 5 кр

= 4,64 МПа. Если считать, что теплопроводность системы изоляции будет равна собственной теплопроводности газа, то среднее значение теплопроводности газообразного

__ Ад

метана [3] составит 1 = 0,03----- . Это по-

I -Ё

зволит получить продолжительность процесса испарения до 14 часов.

Уточним, до какой температуры охладится баллон при полной газификации криопродукта. При толщине стенки в 10 мм масса баллона = 90 кг (V = ж - 27 - 1 - 135 = 11445 т 3,

р = 7,8 , тА = 11445 - 7,8 -10-3 = 89 ёа ).

т

При захолаживании баллона от 300 К до 190 К без внешних теплопритоков потребуется испарить некоторое количество криопродукта, для СПГ это будет

AmNlA •r _ mA • na -(o0-190);

90-0,5-100

AmNiA _----------------_ 17,6 ёй

NIA 255

Масса заправки составляет тд = 12 ёа,

т. е. ктш > тд, и, соответственно, температура баллона будет выше 190 К, а именно:

ґупіп 'Т'

oA “ т0

_ 300 - 12-255

mA -NA

90 0,5

_ 230 Е » -40 °N.

Таким образом, из предварительных оценок следует, что при существующих допустимых уровнях давления заправки и температуры баллона (- 400С) [4] универсальный газовый баллон может иметь допустимую

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

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

Список литературы

1. Топливный баллон./ Довгялло А. И., Лукачев С. В. и др.// Патент № 2163699, Россия, МПК 7Б17С9/02, 99114577/06 заявл. 02.07.1999, опубл. 27.02.2001. Бюл. № 6.

2. Кириллин В. А., Сычев В. В., Шейн-длин А. Е. Техническая термодинамика. - М.: Энергия, 1968.

3. Загорученко В. А., Журавлев А. М. Теплофизические свойства газообразного и жидкого метана. - М.: Изд-во стандартов, 1969.

4. ГОСТ 949-73. Баллоны стальные малого и среднего объема для газов на Р £ 16,6 МПа (200 кгс/см2 ). Технические условия.

PECULIARITIES OF SETTING UP A PROBLEM OF CALCULATING PROCESSES FOR A GENERAL-PURPOSE GAS CYLINDER AS A NON-CLASSICAL TWO-PHASE SYSTEM AND ESTIMATION OF ITS TEMPERATURE CONDITION

© 2006 T. V. Ashikhmina, A. I. Dovgyallo Samara State Aerospace University

The paper describes the processes in a general-purpose gas cylinder different from cryogenic capacities with low pressure levels stored without draining in that heat isolation is inside the vessel. Existing approaches to solving the problem of process calculation for similar two-phase systems are analysed. The results of calculating cylinder temperature condition carried out using a simplified method are presented.