В физике теория динамо предлагает механизм, посредством которого небесное тело, такое как Земля или звезда, генерирует магнитное поле . Теория динамо описывает процесс, посредством которого вращающаяся, конвективная и электропроводящая жидкость может поддерживать магнитное поле в течение астрономических временных масштабов. Считается, что динамо является источником магнитного поля Земли и магнитных полей Меркурия и планет-гигантов .
Когда Уильям Гилберт опубликовал de Magnete в 1600 году, он пришел к выводу, что Земля является магнитной, и предложил первую гипотезу происхождения этого магнетизма: постоянный магнетизм, такой как тот, что обнаружен в магните . В 1822 году Андре-Мари Ампер предположил, что внутренние токи ответственны за магнетизм Земли. [2] В 1919 году Джозеф Лармор предположил, что динамо может генерировать поле. [3] [4] Однако даже после того, как он выдвинул свою гипотезу, некоторые выдающиеся ученые выдвинули альтернативные объяснения. Лауреат Нобелевской премии Патрик Блэкетт провел серию экспериментов в поисках фундаментальной связи между угловым моментом и магнитным моментом , но не нашел ни одной. [5] [6]
Уолтер М. Эльзассер , считающийся «отцом» ныне принятой теории динамо как объяснения магнетизма Земли, предположил, что это магнитное поле является результатом электрических токов, индуцированных в жидком внешнем ядре Земли. Он раскрыл историю магнитного поля Земли, став пионером в изучении магнитной ориентации минералов в горных породах.
Для того, чтобы поддерживать магнитное поле против омического распада (который произошел бы для дипольного поля через 20 000 лет), внешнее ядро должно быть конвективным. Конвекция, вероятно , является некоторой комбинацией тепловой и композиционной конвекции. Мантия контролирует скорость, с которой тепло извлекается из ядра. Источники тепла включают гравитационную энергию, выделяемую при сжатии ядра, гравитационную энергию, выделяемую при отторжении легких элементов (вероятно, серы , кислорода или кремния ) на внутренней границе ядра по мере его роста, скрытую теплоту кристаллизации на внутренней границе ядра и радиоактивность калия , урана и тория . [7]
На заре 21-го века численное моделирование магнитного поля Земли не было успешно продемонстрировано. Первоначальные модели были сосредоточены на генерации поля путем конвекции в жидком внешнем ядре планеты. Было возможно показать генерацию сильного поля, подобного земному, когда модель предполагала однородную температуру ядра-поверхности и исключительно высокую вязкость для жидкости ядра. Вычисления, которые включали более реалистичные значения параметров, давали магнитные поля, которые были менее земными, но указывали, что уточнения модели [ какие? ] могут в конечном итоге привести к точной аналитической модели. Небольшие изменения температуры ядра-поверхности в диапазоне нескольких милликельвинов приводят к значительному увеличению конвективного потока и создают более реалистичные магнитные поля. [8] [9]
Теория динамо описывает процесс, посредством которого вращающаяся, конвективная и электропроводящая жидкость действует для поддержания магнитного поля. Эта теория используется для объяснения наличия аномально долгоживущих магнитных полей в астрофизических телах. Проводящая жидкость в геодинамо — это жидкое железо во внешнем ядре, а в солнечном динамо — ионизированный газ в тахоклине . Теория динамо астрофизических тел использует уравнения магнитогидродинамики для исследования того, как жидкость может непрерывно восстанавливать магнитное поле. [10]
Когда-то считалось, что диполь , который составляет большую часть магнитного поля Земли и смещен вдоль оси вращения на 11,3 градуса, был вызван постоянным намагничиванием материалов в Земле. Это означает, что теория динамо изначально использовалась для объяснения магнитного поля Солнца в его связи с полем Земли. Однако эта гипотеза, первоначально предложенная Джозефом Лармором в 1919 году, была изменена из-за обширных исследований магнитных вековых вариаций , палеомагнетизма (включая инверсию полярности ), сейсмологии и обилия элементов в Солнечной системе. Кроме того, применение теорий Карла Фридриха Гаусса к магнитным наблюдениям показало, что магнитное поле Земли имеет внутреннее, а не внешнее происхождение.
Для работы динамо-машины необходимы три условия:
В случае Земли магнитное поле индуцируется и постоянно поддерживается конвекцией жидкого железа во внешнем ядре. Требованием для индукции поля является вращающаяся жидкость. Вращение во внешнем ядре обеспечивается эффектом Кориолиса , вызванным вращением Земли. Сила Кориолиса стремится организовать движения жидкости и электрические токи в столбы (см. также столбцы Тейлора ), выровненные с осью вращения. Индукция или генерация магнитного поля описывается уравнением индукции : где u — скорость, B — магнитное поле, t — время, а — магнитная диффузия с электропроводностью и проницаемостью . Отношение второго члена в правой части к первому члену дает магнитное число Рейнольдса , безразмерное отношение адвекции магнитного поля к диффузии.
Приливные силы между небесными телами, вращающимися по орбите, вызывают трение, которое нагревает их внутреннюю часть. Это известно как приливной нагрев, и это помогает поддерживать внутреннюю часть в жидком состоянии. Жидкая внутренняя часть, которая может проводить электричество, необходима для создания динамо. Энцелад Сатурна и Ио Юпитера имеют достаточно приливного нагрева, чтобы разжижить свои внутренние ядра, но они не могут создать динамо, потому что не могут проводить электричество. [12] [13] Меркурий, несмотря на свой небольшой размер, имеет магнитное поле, потому что у него есть проводящее жидкое ядро, созданное его железным составом и трением, возникающим в результате его высокоэллиптической орбиты. [14] Предполагается, что у Луны когда-то было магнитное поле, основанное на доказательствах из намагниченных лунных пород, из-за ее кратковременного близкого расстояния к Земле, создающего приливной нагрев. [15] Орбита и вращение планеты помогают обеспечить жидкое ядро и дополняют кинетическую энергию, которая поддерживает действие динамо.
В кинематической теории динамо поле скорости предписано , а не является динамической переменной: Модель не предусматривает искажения потока в ответ на магнитное поле. Этот метод не может обеспечить переменное во времени поведение полностью нелинейного хаотического динамо, но может быть использован для изучения того, как напряженность магнитного поля изменяется в зависимости от структуры потока и скорости.
Используя уравнения Максвелла одновременно с ротором закона Ома , можно вывести то, что по сути является линейным уравнением собственных значений для магнитных полей ( B ), что можно сделать, предположив, что магнитное поле не зависит от поля скорости. Получается критическое магнитное число Рейнольдса , выше которого сила потока достаточна для усиления наложенного магнитного поля, а ниже которого магнитное поле рассеивается.
Наиболее функциональная особенность теории кинематического динамо заключается в том, что ее можно использовать для проверки того, способно ли поле скорости к действию динамо. Экспериментально прикладывая определенное поле скорости к небольшому магнитному полю, можно наблюдать, имеет ли магнитное поле тенденцию к росту (или нет) в ответ на приложенный поток. Если магнитное поле растет, то система либо способна к действию динамо, либо является динамо, но если магнитное поле не растет, то ее просто называют «не динамо».
Аналогичный метод, называемый мембранной парадигмой, представляет собой способ изучения черных дыр , позволяющий описать вещество вблизи их поверхности на языке теории динамо.
Кинематическое динамо можно также рассматривать как явление спонтанного нарушения топологической суперсимметрии соответствующего стохастического дифференциального уравнения, связанного с потоком фоновой материи. [16] В рамках стохастической суперсимметричной теории эта суперсимметрия является внутренним свойством всех стохастических дифференциальных уравнений , ее интерпретация заключается в том, что фазовое пространство модели сохраняет непрерывность посредством непрерывных потоков времени. Когда непрерывность этого потока спонтанно нарушается, система оказывается в стохастическом состоянии детерминированного хаоса . [17] Другими словами, кинематическое динамо возникает из-за хаотического потока в базовой фоновой материи.
Кинематическое приближение становится недействительным, когда магнитное поле становится достаточно сильным, чтобы повлиять на движения жидкости. В этом случае поле скорости начинает подвергаться влиянию силы Лоренца , и поэтому уравнение индукции больше не является линейным в магнитном поле. В большинстве случаев это приводит к гашению амплитуды динамо. Такие динамо иногда также называют гидромагнитными динамо . [18] Практически все динамо в астрофизике и геофизике являются гидромагнитными динамо.
Основная идея теории заключается в том, что любое малое магнитное поле, существующее во внешнем ядре, создает токи в движущейся жидкости из-за силы Лоренца. Эти токи создают дополнительное магнитное поле из-за закона Ампера . При движении жидкости токи переносятся таким образом, что магнитное поле становится сильнее (пока отрицательно [19] ). Таким образом, «зародышевое» магнитное поле может становиться все сильнее и сильнее, пока не достигнет некоторого значения, связанного с существующими немагнитными силами.
Для моделирования полностью нелинейных динамо используются числовые модели. Используются следующие уравнения:
Эти уравнения затем обезразмериваются, вводя безразмерные параметры, где R a - число Рэлея , E - число Экмана , P r и P m - число Прандтля и магнитное число Прандтля . Масштабирование магнитного поля часто осуществляется в единицах числа Эльзассера
Скалярное произведение приведенной выше формы уравнения Навье-Стокса с дает скорость увеличения плотности кинетической энергии, , в левой части. Последний член в правой части тогда , локальный вклад в кинетическую энергию из-за силы Лоренца .
Скалярное произведение уравнения индукции с дает скорость увеличения плотности магнитной энергии, , в левой части. Последний член в правой части равен Тогда Поскольку уравнение является интегрируемым по объему, этот член эквивалентен с точностью до граничного члена (и с двойным использованием тождества скалярного тройного произведения ) (где использовалось одно из уравнений Максвелла). Это локальный вклад в магнитную энергию из-за движения жидкости.
Таким образом, этот термин — скорость преобразования кинетической энергии в магнитную энергию. Она должна быть неотрицательной хотя бы в части объема, чтобы динамо создавало магнитное поле. [19]
Из приведенной выше диаграммы неясно, почему этот член должен быть положительным. Простой аргумент может быть основан на рассмотрении чистых эффектов. Чтобы создать магнитное поле, чистый электрический ток должен обернуться вокруг оси вращения планеты. В этом случае, чтобы член был положительным, чистый поток проводящей материи должен быть направлен к оси вращения. Диаграмма показывает только чистый поток от полюсов к экватору. Однако сохранение массы требует дополнительного потока от экватора к полюсам. Если бы этот поток был вдоль оси вращения, это означало бы, что циркуляция была бы завершена потоком от тех, которые показаны, к оси вращения, создавая желаемый эффект.
Приведенная выше формула для скорости преобразования кинетической энергии в магнитную энергию эквивалентна скорости работы, совершаемой силой на внешнем веществе ядра, скорость которой . Эта работа является результатом немагнитных сил, действующих на жидкость.
Из них гравитационная сила и центробежная сила являются консервативными и, следовательно, не оказывают общего влияния на движение жидкости в замкнутых контурах. Число Экмана (определенное выше), которое является отношением между двумя оставшимися силами, а именно вязкостью и силой Кориолиса, очень мало внутри внешнего ядра Земли, поскольку его вязкость мала (1,2–1,5 ×10−2 паскаль -секунда [20] ) из-за его текучести.
Таким образом, основной усредненный по времени вклад в работу вносит сила Кориолиса, величина которой, хотя и связана с этой величиной, связана лишь косвенно и в общем случае не является одинаковой локально (то есть они влияют друг на друга, но не в одном и том же месте и времени).
Плотность тока J сама по себе является результатом магнитного поля согласно закону Ома . Опять же, из-за движения материи и тока это не обязательно поле в том же месте и в то же время. Однако эти соотношения все еще могут быть использованы для выведения порядков величин рассматриваемых величин.
С точки зрения порядка величины и , давая или:
Точное соотношение между обеими сторонами равно квадратному корню из числа Эльзассера .
Обратите внимание, что направление магнитного поля не может быть выведено из этого приближения (по крайней мере, его знак), поскольку оно представляется квадратичным и, действительно, иногда перевернутым , хотя в целом оно лежит на той же оси, что и .
Для внешнего ядра Земли ρ приблизительно равно 10 4 кг/м 3 , [20] Ω = 2 π /день = 7,3×10 −5 /секунду, а σ приблизительно равно 10 7 Ω −1 м −1 . [21] Это дает 2,7×10 −4 Тесла .
Магнитное поле магнитного диполя имеет обратную кубическую зависимость от расстояния, поэтому его порядок величины на поверхности Земли можно приблизительно оценить, умножив приведенный выше результат на ( R внешнего ядра ⁄ R Земли ) 3 = ( 2890 ⁄ 6370 ) 3 = 0,093, что дает 2,5×10−5 Тесла , что недалеко от измеренного значения 3×10−5 Тесла на экваторе .
В целом, модели геодинамо пытаются создавать магнитные поля, согласующиеся с наблюдаемыми данными при определенных условиях и уравнениях, как упоминалось в разделах выше. Успешная реализация уравнений магнитогидродинамики имела особое значение, поскольку они подтолкнули модели динамо к самосогласованности. Хотя модели геодинамо особенно распространены, модели динамо не обязательно ограничиваются геодинамо; солнечные и общие модели динамо также представляют интерес. Изучение моделей динамо полезно в области геофизики, поскольку это позволяет определить, как различные механизмы формируют магнитные поля, подобные тем, которые создаются астрофизическими телами, такими как Земля, и как они заставляют магнитные поля проявлять определенные особенности, такие как смена полюсов.
Уравнения, используемые в числовых моделях динамо, очень сложны. Десятилетиями теоретики были ограничены двумерными кинематическими моделями динамо , описанными выше, в которых движение жидкости выбирается заранее, а влияние на магнитное поле рассчитывается. Прогресс от линейных к нелинейным, трехмерным моделям динамо был в значительной степени затруднен поиском решений магнитогидродинамических уравнений, которые устраняют необходимость во многих предположениях, сделанных в кинематических моделях, и допускают самосогласованность.
Первые самосогласованные модели динамо, определяющие как движение жидкости, так и магнитное поле, были разработаны двумя группами в 1995 году: одной в Японии [22] и одной в Соединенных Штатах. [23] [24] Последняя была создана как модель относительно геодинамо и привлекла значительное внимание, поскольку она успешно воспроизводила некоторые характеристики поля Земли. [19] После этого прорыва произошел большой всплеск в разработке разумных трехмерных моделей динамо. [19]
Хотя сейчас существует много самосогласованных моделей, между моделями существуют значительные различия, как в результатах, которые они производят, так и в способе их разработки. [19] Учитывая сложность разработки модели геодинамо, существует много мест, где могут возникнуть расхождения, например, при принятии предположений, касающихся механизмов, которые обеспечивают энергию для динамо, при выборе значений для параметров, используемых в уравнениях, или при нормализации уравнений. Несмотря на множество различий, которые могут возникнуть, большинство моделей имеют общие черты, такие как четкие осевые диполи. Во многих из этих моделей также были успешно воссозданы такие явления, как вековые вариации и геомагнитные инверсии полярности . [19]
Многие наблюдения могут быть сделаны с помощью моделей динамо. Модели могут быть использованы для оценки того, как магнитные поля меняются со временем, и могут быть сравнены с наблюдаемыми палеомагнитными данными, чтобы найти сходства между моделью и Землей. Однако из-за неопределенности палеомагнитных наблюдений сравнения могут быть не совсем достоверными или полезными. [19] Упрощенные модели геодинамо показали взаимосвязи между числом динамо (определяемым дисперсией скоростей вращения во внешнем ядре и зеркально-асимметричной конвекцией (например, когда конвекция благоприятствует одному направлению на севере, а другому на юге)) и инверсиями магнитных полюсов, а также обнаружили сходство между геодинамо и динамо Солнца. [19] Во многих моделях кажется, что магнитные поля имеют несколько случайные величины, которые следуют нормальной тенденции, которая в среднем равна нулю. [ 19] В дополнение к этим наблюдениям, общие наблюдения о механизмах, приводящих в действие геодинамо, могут быть сделаны на основе того, насколько точно модель отражает фактические данные, собранные с Земли.
Сложность моделирования динамо настолько велика, что модели геодинамо ограничены нынешней мощностью суперкомпьютеров , в частности потому, что расчет чисел Экмана и Рэлея внешнего ядра чрезвычайно сложен и требует огромного количества вычислений.
Было предложено много улучшений в моделировании динамо с момента самосогласованного прорыва в 1995 году. Одним из предложений по изучению сложных изменений магнитного поля является применение спектральных методов для упрощения вычислений. [25] В конечном счете, пока не будут достигнуты значительные улучшения в вычислительной мощности, методы вычисления реалистичных моделей динамо должны быть сделаны более эффективными, поэтому внесение улучшений в методы вычисления модели имеет большое значение для продвижения численного моделирования динамо.