Тороидальная система координат

Материал из testwiki
Перейти к навигации Перейти к поиску

Тороидальная система координат — ортогональная система координат в пространстве, координатными поверхностями которой являются торы, сферы и полуплоскости. Данная система координат может быть получена посредством вращения двумерной биполярной системы координат вокруг оси, равноудалённой от фокусов биполярной системы.

Определение

Тороидальная система координат (α,β,φ) определяется посредством формул перехода из этих координат в декартовы координаты:

x=cshαcos⁡φchα−cos⁡βy=cshαsin⁡φchα−cos⁡βz=csin⁡βchα−cos⁡β,

где c>0 — масштабный множитель и радиус окружности x2+y2=c2,z=0 в которую вырождается тороидальная координатная поверхность α=α0=const при α0→∞. Пределы изменения координаты 0⩽α<∞. Обращаясь в бесконечность на указанной окружности, она стремится к нулю на бесконечности, а также в любой точке оси x=0,y=0. Две другие координаты являются циклическими с периодом 2π, например можно выбрать −π<β⩽π,−π<φ⩽π

Формулы перехода из тороидальных координат (α(z,r),β(z,r),φ) в цилиндрические координаты (z,r,φ):

r=cshαchα−cos⁡βz=csin⁡βchα−cos⁡β.

Для обратного преобразования при известных цилиндрических координатах точки (z,r,φ) вычисляют значения R±=r2+z2±2cr — максимальное и минимальное расстояние от данной точки до окружности r=c,z=0, через которые затем выражаются

κ2(α)=1−e−2α=4rcR+2,κ′=e−α=R−R+,sin⁡β=2zcR−R+

Альтернативное определение

В русскоязычной литературе тороидальными могут называться и более простые координаты (ρ,β,φ), такие, что:

z=ρsinβr=c+ρcosβ

(в англоязычной литературе такие координаты называют Шаблон:Lang-en, а не Шаблон:Lang-en). В этом случае циклические координаты β,φ называют полоидальным и тороидальным углами соответственно. В приложении к расчётам тородальных плазменных конфигураций, таких как токамак, помимо этих терминов ещё используется термин „магнитная ось“ для окружности r=c,z=0, на которой ρ=0. Вблизи магнитной оси координаты β для обеих систем приближенно совпадают, а координаты ρ≡R− и u связываются между собой соотношением: 2κ′(u)≈ρ/c. Могут также вводиться криволинейные потоковые координатыШаблон:Sfn, в которых координатными поверхностями являются топологически тороидальные магнитные поверхности (на которых давление плазмы постоянно, а нормальная компонента магнитного поля равна нулю. В этом случае являющаяся аналогом переменных α или ρ „потоковая“ координата служит только „меткой“ магнитной поверхности и её числовое значение несущественно.

Свойства

Координатные поверхности

α=const — торы

(x2+y2−ccthα)2+z2=(cshα)2,

β=const — сферы

(z−cctgβ)2+x2+y2=(csin⁡β)2,

φ=const — полуплоскости

xcos⁡φ=ysin⁡φ.

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

gij=(c2(chα−cos⁡β)2000c2(chα−cos⁡β)2000c2sh2α(chα−cos⁡β)2),gij=((chα−cos⁡β)2c2000(chα−cos⁡β)2c2000(chα−cos⁡β)2c2sh2α).

Он является диагональным, так как тороидальная система координат является ортогональной.

  • Квадрат линейного элемента:
ds2=c2(chα−cos⁡β)2(dα2+dβ2+sh2αdφ2).
  • Квадрат элемента площади:
dS2=c4(chα−cos⁡β)4((dαdβ)2+sh2α(dαdφ)2+sh2α(dβdφ)2).
  • Элемент объёма:
dV=c3shα(chα−cos⁡β)3dαdβdφ.
hα=hβ=cchα−cos⁡β,hφ=cshαchα−cos⁡β.
∂(x,y,z)∂(α,β,φ)=c3shα(chα−cos⁡β)3.
Γij1=(0−sin⁡βchα−cos⁡β0−sin⁡βchα−cos⁡βshαchα−cos⁡β000shα(chαcos⁡β−1)chα−cos⁡β),
Γij2=(sin⁡βchα−cos⁡β−shαchα−cos⁡β0−shαchα−cos⁡β0000sh2αsin⁡βchα−cos⁡β),
Γij3=(00−1(chα−cos⁡β)sh2α00−sin⁡βchα−cos⁡β−1(chα−cos⁡β)sh2α−sin⁡βchα−cos⁡β0).

Вид дифференциальных операторов в тороидальных координатах

  • Градиент скалярной функции в тороидальных координатах задается следующим выражением:
grad⁡U(α,β,φ)=chα−cos⁡βc(∂U∂αe→α+∂U∂βe→β+1shα∂U∂φe→φ).
 div⁡𝐅=(chα−cos⁡β)2c2shα(∂Fα∂α+∂Fβ∂β+shα∂Fφ∂φ)−chα−cos⁡βc2shα(Fαshα−Fβsin⁡β)
Δu=(chα−cos⁡β)3c2shα(∂∂α(shαchα−cos⁡β∂u∂α)+∂∂β(shαchα−cos⁡β∂u∂β)+1(chα−cos⁡β)shα∂2u∂φ2)

Дифференциальные уравнения в тороидальных координатах

Уравнение Лапласа в тороидальных координатах имеет вид:

(∂∂α(shαchα−cos⁡β∂u∂α)+∂∂β(shαchα−cos⁡β∂u∂β)+1(chα−cos⁡β)shα∂2u∂φ2)=0

Решение удобно искать в виде:

u=v2chα−2cos⁡β,

тогда уравнение для функции v:

vαα+vββ+vαcthα+14v+1sh2αvφφ=0.

После чего можно разделить переменные:

v=A(α)B(β)Φ(φ).

В результате получится система:

{A″+cthαA′+(14−kφ2sh2α−kβ2)A=0B″+kβ2B=0Φ″+kφ2Φ=0

В случае уравнения Гельмгольца в тороидальных координатах переменные не делятся.

Примечания

Шаблон:Примечания

Литература

Ссылки

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