Товары в корзине: 0 шт Оформить заказ
Стр. 1 

43 страницы

349.00 ₽

Купить Р 401-80 — бумажный документ с голограммой и синими печатями. подробнее

Распространяем нормативную документацию с 1999 года. Пробиваем чеки, платим налоги, принимаем к оплате все законные формы платежей без дополнительных процентов. Наши клиенты защищены Законом. ООО "ЦНТИ Нормоконтроль"

Наши цены ниже, чем в других местах, потому что мы работаем напрямую с поставщиками документов.

Способы доставки

  • Срочная курьерская доставка (1-3 дня)
  • Курьерская доставка (7 дней)
  • Самовывоз из московского офиса
  • Почта РФ

В Руководстве изложена методика расчета газодинамических параметров при заполнении газопроводов диаметрами 1020, 1220 и 1420 мм природным газом для испытаний, приведены программа ASTRA для расчета на ЭВМ БЭСМ-6 процесса заполнения газопровода и приближенный аналитический метод расчета.

 Скачать PDF

Оглавление

1 Общие положения

2 Методика расчета газодинамических параметров при заполнении магистральных газопроводов газом для испытаний

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

Приложения

Литература

 
Дата введения03.07.1980
Добавлен в базу01.09.2013
Актуализация01.01.2021

Этот документ находится в:

Организации:

03.07.1980УтвержденВНИИСТ
РазработанВНИИСТ
ИзданРотапринт ВНИИСТа1981 г.
Стр. 1
стр. 1
Стр. 2
стр. 2
Стр. 3
стр. 3
Стр. 4
стр. 4
Стр. 5
стр. 5
Стр. 6
стр. 6
Стр. 7
стр. 7
Стр. 8
стр. 8
Стр. 9
стр. 9
Стр. 10
стр. 10
Стр. 11
стр. 11
Стр. 12
стр. 12
Стр. 13
стр. 13
Стр. 14
стр. 14
Стр. 15
стр. 15
Стр. 16
стр. 16
Стр. 17
стр. 17
Стр. 18
стр. 18
Стр. 19
стр. 19
Стр. 20
стр. 20
Стр. 21
стр. 21
Стр. 22
стр. 22
Стр. 23
стр. 23
Стр. 24
стр. 24
Стр. 25
стр. 25
Стр. 26
стр. 26
Стр. 27
стр. 27
Стр. 28
стр. 28
Стр. 29
стр. 29
Стр. 30
стр. 30

МИНИСТЕРСТВО СТРОИТЕЛЬСТВА ПРЕДПРИЯТИЙ НЕФТЯНОЙ И ГАЗОВОЙ ПРОМЫШЛЕННОСТИ

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

ВНИИСТ

рип

ПО РАСЧЕТУ ОСНОВНЫХ ПАРАМЕТРОВ ЗАПОЛНЕНИЯ ГАЗОПРОВОДОВ ГАЗОМ ДЛЯ ИСПЫТАНИЯ

Р 401-80

Москва 1981

МИНИСТЕРСТВО СТРОИТЕЛЬСТВА ПРЕДПРИЯТИЙ НЕФТЯНОЙ И ГАЗОВОЙ ПРОМЫШЛЕННОСТИ

ВСЕСОЮЗНЫЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ПО, СТРОИТЕЛЬСТВУ МАГИСТРАЛЬНЫХ ТРУБОПРОВОДОВ

ВНИИСТ

ПО РАСЧЕТУ ОСНОВНЫХ ПАРАМЕТРОВ ЗАПОЛНЕНИЯ ГАЗОПРОВОДОВ ГАЗОМ ДЛЯ ИСПЫТАНИЯ

Р 401-80

Москва 1981

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

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

* Окончание расчета соответствует моменту, когда по всей длине трубопровода давление газа будет равно испытательному давлению с учетом погрешности измерительной аппаратуры, т.е.у:

P(o,L) - P(L,t) = 0,005 фит    аэ)

2.10.    Расчет газодинамических параметров (давления, температуры, расхода) в процессе заполнения участка магистрального газопровода га&м для испытания выполняется по программе

РЕАЛИЗАЦИЯ МЕТОДА НА ЭВМ

2.11.    Алгоритм программы Я$ТМ и подпрограмм PRQG0H Fi , RLRM- реализован в виде программы на алгоритмическом языке ФОРТРАН для ЭВМ БЭСМ-б (прил.2).

В программе fiSTRR приняты следующие обозначения, приведенные в табл.1.

Таблица I

Сшволы ал- ] Идентификаторы ь Наименование параметров и горитма ; программы |    массивов

---:-■-у------------ '    ---

L

Й1

Протяженность заполняемого участка

D

В

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

Рисп

F1MAX

Испытательное давление

ТС

Критическая температура газа

PC

Критическое давление газа


W'.f—'....... '}    ■    "■     -    —

Символ ал-(Идентификаторы ; Наименование параметров и мае-горитма ; программы ]    сивов

Ггр

TGR

Температура грунта

Р0

PCEP

Давление в источнике испытательной среды

Гс

TCEP

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

Р»

PH

Начальное давление газа в заполняемом трубопроводе

Тн

TH

Начальная температура газа в заполняемом участке газопровода

GH

Начальный расход газа в заполняемом трубопроводе

Pi

SR

Коэффициент Джоуля-Томсона

К

RH

Коэффициент теплопередачи

X

RLRM

Коэффициент гидравлического сопро

тивления

г

2

Коэффициент сжимаемости газа

об

RIFP

Скорость подьема давления

9

Q.

Плотность теплового потока через стенку трубы

Р(о,ь)

Fl(x)

Граничное условие для давления

Г(0, Ъ)

F2 (X)

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

G(L}b)

GK(N)

Граничное условие для расхода

/dz\

\ЭР/т

ZP

Вспомогательные функции

Г

(bZ\

2T

\5t7p

Re

RE

Re,

R El

"Числа Рейнольдса

Rek

REZ

>

N~1

N-l

Число интервалов разбиения об

k

H

ласти интегрирования

Шаг по координате при численном

T

TRU

THUN

интегрировании

V

w Шаги по времени при численном интегрировании

II

Символ ал- [идентификаторы ] “Наименование параметров и иас-горитма I программы ;    сивов


1 3 N7 NP NPX


N


Рабочие переменные целого типа


У


2.12.    В программе RSTRA (см.прид.2) приняты следующие значения термодинамических величин (газ - метан):

R = э15 Дд/кг*К'- газовая постоянная;

Ср = 2930 Дж/кг-К - теплоемкость газа при постоянном давлении;

Ркр= 44,5 ■Ю"’ Н/м2 - критическое давление газа;

Ткр= 190,7 К - критическая температура газа;

Dj, = 4,08*10“^ м2*К/Н - коэффициент Джоуля-Томсона ,

2.13.    При заполнении трубопровода природным газом с другими термодинамическими характеристиками необходимо изменить коэффициенты Al, А2, АЗ, А4, АЗ в соответствии с формулами:


AI =

2 кЛ

srfifd5

>

А2 = ■

51 д-£г

4R

>

It

4ft2

71^ СрДг

»

А4 =

сРр

1

А5 =

R


(20)

(21)

(22)

(23)

(24)


2.14. При реализации программы расчета газодинамических параметров на ЭВМ на печать выводим значения давления (РР),


12


температуры (ТТ) и расхода (    )    газа    внутри заполняемого

газопровода в дискретных точках-узлах по длине трубопровода во времени (TIME).

Для каждого момента времени (TIME) значения найденных функций Р, Т и & печатаются построчно (прил.З). Количество выдаваемых на печать точек определяется значением операторов NPX (по координате х ) и ' nP (по времени).

2.15. Ра§мв£|Ость выводимых на пеяать значений:

РР - давление, кгс/см2;

ТТ - температура, К;

GG- расход, м3/ч;

TIME - время, ч.

2.16. Первое число в строке (см.прил.З), определяющее величину данного параметра (РР), (ТТ), (OG ), соответствует сечению х = 0.

для определения расстояния между следующим сечением, в котором выдаются на печать рассчитанные параметры,и началом заполняемого участка ( X s 0) следует воспользоваться формулой

(25)

^ L. NPX , N- 1

где п - порядковый номер числа в соответствующей строке СРРЗт (ТТ), (вЪ);

L - протяженность участка;

(N-l) - количество интервалов разбиения области интегри-у рования;

/ХРХ - количество точек (узлов), через которое печатаются рассчитываемые параметры ( NPX = I - параметры печатаются во всех точках, NPX = 2 - через две, начиная с точки X =0, и т.д.).

3. АНАЛИТИЧЕСКИЙ МЕТОД РАСЧЕТА ПРОЦЕССА ЗАПОЛНЕНИЙ УЧАСТКА МАГИСТРАЛЬНОГО ТРУБОПРОВОДА ГАЗОМ ДЛй ИСПЫТАНИИ '

3*1. Изотермический процесс заполнения замкнутого участка трубопровода постоянного диаметра природным газом до испытательного давления описывает система уравнений (26) - (28)* которая вытекает из системы уравнений (1)-(4) при Т * const

bfi

It


4“


9(jpw) 3 x


= 0;


уравнение движения:


4-


Эх


(p\Pw2)=-/>



(26)

(27)


уравнение состояния газа:

P=J>CZ


(28)


Ё уравнениях (26) - (2d) приняты обозначения:

fi - средняя по сечению трубы плотность газа, кг/м3;

W - средняя по сечению трубы скорость движения газа,м/с,; р - среднее по сечению трубы давление газе, Н/м^;

X - коэффициент гидравлического сопротивления;

D - диаметр заполняемого трубопровода^; с - скорость звука в газа, м/с.


3.2.    Аналитические расчетные формулы получим, приняв для решения системы уравнений (26) - (28) следуюи^е допущения:

пренебрегаем конвективным членом ^    ;

Эх

пренебрегаем изменением массовой скорости (в^Л во времени, т.е.    ^    о    ;

принимаем коэффициент гидравлического сопротивления постоянным X = const.

3.3.    В качестве начальных условий для рассматриваемой системы уравнений принимаем постоянное давление по длине трубопровода


Р М)-Рн .


(29)


3.4. Граничные условия, отражающие процесс подъема давления по линейному закону на левом конце трубопровода (Х = 0) и перекрытия правого конца при заполнении газопровода газом/ приняты следующими (рис.З):


14


P(o,t)=


рн

p

r «сл


при


(30)

ПР* t & t, ;    /


a p


0,


(31)


X=L


где ос - скорость подъема давления во входном сечении заполняемого учаотка трубопровода

а -Р х _ иоп. 'н

о- оС



Рис.З. График изменения давления газа (Р) при заполнении участка магистрального газопровода газом для испытания:

1 - изменение давления в сечении * & 0;

2 - изменение давления в сечении х = ^


15


3.5. Изменение давления газа в заполненной трубопроводе подчиняется соотношению, которое было найдено в результате решения линеаризованной системы уравнений (26) - (28) при соответствующих начальных и граничных условиях:


где


i

а *.•«:'

(33)

^ S срС. Л» L

"" ззсгс2Д'

(34)

.. _ к(гп.-и) .

(35)

" 2L

+. — .^.ЧИ1 ~ Р Н 0 od

(36)

РЭфф * з.эг-ю5 н/“2

- эффективное давление, значение которого подучено из численных расчетов газодинамических параметров при заполнении трубопроводов газом на ЭВМ.


3,6. С достаточной степенью точности для приближенного аналитического описания изменения давления газа внутри заполняемого участка можно ограничиться первым членом разложения ( И = О):


=<


"    2oCs6a/~)    Г    ,    г,Л

[exp(^t0)-lJ'eyf>(-(u<Bi) при t =* t0.


при 0<=t«C t :


(3?)


Изменение давления газа в конце заполняемого участка определяем по формуле:


rti


P(L th


(38)


■exp -


при


'о •


3.7. Задаваясь допустимой погрешностью измерения давле -ния в конце трубопровода (£,) (см.рис.3)


(«)


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

r3tjn?£3P*


ZtLKu

L    __

зхЧг1/Р.



3&, сгЛ Х3Р;



з0><р


3.8. Пример приближенного аналитического расчета газо -динамических параметров в процессе заполнения участка магистрального трубопровода природным газом для испытания приведен в прил.4.


Приложение I

АЛГОРИТМ ПРОГРАММЫ R$TRfi И ПОДПРОГРАММ PRO&ON, W, RLRN, АЛГОРИТМ ПРОГРАММЫ ASTM

0-'


I

Описание массивов.

2

Описание блоков

COMM ON

3

Описание оператора-функции граничного условия для температуры F2 (х)

4

Ввод постоянных ИСХОДНЫХ данных

5

Ввод параметров явной конечно-разностной схемы

К, N^HPxjhUJAVN, NT

6

Ввод начальных условий рн,тн,С>н

7

Печать начальных условий

ТТ

8

Начало цикла по времени

9

Начало цикла по координате (I) Г=2**-ы1

Т"

К

Обращение к подпрограмме ALAM (&&) Вычисление а1

-0


■ II

вычисление 2,2Т, HP

12

Вычисление 2Е1

13

Вычисление

ФМК’(те/?-Г(1))




©-



15

Вычисление Р«(ОвК*)*])Ри ТН(Г)-Т(1)ЗД

16

Вычисление Pnfi) «Pflj TN(r)=T(IJ

17

continue



Вычисление в граничных

19

узлах PNpKTH(f|,

РН/Н)х th(n), Gn(n) . . Обращение к п/л FIOO

—*1

20

Начало цикла по координате^ (Г)

-кв>

21


УДК 622.691.47.001.2

В настоящей Руководстве изложена методика расчета газодинамических параметров при заполнении газопроводов диаметрами 1020, 1220 и 1420 мм природ -ним газом для испытаний, приведены программа RStrr для расчета на ЭВМ БЭСМ-б процесса заполнения газопровода и приближенный аналитический метод расчета.

Руководство разработано на основе теоретиче -ских исследований и натурных наблюдений при испытаниях магистральных газопроводов лабораторией надежности конструкций и методов испытаний магистральных трубопроводов ВНИИСТа совместно с лабораторией экономического анализа проектных решений в нефтегазовом строительстве НИПИЭСУнефтегазстроя.

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

Руководство составили инженеры Ю.В.Колотилов, Е.М.Климовский, А.Й.Тоут (ВНИИСТ), кацд.техн.наук БЛ.Кривонеив, канд.хим.наук Е.Я.Клименко®, канд. техн.наук'В.И.Агапкин, инж.М.Я.Куракина (НЙПИЭСУ-нефтегазстрой). В составлении Руководства принимали участие д-р техн.наук П.И.Тугунов (Уфимский нефтяной институт), канд.техн.наук В.В.Постников (Главсибтру-бопроводстрой).

(п) Всесоюзный научно-исследовательский институт по строи-w тельству магистральных трубопроводов (ВНИИСТ), 1981

©


бо)

21


Обращение к п/п ^лМ(бб)


Вычисление AI


22


Вычисление 22



37

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

38

Обращение к п/п р1(х)

39

Обращение к п/п PROG 0 а/

40

Получение массива давления

41


ВНИИСТ

ниливс;



Руководство


по расчету основных параметров заполнения газопроводов газом


для испытания


Р 401-80


I. ОБЩИЕ ПОЛОЖЕНИЯ

1.1.    Руководство разработано в развитие глав СНиП П-45-75 "Магистральные трубопроводы. Нормы проектирования" и

СНиП Ш-д. 10-72 "Магистральные трубопроводы. Правила производства и приемки работ".

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

1.3.    Одним из этапов комплекса работ по пневматическому испытанию является этап заполнения газопровода природным газом до создания в нем давления, необходимого для испытания газопровода на прочность.

1.4.    Источником природного газа для испытания магистрального газопровода служат:

месторождения газа, примыкающие к газопроводу или расположенные непосредственно около него;

действующий газопровод, пересекающий строящуюся магист -раль или проходящий непосредственно около нее.

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

1.6.    При пневматическом испытании газопроводов скорость подъема давления в них регламентирована СНиП 1-Д. 10-72.

1.7.    При испытании на прочность и герметичность для измерения давления применяют манометры класса точности не ниже 1,5.

Внесено ЛНИТ ВНИИСТа ! Утверждено ВНИИСТом ! Разработано и ЛЭАПР НШШЭСУнеф-    3    июля    1980 г.    впервые

тегазстроя    i    __

* - ______ * **


3


Рие.1. Принципиальная -схема наполнения участка магистрального тру -бопровода природным газом для ис -питания:

a-от промысла: б-иэ действующего газопровода; 1-источник газа; П-за-полняемый участок

2. МЕТОДИКА РАСЧЕТА ГАЗОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ПРИ ЗАПОЛНЕНИИ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ ГАЗОМ ДЛЯ ИСПЫТАНИЙ

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

2.1. Система, описывающая нестационарное неиэотермическое течение газа внутри заполняемого газопровода (постоянного диаметра J) ) имеет следующий вид:

*

4

уравнение неразрывности:

cd

Э.я , 9QwJ .

9t 9 х “ и'

(2)


уравнение импульсов:

уравнение энергии:

д

J \лЛ

0

л ( Р м1)

и

rrzj

+Н

ГК *2j

уравнение состояния:

(4)

= в/?Т

Р_

f

В уравнениях (I) - (4) приняты обозначения:

О

Р - среднее по сечению трубы давление газа, Н/м ;

jo - средняя по сечению трубы плотность газа, кг/м3;

w - средняя по сечению трубы скорость газа, м/с;

Л - коэффициент гидравлического сопротивления;

д - плотность теплового потока через утенку трубы,Вт/м2;

Т - средняя по'сечению трубы температура газа, К;

и - удельная внутренняя энергия газа, м22;

ъ    -    коэффициент сжимаемости газа;

R    -    газовая постоянная, Дж/кг.К;

Я - геометрическая высота центра тяжести рассматриваемого элемента объема газа, м;

р

д - ускорение свободного падения, м/с4*;

^    - диаметр заполняемого трубопровода, м;

Jf - коэффициент теплопроводности газа, Вт/м*К; х    -    координата вдоль оси трубы, м;

t    -    время, с.

5

2.2. Для получения искомых параметров Сдавления, темпег ратуры и расхода) при заполнении газопровода газом используем нелинейные дифференциальные уравнения (I) - (4), при этом принимаем следующие допущения:

теплопередача вдоль оси трубопровода за счет теплопроводности газа отсутствует;

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

2.3. Закон теплообмена между стенкой трубы и грунтом принимаем в форме Ньютона;

$~-к(Т-Тгр) ,    (5)

где Тгр - температура грунта, К;

к - коэффициент теплопередачи от газа в грунт, Вт/м2.К.

Р Т

нр '


Я 71


U 0,07


н


(6)


2.4.*' Для определения коэффициента сжимаемости реального газа используем аналитическую зависимость 2=*Е(Р,Т) типа Бертдо, справедливую в закритической области температур и при умеренных давлениях (до 98,1*10'’ Н/м2):

где Ркр - критическое давление газа, Н/м2; Ткр - критическая температура газа, К.

2.5. Коэффициент гидравлического сопротивления X как

функции числа Рейнольдса, диаметра и гидравлической.верахова-

тости определяем по формулам:

II

^5

при Re 200Q; (?)

Л = 0,0025‘(Re ft

при 2000< Re«4ООО; (8)

Л = 0„ЗШ •^Re'j0,25

при 4000 <■ Re ^ Ret; (9)

при Re,cRes£ Rez; (10)

при Re , (II)

40

R 9 7Г-Т~гг 9

51 j/ll)

(12)

6

Re


t


ss.s

у, y/7


(13)


где Re

Re1}Rea

К.



66S-765-U^A

*■*- —жг—~ >

число Рейнольдса; ^ переходные числа Рейнольдса; эквивалентная'гидравлическая шероховатость, и; коэффициент динамической вязкости, Н»с/м2; расход газа, Н/с.


(14)


2.6. Чтобы найти однозначное решение системы уравнений (I) - (4), описывающих процесс заполнения газопровода газом для испытаний, необходимо задать начальные и граничные условия.

В качестве начальных условий в практических расчетах задаем распределение давления (Р), температуры (Т) и расхода газа (Ь) по длине газопровода в начальный момент времени:'


Р(х>0) —Рн >    Т(х50)=Тн    ,    &(х,0)«&я    .    (15)


В качестве граничных условий задаем изменение давления и температуры газа в начале заполняемого участка и расхода в конце заполняемого участка в зависимости от времени.

На рис.2 приведен график изменения давления в начале ( х = 0 ) заполняемого участка трубопровода. Аналитическую зависимость изменения давления газа во времени принимаем .следующей:


Рн Lt

При

р(°>*) Рим J

р

г исп

при

P(0,t) ^ рисп ,

(16)


где ос


‘ UCfl


скорость подъема давления в сечении х - 0, Н/с*м£ испытательное давление, Н/м2: начальное давление газа, Н/м •


Граничное условие для температуры газа в сечении х = О принимаем в виде:


7


П0,ъ)~ r0~J)i [p0-P(0,t)J ,

(17)

2

где D- - коэффициент Джоуля—Томсона, м.К/Н;


Рв - давление гага в источнике испытательной среды,Н/м2 Т0 - температура газа в источнике испытательной среды,К.

Рис,2. График изменения давления в начале заполняемого трубопровода ( * = 0):

I- подъем давления до Роем ; 2-осмотр трубопровода; 3 - польем давления до Р^

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

(18)

G (д ,t) =0 .

2.7. В качестве исходных данных для расчета газолинами чесних параметров зададим следующие величины:

L - протяженность заполняемого участка, м;

])    -    диаметр заполняемого газопровода, м;

РКр - критическое давление газа, Н/м2;

-    критическую температуру гага, К;

*:

л-

К

R

р““

'Грр

Л,

Т0

-    эквивалентную гидравлическую шероховатость, и:

-    коэффициент динамической вязкости гаэа, Н*с/м2:

-    коэффициент теплопередачи от газа в грунт,Br/ir.K;

-    газовую постоянную, Дж/кг*К;

-    теплоемкость газа при постоянном давлении, Лк/кг*К;

-    коэффициент Джоуля-Томсона . м2*К/Н»

-    испытательное давление, Н/м2; начальное давление газа в заполняемом участке, Н/м2

-    температуру грунта, К;

-    давление газа в источнике испытательной среды, Н/м2

-    температуру газа в источнике испытательной среды,К;

-    скорость подъема давления в сечении х = 0, Н/м2.с.

Сечение я = 0 принимаем за начало заполняемого трубопровода, т.е. сечение & = 0 соответствует месту подключения источника испытательной средь*, а сечение х = L является концом заполняемого участка.

2.3. Дифференциальные уравнения (I) - (4), описывающие газодинамические процессы, происходящие внутри заполняемого участка газопровода, аппроксимируются замкнутой системой алгебраических соотношений, которая решается при помощи комбина -ций явной и неявной схем метода сеток*

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

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

2.9. Алгоритм расчета (прил.1) газодинамических парамет-

Э