For investigating efficiently the stagnation kinetic-process of Z-pinch,we develop a novel modified electrostatic implicit particle-in-cell algorithm in radial one-dimension for Z-pinch simulation in which a small-ang...For investigating efficiently the stagnation kinetic-process of Z-pinch,we develop a novel modified electrostatic implicit particle-in-cell algorithm in radial one-dimension for Z-pinch simulation in which a small-angle cumulative binary collision algorithm is used.In our algorithm,the electric field in z-direction is solved by a parallel electrode-plate model,the azimuthal magnetic field is obtained by Ampere’s law,and the term for charged particle gyromotion is approximated by the cross product of the averaged velocity and magnetic field.In simulation results of 2 MA deuterium plasma shell Zpinch,the mass-center implosion trajectory agrees generally with that obtained by one-dimensional MHD simulation,and the plasma current also closely aligns with the external current.The phase space diagrams and radial-velocity probability distributions of ions and electrons are obtained.The main kinetic characteristic of electron motion is thermal equilibrium and oscillation,which should be oscillated around the ions,while that of ion motion is implosion inwards.In the region of stagnation radius,the radial-velocity probability distribution of ions transits from the non-equilibrium to equilibrium state with the current increasing,while of electrons is basically the equilibrium state.When the initial ion density and current peak are not high enough,the ions may not reach their thermal equilibrium state through collisions even in its stagnation phase.展开更多
基金the graduated students Zhixing Feng,Xiaoqiang Zhang,and Deli Fang for their excellent works to develop the PIC simulation codes of Z-pinch.This research was partly supported by the National Natural Science Foundation of China(Grant Nos.11675025 and 11135007)the Innovation Project of China Academy of Engineering Physics(Grant No.CX2019030).
文摘For investigating efficiently the stagnation kinetic-process of Z-pinch,we develop a novel modified electrostatic implicit particle-in-cell algorithm in radial one-dimension for Z-pinch simulation in which a small-angle cumulative binary collision algorithm is used.In our algorithm,the electric field in z-direction is solved by a parallel electrode-plate model,the azimuthal magnetic field is obtained by Ampere’s law,and the term for charged particle gyromotion is approximated by the cross product of the averaged velocity and magnetic field.In simulation results of 2 MA deuterium plasma shell Zpinch,the mass-center implosion trajectory agrees generally with that obtained by one-dimensional MHD simulation,and the plasma current also closely aligns with the external current.The phase space diagrams and radial-velocity probability distributions of ions and electrons are obtained.The main kinetic characteristic of electron motion is thermal equilibrium and oscillation,which should be oscillated around the ions,while that of ion motion is implosion inwards.In the region of stagnation radius,the radial-velocity probability distribution of ions transits from the non-equilibrium to equilibrium state with the current increasing,while of electrons is basically the equilibrium state.When the initial ion density and current peak are not high enough,the ions may not reach their thermal equilibrium state through collisions even in its stagnation phase.