ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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/:.*ЧмС,я
науч.чо-то.хнл .аская
бполнотеьа СССР
скэс:;-пдг.р
читального зада