Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/~math/magngidrwaves/node6.html
Дата изменения: Fri Apr 25 19:23:49 2003
Дата индексирования: Mon Oct 1 23:00:24 2012
Кодировка: koi8-r
Поисковые слова: protoplanetary disk

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

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

$$div \vec{H}=0$$ (38)

$$\vec{E}=-\frac{1}{c}[\vec{v} \vec{H}]+\frac{c}{4\pi\sigma}rot \vec{H}$$ (39)

$$\frac{\partial{\vec{H}}}{\partial t}=rot [\vec{v} \vec{H}]+\frac{c^2}{4\pi\sigma}\Delta\vec{H}$$ (40)

Также,напоминая, используя уравнение непрерывности, из уравнения (40) можно получить другое просто полезное уравнение:

$$\frac{d}{dt}\left(\frac{\vec{H}}{\rho}\right)=\left(\frac{\vec{H}}{\rho} \vec{\nabla}\right)\vec{v} +\frac{c^2}{4\pi\sigma}\Delta\vec{H}$$ (41)

Наконец, можно написать выражение для величины объемной силы действующей на жидкость - это сила Лоренца. Собственно, так как объемных зарядов нет, то объемная сила

$$\vec{f}=\frac{1}{c}[\vec{j} \vec{H}]$$ (42)

Соотвественно так как

$$j=\sigma \vec{E'}=\sigma (\vec{E}+\frac{1}{c}[\vec{v} \vec{H}])$$ (43)

подставляя $\vec{E}$ из уравнения (39) получим, что

$$\vec{f}=-\frac{1}{4 \pi}[\vec{H} rot \vec{H}]$$ (44)

Теперь немного снова об уравнениях гидродинамики. Если рассматривать не идеальную жидкость, а жидкость с вязкостью, то мы будем иметь не уравнение Эйлера, а уравнение Навье-Стокса. Напомним откуда оно собственно берется: При введении в рассмотрение вязких сил тензор натяжений перестает быть диагональным, и к тензору натяжений, в первом приближении, добавляются члены линейные по градиенту скорости, и соотвественно учитывая симметричность этого тензора оказывается, что общий вид вязкого тензора напряжений зависит от двух констант, т. н. первой и второй вязкостей.

$$p_{ik}=-p \delta_{ik}+\eta \left(\frac{\partial v_k}{\partial\... ...s}{\partial x_s}\right) +\zeta \delta_{ik}\frac{\partial v_s}{\partial x_s}$$ (45)

Для такого тензора напряжений вместо уравнения Эйлера получается соотвественно следующее уравнение (Навье-Стокса):

$$\rho\frac{d\vec{v}}{dt}=-\vec{\nabla} p+f +\eta\Delta\vec{v}+(\zeta+\frac{\eta}{3})\vec{\nabla}(\vec{\nabla}\vec{v})$$ (46)

Для рассмотрения магнитогидродинамической волны в жидкости нам надо рассмотреть малые добавки к стационарным значениям параметров в жидкости, таких как H, $\rho$, $\vec{v}$, p.

$$\vec{H}(\vec{r},t)=\vec{H_0}+\vec{\delta H}(\vec{r},t) p(\vec{r},t)=p_0+\delta p(\vec{r},t) \rho(\vec{r},t)=\rho_0+\delta\rho(\vec{r},t) \vec{v}(\vec{r},t) = \vec{v}_0 + \vec{\delta v}(\vec{r},t)$$ (47)

Причем мы будем рассматривать волны в жидкости при условии отсутствия макроскопических потоков в жидкости. Таким образом стационарное значение значение скорости принято равным нулю $\vec{v}_0=\vec{0}$. Для других параметров условие малости возмущений запишется следующим образом

$$\delta\rho\ll\rho_0 \delta \vec{H}\ll \vec{H}_0 \delta p \ll p_0$$ (48)

Макроскопически магнитное поле будем считать однородным, жидкость однородной. Для получения уравнений распространения возмущений в среде нам очевидно должно понадобится уравнение состояния, вообще говоря выглядящее вот так $p=p(\rho,s)$, где s - энтропия.

Дальше, можно собственно, либо сказать, что рассматриваем например адиабатический процесс и тогда $p=p(\rho)=p_0 \left(\frac{\rho}{\rho_0}\right)^\gamma$. Можно же сказать, что так как мы рассматриваем малые возмущения, то учитывая то, что $\rho T\frac{ds}{dt}=\mathcal{D}$, где $\mathcal{D}$ - диссипативная функция, которая квадратична по скоростям - т.е. в нашем приближении ей можно пренебречь, а следовательно считать s постоянной, и таким образом иметь опять же уравнение состояния в виде $p=p(\rho)$ Учитывая то, что мы имеем уравнение состояния $p=p(\rho)$ , то раскладывая получим:

$$\delta p=u_0^2 \delta \rho u_0^2 = \left.\left(\frac{\partial p}{\partial \rho}\right)\right\vert _{\rho=\rho_0}$$ (49)

где через $u_0$ обозначена скорость звука в жидкости. Теперь, для дальнейших вычислений заменим субстанциональные производные на эйлеровы производные. Дело в том ,что, как известно, субстанциональная производная выражается через Эйлерову производную следующим образом: $\frac{d}{dt}=\frac{\partial}{\partial t} + (\vec{v} \nabla)$. А в предположении малых возмущении мы просто получаем равенство субстанциональной и Эйлеровой производных. Наконец напишем несколько очевидных равенств верных с точностью до бесконечно малых более высокого порядка

$$[\vec{H} rot \vec{H}]=[H_0 rot \vec{\delta H}] [\vec{v}\vec{H}]=[\vec{\delta v} \vec{H_0}]$$ (50)

Соответственно, тогда основные уравнения запишутся следующим образом:

$$\rho_0\frac{\partial\vec{v}}{\partial t}= -u_0^2 \vec{\nabla}\del... ...}]+\eta\Delta\vec{v}+ (\zeta+\frac{\eta}{3})\vec{\nabla}(\vec{\nabla} \vec{v})$$ (51)

$$\frac{\partial{\vec{\delta H}}}{\partial t}= rot [\delta\vec{v} \vec{H}_0]+ \frac{c^2}{4 \pi \sigma}\Delta \vec{\delta H}$$ (52)

$$\frac{\partial\delta \rho}{\partial t}=-\rho_0 (\vec{\nabla} \vec{v})$$ (53)

$$\vec{E}=-\frac{1}{c}[\vec{\delta v} \vec{H}_0]+\frac{c}{4\pi\sigma}rot \vec{\delta H}$$ (54)

Наша задача сейчас - получить волновое уравнение для возмущений в жидкости. Если такое волновое уравнение будет получено, то это будет автоматически означать существование волн. Однако получить его в самом общем виде сложно, если вообще возможно. Поэтому оно будет получено собственно для плоских волн. Но для начала преобразуем уравнения (51)- (53). Как всегда бывает, для получения волнового уравнения надо к одному уравнению применить дифференциальный оператор по координатам, а другое уравнение продифференцировать по времени. Так и сделаем. Уравнение (51) продифференцируем по времени, а к уравнению (53) применим оператор градиента и после этого подставим в продифференцированное (51) и поделим на $\rho_0$. Тогда получим:

$$\ddot{\vec{v}}-u_0^2 \vec{\nabla}(\vec{\nabla} \vec{\delta v})+ ... ...\frac{1}{\rho_0}(\zeta+\frac{\eta}{3})\vec{\nabla}(\vec{\nabla} \dot{\vec{v}})$$ (55)

Теперь собственно рассмотрим плоские волны:

$$\delta \rho(\vec{r},t) = \rho e^{-i\omega t+i\vec{k}\vec{r}} \delta \vec{H}(\vec{r},t) = \vec{h} e^{-i\omega t+i\vec{k}\vec{r}} \delta \vec{v}(\vec{r},t) = \vec{v} e^{-i\omega t+i\vec{k}\vec{r}}$$ (56)

Здесь, $\rho, \vec{h}, \vec{v}$ - постоянные величины , являющиеся амплитудами соответствующих колебаний. Теперь, учитывая то, что

$$rot(\vec{A} b)=b rot \vec{A} + \vec{A} \times grad b div(\vec{A} b)=b div \vec{A} + (\vec{A} grad b)$$ (57)

$$grad (e^{-i\omega t+i\vec{k}\vec{r}})=i \vec{k} e^{-i\omega t+i\vec{k}\vec{r}} \frac{\partial}{\partial t} e^{-i\omega t+i\vec{k}\vec{r}} = -i \omega e^{-i\omega t+i\vec{k}\vec{r}}$$ (58)

уравнение (55) преобразуется в

$$-\omega^2 \vec{v}+u_0^2 \vec{k} (\vec{k} \vec{v})+ \frac{\ome... ...i \omega}{\rho_0}\left(\zeta+\frac{\eta}{3}\right)\vec{k} (\vec{k} \vec{v})=0$$ (59)

Чтобы получить волновое уравнение нам надо избавиться в последнем уравнении от еще одной независимой переменной и оставить только скорость. Для этого подставляя уравнения (56) в уравнение (52). Тогда получим:

$$\left(-i \omega+\frac{c^2 k^2}{4 \pi \sigma}\right)\vec{h}= i [\vec{k}[\vec{v} \vec{H}_0]] \Rightarrow \vec{h}=-\frac{1}{\omega} \frac{1}{1+i \frac{c^2 k^2}{4 \pi \sigma \omega}} \left[\vec{k}\left[\vec{v} \vec{H}_0\right]\right]$$ (60)

Далее, подставим последнее выражение в (59), обозначив при этом $\vec{V}_a=\frac{1}{\sqrt{4 \pi \rho_0}}\vec{H}_0$ (по размерности это скорость) получим, приводя подобные члены:

$$-\omega^2\left(1+\frac{i \eta}{\omega \rho_0}\vec{k}^2\right)\v... ...left[\vec{k}\left[\vec{k}\left[\vec{v} \vec{V}_a\right]\right]\right]\right]=0$$ (61)

Так собственно мы получили собственно фактически волновое уравнение для скорости в достаточно общих предположениях (малые возмущения в вязкой, неидеально проводящей жидкости в однородном магнитном поле). Получим также несколько вспомогательных уравнений. Подставляя  (56) в (53) получим:

$$\delta\rho=\frac{\rho_0}{\omega}(\vec{k} \vec{v})$$ (62)

Подставляя же (56) и (60) в (54) получим:

$$\vec{E}=-\frac{1}{c}[\vec{v} \vec{H}_0]-\frac{i c}{4 \pi \sig... ...ega}} \left[\vec{k}\left[\vec{k}\left[\vec{v} \vec{H}_0\right]\right]\right]$$ (63)

Теперь же наконец можно все-таки перейти к тем приближениям, о которых говорилось до этого: т. е. приближения идеальной и идеально проводящей жидкости. Эти приближения естественно сводятся к следующему: $\eta \to 0, \zeta \to 0$ и $\sigma \to \infty$ Тогда наш основной набор уравнений будет выглядеть следующим образом:

$$-\omega^2\vec{v}+u_0^2\vec{k}(\vec{k} \vec{v})- \left[\vec{V}_a\left[\vec{k}\left[\vec{k}\left[\vec{v} \vec{V}_a\right]\right]\right]\right]=0$$ (64)

$$\vec{h}=-\frac{1}{\omega}[\vec{k}[\vec{v} \vec{H}_0]] \vec{E}=-\frac{1}{c}[\vec{v} \vec{H}_0] \rho=\frac{\rho_0}{\omega}(\vec{k} \vec{v})$$ (65)

Раскроем векторное произведение

$$\left[\vec{V}_a\left[\vec{k}\left[\vec{k}\left[\vec{v} \vec{V}_a... ...ight]- (\vec{k} \vec{v})\left[\vec{V}_a\left[\vec{k} \vec{V}_a\right]\right]=$$ (66)

$$= (\vec{k} \vec{V}_a)(\vec{k}(\vec{V}_a \vec{v})-\vec{v}(\vec{k... ...V}_a)) -(\vec{k} \vec{v})(\vec{k} \vec{V}_a^2-\vec{V}_a(\vec{k} \vec{V}_a))=$$ (67)

$$= (\vec{k} \vec{V}_a)(\vec{k}(\vec{V}_a \vec{v})-\vec{v}(\vec{k} \vec{V}_a)+ \vec{V}_a(\vec{k} \vec{v}))-(\vec{k} \vec{v})\vec{k}\vec{V}_a^2$$

Тогда уравнение (64) запишется следующим образом:

$$-\omega^2\vec{v}+(u_0^2+V_a^2)\vec{k}(\vec{k} \vec{v})+\left( -\... ...}(\vec{k} \vec{V}_a)- \vec{V}_a(\vec{k} \vec{v})\right)(\vec{k} \vec{V}_a)=0$$ (68)

И называется это уравнение Основным спектральным уравнением. Теперь перейдем собственно к его решению. Очевидно, так как уравнение векторное, то мы имеем на самом деле однородную линейную систему из трех уравнений, которая имеет решение только при некотором условиях на коэффициенты (эти условия на самом деле будут дисперсионными уравнениями для соответствующих волн (т.е. соотношениями между волновым вектором и частотой)). Распишем для получения этих условий уравнение (68) по осям. Причем для этого ориентируем систему координат так, чтобы вектор $\vec{k}$ был направлен вдоль оси Y, а вектор $\vec{H}_0$ лежал в плоскости YZ (очевидно, что мы всегда можем так выбрать систему координат). Итак тогда:

$$\vec{k}=k \vec{e}_2 \vec{H}_0=H_y \vec{e}_2+H_z \vec{e}_3$$ (69)

Итак, уравнение (68) запишется:

$$(-\omega^2+k^2V_y^2)\vec{v}+(u_0^2+V^2)k^2 v_y\vec{e}_2 -k^2 V_y(v_y V_y+v_z V_z)\vec{e}_2 -k^2V_yv_y(V_y\vec{e}_2+V_z\vec{e}_3)=0$$ (70)

Распишем последнее уравнение по координатам:

$$\begin{align*}\begin{cases}(-\omega^2+k^2 V_y^2)v_x=0 \ (-\omega^2+k^2 V_y^2+... ...2 V_yV_zv_z=0\ -k^2V_yV_zv_y+(-\omega^2+k^2V_y^2)v_z=0 \end{cases}\end{align*}$$ (71)

Видно что в системе уравнения делятся, т.е. есть уравнение содержащее только $v_x$, а другая пара уравнений содержит $v_y$ с $v_z$. Это значит на самом деле,что в жидкости могут распространяться независимо!!! друг от друга два типа возмущений. Рассмотрим первый тип возмущений, происходящий из первого уравнения системы. Он соответствует $v_x\ne0$. Для этого типа волн (называемых Альвеновскими) мы получаем дисперсионное уравнение:

$$\omega^2=k^2 V_y^2$$ (72)

или более общо:

$$\omega^2=(\vec{k} \vec{V})^2$$ (73)

А отсюда

$$v_{gr}=\frac{\partial\omega}{\partial\vec{k}}=\vec{V}_a=\frac{1}{\sqrt{4 \pi \rho_0}}\vec{H} v_{ph}=\frac{\omega}{k}=\frac{1}{\sqrt{4 \pi \rho_0}}H_y$$ (74)

Как нетрудно видеть теперь мы получили поперечную волну, причем изменений плотности в ней нет и она следовательно возможна как в сжимаемой так и в несжимаемой жидкости. Энергия в этой волне ( в соответствии с направлением групповой скорости) распространяется вдоль магнитного поля.

Рассмотрим теперь последние два уравнения из (71) Опять же распишем условие существования нетривиального решения у двух линейных уравнений - равенство нулю определителя. Получится следующее уравнение:

$$\omega^4-\omega^2k^2(u_0^2+V^2)+k^4((u_0^2+V_z^2)V_y^2-V_y^2V_z^2)=0$$ (75)

Решая это биквадратное уравнение, получим:

$$\omega=k \sqrt{\frac{1}{2} \left(u_0^2+V_a^2 \pm\sqrt{\left(u_0^2+V_a^2\right)^2-4 V_y^2 u_0^2}\right)}$$ (76)

$$\omega=k \sqrt{\frac{1}{2} \left(u_0^2+\frac{\vec{H}_{0}^2}{4 ... ...H}_{0}^2}{4 \pi \rho_0}\right)^2-\frac{H_{0y}^2}{\pi \rho_0} u_0^2}\right)}$$ (77)

Эти волны уже не являются ни чисто продольными, ни чисто поперечными, хотя и называются магнитозвуковыми (медленной магнитозвуковой ,если в последней формуле выбран знак минус, и быстрой, если наоборот). Называются они магнитозвуковыми, так как в приближении слабого магнитного поля $V_a^2«u_0^2$ видно что скорость быстрой магнитозвуковой волны будет равна скорости звука, и из уравнения (68) видно, что волна получится продольная.