摘要
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日(中国期刊网平台首次上网日期,不代表论文的发表时间)