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

книги из ГПНТБ / Чугаевский, Ю. В. Элементы теории нелинейных и быстропеременных волновых процессов

.pdf
Скачиваний:
12
Добавлен:
19.10.2023
Размер:
7.46 Mб
Скачать

)ид решений этого уравнения зависит от знака коэффициента

при

зтором

члене. В случае положительного

множителя

 

с а и и>0 )

решением будет

"медленный" солифон (приложение,

а)

)

 

 

 

 

x- c t\

 

 

 

 

 

Л

 

 

 

 

 

 

(69)

 

 

‘ М

 

=

 

 

 

 

 

3

 

 

 

 

 

 

Vfc =

3 0г - с 2)

 

 

 

 

 

 

 

 

 

 

 

 

(70)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ширина которого, как видно, изменяется от 0 до

 

 

при

изме­

нении

скорости

с

от 0

до

co = u)0J'/y.

 

 

 

 

 

 

 

 

"Быстрый"

солифон

( р с г

>

)

запишется

как

 

 

 

V = vc seek

 

 

*0 =

3

(С *-С о)

 

h =

 

2 .С

 

 

(71)

 

 

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\ ^ = COf

 

 

Лпределе

ejg —

о

"быстрый" солифон,

естественно,

совпа-

дает с

"коротким".

 

 

 

 

 

На ри с.2

показаны

 

про­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

фили

"медленных"

солифонов

 

 

 

 

 

 

 

 

для некоторых значений пара­

 

 

 

 

 

 

 

 

метров

v

h . Профили

 

"бы­

 

 

 

 

 

 

 

 

стрых"

солифонов

получаются

 

 

 

 

 

 

 

 

смещением приведенных кривых

 

 

 

 

 

 

 

по

вертикали

вверх на

 

вели­

 

 

 

 

 

 

 

 

чину

J vQ }

т .е .

будут

рас­

 

 

 

 

 

 

 

положены полностью в облаоти

 

 

 

 

 

 

 

положительных окоростей. Дви­

 

 

 

 

 

 

 

 

жение среды пооле "быстрого"

 

 

 

 

 

 

 

 

содифона

и перед

ним практи­

 

 

 

 

 

 

 

 

чески отсутствует. Медленный

 

 

 

 

 

 

 

 

же солифон имеет отруктуру -

 

 

Р и с .

2

 

 

 

внутреннюю область ( 0

^ v

 

 

 

 

 

6

 

)

положительных

скоро-

 

 

 

 

 

 

 

 

 

стей

и

внешвюю

область

{O&vz-,

-3 v(.

)

 

отрицательных,

 

 

 

 

Перед возмущением и после него жидкость движется

с

постоянной

скоростью

 

О 1/

 

 

 

солифон возможен

 

 

только

V = ----- ,

т .е . медленный

 

 

в набегающем потоке.

20

д) Влияние вязкости и релаксационных потерь, "плавниковая" деформация

Дифференцирование уравнения состояния с диссипативным чле-

ном

 

д*Р , 'f l '

32р

 

г др

 

 

Ж

 

 

 

 

 

 

 

 

 

 

 

 

~ ^Г х

3

d*dt

°

дх

 

 

дх

 

 

 

 

 

 

 

 

 

 

и подстановка

отсюда производной

дс

дх

в

уравнение

непре

рывности

 

(49)

даст

 

 

 

 

 

 

 

 

 

'°дх

з

 

 

 

 

11

+f °

 

д \

=

 

о.

(72)

 

дх2 + у дх

3xed t

 

 

 

Подставляя сюда значение

др

 

из

уравнения

движения

 

 

 

 

 

 

 

ж =

 

dv

 

d v \

 

 

d zv

( 73)

 

 

 

Ро ( d

t * V

д х )

 

~дхт

 

 

 

 

дх

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ш? f

 

d t

У

 

хм

dv + Р '

d*v + &

d*v

;

3 ^ - О

(74)

 

 

 

 

 

 

Ж * 7 ' ' 1 к г

А ~ г к г

 

Л '# г

1741

что после дифференцирования по времени дает уравнение

 

си' 111

 

d 2v

 

v ( dv

dv

,,

dzv )

,

d^v

 

 

 

дх?

 

-Y J F

r \ lt

 

+ ^ w ) + g ^ d F

 

 

 

 

 

д \

=

о,

 

 

 

 

 

 

 

 

 

(76)

 

 

dxzd t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оторое отличается

от

уравнения

(53) последним

членом,

у ч и т -

1ающим вязкие и релакоапионше потери, В случае отационарных волн отсюда подучаем

 

 

 

(76)

или после понижения поря.дка

 

 

 

 

 

,,/

 

v +\>v +ifvz+ <f>v = о ,

=

9

• ( 77)

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

 

v = ? 0 ( t ) T ( t ) ,

 

 

(78)

считая T ( t )

медленно меняющейся

амплитудой,

a v0 ( t )

ре­

шением уравнения (7?) без диссипативного члена,

т ,е ,- нулевым

приближением,

 

 

 

 

 

+ 9'*° +<f%z = 0.

 

 

(79)

Подставляя

(78) в (77) и имея в

виду (79 ),

получим

 

 

Т - 1 -

 

 

(80)

Если, в частности, принять в качестве нулевого приближения "бы­ стрый" солифон (71 ), то функция (80) будет иметь вид

 

Т= / -

sh ( - I х

) .

(80а)

 

 

h<fvc

\

h

J

 

 

Поскольку sh

- функция нечетная, влияние вязкости

и ре­

лаксации приведет к несимметричному искажению геометрии

соли-

фона:

задний его фронт будет

прогибаться, а передний ,

наобо­

р о т ,-

выпучиваться.

Профиль

уединенной

волны., будет деформиро­

ваться в сторону образования характерной формы "акульего плав­ ника" (р и с .З .а ).

Следует заметить, что деформация "плавникового" типа мо­

жет быть вызвана также наличием в системе нечетной, например,ку­

бической нелинейности (см ., например, ри с.21, кривая J3 ).

22

е)

Диш ерсионное "заострение"

профиля

волны

 

Учет в разложении функции с( к)

членов

четвертого

поряд­

на иоЛ

 

 

 

 

 

 

 

с(А) =* с0 -Ъ кг + d k 4

,

d=

const

 

 

приводит

к пополнению уравнения (64) производной пятого поряд­

ка

 

 

 

 

, ~dsv■

 

 

 

dv

dv

£

d3v

 

(64а)

 

dt

дх

 

дх 3

дх 5

 

 

 

 

 

Квазистационарному решению получающегося отсвда

уравне-

ния

dvm+ b v " -cv + -^- = 0

.ложно придать вид

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

vg

следует брать

длинный

солифок (66)

(ри с.3 ,6 , кривая

оС ) .

 

 

 

 

Из поведения функции

v"'

(кривая J3

) видно,

что

высшие

дишерсионные

члены

заостряют

вершину

солифона (кривая/ ), де­

формируя его

геометрию в сторону так

называемого

т р о х о —

и д а л ь н о г о

профиля Герстнера

(кривая & ) .

Гравитаци­

онные волны такого

типа без

труда можно наблюдать на

глубокой

воде в естественных условиях (подробнее о трохоидальных волнах

см. в

[2 2 ] ) .

 

 

 

ж) 0 числах Рейнольдса и Маха

 

В связи с тем, что в

нелинейной акустике

диспергирующей

жидкости

оказываются возможными стационарные волны, вообще го­

воря, любой скорости распространения, возникает

необходимость

ввести

определенные коррективы в понимание критериев Re и М.

Как

видно из уравнения (7 7 ),

критерием" недиссипативного

приближения является неравенство

 

 

 

 

( v v | <<

\fv\ ,

 

 

которое

 

удобно переписать,

введя

численный критерий

<.3

 

 

 

Ц ч>о-Гс ‘

|1- « l

Re >> /

 

 

 

 

 

(81)

 

 

 

/Г= -

 

 

«К

 

 

 

 

 

 

 

 

 

he

 

 

 

 

 

 

 

 

 

 

 

где

Re =

-

число Рейнольдса

для данного процесса, х =

 

 

и,

кроме

того, принято простоты ради

j v

 

- -

J3'

Для

 

 

 

3

весьма

медленных волн, следовательно,

 

даже при

 

Re<<1

чис­

ло

К

может значительно

превышать I

и допускать,

таким

обра­

зом, пренебрежение диссипацией. Так,

при

х

~

Ю- ^

и

R e ~

~

ХСГ^

-t

10“ ^ число

К

имеет порядок

Хс£

*

ХО^,

допускающий

бездиссипативное

приближение.

 

 

 

 

М =

-jr-

под с

 

Аналогично,

при

определении

числа

Маха

следует понимать не скорость линейного бездисперсионного

зву­

ка,

но

скорость

рассматриваемого

нелинейного

стационарного

волнового процесса.

 

 

 

 

 

 

 

 

 

 

 

 

§

4. Кноидальные волны, "бисолифоны",

 

 

 

 

 

 

 

 

 

"волновая

сверхпроводимость"

 

 

 

 

 

 

 

 

Уравнения ангармонического вибратора с квадратичной нели­

нейностью,

рассмотренные

в предыдущем параграфе,

 

имеют

перио­

дические решения, дающие стационарные периодические волны,опи­

сываемые эллиптическими функциями. Впервые такие решения были

получены, по-видимому, Кортевегом и де Вризом для гравитацион­

ных волн на мелкой воде и были названы ими

к н о

и д е л ь ­

н ы м и ( сп -

эллиптический

косинус).

 

 

Рассмотрим

уравнение (68)

при ^ с 2 < о-'2

, т .е .

медленный

волновой процесс. Если амплитуда волны меньше

амплитуды соли-

фона

 

 

 

 

3 ( с % - с г )

(82)

^=

то решение уравнения (68) имеет вид (см.приложение, а) )

v = vo -iv 0-vt)s n * (£ z j± L , 9 ) = Vf +(v«0-v ,)c/7 z (lil£ £ A ,

(83)

Здесь v

-

ордината

максимума волны -(горба),

vf-

ордината

минимума

(шадины),

(vg - vf )

- следовательно,

амплитуда вол­

ны. Модуль

^ имеет

значение

 

 

 

 

 

 

9

vО

-V2

 

(84)

 

 

 

 

 

 

24

где

 

 

 

 

 

 

 

 

 

 

 

Ь

= -

^ r~ ±

J

\J(Vc +v0) 2 -4V0 (vc + V0)

(85)

 

 

 

 

 

 

WB> V >

V, > V2 .

 

 

Длина

волны

h

определяется

выражением

 

 

 

 

 

 

h =

Z 'JJc

/ Zf(v0-v2) '

 

( 86)

 

Если амплитуда

кноидальной

 

 

волны

близка

к амплитуде

солифо-

 

 

на, волновой процесс представляет

 

 

собой

цепочку периодически

сле-г

 

 

дующих друг

за

другом

импульсов,

 

 

форма

которых близка

к

солифону

 

 

(р и с .4,

кривая

а ).

Чем ближе

ам­

Р и с . 4

 

плитуда процесса к солифовной.тем

 

 

 

дальше импульсы друг

от друга. В пределе расстояние‘ между

им­

пульсами увеличивается до бесконечности и кноидальная волна вы­

рождается в уединенную (кривая с ) .

В пределе малых амплитуд влияние нелинейности становится

пренебрежимо малым и процесс вырождается в монохроматический.

 

Обе эти тенденции в нелинейной эволшди кноидальных

вол',

нетрудно проследить

аналитически.

Действительно, при

прибли

кении

амплитуды кноидальной

волны

(83) к солифоннрй ( v0~ v, )

vc

имеет место

v' — •vz— »—

Пользуясь

фор­

мулой

вырождения для

эллиптического синуса

 

 

 

 

 

 

 

 

 

(8'?)

из (83)

находим

 

 

 

 

 

 

V ( x - c t ) /

 

 

 

 

(8-3;

 

 

h

 

 

 

 

 

При уменьшении амплитуды' волны модуль стремится

к

нулю

 

 

0 ) и формулы

 

 

 

 

 

 

 

Sin U’

Cn{U’Z)^ - ^ o ^ CaSU'

 

(89)

дают гармоническую волну (кривая Ь ):

Зак.240

25

V (x - ct) = v0 c os

,

Zc

(90)

к = 'fwJFifF

 

 

 

На р и с.5 изображены профили кноидальиых волн для

некоторых

значений параметров ( Vp- V,), h .

 

 

(v, v ) по-

Вид решений уравнения (68)

 

на фазовой плоскости

Р и о_ 5 Р и с . 6

казал на р и с.6 . Фазовая картина решений при отрицательном коэф­ фициенте у второго члена этого уравнения (быстрый пооцесс) сме­

стится на

величину 2 vc / 3 .

 

 

Учет

даооипацм приводит к уравнению (77),

численные реше­

ния которого для следующего вабора коэффициентов; V = Orf;

f = fi

0.4-

приведены на рио.7 . Кривые й ,Ь t C,

как видно,

пред­

ставляют затухающие колебания квавикноидального тала, которые о

уменьшением дисперсии принимают форму ударной волны

с

осцилли­

рующей структурой [16,44,46], вырождаясь

в пределе

V

О в

апериодический скачок. Между периодическим

с и неустойчивым d

решениями лежит другая предельная

апериодическая волна, пример­

ный профиль которой показан на р и с.8 (кривая а

) . Для сравне­

ния здесь же приведена уединенная

волна ( Ъ ),

соответствующая

бездиссяпативиому приближению. Помимо плавниковой деформации от солифона это решение отличается "расщепленностью" асимптотик:

йб

 

 

Р и G. 7

vpr — o

o

Следует, однако, иметь в виду, что

образование расщепленного солифона (или, в общем, расщепленного солитона) имеет место на бесконечно большом расстоянии от излу­

чателя.

Б связи

с

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

рас­

сматривать квазистационарную волну

б и с о л и т о н н о г о

типа -

спаренную комбинацию горба

 

и впадины

(р и с.8, кривая

с ,

фазовая

картина на р и с.9 ). Хотя

"бисолитон”

уже не является,

строго говоря, стационарной волной, однако его отклонение

от

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

может

быть весьма

мало и сказываться поэтому

толь-

V

Р и с . 8

Р и с . 9

27

кс) аа очень больших расстояниях.

В гемо-

и нейродинамике

(гл.Ш

а У1), где проводящие системы не

являются

в этом смысле

без­

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

вполне оправдано.

 

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

§ 5. Простая волна и интенсивный солифон

Выход за пределы солифонной сепаратрисы во внешнюю область

незамкнутых фазовых траекторий,

дающих неустойчивые решения

уравнения (68 ), сопровождается

нарушением стационарности уеди­

ненной волны из-за нарастания нелинейных тенденций,

что приво­

дит к деформации ее профиля.

 

 

Нелинейная деформация волновой геометрии была рассмотрена,

как известно, еще Риманом [ 4б]

на примере простой

волны - ин­

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

жидкости.

Приведем римановское решение о целью сохранить целы ость изло­ жения.

Очевидные "манипуляции" позволяют представить уравнения Эйлера и непрерывности в виде

 

 

 

 

(92)

Умножив

первое уравнение

на

 

и вычтя иэ него вто­

рое, получим для местной скорости

волны с ( р ) :

 

 

 

 

(93)

Подставив это

значение скорости

в

(91),

найдем

 

 

 

 

(94)

28

Уравнение (94) имеет очевидные решения

(95)

представляющие бегущие в противоположных направлениях ( ± л ; прогрессивные волны, которые и именуются "простыми". Как видно,

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

ct = ± с + v .

(96)

Отбросив в уравнении (9) члены, содержащие у3” и р , не­ трудно получить для местной скорости выражение

~ Со \ 1

(97)

3/Ъ

 

Если же, следуя традиции, воспользоваться адиабатическим уравнением состояния Пуассона - Тэта

- - P J 4 - Y ,

(98)

ГР >

 

и

то из (93) вытекает

 

г-t

С“ ± с.

Г-1

(99)

А

= * Со + ~ 2

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

 

 

с. = * сл +

-Ш - v ,

 

 

 

( 100)

то

нетрудно видеть, что

точки профиля волны, для которых

ср и

v

совпадают по знаку,

перемещаются с бо'лыпей скоростью,

чем

нули ( у =

0 ), а точки,

для которых эти знаки противоположны,-

с

меньшей.

Это обстоятельство приводит к одному из

важнейших

эффектов самовоздействия

нелинейных волн -

укручению

передне­

го фронта волны, следствием чего становится

образование разры­

ва

(ударной волны). Если речь идет

о первоначально синусоидаль­

ной волнеу

то нелинейная

деформация

ее профиля приводит к

пи­

лообразной геометрии, которая экспериментально наблюдалась ря­ дом исследователей [47,48].

Имея в виду (100) и переопределяя колебательную скорость

Соседние файлы в папке книги из ГПНТБ