The structure of chaotic magnetic field lines in a tokamak withexternal nonsymmetric magnetic perturbations

We consider the effects of external nonsymmetric magnetostatic perturbations caused by resonant helical windings and a chaotic magnetic limiter on the plasma confined in a tokamak. The main purpose of both types of perturbation is to create a region in which field lines are chaotic in the Lagrangian sense: two initially nearby field lines diverge exponentially through many turns around the tokamak. The equilibrium field is obtained from the equations of magneto-hydrodynamic equilibrium written down in a polar toroidal coordinate system. The magnetic fields generated by the resonant helical windings and the chaotic magnetic limiter are obtained through an analytical solution of Laplace equation. The magnetic field line equations are integrated to give a Hamiltonian mapping of field lines that we use to characterize the structure of chaotic field lines. In the case of resonant windings, we obtained the map by both numerical integration and a Hamiltonian formulation. For a chaotic limiter, we analytically derived a symplectic map by using a Hamiltonian formulation     

托卡马克中漂移不稳定性的回旋动理学二维本征模方程

从弗拉索夫方程出发,导出了托卡马克等离子体中漂移不稳定性的回旋动理学二维本征模积分方程组。该方程组保留了离子的动理学效应,包括沿磁场的运动、磁场梯度和曲率漂移以及有限拉莫半径效应。与传统的采用气球模表象得到的一维回旋动理学方程(其只能给出不稳定模沿磁场线的结构)不同,该方程组不仅能给出托卡马克等离子体中漂移不稳定模的径向结构,同时还考虑了由离子的环形性漂移引起的相邻极向模之间的线性耦合,进而得到模的极向结构。该结果为相应的数值模拟研究提供了理论基础。     

A review on the solution of Grad–Shafranov equation in the cylindrical coordinates based on the Chebyshev collocation technique

Equilibrium reconstruction consists of identifying, from experimental measurements, a distribution of the plasma current density that satisfies the pressure balance constraint. Numerous methods exist to solve the Grad–Shafranov equation, describing the equilibrium of plasma confined by an axisymmetric magnetic field. In this paper, we have proposed a new numerical solution to the Grad–Shafranov equation (an axisymmetric, magnetic field transformed in cylindrical coordinates solved with the Chebyshev collocation method) when the source term (current density function) on the right-hand side is linear. The Chebyshev collocation method is a method for computing highly accurate numerical solutions of differential equations. We describe a circular cross-section of the tokamak and present numerical result of magnetic surfaces on the IR-T1 tokamak and then compare the results with an analytical solution.     

Internal magnetic field measurement in tokamak plasmas using a Zeeman polarimeter

In a tokamak plasma, the poloidal magnetic field profile closely depends on the current density profile. We can deduce the internal magnetic field from the analysis of circular polarization of the spectral lines emitted by the plasma. The theory of the measurement and a detailed design of the Zeeman polarimeter constructed to measure the poloidal field profile in the ADITYA tokamak are presented. The Fabry-Perot which we have employed in our design, with photodiode arrays followed by lock-in detection of the polarization signal, allows the measurement of the fractional circular polarization. In this system He-II line with wavelength 4686 Å is adopted as the monitoring spectral line. The line emission used in the present measurement is not well localized in the plasma, necessiating the use of a spatial inversion procedure to obtain the local values of the field.     

Spectral analysis of the precessional switching of the magnetization in an isotropic thin film

We study how a magnetic field step triggers the precessional switching of the magnetization in an isotropic thin film. Using a variational approach, we make an analytical estimate of the switching frequency. We compare it to more general analytical models, and to the results obtained numerically by direct integration of the equations of motion. We show that the periodic motions of the three magnetization components can be described satisfactorily with truncated Fourier expansions, indicating a relatively high spectral purity of the magnetic response. Our analytical expressions are simple enough to be physically transparent at first sight, in contrast to the results of the more elaborate models that treat also anisotropy.     

Onset of chaotic field line trajectories in tokamaks

Summary The onset of chaotic field line trajectories on tokamak plasmas perturbed by external resonant helical currents has been investigated. The spatial distribution of field lines is obtained integrating numerically the differential equations for the perturbed field lines and investigating the Poincaré maps characteristic of these fields. Several algorithms have been computed to identify the transition to large-scale stochasticity for this quasi-integrable dynamical system, which, for a given m.h.d. equilibrium, is determined by the enhancement of the coupling between the helical perturbation and the toroidicity, as the helical current increases. The authors of this paper have agreed to not receive the proofs for correction. Work partially supported by CNPq and FINEP.     

Plasma confinement in tokamaks with robust torus

We present a non-linear symplectic map that describes the alterations of the magnetic field lines inside the tokamak plasma due to the presence of a robust torus (RT) at the plasma edge. This RT prevents the magnetic field lines from reaching the tokamak wall and reduces, in its vicinity, the islands and invariant curve destruction due to resonant perturbations. The map describes the equilibrium magnetic field lines perturbed by resonances created by ergodic magnetic limiters (EMLs). We present the results obtained for twist and non-twist mappings derived for monotonic and non-monotonic plasma current density radial profiles, respectively. Our results indicate that the RT implementation would decrease the field line transport at the tokamak plasma edge.     

Magnetohydrodynamic equilibrium interpretation of low-aspect-ratio tokamak plasma

It is very important to calculate the equilibrium magnetic field configuration exactly and to estimate the parameters of the device needed for the equilibrium and stability in the design and operate the low-aspect-ratio tokamak, one of the axis-symmetrical torus devices which is widely used in the research of nuclear fusion using plasma. In the previous researches, there were the theoretical and numerical interpretation methods for the high-aspect-ratio tokamak plasma. But for the low-aspect-ratio tokamak plasma, numerical methods such as the finite element method and iteration method are usually used. In this research, we interpreted and examined the features of magnetohydrodynamic equilibrium of the low-aspect-ratio tokamak plasma in combination with the Green function method and finite element method which has high numerical value accuracy.     

Bloch方程的精细时程积分及其在射频脉冲设计中的应用

射频脉冲的频率选择性直接影响磁共振成像的质量，而射频脉冲的优化设计又归结为对Bloch方程的求解.尽管在某些情况下Bloch方程存在解析解，但由于其缺乏通用性而且形式上过于复杂而难于得到实用.本文提出一种Bloch方程的精细时程积分算法，并结合全局优化算法给出一个完整的射频脉冲设计方案.精细积分算法具有高效、高精度的特点，对于射频脉冲的设计很有裨益.数值算例表明，设计所得的射频脉冲具有较好的频率选择性. 关键词： 磁共振成像 射频脉冲 Bloch方程 精细时程积分     

Mean Field Analysis of Large-Scale Interacting Populations of Stochastic Conductance-Based Spiking Neurons Using the Klimontovich Method

We investigate the dynamics of large-scale interacting neural populations, composed of conductance based, spiking model neurons with modifiable synaptic connection strengths, which are possibly also subjected to external noisy currents. The network dynamics is controlled by a set of neural population probability distributions (PPD) which are constructed along the same lines as in the Klimontovich approach to the kinetic theory of plasmas. An exact non-closed, nonlinear, system of integro-partial differential equations is derived for the PPDs. As is customary, a closing procedure leads to a mean field limit. The equations we have obtained are of the same type as those which have been recently derived using rigorous techniques of probability theory. The numerical solutions of these so called McKean–Vlasov–Fokker–Planck equations, which are only valid in the limit of infinite size networks, actually shows that the statistical measures as obtained from PPDs are in good agreement with those obtained through direct integration of the stochastic dynamical system for large but finite size networks. Although numerical solutions have been obtained for networks of Fitzhugh–Nagumo model neurons, which are often used to approximate Hodgkin–Huxley model neurons, the theory can be readily applied to networks of general conductance-based model neurons of arbitrary dimension.     

Magnetic Field in a Screw Flow with Fluctuations

We consider the influence of fluctuations in a screw flow of a conducting liquid on the effect of magnetic field self-excitation; the solution of this problem is important for experimental realization of a turbulent dynamo. We propose a theoretical approach based on the solution of averaged equations obtained in the limit of a short correlation time. The applicability of this approach is confirmed by direct numerical simulation of the initial equations. We demonstrate the influence of the correlation of fluctuations on the dynamo effect threshold. It is shown that the solution of the mean-field equations differs from the solution based on direct numerical simulation for a finite correlation time. The advantages and disadvantages of the two approaches are estimates, as well as the importance of the discovered difference in the context of problems of magnetic field self-excitation. The influence of helicity and intermittency on the type of the solution is considered.     

Locating Cantori for Symmetric Tokamap and Symmetric Ergodic Magnetic Limiter Map Using Mean-Energy Error Criterion

We use a method based on the conservation of energy, the mean-energy error criterion, to approximately locate the place of a cantorus by locating the series of its convergents. The mean-energy error curve has nearly stationary parts in the vicinity of elliptic (minimax) orbits, the so-called magnetic islands. Stable minimax orbits converge to orbits homoclinic to a cantorus. By tracing the island series, we limit the cantorus to a narrow region. A near-critical perturbation parameter is used so that, while the cantorus may be destabilized, its high-order minimax orbits remain intact. As illustrations, we consider two symplectic maps, systematically derived from the Hamilton–Jacobi equation and Jacobi’s theorem, in the context of the magnetically confined plasmas in a tokamak: a symmetric tokamap realistically reproduces the main features of a tokamak, and a symmetric ergodic magnetic limiter (EML) map is defined to describe the action of EML rings on the magnetic field lines in the tokamak.     

CT-6B托卡马克的锯齿模拟

一、基本方程与计算方法 1.锯齿模式锯齿上升部分可做如下解释,焦耳加热使电流分布变尖,q减小,致使m=1的撕裂模变得不稳定,磁岛将增长。r_s是q=1处的半径,由下式给出 q(r_s)=1. (1) 磁岛宽度由给出。式中W_o为初值,Y是m=1/n=1模的线性增长率     

Tokamak plasma equilibrium analysis based on the relaxation method with a specified magnetic axis

In this contribution, we have presented two techniques for the determination of plasma equilibrium position in IR-T1 tokamak: relaxation and optical methods. An analysis method of tokamak plasma equilibrium by a relaxation method with a specified magnetic axis is presented. The degrees of freedom due to designated positions of the magnetic axis are possible by using poloidal field coil currents. Stable steady-state tokamak plasma equilibria are calculated along with the magnetohydrodynamic potential energy. The plasma generates a plasma current which partially or fully cancels the magnetic field from the poloidal field coils. For low-temperature plasmas, the plasma current distribution is centrally peaked; for high-temperature plasmas, the plasma current has a hole. A centrally peaked current distribution in a low-temperature plasma is evolved into a current distribution with a hole by increasing the plasma pressure by Ohmic heating, radio frequency heating, or by neutral beam injection heating. In the second technique, an image-processing technique was used for the output signal of the charge coupled device camera and plasma emission intensity profile and then the plasma position was obtained. Results are compared and discussed.     

Non-perturbative renormalization group calculation of the scalar self-energy

We present the first numerical application of a method that we have recently proposed to solve the Non Perturbative Renormalization Group equations and obtain the n-point functions for arbitrary external momenta. This method leads to flow equations for the n-point functions which are also differential equations with respect to a constant background field. This makes them, a priori, difficult to solve. However, we demonstrate in this paper that, within a simple approximation which turns out to be quite accurate, the solution of these flow equations is not more complicated than that of the flow equations obtained in the derivative expansion. Thus, with a numerical effort comparable to that involved in the derivative expansion, we can get the full momentum dependence of the n-point functions. The method is applied, in its leading order, to the calculation of the self-energy in a 3-dimensional scalar field theory, at criticality. Accurate results are obtained over the entire range of momenta.     

Simplified method for describing plasma temperature evolution in a tokamak with a growing magnetic island

A new simplified method is proposed for solving the heat conduction equation in a tokamak plasma with an unsteady magnetic island with allowance for heat transfer along the magnetic field. The method, whose efficiency is illustrated by a simple example, allows one to reduce the three-dimensional problem of plasma temperature evolution in the presence of several magnetic islands to a set of one-dimensional equations.     

修正的频域有限差分法在二维金属光子晶体分析中的应用

与介质光子晶体相比,金属光子晶体的带隙特性在毫米波和亚毫米波波段有着重要的应用价值.基于Yee网格的频域有限差分法推导得出的本征模方程,求解后能方便而又可靠地得出介质光子晶体的带隙图和场分布.但由于金属与介质的本质差异,该方法不能直接应用于金属光子晶体.文中引入了金属表面边界条件,推导了二维金属周期结构的光子带隙本征模方程.通过数值计算,得出了不同晶格结构（正方/三角格子）下两种模式（TE/TM）的全禁带特性,并与介质周期结构的禁带特性进行对比,分析了金属周期结构在模式选择和器件集成方面的优点.     

修正的频域有限差分法在二维金属光子晶体分析中的应用

与介质光子晶体相比,金属光子晶体的带隙特性在毫米波和亚毫米波波段有着重要的应用价值.基于Yee网格的频域有限差分法推导得出的本征模方程,求解后能方便而又可靠地得出介质光子晶体的带隙图和场分布.但由于金属与介质的本质差异,该方法不能直接应用于金属光子晶体.文中引入了金属表面边界条件,推导了二维金属周期结构的光子带隙本征模方程.通过数值计算,得出了不同晶格结构（正方/三角格子）下两种模式（TE/TM）的全禁带特性,并与介质周期结构的禁带特性进行对比,分析了金属周期结构在模式选择和器件集成方面的优点. 关键词： 金属光子晶体 频域有限差分法 全禁带     

Analytical magnetic distributions and Stoner–Wohlfarth theory for near-uniformly magnetised rectangular nanomagnets

The micromagnetic equations which govern the magnetisation distribution have been studied analytically for nano-sized magnets. By exploiting the fact that in such magnets the magnetisation is near-uniform these equations can be linearised. The problem then reduces to the solution of two uncoupled Poisson equations with given source terms subject to certain surface boundary conditions. These source terms are particularly simple and have analytical forms when one restricts the analysis to rectangular nanomagnets. Solutions to the Poisson equations have been found analytically by expansion in Fourier series which automatically satisfies the boundary conditions. The total energies obtained using these Fourier solutions are compared with the energies obtained from numerical solutions to the full non-linear micromagnetic equations in the case of square cross-section nanomagnets (the only case that has been considered numerically) and shown to be in agreement to the order of a few percent. These solutions may now be used to calculate the total energy in the presence of an external magnetic field and in this way extend the Stoner–Wohlfarth theory to near-uniformly magnetised rectangular bodies.