GTC–GyrokineticToroidalCodeTrainingcoursePartIIhor.pptVIP

GTC–GyrokineticToroidalCodeTrainingcoursePartIIhor.ppt

  1. 1、本文档共28页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
GTC–GyrokineticToroidalCodeTrainingcoursePartIIhor.ppt

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

文档评论(0)

ailuojue + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档