Научная статья на тему 'Математическая модель и метод расчета нестационарных режимов в линейных участках магистральных газопроводов'

Математическая модель и метод расчета нестационарных режимов в линейных участках магистральных газопроводов Текст научной статьи по специальности «Математика»

CC BY
657
176
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по математике, автор научной работы — Тевяшев Андрей Дмитриевич, Гусарова Ирина Григорьевна, Каминская Анна Владимировна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Тевяшев Андрей Дмитриевич, Гусарова Ирина Григорьевна, Каминская Анна Владимировна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Mathematical model and method of non-stationary mode in linear pipeline portion of gas-main pipeline

In the article is taken comparative analysis of two implicit finite-difference schemes in modeling of non-stationary nonisothermic modes of natural gas transport in multistrand line sections gas-main pipeline. The formulated problem of calculation non-stationary nonisothermic modes of natural gas transport in multistrand line sections gas-main pipeline and presented the algorithms their calculation. For choice optimal finite-difference scheme take the calculation of model problem by following each of the templates of computational grids, both the pipeline section and linear section of multistrand line sections gas-main pipeline.

Текст научной работы на тему «Математическая модель и метод расчета нестационарных режимов в линейных участках магистральных газопроводов»

Матрица Pr(t) в (11) является трёхдиагональной.

Установлено, что замена

Y =

ґ X

Z Pr(t)e r

v r=0

Z

переводит уравнение (3) в дифференциальное уравнение

е hZ' =

(

\

Е Br(t)e

Vr=0 J

Z

Матрица Br(t) этой системы имеет блочно-диагональный вид (12).

Таким образом, исходная система, состоящая из n уравнений, может быть расщеплена на три подсистемы меньших порядков: k, l, m .

Аналогичные преобразования могут быть проведены и на промежутке [to + 5,ъ].

При достаточно малом 5 решения исследуемой системы можно заменить их линейными комбинациями.

4. Заключение

Научная новизна состоит в следующем. Рассмотрен процесс формирования активной смеси, в которой происходит образование порошковой массы. Описан

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

УДК 004.942:621.6:622.691

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

ТЕВЯШЕВ А.Д, ГУСАРОВА И.Г., КАМИНСКАЯА.В.

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

1. Введение

Современные методы исследования, проектирования и создания сложных систем неразрывно связаны с разработкой, реализацией на ПЭВМ и изучением их моделей. Исследование динамических свойств и характеристик таких систем часто проводится на их моделях, представленных системами дифференциаль-

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

Описана система дифференциальных уравнений, описывающих процесс формирования смеси. Установлено, что её можно расщепить на системы меньших порядков.

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

Литература: 1. Гора Н.Н. Уравнения процесса формирования многокомпонентных смесей // АСУ и приборы автоматики. 2006. Вып. 133. 2. Грановский М.Г., Лавров И.С., Смирнов О.В. Электрообработка жидкостей. Л.: Химия. 1976. 216с. 3. Ходаков Г.С. Физика измельчения. М.: Наука, 1972. 307с. 4. Вазов В. Асимптотические разложения решений обыкновенных дифференциальных уравнений. М.: Мир, 1968. 5. Sibuya Y. Sur un systeme des equations differentielles ordinaires lineaires a coefficients periodiques et contenant des parameters, J. Fac. Sci., Univ. Tokyo, (1), 7 (1954). Р. 229-2414. 6. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: В 10 т. Т.6: Гидродинамика. М.: Наука, 1986. 736с.

Поступила в редколлегию 04.04.2007

Рецензент: д-р техн. наук, проф. Кривуля Г.Ф.

Вовк Александр Владимирович, стажёр-исследователь кафедры прикладной математики ХНУРЭ. Научные интересы: теория вероятностей, случайные процессы. Адрес: Украина, 61100, Харьков, пр. Маршала Жукова, 45, кв. 16, тел. 716-16-88.

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

Однако не вызывает никаких сомнений, что разные задачи предъявляют разные требования к сеткам. Но лишь в отдельных случаях эти требования удается четко сформулировать. В остальных же случаях приходится сравнивать сетки на тестах. Зачастую удачный выбор параметров сетки повышает точность не менее, чем усовершенствование разностной схемы.

144

РИ, 2007, № 2

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

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

Проведем сравнительный анализ двух неявных конечно -разностных схем при моделировании нестационарных неизотермических режимов транспорта природного газа в многониточных линейных участках (ЛУ) магистральных газопроводов (МГ). Первая определена на четырехточечном шаблоне с разностным оператором второго порядка аппроксимации по пространственной переменной и первым порядком - по временной переменной, а вторая - на пятиточечном шаблоне с разностными операторами второго порядка аппроксимации по простр анственной и временной переменным. Сформулируем задачи расчета ННРТГ по ЛУ МГ и приведем алгоритмы их решений. Осуществим программную реализацию рассмотренных алгоритмов по двум выбранным схемам. Для выбор а оптимальной конечно -разностной схемы про -ведем тестрование программы как на участке трубопровода, так и на ЛУ многониточного МГ.

Для этого нужно решить следующие задачи: разработка алгоритма расчета ННРТГ по многониточным ЛУ МГ на основе математической модели и проведение сравнительного анализа выбраных конечно-разностных схем с помощью разработанного программного продукта.

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

Математическая модель нестационарных неизотермических режимов транспорта газа по участку трубопровода (УТ) магистрального газопровода рассмотрена в [1]. В настоящей статье приведена модель для УТ без промежуточных выкладок.

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

д W д t

W2 д P

+ (1 - aST——)— + 2 aST

p2 д x

W d W

P d x

+

+ PST

w|w|

p

+

g P dh aS T dx

(1)

5 P ^ ST 5 W 0 ----h aST---= 0,

d t d x

5T __W 5T „ T2 5 W

----ъ ayST----h aS(y -1)-----ь

d t P d x P d x

4K T WT dh

+---(Y-1) — (T- Tro) + g(y-1)-= 0,

D ' P ^ P dx

(2)

(3)

zgR где a=-e— S

P

_la Y_ Cp 2D ’ Cp - zgR

S - площадь попе-

речного сечения трубы, м2; Cp - удельная теплоемкость газа, Дж/кг- К ; z - коэффициент сжимаемости газа; R - газовая постоянная; W(x,t)T(x, t), P(x,t) -

удельный массовый расход, кг /(м 2с); температура, K ; давление (Па) газа, t, x - временная и пространственная переменные, l - коэффициент гидравлического сопротивления, D - диаметр трубы, м; K -коэффициент теплопередачи от трубы к грунту, Вт/м 2 • K ; Тф - температура грунта, K ; h - глубина залегания трубы, м; g - ускорение свободного падения, м/с2.

Модель структуры многониточного ЛУ МГ представляется ориентированным графом G(V,M), где V -множество узлов графа, М - множество дуг графа соответствующего УТ. Множество узлов V = V1 U V2 U V3 , где IV = v, = v1 - количество

входов ЛУ МГ; |У^| = V2 - количество промежуточных узлов; | V31 = V3 - количество выходов ЛУ МГ.

Систему уравнений (1) - (3) для j-го УТ запишем в матричной форме:

дф-. дф-.

—j + В—j = Фj є M

д t д x j

(4)

Wj 2 aTjSj—- j Jpj wj2 p2

1 - aTjSj 0

где Bj = aTjSj t2 a(Y“1)Sj-pj- pj 0 0 Wj aYTjSj —-' j j p. j

0

Ф.=

- PtjSj

WjWj

g Pj dh

a TjSj dx

4K T - — (Y - 1)tt(T- In) - g(Y-1)

T

Dj Pj

гр

TjWjdh

Pj dx

P

0

j (Wj,Pj,Tj).

Для того чтобы система уравнений (4) была разрешимой, необходимо задать граничные условия для узлов, соответствующих входам и выходам ЛУ МГ.

РИ, 2007, № 2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

145

Рассматриваются граничные узлы 1-го и 2-го типов. Для узлов 1-го типа задано давление, как функция времени, для узлов 2-го типа задан расход, как функция времени.

Граничные условия для m-го выходного и входного узлов имеют вид:

Gm (t) = Gm (t) (узел II типа)

или Руз (t) = Pm(t) (узел I типа), кроме того, на входах

задана температура поступающего газа T? (t) = Tm (t).

Условия согласования параметров газового потока для m-го промежуточного узла принимают следующий вид:

для расходов:

^ Gj(x++ ,t) = £Gi(x+ ,t),mєV2 (5)

j^V? ieV,m , ( )

для давлений:

Pj(x++ ,t)=pi(x+,t)j vmд є Vm, (6)

для температур:

Z ((Gj(x++ ,t))+-Tj(x++ ,t)) + E((Gi(x+ ,t))--Ti(x+ ,t)) =

jeVm

= T? • ( £ ((Gj(x++ ,t)) ++ Z(Gi(x + ,t))-),

jeVrn ieVm

кроме того,

если Gj (x++, t)<0, то Tm(x++ ,t) = TCP(t), jє Vm, если Gi (x+, t)>0, то Tm(x+,t) = T“(t), i є v?,

, |a,a > 0

где (a) =( A, (a) 10, a < 0

- a,a < 0 0,a > 0 ’

+

x

++

x

- на-

чальная и конечная координата соответствующего участка; V - множество узлов сети; V? ,Vm - множество индексов дуг, входящих и выходящих из m-го узла сети; G(x, t) T(x, t), P(x, t) - массовый расход, кг/с, давление (Па) и температура (к) для j-го участка; T? (t) - средняя температура вытекающего из m-го узла газа (к ).

Задается также начальное распределение расходов, давлений и температур для ЛУ: Wj (xj, 0) = Wj (xj), pj(xj. 0) = p"(xj). Tj(xj. °) = T"(xj). где

xj e[xj,x|+ ], vjeМ .

3. Метод решения

Для решения систем дифференциальных уравнений математической модели (4)-(7) использован метод сеток. Аппроксимация уравнений системы проведена по двум неявным конечно-разностным схемам. Первая определена на четырехточечном шаблоне с разностным оператором второго порядка аппроксимации по пространственной переменной и первым порядком

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

а) Пятиточечный шаблон

Для получения численного решения система (4) аппроксимируется разностными уравнениями с использованием неявной конечно-разностной схемы, определенной на пятиточечном шаблоне (см.рис. 1). Разобьем отрезок [x|,x|+ ] на Nj частей для каждогоj-го участка трубопровода. Разобьем временной интер-

T

вал [0, T] на k = 0, [—] временных слоев.

дф

Ж

k-1 -

k-2

4-

n-1

n+1

дф

d x

k

n

k

k

n

n

Рис. 1. Пятиточечный шаблон

В этом случае производные заменяются следующими соотношениями [2]:

дф д t

k = - 4фП-1+dk ~2

n 2 At

n = 0,N

(8)

дф

d x

k

dk -do

—x

kk Pn+1 _ Vn-1

n = 0,

2 —x

_ vN-1

—x

n = 1,N-1, n = N.

(9)

С учетом (8) - (9) система уравнений для j-го участка трубопровода принимает вид:

1 Bk ф^_1 +-3- dk ж^вП ^nk+1 =

2 Ax

2 At

2 Ax

j

= Фk +—фП 1 — ф£ 2,n = 1,Nj-1,

n At n 2 At n j

3 Jc 1 Bk ..k , 1 Bk ,.k _

2—T№" —xB0 n + A> я

= ф 0+± №k-2

0 At 2 At 0 ’

(10)

(11)

n

146

РИ, 2007, № 2

3 Лk , 1 Bk к_____1 Bk к _

2At ^Nj + AxBNj ^Nj -1

Выражения для невязок для любой внутренней точки

j-го участка трубопровода (n = 1, Nj -1) имеют вид:

- Фк _ 2 к-1 _ 1 к-2

ФNj At ^Ni 2 At ^Ni

j At ^Nj 2At ^Nj ' (12) vk’r,j = (—Wk,r,j-i-W*-1,j+-^Wnk-2,j) + —

Tn,i -ІЛ+ 11 Л+ 11 -ІЛ+ 11 -ІЛ

2 At

At

2 At

2 Ax

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2 At

Таким образом, получили систему нелинейных алгебраических уравнений, которая содержит 3(N+1) 2 aSTk,r,jWk.r,j aSTk,r, j (Wk,r,j )2

уравнений и 3(N+1) переменных. Данная система х[---------------“—W—(Wi+:[j _ W^k^rfj) + (1------—W;---------)х

решается методом Ньютона, путем линеаризации системы (10) - (12). Линейная система уравнений решается относительно векторов поправок к неизвестным

Pk,r,j

(Pnk>rJ)2

Wnk’r,j Wnk’r,j

k,r,j _ Pk,r,j)] ^ RCTk,r,j . n n

методом Гаусса с выбором главного элемента. Для ее ^Pn+i _ P„li )]+ PSTj1

решения необходимо вычислить векторы невязок и элементы матриц Якоби для каждого j-го участка трубопровода:

____________ g Pik,r,j dh

Pk,r,j aSTk,r,j dx’

k,r, j 3 k,r, j 1 t-, k,r, j k,r, j ,

¥0 = -AxB0 ^ +

^_^Bk,r,j k,rj _ Фk,r,j ^k-1, j + _^_ k-2J

+ AxBo 91 Ф0 At90 + 2 At90 ’

wk,r, j -_j fok,r, j +_!_ +

Vn “ 2 AxBn ^-1 + 2 At ^ +

,,,tr j - ( 3 Pk,r, j_^Pk-1, j і 1 Pk-2, j) ,

Vn,2 - (2AtPn AtPn + 2 AtPn ) +

+ ^ aSTk,r,j(Wnk5j -

k,r,j _ ( 3 Tk,r,j —2Tk_1,j и_—Tk_2,j) и__!i__х

Vn,3 - (2AtTn AtTn + 2AtTn )

aSTnk,r,j

-В^1 рЦ -фП^--^ +—ak“4n = 1

2 Ax

At

Vn

2 At

%

,Nj-l

хКг -1) j (Wnk;r1 j - w^"j)+yW^^ j (Tnk+r1 j - Tnk!r1 j)] +

izk;r,j =—^N^j +—bN" J ^k:r,j

^Nj 2Ak Nj Ax Nj ^Nj

1 Bk,r,j ™k,r,j _Фk,r,J_ 2 ^k-1,J , 1 k-2,j

-----B^ ’-1 - ф1^

AxBnj ^V1 Ф Ni

<P-

j At ' Nj 2 At ' Nj

V-

Nj

4K Tk,r,j k ■ Tk,r,jWk,r,j dh

+1T(Y - «Tin1"-ТФ)+g<Y-')^L?i;ij.

nn Выражения для определения невязок в граничной точке N j имеют вид:

Запишем выражения для невязок Wj = ( 3 Wk,r,j__2Wk-1,j + 1 Wk-2,j4^J_ Vn = (vn,1, Vn,2,¥11,3) в развернутой форме. Ni,1 (т A Ni A Ni 2 A Ni

Выражения для определения невязок в точке 0 имеют вид:

Vin j = (— W0k’r, j - — W0k_1, j + — W0k-2, j) +

У0Д 2 At 0 At 0 2 At 0

2 aSTkT’r,jWk,r,j

,, aSTk,r,j(Wk,r,j)2

Nj Nj k_ri kri aS1Nj (WNj )

j------(Wk’r,j - Wk,r,j ) + (1----------------------У.

Г (WnJ WnJ-1) + (1 (Pk,r,j)2 *

Nj (PNj )

P

,Wt>.i Wk«.+(1 saq^^'-1)2

-k-ri V 1 0 3 V (Pk,j)2

+— r-

AxL PkrJ

P0

)х х(р;т"pNh)]+ps

Tk,r’j Wk,r’j

TNj WNj

Wk,r, j WNj

P

k,j,l

N

Pk,r, j

g Nj dh aS t k,r’j dx

Nj

х (P1k’r,j - P0k’r,j)] + PS

Tk,r’j Wk^r’j 10 vv0

W0k’r,j

g PQ-^dh,

aS Tk,r’j dx 0

Pkr,j

0

Wk’r’j - ( 3 Pk,r, j __^Pk 1,j j_^Pk 2,j) +

V0,2 ■ (2AtP0 AtP0 + 2AtP0 ) +

+—aST0k,r’j (W1k,r’j - W0k,r’j),

wk,r,j = ^^^Pk,r,j —^Pk 1,j^^Pk 2,j) +

V Nj’2 (2 AtPkJ AtPkJ + 2 AtPkJ ) +

+ — aST,k’r’j (W^^r’j - W^^r’j,),

Ax NJ Nj Nj -1

k,r, j _ ( 3 Tk,rJ j _ 2 Tk—1j j ^ 1 Tk_2>j

aST

aST,

k,r, j 0

Ax 0

,,,k,r,j _ ( 3 Tk,r,j — 2 Tk~1,j + ^^_Tk—2’j) +_

V0,3 ■ (2At T AtT0 + 2At T ) + AxPk,r,j^

х: [(Y -1) T0k-r, j (W1k’r’j - W0k’r,j) + yW0k’r’j (T1k’r’j - T^j)] +

-1---T

VNi,3 ^a,tN| At T

j,3 '2At NJ At NJ 2At NJ

^ ^,j) TN )

k,r, j

Nj

AxP

k,r, j j

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

т5і

N

х [(y -1) T-Nj"’j (W^j - wkk,jr:_-i)+YWjj’j (Tjj>j - TNjr1-1)]+

T

k,r, j

Nj

Tk’r’jwk’r’j

4K N| kr ■ Nj N| dh

+ 1T(Y - l)‘Pktr(TN' Тгр) - g(Y-'>

Nj

Nj

4K Tk,r,j Tk,r,j Wk,r,j

+ 4k(Y -1) T— (Tk’r’j - Тр) + g(Y-1) 0 kr0

D ■ >k,r,j ^ - Pk,r,j dx'

РИ, 2007, №

02

147

д у

дф

k, r, j

в начальной

J о

Элементы матрицы Якоби точке 0 имеют вид (индексы k, r, j опущены):

5¥<u _ 3 .2^_ 1Щ)_(p__p^Wo, +

5у„д _ aSW^[i(wn+i_ wn-1) _ Pn+1 Pn-1

d T

1 AxP„ - ' Pn

, PSWn|Wn| gPi dh

Wn] +

d W 2 At AxP,

0

2 PST0W0I

P0

d ¥0,1 = - 2 aST0W0(W1 - W0) _ +

dP Ax[ P2 ( +

+ aST0W02(P0- 2P1) _ PSTqWqIwQ + _ g dh

Pn

d ¥ n,2

d W

5 ¥n,2 _

5 P "

9 ¥ n,2 aS

aSTn2 dx

= 0,

3

2 At’

Pq3

P02

aST0 dx

d T 2 Ax

9 ¥ n,3 _ yaS Tn

(Wn+1 - Wn_1),

5 ¥0,1 aWoS [2(W W) P1 - pqW] , -[2(W1 - Wo)----------------------Wo] +

5 T

AxP0 + PSw0|Wq|

+ Pq

9 ¥ 0,2 _ d W ”

9 ¥ 0,2

P0

g Po dh aS T02 dx

aT0S

d W 9 ¥n,3 - aSTn 5 P 2Axp2

2 AxPn(Tn« - T-1'+ З^1^-

[(Y-1)T0(W..1-W,,-1)+YW,(T„+1-T„_1)]-

4K T T W dh

- — (y- 1)-2(Tn-Trp)-g(Y-1) n n

D

Pn2 dx

Ax

3

d ¥ n,3 5 T

3

aS

2 At 2 AxPn

-[2(Y-1)Tn(Wn+1 - Wn_1) +

d P 2 At ’

9 ¥0,2 _ aS(W1 - W0)

d T

Ax

9¥0,3 aSTo [(1 T | i ( 1)Tq dh

W=AP0[(1 -2 Y)To+YT1]+g(Y~ «iOd?

= ^aST°- fo-1) To (W1 - Wo)+YWo (T1 - To J -

"i ^

d P

AxP

0

_ iDK(Y- 1)To(To 2Trp} - g(Y- 1)TqWq dh

P02

P02 dx

4K 2Tn - Trp W dh

+ YWn(Tn+1-Tn_1)] +—(y -1) np rp + g(y -1)^— D Pn Pn dx

б) Четырехточечный шаблон

Для получения численного решения система (4) аппроксимируется разностными уравнениями с использованием неявной конечно-разностной схемы, определенной на четырехточечном шаблоне (см. рис. 2). Разобьем отрезок [x|,xj+ ] на Nj частей для каждого j-го участка трубопровода. Разобьем временной

T

интервал [0,T] на k = 0, ["Af] временных слоев.

9¥03 3 aS

= ТГ7 + t-^[2(y -1)To(W1 - Wo) + yWq(T1 - 2Tq)] + и T 2 At AxPq

4K 142T0_^'гр , W0 dh

D P02 P02 dx

дф d t

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

k

n

k t

Элементы матрицы Якоби

d ¥

дф

"|k,r, j

для внутренних

k-1

n-1

—I---

n+1

точек (n = 1, Nj -1) имеют вид:

5¥n,i = +TSTn[(Wn+1 - Wn_1) - Pn+1~ Pn-1

дф d x

5 W 2 At AxP,

Wn] +

+ 2 PS

Tn|Wn|

Pn

Рис. 2. Четырехточечный шаблон

Тогда производные заменяются следующими соотношениями [2]:

3¥„л = oSVW. +1 _ Wn_1)- Pn+1р Pn-1 Wn]-

2

дф

d t

k mk k-1

_ Pn ~Pn

At

n = 0, N

(13)

5 P

AxPn2

Pn

PSTnWn|Wn| + g dh

aSTn dx

дф

d x

148

~Po

Ax

kk Pn+1 ~ Pn-1

n = 0,

2 Ax kk T’N - ^N-1

Ax

n = 1,N-1, n = N.

(14)

РИ, 2007, № 2

0

k

n

n

n

n

k

P

n

С учетом (13) - (14) уравнение (4) для для j-го участка трубопровода записывается следующим образом:

1 1 Bk k , к фк , 1 к-1

— Фо в0 <Р0 +—в0 91 - ф0 + — Фо , (15)

At

Ax

Ax

At

lAxBn ^”k->+ At+ 2AXBk 'и*1=

1

= ФП+A: %k'\n = 1,Nj-1, At J

—*kT +—bn ®n

At Nj Ax Nj Nj

(16)

—bN ®N і = фN +—<pN

Ax NJ Nj 1 Nj At N

к-1

At^j • (17)

Аналогичные выкладки, проведенные для пятиточечного шаблона, справедливы и для разностной схемы, определенной на четырехточечном шаблоне.

Как и в пункте а) система уравнений (15) - (17) решена методом Ньютона. В результате линеаризации получаем линейную систему, которая была решена методом Гаусса с выбором главного элемента, относительно векторов поправок к неизвестным.

4. Результаты расчетов

Приведем сопоставление результатов расчетов сети по двум конечно-разностным схемам.

Рассматриваемая сеть, состоящая из ЛУ МГ, имеет следующие характеристики: суммарная длина всех труб - 170 км, внутренние диаметры труб - 1400 мм, толщина стенок h=10 мм, шероховатость труб -0,03мм, удельная теплоемкость газа Cp = 0.655952 ккал/(кг0С), относительная плотность газа по воздуху A =0,604707, температура грунта на глубине заложения газопровода trp=10°C, коэффициент теплопередачи от газа к грунту K = 1,4 ккал/ (м2ч0С). Расчетная схема представлена на рис. 3.

® о ©

©'

19

©

Входы © 10 3 6 © © 11

ЛУ 0 2© 4708 9 12 ©■

©

14 © 16 ©

©

. Выходы ЛУ

©

5

20

17

115км

35км

20 км

Рис. 3. Расчетная схема Л У МГ

Задача решалась на пространственно-временной сетке с шагами Ax = 20 км и At = 30 мин на период времени 24ч. В качестве начального условия принимаем стационарное течение газа по рассматриваемой ЛЧ с давлением на входах равным 8,398МПа, температурой - 40 0 C и суммарным коммерческим расходом -306,798 млн.м3 /сут. Ниже приводятся графики зависимости параметров газового потока от времени для разных шаблонов расчетов.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

На рис. 4 - 6 представлен переходной процесс, моделирующий заданное по определенному ступенчатому временному закону увеличение и последующее восстановление выходного расхода в 16 и 14 узлах линейной части и постоянным расходом в 17, 10, 18, 15 узлах рассматриваемой сети. Граничные условия представлены в таблице.

Узлы сети Значения в узлах сети

1 P(t) = 8,398 МПа T(t) = 400C

2 P(t) = 8,398 МПа T(t) = 400C

3 P(t) = 8,398 МПа T(t) = 400C

10 q10(t) = 15 млн м3 в сут. t > 0мин.

15 q15(t)=10 млн м3 всут, t>0мин.

16 q6(t) = ■ 15млн м3всутt < 180мин. 50млн м3всут.,180< t < 390мин. 15млнм3 всутt>390мин.

17 9 (15млнм3 всут.Д< 120мин. q9(t)=1 [25 млнм3 всут.Д> 120мин

18 qu(t) = 10млн м3 в сут., t > 0мин.

5. Выводы

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

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

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

РИ, 2007, № 2

149

--------2 узел (четырехточечный шаблон)

5 узел (четырехточечный шаблон) узел (четырехточечный шаблон)

--------15 узел (четырехточечный швблон)

........2 узел (пятиточечный швблон)

- - - '5 узел (пятиточечный шаблон)

™ ™ '16 узел (пятиточе^ый шаблон)

........15 узел (пятиточе^ый шаблон)

200 400 600 800 1000 1200 1400

Время, мин.

5

га

d

Рис. 4. График зависимости расхода от времени в узлах сети ЛУ МГ

2 узел 5 узел 16 узел

15 узел

--------2 узел (четырехточечный шаблон)

5 узел (четырехточечный шаблон) ^^^^*16 узел (четырехточечный шаблон) --------15 узел (четырехточечный шаблон)

■ - ----2 узел (пятиточечный шаблон)

- - - -5 узел (пятиточечный шаблон)

■ ■ '16 узел (пятиточечный шаблон)

- - ----15 узел (пятиточечный шаблон)

600 800 Время, мин.

Рис. 5. График зависимости давления от времени в узлах сети ЛУ МГ

45

40

со 35

30

25

20

2 узел в 5 узел

-2 узел (четырехточечный шаблон) -5 узел (четырехточечный шаблон) •16 узел (четырехточечный шаблон) -15 у зел (четырехточечный шаблон) •2 узел (пятиточечный шаблон)

■5 узел (пятиточечный шаблон)

■ 16 узел (пятиточечный шаблон)

■ 15 узел (пятиточечный шаблон)

200 400

600 800 Время, мин.

1000 1200 1400

Рис. 6. График зависимости температуры от времени в узлах сети ЛУ МГ

0

86

82

80

78

0

200

400

1000

1200

1400

16 узел

15 узел

0

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

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

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

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

Литература: 1. Тевяшев А.Д., Гусарова И.Г., Чуркина А.В. Эффективный метод и алгоритм расчета нестационарных неизотермических режимов транспорта газа в газотранспортной сети произвольной структуры // Восточно-Европейский журнал передовых технологий. 2006. 2/3 (20). С. 45-52. 2. Годунов С.К., Рябенький В. С. Разностные схемы (введение в теорию). Учеб. пособие. М: Наука, 1973. 400с.

Поступила в редколлегию 24.05.2007

Рецензент: д-р техн. наук, проф. Левыкин В.М.

Тевяшев Андрей Дмитриевич, академик УНГА, д-р техн. наук, профессор, зав. кафедрой прикладной математики ХНУРЭ. Научные интересы: теория стохастических моделей. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел.: (057) 702-14-36.

Гусарова Ирина Григорьевна, канд. техн. наук, доцент кафедры прикладной математики ХНУРЭ. Научные интересы: математическое моделирование и управление систем с распределенными параметрами. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел.: (057) 702-14-36.

Каминская Анна Владимировна, аспирантка кафедры прикладной математики ХНУРЭ. Научные интересы: численное моделирование трубопроводных систем. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел.: (057) 702-14-36.

РИ, 2007, № 2

150

i Надоели баннеры? Вы всегда можете отключить рекламу.