Морозов Введение в теорию горячей плазмы Ч.1 2011
.pdf∂v |
|
∂v |
v cs 0. |
(VIII.2.8) |
|
∂t |
∂x |
||||
|
|
|
Это уравнение очень похоже на уравнение (VIII.1.2), описывающее пучок невзаимодействующих частиц.
Перейдем в систему отсчета, движущуюся со
скоростью v cs или v −cs; v′ v cs. Тогда вместо (VIII.2.8) имеем:
∂v′ |
v′ |
∂v′ |
0. |
(VIII.2.9) |
|
∂t |
∂x |
||||
|
|
|
Теперь оно в точности совпадает с уравнением (VIII.1.2), описывающим поток невзаимодействующих частиц. Следовательно, на нелинейной стадии также должно происходить укручение фронта ионно-звуковой волны, переходящее в ударную волну при ∂∂vx → . На
фронте такой волны, как и на фронте всякой ударной волны, происходит диссипация механической энергии вследствие вязкости.
VIII.3. Нелинейные МГД-волны
Рассмотрим поперечную МГД-волну cos 0 , т.е. магнитный звук, в приближении 8H2P 1.
Уравнения, описывающие такую волну, имеют вид:
181
∂v |
v |
∂v |
|
|
|
|
|
|
|
1 |
|
∂ |
|
H2 |
0, |
(VIII.3.1) |
||||||
∂t |
∂x |
|
|
|
min |
∂x 8 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
∂H |
|
|
|
|
∂ |
|
|
v H |
0, |
|
(VIII.3.2) |
|||||||||
|
|
|
∂t |
∂x |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
∂n |
|
|
|
|
∂ |
|
nv 0. |
|
|
(VIII.3.3) |
|||||||||
|
|
|
∂t |
|
∂x |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из последних двух уравнений видно, что
H2 |
n2. |
(VIII.3.4) |
|
8 |
|||
|
|
Далее можно проделать все те же выкладки, что и для ионного звука. Следовательно, также должно наблюдаться укручение и "опрокидывание" фронта.
Что касается альфвеновской волны, для нее v k, поэтому v, v 0, и при больших может происходить распространение нелинейной альфвеновской волны без искажения. Однако при1 альфвеновская волна также опрокидывается из-за перекоса силовых линий.
VIII.4 Волны на мелкой воде
Рассмотрим еще один вид нелинейных волн, имеющих важное значение не только в физике плазмы, но и в геофизике – это волны на мелкой воде.
Пусть на твердую поверхность налит слой воды
182
толщиной h0. Пусть также по этому слою воды распространяется одномерная волна в направлении x, причем длина волны много больше глубины слоя h0, h0 (см. рис.26)
Ввиду малой толщины слоя можно считать, что жидкость смещается только в направлениях x, а скорость этого смещения не зависит от высоты. Жидкость будем считать несжимаемой, const. Давление жидкости определяется гидростатическим законом:
P g h − h0 . |
(VIII.4.1) |
Уравнения гидродинамики, описывающее такую волну, принимают вид:
∂v |
v |
∂v |
g |
∂h |
0. |
(VIII.4.2) |
|
∂t |
∂x |
∂x |
|||||
|
|
|
|
Уравнение непрерывности теперь перепишется так:
∂h |
v |
∂ |
vh 0. |
(VIII.4.3) |
|
∂t |
∂x |
||||
|
|
|
Рис 26. Волна на мелкой воде
183
Будем считать, что высота слоя и скорость связаны однозначно: h h v . Тогда вместо уравнений (VIII.4.2) и (VIII.4.3) можно написать:
|
|
∂v |
v |
∂v |
g dh |
∂v |
0, |
(VIII.4.4) |
||||||
|
|
|
|
|
|
|
|
|||||||
|
|
∂t |
|
∂x |
dv ∂x |
|
|
|||||||
∂h |
|
|
∂v |
v |
∂v |
h |
∂v |
|
0. |
(VIII.4.5) |
||||
∂t |
∂t |
∂x |
|
|||||||||||
|
|
|
|
|
∂x |
|
|
Умножая (VIII.4.5) на dhdv и вычитая из него
(VIII.4.4), получаем:
h ∂∂vx g dhdv 2 ∂∂vx ,
или
∂∂hv hg
и
∂v |
v |
gh |
∂v |
0. |
(VIII.4.6) |
|
∂t |
∂x |
|||||
|
|
|
|
Таким образом, и в этом случае получаем результат, аналогичный результату для ионного звука.
Полагая h h′ h0, где h′ h0, получаем:
h ≈ h0 |
dh |
|
h0 |
|
|
dv v h0 |
v |
g |
. |
(VIII.4.7) |
|
С другой стороны |
|
|
|
|
|
184
|
g h |
≈ |
g h0 |
|
1 |
|
|
|
g |
|
|
dh |
|
1 |
v , |
|
(VIII.4.8) |
||||||
|
2 |
|
|
h0 |
|
dv |
2 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
и вместо (VIII.4.6) получаем: |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
∂v |
|
3 |
v c0 |
∂v |
|
0 , |
|
|
|
(VIII.4.9) |
||||||||||
|
|
|
∂t |
|
∂x |
|
|
|
|
||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где |
c0 |
|
gh0 . |
Переходя |
|
в |
|
систему, |
|||||||||||||||
движущуюся |
со скоростью |
c0 |
|
и |
|
вводя u |
3 |
v, |
|||||||||||||||
|
|
2 |
|||||||||||||||||||||
снова получаем уравнение |
|
∂u |
u |
|
∂u |
|
0 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
∂t |
|
∂x |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cо всеми вытекающими из него следствиями. Заметим, что во всех рассмотренных случаях не
учитывалась дисперсия, т.е. ∂∂k const. Учет даже
слабой дисперсии приводит к важным физическим результатам.
VIII.5. Волны в средах со слабой дисперсией. Уравнение Кортевега-де Вриза. Солитоны
Учтем малую дисперсию при распространении ионного звука. В линейной задаче для этого типа волн был получен следующий закон дисперсии:
k2r2
2 2 De .
pi 1 k2r2De
При krDe 1 дисперсия отсутствовала:≈ pikrDe, и групповая скорость не зависела от k,
185
т.е. vqr pirDe mTei .
Однако учет следующего малого члена в разложении выражения (VIII.5.1) приводит к появлению дисперсии:
|
|
|
1 |
|
|
2 |
|
2 |
|
|
||
≈ pikrDe 1 − |
|
|
|
k |
|
|
rDe |
, |
(VIII.5.2) |
|||
2 |
|
|
||||||||||
vph |
Te |
|
1 |
|
|
|
2 |
2 |
|
|
||
mi |
1 − |
|
|
k |
|
|
/k0 |
. |
(VIII.5.3) |
|||
|
2 |
|
|
Здесь k02 |
1 |
. Величина k0 определяет длину, |
2 |
||
|
3rDe |
на которой изменение групповой скорости оказывается порядка ей самой.
Учтем малую дисперсию в нелинейном уравнении для простых волн. Для этого перейдем в
систему координат, движущуюся |
со скоростью |
||||
c0 pikrDe. |
Линейное |
|
дисперсионное |
||
соотношение (VIII.5.2) преобразуется в следующее: |
|||||
|
|
c0k3 |
|
0. |
(VIII.5.4) |
|
2k02 |
||||
|
|
|
|
Выполняя обратное фурье-преобразование,
получаем линейное уравнение с учетом дисперсии: |
|
||||||
|
∂u |
|
cs |
|
∂3u |
0. |
(VIII.5.5) |
|
∂t |
2k02 |
|
∂x3 |
|||
|
|
|
|
|
Но при конечной амплитуде возмущения
уравнение должно содержать член u ∂∂ux . Таким
образом, полное уравнение должно содержать и нелинейный член, и член, описывающий
186
дисперсию: |
|
|
|
|
|
|
|
|
|
|
∂u |
u |
∂u |
|
cs |
|
∂3u |
0. |
(VIII.5.6) |
|
∂t |
∂x |
2k02 |
|
∂x3 |
||||
|
|
|
|
|
|
Это уравнение носит название "Кортевега-де Вриза" и было получено в 1895 г.
Уравнение Кортевега-де Вриза допускает новый тип решений – уединенную волну (солитон).
Для того, чтобы найти такое решение, перейдем в систему отсчета, движущуюся с некоторой, пока неопределенной, скоростью c и введем новую
переменную x − ct. |
Предположим также, |
что |
|||||||||||||||
скорость u зависит только от |
|
одной |
переменной |
||||||||||||||
: u u , т.е. в новой системе координат имеем |
|||||||||||||||||
уравнение: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
du |
u |
du |
|
|
cs |
|
∂3u |
0. |
(VIII.5.7) |
||||||
|
d |
d |
|
2k02 |
|
∂3 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
Это |
уравнение |
|
|
|
можно |
|
один |
раз |
|||||||||
проинтегрировать по : |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
cs |
|
|
∂2u |
|
a cu − |
|
u2 |
. |
(VIII.5.8) |
|||||
|
|
|
2k02 |
|
∂ 2 |
|
2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Здесь a – постоянная интегрирования. Прибавив к скорости u постоянную величину, можно добиться того, чтобы a обратилась в ноль.
Полученное уравнение в точности совпадает с уравнением, описывающим движение классической материальной точки в поле сил с потенциалом
187
w − |
cu |
2 |
|
u3 |
. |
(VIII.5.9) |
|
|
6 |
||||
2 |
|
|
|
|
Потенцитал представлен на рис.27.
Масса частицы при этом равна m
переменная соответствует времени, а переменная u – отклонению от точки положения равновесия. Точка u 2c соответствует положению равновесия, а малые колебания вблизи этого положения являются линейными:
u 2c u, |
|
|
2c |
x − ct . |
(VIII.5.10) |
u 2c u0 exp ik0 c0 |
Величина u0 здесь соответствует амплитуде колебаний, а решение (VIII.8.5) – периодической волне, распространяющейся со скоростью c. С ростом амплитуды колебания все более отличаются от гармонических. Если u0 достигает величины 3c, решение переходит в уединенную волну, показанную на рис. 28.
Решение при этом имеет вид:
u |
|
u0 |
|
. |
(VIII.5.11) |
|
ch2 |
|
−ct |
|
|||
|
|
|
|
|||
|
|
|
Такое решение называется уединенной волной
(solitary wave), или солитоном.
188
Рис. 27. Эффективный потенциал
Рис. 28. Солитон
Подставим решение (VIII.5.11) в уравнение
189
(VIII.5.8), введя переменную |
y |
x−ct |
. |
Это |
||
|
||||||
уравнение удовлетворяется, если u0 |
3c, т.е. |
|
||||
2u0 |
6cs |
const. |
|
|
(VIII.5.12) |
|
|
|
|
||||
|
k02 |
|
|
|
|
Выбрав амплитуду и параметр 2u0, можно найти ширину солитона.
В среде с отрицательной дисперсией групповая скорость vgr v0gr 1 − 2kk22 . Неопределенная ранее
0
скорость (в системе отсчета, двигающейся со скоростью c0) равна
c |
u0 |
0 , |
(VIII.5.13) |
|
3 |
||||
|
|
|
т.е. солитон распространяется со скоростью больше скорости звука.
Полученные решения имеют интересные приложения в геофизике. С ними связывают иногда встречающиеся в океане одинокие "волны-убийцы", появляющиеся неожиданно и имеющие громадную высоту.
190