Steepest‐descent moment method for three‐dimensional magnetohydrodynamic equilibria
暂无分享,去创建一个
An energy principle is used to obtain the solution of the magnetohydrodynamic (MHD) equilibrium equation J×B−∇p=0 for nested magnetic flux surfaces that are expressed in the inverse coordinate representation x=x(ρ, θ, ζ). Here, θ are ζ are poloidal and toroidal flux coordinate angles, respectively, and p=p(ρ) labels a magnetic surface. Ordinary differential equations in ρ are obtained for the Fourier amplitudes (moments) in the doubly periodic spectral decomposition of x. A steepest‐descent iteration is developed for efficiently solving these nonlinear, coupled moment equations. The existence of a positive‐definite energy functional guarantees the monotonic convergence of this iteration toward an equilibrium solution (in the absence of magnetic island formation). A renormalization parameter λ is introduced to ensure the rapid convergence of the Fourier series for x, while simultaneously satisfying the MHD requirement that magnetic field lines are straight in flux coordinates. A descent iteration is also d...
[1] P. Garabedian,et al. A computational method in plasma physics , 1978 .
[2] G. Hedstrom,et al. Numerical Solution of Partial Differential Equations , 1966 .
[3] Philip Rabinowitz,et al. Methods of Numerical Integration , 1985 .
[4] James M. Ortega,et al. Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.