CC BY
10
УДК 532.591
С.И. ПЕРЕГУДИН, д-р физ. -мат. наук, профессор, peregudinsi@yandex. ru Санкт-Петербургский государственный горный институт (технический университет) С.Е. ХОЛОДОВА, канд. физ. -мат. наук, доцент, kholodovase@yandex. ru Санкт-Петербургский государственный университет
S.I. PEREGUDIN, Dr. in рhys. and math. sc., professor, [email protected] Saint Peterburg State Mining Institute (Technical University)
S.E. KHOLODOVA, PhD in рhys. and math. sc., associate professor, [email protected] Saint Petersburg State University
ВОЛНОВЫЕ ДВИЖЕНИЯ, ВЫЗВАННЫЕ КОЛЕБАНИЯМИ ПЛОСКОЙ СТЕНКИ
Исследование посвящено построению аналитических, в частности, точных решений задачи о движении электропроводной идеальной вращающейся жидкости, моделирующей волновые движения в жидком ядре Земли. Рассматривается задача об излучении волн во вращающуюся несжимаемую электропроводную жидкость плоской стенкой, совершающей, начиная с начального момента, гармонические колебания. Такая постановка задачи может служить первым шагом к исследованию влияния неоднородности жидкой среды на генерацию и поддержание магнитного поля Земли.
Ключевые слова: идеальная электропроводная вращающаяся жидкость, волновые движения, магнитное поле Земли.
THE WAVE MOTIONS CAUSED BY OSCILLATIONS
OF A FLAT WALL
The research is devoted construction analytical, in particular, exact solution of a problem on movement of the electrically conducting ideal rotating liquid modeling wave movements in a liquid kernel of the Earth. Namely, the problem about radiation of waves in rotating incompressible electrically conducting liquid by the flat wall making, since the initial moment, harmonious fluctuations is considered. Such statement of a problem can serve as the first step to research of influence of heterogeneity of the liquid environment on generation and maintenance of a magnetic field of the Earth.
Key words: the electrically conducting ideal rotating liquid, wave motions, the magnetic field of the Earth.
Земное ядро оказывает существенное влияние на многие геофизические явления и процессы, происходящие и происходившие в Земле. Согласно общепризнанным представлениям, магнитное поле Земли возбуждается движением в жидкой части земного ядра, но детали этого процесса все еще не выяснены.
Изучение динамики глубоких недр Земли является важнейшей и вместе с тем труднейшей проблемой, которая решается, прежде всего, методами геофизики с использованием сейсмических и гравитационных методов. Непосредственно проникнуть в земные недра трудно, в связи с чем большое значение приобретают теоретические исследования.
Вопросам о движении электропроводной жидкости посвящено множество работ. Основополагающие из них выполнены У. Эльзассером [8, 9], Э. Буллардом [5], которым удалось доказать существование стационарных и нестационарных решений уравнения магнитной индукции при задании некоторого специального вида поля скоростей.
Вопросам о крупномасштабных движениях электропроводной жидкости посвящены работы [6, 7, 10-12], в которых рассматривалась модель, построенная в приближении быстрого вращения. В рамках этой теории в уравнении движения силой инерции пренебрегается и инерциальные, альфвенов-ские волны и волны Россби отфильтровываются. Кроме того, в пределе быстрого вращения скорость V определяется с точностью до слагаемого, представляющего собой геострофическую скорость. Последнее обстоятельство обусловлено тем, что геострофическая скорость не удовлетворяет магни-тострофическому уравнению. Для преодоления указанных трудностей привлекаются вязкие силы и пренебрегается вязкостью, когда это допустимо.
В работах [6, 7] анализировалась соответствующая задача в слое, заключенном между плоскостями г = 0 е г = ё . Здесь мы предполагаем, что границы слоя представляют собой поверхности, изменяющиеся в пространстве и во времени. Кроме того, в уравнении движения учитываются инерционные силы.
Будем использовать прямоугольную декартову систему координат Охуг. Под объемной силой понимается вектор g, перпендикулярный поверхности г = 0 и направленный в сторону, противоположную вертикальной оси. Ось вращения жидкости совпадает с осью г, т.е. м> = кю, где м> - угловая скорость вращения Земли.
Рассмотрим вращающийся слой электропроводной идеальной несжимаемой жидкости, ограниченный сверху твердой непроницаемой поверхностью г = -Z(х, у) мантии Земли, а снизу поверхностью г = (х, у, t) твердого ядра Земли. Сис-
тема уравнений для описания движения невязкой электропроводящей несжимаемой вращающейся с угловой скоростью п> жидкости в магнитной гидродинамике в переменных Эйлера имеет вид [1-3]
дУ / \ Ур
div V = 0;--ь (у • У )у =---2^ х
дt р
х V - gz +——rot B х B ; (1)
цр
ЯП
— = rot(vXБ); divB = 0,
dt
где В - вектор магнитной индукции; v -скорость жидкости в системе координат, вращающейся с угловой скоростью w; p -давление; р - плотность; g - ускорение свободного падения; ц - магнитная проницаемость, предполагается постоянной.
В статье [4] была проведена редукция нелинейной векторной системы уравнений в частных производных, моделирующей возмущение в слое идеальной электропроводной жидкости, к скалярному уравнению для функции С, имеющему вид
(
\
а д , д
К^Г + bQy — дх ду
УдбА2С „ 2
2Ъ + 2а2-^ + а
дчб
д^
дч4
, 4 д2А2^
л
дч2
дх
+b
Q у
дУ,
ЦР
д^С „ 2 д^
дч4
- + 2а
дч2
4
+ а4^
bQQ - b(e) z Q zQ
(^P)2 HQ '
(2)
и доказано утверждение, что любое решение V(х, у, t), Ь(х, у, t), х, у, t) задачи о малых возмущениях в слое идеальной несжимаемой однородной электропроводной вращающейся жидкости, удовлетворяющее необходимым условиям гладкости, предста-вимо в виде
Ь(х, у, t) = црД (р2 + а2)~;
3
д
д
b
ISSN 0135-3500. Записки Горного института. Т.187
г(х, у, t) = - Р( (р,2 +а2 }г ;
ё
~ = Р1г; Ь = ДЬ ;
^ 1[ Р а
^У ^-а Р
~ Л гх + РК
Vг у + РА у
ГР
V ~ У
= РР, (р,2 + а2 }х
^црР12 - Р2 ацрР, ^ - ацрР, црР,2 - Р2
~(х, у,,) = Р, (р,2^2 + (аР2 }2)Ъ(х, у,,);
F = црР, (р(2 + а2)-Р2
0-
Р=. А+л А,
дх
"0 у
ду
где функция С = Р(Ъ - решение уравнения (2); Ь - вектор магнитной индукции поля, обусловленного волновым движением жидкости; Ь0 - вектор магнитной индукции невозмущенного однородного постоянного магнитного поля; Н0(х, у) - толщина жидкого слоя в состоянии покоя; Ъг0 (х, у,,) и Ь(0 - компоненты магнитного поля на границах слоя.
Верно также и обратное утверждение, что любое решение уравнения (2) порождает решение системы уравнений (1), моделирующей малые возмущения в тонком слое идеальной несжимаемой однородной электропроводной вращающейся жидкости, если построенные по приведенным выше формулам функции V(х, у,,), Ь(х, у,,), г(х, у,,) удовлетворяют в рассматриваемой области условиям гладкости.
Исследуем волновые движения малой
амплитуды в полуплоскости Я + = = ^ = (х, у) е R2 : у > о}, ограниченной горизонтальной стенкой у = 0 . Направим ось Оу параллельно Ь0 и, полагая значение поля на вертикальных границах слоя и функцию глубины слоя не зависящими от х и
введя обозначение = С , получим мате-
ду
матическую постановку данной задачи в виде начально-краевой задачи относительно функции С:
Ъ
0 у
6А2~ + 2а2дМ + а4 д2А
д,ь
(
дг4
д,4
д 2А 2 С
д,2
Ъу д2
цр ду2
д4А2С 2 д А2 4 -Р + 2а2-р + а4 С
42
д/2
(е)
ъ2 0 - Ъ
2 0 2 0
(цр}2 Н 0
— С(х, у, 0) = 0, к = 0,5, С(х, у,,) д,
= С(х, у,,) = 0(t) ехр^/'ю,), (4)
у=0
где 0(,) - функция Хевисайда.
Решение поставленной задачи зависит от одной пространственной координаты у и
от времени ,, т.е. С = С(у,,). В дальнейшем изложении знак тильды опущен.
С помощью преобразования Лапласа по времени ,, получим решение задачи (-)-(4) в явном виде:
(
С( у,,) = 0
, -
л/^Р
л (
Ъ
у
0 у
ехр
(
- /ю
,у
Ъ
0 у УУ
Из этого выражения следует, что гармонические колебания плоской стенки излучают в электропроводную жидкость горизонтально распространяющуюся нестационарную плоскую волну с волновым
числом к = Юд/цр / Ъ0у.
Изучим поведение решения С(у,,) при , ^ да. Так как
(
lim С( х,,) ехр^/'ю,) = ехр
. л/Цр
/ю--
Л
Ъ
у
0 у
то в рассматриваемой задаче существует режим установившихся колебаний с предельной амплитудой, описывающей установившуюся плоскую волну. Таким образом, проведенный анализ позволяет судить о предельном поведении решения при больших временах.
х
X
к
Г
х
X
ЛИТЕРАТУРА
1. Алешков Ю.З. Математическое моделирование физических процессов. СПб, 2001. 264 с.
2. ГунькоЮ.Ф. Электромагнитная газодинамика плазмы / Ю.Ф.Гунько, А.В.Норин, Б.В.Филиппов. СПб, 2003. 176 с.
3. Ландау Л.Д. Электродинамика сплошных сред / Л.Д.Ландау, Е.М.Лифшиц // Теоретическая физика. М., 1992. Т.8. 664 с.
4. Холодова С.Е. Динамика вращающегося слоя идеальной электропроводной несжимаемой жидкости // Журнал вычислительной математики и математической физики. 2008. Т.48. № 5. С.151-167.
5. Bullard E.C. The magnetic field withing the Earth // Proc. Roy. Soc. London Ser. A. 1949. Vol.197. P.433-453.
6. Bus.se F.H. A model of the geodynamo // Geophys. J. Roy. Asron. Soc. 1975. Vol.42. P.437-459.
7. Busse F.H. Generation of planetary by convection // Phys. Earth Planet. Inter. 1976. Vol.12. P.350-358.
8. Elzasser W.M. Induction effects in terrestrial magnetism. Part I: Theory // Phys. Rev. 1946. Vol.69. P.106.
9. Elzasser W.M. Induction effects in terrestrial magnetism, Part II: The secular variation // Phys. Rev. 1946. Vol.70. P.202.
10. Zhang K.-K., Busse F.H. Finite amplitude convection and magnetic field generation in a rotating spherical shell // Geophys. Astrophys. Fluid Dyn. 1988. Vol.44. P.33-54.
11. Zhang K.-K., Busse F.H. Convection driven mag-netohydrodynamic dynamos in rotating spherical shell // Geophys. Astrophys. Fluid Dyn. 1989. Vol.49. P.97-116.
12. Zhang K.-K., Busse F.H. Generation of magnetic fields by convection in a rotating spherical fluid shell of infinite Prandtl number // Phys. Earth Planet. Inter. 1990. Vol.59. P.208-222.
REFERENCES
1. Aleshkov Yu.Z. Mathematical modeling of physical processes. Saint Petersburg, 2001. 264 p.
2. Gun'ko Yu.P., Norin A.V., Philippov B.V. Electromagnetic gas-dynamic of the plasma. Saint Petersburg, 2003. 176 p.
3. Landau L.D., Liphshits E.M. Electrodynamics of the continuous mediums // Theoretical Physics. Moscow, 1992. Vol.8. 664 p.
4. Kholodova S.E. Dynamics of a rotating layer of an ideal electrically conducting incompressible fluid // Journal of computational mathematics and mathematical phisics. 2008. Vol.48. N 5. P.151-167.
5. Bullard E.C. The magnetic field withing the Earth // Proc. Roy. Soc. London Ser. A. 1949. Vol.197. P.433-453.
6. Busse F.H. A model of the geodynamo // Geophys. J. Roy. Asron. Soc. 1975. Vol.42. P.437-459.
7. Busse F.H. Generation of planetary by convection // Phys. Earth Planet. Inter. 1976. Vol.12. P.350-358.
8. Elzasser W.M. Induction effects in terrestrial magnetism. Part I: Theory // Phys. Rev. 1946. Vol.69. P.106.
9. Elzasser W.M. Induction effects in terrestrial magnetism, Part II: The secular variation // Phys. Rev. 1946. Vol.70. P.202.
10. Zhang K.-K., Busse F.H. Finite amplitude convection and magnetic field generation in a rotating spherical shell // Geophys. Astrophys. Fluid Dyn. 1988. Vol.44. P.33-54.
11. Zhang K.-K., Busse F.H. Convection driven mag-netohydrodynamic dynamos in rotating spherical shell // Geophys. Astrophys. Fluid Dyn. 1989. Vol.49. P.97-116.
12. Zhang K.-K., Busse F.H. Generation of magnetic fields by convection in a rotating spherical fluid shell of infinite Prandtl number // Phys. Earth Planet. Inter. 1990. Vol.59. P.208-222.
116_
ISSN 0135-3500. Записки Горного института. Т.187
