ПЛАНЕТНЫЕ КОЛЬЦА КАК РЕЛИКТЫ ПЛАЗМЕННЫХ ПРАКОЛЕЦ Борис И. Рабинович ИКИ РАН Семинар «Механика, Управление, Информатика», 27 января 2005 Институт Космических Исследований Российской Академии Наук

Аннотация Проводится анализ положения и структуры систем колец больших планет, основанный на решении задачи динамики и устойчивости вращающихся колец, состоявших из замагниченной плазмы. Последние рассматриваются при этом как реликты плазменных праколец. Приводятся численные оценки спектра собственных значений безразмерного параметра Ω, определяющего элитные пракольца. Теория прилагается к кольцам Юпитера, Сатурна и Урана и к плазменным торам Юпитера и Сатурна (J-тор и S-тор). 1

PLANETS RINGS AS RELICS OF ANCIENT PLASMA RINGS The analysis is fulfilled of great planets rings systems position and structure based on solution of dynamic and stability problems of the rotating magnetized plasma rings. Planets rings systems are considered as relics of ancient plasma rings. The numerical estimation of the non- dimensional parameters Ω eigen-values spectrum, determining the elite rings are presented. The theory has been applied to Jupiter, Saturn and Uranus rings and to the Jovian torus (Io- or J-torus) and Saturn torus (S-torus). 2

СОДЕРЖАНИЕ Вводные замечания Уравнения движения Магнито-гравитационные волны Магнито-гироскопические волны Приложение к планетным кольцам Приложение к плазменным торам Вместо заключения Литература 3

ВВОДНЫЕ ЗАМЕЧАНИЯ 4

Основное отличие Солнечной системы состоит, по-видимому, в том, что в прежние времена плотность плазмы, окружающей центральное тело, было гораздо выше… В результате могло реализоваться состояние частичной коротации свободно вращающейся плазмы. Строение пояса астероидов и колец Сатурна объясняется этим явлением …» Х. Альвен 5

«…Устойчивость, явление принципиально общее, как-то должна, по-видимому, проявляться в законах природы». «... Квантовани е устойчивых орбит динамики, подобно тому, как в квантовой механике квантуются орбиты электронов. Иначе говоря, устойчивы могут быть траектории не при любых начальных данных, а лишь при некоторых, исключительных». Н.Г. Четаев 6

Автор благодарен Анне Васильевне Калининой, расчеты которой подтвердили правильность главной идеи, Роберту Ефимовичу Ламперу за содействие в публикации предварительного сообщения в СибНИА им. Чаплыгина, Равилю Равильевичу Назирову за внимание и интерес к этой работе. Автор хотел бы особо отметить творческий вклад Виктории Прохоренко в подготовку этой презентации, который трудно переоценить. 7

I. Уравнения движения 8

Модель плазменного кольца d / b

Общие уравнения V/ t + 2 ( 0 V) + 1/ grad p – grad U = = ( 0 rot H H + P)/ ; H/ t + H = rot (V H); div V = 0; div H = 0. V– поле скоростей; H – напряженность магнитного поля; U – потенциал поля массовых сил невозмущенного движения (гравитационных и центробежных); – массовая плотность среды; P – векторное поле внешних сил, действующих в возмущенном движении; (1.2) 10

Уравнения невозмущенного движения A = 1 / Al = V G / a; V G = V G0 (R 0 / r 0 ) 1/2 ; A 2 = A 0 2 (r 0 / R 0 ) 5 ; a = a 0 (R 0 / r 0 ) 3 ; A 0 = 1 / Al 0 = V G0 / a 0 ; V G0 = (g 0 R 0 ) 1/2 ; a 0 = H 0 ( 0 / ) 1/2. V M 2 = V G 2 - 3a 2 ; V M = 0 r 0. (1.4) (1.3) 11

Уравнения возмущенного движения (1) V = V 0 + v; H = H 0 + h; v

Уравнения возмущенного движения (2) rot E / t + 2 ( 0 rot E) 0 / (H 0, ) rot F = 0; rot F / t + (V 0, ) rot F (H 0, ) rot E = -2 H 0 / t; rot r E L j = 0; rot r F L j = 0, (j = 1, 2). (1.7) (1.8) = 0; (p/ U ) = 0; 2 / t ( 0 / t) + (p/ U ) = 0; r (p/ ) L j = 0, (j = 1, 2). 13

II. Магнито-гравитационные волны 14

Свободные колебания (1) M [q mj ( mj, nj ) dS - 2p mj ( 0 x jm, nj ) dS + m=1 S * S * + ( (p U ) nj ) dS] = 0; n = 1, 2, …, M; j = 1,2. S * mj (r, ) = mj (r) cos m ; mj (r, ) = mj (r) sin m ; mj (r) = A mj r m + B mj r -m ; m = 1,2,…; j = 1,2. M = [p mj (t) mj + q mj (t) mj ]. m=1 (2.1) (2.2) 15

Свободные колебания (2) s nj – 2 i o nj s nj + 2 j s nj = 0; s nj = p nj + iq nj ; n = 1, 2,, M; j = 1, 2. (2.3) (2.4) 16

Вынужденные колебания (1) w = r L j = j 0 R 0 cos( - t); h = i x h 0 cos( - t); = - i x ; h 0 = H 0 0 (R 0 / r) 3. (2.5) j 0 = ( j 2 ( + ) - j 2 )/( j ) = = j0 [( j r j ) 2 +V G 2 (1+Al 2 )]/(- 2 + j 2 )r j 2 ; j0 = - = -3 Al 2 (1 – 3 Al 2 ). (2.6) j 2 = - 2 U r 2 L j = 3[(8a 2 – V G 2 )/r 2 ] L j ; j 2 = 2 U G r 2 L j = 2g 0 R 0 2 /r j 3 ; j 2 = - 2 U M r 2 L j = 21a 0 2 R 0 6 /r 8 j. (2.7) 17

Вынужденные колебания (2) 18

III. Магнито-гироскопические волны 19

Общие уравнения rot E / t + 2 ( 0 rot E) 0 / (H 0, ) rot F = 0; rot F / t + (V 0, ) rot F (H 0, ) rot E = -2 H 0 / t; rot r E L j = 0; rot r F L j = 0, (j = 1, 2). (3.1) v* = rot E; h* = rot F. (3.2) r * = r/r 0 ; b * = b/r o

Общее решение N v * = (t) rot i x C (r *, ) + (t) rot i x S (r *, ) ; =1 N h* = (t) rot i x C (r *, ) + (t) rot i x S (r *, ). =1 (3.4) C (r *, ) = sin (r * ) cos ; S (r *, ) = sin (r * ) sin ; (3.5) (r * ) = b * (r * – 1 + b * /2), ζ = + i ; = + i. (3.6) 21

Уравнения Б. Г. ' i d g e -i t + g e i t ( 0 / ) ' - ie - g e -i t - g +2 1 e i t +2 =0; = 1,3,5,, N - 1; g = 0 при = 1; g +2 = 0 при = N – 1. (3.7) ( 0 / ) ' - i e - g 1, e -i t -1 - g +1, e i t +1 = = - q e -i t ; +1 ' +1 i d g +1 e -i t + g +1 2 e i t 2 = 1,3,5,, N - 1; g -1 = 0 при = 1; g +1 2 = 0 при = N – 1; q = 0 при 1. (3.8) 22

Приведение к автономной системе = e i t ; * = e i t. (3.9) = (a 0 R 0 ); = ( 0 ) 1/2 /( a 0 R 0 ); = a 0 R 0 3 t / r 0 4 ; = r 0 2 / (a r 0 ); * 0 = 0. (3.10) = ; a 0 = H 0 ( 0 ) 1/2 ; a = a 0 R 0 3 / r 0 3. (3.11) 23

Уравнения автономной системы ' + i d * + g *, - + g * ' i( ) e* +1 - g * + - g * при +2 при. (3.12) ' i e* - g * -, -1 - g * q * ; ' +1 + i( ) d * ' +1 + g * +, + g * 2 2 = 0, -1 при +2 = 0 при ; q * = 0 при. (3.13) 24

Коэффициентная матрица G 1 при N = 8 id 1 * g 12 * g 12 * -2ie 2 * -g 32 * g 32 * 3id 3 * g 34 * g 34 * -4ie 4 * -g 54 * g 54 * 5id 5 * g 56 * g 56 * -6ie 6 * -g 76 * g 76 * 7id 7 * g 78 * g 78 * -8ie 8 * 25

Уравнения в области Лапласовых изображений. Асимптотическое приближение (p - i ) ( ) ( ) p ; p i + ( ) ( + 2)/2) 2 p. (3.14) (p i ) μ (( - 1) ) (( +1) ) p 0 - q* ; [p i ] ( ) 2 p 0 = 0 q* = 0 при. (3.15) g * = μ/2; g = μ/2 при π/b* или μ. 26

Вынужденные колебания i μ 0 (( - 1) ) 0 (( +1) ) 0 - q* ; i 0 ( ) = 0 0 N ; q* = 0 при. (3.16) -2iΩη 1 0 – ξ 2 0 = - q 1 Ω; -2iΩ ξ η 1 0 = 0. (3.17) 27

Элитные значения параметра (1) e -i t * e -i t a 0 R 0 a 0 R 0 ( 0 ) -1/2 ; (nm) const (nm) (nm) const (nm). m nm = 2 r nm. e i ; e i ; r a = r 4 a 0 R 0 3 ; nm - m nm = 0. (3.18) (3.19) (3.20) (3.21) 28

Периодическое движение на торе p =, q =. 29

Учет диссипативных сил ' i d + + g μ e -i t + g e i t ( 0 / ) ' - ie g e -i t - g +2 e i t +2 = 0. (3.7a) 30

Спектр элитных значений (1) 31

Спектр элитных значений (2) 32 m

Аналитическое выражение для спектра критических значений det G j = ; j = 1, 2; [(μ+2)/(μ+1)] +3 – 2(λ 2 +1) +1 + [μ/(μ+1)] -1 = 0; μ = 1, 3, …, N – 1; 0 = 0; N+2 = 0. z k ch α z k+1 + z k = 0; ch α = k = 0, 1, 2,…, K – 1; K = N/2; z 0 = 0; z K+1 = 0. s - sin ( s / ) s. -3 / < s < -. (3.25) (3.22) (3.23) (3.24) 33

Аналитическое выражение для спектра элитных значений s ( ) sin s ) s + 1, / < s < ; ; 2(s - 1) + для (3.14) m = 2s для (3.15) (3.27) (3.26) 34

Механическая модель s ( d ) 1 s ; s = 1, 2,,, (3.28) s (m ) -1/2 sin s 35 M = M 0 m j ; m = /(K+1); = 1, 2, …, K. T – сила натяжения струны; d – расстояние между массами M

Спектр элитных значений (3) 36

Целые квантовые числа и уравнение Шредингера kn = (n/ ) 1/2 sin ( k / 2n); k = 0, l; l = 1, 2, ; n = l + ; = 1, 2. (3.29) (3.30) В решении уравнения Шредингера, соответствующем атому водорода, числа n, l, k совпадают со следующими квантовыми числами: «Главным квантовым числом" n, «Квантовым числом полного кинетического момента" l, «Магнитным квантовым числом" k. n = K + 1; l = K ;k = 0, s. 37

Спиновые возмущения и полуцелые квантовые числа n j = (n/ ) 1/2 sin (j / 2n); j = 1/2, 3/2,, ; = l - 1/2; n = + ; (3.31) = + ½. Вихревое кольцо 38

ES-область 39 Al 0 > (1/3) 1/2 ; r 01 /R 0 < r 0 /R 0 < r 02 /R 0, r 01 /R 0 = 3 1/5 Al 0 2/5 ; r 02 /R 0 = 8 1/5 Al 0 2/5. = r 0 2 / (a r 0 ); max /. (3.32) (3.33) r 01 /R 0 r 02 /R 0

IV. Приложение к планетным кольцам 40

Системы колец Сатурна и Урана 41 Кольца Сатурна (Cassini) Кольца Урана (Voyager 2) Кольца Сатурна: Фрагмент кольца А (Cassini); Кольца B и C (Voyager 2)

Voyager и кольца Сатурна (1) Руководитель группы отображения миссий Voyager Брэд Смит: «Сотни неожиданных колечек внутри колец появляются перед электронными глазами ` Вояджера `. Деление Кассини, предположительно свободная зона между внешним кольцом А и средним кольцом В, ожила, по меньшей мере, тремя дюжинами малых колечек. Удивительные спицы протянулись в радиальном направлении поперек кольца В». «Большинство астрономов полагают, что с тех пор как кольца еще формировались, этот процесс всегда является рутинным следствием законов небесной механики, особенно - механизма, называемого резонансом. До Вояджера такие резонансы предполагались ответственными за те локальные структуры, которые демонстрируют кольца. Но теперь мониторы JPL (Jet Propulsion Laboratory) видят больше структур не только в кольцах, но и в делении Кассини, чем может объяснить любая симфония резонансов». 42

Voyager и кольца Сатурна (2) «Мы думали, что увидели все, что можно было увидеть, но в этом странном мире колец Сатурна причудливое стало банальностью; когда мы посмотрели сегодня на кольцо, вот, что мы увидели»: То, что показал Смит, - это картина кольца F, расщепленного на три пряди, две из которых оказались переплетенными, как косы. «'Косы' бросают вызов законам небесной механики по многим причинам, - сказал он, - но, очевидно, эти кольца ведут себя правильно (these rings are doing the right thing); я полагаю, что мы просто не очень хорошо понимаем законы…» 43

Cassini и кольца Сатурна Руководитель группы отображения миссий Cassini Каролин Порко: «Стало немедленно очевидно после изучения 61 изображения, что мы обнаружили целую коллекцию новых явлений, которые никогда не удавалось увидеть раньше, а некоторые из которых даже предвидеть. Мы видели много узких колечек, отделенных от широких лент, разрезанных на множество казавшихся некогерентными структур. Но по большей части мы видели волны, много волн…» «Волны являются точно идентифицируемыми, но странными структурами, и мы не имеем сколько-нибудь полного удовлетворительного объяснения того, почему они выглядят именно так. Соответствующие теории `спотыкались` много лет тому назад, когда отсутствовали достаточно точные данные, на которых они могли бы базироваться. Сейчас, вернувшись к Сатурну, мы имеем щедрую информацию, которая, конечно, будет стимулировать изучение формирования и распространения этих волн…» Гамлет: «И в небе и в земле сокрыто больше, чем снится вашей мудрости, Горацио» 44

Гипотеза о стратификации эволюционно зрелого плазменного пракольца Эволюционно зрелое вращающееся пракольцо, состоящее из замагниченной плазмы, оказывается расщепившееся на большое количество элитных колец, определяемых элитными собственными значениями фундаментального безразмерного параметра Ω. Дальнейшая стратификация пракольца и наличие дополнительных возмущений приводит к появлению уплотнений, становящихся центрами аккреции. В этих центрах начинают формироваться, благодаря известному механизму гравитационной неустойчивости, твердые частицы, что, в конечном счете, приводит к картине наблюдаемой сейчас. С этой точки зрения современные планетные кольца представляют собой реликты наиболее "жизнеспособных", эволюционно зрелых плазменных праколец. 45

Эволюция элитных колец 46

Схема плазменного пракольца 47

ES-область для колец Юпитера, Сатурна и Урана 48

Относительная плотность плазмы праколец Сатурна 49

V. Приложение к планетным торам 50

Фотографии J-тора 51

Магнитное поле Юпитера и схема J-тора X Y OGOG 52

Может ли существовать J-тор? Мнение компетентных специалистов «Этот тор является динамически неустойчивым (?!), поскольку направленная к периферии центробежная сила, обусловленная коротацией, значительно превосходит направленную к центру гравитационную силу Юпитера. Эта неустойчивость создает отток плазмы из тора, что должно, в конечном счете, приводить к ее ускользанию в магнитосферу планеты. И хотя такой отток является неизбежным, детальная природа его неизвестна и поэтому является предметом активного обсуждения. Результаты наблюдений, проведенных космическим аппаратом Galileo, возможно, внесут ясность в этот вопрос». Hill, T.W. and A.J. Dessler. Space Physics and Astronomy Converge in Exploration of Jupiter's Magnetosphere // EOS Trans. AGU P Скорее нет, чем да. Судите сами: 53

54 Может ли существовать J-тор? Мнение автора Скорее да, чем нет. Судите сами: Область существования и устойчивости J-тора и пракольца Юпитера (ES-область) Тор Пракольцо

Характерные скорости Абсолютные скорости, км/с 1V Io VJVJ aVMVM Относительные скорости, км/с 1 V Io - V M V J - V Io V J - V M (V J -V Io ) - V M

Измеренные и расчетные значения параметров J-тора (1) Размеры тора R 1 /R 0 ; R 2 /R 0 Измерение Расчет 1 Расчет 2 5.0; 6.2– Амплитуда коротационной волны 20 /R 0 Измерение Расчет 1 Расчет – 0.05–

Измеренные и расчетные значения параметров J-тора (2) Скорость относительноИо км/с Измерение Расчет 1 Расчет Cекундныйприход массы кГ/с Измерение Расчет 1 Расчет Произв. 57

Схема МГД- элемента 58

ES-область для S-тора 59

ВМЕСТО ЗАКЛЮЧЕНИЯ Вернемся в заключение к проблеме превращения совокупности плазменных элитных колец, на которые расщепилось эволюционно зрелое пракольцо, в систему колец, наблюдаемых у каждой из больших планет, состоящее из твердых элементов. Возникает принципиальный вопрос, заключающийся в следующем. Как получилось, что элитные кольца не только «запомнили» до мельчайших деталей свою конфигурацию в момент «рождения», но и сумели сохранить ее до наших дней? 60

Ответ на этот вопрос, по-видимому, дает теорема В.И. Арнольда о вековой устойчивости Солнечной системы. Каждое из элитных плазменных колец превратилось в процессе эволюции в соответствии с гипотезой, сформулированной выше, в пылевое кольцо, состоящее из твердых частиц. В дальнейшем частицы все больше укрупнялись благодаря гравитационной неустойчивости. Однако каждое элитное кольцо сохранило при этом свою индивидуальность и заняло определенное место в пространстве. В конечном счете, сформировалась система многочисленных узких колечек, состоящих из твердых элементов. Эти элементы должны неограниченно долго существовать на тех орбитах, на которых они «родились», за исключением тех орбит, которые были скорректированы за счет резонансных эффектов. 61

Литература 1.Рабинович Б.И. Магнитогидродинамика вращающихся вихревых колец из замагниченной плазмы // ДАН Т С Рабинович Б.И. Плазменное кольцо, вращающееся в гравитационно-магнитном поле. Вопросы устойчивости // ДАН Т С Рабинович Б.И., Прохоренко В.И. Космический аппарат с жидкостью, стабилизированный вращением, плазменный тор и задача Альвена // Полет С Рабинович Б.И, Прохоренко В.И. Задача Альвена и планетные пракольца. Проблема частичной коротации. Препринт ИКИ РАН Пр С Рабинович Б.И, Прохоренко В.И. Задача Альвена и планетные пракольца. Проблема квантования и устойчивости. Препринт ИКИ РАН Пр С