Диполярная система координат

Диполя́рная, или дипо́льная, систе́ма координа́т — трёхмерная криволинейная ортогональная система координат, основанная на точечном (центральном) диполе, точнее, на его инвариантах преобразования координат.

Определение

Точное рассчитанное изображение силовых (красных) и эквипотенциальных (зеленых) линий диполя на плоскости, проходящей через его ось

Пропорционально показаны Земля (синяя) радиуса R с апексом дальней силовой линии 10R и границами термосферы и экзосферы

В диполярной системе координат, привязанной к точечному диполю, каждая точка пространства определяется тремя числами. При этом при фиксации одной из координат получается эквипотенциальная поверхность, а при фиксации двух других — силовая линия. Силовые линии перпендикулярны эквипотенциальным поверхностям. Диполярная система координат имеет вращательную (аксиальную) симметрию относительно оси диполя.

На рисунке справа (рассчитанном на компьютере) на плоскости, проходящей через ось диполя, показаны его силовые линии (красные), а также сечения эквипотенциальных поверхностей этой плоскостью (зеленые). Сам диполь находится в центре рисунка. Рисунок имеет две оси симметрии, горизонтальную и вертикальную, показанные прямыми линиями. Вертикальная прямая на рисунке является осью диполя. Силовые линии нарисованы красным цветом, они более вытянуты, расположены слева и справа от диполя, а линии зеленого цвета, более округлые, расположенные сверху и снизу диполя, суть сечения эквипотенциальных поверхностей («эквипотенциальные линии»). Координатные линии диполярной системы координат в трехмерном пространстве получаются вращением этого рисунка вокруг вертикальной оси.

Диполярная система координат широко используется при математическом моделировании дипольных систем. Причем обозначения координат, их порядок и направление не устоялись и могут меняться.

Определение диполярных координат через координаты других систем

Центры систем координат совпадают, также они соответственно ориентированы относительно друг друга: оси систем и долготы совпадают.

Сферическая система координат

Координатные составляющие $(\alpha ,\beta ,\lambda )$ диполярной системы, моделирующей магнитный диполь, определяются через сферические координаты $(r,\theta ,\phi )$ так:

$\alpha =\frac{r}{{\sin }^2\theta },\beta =\frac{\cos \theta }{r^2},\lambda =\phi .$

В соответствии с терминологией сферической системы координат здесь $r$ — расстояние до начала координат (радиальное расстояние), $\theta$ — зенитный, или полярный, угол, или наклонение, или коширота, $\phi$ — азимутальный угол. Уравнение $\beta =\text{const}$ определяет эквипотенциальную поверхность магнитного поля, а система уравнений $\alpha =\text{const},\lambda =\text{const}$ — силовую линию.

Значения диполярных координат имеют следующие ограничения:

$0⩽\alpha <+\mathrm{\infty }$, $-\mathrm{\infty }<\beta <+\mathrm{\infty }$, $0^{\circ }⩽\lambda <360^{\circ }$,

причем координаты $\beta$ и $\lambda$ (а также $\theta$ и $\phi$) не определены при $\alpha =0$, а координата $\lambda$ (и $\phi$) не определена еще и при $\theta =0^{\circ }$ и $\theta =180^{\circ }$.

Переход от компонент некоторого вектора $𝐀$ в сферических координатах к компонентам в диполярной системе осуществляется по формулам

$\left(\begin{array}{c}{A}_{\alpha }\\ A_{\beta }\\ A_{\lambda }\end{array}\right)=\frac{1}{\sqrt{\delta }}\left(\begin{array}{ccc}\sin \theta & -2\cos \theta & 0\\ -2\cos \theta & -\sin \theta & 0\\ 0& 0& \sqrt{\delta }\end{array}\right)\left(\begin{array}{c}{A}_r\\ A_{\theta }\\ A_{\phi }\end{array}\right),$

где

$\delta =1+3{\cos }^2\theta =1+3{\beta }^2{r}^4.$

Пусть ${𝐞}_{\alpha }$, ${𝐞}_{\beta }$, ${𝐞}_{\lambda }$ — координатные орты в этой диполярной системе координат. Тогда

${𝐞}_{\beta }\times {𝐞}_{\alpha }={𝐞}_{\lambda }$, ${𝐞}_{\alpha }\times {𝐞}_{\lambda }={𝐞}_{\beta }$, ${𝐞}_{\lambda }\times {𝐞}_{\beta }={𝐞}_{\alpha }$,

т. е. так определенная диполярная система координат является, по правилу буравчика, левой.

Выразить однозначно $(r,\theta ,\phi )$ через $(\alpha ,\beta ,\lambda )$ не удается, например, уравнения для определения $r,\theta$ такие:

$\alpha {\beta }^2{r}^4+r-\alpha =0,{\alpha }^2\beta {\sin }^4\theta -\cos \theta =0.$

Иногда используют безразмерное расстояние $\tilde{r}=\frac{r}{r_c}$, где $r_c$ — некоторое фиксированное расстояние, следующим образом:

$\tilde{\alpha }=\frac{{\sin }^2\theta }{\tilde{r}},\tilde{\beta }=\frac{\cos \theta }{{\tilde{r}}^2},\lambda =\phi .$

Тогда

${𝐞}_{\tilde{\alpha }}\times {𝐞}_{\tilde{\beta }}={𝐞}_{\lambda }$, ${𝐞}_{\tilde{\beta }}\times {𝐞}_{\lambda }={𝐞}_{\tilde{\alpha }}$, ${𝐞}_{\lambda }\times {𝐞}_{\tilde{\alpha }}={𝐞}_{\tilde{\beta }}$,

т. е. так определенная диполярная система координат является, по правилу буравчика, правой.

Декартова система координат

Координатные составляющие $(\alpha ,\beta ,\lambda )$ диполярной системы, моделирующей магнитный диполь, определяются через декартовы координаты $(x,y,z)$ и радиальное расстояние $r=\sqrt{x^2+y^2+z^2}$ так:

$\alpha =\frac{r^3}{x^2+y^2},\beta =\frac{z}{r^3},\lambda =\mathrm{arctg}\frac{y}{x}.$

Выразить однозначно $(x,y,z)$ через $(\alpha ,\beta ,\lambda )$ не удается:

$x=\frac{\sqrt{r^3}\cos \lambda }{\sqrt{\alpha }},y=\frac{\sqrt{r^3}\sin \lambda }{\sqrt{\alpha }},z=\beta r^3.$

Дифференциальные характеристики

Первые и вторые производные

Матрица Якоби перехода от декартовых к диполярным координатам имет вид:

$\frac{D(x,y,z)}{D(\alpha ,\beta ,\lambda )}=\frac{1}{\delta }\left(\begin{array}{ccc}{\sin }^3\theta (1-3{\cos }^2\theta )\cos \lambda & -3r^3\sin \theta \cos \theta \cos \lambda & -\delta r\sin \theta \sin \lambda \\ {\sin }^3\theta (1-3{\cos }^2\theta )\sin \lambda & -3r^3\sin \theta \cos \theta \sin \lambda & \delta r\sin \theta \cos \lambda \\ 3{\sin }^4\theta \cos \theta & r^3(1-3{\cos }^2\theta )& 0\end{array}\right).$

Матрица Якоби перехода от сферических к диполярным координатам имет вид:

$\frac{D(r,\theta ,\phi )}{D(\alpha ,\beta ,\lambda )}=\frac{1}{\delta }\left(\begin{array}{ccc}{\sin }^4\theta & -2r^3\cos \theta & 0\\ -(2{\sin }^3\theta \cos \theta )/r& -r^2\sin \theta & 0\\ 0& 0& \delta \end{array}\right).$

Коэффициенты Ламэ

$H_{\alpha }=\frac{{\sin }^3\theta }{\sqrt{\delta }},H_{\beta }=\frac{r^3}{\sqrt{\delta }},H_{\lambda }=r\sin \theta ,H_{\beta }={\alpha }^2{H}_{\alpha }{H}_{\lambda }.$

Вторые производные:

$\frac{{\partial }^{2}r}{\partial {\alpha }^2}=-\frac{4{\sin }^6\theta {\cos }^2\theta (5+3{\cos }^2\theta )}{r{\delta }^3},$

$\frac{{\partial }^{2}r}{\partial {\beta }^2}=\frac{2r^5(15{\cos }^4\theta +10{\cos }^2\theta -1)}{{\delta }^3},$

$\frac{{\partial }^2\theta }{\partial {\alpha }^2}=\frac{2{\sin }^5\theta \cos \theta (9{\cos }^4\theta +16{\cos }^2\theta -1)}{r^2{\delta }^3},$

$\frac{{\partial }^2\theta }{\partial {\beta }^2}=\frac{r^4\sin \theta \cos \theta (9{\cos }^2\theta +11)}{{\delta }^3}.$

Пусть $f$ — некоторая скалярная функция. Ее первые производные в диполярных и сферических координатах связаны:

$\delta \frac{\partial f}{\partial \alpha }={\sin }^4\theta \frac{\partial f}{\partial r}-2\frac{{\sin }^3\theta \cos \theta }{r}\frac{\partial f}{\partial \theta },$

$\delta \frac{\partial f}{\partial \beta }=-2r^3\cos \theta \frac{\partial f}{\partial r}-r^2\sin \theta \frac{\partial f}{\partial \theta },$

$\frac{\partial f}{\partial \lambda }=\frac{\partial f}{\partial \phi },$

или

$\frac{\partial f}{\partial r}=\frac{1}{{\sin }^2\theta }\frac{\partial f}{\partial \alpha }-\frac{2\cos \theta }{r^3}\frac{\partial f}{\partial \beta },$

$\frac{\partial f}{\partial \theta }=-2r\frac{\cos \theta }{{\sin }^3\theta }\frac{\partial f}{\partial \alpha }-\frac{\sin \theta }{r^2}\frac{\partial f}{\partial \beta },$

$\frac{\partial f}{\partial \phi }=\frac{\partial f}{\partial \lambda }.$

Ее оператор Лапласа равен

$\mathrm{\Delta }f=\frac{4}{r{\sin }^4\theta }\frac{\partial f}{\partial \alpha }+\frac{\delta }{{\sin }^6\theta }\frac{{\partial }^{2}f}{\partial {\alpha }^2}+\frac{\delta }{r^6}\frac{{\partial }^{2}f}{\partial {\beta }^2}+\frac{1}{r^2{\sin }^2\theta }\frac{{\partial }^{2}f}{\partial {\lambda }^2}.$

Векторные операции

Координаты векторных дифференциальных операторов в диполярной системе следующие:

$\mathrm{grad}f=\frac{\sqrt{\delta }}{{\sin }^3\theta }\frac{\partial f}{\partial \alpha }{𝐞}_{\alpha }+\frac{\sqrt{\delta }}{r^3}\frac{\partial f}{\partial \beta }{𝐞}_{\beta }+\frac{1}{r\sin \theta }\frac{\partial f}{\partial \lambda }{𝐞}_{\lambda },$

$\mathrm{div}𝐀=\frac{\sqrt{\delta }}{{\sin }^3\theta }\frac{\partial A_{\alpha }}{\partial \alpha }+4\frac{1-3{\cos }^4\theta }{{\sqrt{\delta }}^{3}r\sin \theta }{A}_{\alpha }+\frac{\sqrt{\delta }}{r^3}\frac{\partial A_{\beta }}{\partial \beta }-$

$-3\frac{\cos \theta (3+5{\cos }^2\theta )}{{\sqrt{\delta }}^{3}r}{A}_{\beta }+\frac{1}{r\sin \theta }\frac{\partial A_{\lambda }}{\partial \lambda },$

$\mathrm{rot}𝐀=(\frac{1}{r\sin \theta }\frac{\partial A_{\beta }}{\partial \lambda }-\frac{\sqrt{\delta }}{r^3}\frac{\partial A_{\lambda }}{\partial \beta }+\frac{3\cos \theta }{\sqrt{\delta }r}{A}_{\lambda }){𝐞}_{\alpha }+$

$+(\frac{\sqrt{\delta }}{{\sin }^3\theta }\frac{\partial A_{\lambda }}{\partial \alpha }-\frac{1}{r\sin \theta }\frac{\partial A_{\alpha }}{\partial \lambda }+\frac{1-3{\cos }^2\theta }{\sqrt{\delta }r\sin \theta }{A}_{\lambda }){𝐞}_{\beta }+$

$+(\frac{\sqrt{\delta }}{r^3}\frac{\partial A_{\alpha }}{\partial \beta }-\frac{\sqrt{\delta }}{{\sin }^3\theta }\frac{\partial A_{\beta }}{\partial \alpha }-6\frac{\cos \theta (1+{\cos }^2\theta )}{{\sqrt{\delta }}^{3}r}{A}_{\alpha }-$

$-3\frac{\sin \theta (1+{\cos }^2\theta )}{{\sqrt{\delta }}^{3}r}{A}_{\beta }){𝐞}_{\lambda }.$

Математическое моделирование Земли

Для описания поведения заряженных частиц в магнитном поле Земли наиболее удобна (гораздо удобнее, чем сферическая геомагнитная система координат) диполярная система координат.

Дипольная система координат $(\alpha ,\beta ,\lambda )$ при моделировании Земли строится следующим образом:

• её начало помещено в центр Земли;
• дипольная ось направлена по оси магнитного диполя Земли (геомагнитной оси), проходящей через магнитные полюса;
• координата $\alpha$ меняется перпендикулярно вектору магнитного поля Земли $𝐁$ в плоскости геомагнитного меридиана;
• координата $\beta$ меняется вдоль дипольной силовой линии, т. е. её направление совпадает с направлением вектора $𝐁$;
• координата $\lambda$ меняется в направлении, перпендикулярном первым двум.

Теоретически диполярная система координат может записываться и как левая система координат, когда координатный орт ${𝐞}_{\alpha }$ направлен от центра Земли, например, так:

$\alpha =\frac{r}{{\sin }^2\theta },\beta =\frac{\cos \theta }{r^2},\lambda =\phi ,$

и как правая система координат, когда координатный орт ${𝐞}_{\alpha }$ направлен к центру Земли, например, так:

$\tilde{\alpha }=\frac{{\sin }^2\theta }{\tilde{r}},\tilde{\beta }=\frac{\cos \theta }{{\tilde{r}}^2},\lambda =\phi ,$

где $\tilde{r}=\frac{r}{R_E}$, $R_E$ — радиус Земли.

В соответствии с терминологией сферической системы координат здесь $r$ — расстояние до начала координат (радиальное расстояние), $\theta$ — зенитный, или полярный, угол, или наклонение, или коширота, $\phi$ — азимутальный угол.

Координаты системы имеют следующий физический смысл:

• $\alpha$ и $L=\frac{1}{\tilde{\alpha }}$ — соответственно размерное (в м или км) и безразмерное (в радиусах Земли $R_E$) геоцентрическое расстояние до вершины силовой линии, $L$ — параметр Мак-Илвейна;
• $\beta$ и $\tilde{\beta }$ — соответственно размерный и безразмерный геомагнитный потенциал;
• $\lambda$ — геомагнитная долгота.

Примечания

Шаблон:Навигационная таблица