Status of Continuum Edge Gyrokinetic Code Physics Development
Description
We are developing an edge gyro-kinetic continuum simulation code to study the boundary plasma over a region extending from inside the H-mode pedestal across the separatrix to the divertor plates. A 4-D (ψ, θ, ε, μ) version of this code is presently being implemented, en route to a full 5-D version. A set of gyrokinetic equations[1] are discretized on computational grid which incorporates X-point divertor geometry. The present implementation is a Method of Lines approach where the phase-space derivatives are discretized with finite differences and implicit backwards differencing formulas are used to advance the system in time. A fourth order upwinding algorithm is used for particle cross-field drifts, parallel streaming, and acceleration. Boundary conditions at conducting material surfaces are implemented on the plasma side of the sheath. The Poisson-like equation is solved using GMRES with multi-grid preconditioner from HYPRE. A nonlinear Fokker-Planck collision operator from STELLA[2] in (ν#parallel#,ν#perpendicular#) has been streamlined and integrated into the gyro-kinetic package using the same implicit Newton-Krylov solver and interpolating F and dF/dt|coll to/from (ε, μ) space. With our 4D code we compute the ion thermal flux, ion parallel velocity, self-consistent electric field, and geo-acoustic oscillations, which we compare with standard neoclassical theory for core plasma parameters; and we study the transition from collisional to collisionless end-loss. In the real X-point geometry, we find that the particles are trapped near outside midplane and in the X-point regions due to the magnetic configurations. The sizes of banana orbits are comparable to the pedestal width and/or the SOL width for energetic trapped particles. The effect of the real X-point geometry and edge plasma conditions on standard neoclassical theory will be evaluated, including a comparison of our 4D code with other kinetic neoclassical codes (such as NCLASS[3] and XGC[4]) and experiments
Availability note (English)
Available from OSTI as DE00878209; PURL: https://www.osti.gov/servlets/purl/878209-b4xqBQ/Additional details
Identifiers
Publishing Information
- Imprint Pagination
- 3 p.
- Report number
- UCRL-PROC--212647
Conference
- Title
- 10. IAEA Technical Meeting on H-Mode Physics and Transport Barriers
- Dates
- 28-30 Sep 2005
- Place
- St. Petersburg (Russian Federation)
INIS
- Country of Publication
- United States
- Country of Input or Organization
- United States
- INIS RN
- 37071036
- Subject category
- S70: PLASMA PHYSICS AND FUSION TECHNOLOGY;
- Resource subtype / Literary indicator
- Non-conventional Literature
- Descriptors DEI
- ACCELERATION; ALGORITHMS; BOUNDARY CONDITIONS; COMPUTERIZED SIMULATION; DIVERTORS; ELECTRIC FIELDS; FOKKER-PLANCK EQUATION; GEOMETRY; H-MODE PLASMA CONFINEMENT; N CODES; NEOCLASSICAL TRANSPORT THEORY; OSCILLATIONS; PHASE SPACE; PLASMA SCRAPE-OFF LAYER; PLASMA SIMULATION
- Descriptors DEC
- BOUNDARY LAYERS; CHARGED-PARTICLE TRANSPORT THEORY; COMPUTER CODES; CONFINEMENT; DIFFERENTIAL EQUATIONS; EQUATIONS; LAYERS; MAGNETIC CONFINEMENT; MATHEMATICAL LOGIC; MATHEMATICAL SPACE; MATHEMATICS; PARTIAL DIFFERENTIAL EQUATIONS; PLASMA CONFINEMENT; SIMULATION; SPACE; TRANSPORT THEORY
Optional Information
- Contract/Grant/Project number
- W--7405-ENG-48
- Notes
- PDF-FILE: 3 ; SIZE: 49.9 KBYTES
- Funding organization
- US Department of Energy (United States)