Динамика роста давления жидкости в замкнутом объеме при ее нагревании

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

МАШИНОСТРОЕНИЕ
В. Ш. Шагапов, Ю. А. Юмагулова
УДК 532. 58
ДИНАМИКА РОСТА ДАВЛЕНИЯ ЖИДКОСТИ В ЗАМКНУТОМ ОБЪЕМЕ ПРИ ЕЕ НАГРЕВАНИИ
Предложена и исследована модель повышения давления воды в замкнутом объеме при ее нагревании. Для случая плоской, радиально-симметричной и сферически-симметричной постановки задачи получены решения, описывающие рост давления вследствие термического расширения и фазовых переходов. Проанализировано влияние начального состояния жидкости, геометрического размера замкнутого объема на темп изменения давления. Повышение давления- термическое расширение- парообразование
ВВЕДЕНИЕ
В настоящее время представляют интерес исследования процессов теплопереноса в емкостях различной геометрической конфигурации. Это, в первую очередь, цилиндр и сфера. Такая геометрия обусловлена преимущественно потребностями ракетно-космической и авиационной промышленности, так как топливные баки и емкости для хранения криогенных жидкостей имеют чаще всего именно такую форму. Следует также отметить, что горючие жидкости, хранящиеся в замкнутых емкостях, во время пожара представляют собой опасность и могут привести к взрыву.
Нагревание жидкости в замкнутом объеме сопровождается термическим расширением и фазовыми переходами (испарением), которое может способствовать разрушению конструкции в результате сильного повышения давления. В ряде случаев причина возникновения подобных ситуаций может быть обусловлена только объемным расширением жидкости при ее нагреве. Но особенно сильный рост давления (тепловой удар) может реализоваться за счет фазовых переходов при больших тепловых нагрузках. Поэтому при конструировании и расчете надежности прочности различных технических систем, таких как гидравлические машины и другие устройства, которые работают в переменных температурных условиях, необходимо учитывать термическое расширение используемой жидкости и возможные фазовые переходы.
ПОВЫШЕНИЕ ДАВЛЕНИЯ ЖИДКОСТИ В ЗАМКНУТОМ ОБЪЕМЕ ВСЛЕДСТВИЕ ЕЕ ТЕРМИЧЕСКОГО РАСШИРЕНИЯ
Рассмотрим задачу о нагревании жидкости в виде плоского слоя толщиной 2Ь либо в виде цилиндрической или сферической области радиусом Ь, через боковые стенки. Нагревание происходит с поверхности г = Ь за счет поддержания заданной температуры. При отсутствии кипения жидкости повышение давления будет происходить только за счет термического расширения по всему объему (0 & lt- г & lt- Ь), где г = 0 -координата оси или центра симметрии.
Для теоретического описания процесса повышения давления в объеме при нагреве запишем уравнения неразрывности, теплопроводности и состояния жидкости в линейном приближении [1, 2]:
Р, С,
1, э
)= 0 ЭТ,
(1)
(2)
Контактная информация: 8−960−384−43−52
дг гп дг
эт
дг гп дг ^ дг
р, = р, о (1 — а (Т) (Т — То)+ а (р) (р — ро)), (3) где р, Х, с, и, Т, р — плотность, теплопроводность, теплоемкость, скорость, температура и давление жидкости, нижний индекс 0 здесь и далее соответствует начальным значениям параметров жидкости, а (Т) — коэффициент термического расширения жидкости, а (р) =
= V (Р 1С) — коэффициент сжимаемости, определяемый скоростью звука в жидкости С, п = 0, 1 и 2 соответствует плоской, радиальной и сферической геометрии задачи.
Пусть в исходном состоянии (г & lt- 0) температура жидкости равна Т, 0, а давление р0. С некоторого момента времени г = 0 на границе г = Ь поддерживается постоянная температура Те, ко-
торая выше начальной температуры Т, 0. При этом на оси симметрии емкости (г = 0) выполняется условие отсутствия тепловых потоков (дТ, / дг = 0).
При гипотезе однородности давления (условие гомобаричности др / дг = 0) из уравнения неразрывности (1) и теплопроводности (2) с учетом линейного уравнения состояния (3) можно получить
а (р) г
-л,(Т,(Т).
= а '-п
э
ёг ' дг I дг) дг
где У (Т-1 = 1,/(р, с,) — температуропроводность жидкости.
Учитывая, что скорость жидкости на границах г = 0 и г = Ь равна нулю, интегрирование уравнения (4) по координате от 0 до Ь приводит к следующему дифференциальному уравнению для изменения давления жидкости:
= (п + 1) а'-V
(Т)Л,(Т) 1
а (р)Ь
дг
(5)
Для задания теплового потока на границе г = = Ь необходимо найти решение уравнения теплопроводности (2), удовлетворяющее условиям Т = Т, 0 (г = 0, 0 & lt- г & lt- Ь), дТ, / дг = 0 (г & gt- 0, г = 0), Т = Те (г & gt- 0, г = Ь).
Решение этой задачи известно [3], подставляя в уравнение (5), для п = 0, 1 и 2 соответственно получим:
р = р0 +
^ 1
Х Е (о +1)2
т=0 (2т + 1)
8а (Т) (Те — Т0)
а (р) л2
х
1 — е
-п (Т ^ (2т+1)2 л: 2г/ (4Ь2)
р = р0 +
4а (Т) (Те — Т0)? ^
(р)а2.
а (р Ь
'-л
=1 а"
1 — е-
р = р0 +
х Е
т=1 т
6а (Т) (Те — Т0)
а (р) л2
х
1 — е
-V т л г,
: л2г/б2
(6)
(7)
(8)
Следует отметить, что из решений (6)-(8) можно получить предварительные предельные оценки роста давления р (М) предельным переходом при г ^ да или из уравнения состояния (3), полагая Т, = Те и р, = рю:
ЛМ)
а
= р0 + ^7(Те -Т0).
(Т)
а
(р)
(9)
В численных расчетах применительно к воде использовались следующие физические параметры [4]: а (г) = 10−4 К-1, а (р) = 10−9 Па-1,
V
(Т)
= 1.3 10 7 м2/с. Для начального состояния приняты значения температуры и давления воды Тю = 293 К, р0 = 0,1 МПа, температура на границе равна Те = 353 К.
На рис. 1 представлены зависимости давления от времени, иллюстрирующие влияние линейных размеров замкнутого объема. Здесь и в дальнейшем сплошная, штрихпунктирная и точечная линии соответствуют случаям плоской, радиальной и сферической геометрии задачи. Штриховая линия, выполненная по формуле (9), соответствует величине максимального давления жидкости в замкнутом объеме.
Рис. 1. Динамика роста давления вследствие термического расширения воды при различных значениях полуширины или радиуса Ь
Видно, что при нагревании воды в замкнутой емкости, которая препятствует ее расширению, появляется давление на стенки, которое может привести к ее разрушению. Величина максимального давления не зависит от геометрического размера емкости, а зависит от начальных температур границы и воды. Геометрической размер емкости влияет на темп повышения давления.
ПОВЫШЕНИЕ ДАВЛЕНИЯ
ЖИДКОСТИ ВСЛЕДСТВИЕ ПАРООБРАЗОВАНИЯ
Рассмотрим случай, когда при нагревании емкости с жидкостью через боковые стенки может происходить вскипание жидкости вблизи границы. В этом случае появляется слой пара (г («) & lt- г & lt- Ь), в котором можно записать уравнение теплопроводности в виде
(дТъ дТъ I д (п дТв I пп.
Р-с{1Г+ь--Т)=I/ -г)¦ (|0)
и
Ь
г
т
где р», Х", Сщ, ищ, Т-о — плотность, теплопроводность, теплоемкость, скорость и температура пара.
Для пара примем уравнение Менделеева-Клапейрона
рщ=риКТи, (11)
где ри — давление пара, — газовая постоянная.
На границе раздела пар-жидкость (г = г^) температура равна равновесной температуре фазовых переходов для давления р: Г =
= Тг =Г (р). Кроме того, на этой границе выполняются условия баланса массы и тепла:
(12) (13)
к- г,))=р, (и — г,)),
РДЧ
. ЭТ,, ЭТ _
1'- аТ& quot--1^
= РиГ (=
где рщ — средняя плотность пара, , — удельная теплота парообразования.
Начальное состояние (г & lt- 0) характеризуется температурой жидкости Т, 0 (0 & lt- г & lt- Ь), давлением р0, а толщина слоя пара равна нулю (5 = Ь -- г^) = 0). С некоторого момента времени г = 0 на границе г = Ь поддерживается постоянная температура Те, которая выше температуры кипения жидкости при давлении р0, в центре емкости (г = 0) тепловые потоки отсутствуют (дТ,/дг = 0).
Аналогично предыдущему случаю уравнение для изменения давления жидкости будет иметь вид (4). Проинтегрировав это уравнение по координате от 0 до г^) с учетом малости слоя пара (5 & lt-<- Ь) и его плотности (рщ & lt-<- р,), а также используя условие баланса массы (12), получим
'-(р)Ь. = а (Т-V (Т)| ^ | +5. (14)
а
(п + 1)
щ +6.
Эг І
Согласно уравнению (14) повышение давления жидкости происходит за счет термического расширения и сжатия слоем пара.
С учетом предположения, что толщина парового слоя мала (5 & lt-<- Ь), для распределения температуры из (10) в паровом слое примем линейный закон:
(Т (р)-Те)(ь — г)
Т= Т ±
6
(15)
В уравнении (14) неизвестной величиной является толщина парового слоя 5, которая может быть определена из условия баланса тепла (13) с учетом (15) в виде:
6 = ^ Ри*
(
11 Iі
Эг
Те — Т (р) 6
(16)
Для дальнейшего примем некоторые упрощения. В уравнениях (14) и (16) будем пренебрегать первым слагаемым правой части, отвечающим за термическое расширение. Это означает, что повышение давления жидкости будет происходить только в основном из-за сжатия вследствие образования пара. Тогда система уравнений (14) и (16) относительно неизвестных 5 и р преобразуется к виду
а (р)Ьёр = (п +1)5, 5=1и (Те_- Т (р)). (17) ёг ри, 5
Отсюда с учетом уравнения состояния пара (11) и выражения Тх (р) = Т*п-1(р* / рж) из [2] для зависимости равновесной температуры фазовых переходов от давления (Т* =, / Лщ, р* - эмпирические параметры, определяемые на основе табличных данных), для закона изменения давления жидкости в емкости получим уравнение
(п + 1)2
Ж (а (р)Ь)2 р (р — Ро)1п (р./ р)
Т --
Т.
1п (р./ р) (18)
Для начального этапа закона роста давления, когда (р — р0) & lt-<- р0, из (18) имеем следующее решение:
р = р0+ (п±1). (19)
а (р)Ь
р01п (р•/ р0)
Формула (19) необходима для задания начальных условий Коши при дальнейшем численном продолжении решения уравнения (18). При этом закон роста толщины слоя пара на фоне полученного решения для давления определяется интегрированием второго уравнения
(17).
В численных расчетах для температуры границы, начальной температуры воды и давления использованы значения: Те = 673 К, Т, 0 = 293 К, р0 = 0,1 МПа. Для теплофизических параметров воды и пара приняты величины: Х =
= 0,024 Вт/(мК), Т* = 5,05 ¦ 103 К, р* = 7,67 х х 1010 Па.
На рис. 2 представлены графики зависимости давления в объеме (а) и толщины слоя пара 5 (б) от времени при различных значениях полуширины или радиуса емкости Ь.
Видно, что при сильном перегреве вблизи границы емкости сразу начнется интенсивное парообразование, которое в свою очередь приведет к быстрому повышению давления жидкости во всем объеме. При этом на начальном этапе давление достигает определенного значения и в дальнейшем меняется незначительно. Это
связано с тем, что равновесная температура Тц{р) фазовых переходов, соответствующая значению давления в объеме р становится равной температуре границы нагревания Т. е.
р, МПа
О 36?. мин
б
Рис. 2. Динамика роста давления в объеме вследствие парообразования (а) и толщины слоя пара (б) при различных значениях Ь
Таким образом, процесс парообразования вблизи границы замкнутой емкости продолжается некоторое время и впоследствии прекращается.
На рис. 3 представлены профили температуры вблизи границы по координате х = Ь — г в моменты времени г = 1 (а) и 10 мин (б) соответственно. Точки излома на кривых соответствуют границе фазовых переходов.
Видно, что с течением времени область пара уменьшается, а прогретая область воды увеличивается, что приводит к тому, что через некоторое время на рост давления в объеме термическое расширение воды будет оказывать большее влияние.
т. к
600 —
400
«я X, м
а
т, к 600 —
400 —
0,00 0 01 с 0,02 Х. м 0,03
б
Рис. 3. Профили температуры вблизи границы
ЗАКЛЮЧЕНИЕ
Для случаев плоской, радиальной и сферической геометрии получены аналитические и численные решения, описывающие повышение давления жидкости в замкнутом объеме за счет термического расширения и возможного парообразования.
Показано, что нагревание жидкости на границе замкнутого объема приводит к повышению ее давления, которое определяется интенсивностью нагревания. При этом возможны два режима повышения давления. При первом режиме, когда температура границы нагревания меньше температуры кипения для начального давления жидкости, повышение давления происходит за счет ее термического расширения, темп его зависит от условий нагревания и геометрии замкнутого объема. Интенсивный рост давления наблюдается при втором режиме, когда появляется область парового слоя. В этом случае процесс парообразования длится в течение конечного промежутка времени, а затем паровой слой схлопывается. Это связано с достижением давления значения, при котором температура насыщения становится равной температуре границы нагревания и дальнейшее повышение давления возможно только за счет термического расширения жидкости.
СПИСОК ЛИТЕРАТУРЫ
1. Шагапов В. Ш., Ильясов У. Р., Насыро-
ва Л. А. Тепловой удар в пористой среде, насыщенной жидкостью // Теплофизика и аэромеханика, 2003. Т. 10, № 3. С. 411−422.
2. Нигматулин Р. И. Динамика многофазных сред. М.: Наука, 1987. Ч. 1. 464 с. Ч. 2. 360 с.
3. Карслоу Г., Егер Д. Теплопроводность твердых тел. М.: Наука, 1964. 487 с.
4. Варгафтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей. М.: Наука, 1972. 720 с.
ОБ АВТОРАХ
Шагапов Владислав Шайхулагзамович, академик АН РБ. Д-р физ. -мат. наук.
Юмагулова Юлия Александровна, асп. Бирск. филиала Башкирск. гос. ун-та.

ПоказатьСвернуть
Заполнить форму текущей работой