摘要
Inordertoimplementlarge-scaleandhigh-betatokamaksimulation,anewalgorithmoftheelectromagneticgyrokineticPIC（particle-in-cell）codewasproposedandinstalledontheGpic-MHDcode[GyrokineticPICcodeformagnetohydrodynamic（MHD）simulation].Inthenewalgorithm,thevorticityequationandthegeneralizedOhm’slawalongthemagneticfieldarederivedfromthebasicequationsofthegyrokineticVlasov,Poisson,andAmperesystemandareusedtodescribethespatio-temporalevolutionofthefieldquantitiesoftheelectrostaticpotentialφandthelongitudinalcomponentofthevectorpotentialAz.Thebasicalgorithmisequivalenttosolvingthereduced-MHD-typeequationswithkineticcorrections,inwhichMHDphysicsrelatedtoAlfvenmodesarewelldescribed.Theestimationofperturbedelectronpressurefromparticledynamicsisdominant,whiletheeffectsofothermomentsarenegligible.Anotheradvantageofthealgorithmisthatthelongitudinalinducedelectricfield,Etz=-■z/■t,isexplicitlyestimatedbythegeneralizedOhm’slawandusedintheequationsofmotion.Furthermore,theparticlevelocitiesalongthemagneticfieldareused（vz-formulation）insteadofgeneralizedmomentums（pz-formulation）,hencethereisnoproblemof’cancellation’,whichwouldotherwiseappearwhenAzisestimatedfromtheAmpere’slawinthepz-formulation.ThesuccessfulsimulationofthecollisionlessinternalkinkmodebythenewGpic-MHDwithrealisticvaluesofthelarge-scaleandhigh-betatokamaksrevealedtheusefulnessofthenewalgorithm.
出版日期
2011年05月15日（中国期刊网平台首次上网日期，不代表论文的发表时间）
