Файл: Подводные и подземные взрывы сб. ст.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 84

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

МОЩНЫЕ ПОДВОДНЫЕ ВЗРЫВЫ

13

Условия на ударной волне даются соотношениями Рэн­ кина — Погонно:

 

p0u = Ps(U — us),

(6а)

 

Pot/2 + Ро= РЛU — us)24- ft,

(66)

^+£+*°=-^Чг }1 + £ + ‘!--

<вв)

Здесь индекс

0 относится

к невозмущенным

условиям,

а индексом s

обозначается

состояние на скачке; U

скорость ударной волны. На самой ранней стадии взры­ ва давление на ударной волне весьма высоко, поэтому мы сделаем предположение о сильной ударной волне [6], т. е. давление и энергия в невозмущенной среде будут полагаться равными нулю. Таким образом, в уравнениях (66) и (6в) будем пренебрегать членами, содержащими Ро и во. Тогда единственными размерными постоянными величинами, которые появляются в задаче о гомотермическом взрыве, будут невозмущенная плотность р0 и энергия взрыва Ео.

Задача об одномерном неустановившемся течении при наличии только двух независимых размерных по­ стоянных величин является автомодельной задачей, и в ней количество независимых переменных может быть сведено к одной-единственной безразмерной переменной [6]. Для сферически симметричного случая безразмерная автомодельная переменная обычно выражается так:

где р* — характерная плотность, Е — постоянная, имею­ щая размерность энергии, и R — радиус ударной волны. Второе равенство в (7) позволяет определить константу пропорциональности а между величинами Е и До и по­ казывает также, что Xs = 1.

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

форме [8]:

 

p - t y i ? ) ф(р/р.)>

(8)


14

К. А. КОТ

где ср(р/р*) — произвольная функция приведенной'плот­ ности р/р*, а ф(Г)— функция температуры, которая в предположении гомотермичности зависит только от вре­ мени. Функция ф(7') должна иметь размерность давле­ ния. Из соображений размерности следует, что функция

ф(Г) пропорциональна t~°!\ Ее можно выразить через скорость ударной волны:

4>(D = /tp.£/2.

(9)

Здесь К — безразмерная постоянная, которая опреде­ ляется из решения задачи. Используя условие гомотер­ мичности, уравнение состояния и уравнение (9), можно исключить давление из дифференциальных уравнений и граничных условий задачи. Таким образом,

-§ f= 4 f Ж Л ф(р/р.)] = к и \ ' ,

причем штрих означает дифференцирование по аргумен­

ту р/р*.

Чтобы исключить из уравнений задачи входящую в них явно постоянную К, можно ввести новую автомо­ дельную переменную [8]:

 

Z =

г

 

( 10)

Поскольку

= 1, то получим

Zs = l / Y К- Используя

уравнение (10), выражаем частные производные через автомодельную переменную Z:

д

__ Zs d

д

UZ d

 

, . . .

dr

~ R dZ

dt ~

R

dZ

-

Кроме того, все зависимые переменные можно опреде­ лить как некоторые функции от автомодельной перемен­ ной Z:

и-

 

 

p =

p.G(Z),

ф(Г) =

p.t/2

p =

( 12)

p.C/’tf ( Z )== -*k££-<p(G).

 

z2 ’

 

 

Подставляя

(12)

в дифференциальные уравнения задачи

и используя

выражения

(11), получаем обыкновенные


МОЩНЫЕ ПОДВОДНЫЕ ВЗРЫВЫ

15

дифференциальные уравнения по автомодельной пере­ менной Z. После преобразования выведем систему урав­ нений, каждое из которых содержит только одну произ­ водную:

dF

F

^ Z ( Z - F ) - 2q/ (G)

dZ ~

Z

Ф' (G) - (Z - F)2

I dG __ F | - Z - 2 (Z - f)

G dZ

Z

Ф' (G) ~ (Z ~ F)2 1

п/ (Гг) -

dtp

dtp

dG

d(P/P.) •

 

(13a)

(136)

Аналогичное введение автомодельной переменной в гра­ ничные условия приводит к следующим выражениям:

G0 =

Gs (> - £ ) ■

 

G0 = Gs('i

' ■ ) * + * ' ? > .

o

1

ZsJ

Zs

О

 

 

'

 

 

 

(14a)

(146)

(14b)

Таким образом, решение зависит от вида функции cp(G), и, вообще говоря, полученная система уравнений не мо­ жет быть исследована непосредственно из-за ее неавто­ номной формы.

Для уравнения состояния достаточно общего вида

р = ^(Т){арч ± Ь),

которое соответствует термодинамическим данным для воды при высокой температуре, можно получить авто­ номную систему уравнений. Здесь ф(Г) по-прежнему является произвольной функцией от температуры. Част­ ную функцию плотности, которая лучше описывает тер­ модинамические свойства воды, возьмем в виде <p(G) =

= Gv — 1 =

(р/р*)у— 1, и тогда ее производная

равна

qp'(G) -

где Р =

Y — 1-

 

Проведем теперь преобразование переменных, а

именно

 

 

 

 

V __ ^

Т17 __ Cp/ (G)

/1


16

 

К. А. КОТ

 

Подстановка

этих

переменных и их производных

в. урав­

нения (13)

после

упрощений даст следующие

соотно­

шения:

 

 

 

z dr

+

(16а)

У

dZ

W - (1 - У)2

 

Z 6 V

р ф - у Н и г - и - т

(166)

У

dZ

W ~ ( \ - У)2

 

Эта система уравнений является автономной и путем простого деления первого уравнения на второе может быть сведена к одному обыкновенному дифференциаль­ ному уравнению

dY

У

2У2 — 7У + 5 — 6W

(17)

dW

W

4У2 ( 1 + Р ) — У (8 + р) + 4 (1 — W)

 

Граничные условия на ударной волне нельзя непо­ средственно выразить через преобразованные перемен­ ные, потому что они содержат параметр Zs, который заранее неизвестен и должен быть определен из реше­ ния задачи. Однако для любого выбранного значения Zs можно найти величины У« и IFS по формулам преобра­ зования. Следовательно, геометрическое место всех воз­ можных состояний на скачке в плоскости У, W легко определяется. Заметим, что при конечных, но ненулевых значениях плотности величина W становится бесконечно большой, когда Z равно нулю. Величина У при Z = О равна dF/dZ, и дальнейшее исследование показывает, что У = 0.

Уравнение (17)

имеет особенности там, где знамена­

тель обращается

в нуль, т. е. когда W = 0, или

4У2(1 + р) — У(8 +

р) + 4(1 — W) = 0 . Из определения

№ видно, что первое из вышеуказанных условий требует, чтобы либо Z = оо (бесконечно большой радиус), либо величина плотности обращалась в нуль или в бесконеч­ ность в зависимости от значения (3. Все эти условия фи­ зически неприемлемы. Следовательно, надо рассматри­ вать второе из вышеуказанных условий. Чтобы произ­ водная dY/dW оставалась конечной, когда знаменатель дифференциального уравнения обращается в нуль, не-


МОЩНЫЕ ПОДВОДНЫЕ ВЗРЫВЫ

17

обходимо, чтобы одновременно и числитель проходил через нуль. Это требование приводит к следующим усло­ виям совместности: Y = 0, или 2 Уг — 7К+ 5 — 6W — 0. При этих условиях существуют три особые точки. При

частном выборе величины

р = —0,05 одна

из этих осо­

бых точек получается при

отрицательном

значении W,

что соответствует либо отрицательной плотности, либо отрицательной величине Z. Значит, в области, имеющей физический смысл, будут расположены две особые точ­

ки, а именно Yx = 0, W\ =

1 и Y2

=

'A.

= 9Аб-

Дифференциальное уравнение

(17) имеет вид

dY

Р (W, Y)

 

 

dW ~

Q (W, Y)

 

Из теории таких уравнений [11] известно, что поведение данного уравнения в окрестности изолированной особой точки идентично поведению уравнения, имеющего лине­ аризованные формы числителя Р и знаменателя Q, по­ лученные при их разложении в ряд Тейлора, т. е.

dY _

a( W - W ) + b(Y - Y )

d\V

c ( W - W ) + d ( Y - 7 )

Здесь через W, У обозначены координаты особой точки. Поведение интегральных кривых такого дифференциаль­ ного уравнения в окрестности особенности определяется корнями б характеристического уравнения, которое за­ писывается так:

с — б

d

= б2 — б + с) {ad be) = 0.

а

Ь — б

 

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

Исследование линеаризованной формы уравнения (17) в фазовой плоскости показы! лрт, чтп первая особая

Гсс. П'<|5/:.*ЧмС,я

науч.чо-то.хнл .аская

бполнотеьа СССР

скэс:;-пдг.р

читального зада