GTC–GyrokineticToroidalCodeTrainingcoursePartIIhor.pptVIP

GTC – Gyrokinetic Toroidal Code Training course: Part I Ihor Holod University of California, Irvine;Outline Introduction Gyrokinetic approach GK equations of motion. GK Poisson’s equation. Delta-f method Fluid-kinetic hybrid electron model: continuity equation, Ampere’s law, drift-kinetic equation Tokamak geometry. Magnetic coordinates Particle-in-cell approach. General structure of GTC. Grid. Parallel computing. Setup and loading;Introduction;Gyrokinetic approach;Charged particle motion in constant uniform magnetic filed;Guiding center motion in the inhomogeneous magnetic filed; Gyrocenter equations of motion [Brizard and Hahm, RMP‘07]:;Gyrokinetic Vlasov-Maxwell’s system of equations;Particle-in-cell method;Delta-f method;Fluid-kinetic hybrid electron model - I;Fluid-kinetic hybrid electron model - II;Fluid-kinetic hybrid electron model - III;Tokamak geometry;Magnetic coordinates - definition;Magnetic coordinates - example;GTC grids;Simulation grid;Parallel computing;MPI-Parallelization;OpenMP-Parallelization;;main.F90;main.F90 (cont); if(numereq==1)then ! Cyclone case psiw=0.0375 !poloidal flux at wall ! q and zeff profile is parabolic: q=q1+q2*psi/psiw+q3*(psi/psiw)^2 q=(/0.82,1.1,1.0/) !cyclone case er=(/0.0,0.0,0.0/) itemp0=1.0 ! on-axis thermal ion temperature, unit=T_e0 ftemp0=2.0 ! on-axis fast ion temperature, unit=T_e0 fden0=1.0e-5 ! on-axis fast ion density, unit=n_e0 ! density and temperature profiles are hyperbolic: ! ne=1.0+ne1*(tanh(ne2-(psi/psiw)/ne3)-1.0) ne=(/0.205,0.30,0.4/) !Cyclone base case te=(/0.415,0.18,0.4/) ti=(/0.415,0.18,0.4/) ! ITG tf=ti ! fast ion temperature profile nf=ne ! fast ion density profile;! allocate memory for equilibrium spline array allocate (stpp(lsp),mipp(lsp),mapp(lsp),nepp(3,lsp),tepp(3,lsp), nipp(3,lsp),tipp(3,lsp),zepp(3,lsp),ropp(3,lsp),erpp(3,lsp), qpsi(3,lsp),gpsi(3,lsp),ppsi(3,lsp),rpsi(3,lsp),torpsi(3,lsp), bsp(9,lsp,lst),xsp(9,lsp,l