162 research outputs found
Two dimensional code for modeling of high ione cyclotron harmonic fast wave heating and current drive
A two dimensional numerical code for computation of the electromagnetic field of a fast magnetosonic wave in a tokamak at high harmonics of the ion cyclotron frequency has been developed. The code computes the finite differ-ence solution of Maxwell’s equations for separate toroidal harmonics making use of the toroidal symmetry of tokamak plasmas. The proper boundary conditions are prescribed at the realistic tokamak vessel. The currents in the RF antenna are specified externally and then used in Ampere’s law. The main poloidal tokamak magnetic field and the "kinetic" part of the dielectric permeability tensor are treated iteratively. The code has been verified against known analytical solutions and first calculations of current drive in the spherical torus are presented.Разработан двумерный код для расчёта электромагнитного поля быстрой магнитозвуковой волны в токамаках в диапазоне высоких гармоник ионной циклотронной частоты. Используя тороидальную симметрию плазмы токамака, методом конечных разностей код находит решение уравнений Максвелла для отдельной тороидальной гармоники. Граничные условия задаются на реальной стенке камеры токамака. ВЧ-антенна моделируется введением «внешнего» тока в закон Ампера. Полоидальное удерживающее магнитное поле и вклад «тепловой» части тензора диэлектрической проницаемости плазмы включены в расчёты с помощью итераций. Код был верифицирован путём сравнения расчётов с известным аналитическим решением. Представлены первые расчёты тока увлечения в сферическом торе.Розроблено двовимірний код для розрахунку електромагнітного поля швидкої магнітозвукової хвилі в токамаках у діапазоні високих гармонік циклотронної частоти іонів. Використовуючи тороїдальну симетрію плазми токамаків, код знаходить рішення рівнянь Максвела методом скінчених різниць для окремої тороїдальної гармоніки. Граничні умови задаються на реальній стінці камери токамака. ВЧ-антена моделюється введенням «зовнішнього» струму в закон Ампера. Полоїдальне магнітне поле, що утримує, і внесок від «теплової» частини тензора діелектричної проникливості плазми включені в розрахунки за допомогою ітерацій. Код було перевірено шляхом порівняння розрахунків з відомим аналітичним рішенням. Наведено перші розрахунки струму захоплення в сферичному торі
Integration of the guiding-center equations in toroidal fields utilizing a local linearization approach
Computation of neoclassical toroidal viscosity with the account of non-standard orbits in a tokamak
Calculations of on effective ripple for a stellarator magnetic field computed by the HINT2 code
Improvement of collisionless particle confinement in a non-quasi-symmetric stellarator vacuum magnetic field
On the magnetic nature of electron transport barriers in tokamaks
The formation of internal transport barriers in the vicinity of rational magnetic surfaces in tokamaks with braided magnetic fields is studied for a simplified model of the perturbed magnetic field with a broad spatial spectrum and a monotonous shear profile. The island overlap criterion is used to derive a condition for barrier formation. This condition links the amplitude and the spectral width of the perturbation with the shear parameter. Numerical experiments with the MHD Monte-Carlo code E3D, where the problem of plasma heat conductivity is solved in 3D, confirm this formation of transport barriers in the case of a monotonous shear profile. Assuming that experimentally observed electron internal transport barriers are the result of local reduction of electron heat transport due to the magnetic field braiding, the amplitude and spectral width of magnetic perturbations are estimated for the tokamak RTP
Kinetic estimate of the shielding of resonant magnetic field perturbations by the plasma in DIII-D
Quasi-geometric integration of guiding-center orbits in piecewise linear toroidal fields
A numerical integration method for guiding-center orbits of charged particles
in toroidal fusion devices with three-dimensional field geometry is described.
Here, high order interpolation of electromagnetic fields in space is replaced
by a special linear interpolation, leading to locally linear Hamiltonian
equations of motion with piecewise constant coefficients. This approach reduces
computational effort and noise sensitivity while the conservation of total
energy, magnetic moment and phase space volume is retained. The underlying
formulation treats motion in piecewise linear fields exactly and thus preserves
the non-canonical symplectic form. The algorithm itself is only quasi-geometric
due to a series expansion in the orbit parameter. For practical purposes an
expansion to the fourth order retains geometric properties down to computer
accuracy in typical examples. When applied to collisionless guiding-center
orbits in an axisymmetric tokamak and a realistic three-dimensional stellarator
configuration, the method demonstrates stable long-term orbit dynamics
conserving invariants. In Monte Carlo evaluation of transport coefficients, the
computational efficiency of quasi-geometric integration is an order of
magnitude higher than with a standard fourth order Runge-Kutta integrator.Comment: 38 pages, 11 figure
- …