The processes of electric ion extraction from plasma induced by pulse lasers are simulated by particle-in-cell(PIC)method and hybrid-PIC method.A new calculation scheme named preprocessing hybrid-PIC is presented beca...The processes of electric ion extraction from plasma induced by pulse lasers are simulated by particle-in-cell(PIC)method and hybrid-PIC method.A new calculation scheme named preprocessing hybrid-PIC is presented because neither of the two methods above is omnipotent,especially under the circumstance of high initial plasma density.The new scheme provides credible results with less computational consumption than PIC method in both one-and two-dimensional simulations,except forΠ-type electrode configuration.The simulation results show that the M-type performs best in all electrode configurations in both high-density and low-density plasma conditions.展开更多
Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging ...Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.展开更多
A two-and three-dimensional velocity space axisymmetric hybrid-PIC model of Hall thruster discharge called Hybrid2D has been developed.The particle-in-cell(PIC) method was used for neutrals and ions(heavy species),and...A two-and three-dimensional velocity space axisymmetric hybrid-PIC model of Hall thruster discharge called Hybrid2D has been developed.The particle-in-cell(PIC) method was used for neutrals and ions(heavy species),and fluid dynamics on a magnetic field-aligned(MFA) mesh was used for electrons.A time-saving method for heavy species moment interpolation on a MFA mesh was developed.The method comprises using regular rectangle and irregular triangle meshes,connected to each other on a pre-processing stage.The electron fluid model takes into account neither inertia terms nor viscous terms and includes an electron temperature equation with a heat flux term.The developed model was used to calculate all heavy species moments up to the third one in a stationary case.The analysis of the viscosity and the heat flux impact on the force and energy balance has shown that for the calculated geometry of the Hall thruster,the viscosity and the heat flux terms have the same magnitude as the other terms and could not be omitted.Also,it was shown that the heat flux is not proportional to the temperature gradient and,consequently,the highest moments should be calculated to close the neutral fluid equation system.At the same time,ions can only be modeled as a cold non-viscous fluid when the sole aim of modeling is the calculation of the operating parameters or distribution of the local parameters along the centerline of the discharge channel.This is because the magnitude of the viscosity and the temperature gradient terms are negligible at the centerline.However,when a simulation’s focus is either on the radial divergence of the plume or on magnetic pole erosion,three components of the ion temperature should be taken into consideration.The non-diagonal terms of ion pressure tensor have a lower impact than the diagonal terms.According to the study,a zero heat flux condition could be used to close the ion equation system in calculated geometry.展开更多
电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte C...电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。展开更多
Double tearing modes(DTMs),induced by double current sheet configurations or two neighboring rational surfaces with the same safety factor in tokamaks,are widely observed in solar,space,and fusion plasmas.In this pape...Double tearing modes(DTMs),induced by double current sheet configurations or two neighboring rational surfaces with the same safety factor in tokamaks,are widely observed in solar,space,and fusion plasmas.In this paper,the evolution of DTMs without a guide field is investigated numerically using a hybrid model(electron fluid+ion PIC).The overall evolution processes of DTMs are qualitatively consistent with previous works using other models.The particle dynamics during the evolution of DTMs is analyzed in detail.Behaviors of ions and electrons present different characteristics around the reconnection region which gives rise to Hall effects producing the out-of-plane quadrupole magnetic field.In the explosive reconnection process with interactions between two DTMs islands,the asymmetric drive and the thin current layer feature lead to the emergence of secondary magnetic islands which develop with the late evolution of the DTMs.展开更多
为了能够将脉冲等离子体推力器成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法一体化模拟实验室PPT羽流。验证计算...为了能够将脉冲等离子体推力器成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法一体化模拟实验室PPT羽流。验证计算显示该模型具有一体化模拟脉冲等离子体推力器羽流的能力。对不同初始放电能量下的羽流场进行模拟,给出了离子、中性粒子、电子温度、轴线上质量流率和出口平面返流质量流率的变化情况。计算结果显示高放电能量下返流量更大,同时中性粒子在返流中所占比例也越大。展开更多
The parametric instability related to ion motion and the resulting cross-beam energy transfer are important aspects in the physics of inertial confinement fusion.The numerical simulation of the above physical problems...The parametric instability related to ion motion and the resulting cross-beam energy transfer are important aspects in the physics of inertial confinement fusion.The numerical simulation of the above physical problems still faces great technical challenges.This paper introduces a 2D hybrid-kinetic particle-in-cell(PIC)code,CBETor.In this code,the motion of ions is described by the kinetic method,the motion of electrons is described by the simplified fluid method and the propagation of laser in plasma is described by solving the wave equation.We use CBETor and the popular fully kinetic PIC code EPOCH to simulate the stimulated Brillouin scattering and cross-beam energy transfer process,respectively.The physical images are in good agreement,but CBETor can significantly reduce the amount of calculation.With the premise of correctly simulating the ion dynamics,our hybrid-kinetic code can effectively suppress the noise of numerical simulation and significantly expand the simulation scale of physical problems.CBETor is very suitable for simulating the physical process dominated by ion motion in the interaction of medium intensity laser and underdense plasma.展开更多
目前微小卫星正在积极的发展中,脉冲等离子体推力器是其推进系统的一个重要发展方向,为了能够将PPT成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型用于DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混...目前微小卫星正在积极的发展中,脉冲等离子体推力器是其推进系统的一个重要发展方向,为了能够将PPT成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型用于DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法模拟PPT羽流的入口条件计算,一体化模拟实验室PPT羽流,对不同电容情况下的羽流场进行模拟,并与实验结果进行了比较。计算结果显示高电容下带来更高的质量流量,更高的中性粒子的含量,同时返流的影响域更广。在推力器入口附近,CEX碰撞与一般碰撞形式共同存在,且频率很高,在羽流外围,CEX碰撞成为碰撞的主要形式。展开更多
将三维MHD双温入口模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法,模拟实验室PPT样机羽流。验证计算显示,该模型具有模拟脉冲等离子体推力器羽流的能力。对不同初始电压和电...将三维MHD双温入口模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法,模拟实验室PPT样机羽流。验证计算显示,该模型具有模拟脉冲等离子体推力器羽流的能力。对不同初始电压和电容下的羽流场进行了模拟,给出了出口平面返流质量流率的变化情况和20μs羽流中离子和CEX离子的分布情况。计算结果显示,高能量状态对应高质量流率,高质量流率对应高动量,高动量离子和中性粒子对航天器撞击会造成更强影响。展开更多
文摘The processes of electric ion extraction from plasma induced by pulse lasers are simulated by particle-in-cell(PIC)method and hybrid-PIC method.A new calculation scheme named preprocessing hybrid-PIC is presented because neither of the two methods above is omnipotent,especially under the circumstance of high initial plasma density.The new scheme provides credible results with less computational consumption than PIC method in both one-and two-dimensional simulations,except forΠ-type electrode configuration.The simulation results show that the M-type performs best in all electrode configurations in both high-density and low-density plasma conditions.
基金supported by the National Natural Science Foundation of China(NNSFC)grants 42074202,42274196Strategic Priority Research Program of Chinese Academy of Sciences grant XDB41000000ISSI-BJ International Team Interaction between magnetic reconnection and turbulence:From the Sun to the Earth。
文摘Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.
文摘A two-and three-dimensional velocity space axisymmetric hybrid-PIC model of Hall thruster discharge called Hybrid2D has been developed.The particle-in-cell(PIC) method was used for neutrals and ions(heavy species),and fluid dynamics on a magnetic field-aligned(MFA) mesh was used for electrons.A time-saving method for heavy species moment interpolation on a MFA mesh was developed.The method comprises using regular rectangle and irregular triangle meshes,connected to each other on a pre-processing stage.The electron fluid model takes into account neither inertia terms nor viscous terms and includes an electron temperature equation with a heat flux term.The developed model was used to calculate all heavy species moments up to the third one in a stationary case.The analysis of the viscosity and the heat flux impact on the force and energy balance has shown that for the calculated geometry of the Hall thruster,the viscosity and the heat flux terms have the same magnitude as the other terms and could not be omitted.Also,it was shown that the heat flux is not proportional to the temperature gradient and,consequently,the highest moments should be calculated to close the neutral fluid equation system.At the same time,ions can only be modeled as a cold non-viscous fluid when the sole aim of modeling is the calculation of the operating parameters or distribution of the local parameters along the centerline of the discharge channel.This is because the magnitude of the viscosity and the temperature gradient terms are negligible at the centerline.However,when a simulation’s focus is either on the radial divergence of the plume or on magnetic pole erosion,three components of the ion temperature should be taken into consideration.The non-diagonal terms of ion pressure tensor have a lower impact than the diagonal terms.According to the study,a zero heat flux condition could be used to close the ion equation system in calculated geometry.
文摘电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。
基金supported by National Natural Science Foundation of China(Nos.11705039,11975087,and 42011530086)the China Postdoctoral Science Foundation(No.2018M631918)the Heilongjiang Postdoctoral Foundation(No.184744)。
文摘Double tearing modes(DTMs),induced by double current sheet configurations or two neighboring rational surfaces with the same safety factor in tokamaks,are widely observed in solar,space,and fusion plasmas.In this paper,the evolution of DTMs without a guide field is investigated numerically using a hybrid model(electron fluid+ion PIC).The overall evolution processes of DTMs are qualitatively consistent with previous works using other models.The particle dynamics during the evolution of DTMs is analyzed in detail.Behaviors of ions and electrons present different characteristics around the reconnection region which gives rise to Hall effects producing the out-of-plane quadrupole magnetic field.In the explosive reconnection process with interactions between two DTMs islands,the asymmetric drive and the thin current layer feature lead to the emergence of secondary magnetic islands which develop with the late evolution of the DTMs.
文摘为了能够将脉冲等离子体推力器成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法一体化模拟实验室PPT羽流。验证计算显示该模型具有一体化模拟脉冲等离子体推力器羽流的能力。对不同初始放电能量下的羽流场进行模拟,给出了离子、中性粒子、电子温度、轴线上质量流率和出口平面返流质量流率的变化情况。计算结果显示高放电能量下返流量更大,同时中性粒子在返流中所占比例也越大。
基金supported by National Natural Science Foundation of China (Nos. 11774430, 11875091, 12075157 and 11975062)
文摘The parametric instability related to ion motion and the resulting cross-beam energy transfer are important aspects in the physics of inertial confinement fusion.The numerical simulation of the above physical problems still faces great technical challenges.This paper introduces a 2D hybrid-kinetic particle-in-cell(PIC)code,CBETor.In this code,the motion of ions is described by the kinetic method,the motion of electrons is described by the simplified fluid method and the propagation of laser in plasma is described by solving the wave equation.We use CBETor and the popular fully kinetic PIC code EPOCH to simulate the stimulated Brillouin scattering and cross-beam energy transfer process,respectively.The physical images are in good agreement,but CBETor can significantly reduce the amount of calculation.With the premise of correctly simulating the ion dynamics,our hybrid-kinetic code can effectively suppress the noise of numerical simulation and significantly expand the simulation scale of physical problems.CBETor is very suitable for simulating the physical process dominated by ion motion in the interaction of medium intensity laser and underdense plasma.
文摘目前微小卫星正在积极的发展中,脉冲等离子体推力器是其推进系统的一个重要发展方向,为了能够将PPT成功地运用于空间,需对其羽流进行研究。将一维MHD双温放电模型用于DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法模拟PPT羽流的入口条件计算,一体化模拟实验室PPT羽流,对不同电容情况下的羽流场进行模拟,并与实验结果进行了比较。计算结果显示高电容下带来更高的质量流量,更高的中性粒子的含量,同时返流的影响域更广。在推力器入口附近,CEX碰撞与一般碰撞形式共同存在,且频率很高,在羽流外围,CEX碰撞成为碰撞的主要形式。
文摘将三维MHD双温入口模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法,模拟实验室PPT样机羽流。验证计算显示,该模型具有模拟脉冲等离子体推力器羽流的能力。对不同初始电压和电容下的羽流场进行了模拟,给出了出口平面返流质量流率的变化情况和20μs羽流中离子和CEX离子的分布情况。计算结果显示,高能量状态对应高质量流率,高质量流率对应高动量,高动量离子和中性粒子对航天器撞击会造成更强影响。