FINITE ELEMENT METHODS TO CALCULATE IDEAL MAGNETOHYDRODYNAMIC FLOW EQUILIBRIA IN TOKAMAKS
The MHD equations which describe a stationary equilibrium in axisymmetric geometry lead to a quasilinear partial differential equation for the flux PSI, one nonlinear and three linear algebraic equations for the mass density rho, the poloidal current J and the two tangential velocity components nu-t and nu-p, respectively. These equations are solved numerically by a continuation method, by a Picard iteration and finite element approaches using conforming and non-conforming (hybrid) finite elements. The nodal points are redistributed iteratively such that they fall on initially prescribed constant flux surfaces. For different coordinate systems the different finite element approaches are applied to a static case and compared. The influence of flow is shown by applying our fixed foundary code CLIO to a particular plasma configuration (PDX Tokamak).
Bibliographic Reference: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, VOL. 52 (1985), PP. 675-682
Record Number: 1989124054500 / Last updated on: 1987-01-01
Available languages: en