CD-EXCALIBUR-FMS0058-M2.8.2_TechnicalPhysicsModelSelection ========================================================== .. meta:: :description: technical note :keywords: ExCALIBUR,Technical,report,on,Physics,model,selection,M2.8.2,The,report,describes,work,for,ExCALIBUR,project,NEPTUNE,at,Milestone,2.8.2.,This,report,provides,a,digest,of,the,reports,produced,by,the,Oxford,group,regarding,the,develop-,ment,of,a,drift,kinetic,model,for,plasma,in,the,tokamak,edge.,There,are,eight,reports,to,date,,five,of,which,concern,the,development,of,the,model:,“Physics,in,the,edge,of,fusion,devices,(2047357-,TN-03-02)”,[1];,“1D,drift,kinetic,models,with,periodic,boundary,conditions,(2047357-TN-01-02),[2];,“1D,drift,kinetic,models,with,wall,boundary,conditions,(2047357-TN-05),[3]”;,“2D,drift,kinetic,model,with,wall,boundary,conditions,(2047357-TN-07-01)”,[4];,and,“2D,drift,kinetic,model,with,periodic,boundary,conditions,(2047357-TN-09-01)”,[5].,The,remaining,three,reports,concern,the,numerical,implementation,of,the,model.,The,five,theory,reports,are,appended,to,this,report,in,their,logical,reading,order:,TN-03,,TN-01,,TN-05,,TN-07,,TN-09.,This,report,summarizes,the,five,reports,on,theory,development,,discussing,the,state,of,the,field,before,this,work,in,Section,1,,and,the,development,of,the,model,in,Section,2.,The,work,is,sum-,marized,and,future,work,is,discussed,in,Section,3.,1,Background,The,first,report,[1],describes,the,important,physical,regimes,in,the,tokamak,edge,plasma,and,reviews,the,state-of-the-art,of,edge,plasma,codes.,The,report,focusses,on,known,problems,of,fluid,models,,but,emphasises,the,importance,of,fluid,models,and,states,that,research,with,fluid,model,should,proceed,in,parallel,to,research,on,kinetic,models.,In,particular,,one,challenge,it,notes,for,kinetic,modelling,is,to,match,to,fluid,models,efficiently,,and,to,determine,which,spatial,regions,may,be,treated,using,fluid,models.,The,report,contains,four,topic,sections:,a,discussion,of,characteristic,time,and,length,scales,in,the,edge;,an,overview,of,the,existing,fluid,models,for,the,edge;,a,list,physical,phenomena,that,cannot,be,captured,by,fluid,models,(and,an,explanation,of,how,these,phenomena,are,currently,treated,by,fluid,models);,and,a,discussion,of,the,proposed,work,to,develop,a,complete,set,of,kinetic,equations,for,the,edge.,There,are,two,key,limitations,of,using,fluid,codes,to,model,edge,plasma.,The,first,is,that,the,calculated,diffusion,heat,flux,is,accurate,only,in,a,very,collisional,plasma.,This,is,because,the,heat,flux,is,predominantly,due,to,energetic,particles,which,are,much,less,collisional,than,thermal,particles.,Fluid,models,overestimate,the,number,of,energetic,particles,present,in,hot,plasmas,,while,underestimating,the,number,of,energetic,particles,present,in,cold,plasmas.,This,results,in,inaccurate,heat,fluxes.,The,second,limitation,concerns,boundary,conditions.,At,the,wall,,ions,and,electron,recombine,,meaning,that,near,the,wall,,the,distribution,function,vanishes,for,velocities,corresponding,to,motion,into,the,wall.,Thus,the,distribution,function,of,charged,particles,is,highly,non-Maxwellian,near,the,wall,and,must,be,treated,kinetically.,Together,,these,two,effects,result,in,an,incorrect,description,of,energetic,electrons,,which,in,turn,yields,incorrect,heat,transport,and,thus,incorrect,density,and,temperature,profiles.,Such,effects,have,been,seen,in,1D,kinetic,simulations,[6,,7].,Kinetic,models,for,neutrals,have,long,been,used,by,the,community,in,codes,such,as,EIRENE,[8],or,DEGAS,[9].,However,,energetic,electrons,also,mediate,important,collisional,processes,in,neu-,tral,models,,like,radiation,,ionization,and,recombination,,meaning,that,fluid-plasma,simulations,of,models,containing,neutrals,can,be,inaccurate,,even,when,kinetic,models,are,used,for,the,neutrals.,The,community,is,also,starting,to,develop,edge,code,that,treat,charged,particles,kinetically,,such,as,XGC,[10],,GKEYLL,[11],and,COGENT,[12].,However,,all,these,models,are,based,on,δf,gy-,rokinetics,[13,,14],,originally,developed,for,the,tokamak,core,,and,which,assume,that,turbulence,structures,are,larger,than,the,ion,gyroradius,,and,that,density,and,temperature,are,slowly,varying.,These,assumption,no,longer,hold,at,the,edge,,and,it,is,necessary,instead,to,use,“full,f,”,gyroki-,netics,[15,,16],where,the,distribution,function,is,no,longer,taken,to,be,a,small,perturbation,about,a,Maxwellian.,Full,f,gyrokinetics,is,more,general,formulation,,and,by,imposing,appropriate,addi-,tional,orderings,,one,may,recover,either,δf,gyrokinetics,for,the,core,,or,a,regime,appropriate,for,describing,the,edge,plasma.,The,latter,is,essentially,drift,kinetics,but,with,corrections,due,to,finite,gyroradius,effects.,These,finite,gyroradius,effects,are,important,in,two,places.,Firstly,,damping,terms,appear,in,the,kinetic,equation,which,stabilize,the,small,wavelength,instabilities,that,appear,in,drift,kinetics.,Secondly,,a,“polarization,density”,term,appears,in,the,field,equations,(quasineu-,trality,and,Amp,´ere’s,law).,This,term,is,often,small,,but,can,become,important,to,determine,the,electromagnetic,fields,in,certain,scenarios,,such,as,in,shear,Alfv,´en,waves.,While,implementing,the,damping,terms,is,relatively,straightforward,in,PIC,codes,,doing,so,has,not,been,attempted,in,full,f,continuum,codes.,Moreover,,the,correction,terms,in,the,field,equations,can,lead,to,significant,complications.,For,example,,to,include,these,corrections,while,also,ensuring,the,scheme,conserves,energy,and,momentum,would,lead,to,the,need,to,solve,a,nonlinear,equation,at,every,point,in,velocity,space.,In,addition,to,these,problems,,there,are,other,issues,a,model,must,address.,There,is,as,yet,no,gyrokinetic,formulation,for,the,wall,boundary,conditions,(only,a,much,less,general,formulation,for,ions,on,magnetic,field,lines,that,intersect,the,wall,at,a,shallow,angle,[17]).,It,also,remains,to,develop,a,model,for,collisions,between,charged,particles,and,neutrals,,which,is,challenging,in,a,gyrokinetic,framework,as,gyrophase,can,be,important,in,the,process.,Finally,,gyrokinetics,assumes,fluctuations,have,a,characteristic,scale,of,the,ion,gyroradius.,In,order,to,describe,Edge,Localized,Modes,–,significant,disruptions,in,the,plasma,at,the,tokamak,edge,–,it,is,necessary,to,extend,gyrokinetics,to,described,fluctuations,on,larger,MHD,scales.,2,1.1,Workplan,It,is,therefore,established,that,accurate,modelling,of,the,tokamak,edge,requires,the,inclusion,of,kinetic,effects.,A,tractable,model,will,also,use,gyroaveraging,to,remove,the,fastest,timescales,from,the,problem.,There,is,however,,no,existing,systematic,derivation,containing,both,of,these,for,plasma,in,the,tokamak,edge.,Making,such,a,derivation,is,the,main,thrust,of,this,work,package.,The,main,difficulty,is,performing,the,derivation,while,ensuring,that,the,equations,contain,all,the,information,required,to,completely,determine,the,system.,One,specific,problem,is,that,the,elec-,trostatic,potential,must,be,determined,from,the,quasineutrality,condition.,However,,even,small,deviations,from,charge,neutrality,cause,large,variations,in,the,electrostatic,potential,,leading,to,numerical,instabilities.,Thus,the,derived,system,must,be,amenable,to,implicit,time,stepping,,as,solving,quasineutrality,for,the,potential,explicitly,strongly,limits,possible,timesteps,[18].,The,derivation,introduces,a,novel,“moment,kinetic”,approach,,which,evolves,the,fluid,moments,of,the,distribution,function,separately,from,the,rest,of,the,distribution,function.,As,well,as,allowing,efficient,matching,between,fluid,and,kinetic,models,,this,particular,decomposition,proves,vital,for,the,derivation,of,a,tractable,system,in,2D.,The,work,plan,is,to,derive,this,model,,first,in,1D,with,periodic,and,wall,boundary,conditions,ap-,propriate,for,closed,and,open,field,line,regions,of,the,edge,respectively.,The,model,will,then,be,extended,to,3D,,but,using,an,axisymmetric,helical,field.,This,will,provide,a,model,that,is,improves,upon,models,implemented,in,existing,continuum,edge,codes.,The,model,is,limited,by,being,es-,sentially,2D,,and,for,using,a,simplified,helical,field,,rather,than,a,realistic,tokamak,field.,Both,of,these,are,candidates,for,future,work,,as,discussed,in,Section,3.1.,2,Theory,development,2.1,1D,drift,kinetics,with,periodic,boundary,conditions,The,first,report,“1D,drift,kinetic,models,with,periodic,boundary,conditions,(2047357-TN-01-02),[2],develops,a,proof-of-concept,for,the,moment,kinetic,approach.,It,uses,a,simple,1D,drift,kinetic,model,with,periodic,boundary,conditions,,appropriate,for,the,edge,region,inside,the,separatrix,that,has,closed,field,lines.,In,the,standard,kinetic,approaches,,one,solves,for,perturbations,about,a,Maxwellian,that,has,a,single,global,reference,density,,bulk,velocity,and,thermal,velocity.,In,the,moment-kinetic,approach,,one,instead,perturbs,about,a,Maxwellian,where,the,density,,bulk,velocity,and,thermal,velocity,are,themselves,variables,to,be,evolved,using,self-consistent,fluid,equations.,Thus,one,replaces,a,single,kinetic,equation,for,the,distribution,function,with,a,set,of,a,fluid,equations,and,a,kinetic,equation,for,a,modified,distribution,function.,The,main,motivation,for,this,approach,is,the,hope,that,it,would,allow,the,electrostatic,potential,to,be,derived,through,a,vorticity,equation,(using,fluid,variables),rather,than,through,a,field,solve,(requiring,kinetic,variables).,This,simplifies,the,solution,procedure,,but,requires,further,theory,de-,velopment.,A,secondary,benefit,of,splitting,out,the,moments,is,that,it,helps,the,numerics,when,3,dealing,with,dramatic,change,in,thermal,velocity,from,pedestal,to,divertor.,This,approach,normal-,izes,the,velocity,variable,in,the,modified,distribution,function,to,the,local,thermal,velocity,,meaning,no,regridding,is,necessary,for,different,regimes,in,the,plasma.,Finally,,this,approach,also,allows,one,to,develop,an,adaptive,model,which,would,automatically,switch,from,solving,the,fluid+modified,kinetic,equations,to,only,fluid,equations,when,the,distribution,function,was,sufficiently,close,to,Maxwellian.,This,may,prove,to,be,a,very,natural,and,computationally,tight,way,to,couple,a,fluid,model,to,a,kinetic,model.,2.2,1D,drift,kinetics,with,wall,boundary,conditions,In,the,second,report,[3],,the,1+1D,drift,kinetic,model,is,extended,so,that,it,can,treat,the,wall,bound-,ary,conditions,of,the,open,field,line,region,of,the,edge.,The,wall,is,assumed,to,be,perpendicular,to,the,magnetic,field,lines,so,that,simplified,boundary,conditions,,the,logical,sheath,boundary,conditions,[19],can,be,used.,The,logical,boundary,conditions,account,for,a,thin,sheath,of,non-neutral,plasma,of,width,order,Debye,length,in,front,of,the,wall.,There,is,a,potential,drop,across,this,sheath,which,repels,electrons,–,otherwise,electrons,(with,their,lower,mass,and,higher,thermal,velocity),would,flow,into,the,wall,at,a,greater,rate,than,ions,,violating,quasineutrality.,Requiring,current,into,the,sheath,to,vanish,gives,a,condition,for,the,electrostatic,potential,on,the,boundaries.,This,,in,conjunction,with,the,electron,parallel,momentum,equation,,can,be,used,to,find,the,electrostatic,potential,across,the,whole,domain.,As,before,,the,low,moments,of,the,distribution,function,are,treated,separately,so,that,for,ions,and,neutrals,three,fluid,equations,and,one,kinetic,equation,for,the,modified,distribution,function,must,be,evolved.,For,electrons,,two,fluid,equations,and,a,kinetic,equation,must,be,evolved,,with,the,third,fluid,equation,for,parallel,momentum,being,used,to,determine,the,electrostatic,potential.,2.3,2D,drift,kinetic,model,with,wall,boundary,conditions,The,third,report,[4],extends,the,previous,model,with,wall,boundary,conditions,to,a,2D,model,with,an,axisymmetric,helical,magnetic,field.,While,being,simpler,,the,helical,field,has,similarities,to,the,magnetic,field,in,the,edge,of,a,tokamak.,As,the,magnetic,field,is,now,not,perpendicular,to,the,walls,,it,is,in,principle,necessary,to,consider,the,magnetic,presheath,[20].,However,,the,boundary,conditions,that,this,would,entail,in,drift,kinet-,ics,are,both,complicated,and,an,area,of,open,research,[17,,21].,To,avoid,this,,it,is,assumed,that,the,electron,gyroradius,is,much,smaller,than,the,Debye,length,,in,which,case,boundary,conditions,similar,to,the,logical,boundary,conditions,[19],can,again,be,applied.,In,this,configuration,,the,ion,and,electron,bulk,properties,only,vary,along,field,lines,,so,as,before,,the,system,comprises,three,fluid,equations,and,a,kinetic,equations,for,ions,,and,two,fluid,equa-,tions,and,a,kinetic,equation,for,electrons.,As,before,,the,parallel,electron,momentum,equation,is,used,along,with,boundary,conditions,to,determine,the,electrostatic,potential.,Now,however,,the,distribution,functions,depend,on,two,spatial,and,two,velocity,space,coordinates.,The,velocity,coor-,4,dinates,are,parallel,and,perpendicular,to,the,magnetic,field,,with,the,third,direction,,the,gyrophase,,removed,by,gyroaveraging.,As,neutrals,are,unaffected,by,the,magnetic,field,,the,problem,for,neutrals,is,now,essentially,in,2+3D.,There,are,thus,five,fluid,equations,to,solve,(for,density,and,thermal,velocity,,plus,three,bulk,velocity,components).,The,kinetic,equation,for,neutrals,now,also,depends,on,two,spatial,directions,and,three,velocity,space,directions.,2.4,2D,drift,kinetic,model,with,periodic,boundary,conditions,The,fourth,report,[5],derives,a,2D,model,with,periodic,boundary,conditions,relevant,for,closed,flux,surfaces,,to,complement,the,wall,boundary,conditions,derived,in,a,previous,report.,That,report,described,how,to,obtain,the,electrostatic,potential,and,the,parallel,electron,bulk,velocity,with,wall,boundary,conditions.,The,same,prescription,with,periodic,boundary,conditions,only,allows,one,to,find,the,electron,velocity;,to,obtain,the,electrostatic,potential,,higher-order,terms,must,be,included,in,the,current,conservation,equations.,The,main,work,of,this,report,is,deriving,these,higher-order,corrections,,and,consequently,a,final,,additional,equation,(for,current,conservation,at,higher,order),to,solve,the,electrostatic,potential.,The,derivation,of,this,additional,equation,is,only,possible,in,the,moment-kinetic,framework,,which,evolves,the,ion,and,electron,densities,independently,of,their,normalized,distribution,functions.,With,unnormalized,distribution,functions,,the,corresponding,higher-order,current,conservation,equation,is,not,consistent,with,the,density,found,in,the,lower-order,kinetic,equations.,3,Summary,This,body,of,work,has,rigorously,and,systematically,derived,2D,drift,kinetic,equations,,appropriate,for,the,kinetic,treatment,of,ions,,electrons,and,neutral,in,the,tokamak,edge.,There,is,one,more,report,on,theory,development,to,be,submitted,under,this,grant,,“2D,drift,kinetics,in,a,helical,field,with,a,‘separatrix”’,(which,is,currently,under,review),and,which,combines,the,2D,drift,kinetic,models,for,the,open,and,closed,field,line,regions,into,a,single,model.,This,will,complete,the,theoretical,work,under,the,current,grant.,The,theoretical,developments,have,been,accompanied,by,the,numerical,implementation,of,these,equations,,using,a,Chebyshev,spectral,finite,element,implementation,in,Julia.,The,code,develop-,ments,are,described,in,separate,reports,[22,,23,,24],,and,the,software,is,publicly,available,[25].,At,present,,this,code,only,implements,a,subset,of,the,models,presented,here;,future,reports,will,discuss,the,software,developments,when,the,code,is,more,mature.,3.1,Future,directions,There,are,a,number,of,directions,for,future,work.,The,first,is,to,extend,the,model,into,3D,with,realistic,tokamak,geometry.,The,main,difficulty,here,is,the,inclusion,of,drifts,in,the,direction,per-,pendicular,to,the,flux,surfaces.,These,lead,to,finite-width,particle,orbits,which,can,significantly,5,alter,the,physics,in,quiescent,2D,[26],and,turbulent,3D,[27],plasmas.,The,rigorous,derivation,of,these,drifts,requires,further,work.,The,resulting,3D,model,is,also,expected,to,have,instabilities,on,the,grid,scale.,While,the,model,could,be,stabilized,using,numerical,dissipation,on,the,grid,scale,,further,work,is,required,to,prop-,erly,derive,the,finite,ion,gyroradius,effects,that,determine,behaviour,of,fluctuations,at,this,scale.,This,is,particularly,important,for,studying,high,confinement,mode,(H-mode),where,the,stabiliza-,tion,of,large-scale,turbulence,is,not,well,understood.,In,H-mode,,turbulence,at,the,ion,gyroscale,and,smaller,is,important,[27],,and,so,the,accurate,capture,of,finite,ion,gyroradius,effects,is,vital.,Including,these,finite,gyroradius,effects,has,the,potential,to,be,both,theoretically,and,numerically,challenging.,The,wall,boundary,conditions,in,the,model,are,also,simplified,conditions,which,assume,the,mag-,netic,field,is,perpendicular,to,the,wall.,In,the,edge,,this,is,often,not,the,case,,with,magnetic,field,lines,meeting,the,wall,at,very,shallow,angles.,The,boundary,conditions,are,formally,valid,in,the,limit,of,electron,gyroradius,being,much,smaller,than,the,Debye,length,,but,relaxing,this,condition,requires,further,study.,This,is,currently,being,pursued,by,Parra,and,coworkers,outside,ExCAL-,IBUR.,Finally,,it,would,be,beneficial,to,extend,the,model,to,treat,fluctuations,in,the,magnetic,field,,as,well,as,electric,field.,While,electrostatic,models,are,often,a,good,approximation,for,edge,plasmas,,including,electromagnetic,effects,would,allow,this,model,to,connect,with,the,MHD,regime,,and,so,to,study,large-scale,plasma,disturbances.,6,Acknowledgement,The,support,of,the,UK,Meteorological,Office,and,Strategic,Priorities,Fund,is,acknowledged.,References,[1],F.,I.,Parra,,M.,Barnes,,and,M.,Hardman.,Physics,in,the,edge,of,fusion,devices.,Tech-,nical,Report,2047357-TN-03-02,,UKAEA,Project,Neptune,,2021.,https://github.com/,ExCALIBUR-NEPTUNE/Documents/blob/main/reports/2047357/TN-03.pdf.,[2],F.I,Parra,,M.,Barnes,,and,M.R.,Hardman.,1D,drift,kinetic,models,with,periodic,boundary,conditions.,Technical,Report,2047357-TN-01-02,,UKAEA,Project,Neptune,,2021.,https:,//github.com/ExCALIBUR-NEPTUNE/Documents/blob/main/reports/2047357/TN-01.pdf.,[3],F.,I.,Parra,,M.,Barnes,,and,M.,Hardman.,1D,drift,kinetic,models,with,wall,boundary,conditions.,Technical,Report,2047357-TN-05-01,,UKAEA,Project,Neptune,,2021.,https://github.com/,ExCALIBUR-NEPTUNE/Documents/blob/main/reports/2047357/TN-05.pdf.,[4],F.,I.,Parra,,M.,Barnes,,and,M.,Hardman.,2D,drift,kinetic,models,with,wall,boundary,conditions.,Technical,Report,2047357-TN-07-01,,UKAEA,Project,Neptune,,2021.,https://github.com/,ExCALIBUR-NEPTUNE/Documents/blob/main/reports/2047357/TN-07.pdf.,[5],F.,I.,Parra,,M.,Barnes,,and,M.,Hardman.,2D,drift,kinetic,models,with,periodic,boundary,conditions.,Technical,Report,2047357-TN-09-01,,UKAEA,Project,Neptune,,2021.,https:,//github.com/ExCALIBUR-NEPTUNE/Documents/blob/main/reports/2047357/TN-09.pdf.,[6],D.,Tskhakaya,,S.,Jachmich,,T.,Eich,,W.,Fundamenski,,and,JET,EFDA,contributors.,Interpre-,tation,of,divertor,Langmuir,probe,measurements,during,the,ELMs,at,JET.,Journal,of,Nuclear,Materials,,415(1):S860–S864,,2011.,[7],A.,V.,Chankin,and,D.,P.,Coster.,On,the,locality,of,parallel,transport,of,heat,carrying,electrons,in,the,SOL.,Journal,of,Nuclear,Materials,,463:498–501,,2015.,[8],D.,Reiter,,M.,Baelmans,,and,P.,Boerner.,The,EIRENE,and,B2-EIRENE,codes.,Fusion,Science,and,Technology,,47(2):172–186,,2005.,[9],D.,Stotler,and,C.,Karney.,Neutral,gas,transport,modeling,with,DEGAS,2.,Contributions,to,Plasma,Physics,,34(2-3):392–397,,1994.,[10],S.,Ku,,R.,Hager,,C.-S.,Chang,,J.,M.,Kwon,,and,S.,E.,Parker.,A,new,hybrid-Lagrangian,nu-,merical,scheme,for,gyrokinetic,simulation,of,tokamak,edge,plasma.,Journal,of,Computational,Physics,,315:467–475,,2016.,[11],A.,H.,Hakim,,N.,R.,Mandell,,T.,N.,Bernard,,M.,Francisquez,,G.,W.,Hammett,,and,E.,L.,Shi.,Continuum,electromagnetic,gyrokinetic,simulations,of,turbulence,in,the,tokamak,scrape-off,layer,and,laboratory,devices.,Physics,of,Plasmas,,27(4):042304,,2020.,7,[12],M.,A.,Dorf,,M.,R.,Dorr,,J.,A.,Hittinger,,R.,H.,Cohen,,and,T.,D.,Rognlien.,Continuum,kinetic,modeling,of,the,tokamak,plasma,edge.,Physics,of,Plasmas,,23(5):056102,,2016.,[13],P.,J.,Catto.,Linearized,gyro-kinetics.,Plasma,Physics,,20(7):719,,1978.,[14],E.,A.,Frieman,and,L.,Chen.,Nonlinear,gyrokinetic,equations,for,low-frequency,electromag-,netic,waves,in,general,plasma,equilibria.,The,Physics,of,Fluids,,25(3):502–508,,1982.,[15],F.,I.,Parra,and,P.,J.,Catto.,Limitations,of,gyrokinetics,on,transport,time,scales.,Plasma,Physics,and,Controlled,Fusion,,50(6):065014,,2008.,[16],A.,J.,Brizard,and,T.,S.,Hahm.,Foundations,of,nonlinear,gyrokinetic,theory.,Reviews,of,modern,physics,,79(2):421–468,,2007.,[17],A.,Geraldini,,F.I.,Parra,,and,F.,Militello.,Dependence,on,ion,temperature,of,shallow-angle,magnetic,presheaths,with,adiabatic,electrons.,Journal,of,Plasma,Physics,,85(6),,2019.,[18],W.W.,Lee.,Gyrokinetic,particle,simulation,model.,Journal,of,Computational,Physics,,72(1):243–269,,1987.,[19],S.,E.,Parker,,R.,J.,Procassini,,C.,K.,Birdsall,,and,B.,I.,Cohen.,A,suitable,boundary,condition,for,bounded,plasma,simulation,without,sheath,resolution.,Journal,of,Computational,Physics,,104(1):41–49,,1993.,[20],R.,Chodura.,Plasma-wall,transition,in,an,oblique,magnetic,field.,The,Physics,of,Fluids,,25(9):1628–1633,,1982.,[21],A.,Geraldini.,Large,gyro-orbit,model,of,ion,velocity,distribution,in,plasma,near,a,wall,in,a,grazing-angle,magnetic,field.,Journal,of,Plasma,Physics,,87(1),,2021.,[22],M.,Barnes,,F.,I.,Parra,,and,M.,R.,Hardman.,Numerical,study,of,1D,drift,kinetic,models,with,periodic,boundary,conditions.,Technical,Report,2047357-TN-02-02,,UKAEA,Project,Nep-,tune,,2021.,https://github.com/ExCALIBUR-NEPTUNE/Documents/blob/main/reports/,2047357/TN-02.pdf.,[23],M.,Barnes,,F.,I.,Parra,,M.,R.,Hardman,,and,J.,Omotani.,Numerical,study,of,1+1D,,moment-,based,drift,kinetic,models,with,periodic,boundary,conditions.,Technical,Report,2047357-TN-,04,,UKAEA,Project,Neptune,,2021.,https://github.com/ExCALIBUR-NEPTUNE/Documents/,blob/main/reports/2047357/TN-04.pdf.,[24],M.,Barnes,,F.,I.,Parra,,M.,R.,Hardman,,and,J.,Omotani.,Numerical,study,of,1+1D,drift,kinetic,model,with,wall,boundary,conditions.,Technical,Report,2047357-TN-08,,UKAEA,Project,Neptune,,2021.,https://github.com/ExCALIBUR-NEPTUNE/Documents/blob/main/,reports/2047357/TN-08.pdf.,[25],M.,A.,Barnes,and,J.,T.,Omotani.,moment,kinetics,GitHub,repository.,https://github.com/,mabarnes/moment_kinetics,,2021.,Accessed:,December,2021.,[26],G.,Kagan,and,P.,J.,Catto.,Arbitrary,poloidal,gyroradius,effects,in,tokamak,pedestals,and,transport,barriers.,Plasma,Physics,and,Controlled,Fusion,,50(8):085010,,2008.,8,[27],J.,F.,Parisi,,F.,I.,Parra,,C.,M.,Roach,,C.,Giroud,,W.,Dorland,,D.,R.,Hatch,,M.,Barnes,,J.,C.,Hillesheim,,N.,Aiba,,J.,Ball,,et,al.,Toroidal,and,slab,ETG,instability,dominance,in,the,linear,spectrum,of,JET-ILW,pedestals.,Nuclear,Fusion,,60(12):126045,,2020.,UKAEA,REFERENCE,AND,APPROVAL,SHEET,Client,Reference:,UKAEA,Reference:,CD/EXCALIBUR-FMS/0058,Issue:,Date:,1.0,9,December,2021,Project,Name:,ExCALIBUR,Fusion,Modelling,System,Prepared,By:,Name,and,Department,Joseph,Parker,Wayne,Arter,Signature,N/A,N/A,Date,9,December,2021,9,December,2021,BD,Reviewed,By:,Rob,Akers,9,December,2021,Advanced,Dept.,Manager,Computing,Approved,By:,Rob,Akers,9,December,2021,Advanced,Dept.,Manager,Computing,9,REPORT,2047357-TN-03-02,M1.2,1,Physics,in,the,edge,of,fusion,devices,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Rudolf,Peierls,Centre,for,Theoretical,Physics,,University,of,Oxford,,Oxford,OX1,3PU,,UK,(This,version,is,of,23,April,2021),1.,Introduction,We,present,a,brief,overview,of,the,current,state,of,fusion,device,edge,modeling.,This,re-,port,is,written,for,ExCALIBUR,NEPTUNE,contract,T/NA085/20.,This,contract,stated,mission,is,the,development,of,kinetic,models,for,the,edge.,Thus,,this,report,will,focus,on,the,known,problems,of,fluid,models,without,making,much,emphasis,on,their,many,suc-,cesses.,This,lack,of,praise,for,fluid,models,is,not,the,objective,of,this,report,and,hence,we,would,like,to,start,by,reassuring,the,reader,that,we,believe,current,fluid,codes,have,much,to,offer,and,need,to,be,pursued,in,parallel,to,kinetic,modeling.,One,of,the,challenges,for,kinetic,modeling,is,to,devise,methods,to,match,with,fluid,simulations,efficiently.,With,this,capability,,the,edge,can,be,split,into,non-overlapping,spatial,regions,with,different,levels,of,kinetic,sophistication.,The,remainder,of,this,reports,is,organized,as,follows.,In,section,2,we,explain,what,the,characteristic,time,and,length,scales,are,in,the,edge.,In,section,3,,we,give,an,overview,of,the,existing,fluid,models,for,the,edge.,In,section,4,we,list,physical,phenomena,that,cannot,be,captured,by,fluid,models,,and,we,explain,how,these,phenomena,have,been,addressed,so,far.,Finally,,we,discuss,our,proposal,to,develop,a,complete,set,of,kinetic,equations,for,the,edge,in,section,5.,2.,The,edge,In,this,report,,we,call,the,edge,the,plasma,that,surrounds,the,separatrix.,The,separatrix,(represented,as,a,red,line,in,figure,1),is,the,flux,surface,that,separates,the,region,where,magnetic,field,lines,are,in,contact,with,the,walls,of,the,vessel,from,the,region,in,which,magnetic,field,lines,form,nested,toroidal,flux,surfaces.,This,means,that,the,edge,includes,both,open,field,lines,(those,who,are,in,contact,with,walls),and,closed,field,lines,(those,that,form,nested,toroidal,flux,surfaces).,We,need,to,distinguish,the,core,,where,the,fusion,reactions,are,supposed,to,take,place,,from,the,edge.,In,this,report,,we,will,use,a,theoretical,criterion,to,distinguish,one,from,the,other:,whether,the,characteristic,transit,times,in,the,directions,parallel,and,perpendicular,to,the,magnetic,field,are,comparable,or,not.,In,the,core,,both,ions,and,electron,closely,follow,magnetic,field,lines,and,can,travel,around,the,device,many,times,before,collisions,,turbulent,fluctuations,or,other,effects,drive,them,away,from,the,magnetic,field,line,on,which,they,started.,As,a,result,,density,and,temperature,are,almost,constant,along,magnetic,field,lines,,and,one,only,worries,about,small,fluctuations,around,the,mostly,quiescent,profiles,of,density,and,temperature.,This,is,,of,course,,a,highly,idealized,situation,that,ignores,violent,events,taking,place,in,the,core,,such,as,sawteeth,(Hastie,1997),,but,it,is,a,useful,one.,In,contrast,,in,the,edge,,charged,particles,that,were,well-confined,in,the,core,region,cross,the,separatrix,and,eventually,reach,the,wall,by,following,magnetic,field,lines.,The,displacements,of,ions,and,electrons,away,from,the,magnetic,field,line,in,which,they,2,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Figure,1.,Sketch,of,a,tokamak,edge,on,a,plane,that,contains,the,axis,of,symmetry,(represented,here,by,a,dash-dot,vertical,line).,The,edge,region,is,in,blue,,and,the,separatrix,is,the,red,line.,started,are,still,small,compared,to,the,size,of,the,device,,but,they,are,not,small,compared,with,the,characteristic,lengths,of,the,density,and,temperature,in,the,edge,region.,If,the,characteristic,lengths,of,density,and,temperature,along,and,across,magnetic,field,lines,are,L(cid:107),and,L⊥,,the,edge,is,characterized,by,|v(cid:107)|,L(cid:107),∼,|vd|,L⊥,,,(2.1),where,v(cid:107),is,the,characteristic,parallel,velocity,of,the,particles,,and,vd,their,small,drift,perpendicular,to,the,magnetic,field.,Since,|v(cid:107)|,(cid:29),|vd|,,L⊥,is,much,smaller,than,L(cid:107).,In,the,core,,for,comparison,,we,find,that,L(cid:107),and,L⊥,are,of,similar,order,and,comparable,to,the,machine,size,,giving,|v(cid:107)|,L(cid:107),(cid:29),|vd|,L⊥,.,(2.2),There,is,another,aspect,that,makes,the,edge,very,different,from,the,core.,The,tem-,perature,of,both,electrons,and,ions,is,kept,low,near,the,wall,because,the,wall,is,a,very,effective,sink,of,energy,(wall,materials,that,prevent,slow,Hydrogen,particles,from,going,back,into,the,plasma,are,notable,exceptions,where,the,temperature,of,the,plasma,can,be,large,near,the,wall;,see,,for,example,,Schmitt,et,al.,(2015),for,liquid,Lithium,divertors,,or,Jackson,et,al.,(1991),for,boronized,walls).,At,low,plasma,temperatures,,the,plasma,is,partially,ionized,and,collisions,between,the,charged,particles,in,the,plasma,and,neutrals,become,important.,In,extreme,limits,,detachment,occurs,,that,is,,the,plasma,tempera-,ture,decreases,sufficiently,due,to,radiation,that,the,plasma,recombines,and,a,cushion,of,neutrals,appears,in,front,of,the,walls,protecting,them,(Krasheninnikov,&,Kukushkin,2017).,In,addition,to,limiting,the,temperature,,the,presence,of,the,wall,controls,the,size,of,the,electric,field,and,the,flows,in,the,open,field,line,region.,We,will,discuss,these,effects,in,more,detail,in,section,3.,We,finish,this,section,by,calculating,a,few,characteristic,time,and,length,scales.,In,Militello,&,Fundamenski,(2011),,one,can,find,a,concise,summary,of,typical,values,of,plasma,characteristics,in,the,edge,of,tokamaks.,In,current,tokamaks,,the,magnitude,of,Closedfieldline,regionOpen,field,line,regionWalls,Physics,in,the,edge,of,fusion,devices,3,the,magnetic,field,B,is,typically,2,T,,and,the,characteristic,size,of,the,device,is,a,few,meters.,The,plasma,temperature,T,ranges,from,1,keV,in,the,closed,magnetic,field,line,part,of,the,edge,to,100,eV,at,the,separatrix,down,to,10,eV,near,the,wall.,In,detached,plasmas,,the,temperature,drops,to,1,eV,near,the,wall.,The,electron,density,ne,ranges,from,1019,m−3,in,the,open,field,line,region,to,1020,m−3,in,the,closed,field,line,region.,Neutral,density,nn,is,usually,smaller,than,ne,and,it,ranges,from,1015,m−3,in,the,closed,field,line,region,to,1019,m−3,in,the,open,field,line,region,(Colchin,et,al.,2000;,Scotti,et,al.,2021).,With,these,quantities,,we,calculate,several,characteristic,frequencies,of,interest:,•,The,gyrofrequencies,of,both,Deuterium,and,electrons,,Ωi,:=,eB/mD,and,Ωe,:=,eB/me,,are,the,characteristic,frequencies,of,the,nearly,circular,motion,of,the,charged,particles,around,magnetic,field,lines.,Here,e,is,the,proton,charge,,and,mD,and,me,are,the,Deuterium,and,electron,masses.,•,The,transit,frequencies,for,Deuterium,ions,and,electrons,,vtD/L(cid:107),and,vte/L(cid:107),,are,the,inverse,of,the,time,that,it,takes,charged,particles,to,move,along,a,magnetic,field,line,from,one,wall,to,another,in,the,open,field,line,region,,and,the,inverse,of,the,time,that,it,takes,charged,particles,to,sample,a,flux,surface,in,the,closed,field,line,region.,Here,vtD,:=,(cid:112)2T,/mD,and,vte,:=,(cid:112)2T,/me,are,the,Deuterium,and,electron,thermal,speeds,,and,L(cid:107),is,the,characteristic,length,of,magnetic,field,lines,in,the,edge,,which,we,take,to,be,10,m.,•,The,transit,frequency,for,Deuterium,neutrals,,vtD/L⊥,,is,the,inverse,of,the,time,that,it,takes,neutral,atoms,to,cross,the,edge,region.,The,characteristic,scale,of,variation,of,density,and,temperature,across,magnetic,field,lines,,L⊥,,is,of,the,order,of,5,cm,or,larger,in,the,closed,field,line,region,(Sugihara,et,al.,2000),,and,of,the,order,of,1,cm,or,larger,in,the,open,field,line,region,(Goldston,2012).,•,The,collision,frequencies,and,νii,:=,√,4,π,3,e4ne,ln,Λ,(4π(cid:15)0)2m1/2,D,T,3/2,νep,:=,√,4,2π,3,e4ne,ln,Λ,(4π(cid:15)0)2m1/2,e,T,3/2,(2.3),(2.4),describe,how,often,Deuterium,ions,collide,with,each,other,or,electrons,collide,with,other,electrons,or,Deuterium,ions,,respectively.,Here,ln,Λ,≈,15,is,the,Coulomb,logarithm,and,(cid:15)0,2,difference,between,νii,and,νep,is,a,convention,is,the,vacuum,permittivity.,The,factor,of,introduced,by,Braginskii,(Braginskii,1958).,√,•,We,also,need,the,collision,frequencies,that,describe,how,often,Deuterium,ions,collide,with,Deuterium,neutral,atoms,,νin,:=,nnvtDσin,,and,how,often,Deuterium,neutral,atoms,collide,with,Deuterium,ions,,νni,:=,nivtDσin.,Here,,the,ion-neutral,cross,section,σin,is,of,order,10−18,m2,(Lindsay,&,Stebbings,2005).,Similarly,,we,need,the,collision,frequencies,that,describes,how,often,electrons,collide,with,Deuterium,neutral,atoms,,νen,:=,nnvteσen,,and,how,often,neutrals,are,ionized,,νion,:=,nevteσion.,Here,,the,electron-neutral,collision,cross,section,σen,is,of,order,10−19,m2,(Brackmann,et,al.,1958),and,the,ionization,collision,cross,section,σion,is,of,order,10−20,m2,(Zel’dovich,&,Raizer,2013).,All,these,frequencies,are,shown,in,table,1.,It,is,clear,that,the,collision,frequencies,are,the,ones,that,change,the,most,across,the,edge.,For,ions,and,electrons,,collisions,become,very,important,in,the,cooler,plasma,of,the,open,field,line,region,,but,are,infrequent,in,the,closed,field,line,region.,Conversely,,for,neutrals,,collision,are,more,frequent,in,the,closed,field,line,region,than,in,the,open,field,line,region.,We,also,calculate,a,few,length,scales,of,interest:,4,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Ion,freq.,[104,s−1],Ωi,Region,Closed,lines,9400,9400,Open,lines,vtD/L(cid:107),3.1,0.31,νin,νii,0.16,0.031,16,31,Neutral,freq.,[104,s−1],vtD/L⊥,νni,3100,31,νion,1800,18,610,310,Electron,freq.,[106,s−1],Ωe,νen,vte/L(cid:107),νep,320000,320000,0.13,0.0018,13,1.8,0.18,1.8,Table,1.,Characteristic,frequencies,for,the,closed,field,line,region,(B,=,2,T,,ne,=,1020,m−3,,nn,=,1015,m−3,,T,=,1,keV),and,the,open,field,line,region,(B,=,2,T,,ne,=,1019,m−3,,nn,=,1019,m−3,,T,=,10,eV).,Parallel,lengths,[1,m],L(cid:107),λii,Region,Closed,lines,10,190,1000,140,10000,Open,lines,λen,L⊥,5,1,λin,λep,10,0.19,0.14,0.1,1,Perpendicular,lengths,[1,cm],ρe,ρi,0.33,0.0056,0.033,0.00056,λni,1,10,λion,1.7,17,Table,2.,Characteristic,lengths,for,the,closed,field,line,region,(B,=,2,T,,ne,=,1020,m−3,,nn,=,1015,m−3,,T,=,1,keV),and,the,open,field,line,region,(B,=,2,T,,ne,=,1019,m−3,,nn,=,1019,m−3,,T,=,10,eV).,•,The,mean,free,paths,λii,:=,vtD/νii,and,λin,:=,vtD/νin,are,the,distances,that,a,Deu-,terium,ion,can,travel,before,colliding,with,another,Deuterium,ion,or,with,a,Deuterium,atom,,respectively.,Similarly,,the,mean,free,paths,λni,:=,vtD/νni,and,λion,:=,vtD/νion,are,the,distances,that,a,Deuterium,neutral,atom,can,move,before,colliding,with,an,ion,or,getting,ionized,,respectively,,and,λep,:=,vte/νep,and,λen,:=,vte/νen,are,the,distances,that,an,electron,can,move,before,colliding,with,another,charged,particle,or,with,a,neutral,,respectively.,•,The,Deuterium,and,electron,gyroradii,,ρi,:=,vtD/Ωi,and,ρe,:=,vte/Ωe,,are,the,characteristic,size,of,the,gyration,of,charged,particles,around,magnetic,field,lines.,All,these,lengths,are,given,in,table,2.,Unsurprisingly,,we,see,that,the,mean,free,paths,,inversely,proportional,to,the,collision,frequencies,,are,the,characteristic,lengths,that,change,the,most,across,the,edge.,In,the,direction,perpendicular,to,the,magnetic,field,,we,see,that,the,gyroradii,are,small,compared,to,the,characteristic,lengths.,3.,Drift-ordered,fluid,models,Due,to,the,large,collision,frequencies,in,the,open,field,line,region,,the,bulk,of,the,ion,and,electron,distribution,functions,is,Maxwellian.,For,this,reason,,much,of,the,edge,modeling,has,been,based,on,plasma,fluid,equations,derived,in,the,limit,of,large,Coulomb,collisions,(Braginskii,1958).,Since,the,most,readily,available,fluid,equations,for,magnetized,plasmas,assume,that,there,are,no,neutrals,,most,fluid,models,ignore,neutrals,,with,some,notable,exceptions.,The,main,difference,between,the,fluid,equations,used,in,edge,models,and,the,usual,fluid,equations,is,that,diffusivities,are,anisotropic.,As,demonstrated,by,table,1,,charged,parti-,cles,gyrate,around,magnetic,field,lines,many,times,before,having,a,collision.,Thus,,parti-,cles,barely,move,across,magnetic,field,lines,and,the,diffusivity,across,magnetic,field,lines,is,much,smaller,than,the,diffusivity,along,them.,To,capture,the,effect,of,this,anisotropy,correctly,and,efficiently,,one,has,to,either,use,flux,coordinates,that,follow,magnetic,field,lines,(Beer,et,al.,1995),or,employ,appropriate,discretizations,(Hariri,&,Ottaviani,2013).,Physics,in,the,edge,of,fusion,devices,5,The,problem,with,following,magnetic,field,lines,is,that,it,is,difficult,to,find,a,grid,that,both,aligns,with,the,magnetic,field,lines,and,extends,to,the,walls,of,the,vessel.,Unfor-,tunately,,the,wall,geometry,is,important,because,it,determines,how,far,neutrals,leaving,the,wall,travel,into,the,plasma,(Wiesen,et,al.,2018).,Finding,grids,that,adjust,to,the,wall,geometry,has,become,a,problem,of,great,interest,in,recent,years,,and,the,community,has,been,trying,exotic,methods,to,address,it,(Isoardi,et,al.,2010;,Bufferand,et,al.,2019).,An,important,aspect,of,the,edge,fluid,models,is,the,size,of,the,flow.,The,potential,difference,between,the,wall,and,the,plasma,is,controlled,by,the,non-neutral,Debye,sheath,that,forms,around,the,walls,(Riemann,1991).,Typically,,the,potential,drop,between,the,wall,and,the,plasma,has,to,be,several,times,the,electron,temperature,because,otherwise,a,large,electron,current,flows,into,the,wall,,breaking,the,neutrality,of,the,plasma.,If,one,assumes,that,the,wall,in,contact,with,the,plasma,is,a,conductor,and,hence,the,potential,is,constant,across,its,volume,,the,potential,differences,within,the,plasma,are,restricted,to,be,of,the,order,of,the,electron,temperature,,φ,∼,T,/e.,This,gives,an,electric,field,E,=,−∇φ,of,the,order,of,T,/eL⊥.,The,perpendicular,velocity,of,the,fluid,is,determined,by,the,balance,between,the,electric,and,magnetic,forces,,This,equation,gives,a,perpendicular,velocity,of,order,u⊥,×,B,∼,E.,u⊥,∼,|E|,B,∼,ρi,L⊥,vtD.,(3.1),(3.2),Thus,,according,to,table,2,,u⊥,is,significantly,smaller,than,vtD.,The,fluid,equations,obtained,with,the,ordering,(3.2),are,known,as,drift-ordered,equations,because,the,flow,is,of,the,same,order,as,the,slow,particle,drifts,–,the,other,possible,ordering,is,the,high,flow,ordering,that,assumes,that,the,perpendicular,velocity,is,sonic.,Importantly,,for,fluid,velocities,of,the,size,given,in,equation,(3.2),,one,needs,to,keep,terms,that,depend,on,the,gradient,of,the,temperature,and,the,pressure,in,the,stress,tensor,to,be,completely,consistent,(Mikhailovskii,&,Tsypin,1971;,Simakov,&,Catto,2003;,Catto,&,Simakov,2004).,The,system,of,drift-ordered,fluid,equations,is,usually,comprised,of,•,one,continuity,equation,per,ion,species,(the,electron,density,need,not,be,calculated,because,it,is,determined,by,quasineutrality),,•,one,conservation,equation,for,the,component,of,the,plasma,momentum,parallel,to,the,magnetic,field,,•,a,vorticity,equation,that,determines,the,electrostatic,potential,,•,in,electromagnetic,models,,Amp´ere’s,law,,and,•,one,conservation,equation,for,the,energy,of,all,ion,species,and,another,one,for,the,energy,of,the,electrons.,Note,that,there,is,one,single,conservation,equation,for,the,whole,plasma,parallel,momen-,tum,and,one,single,conservation,equation,for,the,energy,of,all,the,ion,species,,and,not,several,conservation,equations,,one,per,ion,species.,The,reason,why,all,ion,species,must,be,considered,as,one,in,these,fluid,equations,is,that,,within,the,large,collision,frequency,approximation,,all,ion,species,have,the,same,temperature,and,all,charged,species,have,the,same,average,flow.,The,electron,temperature,can,be,different,from,the,ion,temperature,due,to,the,mass,difference,between,the,two,species.,The,temperature,differences,between,the,different,ion,species,are,of,the,same,order,as,kinetic,effects,that,are,neglected.,The,differences,between,the,parallel,flows,are,calculated,and,used,in,the,vorticity,equation,,where,they,are,needed,because,the,electric,current,enters,in,the,Lorenz,force.,As,we,mentioned,at,the,start,,the,drift-ordered,fluid,equations,usually,implemented,6,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,in,edge,codes,(Mikhailovskii,&,Tsypin,1971;,Simakov,&,Catto,2003;,Catto,&,Simakov,2004),ignore,neutrals.,There,have,been,attempts,to,include,neutrals,in,the,limit,where,both,ions,and,neutrals,are,sufficiently,collisional,that,both,species,can,be,treated,as,fluids,(Hazeltine,et,al.,1992;,Catto,1994;,Helander,et,al.,1994).,The,advantage,of,these,models,is,that,they,account,for,the,effect,that,ion-neutral,collisions,have,on,the,diffusivities,,but,it,is,generally,believed,that,these,models,are,insufficient,because,neutrals,have,long,mean,free,paths,in,large,regions,of,the,edge,(see,tables,1,and,2).,There,have,been,some,attempts,at,simplified,kinetic,treatments,of,the,neutrals,(Wersal,&,Ricci,2015),and,there,exist,sophisticated,Monte-Carlo,approaches,with,a,large,collection,of,collisions,such,as,EIRENE,(Reiter,et,al.,2005),or,DEGAS,(Stotler,&,Karney,1994).,Note,that,these,Monte-,Carlo,treatments,are,not,completely,self-consistent,because,they,assume,that,the,ion,and,electron,distribution,functions,are,Maxwellian,and,,depending,on,the,version,of,the,code,,they,ignore,elastic,collisions,and,they,average,over,the,dependence,of,the,differential,cross,section,on,the,scattering,angle.,Moreover,,these,neutral,kinetic,codes,are,sensitive,to,the,parameters,that,determine,how,particles,interact,with,the,wall,(Chankin,et,al.,2021),,and,these,parameters,are,not,well,known.,There,are,two,types,of,problems,in,which,the,drift-ordered,fluid,equations,are,used:,2D,profiles,and,turbulence.,The,objective,of,2D,fluid,solvers,such,as,SOLPS,(Wiesen,et,al.,2015),,SolEdge2D,(Bufferand,et,al.,2015),,UEDGE,(Rognlien,et,al.,2007),or,EDGE2D,(Simonini,et,al.,1994),is,to,determine,the,toroidally,averaged,density,and,temperature,profiles,in,the,edge.,These,codes,cannot,model,turbulent,fluctuation,because,they,are,missing,the,third,dimension.,For,this,reason,,instead,of,the,collisional,perpendicular,diffusion,coefficients,,these,codes,use,enhanced,perpendicular,diffusion,coefficients,that,are,chosen,to,match,the,experimental,observations.,These,2D,codes,are,meant,to,provide,detailed,understanding,of,transport,along,magnetic,field,lines,and,of,neutrals,,as,some,of,them,are,coupled,to,Monte-Carlo,neutral,codes,such,as,EIRENE,(Reiter,et,al.,2005),or,DEGAS,(Stotler,&,Karney,1994).,Turbulence,codes,such,as,GBS,(Halpern,et,al.,2016),,TOKAM3X,(Tamain,et,al.,2016),,Hermes,(Dudson,&,Leddy,2017),or,GRILLIX,(Stegmeir,et,al.,2018),are,3D,fluid,codes.,Originally,,fluid,turbulence,codes,assumed,that,ions,were,much,colder,than,electrons,and,that,the,turbulent,fluctuations,were,small,compared,to,an,almost,constant,background,(Zeiler,et,al.,1996).,However,,it,was,soon,realized,that,this,treatment,is,not,appropriate,for,the,edge.,Ions,are,not,cold,and,the,density,and,temperature,profiles,cannot,be,easily,split,into,a,slowly,varying,piece,plus,small,fluctuations,due,to,the,presence,of,the,wall.,The,same,wall,boundary,conditions,that,constrain,the,fluid,velocity,perpendicular,to,the,magnetic,field,to,be,subsonic,,as,shown,in,equation,(3.2),,require,that,the,fluid,velocity,parallel,to,the,magnetic,field,be,sonic,near,the,wall,(Chodura,1982).,The,pressure,drops,along,magnetic,field,lines,have,to,be,significant,to,ensure,that,the,parallel,flow,is,accelerated,sufficiently.,This,is,incompatible,with,the,assumption,that,the,fluctuations,are,small.,Thus,,most,current,edge,turbulence,codes,have,tried,to,lift,the,assumption,of,small,fluctuations,,and,they,have,done,so,by,differing,degrees,depending,on,the,code,or,the,version,of,the,code,in,use.,Overall,,fluid,codes,are,maturing,,and,although,there,is,still,work,to,be,done,,there,is,starting,to,be,a,consensus,on,the,physics,that,they,must,include.,The,same,cannot,be,said,about,kinetic,effects,in,the,edge.,4.,Kinetic,effects,Tables,1,and,2,are,evidence,that,kinetic,effects,must,be,taken,into,account:,in,the,closed,magnetic,field,line,region,and,in,part,of,the,open,field,line,region,,ions,and,electrons,can,Physics,in,the,edge,of,fusion,devices,7,travel,a,distance,of,the,size,of,the,device,without,suffering,a,single,collision.,Moreover,,kinetic,effects,are,surprisingly,important,even,in,the,regions,where,the,collision,frequency,is,high.,There,are,two,reasons,for,this.,•,The,diffusive,heat,flux,calculated,in,the,fluid,limit,is,only,accurate,if,νep,(cid:38),100,vte/L(cid:107),because,the,heat,flux,is,dominated,by,energetic,particles,that,collide,much,less,than,the,thermal,particles,(Gurevich,&,Istomin,1979;,Gray,&,Kilkenny,1980).,The,fluid,solution,assumes,that,there,are,too,many,energetic,particles,in,the,hot,plasma,regions,,and,predicts,too,few,energetic,particles,in,the,cooler,regions.,This,is,a,well,known,issue,,and,in,2D,fluid,simulations,,it,is,resolved,by,imposing,an,upper,bound,for,the,electron,heat,flux,that,is,of,order,neT,vte.,•,The,other,region,where,kinetic,effects,are,important,is,near,the,wall.,The,ions,and,electrons,that,reach,the,wall,recombine,and,do,not,come,back,as,charged,particles.,As,a,result,of,this,recombination,,the,charged,particle,distribution,functions,vanish,for,significant,parts,of,the,velocity,space,and,cannot,be,approximated,by,Maxwellians,,as,one,would,need,to,be,the,case,to,be,able,to,use,fluid,equations.,Thus,,the,treatment,of,charged,particles,near,the,wall,has,to,be,kinetic,(Loizu,et,al.,2011;,Geraldini,et,al.,2018).,The,combination,of,these,two,issues,can,manifest,in,a,population,of,energetic,electrons,that,cannot,be,predicted,by,fluid,models.,In,turn,,these,energetic,electrons,affect,the,density,and,temperature,by,suppressing,or,enhancing,heat,transport,and,by,changing,the,potential,difference,across,the,Debye,sheaths.,Evidence,of,these,effects,has,been,found,in,direct,kinetic,simulations,of,1D,problems,(Tskhakaya,et,al.,2011;,Chankin,&,Coster,2015).,The,fact,that,kinetic,effects,are,important,has,been,recognized,for,neutrals,,leading,to,kinetic,codes,for,them,(Reiter,et,al.,2005;,Stotler,&,Karney,1994).,However,,as,we,have,pointed,out,before,,these,codes,ignore,any,possible,non-Maxwellian,features,in,the,ion,and,electron,distribution,functions,,are,sensitive,to,the,parameters,that,determine,how,particles,interact,with,the,wall,and,,depending,on,the,version,,do,not,include,many,elastic,collisions,or,detailed,collision,physics.,The,non-Maxwellian,features,of,the,electron,dis-,tribution,function,are,particularly,important,as,electrons,mediate,many,of,the,processes,considered,important,in,these,neutral,models:,radiation,,ionization,,recombination,,etc.,For,ions,and,electrons,,the,community,is,starting,to,construct,edge,kinetic,codes,such,as,XGC,(Ku,et,al.,2016),,GKEYLL,(Hakim,et,al.,2020),or,COGENT,(Dorf,et,al.,2016).,These,codes,are,broadly,based,on,the,same,gyroaveraged,kinetic,models,used,in,δf,gyrokinetic,codes,such,as,GS2,(Kotschenreuther,et,al.,1995;,Dorland,et,al.,2000),,GENE,(Dannert,&,Jenko,2005),or,stella,(Barnes,et,al.,2019),,but,are,very,different,from,them,because,of,the,edge,particularities,,as,we,proceed,to,explain.,The,idea,behind,gyroaveraged,kinetic,models,is,to,average,over,the,fast,gyrofrequency,time,scale,to,avoid,a,cripplingly,small,time,step.,All,gyroaveraged,kinetic,models,are,based,on,an,asymptotic,expansion,in,ρi/L⊥,(cid:28),1,(see,table,2,for,values,of,ρi,and,L⊥).,The,simplest,possible,approach,is,drift,kinetics,(Hazeltine,1973),that,assumes,that,the,size,of,all,turbulent,structures,is,much,larger,than,ρi.,Unfortunately,,in,the,presence,of,temperature,and,density,gradients,,drift,kinetics,develops,instabilities,at,the,grid,scale.,These,instabilities,can,be,stabilized,by,numerical,dissipation,,and,in,the,real,word,,they,are,stabilized,by,finite,gyroradius,effects,that,drift,kinetics,neglects.,Gyrokinetics,(Catto,1978;,Frieman,&,Chen,1982),was,developed,to,solve,this,problem,in,the,core,of,the,tokamak.,Initially,,δf,gyrokinetics,assumed,that,turbulent,fluctuations,had,a,characteristic,size,of,the,order,of,ρi,and,their,amplitude,was,small,by,a,factor,of,ρi/L⊥,(cid:28),1.,This,ensures,that,the,gradients,of,the,fluctuations,are,comparable,to,the,background,gradient,and,not,larger.,However,,as,we,have,explained,in,section,3,,in,the,edge,it,is,not,possible,to,assume,that,density,and,temperature,are,slowly,varying,quantities,plus,small,fluctuations.,This,fact,8,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,was,part,of,the,justification,to,develop,what,is,know,as,full,f,gyrokinetics,,in,which,the,distribution,function,is,not,assumed,to,be,composed,of,a,slowly,varying,piece,plus,small,fluctuations.,See,Parra,&,Catto,(2008),and,Brizard,&,Hahm,(2007),for,two,different,types,of,derivations,of,full,f,gyrokinetics.,Overall,,full,f,gyrokinetics,can,be,seen,as,a,mixture,of,drift,kinetics,and,the,original,δf,gyrokinetics.,The,distribution,function,is,allowed,to,have,wavelengths,that,range,from,L⊥,to,ρi,,but,the,size,of,these,different,Fourier,components,has,to,be,sufficiently,small,that,the,gradient,of,the,distribution,function,is,never,larger,than,1/L⊥.,Assuming,a,core,ordering,(2.2),,one,naturally,recovers,δf,gyrokinetics,,showing,that,the,distribution,function,has,to,be,a,slowly,varying,piece,plus,small,fluctuations,(Parra,&,Catto,2010;,Calvo,&,Parra,2012).,In,the,edge,ordering,(2.1),,full,f,gyrokinetics,fundamentally,becomes,drift,kinetics,with,some,small,corrections.,These,corrections,are,of,two,types.,•,The,electromagnetic,fields,appear,in,the,kinetic,equation,averaged,over,circular,gyro-,orbits.,Thus,,electromagnetic,fluctuations,that,have,characteristic,lengths,much,smaller,than,the,gyroradius,are,averaged,over,and,do,not,drive,fluctuations,in,the,distribution,function.,Since,these,short,wavelength,electromagnetic,fluctuations,can,only,survive,if,there,are,charge,and,current,fluctuations,of,similar,wavelength,,which,they,are,not,able,to,drive,,they,eventually,damp,and,disappear.,•,The,field,equations,,quasineutrality,and,Amp´ere’s,law,,contain,densities,and,currents,that,one,calculates,from,the,ion,and,electron,distribution,functions.,As,a,result,of,the,gyrokinetic,expansion,,these,distribution,functions,have,finite,gyroradius,corrections,that,give,terms,that,are,small,in,ρi/L⊥,(cid:28),1.,One,of,these,small,terms,in,particular,,the,polarization,density,in,the,quasineutrality,equation,,is,important,because,,despite,its,small,size,,it,can,determine,the,electric,field,in,different,situations.,For,example,,for,shear,Alfven,waves,,the,polarization,density,is,small,,but,so,are,the,rest,of,the,contributions,to,the,density,,so,in,the,end,a,balance,between,the,polarization,density,and,another,term,determines,the,fluctuating,electric,field.,Another,example,is,the,component,of,the,electric,field,perpendicular,to,the,flux,surfaces,,which,is,also,determined,by,a,balance,between,the,polarization,density,and,other,terms,(Parra,&,Catto,2009).,The,small,finite,gyroradius,terms,in,the,full,f,gyrokinetic,equations,are,important,to,stabilize,the,short,wavelength,instabilities,,and,some,of,them,(e.g.,the,polarization,density),can,also,be,important,for,certain,aspects,of,the,physics.,Keeping,the,finite,gyro-,radius,effects,in,the,kinetic,equation,is,relatively,straightforward,in,PIC,codes,,although,one,has,to,be,careful,with,the,accuracy,of,the,average,(Guadagni,&,Cerfon,2017).,We,are,not,aware,of,any,full,f,edge,gyrokinetic,code,that,retains,finite,gyroradius,effects,in,the,kinetic,equation.,The,corrections,to,the,field,equations,are,much,more,difficult,to,retain,–,formulations,that,explicitly,try,to,conserve,energy,and,momentum,exactly,require,solving,nonlinear,equations,for,every,element,in,velocity,space,,for,instance.,This,has,driven,the,community,towards,simplifying,these,terms.,As,a,result,of,these,simplifications,,several,edge,codes,solve,drift,kinetics,with,some,ad,hoc,additions,to,the,field,equations,,such,as,a,simplified,polarization,density.,In,addition,to,these,fundamental,issues,,edge,drift,kinetics,and,gyrokinetics,need,to,address,other,problems.,As,gyrokinetics,was,devised,to,model,turbulent,fluctuations,with,a,spatial,size,of,the,order,of,ρi,in,tokamak,cores,,most,available,models,do,not,include,features,that,are,important,for,edge,physics.,•,There,is,no,gyrokinetic,formulation,for,the,wall,boundary,conditions.,Only,recently,one,such,formulation,was,developed,for,drift,kinetic,ions,in,magnetic,fields,that,reach,the,wall,with,a,grazing,angle,(Geraldini,et,al.,2018).,This,work,has,to,be,generalized,to,Physics,in,the,edge,of,fusion,devices,9,electrons,,more,general,angles,between,the,wall,and,the,magnetic,field,,and,eventually,to,gyrokinetics.,•,Gyrokinetics,does,not,usually,include,collisions,with,neutrals,,and,such,collisions,are,important,because,neutrals,break,the,symmetry,introduced,by,the,fast,charged,par-,ticle,gyration,around,magnetic,field,lines.,In,other,words,,the,dependence,on,gyrophase,,neglected,in,gyrokinetics,,can,become,important.,•,Gyrokinetics,was,developed,for,fluctuations,with,characteristic,scales,of,the,order,of,the,ion,gyroradius.,At,these,scales,,the,magnetic,field,fluctuations,have,a,very,constrained,form.,For,this,reason,,the,extension,of,gyrokinetics,to,include,MagnetoHydroDynamic,(MHD),modes,is,non-trivial,and,an,active,area,of,research.,See,Zheng,et,al.,(2007),for,a,theoretical,treatment,,and,Collar,et,al.,(2020),for,recent,numerical,work,in,this,area.,A,gyrokinetic,model,of,the,edge,should,be,able,to,reproduce,MHD,results,,as,MHD,modes,are,believed,to,be,the,main,drive,of,the,eruptions,known,as,Edge,Localized,Modes,(ELMs),(Ham,et,al.,2020).,5.,Discussion,A,complete,edge,description,requires,kinetic,effects.,Current,attempts,to,model,kinetic,effects,in,the,edge,rely,heavily,on,gyroaveraged,kinetic,models,that,average,over,the,very,fast,gyrofrequency,timescale.,In,our,opinion,,there,does,not,exist,a,systematic,procedure,to,choose,the,relevant,finite,gyroradius,effects,to,be,kept,in,the,kinetic,and,field,equations.,The,first,objective,of,a,kinetic,modeling,effort,for,the,edge,must,be,to,establish,the,finite,gyroradius,effects,to,be,kept,in,the,equations.,To,do,so,,for,contract,T/NA085/20,,we,have,proposed,as,a,first,attempt,to,construct,a,drift,kinetic,model.,This,drift,kinetic,model,will,be,unstable,at,grid,scales,,but,hyperviscosity,might,be,enough,to,stabilize,such,scales,if,they,do,not,contribute,much,to,transport,(as,one,would,expect,due,to,their,small,size).,For,the,contract,work,,we,will,use,the,model,only,in,1D,and,2D,configurations,that,cannot,develop,these,grid,scale,instabilities.,We,will,then,determine,analytically,which,finite,gyroradius,effects,must,be,kept,in,the,equations,to,determine,every,part,of,the,problem,,and,in,particular,the,component,of,the,electric,field,perpendicular,to,the,flux,surfaces.,A,well,known,issue,arises,when,keeping,finite,gyroradius,effects,in,the,field,equations.,In,quasineutrality,,the,only,term,that,contains,the,electric,potential,explicitly,is,the,small,finite,gyroradius,correction.,Thus,,unless,an,implicit,time,stepping,algorithm,is,employed,,one,needs,to,solve,for,the,potential,by,inverting,a,small,term,in,the,equation.,This,procedure,limits,the,time,step,size,severely,(Lee,1987;,Barnes,et,al.,2019).,Thus,,in,addition,to,keeping,finite,gyroradius,effects,,we,need,to,make,sure,that,the,kinetic,equations,that,we,obtain,are,amenable,to,implicit,time,stepping,methods.,In,addition,to,studying,finite,gyroradius,effects,,we,will,determine,the,effect,that,collisions,with,neutrals,have,on,the,gyrokinetic,formalism,by,introducing,charge,exchange,collisions,and,ionization,collisions.,The,final,deliverable,of,contract,T/NA085/20,will,be,a,drift,kinetic,model,with,neutrals,and,the,finite,gyroradius,terms,that,are,needed,to,calculate,the,component,of,the,electric,field,that,is,perpendicular,to,the,flux,surfaces.,Both,of,these,features,will,be,improvements,on,the,models,implemented,in,existing,continuum,edge,codes.,By,the,end,of,the,contract,,the,model,will,have,been,tested,in,1D,and,2D,problems,,and,even,though,it,will,have,a,3D,version,,this,3D,version,will,not,have,been,tested,in,the,turbulent,regime.,Moreover,,all,the,work,will,have,been,performed,in,a,helical,field,and,not,in,a,diverted,tokamak,field,because,the,helical,field,is,the,state-of-the-art,for,continuum,edge,kinetic,codes,(only,recently,,in,the,2020,Annual,APS,DPP,meeting,,the,10,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,two,leading,continuum,codes,,GKEYLL,and,COGENT,,have,reported,the,first,runs,with,more,complex,tokamak,geometry).,The,limitations,of,the,model,delivered,at,the,end,of,contract,T/NA085/20,leave,two,obvious,avenues,for,future,work.,•,The,main,theoretical,difficulty,added,by,tokamak,magnetic,fields,is,the,presence,of,a,component,of,the,drifts,in,the,direction,perpendicular,to,the,flux,surface,that,does,not,exist,in,the,case,of,a,helical,field.,This,component,of,the,drifts,leads,to,finite,orbit,widths,and,can,significantly,modify,some,of,the,physics,of,both,quiescent,2D,plasmas,(Kagan,&,Catto,2008;,Landreman,et,al.,2014),and,3D,turbulence,(Parisi,et,al.,2020).,Adding,this,component,of,the,particle,drift,to,the,equations,will,require,extra,work.,•,As,we,have,explained,above,,the,3D,model,is,expected,to,have,grid,scale,instabilities,that,can,be,stabilized,by,hyperviscosity,or,other,numerical,methods,of,dissipation.,To,properly,capture,the,fluctuations,at,the,ion,gyroradius,scale,,one,would,have,to,include,more,precise,finite,gyroradius,effects.,This,is,an,obvious,extension,of,the,work,in,contract,T/NA085/20,,and,it,is,very,important,for,High,confinement,mode,(H-mode),where,the,turbulence,at,scales,larger,than,the,ion,gyroradius,is,stabilized,for,still,unclear,reasons.,In,this,regime,,turbulence,at,scales,of,the,order,of,or,smaller,than,the,ion,gyroradius,is,important,(Hillesheim,et,al.,2016;,Hatch,et,al.,2017;,Parisi,et,al.,2020),,and,the,finite,gyroradius,effects,become,crucial,for,ions.,We,foresee,that,adding,more,detailed,finite,gyroradius,effects,will,be,theoretically,and,numerically,challenging,and,will,require,dedicated,work.,Two,other,aspects,of,an,edge,kinetic,model,are,beyond,the,scope,of,contract,T/NA085/20.,•,To,test,the,effect,of,wall,boundary,conditions,on,drift,kinetics,,we,will,impose,a,simplified,version,of,the,boundary,conditions,that,are,valid,in,the,limit,in,which,the,electron,gyroradius,is,much,smaller,than,the,Debye,length.,However,,more,detailed,boundary,conditions,must,be,found,because,usually,the,electron,gyroradius,is,larger,than,or,comparable,to,the,Debye,length.,This,is,work,that,is,being,pursued,by,one,of,the,authors,of,this,report,(F.I.P.),with,other,sources,of,funding.,•,We,pointed,out,above,that,it,would,be,desirable,to,be,able,to,recover,MHD,modes,with,the,edge,kinetic,model.,In,the,work,for,contract,T/NA085/20,,the,fluctuations,in,the,magnetic,field,will,be,neglected,(this,is,a,good,approximation,for,many,edge,plasmas),,and,hence,it,will,not,be,possible,to,explore,connections,with,MHD.,This,is,an,area,of,research,where,ExCALIBUR,NEPTUNE,could,benefit,from,collaboration,with,the,EPSRC,Programme,Grant,‘Turbulent,Dynamics,of,Tokamak,Plasmas,(TDoTP)’.,Several,PIs,in,the,ExCALIBUR,NEPTUNE,project,also,belong,to,TDoTP.,REFERENCES,Barnes,,M.,,Parra,,F.I.,&,Landreman,,M.,2019,stella:,An,operator-split,,implicit-explicit,δf,-gyrokinetic,code,for,general,magnetic,field,configurations.,J.,Comput.,Phys.,391,,365.,Beer,,M.A.,,Cowley,,S.C.,&,Hammett,,G.W.,1995,Field-aligned,coordinates,for,nonlinear,simulations,of,tokamak,turbulence.,Phys.,Plasmas,2,,2687.,Brackmann,,R.T.,,Fite,,W.L.,&,Neynaber,,R.H.,1958,Collisions,of,Electrons,with,Hydrogen,Atoms.,III.,Elastic,Scattering.,Phys.,Rev.,112,,1157.,Braginskii,,S.I.,1958,Transport,phenomena,in,a,completely,ionized,two-temperature,plasma.,Sov.,Phys.,JETP,6,,358.,Brizard,,A.J.,&,Hahm,,T.S.,2007,Foundations,of,nonlinear,gyrokinetic,theory.,Rev.,Mod.,Phys.,79,,421.,Bufferand,,H.,,Ciraolo,,G.,,Marandet,,Y.,,Bucalossi,,J.,,Ghendrih,,Ph.,,Gunn,,J.,,Mellet,,N.,,Tamain,,P.,,Leybros,,R.,,Fedorczak,,N.,,Schwander,,F.,&,Serre,,E.,2015,Numerical,modelling,for,divertor,design,of,the,WEST,device,with,a,focus,on,plasmawall,interactions.,Nucl.,Fusion,55,,053025.,Bufferand,,H.,,Tamain,,P.,,Baschetti,,S.,,Bucalossi,,J.,,Ciraolo,,G.,,Fedorczak,,N.,,Physics,in,the,edge,of,fusion,devices,11,Ghendrih,,Ph.,,Nespoli,,F.,,Schwander,,F.,&,E.,Serre,and,,Y.,Marandet,2019,Three-dimensional,modelling,of,edge,multi-component,plasma,taking,into,account,realistic,wall,geometry.,Nucl.,Mat.,Energy,18,,82.,Calvo,,I.,&,Parra,,F.I.,2012,Long-wavelength,limit,of,gyrokinetics,in,a,turbulent,tokamak,and,its,intrinsic,ambipolarity.,Plasma,Phys.,Control.,Fusion,54,,115007.,Catto,,P.J.,1978,Linearized,gyro-kinetics.,Plasma,Phys.,20,,719.,Catto,,P.J.,1994,A,short,mean-free,path,,coupled,neutral-ion,transport,description,of,a,toka-,mak,edge,plasma.,Phys.,Plasmas,1,,1936.,Catto,,P.J.,&,Simakov,,A.N.,2004,A,drift,ordered,short,mean,free,path,description,for,magnetized,plasma,allowing,strong,spatial,anisotropy.,Phys.,Plasmas,11,,90.,Chankin,,A.V.,,Corrigan,,G.,&,JET,Contributors,2021,EDGE2D-EIRENE,modeling,of,the,impact,of,wall,materials,on,core,edge,,scrape-off,layer,and,divertor,parameters.,Plasma,Phys.,Control.,Fusion,63,,035010.,Chankin,,A.V.,&,Coster,,D.P.,2015,On,the,locality,of,parallel,transport,of,heat,carrying,electrons,in,the,SOL.,J.,Nucl.,Mater.,463,,498.,Chodura,,R.,1982,Plasma-wall,transition,in,an,oblique,magnetic,field.,Phys.,Fluids,25,,1628.,Colchin,,R.J.,,Maingi,,R.,,Fenstermacher,,M.E.,,Carlstrom,,T.N.,,Isler,,R.C.,,Owen,,L.W.,&,Groebner,,R.J.,2000,Measurement,of,neutral,density,near,the,X,point,in,the,DIII-D,tokamak.,Nucl.,Fusion,40,,175.,Collar,,J.P.M.,,McMillan,,B.F.,,Saarelma,,S.,&,Bottino,,A.,2020,Comparing,electro-,magnetic,instabilities,in,global,gyrokinetic,simulations,with,local,and,MHD,models.,Plasma,Phys.,Control.,Fusion,62,,095005.,Dannert,,T.,&,Jenko,,F.,2005,Gyrokinetic,simulation,of,collisionless,trapped-electron,mode,turbulence.,Phys.,Plasmas,12,,072309.,Dorf,,M.A.,,Dorr,,M.R.,,Hittinger,,J.A.,,Cohen,,R.H.,&,Rognlien,,T.D.,2016,Contin-,uum,kinetic,modeling,of,the,tokamak,plasma,edge.,Phys.,Plasmas,23,,056102.,Dorland,,W.,,Jenko,,F.,,Kotschenreuther,,M.,&,Rogers,,B.N.,2000,Electron,Temper-,ature,Gradient,Turbulence.,Phys.,Rev.,Lett.,85,,5579.,Dudson,,B.D.,&,Leddy,,J.,2017,Hermes:,global,plasma,edge,fluid,turbulence,simulations.,Plasma,Phys.,Control.,Fusion,59,,054010.,Frieman,,E.A.,&,Chen,,L.,1982,Nonlinear,gyrokinetic,equations,for,low-frequency,electro-,magnetic,waves,in,general,plasma,equilibria.,Phys.,Fluids,25,,502.,Geraldini,,A.,,Parra,,F.I.,&,Militello,,F.,2018,Solution,to,a,collisionless,shallow-angle,magnetic,presheath,with,kinetic,ions.,Plasma,Phys.,Control.,Fusion,60,,125002.,Goldston,,R.J.,2012,Heuristic,drift-based,model,of,the,power,scrape-off,width,in,low-gas-puff,H-mode,tokamaks.,Nucl.,Fusion,52,,013009.,Gray,,D.R.,&,Kilkenny,,J.D.,1980,The,measurement,of,ion,acoustic,turbulence,and,reduced,thermal,conductivity,caused,by,a,large,temperature,gradient,in,a,laser,heated,plasma.,Plasma,Phys.,22,,81.,Guadagni,,J.,&,Cerfon,,A.J.,2017,Fast,and,spectrally,accurate,evaluation,of,gyroaverages,in,non-periodic,gyrokinetic,Poisson,simulations.,J.,Plasma,Phys.,83,,905830407.,Gurevich,,A.,&,Istomin,,Ya.N.,1979,Thermal,runaway,and,convective,heat,transport,by,fast,electrons,in,a,plasma.,Sov.,Phys.,JETP,50,,470.,Hakim,,A.H.,,Mandell,,N.R.,,T.N.,Bernard,and,,M.,Francisquez,,Hammett,,G.W.,&,Shi,,E.,L.,2020,Continuum,electromagnetic,gyrokinetic,simulations,of,turbulence,in,the,tokamak,scrape-off,layer,and,laboratory,devices.,Phys.,Plasmas,27,,042304.,Halpern,,F.D.,,Ricci,,P.,,Jolliet,,S.,,Loizu,,J.,,Morales,,J.,,Mosetto,,A.,,Musil,,F.,,Riva,,F.,,Tran,,T.M.,&,Wersal,,C.,2016,The,GBS,code,for,tokamak,scrape-off,layer,simulations.,J.,Comput.,Phys.,315,,388.,Ham,,C.,,Kirk,,A.,,Pamela,,S.,&,Wilson,,H.,2020,Filamentary,plasma,eruptions,and,their,control,on,the,route,to,fusion,energy.,Nature,Rev.,Phys.,2,,159.,Hariri,,F.,&,Ottaviani,,M.,2013,A,flux-coordinate,independent,field-aligned,approach,to,plasma,turbulence,simulations.,Comput.,Phys.,Comm.,184,,2419.,Hastie,,R.J.,1997,Sawtooth,instability,in,tokamak,plasmas.,Astrophys.,Space,Sci.,256,,177.,Hatch,,D.R.,,Kotschenreuther,,M.,,Mahajan,,S.,,Valanju,,P.,&,Liu,,X.,2017,A,gyroki-,netic,perspective,on,the,JET-ILW,pedestal.,Nucl.,Fusion,57,,036020.,Hazeltine,,R.D.,1973,Recursive,derivation,of,drift-kinetic,equation.,Plasma,Phys.,15,,77–80.,12,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Hazeltine,,R.D.,,Calvin,,M.D.,,Valanju,,P.M.,&,Solano,,E.R.,1992,Analytical,calculation,of,neutral,transport,and,its,effect,on,ions.,Nucl.,Fusion,32,,3.,Helander,,P.,,Krasheninnikov,,S.I.,&,Catto,,P.J.,1994,Fluid,equations,for,a,partially,ionized,plasma.,Phys.,Plasmas,1,,3174.,Hillesheim,,J.C.,,Dickinson,,D.,,Roach,,C.M.,,Saarelma,,S.,,Scannell,,R.,,Kirk,,A.,,Crocker,,N.A.,,Peebles,,W.A.,,Meyer,,H.,&,the,MAST,Team,2016,Intermediate-,k,density,and,magnetic,field,fluctuations,during,inter-ELM,pedestal,evolution,in,MAST.,Plasma,Phys.,Control.,Fusion,58,,014020.,Isoardi,,L.,,Chiavassa,,G.,,Ciraolo,,G.,,Haldenwang,,P.,,Serre,,E.,,Ghendrih,,Ph.,,Sarazin,,Y.,,Schwander,,F.,&,Tamain,,P.,2010,Penalization,modeling,of,a,limiter,in,the,Tokamak,edge,plasma.,J.,Comput.,Phys.,229,,2220.,Jackson,,G.L.,,Winter,,J.,,Taylor,,T.S.,,Burrell,,K.H.,,DeBoo,,J.C.,,Greenfield,,C.M.,,Groebner,,R.J.,,Hodapp,,T.,,Holtrop,,K.,,Lazarus,,E.A.,,Lao,,L.L.,,Lipp-,mann,,S.I.,,Osborne,,T.H.,,Petrie,,T.W.,,Phillips,,J.,,James,,R.,,Schissel,,D.P.,,Strait,,E.J.,,Turnbull,,A.D.,,West,,W.P.,&,the,DIII-D,Team,1991,Regime,of,very,high,confinement,in,the,boronized,DIII-D,tokamak.,Phys.,Rev.,Lett.,67,,3098.,Kagan,,G.,&,Catto,,P.J.,2008,Arbitrary,poloidal,gyroradius,effects,in,tokamak,pedestals,and,transport,barriers.,Plasma,Phys.,Control.,Fusion,50,,085010.,Kotschenreuther,,M.,,Rewoldt,,G.,&,Tang,,W.M.,1995,Comparison,of,initial,value,and,eigenvalue,codes,for,kinetic,toroidal,plasma,instabilities.,Comput.,Phys.,Comm.,88,,128.,Krasheninnikov,,S.I.,&,Kukushkin,,A.S.,2017,Physics,of,ultimate,detachment,of,a,tokamak,divertor,plasma.,J.,Plasma,Phys.,83,,155830501.,Ku,,S.,,Hager,,R.,,Chang,,C.S.,,Kwon,,J.M.,&,Parker,,S.E.,2016,A,new,hybrid-Lagrangian,numerical,scheme,for,gyrokinetic,simulation,of,tokamak,edge,plasma.,J.,Comput.,Phys.,315,,467.,Landreman,,M.,,Parra,,F.I.,,Catto,,P.J.,,Ernst,,D.R.,&,Pusztai,,I.,2014,Radially,global,δf,computation,of,neoclassical,phenomena,in,a,tokamak,pedestal.,Plasma,Phys.,Control.,Fusion,56,,045005.,Lee,,W.W.,1987,Gyrokinetic,particle,simulation,model.,J.,Comput.,Phys.,72,,243.,Lindsay,,B.G.,&,Stebbings,,R.F.,2005,Charge,transfer,cross,sections,for,energetic,neutral,atom,data,analysis.,J.,Geophys.,Res.,110,,A12213.,Loizu,,J.,,Ricci,,P.,&,Theiler,,C.,2011,Existence,of,subsonic,plasma,sheaths.,Phys.,Rev.,E,83,,016406.,Mikhailovskii,,A.B.,&,Tsypin,,V.S.,1971,Transport,equations,and,gradient,instabilities,in,a,high,pressure,collisional,plasma.,Plasma,Phys.,13,,785.,Militello,,F.,&,Fundamenski,,W.,2011,Multi-machine,comparison,of,drift,fluid,dimensionless,parameters.,Plasma,Phys.,Control.,Fusion,53,,095002.,Parisi,,J.F.,,Parra,,F.I.,,Roach,,C.M.,,Giroud,,C.,,Dorland,,W.,,Hatch,,D.R.,,Barnes,,M.,,Hillesheim,,J.C.,,Aiba,,N.,&,Ball,,J.,2020,Toroidal,and,slab,ETG,instability,dominance,in,the,linear,spectrum,of,JET-ILW,pedestals.,Nucl.,Fusion,60,,126045.,Parra,,F.I.,&,Catto,,P.J.,2008,Limitations,of,gyrokinetics,on,transport,time,scales.,Plasma,Phys.,Control.,Fusion,50,,065014.,Parra,,F.I.,&,Catto,,P.J.,2009,Vorticity,and,intrinsic,ambipolarity,in,turbulent,tokamaks.,Plasma,Phys.,Control.,Fusion,51,,095008.,Parra,,F.I.,&,Catto,,P.J.,2010,Turbulent,transport,of,toroidal,angular,momentum,in,low,flow,gyrokinetics.,Plasma,Phys.,Control.,Fusion,52,,045004.,Reiter,,D.,,Baelmans,,M.,&,B¨orner,,P.,2005,The,EIRENE,and,B2-EIRENE,Codes.,Fusion,Sci.,Tech.,47,,172.,Riemann,,K.-U.,1991,The,Bohm,criterion,and,sheath,formation.,J.,Phys.,D:,Appl.,Phys.,24,,493.,Rognlien,,T.D.,,Rensink,,M.E.,&,Smith,,G.R.,2007,User,manual,for,the,UEDGE,edge-,plasma,transport,code.,Tech.,Rep.,UCRL-ID-137121.,Lawrence,Livermore,National,Labo-,ratory.,Schmitt,,J.C.,,Bell,,R.E.,,Boyle,,D.P.,,Esposti,,B.,,Kaita,,R.,,Kozub,,T.,,LeBlanc,,B.P.,,Lucia,,M.,,Maingi,,R.,,Majeski,,R.,,Merino,,E.,,Punjabi-Vinoth,,S.,,Tchilin-,gurian,,G.,,Capece,,A.,,Koel,,B.,,Roszell,,J.,,Biewer,,T.,M.,,Gray,,T.,K.,,Kubota,,Physics,in,the,edge,of,fusion,devices,13,S.,,Beiersdorfer,,P.,,Widmann,,K.,&,Tritz,,K.,2015,High,performance,discharges,in,the,Lithium,Tokamak,eXperiment,with,liquid,Lithium,walls.,Phys.,Plasmas,22,,056112.,Scotti,,F.,,Stotler,,D.P.,,Bell,,R.E.,,LeBlanc,,B.P.,,Sabbagh,,S.A.,,Soukhanovskii,,V.A.,,Umansky,,M.V.,&,Zweben,,S.J.,2021,Outer,midplane,neutral,density,measure-,ments,and,H-mode,fueling,studies,in,NSTX-U.,Nucl.,Fusion,61,,036002.,Simakov,,A.N.,&,Catto,,P.J.,2003,Drift-ordered,fluid,equations,for,field-aligned,modes,in,low-β,collisional,plasma,with,equilibrium,pressure,pedestals.,Phys.,Plasmas,10,,4744.,Simonini,,R.,,Corrigan,,G.,,Radford,,G.,,Spence,,J.,&,Taroni,,A.,1994,Models,and,Numerics,in,the,Multi-Fluid,2-D,Edge,Plasma,Code,EDGE2D/U.,Contrib.,Plasma,Phys.,34,,368.,Stegmeir,,A.,,Coster,,D.,,Ross,,A.,,Maj,,O.,,Lackner,,K.,&,Poli,,E.,2018,GRILLIX:,a,3D,turbulence,code,based,on,the,flux-coordinate,independent,approach.,Plasma,Phys.,Control.,Fusion,60,,035005.,Stotler,,D.,&,Karney,,C.,1994,Neutral,Gas,Transport,Modeling,with,DEGAS,2.,Contrib.,Plasma,Phys.,34,,392.,Sugihara,,M.,,Igitkhanov,,Yu.,,Janeschitz,,G.,,Hubbard,,A.E.,,Kamada,,Y.,,Linger-,tat,,J.,,Osborne,,T.H.,&,Suttrop,,W.,2000,A,model,for,H,mode,pedestal,width,scaling,using,the,International,Pedestal,Database.,Nucl.,Fusion,40,,1743.,Tamain,,P.,,Bufferand,,H.,,Ciraolo,,G.,,Colin,,C.,,Galassi,,D.,,Ph.Ghendrih,,Schwan-,der,,F.,&,Serre,,E.,2016,The,TOKAM3X,code,for,edge,turbulence,fluid,simulations,of,tokamak,plasmas,in,versatile,magnetic,geometries.,J.,Comput.,Phys.,321,,606.,Tskhakaya,,D.,,Jachmich,,S.,,Eich,,T.,,Fundamenski,,W.,&,JET,EFDA,contributors,2011,Interpretation,of,divertor,Langmuir,probe,measurements,during,the,ELMs,at,JET.,J.,Nucl.,Mater.,415,,S860.,Wersal,,C.,&,Ricci,,P.,2015,A,first-principles,self-consistent,model,of,plasma,turbulence,and,kinetic,neutral,dynamics,in,the,tokamak,scrape-off,layer.,Nucl.,Fusion,55,,123014.,Wiesen,,S.,,Brezinsek,,S.,,Bonnin,,X.,,Delabie,,E.,,Frassinetti,,L.,,Groth,,M.,,Guille-,maut,,C.,,Harrison,,J.,,Harting,,D.,,Henderson,,S.,,Huber,,A.,,Kruezi,,U.,,Pitts,,R.A.,,Wischmeier,,M.,&,JET,contributors,2018,On,the,role,of,finite,grid,extent,in,SOLPS-ITER,edge,plasma,simulations,for,JET,H-mode,discharges,with,metallic,wall.,Nucl.,Mat.,Energy,17,,174.,Wiesen,,S.,,Reiter,,D.,,Kotov,,V.,,Baelmans,,M.,,Dekeyser,,W.,,Kukushkin,,A.S.,,Lisgo,,S.W.,,Pitts,,R.A.,,Rozhansky,,V.,,Saibene,,G.,,Veselova,,I.,&,Voskoboynikov,,S.,2015,The,new,SOLPS-ITER,code,package.,J.,Nucl.,Mater.,463,,480.,Zeiler,,A.,,Biskamp,,D.,,Drake,,J.F.,&,Guzdar,,P.N.,1996,Three-dimensional,fluid,simu-,lations,of,tokamak,edge,turbulence.,Phys.,Plasmas,3,,2951.,Zel’dovich,,Ya.,B.,&,Raizer,,Yu.,P.,2013,Physics,of,Shock,Waves,and,High-Temperature,Hydrodynamic,Phenomena.,Dover,Publications.,Zheng,,L.J.,,Kotschenreuther,,M.T.,&,Dam,,J.W.,Van,2007,Revisiting,linear,gyrokinetics,to,recover,ideal,magnetohydrodynamics,and,missing,finite,Larmor,radius,effects.,Phys.,Plasmas,14,,072505.,REPORT,2047357-TN-01-02,M1.1,1,1D,drift,kinetic,models,with,periodic,boundary,conditions,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Rudolf,Peierls,Centre,for,Theoretical,Physics,,University,of,Oxford,,Oxford,OX1,3PU,,UK,(This,version,is,of,22,February,2021),1.,Introduction,In,this,report,,we,propose,1D,drift,kinetic,equations,to,test,the,possibility,of,extracting,low,order,moments,from,the,distribution,functions,for,implicit,methods.,The,model,that,we,present,here,has,periodic,boundary,conditions,,adequate,for,the,closed,field,line,region,of,the,edge.,We,will,address,wall,boundary,conditions,for,open,field,lines,in,the,reports,for,milestones,M1.3,,M2.4,and,M2.5.,2.,1D,electrostatic,drift,kinetics,We,consider,a,plasma,with,one,ion,species,with,charge,e,and,mass,mi,,electrons,with,charge,e,and,mass,me,,and,one,species,of,neutrals,with,mass,−,mn,=,mi.,(2.1),The,plasma,is,magnetized,by,a,constant,magnetic,field,B,=,Bˆz,,and,we,assume,that,the,plasma,only,varies,along,magnetic,field,lines.,In,this,case,,the,electric,field,produced,by,the,plasma,is,electrostatic,,E,=,(∂φ/∂z)ˆz.,The,potential,φ(z,,t),depends,on,the,position,along,magnetic,field,lines,z,and,on,time,t.,−,If,we,assume,that,the,gyroradii,are,small,compared,to,the,length,scales,of,interest,,and,that,the,gyrofrequencies,are,much,larger,than,the,frequencies,that,we,want,to,model,,,t),of,the,different,species,s,=,(Hazeltine,1973),,the,distribution,functions,fs(z,,v,(cid:107),i,,e,,n,only,depend,on,the,component,of,the,velocity,parallel,to,the,magnetic,field,v,and,the,magnitude,of,the,velocity,perpendicular,to,the,magnetic,field,v,,,and,are,independent,of,the,direction,of,the,velocity,perpendicular,to,the,magnetic,field.,Thus,,the,distribution,functions,that,in,general,can,depend,on,three,spatial,variables,r,,three,components,of,the,velocity,v,and,the,time,t,depend,only,on,z,,v,and,t,,,,v,,,v,⊥,⊥,(cid:107),fs(r,,v,,t),=,fs(z,,v,(cid:107),(cid:107),,,v,⊥,,,t).,⊥,(2.2),(2.3),The,equations,for,the,distribution,functions,of,the,different,species,are,∂fi,∂t,+,v,∂fi,∂z,−,e,mi,∂φ,∂z,(cid:107),∂fi,∂v,(cid:107),=,Cii[fi],+,Cin[fi,,fn],,∂fe,∂t,+,v,∂fe,∂z,(cid:107),+,e,me,∂φ,∂z,∂fe,∂v,(cid:107),and,=,Cee[fe],+,Cei[fe,,fi],+,Cen[fe,,fn],(2.4),∂fn,∂t,+,v,∂fn,∂z,(cid:107),=,Cni[fn,,fi].,(2.5),Here,we,have,included,ion-ion,and,electron-electron,collisions,,modeled,by,the,Fokker-,Planck,collision,operators,Cii[fi],and,Cee[fe],(Rosenbluth,et,al.,1957),,elastic,electron-ion,2,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,and,electron-neutral,collisions,,modeled,by,the,simplified,Fokker-Planck,collision,operator,Cei[fe,,fi],(Braginskii,1958),and,the,Boltzmann,collision,operator,Cen[fe,,fn],,and,charge-,exchange,collisions,,represented,by,the,simplified,Boltzmann,collision,operators,and,Cin[fi,,fn],=,Cni[fn,,fi],=,(cid:90),(cid:90),−,−,v,Rin(,|,v(cid:48),),[fi(v)fn(v(cid:48)),|,−,−,fi(v(cid:48))fn(v)],d3v(cid:48),(2.6),Rin(,v,|,v(cid:48),|,−,),[fn(v)fi(v(cid:48)),−,fn(v(cid:48))fi(v)],d3v(cid:48),,(2.7),To,simplify,our,equations,,we,assume,that,the,function,Rin,is,constant,(Connor,1977;,Hazeltine,et,al.,1992;,Catto,1994),,finding,and,where,the,densities,are,Cin[fi,,fn],=,Rin,(nnfi,−,−,nifn),Cni[fn,,fi],=,Rin,(nifn,−,−,nnfi),,,ns(z,,t),:=,2π,(cid:90),∞,−∞,dv,(cid:90),∞,(cid:107),0,dv,v,fs(z,,v,(cid:107),⊥,⊥,,,v,⊥,,,t).,(2.8),(2.9),(2.10),Note,that,we,can,neglect,the,effect,of,electron,collisions,on,ions,and,on,neutrals,due,to,the,smallness,of,the,electron,mass,(Braginskii,1958).,We,have,also,neglected,neutral-neutral,collisions,because,,in,current,fusion,devices,,the,neutral,density,is,sufficiently,small,that,the,neutral-neutral,collisions,are,rare.,The,kinetic,equations,will,be,solved,in,the,interval,z,periodic,boundary,conditions,at,z,=,0,and,z,=,L,,∈,[0,,L],,and,we,will,impose,fs(z,=,0,,v,(cid:107),,,v,⊥,,,t),=,fs(z,=,L,,v,(cid:107),,,v,⊥,,,t).,Finally,,the,potential,φ(z,,t),is,determined,by,the,quasineutrality,equation,ni,=,ne.,(2.11),(2.12),To,solve,this,equation,,we,need,to,treat,the,equations,implicitly,as,the,potential,enters,only,via,its,effect,on,∂fi/∂t,and,∂fe/∂t.,The,need,to,use,implicit,methods,is,one,of,the,reasons,why,we,are,trying,to,extract,some,of,the,low,order,moments,from,the,distribution,function,,notably,the,density.,Before,we,treat,the,complete,problem,,we,will,simplify,the,treatment,of,electrons,to,obtain,a,system,of,equations,that,can,be,solved,with,an,explicit,time,advance,so,that,we,can,compare,our,implicit,schemes,with,an,explicit,numerical,method.,Instead,of,solving,for,fe,,we,will,use,a,Maxwell-Boltzmann,response,,ne(z,,t),=,Ne,exp,(cid:18),eφ(z,,t),Te,(cid:19),,,(2.13),where,Ne,and,Te,are,constants,(see,Appendix,A,for,a,derivation,of,the,Maxwell-Boltzmann,response).,Moreover,,the,full,Fokker-Planck,ion-ion,collision,operator,Cii[fi],is,a,com-,plicated,integro-differential,operator,that,we,will,not,implement,in,the,first,versions,of,our,drift,kinetic,code,,so,we,do,not,include,it,in,the,equations,for,now.,Thus,,the,final,,,v,simplified,model,for,fi(z,,v,(cid:107),,,t),and,φ(z,,t),is,given,by,the,equations,,,v,⊥,∂fi,∂t,+,v,(cid:107),⊥,∂fi,∂z,−,,,t),,fn(z,,v,(cid:107),∂fi,∂v,(cid:107),e,mi,∂φ,∂z,=,Rin(nnfi,−,−,nifn),,(2.14),1D,drift,kinetic,models,with,periodic,boundary,conditions,∂fn,∂t,+,v,∂fn,∂z,(cid:107),=,Rin(nifn,−,−,nnfi),and,ni,=,Ne,exp,(cid:19),,,(cid:18),eφ,Te,3,(2.15),(2.16),with,periodic,boundary,conditions,(2.11).,This,system,of,equations,can,be,solved,explic-,itly,because,the,simple,electron,model,allows,one,to,obtain,φ,as,a,function,of,ni.,3.,1D,moment,drift,kinetics,Instead,of,solving,for,fs(z,,v,(cid:107),v3,ts(z,,t),ns(z,,t),Fs(z,,w,(cid:107),,,t),:=,,,w,,,v,⊥,⊥,,,t),,we,solve,for,(cid:16),fs,z,,us,(cid:107),(z,,t),+,vts(z,,t)w,(cid:107),,,vts(z,,t)w,(cid:17),,,t,,,⊥,where,we,have,defined,the,normalized,velocities,w,(z,,v,(cid:107),(cid:107),,,t),:=,v,(z,,t),(cid:107),−,us,(cid:107),vts(z,,t),and,the,average,parallel,velocity,w,⊥,(z,,v,⊥,,,t),:=,v,⊥,vts(z,,t),,,(z,,t),:=,us,(cid:107),2π,ns,(cid:90),∞,−∞,dv,(cid:90),∞,(cid:107),0,dv,v,v,fs(z,,v,(cid:107),(cid:107),⊥,⊥,,,v,⊥,,,t),and,the,thermal,speed,vts(z,,t),:=,(cid:115),4π,3ns,(cid:90),∞,−∞,dv,(cid:90),∞,(cid:107),0,dv,[(v,v,⊥,⊥,(cid:107),−,us,(cid:107),(z,,t))2,+,v2,⊥,]fs(z,,v,(cid:107),,,v,⊥,,,t).,(3.5),According,to,its,definition,,Fs(z,,w,(cid:107),(cid:90),∞,(cid:90),∞,,,w,⊥,,,t),must,satisfy,the,conditions,2π,dw,(cid:107),0,dw,⊥,w,⊥,Fs(z,,w,(cid:107),,,w,⊥,,,t),=,1,,−∞,(cid:90),∞,dw,2π,−∞,(cid:90),∞,(cid:107),0,dw,w,⊥,w,Fs(z,,w,(cid:107),(cid:107),⊥,,,w,⊥,,,t),=,0,(cid:90),∞,2π,(cid:90),∞,0,dw,(cid:107),dw,⊥,w,⊥,(w2,(cid:107),+,w2,⊥,)Fs(z,,w,(cid:107),,,w,⊥,,,t),=,3,2,and,at,every,point,z,and,time,t.,−∞,The,equations,for,ions,become,∂ni,∂t,+,∂,∂z,(cid:0)niui,(cid:107),(cid:1),=,0,,nimi,(cid:18),∂ui,∂t,(cid:107),+,ui,(cid:19),∂ui,(cid:107),∂z,(cid:107),=,∂pi,(cid:107),∂z,−,−,eni,∂φ,∂z,+,ninnmiRin(un,),,ui,(cid:107),(cid:107),−,(3.10),(3.1),(3.2),(3.3),(3.4),(3.6),(3.7),(3.8),(3.9),4,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,3,2,nimivti,(cid:18),∂vti,∂t,+,ui,(cid:107),∂vti,∂z,(cid:19),=,∂qi,(cid:107),∂z,−,∂ui,(cid:107),∂z,pi,(cid:107),−,+,3,4,ninnmiRin(v2,tn,−,v2,ti),+,1,2,ninnmiRin(un,)2,ui,(cid:107),(cid:107),−,and,∂Fi,∂t,+,˙zi,∂Fi,∂z,+,˙w,i,(cid:107),Here,,we,have,defined,the,coefficients,∂Fi,∂w,(cid:107),+,˙w,⊥,i,∂Fi,∂w,⊥,=,˙Fi,+,Cin.,˙zs[Fs](z,,w,(cid:107),,,t),:=,us,(cid:107),+,vtsw,,,(cid:107),(3.11),(3.12),(3.13),˙w,s[Fs](z,,w,(cid:107),(cid:107),,,t),:=,1,nsmsvts,∂ps,(cid:107),∂z,+,2w,(cid:107),3nsmsv2,ts,(cid:20),∂qs,(cid:107),∂z,(cid:18),+,ps,(cid:107),−,3,2,nsmsv2,ts,(cid:21),(cid:19),∂us,∂z,(cid:107),∂vts,∂z,,,w2,(cid:107),−,˙w,⊥,s[Fs](z,,w,(cid:107),,,w,⊥,,,t),:=,2w,3nsmsv2,ts,⊥,(cid:18),∂qs,(cid:107),∂z,+,ps,(cid:107),∂us,(cid:107),∂z,(cid:19),w,w,(cid:107),⊥,−,∂vts,∂z,(3.14),(3.15),and,˙Fs[Fs](z,,w,(cid:107),,,w,⊥,,,t),:=,w,(cid:107),(cid:34),(cid:18),3,2,nsmsv2,ts,−,the,parallel,pressure,∂vts,∂z,−,vts,ns,(cid:18),∂qs,(cid:107),∂z,(cid:19),∂ns,∂z,(cid:18),+,ps,(cid:107),−,1,2,nsmsv2,ts,(cid:19),(cid:35),(cid:19),∂us,∂z,(cid:107),Fs,,(3.16),[Fs](z,,t),:=,2πnsmsv2,ts,ps,(cid:107),(cid:90),∞,−∞,(cid:90),∞,0,dw,(cid:107),dw,⊥,w,⊥,w2,(cid:107),Fs(z,,w,(cid:107),,,w,⊥,,,t),(3.17),the,parallel,heat,flux,[Fs](z,,t),:=,πnsmsv3,ts,qs,(cid:107),(cid:90),∞,−∞,dw,(cid:90),∞,(cid:107),0,dw,⊥,w,w,(cid:107),(w2,(cid:107),+,w2,⊥,⊥,)Fs(z,,w,(cid:107),,,w,⊥,,,t),,(3.18),and,the,modified,charge,exchange,collision,operator,Cin[Fi,,Fn,,nn,,ui,:=,nnRin,(cid:107),−,(cid:107),,,un,(cid:20),Fi,−,∂,∂w,(cid:107),∂,∂w,⊥,,,vti,,vtn](z,,w,(cid:107),(cid:18),Fn,z,,v3,ti,v3,tn,(cid:20)(cid:18),un,ui,(cid:107),(cid:107),−,vti,(cid:18),v2,tn,v2,ti,−,(cid:20),w2,⊥,2,,,t),,,w,ui,⊥,(cid:107),−,vtn,w,(cid:107),2,+,1,+,,,vti,vtn,2(un,un,(cid:107),+,vti,vtn,w,(cid:107),(cid:18),v2,tn,v2,ti,−,2(un,(cid:107),−,3v2,ti,1,+,(cid:19),)2,ui,(cid:107),+,nnRin,+,nnRin,w,⊥,(cid:19)(cid:21),,,t,⊥,w,(cid:19)(cid:19),)2,(cid:107),(cid:21),Fi,ui,(cid:107),−,3v2,ti,(cid:21),Fi,.,(3.19),Note,that,the,differential,terms,in,this,modified,collision,operator,could,have,been,in-,i,and,˙Fi,,but,we,have,decided,to,make,cluded,in,the,definitions,of,the,coefficients,˙w,(cid:107),them,part,of,a,modified,collision,operator,instead,to,separate,the,effect,of,collisions,clearly.,This,split,should,not,be,taken,as,a,suggestion,on,how,to,implement,these,terms,in,a,code.,˙w,i,,⊥,1D,drift,kinetic,models,with,periodic,boundary,conditions,5,The,equations,for,the,neutrals,are,∂nn,∂t,+,∂,∂z,(cid:0)nnun,(cid:107),(cid:1),=,0,,nnmi,(cid:18),∂un,∂t,(cid:107),+,un,(cid:19),∂un,∂z,(cid:107),(cid:107),=,∂pn,∂z,(cid:107),−,+,ninnmiRin(ui,(cid:107),−,un,),,(cid:107),(3.20),(3.21),3,2,nnmivtn,(cid:18),∂vtn,∂t,+,un,∂vtn,∂z,(cid:107),(cid:19),=,∂qn,(cid:107),∂z,−,−,∂un,∂z,pn,(cid:107),3,4,nnnimiRin(v2,ti,−,v2,tn),nnnimiRin(un,(cid:107),−,)2,ui,(cid:107),(3.22),(cid:107),+,+,1,2,and,Cni.,Here,,we,have,defined,the,modified,charge,exchange,collision,operator,+,˙w,(cid:107),+,˙zn,+,˙w,⊥,⊥,n,n,=,˙Fn,+,∂Fn,∂t,∂Fn,∂z,∂Fn,∂w,(cid:107),∂Fn,∂w,(3.23),Cni[Fn,,Fi,,ni,,un,niRin,:=,(cid:107),−,(cid:107),,,ui,(cid:20),Fn,−,∂,∂w,(cid:107),∂,∂w,⊥,,,t),,,w,un,⊥,,,vtn,,vti](z,,w,(cid:107),(cid:18),v3,tn,v3,ti,(cid:20)(cid:18),ui,un,Fi,z,,+,(cid:107),(cid:107),−,vtn,(cid:18),v2,ti,v2,tn,−,(cid:20),w2,⊥,2,(cid:107),−,vti,w,(cid:107),2,1,+,,,vtn,vti,2(un,ui,(cid:107),+,vtn,vti,w,(cid:107),(cid:18),v2,ti,v2,tn,−,2(un,(cid:107),−,3v2,tn,1,+,(cid:19),)2,ui,(cid:107),+,niRin,+,niRin,w,⊥,(cid:19)(cid:21),,,t,⊥,w,(cid:19)(cid:19),)2,(cid:21),Fn,ui,(cid:107),(cid:107),−,3v2,tn,(cid:21),Fn,.,(3.24),Equations,(3.12),and,(3.23),for,Fi,and,Fn,are,constructed,such,that,conditions,(3.6),,(3.7),and,(3.8),are,satisfied,at,all,times,if,they,are,satisfied,at,t,=,0.,4.,Linear,test,One,possible,test,for,the,sets,of,1D,equations,described,above,is,the,evolution,of,small,perturbations,to,a,uniform,Maxwellian,equilibrium.,We,assume,the,following,form,for,the,ion,and,neutral,distribution,functions,,fs(z,,v,(cid:107),,,v,⊥,,,t),=,fM,s(v,,,v,⊥,(cid:107),),+,fs1(v,,,v,⊥,(cid:107),)[exp(ik,z,(cid:107),−,iωt),+,complex,conjugate],,(4.1),where,fM,s(v,,,v,⊥,(cid:107),),=,ns,(cid:19)3/2,(cid:18),mi,2πTh,(cid:32),−,+,v2,⊥,mi(v2,(cid:107),2Th,(cid:33),),.,exp,(4.2),Note,that,both,species,share,the,same,constant,temperature,Th.,To,ensure,that,the,potential,is,small,,we,assume,ni,=,Ne.,Since,the,perturbations,fs1(v,can,be,linearized,to,give,,,v,⊥,(cid:107),),and,φ,are,small,,equations,(2.14),,(2.15),and,(2.16),v,(k,(cid:107),(cid:107),−,ω,−,innRin)fi1,+,iniRinfn1,=,eφ,Th,−,k,v,(cid:107),(cid:107),fM,i,+,iRin(nn1fM,i,−,ni1fM,n),,(4.3),innRinfi1,+,(k,(cid:107),v,ω,iniRin)fn1,=,iRin(ni1fM,n,−,−,(cid:107),−,nn1fM,i),and,ni1,ni,=,eφ,Te,.,(4.4),(4.5),6,Felix,I.,Parra,,Michael,Barnes,and,Michael,Hardman,Figure,1.,Solutions,to,the,dispersion,relation,(4.7):,acoustic,waves,(solid,lines),and,non-prop-,agating,modes,(dashed,lines).,(a),Real,frequency,ωr,:=,Re(ω),and,(b),growth,rate,γ,:=,Im(ω),as,functions,of,the,charge,exchange,collision,frequency,(ni,+,nn)Rin,for,Te/Th,=,1,and,several,values,of,the,parameter,ni/(ni,+,nn).,(c),Real,frequency,ωr,and,(d),growth,rate,γ,as,functions,of,the,the,charge,exchange,collision,frequency,(ni,+,nn)Rin,for,ni/(ni,+,nn),=,1/2,and,several,values,of,the,parameter,Te/Th.,Here,,we,have,defined,the,perturbations,to,the,density,as,ns1,=,2π,(cid:90),∞,−∞,dv,(cid:90),∞,(cid:107),0,dv,fs1.,v,⊥,⊥,(4.6),Solving,for,the,functions,fi1,and,fn1,as,functions,of,ni1,and,nn1,and,then,integrating,fi1,and,fn1,over,velocity,space,,we,find,the,equations,(cid:18),Aii,Ain,Ani,Ann,(cid:19),(cid:18),ni1,nn1,(cid:19),(cid:18),0,0,=,(cid:19),,,(4.7),where,the,elements,of,the,matrix,are,Te,Th,ni,ni,+,nn,Aii,=,1,+,Te,Th,+,ζZ(ζ),+,nn,ni,+,nn,(cid:20)(cid:18),1,+,(cid:19),Te,Th,ζin,−,ζ,(cid:21),Z(ζin),(4.8),Ain,=,−,(cid:26)(cid:20)(cid:18),ni,(ζin,−,ni,+,nn,(cid:19),ζ)Z(ζin),(cid:21),Ani,=,nn,ni,+,nn,−,1,+,Te,Th,ζin,−,ζ,Z(ζin),(cid:27),ζZ(ζ),Te,Th,−,and,Here,,we,have,defined,Ann,=,1,+,ni,ni,+,nn,(ζin,−,ζ)Z(ζin).,(4.9),(4.10),(4.11),with,vth,:=,(cid:112),k,|,2Th/mi,,and,we,have,used,the,plasma,dispersion,function,(Fried,&,Conte,vth,(cid:107)|,(cid:107)|,k,|,ζ,:=,ω,vth,,,ζin,:=,ω,+,i(ni,+,nn)Rin,,,(4.12),00.511.5200.511.513/41/21/4100.511.52-2.5-2-1.5-1-0.5000.511.5211.21.41.61.821/41/212400.511.52-2.5-2-1.5-1-0.5000.511.5200.511.513/41/21/4100.511.52-2.5-2-1.5-1-0.5000.511.5211.21.41.61.821/41/212400.511.52-2.5-2-1.5-1-0.50AAACCnicbVDLSgMxFM34rPU16tJNtAgVoc6IqMuiG5dV7APaIWTStA3NZIYkUyjTWbvxV9y4UMStX+DOvzFtZ6GtBwKHc+7l5hw/4kxpx/m2FhaXlldWc2v59Y3NrW17Z7emwlgSWiUhD2XDx4pyJmhVM81pI5IUBz6ndb9/M/brAyoVC8WDHkbUC3BXsA4jWBsJ2QdFgRg8gQKJY3iPEibS01EftUZwBAco0b0U2QWn5EwA54mbkQLIUEH2V6sdkjigQhOOlWq6TqS9BEvNCKdpvhUrGmHSx13aNFTggCovmURJ4ZFR2rATSvOEhhP190aCA6WGgW8mA6x7atYbi/95zVh3rjwTL4o1FWR6qBNzqEM47gW2maRE86EhmEhm/gpJD0tMtGkvb0pwZyPPk9pZyb0oOXfnhfJ1VkcO7INDUAQuuARlcAsqoAoIeATP4BW8WU/Wi/VufUxHF6xsZw/8gfX5A2NImXk=(ni+nn)Rin/|kk|vthAAACCnicbVDLSgMxFM34rPU16tJNtAgVoc6IqMuiG5dV7APaIWTStA3NZIYkUyjTWbvxV9y4UMStX+DOvzFtZ6GtBwKHc+7l5hw/4kxpx/m2FhaXlldWc2v59Y3NrW17Z7emwlgSWiUhD2XDx4pyJmhVM81pI5IUBz6ndb9/M/brAyoVC8WDHkbUC3BXsA4jWBsJ2QdFgRg8gQKJY3iPEibS01EftUZwBAco0b0U2QWn5EwA54mbkQLIUEH2V6sdkjigQhOOlWq6TqS9BEvNCKdpvhUrGmHSx13aNFTggCovmURJ4ZFR2rATSvOEhhP190aCA6WGgW8mA6x7atYbi/95zVh3rjwTL4o1FWR6qBNzqEM47gW2maRE86EhmEhm/gpJD0tMtGkvb0pwZyPPk9pZyb0oOXfnhfJ1VkcO7INDUAQuuARlcAsqoAoIeATP4BW8WU/Wi/VufUxHF6xsZw/8gfX5A2NImXk=(ni+nn)Rin/|kk|vthAAACCnicbVDLSgMxFM34rPU16tJNtAgVoc6IqMuiG5dV7APaIWTStA3NZIYkUyjTWbvxV9y4UMStX+DOvzFtZ6GtBwKHc+7l5hw/4kxpx/m2FhaXlldWc2v59Y3NrW17Z7emwlgSWiUhD2XDx4pyJmhVM81pI5IUBz6ndb9/M/brAyoVC8WDHkbUC3BXsA4jWBsJ2QdFgRg8gQKJY3iPEibS01EftUZwBAco0b0U2QWn5EwA54mbkQLIUEH2V6sdkjigQhOOlWq6TqS9BEvNCKdpvhUrGmHSx13aNFTggCovmURJ4ZFR2rATSvOEhhP190aCA6WGgW8mA6x7atYbi/95zVh3rjwTL4o1FWR6qBNzqEM47gW2maRE86EhmEhm/gpJD0tMtGkvb0pwZyPPk9pZyb0oOXfnhfJ1VkcO7INDUAQuuARlcAsqoAoIeATP4BW8WU/Wi/VufUxHF6xsZw/8gfX5A2NImXk=(ni+nn)Rin/|kk|vthAAACCnicbVDLSgMxFM34rPU16tJNtAgVoc6IqMuiG5dV7APaIWTStA3NZIYkUyjTWbvxV9y4UMStX+DOvzFtZ6GtBwKHc+7l5hw/4kxpx/m2FhaXlldWc2v59Y3NrW17Z7emwlgSWiUhD2XDx4pyJmhVM81pI5IUBz6ndb9/M/brAyoVC8WDHkbUC3BXsA4jWBsJ2QdFgRg8gQKJY3iPEibS01EftUZwBAco0b0U2QWn5EwA54mbkQLIUEH2V6sdkjigQhOOlWq6TqS9BEvNCKdpvhUrGmHSx13aNFTggCovmURJ4ZFR2rATSvOEhhP190aCA6WGgW8mA6x7atYbi/95zVh3rjwTL4o1FWR6qBNzqEM47gW2maRE86EhmEhm/gpJD0tMtGkvb0pwZyPPk9pZyb0oOXfnhfJ1VkcO7INDUAQuuARlcAsqoAoIeATP4BW8WU/Wi/VufUxHF6xsZw/8gfX5A2NImXk=(ni+nn)Rin/|kk|vth