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

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

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

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

Добавлен: 15.10.2024

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

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

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

УДАРНЫЕ ВОЛНЫ В ВОДЕ С ПУЗЫРЬКАМИ ВОЗДУХА

jg]

ный случай, когда сферические пузырьки воздуха фиксированного радиуса R и постоянной начальной температуры Г2 (после сжатия) окружены водой с температурой Т\. Если пренебречь изменениями коэффициента термодиффузии Da в воздухе, средняя температура в пузырь­ ке в момент времени t будет выражаться в виде [10]

Г = Г, + 6(7\ 7 Tl) £ ~~2~ е~°а<

(14)

Относительное изменение средней температуры воздуха по сравнению с изменением от начального значения Гг до предельного Т\ равно

 

6

оо

1

 

Т2 — Т

V

(15)

т»-т,

я2

Z i

п2 е

 

 

П=1

 

 

В дальнейшем будем предполагать, что охлаждением пу­ зырька можно пренебречь, если относительное охлаж­ дение составляет менее 10%, а если охлаждение составляет более 90%, то можно пренебречь разницей температур пузырьков воздуха и окружающей воды. Со­ ответствующие времена можно найти путем численного решения уравнения (15)

 

 

Uо%~ 0,0009R2/Da,

 

(16)

 

 

t90%~0,18R2/Da

 

(17)

при

температуре 21 °С,

давлении 14,7 фунт/дюйм2,

при

£>а =

2,2-10~4

фут2/с. При сжатии величина Па=Ха/раСа

изменяется, так как меняется плотность ра, а величина

ха/са зависитот температуры приближенно как

Г0,8.

Следовательно, для пузырьков радиусом Ri при нор­

мальной плотности вышеприведенные равенства запи­

шутся в виде

trn

~4(tf?/tf)(7y7f'8,

 

(18)

 

 

 

 

 

£эо%

~ 800(/?i/.ft) ( 7 У Т ) 0’8,

 

(19)

где t в

секундах,

R в футах, Ti =

21 °С.Если

пузы

в начальный

момент

характеризуются

температурой

6 Зак, 741


162 6. Р. ПАРКИН, ф. р. ГИЛМОР, Г. Л. БРОУД

Г, = 210С и давлением Р\ = 14,7 фунт/дюйм2, а затем сжимаются адиабатически до давления рг (при у = 4/3),

то RIR\ = {p2lP\)~'U и Г/Г, = (р2/р,)'/', так что из уравне­ ния (18) получаем

(2°)

Если охлаждение происходит при постоянном давле­ нии р2, то температура пузырька Г стремится к темпе­ ратуре Г,, а радиус R — к его изотермическому значе­

нию Ri (p2/pi)_ /s. Поскольку при 90%-ном охлаждении основное время протекает тогда, когда пузырьки почти холодные, то уравнение (19) дает

 

(21)

Времена охлаждения на 10 и 90%

приведены на

рис. 1 наряду с временами схлопывания,

вычисленными

ранее. Видно, что пузырьки, имеющие в начальный мо­ мент радиус менее 0,01 дюйма, охлаждаются почти до температуры воды через несколько миллисекунд или быстрее. Несомненно, что времена схлопывания для та­ ких пузырьков имеют тот же порядок по величине, так что в точной теории при исследовании поведения пу­ зырьков следует рассматривать как динамику, так и тепловые потоки, но это не изменяет вывода о том, что времена установления динамического и теплового равно­ весия имеют порядок нескольких миллисекунд. Для бо­ лее крупных пузырьков времена охлаждения больше, а для пузырьков, радиус которых в начальный момент имеет порядок 1 дюйма или больше, во многих случаях достаточным является адиабатическое приближение (при условии что при схлопывании пузырьки не распадаются на более мелкие).

Измерения Кэмпбелла и Питчера [3] по ударным вол­ нам малой амплитуды в воде с пузырьками воздуха ра­ диусом порядка 0,005 дюйма показали, что давления достигают своих изотермических значений при временах, меньших времени разрешающей способности аппарату­ ры, которое составляет около 2 мс. Это не противоречит величине ^о% = 0,2 мс, данной на рис. 1,

УДАРНЫЕ ВОЛНЫ В ВОДЕ С ПУЗЫРЬКАМИ ВОЗДУХА

163

Д. Растворение пузырьков воздуха в воде

Даже когда воздушно-водяная смесь первоначально насыщена воздухом, при увеличении давления раствори­ мость воздуха в воде увеличивается. Задача о растворе­ нии сферических пузырьков в полностью ненасыщенной воде при температуре 21 °С и постоянном давлении была решена Эпштейном и Плессетом [11]; их результаты можно записать в виде

^ = 1,5 • 108

R2,

(22)

где td — время в секундах, требуемое для

того, чтобы

пузырек полностью растворился,

a R — радиус в футах

до начала растворения. Ниже будет показано, что время td намного больше, чем времена схлопывания и охлажде­ ния, что оправдывает предположение о постоянстве тем­ пературы и давления, принимаемое в данной работе (а также позволяет не учитывать растворимость возду­ ха при рассмотрении схлопывания и охлаждения пу­ зырьков). Для пузырьков, имеющих в начальный мо­ мент радиус R 1 и давление рь при сжатии до давления р2 уравнение (22) дает

td =\,5-lO*R](p2lPl) - \

(23)

Если вода в начальный момент насыщена воздухом при давлении р\, то уравнение (23) справедливо только при условии pdp\ > 1. При P2IP1 2 вода после сжатия ста­ новится полунасыщенной, а время td увеличивается в два раза [11] по сравнению с величиной, определяемой по уравнению (23).

Из формул (21) и (23) получаем соотношение

 

- ^ - = 2-10ЧРг/р,)"'.

(24)

из которого видно, что времена растворения много боль­ ше времен охлаждения для всех случаев, представляю­ щих интерес. Для пузырьков радиусом больше 10-3 дюйма при давлениях меньше 104 фунт/дюйм2 время раство­ рения, вычисляемое по выражению (23), обычно состав­ ляет больше 10-2 с. Для пузырьков бодьщих размеров и

6*


164

Б. Р. ПАРКИН, Ф. Р. ГИЛМОР, Г. Л. БРОУД

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

Е. Уравнения состояния для воздушно-водяных смесей

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

1)тепловое равновесие без растворения воздуха;

2)изолированные пузырьки (растворения не происходит);

3)тепловое равновесие и равновесие процесса рас­ творения.

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

УДАРНЫЕ ВОЛНЫ В ВОДЕ С ПУЗЫРЬКАМИ ВОЗДУХА

165

скольку в обычных воздушно-водяных смесях отношение масс воздуха и воды мало, то после растворения воз­ духа смесь будет вести себя приближенно, как чистая вода.

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

 

Р

__г_ В

 

(25)

 

Р аТа

Ш ’

 

 

 

 

где р — давление

(одинаковое для

воздуха

и воды);

ра, Та, пг — плотность, температура и молекулярный вес

воздуха соответственно; В — универсальная

газовая по­

стоянная. Плотность воды рш связана

с давлением соот­

ношением

РШ= Р*(1 +Plk),

 

(26)

 

 

где k — модуль

объемного

сжатия,

а р* — плотность

воды при нулевом давлении

(которая незначительно от­

личается от плотности при давлении в 1 атм). При тем­ пературе 21 °С, которая принималась в численных рас­ четах в данной работе, имеем

р* = 62,4 фунт/фут3,

и

k = 3 ■105 фунт/дюйм2

при низких давлениях. Хотя k увеличивается при увели­ чении давления, ошибка в вычислении рю вследствие предположения о постоянстве k в соотношении (26) со­ ставляет только 10% при давлении 104 фунт/дюйм2, в то время как ошибка, вносимая в соотношение на удар­ ной волне, обычно намного меньше, поскольку основное изменение плотности происходит за счет сжатия пу­ зырьков.

Пусть теперь р, есть отношение массы воздуха к мас­ се воды в смеси. Приравнивая полный объем массы смеси l + i-i к объему воздуха плюс объем воды, полу­ чаем

1 + 1* _ _н_ , __L _ J L _ i________ 1

(27)

р

Ра

Pw

Ра "Г Р* (1 + p !h ) '


166 Б. Р. ПАРКИН, Ф. Р. ГИЛМОР. Г. Л. БРОУД

где р — плотность смеси. Исключая ра из уравнений (25) и (27), получаем уравнение состояния для смеси

Р Г1+ и _____ 1

Т В

(28)

Та L рр

|хр* (1 +

p/k) J т '

Для теории теплового равновесия (при отсутствии про­ цесса растворения) и теории изолированных пузырьков уравнение (28) описывает термодинамическое состоя­ ние смеси в каждый момент времени. Если смесь, под­ чиняющаяся теории растворимых пузырьков, внезапно сжимается, то уравнение (28) справедливо только до сжатия. После сжатия термодинамическое состояние описывается уравнением (26).

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

АЕ =

С + ц с0

А Г »

с

 

(29)

1+ ц

1+ Ц АД

где с — удельная теплоемкость воды

(с= 1), cv —удель­

ная теплоемкость

воздуха

при

постоянном объеме, а

АТ — изменение температуры смеси

(одинаковое

для

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

Для теории изолированных пузырьков можно запи­

сать

ЦС0

 

 

АЕ =

А Та + 1+1» *

(30)

1+

где AEw— изменение

энергии воды, рассчитанное

на

единицу массы; AEw можно заменить величиной сД7\», но в любом случае мы имеем дополнительную тер­ модинамическую переменную (Ew или Tw) для смеси по сравнению с обычной жидкостью. Для решения задач о течении жидкости должны быть получены дополни­ тельные соотношения, содержащие термодинамические параметры. Для невязких адиабатических течений изме­ рения энергии водь! обусловлены работой обратимого


УДАРНЫЕ ВОЛНЫ В ВОДЕ С ПУЗЫРЬКАМИ ВОЗДУХА 167

сжатия

АЕШ= - JРа pd( 1/рш) = Pi

_ А .

П п

1

+

р *

L

1

+

p*!k

_ L ____!___________ ! _

Pi/k

■*" 1 + p*/k

i + р ,/ *

р \~ р\

2 р| —Р]

3 Р2 —р] 4 pf —pf

(3 1 )

3 pV h

б- р*£4 "*

Для течений с диссипацией, таких, как течение с ударными волнами и течение в пограничном слое, необходи­ мо найти способ разумного разделения диссипируемой энергии между фазами воздуха и воды. Поскольку вплоть до давлений 104 фунт/дюйм2 ударные волны яв­ ляются довольно «слабыми» и недиссипирующими в воде, в то время как диссипация, необходимая для со­ хранения ударного профиля, по-видимому, имеет место в волнах сжатия, которые отражаются назад и вперед внутри индивидуальных пузырьков, то для таких удар­

ных волн

разумно использовать уравнение (31).

(Течение

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

в данном

отчете.)

Такая аппроксимация справедлива

только в

том случае, когда содержание воздуха в жидкости не очень мало. Если допустить, что ц-э-0 при Та^-оо, то в любом бесконечно малом объеме воздуха количество диссипируемой энергии будет конечным. Так как извест­ но, что диссипируемая энергия является величиной третьего порядка малости по отношению к разности дав­ ления в ударной волне, то не удивительно, что величи­ на ДТа может быть вычислена с большой погрешностью, если только в разложении (31) не учитывать несколько членов ряда; этот факт был случайно обнаружен при проведении излагаемых здесь численных расчетов.

Вывод уравнений (28) — (30) основан на законе идеального газа и предположении о постоянстве удель­ ной теплоемкости воздуха в пузырьках. Эти предполо­ жения, строго говоря, несправедливы при высоких дав­ лениях, рассматриваемых здесь. Однако основным фи­ зическим свойством смеси, которое определяет свойства потока, является изменение объема при сжатии. При