Построение математической модели переходных процессов в газопроводах с учетом неизотермичности и сжимаемости потока
Бесплатный доступ
Предложена математическая модель переходных процессов в газопроводах, в явном виде определяющая распределение температур по длине газопровода. Кратко описан вывод системы дифференциальных уравнений, связывающих основные параметры газового потока. Структура полученных уравнений лучшим образом подходит для последующего применения численных схем.
Газопровод, переходный процесс, математическая модель
Короткий адрес: https://sciup.org/147154342
IDR: 147154342
Текст краткого сообщения Построение математической модели переходных процессов в газопроводах с учетом неизотермичности и сжимаемости потока
В литературе представлено большое количество математических моделей потоков газа в газопроводах [1–4], неизотермичность процессов при этом обычно учитывается в виде уравнений, включающих изменение внутренней энергии, энталь-
пии, энтропии и т. д., но не конкретно температуры газа. В представленной статье из фундаментальных уравнений механики жидкости и газа
Здесь У * - скорость присоединяющихся (отделяющихся) частиц перед их присоединением (после отделения) от основного потока; д - ускорение свободного падения; у – ордината оси х , отсчитываемая от горизонтальной плоскости; А - коэффициент гидравлического трения.
Принимая те же положения о притоке и диаметре трубы, что и в п. 1, имеем:
получена математическая модель переходных
процессов в газопроводе, состоящая из четырех уравнений, одно из которых включает изменение
Э(рУ) Э(р+рУ2) = _(dy ЛЩЩ Х
Эt Эх р(д dx+ 2D ) .
температуры в явном виде.
Широко распространена модель, представленная в [1]. Основные дифференциальные урав-
нения, описывающие неустановившееся движение газа в трубе, выводятся из законов сохранения массы, импульса и энергии. Движение газа предполагается одномерным. При этом площадь сечения трубы ω – переменная величина. Из сопротивлений рассматриваются лишь сопротивления трения газа о стенку трубы.
1. Уравнение потока массы:
^(ры)+-Эх(рУы) _nDm*. (1)
Здесь ρ – плотность газа; ω, D – площадь поперечного сечения и диаметр трубы соответственно; т * -распределенный приток массы через боковую поверхность, приходящийся на единицу площади трубы в единицу времени. В случае транспортировки газа по трубопроводу, очевидно, что никакого притока массы не будет, и для моделирования примем участок трубы постоянного диаметра. Таким образом, уравнение потока массы прини-
Приведем полученное уравнение к более подходящему для моделирования виду, используя базовые правила нахождения производных сложных функций и уравнение потока массы. Опустив вывод, вследствие ограниченности объема статьи, получаем:
р^ + ^+ри^-рГд^У + т к -t Эх Эх к \у dx 2D /
3. Уравнение энергии
Согласно [1] уравнение энергии имеет вид
Э-(рЕ) + Э-рУ(Е + ^) =
Эtvp v Эх^ X р/
4q = D
-
рУд^ + 1^ * (Е + Е ). ' а dx D X р/
мает вид:
ар + a(pu) _ о Эt Эх
ар + и ар + р^_о. Эt Эх м Эх
2. Уравнение потока импульса
Базовая формулировка этого уравнения также принята из [1] и имеет следующий вид:
у (рУш) + "Э[(Р + рУ2)ш] = at Эх
, . dto ( dy A|U|UX ,
(р ■ m . uJ ^л.ф^. . !. (4)
Здесь [1] E – полная энергия единицы массы, равная сумме внутренней u и кинетической энергии У 2 /2:
Е = и + У^.
Теплопередача от газа, текущего в трубе, к внешней среде с известной температурой Т н задается в виде закона Ньютона:
q = ас(Тн - Т).
Как уже отмечалось выше, приток массы через боковую поверхность отсутствует, значит, уравнение (3.1) можно представить в виде
£(р(и+9)+Хи+^+р)_
_ 4q ~ D
-
руд ДУ .
Для дальнейшего преобразования были использованы следующие зависимости:
– из курса термодинамики известно, что удельная энтальпия вычисляется по формуле i = и + р
Краткие сообщения
др д(рП)
- уравнения потока массы — = —дХ" ;
- уравнения потока импульса Р д" + Р" дх =
Р-,^-. ^'"j к дх 2D /
др дх .
Используя перечисленные уравнения и правила нахождения производных, получаем:
р(д1 + и^1) = 4£ + дР + р^№ + идР. (11)
м \д£ дх/ D д1 м 2D дх v 7
Теперь обратимся к основам термодинамики. Известно [1, 2], что приращение энтальпии можно определить по следующей формуле:
Мд^+^р.
Частная производная от энтальпии по температуре при p = const равна теплоемкости тела, т. е. ■^
Также известно выражение для частной производной от энтальпии по давлению при постоянной температуре [2], которое имеет следующий вид:
где V - удельный объем газа: V = р .
Относительное изменение объема тела при нагревании его на dT градусов при постоянном давлении есть коэффициент объемного теплового расширения: ^ = ^ (^) .
V \oi / Р
Подставляя вышеприведенные равенства в уравнение (11) и производя необходимые математические преобразования, получаем:
рс»^."?) = тв-дЕ."^) +
' р Vdt дх/ 1 Vdt дх/
! 4р . W
D Р 2D
.
Для нахождения частных производных давления в правой части уравнения (14) применим уравнение состояния газа. Для реальных газов оно имеет следующий вид:
р = ZpRT,
где R – газовая постоянная, Z – коэффициент сжимаемости.
Продифференцировав (15) по времени и коор-
динате и также используя уравнение неразрывности, приводим уравнение (14) к виду
(рср-рв)£+иЭ + рвтд" =
= eRpT2-^^."^^).4^. р1£". (16)
-
1 ' V д£ дх/ D ' 2D v 7
-
4. Начальные условия
Таким образом, окончательно получаем систему дифференциальных уравнений, описывающих движение газа в газопроводе, состоящую из четырех уравнений (3), (6), (16) и (15).
Перед началом переходного процесса, инициированного тем или иным фактором, поток считается установившимся. То есть для получения уравнений, описывающих начальные условия, достаточно из полученной системы дифференциальных уравнений откинуть производные по времени, входящих в нее величин. Далее, применяя выбранную численную схему, получаем обычные зависимости искомых величин от координаты х.
Вывод. Вышеприведенные уравнения описывают динамические процессы в газопроводах. Для полноты формулировки математической модели к данным уравнениям необходимо присоединить еще граничные условия, указывающие конкретные особенности протекания изучаемых процессов. В связи с многовариантностью граничные условия не приводятся, их формализация не представляет каких-либо затруднений.
Список литературы Построение математической модели переходных процессов в газопроводах с учетом неизотермичности и сжимаемости потока
- Неизотермическое течение газа в трубах/под ред. О.Ф. Васильева. -Новосибирск: Наука, 1978. -128 с.
- Вукалович, М.П. Термодинамика/М.П. Вукалович, И.И. Новиков. -М.: Машиностроение, 1972. -672 с.
- Панферов, В.И. Численное моделирование переходных процессов в газопроводах/В.И. Панферов, А.А. Февралев//Вестник ЮУрГУ. Серия «Строительство». -2008. -Вып. 7. -№ 25(125). -С. 40-45.
- Чарный, И.А. Основы газовой динамики/И.А. Чарный. -М.: Гостоптехиздат, 1961. -196 с.