Жекамухов М.К., Жекамухова И.М. Сход снежных лавин и возникновение воздушных ударных волн

Научная статья

Электронный журнал ИССЛЕДОВАНО в РОССИИ

512

Сход снежных лавин и возникновение воздушных ударных волн.

Жекамухов М. К., Жекамухова И. М. ( Jekamuhova@mail.ru ) Кабардино-Балкарский государственный университет

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

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

В данной работе рассматриваются отдельные аспекты этой проблемы.

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

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

d(mV)аа /.а ,аа чаа cfpSV2

Ч----- = mg (sin \|/-kf cosy)Ч'-Ч--------- ,аа (1)

Электронный журнал ИССЛЕДОВАНО в РОССИИ 513а dVаа /.аа ,а ч pV2 Sа 1 dm

Ч- = g(siny-fcosy)-cf Ч------------- V---- Ч. (2)

atа 2 m m at

Пусть X - длина лавинного снега, h - его средняя высота, Н - ширина захвата

лавины. Тогда отношение

S _аа 1 m psX'

где ps - плотность лавинного снега.

Закон изменения массы лавины можно записать в виде

XPsohOHV,

где pso - плотность снежного покрова на склоне горы, х - коэффициент захвата снега лавиной, ho - мощность снежного покрова. При этом отношение

1 dm Ps0 h0v "X-

m dt Psh?i

Уравнение (2) теперь принимает вид

dVаа (жа fа ч aV2

Ч- = g(siny-fcosy)--------- Ч, (3)

dtа X

pа Ps0h0

где a = Cf----- 1--- Ч.

Ps Psh Ваа дальнейшемаа предположим,аа чтоаа отношениеаа % остаетсяаа постояннымаа при

движении лавины.

Приравниваяаа правуюаа частьаа уравненияаа (3)аа каа нулю,аа найдемаа установившуюся скорость движения лавины:

Vm = J Ч (sin \|/ - f cos \\f) . V a

Введем также некоторое характерное время т, определяемое равенством Vmx = X,

т.е.

g sin \|/ Ч f cos\|/

Электронный журнал ИССЛЕДОВАНО в РОССИИ 514 f = l-V2.аа (4)

где V =------- безразмерная скорость, t =------ безразмерное время.

Решение уравнения (4), удовлетворяющее начальному условию V(0) = 0, имеет вид

V = Vmth-. (5)

Отсюда видно, что скорость V стремится асимптотически к установившемуся значению Vm.

Для закона движения снежной лавины получаем формулу

x = VmxlnchЧ.аа (6)

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

Как известно [1], впереди поршня, толкающего воздух в трубе, профиль волны в каждый момент времени описывается формулой

x = (c0+^ujt+f(u),аа (7)

где и - скорость течения, Со - скорость звука в воздухе, f(u) - произвольная функция, зависящая от скорости, у = 1,4 - показатель адиабаты воздуха.

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

ж- 2С Хаа х = 4^- (8)

а (у +1)' а (у +1)

Электронный журнал ИССЛЕДОВАНО в РОССИИ 515а

VmTlnch^ = fc0+^Vmth^jt + ffvmth^

(9)

tаа z

Обозначим Vm т th Ч = z, тогда t = т arc th -

Подставляя значение t в равенство (9) и учитывая соотношение

, tа 1 1

СПЧ= ,

'1-th

2аа z

11-th

InchЧ= Ч In

т 2

2аа ^

vа V,

будем иметь

Уш1ы

,2аа

Vаа V,

cq н------ z xarcth------- i-f(z).

Отсюда получим

f(u)=-bllln

.2аа

Vаа ^mj

y+1 ]аа , u

eg н----- u т arc th

Равенство (7) теперь принимает вид

х=аа с0 +

у + 1

uаа t-

Vmx

.2^

y+1 ]аа , u

eg н----- u xarcth

Время образования разрыва (ударной волны) может быть найдено из условия [1]

ГдхЛ

и=0

\диД

Оно приводит к равенству

ylt_cox=0|

Таким образом, получаем формулы, определяющие время и место образования ударной волны

2с

2с,

(10)

(у + 1)'

(У + 1)'

а\

x = c0t

Электронный журнал ИССЛЕДОВАНО в РОССИИ

516

где а.\

'ni/

На рис.1 схематически показано распределение скорости воздуха между поршнем и фронтом ударной волны в произвольный момент времени t.

Рис. 1 Распределение скорости газа перед поршнем в момент времени t.

Приведенные выше количественные формулы относятся к движению поршня внутри бесконечно длинной трубы. При движении лавин ситуация совершенно другая -лавина обтекается воздушным потоком. Следовательно, условия возникновения лавин будут более жесткими. Но если даже допустить, что лавина толкает весь воздух, находящийся впереди нее, как это имеет место в случае движения поршня в трубе, то подставляя в формулу (10) значение т = 5,5 с, получающееся при Vm =18 м/с и Х= 100м, получаем t = 84 с, х = 280 км; при максимально возможной скорости движения лавин, равной 50 м/с, t = 29 с, х =9 км.

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

Расчет удара снежной лавины о неподвижное препятствие. При проектировании мостов и линий передач, а также при строительстве противолавинных сооружений в лавиноопасных горных районах необходимо иметь расчетные данные о возможных значениях силы удара лавин о препятствие. Существует ряд эмпирических и полуэмпирических формул, с помощью которых такие расчеты проводятся. Наиболее распространенные из них приводятся, например, в [2]. Однако эти формулы физически обоснованы недостаточно, а расчетные значения силы удара лавин о препятствие по ним частоа плохоа согласуютсяа с экспериментальными данными.а Отсюда следует

Электронный журнал ИССЛЕДОВАНО в РОССИИ 517 Р = Ар\аа (11)

где А - своя постоянная для каждого типа снега, а показатель к меняется в пределах от 7 до 14, в зависимости от типа снега.

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

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

Будем задавать это уравнение в виде

(12)

Р1=Р0

'--Ps<T

VPi-Psl

где р; - плотность льда, psi - плотность снега и pi - давление в ударной волне, параметр к играет роль показателя адиабаты снега. Хорошее совпадение с экспериментальными данными [2] получается, если показатель к принять равным 3,5.

Таким образом, ударная адиабата снежно-воздушной массы представлена в виде достаточно простой аналитической формулы; она справедлива до давлений, сравнимых с пределом прочности кристаллов льда, т.е. до 5 - 10 атм. При более высоких давлениях показатель к будет возрастать.

Электронный журнал ИССЛЕДОВАНО в РОССИИ 518

Рис.2. Кривые зависимости давления в ударной волне от плотности снега.

ps, г/см

Пусть снежная лавина длиной X ударяется о неподвижное препятствие, расположенное перпендикулярно к направлению ее движения. После удара по лавинному снегу бежит отраженная ударная волна; между препятствием и ударной волной возникает застойная зона.

Законы сохранения массы и импульса на фронте ударной волны можно записать в виде

PsODs =Psl(Ds -ul).а (13)

Po+PsoDs^Pi+PsilDs-Ui)2,аа (14)

где Ds - скорость ударной волны в снеге, ui - массовая скорость за фронтом ударной волны.

Из теоремы сложения скоростей при сложном движении тел следует, что массовая скорость равна скорости удара: ui = V.

Присоединяя к системе (13) - (14) уравнение ударной адиабаты (12), получим замкнутую систему уравнений относительно трех неизвестных величин : pi, р si и Ds.

Из равенств (13) и (14) легко получается известная формула С. А. Христиановича, выражающая зависимость между давлением и плотностью за фронтом ударной волны

PsOPslаа Лг2

-V"

PsQ-Psl

Электронный журнал ИССЛЕДОВАНО в РОССИИ 519

Эта формула проста и удобна при оценочных расчетах. Из равенства (13) следует

.Ps(L

Psl

а из равенств (13)и(14) получим

Pl=Po+PsoVDs Подставляя (15) в (16), будем иметь

(16)

PL = 1+P_sOv

Ро Ро

Из равенства (12) выразим р si через рь

PsO Psl

(17)

Psl =Pi "(Pi -Pso)

^Po Подставляя (18) в (17), получим уравнение, связывающее V и pi

(18)

IPsOаа )

fe-i:

)(^1/k-l)

^l/k.

_1+PsO Pi

(19)

.vs

где t,

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

Po Vs^ 10,5 м/с.

Подставляя значения psihVb(15), будем иметь

f пV

vPsOy

PsO

Pi-PsO

^l/k_1 + Ps^

fe-D-wrf

(20)

Скоростьаа распространенияаа звук ваа снежнойаа массе,аа сжатойа ударнойаа волной, определяется формулой

1 k+l

(21)

cs0^

csl

^dpjV VdPsb

где cs0

kpp .Pi" PsO

Электронный журнал ИССЛЕДОВАНО в РОССИИаа 520а Таким образом, все параметры снежно-воздушной массы за фронтом ударной волны выражаются через давление.

На рис.3а приводятся кривые зависимости относительного давленияа у , где

Ар = р^-ро, от скорости удара лавины, построенные по формуле (19) при различных значениях плотности снега до удара.

АР

Ро

25 20 15 10 5

Рис. 3 Зависимость давления в ударной волне от скорости удара лавины о препятствие.

Кривые 1, 2, 3, 4 - для ps = 0.1,

0.2, 0.3,аа 0.4аа кг/м3

соответственно.

50 V, м/с

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

известной формуле р = pso Vа , которой часто пользуются при расчетах противолавинных

сооружений, примерно в 2 - 4 раза в интервале плотностей снега от pso = 200 кг/маа до ps0 = 400 кг/м3.

Н рис.4аа приводятсяаа кривыеаа зависимостиаа Dsаа отаа скоростиаа удар лавины,

построенныеаа поаа формуле

различных

D, м/с

жсо

(20)а при значениях pso.

.и.

ж?м

41J

Рис. 4 Зависимость скорости ударной волны в снеге от скорости удара лавины о препятствие.

V, м/с

Электронный журнал ИССЛЕДОВАНО в РОССИИ

521

Кривые 1, 2, 3, 4 - для ps = 0.1, 0.2, 0.3, 0.4 кг/м соответственно.

Эффект действия лавины на неподвижное препятствие определяется величиной ударного импульса К = p-At, где At - длительность удара лавины о препятствие. Величина

At состоит из двух слагаемых: At = Ati + А12, где Atj = Ay. - промежуток времени, в

течение которого ударная волна, отраженная от препятствия, достигает хвостовой части

лавины, At 2 = 7аа - промежуток времени, в течение которого волна разгрузки бежит по /csl

сжатому снегу. После достижения волной разгрузки поверхности препятствия, давление в

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

снега, можем записать

(

1а 1

Ч + Ч

^Dsа csl.

Api л-----------

(22)

Зная скорость удара V, по графикам 4 могут быть определены соответствующие значения Ds, а затем по формуле (22) и ударный импульс.

Воздушные ударные волны, возникающие при ударе снежной лавины о препятствие. Обширный материал наблюдений за сходом снежных лавин показывает, что сход мощных лавин часто сопровождается возникновением воздушных ударных волн. При этом имеются случаи, когда лесоповалы и другие разрушения происходят не впереди лавины, а в стороне от нее [3]. Эти разрушения уже не связаны с воздушными потоками, образующимися вокруг лавины, а являются следствием действия ударных волн. Однако до настоящего времени механизм образования этих волн остается неясным. Ниже предпринята попытка установления механизма возникновения воздушных ударных волн при ударе лавин о препятствие.

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

На рис.5 схематически показана разгрузка снега в ударной волне.

"Ж"

Электронный журнал ИССЛЕДОВАНО в РОССИИ Pi

Csl

522

Р2

u3 3

Рз

D 4аа ро

Рис. 5 Схема разгрузки сжатой снежно-воздушной массы.

1 - масса снега, сжатого ударной волной; 2 - волна разрежения; 3 - воздушная ударная волна; 4 - невозмущенный газ (воздух); ab - фронт разгрузки снега; cd - передний фронт разгрузки, который играет роль поршня; ef- фронт воздушной ударной волны; стрелками на рисунке показаны направление движения среды и перемещения слабых разрывов, ограничивающих волну разрежения.

Волна разгрузки бежит по сжатому снегу со скоростью csi, а передний фронт волны разрежения движется вправо со скоростью V2, толкая как поршень воздух; зона 3 между cd и ef охвачена воздушной ударной волной, передний фронт которой движется со скоростью D по невозмущенному воздуху.

Скорость течения в зоне разгрузки определяется формулой [4]

dp Pscs

(23)

u2= j

Р2

где p2 - давление на стыке между снежной массой и сжатым ударной волной воздухом,

cs0

k+1 2k

( ~ Y

местная скорость звука в снеге, ps = pj -(pj -pso)

vPoy Подставляя значения cs и ps в (23), получаем

k-1 zk+1dz

k-1 ,k+l

2k po Pi-Po

k +1 pj у kp

PsO

f ~ Y

vPoy

(24)

'1 J

k+1 где z = E, 2k .

При k = 3,5 интеграл в равенстве (24) точно не берется. Однако, если произведем k-1

замену отношения

0,556а на 0,5, то при этом интеграл раскрывается, а ошибка,

k + 1

возникающая при такой замене, не превосходит 10%. Таким образом, можем записать

Электронный журнал ИССЛЕДОВАНО в РОССИИ 523а Давление и скорость за фронтом ударной волны выражаются через число Маха

М = Ч и параметры невозмущенного газа известными формулами [1] со

Рз=Ро

2у у + 1

(

2а У-1

2у

2с0

у + 1

(

^а Mz)

где Со - скорость звука в воздухе.

Подставляяаа значенияаа рзаа иаа изаа ваа (28),аа получимаа трансцендентноеаа уравнение относительно М:

2 ПГаа !Д../-

(

2 с0

у + 1

2а У^1

(26)

2lL

cs0

Po "VT + 1

Обычноа ударныеаа волны,аа возникающиеаа приа сходеа лавин,аа являютсяаа слабыми, поэтому число М можно представить в виде

М = 1 + а,аа (27)

где а л 1.

Подставляя (27) в (26) и ограничиваясь линейными по а членами, получим

^P1V

у + 1аа Vs

(28)

vPoy

8 Уаа c0cs0

Из формулы (31) видно, что а действительно является малой величиной. Значения D, Уз и рз выражаются через а следующими формулами

D = c0(l + a);а u3

4с0а 2у

ЧЧа;а р3=роЧ-

у+1а у+1

2а +

у + 1 ^7

Чтобы получить представление о характерных значениях параметров воздушных ударных волн, возникающих при ударе лавин о препятствие, примем pso = 300 кг/м ,

V = 36 м/с; при этом отношение ^у = 10, a = 0,016, D = 1,016со, из = 9 м/с, а амплитуда

ударной волны Ар = 0,037, р0 = 370 кг/м2.

Для сравнения укажем, что динамическое давление, которое создается при натекании воздушного потока на неподвижное препятствие, составляет при скорости потока V = 20 и V = 50 м/с соответственно 40 и 250 кг/ма .

Электронный журнал ИССЛЕДОВАНО в РОССИИ 524 Выше мы рассматривали прямой удар лавины о плоскую поверхность препятствия. В случае косого удара внутри снежно-воздушной массы будут возникать отраженные ударные волны, наклоненные к поверхности препятствия под некоторым углом, отличным от прямого. При выходе таких волн на боковую поверхность снежной лавины возникают боковые воздушные ударные волны, направления движения которых не совпадают с направлением движения лавины. Такие волны могут возникать, например, при косом ударе лавин о скальные препятствия, при выходе лавин со склона горы на горизонтальную поверхность, т.е. при любых резких поворотах лавин, движущихся с большой скоростью.

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

Обозначения.

kf - коэффициент трения скольжения; \|/ - угол наклона склона горы к горизонту; Cf -коэффициент сопротивления; р - плотность воздуха; S - площадь миделя лавины; V -скорость лавины; m - масса лавины; g -ускорение свободного падения; t - время; h -средняя высота лавины; Н - ширина захвата лавины; X - длина лавины; ps - плотность снега; pso - начальная плотность снега (до удара); psi - плотность снега после удара; х -коэффициент захвата; ho - мощность снежного покрова; а - безразмерный коэффициент; Vm - максимальная скорость лавины; т - характерное время; V - безразмерная скорость; t - безразмерное время; а\ - характерное значение ускорения лавины; х - координата; у -показатель адиабаты воздуха; Со - скорость звука в воздухе; а - ускорение поршня; и -скорость течения воздуха; ро - плотность воздуха в порах снега до удара; pi - плотность воздуха в порах снега после удара; р; - плотность льда; ро - атмосферное давление; pi -давление в сжатом снеге; к - показатель адиабаты снега; Ds - скорость ударной волны в снеге; ui - массовая скорость за фронтом ударной волны в снеге; Vs - постоянная

величина, имеющая размерность скорости; Е, = "у ; К - ударный импульс; cs - местная

скорость звука в снеге; иг - скорость передней границы зоны разрежения снега; из -скорость газа за фронтом воздушной ударной волны; рг - давление на стыке снега и

Электронный журнал ИССЛЕДОВАНО в РОССИИ 525а итература.

андау Л.Д., ЛифшицЕ.М. Гидродинамика. М.: Наука, 1986.
Доловаа М.А.,аа Халкечеваа В.А. Физик снег иаа динамик снежныхаа лавин. Л.: Гидрометеоиздат, 1972.
ФляйгВ. Внимание, лавины!. М.ИЛ, 1960.
Зельдовичаа Я.Б.,аа Райзераа Ю.П. Физик ударныхаа волнаа иаа высокотемпературных гидродинамических течений. М.: Наука. 1966.

Все научные статьи

Blog

Жекамухов М.К., Жекамухова И.М. Сход снежных лавин и возникновение воздушных ударных волн