+7 499 196 95 39
Программа предназначена для расчета эффекта кривизны магнитного поля на устойчивость винтовых мод в токамаке, включая идеальные, резистивные моды и магнитные острова. В общем виде входными данными для расчета являются параметры равновесной конфигурации, включая задание 2D распределения потока полоидального магнитного поля и радиальных профилей давления и тока плазмы. Вклады в уравнение эволюции амплитуд вышеуказанных неустойчивостей рассчитываются путем интегрирования с соответствующими весами вдоль полного полоидального оборота силовой линии магнитного поля. Интегрирование выполняется по методу Рунге-Кутта 4-го порядка. Предусмотрена возможность аналитического задания магнитного поля для вычисления контурных интегралов. Для сравнения результатов, расчетов с имеющимися точными аналитическим решениями, а также для выявления эффектов геодезической кривизны предусмотрен предельный переход к цилиндрической геометрии. FORTRAN