Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ФИЗИЧЕСКИЕ СВОЙСТВА НЕФТЕЙ

.pdf
Скачиваний:
117
Добавлен:
20.08.2019
Размер:
2.32 Mб
Скачать

vk.com/club152685050 | vk.com/id446425943

325

Отсюда находим:

pн. = 3,82 +5,78 105 21000 2 105 6,0 МПа,350

то есть найденное ранее значение давления pн. практически

не изменилось. Следовательно, степень сжатия ε, которую должны обеспечивать нагнетатели 370-18-1, равна

6,04,7 1,28.

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

тателя, учитывая их параллельное соединение:

 

 

 

 

 

 

 

R =

Rвозд.

=

 

287,1

441,7 Дж/(кг К);

 

 

 

 

 

 

 

 

 

0,65

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4,7

 

285

3,668

 

 

 

 

 

 

 

 

 

 

Zв. =

10,4273

 

 

 

 

 

 

 

 

 

0,895 ;

 

 

 

 

 

 

4,75

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

195

 

 

 

 

 

 

 

 

 

 

 

ρст. = ρвозд. ∆ =1,204 0,65 0,783 кг/м3;

 

 

 

 

 

 

 

ρв.

=

 

 

pв.

 

=

 

 

 

4,7 106

 

 

 

41,716 кг/м3;

 

 

 

 

 

 

Zв.RTв.

0,895 441,7 285

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qв. = Qк.ρст.

ρв. =

[(21000 2)

 

350] 106

 

0,783

3

/мин.

 

 

 

 

 

 

 

 

24

60

 

 

41,716

391м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Затем определяем приведенные параметры режима ра-

боты центробежного нагнетателя:

 

 

 

 

 

 

 

n

 

 

=

 

n

 

 

 

Zпр.Rпр.Tпр.

=

 

n

 

0,90 490 288

1,062

n

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0

Zв.RTв.

 

n0

 

0,895 441,7 285

n0

n0

 

пр.

 

 

 

 

 

 

 

 

 

 

 

 

 

(Qв. )

 

 

= Qв.

n0

 

= 391

n0

м3/мин.

 

 

 

 

 

 

 

 

 

 

пр.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку степень сжатия ε уже известна и равна 1,28 , то необходимо, используя приведенные характеристики 370-18-1, рис. 1.13, подобрать значение nn0 так, чтобы точка с координа-

vk.com/club152685050 | vk.com/id446425943

326

тами (Qв. )пр. = 391(nn0 ) и ε =1,28 лежала на характеристи-

ке(nn0 )пр. =1,062 nn0 . Подбор осуществляем методом после-

довательных приближений.

1) Полагаем (nn0 )пр. =1,0 n n0 =1,0 1,062 0,942 ;

(Qв. )

пр.

= 391 0,942 415 м3/мин ε 1,25 (см.

 

 

рис.1.13), что меньше необходимого значения 1,28. Следовательно, (nn0 )пр. нужно увеличить.

2) Полагаем (nn0 )пр. =1,05 n n0 =1,05 1,062 0,989 ;

(Qв. )

пр.

= 391 0,989 395 м3/мин ε 1,28 (см. рис.

 

 

1.13), следовательно, решение найдено.

Имеем: n = 0,989 n0 = 0,989 4800 4750 об/мин.

234. Определим сначала степень ε1 сжатия газа первым

нагнетателем. Имеем:

 

 

 

 

 

 

 

 

R =

8314

 

489 Дж/(кг K);

 

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3,5

 

283 3,668

 

 

 

Zв.1

=10,4273

 

 

 

 

 

 

 

0,951;

4,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

170

 

 

 

 

 

ρв.1

=

 

pв.1

 

=

 

 

 

3,5 106

 

 

26,595 ;

 

Zв.1RTв.1

 

0,951 489

283

 

 

 

 

 

 

 

 

 

 

 

 

n

 

=

 

5300

 

 

0,91 490 288

 

 

0,85 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6150

 

0,951 489 283

n0

пр.

 

 

 

 

(Qв. )пр. = 250 53006150 290 м3/мин.

Используя характеристики Н-300-1,23, представленные на рис. 1.14, находим: ε1 1,19 .

Рассчитываем давление, температуру, см. (132), и расход газа на входе второго нагнетателя:

vk.com/club152685050 | vk.com/id446425943

327

pв.2

 

= pн.1 = ε1 pв.1 =1,19 3,5 = 4,165 МПа;

 

Tв.2

 

 

 

 

 

 

m1

 

 

 

 

 

 

1,271

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ε1

=1,19

1,27

 

1,0377 ;

T

 

= 283 1,0377 294 К;

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tв.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4,165

 

 

294 3,668

 

 

 

Zв.2

 

=10,4273

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,949 ;

 

4,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

170

 

 

 

 

 

ρв.2

=

 

 

 

 

pв.2

 

=

 

 

 

 

4,165 106

 

 

 

30,528кг/м3;

Zв.2RTв.2

 

0,949 489 294

 

 

 

 

 

 

 

 

 

 

 

Qв.2

=

 

ρв.1

 

Qв.1

=

 

26,595

 

250 218 м3/мин;

 

ρв.2

30,528

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

=

5700

 

 

 

 

0,91 490 288

0,90 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6150

 

 

 

0,949 489 294

n0

 

пр.

 

 

 

 

 

 

 

 

 

 

(Qв. )

 

 

= 218

6150

235 м3/мин.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр.

 

 

 

 

 

 

 

 

 

 

 

5700

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя характеристики Н-300-1,23, представленные на

рис. 1.14, находим: ε2

1,23.

 

 

 

 

 

 

 

 

Степень

 

ε

 

 

 

сжатия газа системой двух нагнетателей равна

произведению

ε1 ε2

степеней сжатия отдельных нагнетателей.

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε = ε1ε2 =1,19 1,23 1,46 .

 

 

 

 

 

Используя формулу (132), находим:

 

 

 

Tн.2

 

 

 

 

 

 

 

m1

 

 

 

 

 

 

 

1,271

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 294 1,045 307 К.

 

 

 

= ε2

 

=1,23 1,27

1,045 T

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tв.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

235. Определим сначала степень ε1 сжатия газа первым

нагнетателем. Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

R =

 

 

Rвозд.

=

 

287,1

463,1Дж/(кг K);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,62

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zв.1

 

=10,4273

3,2

 

288 3,668

0,924 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

vk.com/club152685050 | vk.com/id446425943

 

 

 

 

 

 

 

 

328

 

ρв.1

=

 

3,2 106

 

25,967 кг/м3;

0,924 463,1 288

 

 

 

 

 

 

n

 

=

4800

 

0,90

490 288

1,015;

 

 

 

 

 

 

 

 

4800

0,924

463,1 288

n0

пр.

 

 

(Qв.1 )пр. = 500 48004800 = 500 м3/мин.

Используя характеристики нагнетателя 370-18-1, представленные на рис. 1.13, находим: ε1 1,235.

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

pв.2

 

= pн.1 = ε1 pв.1 =1,235 3,2 = 3,952 МПа;

 

 

 

Tв.2

 

 

 

 

 

m1

 

 

 

 

 

1,251

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 288 1,043 300,4 К;

 

 

= ε1

=1,235 1,25

 

1,043

; T

 

 

 

m

 

 

 

 

 

 

 

 

Tв.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zв.2

 

=10,4273

3,952

 

300,4 3,668

 

 

 

 

4,7

 

 

 

 

 

 

0,919 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

ρв.2

=

 

 

 

 

 

3,952 106

 

 

30,912 кг/м3;

 

 

0,919 463,1 300,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qв.2

=

 

ρв.1

 

Qв.1

=

 

25,967

 

500 420 м3/мин;

 

 

 

ρв.2

 

30,912

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

=

n2

 

 

 

0,90 490 288

0,997

n2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0

0,919 463,1 300,4

n0

n0

 

пр.

 

 

 

 

 

 

(Qв.2 )

 

 

=

 

 

 

420

 

м3/мин, где n

0

= 4800 об/мин.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр.

 

 

 

n2 n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

 

суммарная

степень

сжатия ε равна ε1ε2 ,

тоε2 = ε ε1 =1,5 1,235 1,21,

то есть степень сжатия ε2 во

втором нагнетателе известна. Теперь необходимо, используя приведенные характеристики нагнетателя 370-18-1, рис. 1.13, подобрать значение n2 n0 так, чтобы точка с координатами

vk.com/club152685050 | vk.com/id446425943

 

 

 

 

329

 

 

 

 

(Qв.2 )пр. = 420 (n2 n0 )

и ε2

=1,21 лежала

на

характеристи-

ке(n n0 )пр. = 0,997 n2

n0 . Подбор осуществляем методом по-

следовательных приближений.

 

 

 

 

 

Полагаем (n2 n0 )пр. = 0,95 n2 n0 =0,95 0,977 0,972;

тогда (Qв.2 )

пр.

= 420 0,972 432 м3/мин

 

ε

2

1,21 (см.

 

 

 

 

 

 

 

рис. 1.13). Следовательно,

(n n0 )пр. найдено уже в первом

приближении и в последующих приближениях необходи-

мости

 

нет.

 

 

Имеем: n2 = 0,972 n0 = 0,972 4800 4670

об/мин.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определяем мощность N(1)влп. на валу привода первого

ЦБН. Согласно формуле (130), имеем:

 

 

 

 

 

 

 

(1)

 

n1

3

 

N

 

3

 

N

 

N

 

N

 

= ρв.1

 

 

 

 

 

 

= 25,967 (1,0)

 

 

 

25,967

 

 

,

 

 

 

 

 

 

 

 

n0

 

 

ρв. пр.1

 

 

 

ρв. пр.

 

ρв. пр.

 

где значение

(N ρв. )пр.

берется согласно приведенной ха-

рактеристике 370-18-1 при приведенном расходе 500 м3/мин,

рис.

1.13:

(N ρв. )

260 кВт/(кг/м3).

 

Отсюда

находим:

 

 

 

 

 

 

 

пр.

 

 

 

 

 

 

 

 

 

 

N(1) = 25,967 260 6750

кВт и мощность

N(1)влп.

на валу

привода: N(1)влп. = 6750 +100 = 6850 кВт.

 

 

 

 

 

 

 

Определяем мощность N(2)влп. на валу привода второго

ЦБН. Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

n2

3

 

N

 

3

 

N

 

 

N

 

N

 

= ρв.2

 

 

 

 

 

= 30,912 (0,972)

 

 

 

28,4

 

 

,

 

 

 

 

 

 

 

n0

 

 

ρв. пр.

 

 

 

ρв. пр.

 

ρв. пр.

 

где значение

(N ρв. )пр.

берется согласно приведенной ха-

рактеристике 370-18-1, рис. 1.13 при приведенном расходе 432

м3/мин: (N ρв. )

пр.

250 кВт/(кг/м3). Отсюда находим:

 

 

N(2) = 28,4 250 =

7100 кВт и мощность Nвлп. на валу приво-

да: N(2)влп. = 7100

+100 = 7200 кВт.

vk.com/club152685050 | vk.com/id446425943

330

Таким образом, суммарная мощность, потребляемая обоими нагнетателями составляет 7200 + 6850 =14050 кВт или 14,05 МВт.

236. Согласно формуле (137), адиабатическая скорость звука c в газе (то есть волны быстрых колебаний, происходящих без теплообмена) определяется выражением:

c = γ ZRT , гдеγ = CP

CV = 2500 2030 1,232 .

Далее находим: R=

8314

 

467 Дж/(кг К),

17,8

 

 

 

 

 

 

 

4,3

 

288

3,668

Z =10,4273

 

 

 

 

 

0,910 ,

4,8

194

 

 

 

 

c = γ ZRT = 1,232 0,91 467 288 388м/с.

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

Рассчитываем параметры газа:

 

R =

 

Rвозд.

=

287,1 = 486,6 Дж/(кг К),

 

 

 

 

0,59

 

 

 

 

 

 

 

 

 

5,2

 

303

3,668

 

Z =

10,4273

 

 

 

 

0,884 ,

4,55

205

 

 

 

 

 

 

 

 

γ = CP CV = 2400 1913 =1,255,

c =

 

γ ZRT = 1,255 0,884 486,6 303 404 м/с.

238. Определим сначала параметрc γ = ZRT , назы-

ваемый изотермической скоростью звука в газе: c γ = ZRT = 0,9 490 288 356,4 м/с.

Далее по формуле (116) вычисляем коэффициент λ гидравлического сопротивления:

vk.com/club152685050 | vk.com/id446425943

 

 

 

331

 

 

2k

0,2

2 0,05

0,2

λ = 0,067

d

 

= 0,067

 

0,011.

 

 

996

 

 

По формулам (124) и (122) рассчитываем коммерческие расходы газа до (индекс 0) и после (индекс 1) отключения

ГПА. Имеем:

 

 

 

 

 

 

 

 

 

 

 

A =17,02 106

 

9962,6

 

 

85,9 ,

 

 

 

 

288 0,59

 

 

 

 

 

 

0,9

 

 

 

 

(Qк. )

= 85,9

5,12 3,82

26,67 млн.м3/сутки,

 

 

0

 

 

120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Qк. )

= 85,9

4,52 3,82

18,90 млн.м3/сутки.

 

 

1

 

 

120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим среднюю скорость

v0 газа на участке газо-

провода до отключения ГПА:

 

 

ρст. (Qк. )0

 

 

ρср.v0 S = ρст. (Qк. ) v0

=

 

 

 

.

pср.0

(ZRTср. )

πd2

4

 

 

 

0

 

 

 

 

Поскольку

ρст. = ρвозд. ∆ =1,204 0,59 0,710

кг/м3, а pср. ,

согласно (112) составляет:

 

 

 

 

 

 

 

 

2

 

3,82

 

 

 

 

 

 

 

 

 

pср.0 =

3

5,1+

 

 

 

= 4,48 МПа,

 

 

 

 

 

 

 

 

 

 

 

5,1+3,8

 

 

 

 

 

 

 

находим:

 

 

0,710 26,67 106 (24 3600)

 

v0 =

 

 

 

8,0 м/с.

4,48 106

(0,9 490 288) 3,14 0,9962 4

Аналогично вычисляем скорость v1 газа на участке га-

зопровода после отключения ГПА:

 

v1 =

 

 

 

 

ρст. (Qк. )1

 

,

 

pср.1

(ZRTср. )πd2

4

 

 

 

 

2

 

 

 

3,82

 

 

pср.1 =

 

4,5

+

 

= 4,15 МПа,

 

3

4,5 +3,8

 

 

 

 

 

 

 

 

 

vk.com/club152685050 | vk.com/id446425943

 

332

 

v1 =

0,710 18,9 106 (24 3600)

6,0 м/с.

4,15 106 (0,9 490 288) 3,14 0,9962 4

Таким образом, из-за отключения одного из ГПА средняя скорость газа на участке газопровода уменьшилась с 8,0 до 6,0 м/с, поэтому среднюю скорость vср. в рассматривае-

мом переходном процессе можно принять равной 7,0 м/c.

Далее имеем:

 

 

 

 

 

 

 

a

2

=

(c γ)2 d

=

(356,4)2 0,996

1,64

10

6

2

/с.

 

λ vср.

0,011 7

 

м

 

 

 

 

 

 

 

 

 

239. Решение уравнения (140) типа теплопроводности с начальным условием qк. (x,0)= 0 и краевыми условиями:

qк. (0,t)= q = const. и q 0 при x → ∞, имеет на полуоси x > 0 , как известно, см. формулы (71) и (72), следующий вид:

 

2

2

 

x

 

 

2

 

 

 

q

 

(x,t)= q

erf

 

 

 

 

= q

 

1

к.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

t

 

 

 

 

2

x 2a

t

 

−ς2

 

 

 

e

 

dς .

π

 

 

0

 

 

 

 

Понимая под q2 изменение квадрата коммерческого расхода газа, имеем:

[Qк.0 + ∆Qк. (L,t)]2 Qк.0

2

2

L 2a

t

−ς2

 

2

2

=1

 

e

 

dς.

π

 

[Qк.0 + ∆Qк.0 ]

Qк.0

 

0

 

 

 

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

[Qк.0 +0,01 (0,25 Qк2.0 )]2 Qк.02 0,0089 . [Qк.0 +0,25 Qк.0 ] Qк.02

Отсюда получаем уравнение

2150000

e−ς2 dς = 0,9911

π04 1,64 106 t

для определения искомого момента t времени.

vk.com/club152685050 | vk.com/id446425943

333

Используя таблицы значений интеграла вероятностей [11], находим:

2 150000

4 1,64 106 t

 

−ς2

 

 

150000

 

 

 

e

 

dς = 0,955

при

 

1,85 .

π

 

4 1,64 106 t

0

 

 

 

 

 

или t 1002 с (16,7 мин).

240. Для описания переходного процесса на участке газопровода используем уравнение (139), а также выражения (141) и

(142):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 (x,t)

 

2

 

2p2 (x,t)

 

 

& 2

 

2

 

 

γd

p2

 

 

 

t

= a

 

 

x2

 

,

(1),

M

= −S

 

 

λc2

x , (2)

где

 

 

массовый

расход

газа;

 

a

2

 

2

d

(γλvср. );

 

M

 

 

= c

 

2

 

 

 

 

 

 

 

′′

 

 

′′

 

 

 

 

 

c

γ = Zср.RTср. ;

 

 

vср. = 0,5(vср. + vср. );

vср. ,vср.

средние ско-

 

 

 

рости

газа

в

старом

и

 

новом

стационарных

 

режимах;

S= πd2 4 площадь сечения трубопровода.

Вновом стационарном режиме

 

2

2

 

2

2

& 2

 

2

 

2

pк.

2

 

 

pн.

pк.

 

γdS

 

pн.

 

p

 

(x)= pн.

 

L

x , M

=

 

 

 

L

,

 

 

λc2

 

поэтому решение задачи можно представить в виде:

p 2 (x,t)= pн.2 pн.2 pк.2 x +(pн.2 pн.2 ) Ф(x,t), (3) L

где Ф(x,t)неизвестная безразмерная функция, характеризую-

щая отклонение давления в сечении x газопровода в момент времени t от значения в новом стационарном режиме. Очевидно, что

функция Ф(x,t), как и функция p 2 (x,t) удовлетворяет основному уравнению (1).

Кроме того, функция Ф(x,t) удовлетворяет следующим

краевым и начальному условиям:

Вначале участка (при x = 0 ): Ф(0,t)= 0 для всех t > 0 .

Вконце участка (при x = L ): Ф(L,t)= 0 для всех t > 0 .

vk.com/club152685050 | vk.com/id446425943

334

При t = 0 (начальное условие) на участке 0 x L газопровода существовал стационарный режим с давлениями pн. в

начале участка и pк. в его конце. Учитывая это обстоятельство и представление (3), полагаем:

p 2 (x,t)= pн.2

pн.2 pк.2

x +(pн.2 pн.2 ) Ф(x,0)pн.2

pн.

2 pк.2

x,

L

 

L

 

 

 

 

откуда имеем: Ф(x,0)=1xL .

Решение уравнения (1) с полученными краевыми и началь-

ным условиями методом разделения переменных. Согласно это-

му методу, функция Ф(x,t) представляется в виде ряда, каждый

член которого есть произведение функции, зависящей только от t, на функцию, зависящую только от x:

Ф(x,t)= Θn (t) Rn (x). (4)

n =1

Потребуем, чтобы каждый член этого ряда в отдельности удовлетворял уравнению (1). Получим:

 

dΘn (t)

 

Rn (x)= a2

Θn (t)

d2Rn (x)

 

 

 

dt

dx2

 

 

 

 

 

 

 

 

 

 

или, разделив обе части уравнения на произведение Θn (t)Rn (x):

1

 

1

 

 

dΘn (t)

=

1

 

d2Rn (x)

.

 

a2

Θn (t)

 

dt

Rn (x)

dx2

 

 

 

 

 

 

Левая часть этого уравнения зависит только от t, правая – только от x. Такое может быть только в случае, если каждая из этих частей есть константа. Имеем:

1

 

1 dΘn (t)

 

1

 

d2Rn (x)

2

 

 

 

 

 

=

 

 

 

= −µ n = const.

a2

Θn (t) dt

Rn (x)

dx2

 

 

 

 

Для существования решения эта константа (см. учебники по уравнениям математической физики) должна быть отрицатель-

ной, поэтому мы обозначили ее (−µ2n ). Далее имеем: