Based on a semi-resolved CFD-DEM coupling method,this study proposed a method that uses the minimum distance between the fluid grid and the particle boundary as a reference value to determine the degree of influence o...Based on a semi-resolved CFD-DEM coupling method,this study proposed a method that uses the minimum distance between the fluid grid and the particle boundary as a reference value to determine the degree of influence of the target fluid grid on the particle's drag force.A fluidized bed of rod-like particles was chosen as a typical case to investigate the effect of different fluid grid scales on various fluidized bed characteristic parameters.The calculation performance of the semi-resolved and unre-solved CFD-DEM coupling algorithm on key fluidized bed characteristic parameters such as average pressure drop,particle frequency distribution with bed height,and particle orientation distribution were compared.It was found that the semi-resolved CFD-DEM coupling algorithm gradually obtained results with higher consistency with decreasing fluid grid scale for key parameters such as particle frequency distribution with bed height,particle orientation distribution,and time-history mixing index,exhibiting a phenomenon similar to grid independence in fluid simulation.By comparing with experimental results,it was verified that the semi-resolved CFD-DEM coupling algorithm can be applied to simulate multi-granular gas-solid systems with fluid grid scales equivalent to particle scales.This algorithm solves the limitation of fluid grid scale in the unresolved CFD-DEM coupling framework and improves the grid adaptability of the CFD-DEM coupling simulation algorithm.展开更多
The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on th...The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.展开更多
This work investigated the pyrolysis reaction of waste resin in a fluidized bed reactor.It was found that the pyrolysis-generated ash would adhere to the surface of ceramic particles,causing particle agglomeration and...This work investigated the pyrolysis reaction of waste resin in a fluidized bed reactor.It was found that the pyrolysis-generated ash would adhere to the surface of ceramic particles,causing particle agglomeration and defluidization.Adding kaolin could effectively inhibit the particle agglomeration during the fluidized pyrolysis reaction through physical isolation and chemical reaction.On the one hand,kaolin could form a coating layer on the surface of ceramic particles to prevent the adhesion of organic ash generated by the pyrolysis of resin.On the other hand,when a sufficient amount of kaolin(-0.2%(mass))was added,the activated kaolin could fully contact with the Na+ ions generated by the pyrolysis of resin and react to form a high-melting aluminosilicate mineral(nepheline),which could reduce the formation of low-melting-point sodium sulfate and thereby avoid the agglomeration of ceramic particles.展开更多
The existing indoor fusion positioning methods based on Pedestrian Dead Reckoning(PDR)and geomagnetic technology have the problems of large initial position error,low sensor accuracy,and geomagnetic mismatch.In this s...The existing indoor fusion positioning methods based on Pedestrian Dead Reckoning(PDR)and geomagnetic technology have the problems of large initial position error,low sensor accuracy,and geomagnetic mismatch.In this study,a novel indoor fusion positioning approach based on the improved particle filter algorithm by geomagnetic iterative matching is proposed,where Wi-Fi,PDR,and geomagnetic signals are integrated to improve indoor positioning performances.One important contribution is that geomagnetic iterative matching is firstly proposed based on the particle filter algorithm.During the positioning process,an iterative window and a constraint window are introduced to limit the particle generation range and the geomagnetic matching range respectively.The position is corrected several times based on geomagnetic iterative matching in the location correction stage when the pedestrian movement is detected,which made up for the shortage of only one time of geomagnetic correction in the existing particle filter algorithm.In addition,this study also proposes a real-time step detection algorithm based on multi-threshold constraints to judge whether pedestrians are moving,which satisfies the real-time requirement of our fusion positioning approach.Through experimental verification,the average positioning accuracy of the proposed approach reaches 1.59 m,which improves 33.2%compared with the existing particle filter fusion positioning algorithms.展开更多
The traditional standard wet sieving method uses steel sieves with aperture?0.063 mm and can only determine the particle size distribution(PSD)of gravel and sand in general soil.This paper extends the traditional meth...The traditional standard wet sieving method uses steel sieves with aperture?0.063 mm and can only determine the particle size distribution(PSD)of gravel and sand in general soil.This paper extends the traditional method and presents an extended wet sieving method.The extended method uses both the steel sieves and the nylon filter cloth sieves.The apertures of the cloth sieves are smaller than 0.063 mm and equal 0.048 mm,0.038 mm,0.014 mm,0.012 mm,0.0063 mm,0.004 mm,0.003 mm,0.002 mm,and 0.001 mm,respectively.The extended method uses five steps to separate the general soil into many material sub-groups of gravel,sand,silt and clay with known particle size ranges.The complete PSD of the general soil is then calculated from the dry masses of the individual material sub-groups.The extended method is demonstrated with a general soil of completely decomposed granite(CDG)in Hong Kong,China.The silt and clay materials with different particle size ranges are further examined,checked and verified using stereomicroscopic observation,physical and chemical property tests.The results further confirm the correctness of the extended wet sieving method.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
An E//B neutral particle analyzer(NPA)has been designed and is under development at Sichuan University and Southwestern Institute of Physics.The main purpose of the E//B NPA is to measure the distribution function of ...An E//B neutral particle analyzer(NPA)has been designed and is under development at Sichuan University and Southwestern Institute of Physics.The main purpose of the E//B NPA is to measure the distribution function of fast ions in the HL-2A/3 tokamak.The E//B NPA contains three main units,i.e.the stripping unit,the analyzing unit and the detection unit.A gas stripping chamber was adopted as the stripping unit.The results of the simulations and beam tests for the stripping chamber are presented.Parallel electric and magnetic fields provided by a NdFeB permanent magnet and two parallel electric plates were designed and constructed for the analyzing unit.The calibration of the magnetic and electric fields was performed using a 50 kV electron cyclotron resonance ion source(ECRIS)platform.The detection unit consists of 32lutetium-yttrium oxyorthosilicate(LYSO)detector modules arranged in two rows.The response functions ofα,hydrogen ions(H^(+),H_(2)^(+)and H_(3)^(+))andγfor a detector module were measured with^(241)Am,^(137)Cs and^(152)Eu sources together with the 50 kV ECRIS platform.The overall results indicate that the designed E//B NPA device is capable of measuring the intensity of neutral hydrogen and deuteron atoms with energy higher than 20 keV.展开更多
Deformable gel particles(DGPs) possess the capability of deep profile control and flooding. However, the deep migration behavior and plugging mechanism along their path remain unclear. Breakage, an inevitable phenomen...Deformable gel particles(DGPs) possess the capability of deep profile control and flooding. However, the deep migration behavior and plugging mechanism along their path remain unclear. Breakage, an inevitable phenomenon during particle migration, significantly impacts the deep plugging effect. Due to the complexity of the process, few studies have been conducted on this subject. In this paper, we conducted DGP flow experiments using a physical model of a multi-point sandpack under various injection rates and particle sizes. Particle size and concentration tests were performed at each measurement point to investigate the transportation behavior of particles in the deep part of the reservoir. The residual resistance coefficient and concentration changes along the porous media were combined to analyze the plugging performance of DGPs. Furthermore, the particle breakage along their path was revealed by analyzing the changes in particle size along the way. A mathematical model of breakage and concentration changes along the path was established. The results showed that the passage after breakage is a significant migration behavior of particles in porous media. The particles were reduced to less than half of their initial size at the front of the porous media. Breakage is an essential reason for the continuous decreases in particle concentration, size, and residual resistance coefficient. However, the particles can remain in porous media after breakage and play a significant role in deep plugging. Higher injection rates or larger particle sizes resulted in faster breakage along the injection direction, higher degrees of breakage, and faster decreases in residual resistance coefficient along the path. These conditions also led to a weaker deep plugging ability. Smaller particles were more evenly retained along the path, but more particles flowed out of the porous media, resulting in a poor deep plugging effect. The particle size is a function of particle size before injection, transport distance, and different injection parameters(injection rate or the diameter ratio of DGP to throat). Likewise, the particle concentration is a function of initial concentration, transport distance, and different injection parameters. These models can be utilized to optimize particle injection parameters, thereby achieving the goal of fine-tuning oil displacement.展开更多
A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDO...A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.展开更多
It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on b...It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on bubble motion and gas-liquid mass transfer process in a bubble column.Particle image velocimetry(PIV)technology was employed to analyze the flow field and bubble motion behavior,and dynamic oxygen absorption technology was used to measure the gas-liquid volumetric mass transfer coefficient(kLa).The effect of embedding wire mesh,adding solid particles,and wire mesh coupling solid particles on the flow characteristic and kLa were analyzed and compared.The results show that the gas-liquid interface area increases by 33%-72%when using the wire mesh coupling solid particles strategy compared to the gas-liquid two-phase flow,which is superior to the other two strengthening methods.Compared with the system without reinforcement,kLa in the bubble column increased by 0.5-1.8 times with wire mesh coupling solid particles method,which is higher than the sum of kLa increases with inserting wire mesh and adding particles,and the coupling reinforcement mechanism for affecting gas-liquid mass transfer process was discussed to provide a new idea for enhancing gas-liquid mass transfer.展开更多
Realizing the full potential of ultrahigh-intensity lasers for particle and radiation generation will require multi-beam arrangements due to technology limitations.Here,we investigate how to optimize their coupling wi...Realizing the full potential of ultrahigh-intensity lasers for particle and radiation generation will require multi-beam arrangements due to technology limitations.Here,we investigate how to optimize their coupling with solid targets.Experimentally,we show that overlapping two intense lasers in a mirror-like configuration onto a solid with a large preplasma can greatly improve the generation of hot electrons at the target front and ion acceleration at the target backside.The underlying mechanisms are analyzed through multidimensional particle-in-cell simulations,revealing that the self-induced magnetic fields driven by the two laser beams at the target front are susceptible to reconnection,which is one possible mechanism to boost electron energization.In addition,the resistive magnetic field generated during the transport of the hot electrons in the target bulk tends to improve their collimation.Our simulations also indicate that such effects can be further enhanced by overlapping more than two laser beams.展开更多
Accurate vehicle dynamic information plays an important role in vehicle driving safety.However,due to the characteristics of high mobility and multiple controllable degrees of freedom of drive-by-wire chassis vehicles...Accurate vehicle dynamic information plays an important role in vehicle driving safety.However,due to the characteristics of high mobility and multiple controllable degrees of freedom of drive-by-wire chassis vehicles,the current mature application of traditional vehicle state estimation algorithms can not meet the requirements of drive-by-wire chassis vehicle state estimation.This paper proposes a state estimation method for drive-by-wire chassis vehicle based on the dual unscented particle filter algorithm,which make full use of the known advantages of the four-wheel drive torque and steer angle parameters of the drive-by-wire chassis vehicle.In the dual unscented particle filter algorithm,two unscented particle filter transfer information to each other,observe the vehicle state information and the tire force parameter information of the four wheels respectively,which reduce the influence of parameter uncertainty and model parameter changes on the estimation accuracy during driving.The performance with the dual unscented particle filter algorithm,which is analyzed in terms of the time-average square error,is superior of the unscented Kalman filter algorithm.The effectiveness of the algorithm is further verified by driving simulator test.In this paper,a vehicle state estimator based on dual unscented particle filter algorithm was proposed for the first time to improve the estimation accuracy of vehicle parameters and states.展开更多
Sediment accumulation on the bed of open sewers and drains reduces hydraulic efficiency and can cause localized flooding.Slotted invert traps installed underneath the bed of open sewers and drains can eliminate sedime...Sediment accumulation on the bed of open sewers and drains reduces hydraulic efficiency and can cause localized flooding.Slotted invert traps installed underneath the bed of open sewers and drains can eliminate sediment build-up by catching sediment load.Previous three-dimensional(3D)computational studies have examined the particle trapping performance of invert traps of different shapes and depths under varied sediment and flow conditions,considering particles as spheres.For two-dimensional and 3D numerical modeling,researchers assumed the lid geometry to be a thin line and a plane,respectively.In this 3D numerical study,the particle trapping efficiency of a slotted irregular hexagonal invert trap fitted at the flume bottom was examined by incorporating the particle shape factor of non-spherical sewage solid particles and the thicknesses of upstream and downstream lids over the trap in the discrete phase model of the ANSYS Fluent 2020 R1 software.The volume of fluid(VOF)and the realizable k-turbulence models were used to predict the velocity field.The two-dimensional particle image velocimetry(PIV)was used to measure the velocity field inside the invert trap.The results showed that the thicknesses of upstream and downstream lids affected the velocity field and turbulent kinetic energy at all flow depths.The joint impact of the particle shape factor and lid thickness on the trap efficiency was significant.When both the lid thickness and particle shape factor were considered in the numerical modeling,trap efficiencies were underestimated,with relative errors of-8.66%to-0.65%in comparison to the experimental values of Mohsin and Kaushal(2017).They were also lower than the values predicted by Mohsin and Kaushal(2017),which showed an overall overestimation with errors of-2.3%to 17.4%.展开更多
Wind erosion is a geomorphic process in arid and semi-arid areas and has substantial implications for regional climate and desertification.In the Columbia Plateau of northwestern United States,the emissions from fine ...Wind erosion is a geomorphic process in arid and semi-arid areas and has substantial implications for regional climate and desertification.In the Columbia Plateau of northwestern United States,the emissions from fine particles of loessial soils often contribute to the exceedance of inhalable particulate matter(PM)with an aerodynamic diameter of 10μm or less(PM10)according to the air quality standards.However,little is known about the threshold friction velocity(TFV)for particles of different sizes that comprise these soils.In this study,soil samples of two representative soil types(Warden sandy loam and Ritzville silt loam)collected from the Columbia Plateau were sieved to seven particle size fractions,and an experiment was then conducted to determine the relationship between TFV and particle size fraction.The results revealed that soil particle size significantly affected the initiation of soil movement and TFV;TFV ranged 0.304-0.844 and 0.249-0.739 m/s for different particle size fractions of Ritzville silt loam and Warden sandy loam,respectively.PM10 and total suspended particulates(TSP)emissions from a bed of 63-90μm soil particles were markedly higher for Warden sandy loam than for Ritzville silt loam.Together with the lower TFV of Warden sandy loam,dust emissions from fine particles(<100μm in diameter)of Warden sandy loam thus may be a main contributor to dust in the region's atmosphere,since the PM10 emissions from the soil erosion surfaces and its ensuing suspension within the atmosphere constitute an essential process of soil erosion in the Columbia Plateau.Developing and implementing strategic land management practices on sandy loam soils is therefore necessary to control dust emissions in the Columbia Plateau.展开更多
The problems associated with vibrations of viaducts and low-frequency structural noise radiation caused by train excitation continue to increase in importance.A new floating-slab track vibration isolator-non-obstructi...The problems associated with vibrations of viaducts and low-frequency structural noise radiation caused by train excitation continue to increase in importance.A new floating-slab track vibration isolator-non-obstructive particle damping-phononic crystal vibration isolator is proposed herein,which uses the particle damping vibration absorption technology and bandgap vibration control theory.The vibration reduction performance of the NOPD-PCVI was analyzed from the perspective of vibration control.The paper explores the structure-borne noise reduction performance of the NOPD-PCVIs installed on different bridge structures under varying service conditions encountered in practical engineering applications.The load transferred to the bridge is obtained from a coupled train-FST-bridge analytical model considering the different structural parameters of bridges.The vibration responses are obtained using the finite element method,while the structural noise radiation is simulated using the frequency-domain boundary element method.Using the particle swarm optimization algorithm,the parameters of the NOPD-PCVI are optimized so that its frequency bandgap matches the dominant bridge structural noise frequency range.The noise reduction performance of the NOPD-PCVIs is compared to the steel-spring isolation under different service conditions.展开更多
Effects of plasma equilibrium parameters on the alpha particle loss with the toroidal field ripple based on the CFETR steady-state scenario have been numerically investigated by the orbit-following code GYCAVA. It is ...Effects of plasma equilibrium parameters on the alpha particle loss with the toroidal field ripple based on the CFETR steady-state scenario have been numerically investigated by the orbit-following code GYCAVA. It is found that alpha particle losses decrease and loss regions become narrower with the plasma current increasing or with the magnetic field decreasing. It is because the ripple stochastic transport and the ripple well loss of alpha particle are reduced with the safety factor decreasing. Decrease of the plasma density and temperature can reduce alpha particle losses due to enhancement of the slowing-down effect. The direction of the toroidal magnetic field can significantly affect heat loads induced by lost alpha particle. The vertical asymmetry of heat loads induced by the clockwise and counter-clockwise toroidal magnetic fields are due to the fact that the ripple distribution is asymmetric about the mid-plane, which can be explained by the typical orbits of alpha particle. The maximal heat load of alpha particle for the clockwise toroidal magnetic field is much smaller than that for the counter-clockwise one.展开更多
Numerical simulation is employed to investigate the initial state of avalanche in polydisperse particle systems.Nucleation and propagation processes are illustrated for pentadisperse and triadisperse particle systems,...Numerical simulation is employed to investigate the initial state of avalanche in polydisperse particle systems.Nucleation and propagation processes are illustrated for pentadisperse and triadisperse particle systems,respectively.In these processes,particles involved in the avalanche grow slowly in the early stage and explosively in the later stage,which is clearly different from the continuous and steady growth trend in the monodisperse system.By examining the avalanche propagation,the number growth of particles involved in the avalanche and the slope of the number growth,the initial state can be divided into three stages:T1(nucleation stage),T2(propagation stage),T3(overall avalanche stage).We focus on the characteristics of the avalanche in the T2 stage,and find that propagation distances increase almost linearly in both axial and radial directions in polydisperse systems.We also consider the distribution characteristics of the average coordination number and average velocity for the moving particles.The results support that the polydisperse particle systems are more stable in the T2 stage.展开更多
The effects of high-volume slag-fly ash cement with different particle sizes on hydration degree,microstructure and mechanical properties were systematically studied,by means of laser particle size(DLS),X-ray diffract...The effects of high-volume slag-fly ash cement with different particle sizes on hydration degree,microstructure and mechanical properties were systematically studied,by means of laser particle size(DLS),X-ray diffraction (XRD),comprehensive thermal analysis (TG-DTA),scanning electron microscopy(SEM) and mechanical properties tests.The results show that suitable particle size distribution of cementitious material has significantly promoting effects on hydration reaction rate and mechanical properties.Compared with slag without further grinding,the slag after ball milling for 4 h has an obvious improvement in reactivity,which also provides a faster hydration rate and higher compressive strength for the cementitious material.When the slag milled for 1 and 4 h is mixed at a mass ratio of 2:1 (i e,slag with D_(50) of 7.4μm and average size of 9.9μm,and slag with D_(50) value of 2.6μm and average size of 5.3μm),and a certain amount of fly ash is added in,the most obvious improvement of compressive strength of cement is achieved.展开更多
The macroscopic flow behavior and rheological properties of cemented paste backfill(CPB)are highly impacted by the inherent structure of the paste matrix.In this study,the effects of shear-induced forces and proportio...The macroscopic flow behavior and rheological properties of cemented paste backfill(CPB)are highly impacted by the inherent structure of the paste matrix.In this study,the effects of shear-induced forces and proportioning parameters on the microstructure of fresh CPB were studied.The size evolution and distribution of floc/agglomerate/particles of paste were monitored by focused beam reflection measuring(FBRM)technique,and the influencing factors of aggregation and breakage kinetics of CPB were discussed.The results indicate that influenced by both internal and external factors,the paste kinetics evolution covers the dynamic phase and the stable phase.Increasing the mass content or the cement-tailings ratio can accelerate aggregation kinetics,which is advantageous for the rise of average floc size.Besides,the admixture and high shear can improve breaking kinetics,which is beneficial to reduce the average floc size.The chord length resembles a normal distribution somewhat,with a peak value of approximate 20μm.The particle disaggregation con-stant(k_(2))is positively correlated with the agitation rate,and k_(2) is five orders of magnitude greater than the particle aggregation constant(k1).The kinetics model depicts the evolution law of particles over time quantitatively and provides a theoretical foundation for the micromechanics of complicated rheological behavior of paste.展开更多
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.展开更多
基金funded by the National Natural Science Foundation of China (grant No.11972250)the National Key R&D Program of China (grant Nos.22YFE0207000 and 2022YFC3004505).
文摘Based on a semi-resolved CFD-DEM coupling method,this study proposed a method that uses the minimum distance between the fluid grid and the particle boundary as a reference value to determine the degree of influence of the target fluid grid on the particle's drag force.A fluidized bed of rod-like particles was chosen as a typical case to investigate the effect of different fluid grid scales on various fluidized bed characteristic parameters.The calculation performance of the semi-resolved and unre-solved CFD-DEM coupling algorithm on key fluidized bed characteristic parameters such as average pressure drop,particle frequency distribution with bed height,and particle orientation distribution were compared.It was found that the semi-resolved CFD-DEM coupling algorithm gradually obtained results with higher consistency with decreasing fluid grid scale for key parameters such as particle frequency distribution with bed height,particle orientation distribution,and time-history mixing index,exhibiting a phenomenon similar to grid independence in fluid simulation.By comparing with experimental results,it was verified that the semi-resolved CFD-DEM coupling algorithm can be applied to simulate multi-granular gas-solid systems with fluid grid scales equivalent to particle scales.This algorithm solves the limitation of fluid grid scale in the unresolved CFD-DEM coupling framework and improves the grid adaptability of the CFD-DEM coupling simulation algorithm.
基金This work was supported by the National Key R&D Program‘Transportation Infrastructure’project(No.2022YFB2603400).
文摘The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.
基金support and encouragement of the Joint Funds of the National Natural Science Foundation of China(No.U21B2095)the Major Research Project of National Natural Science Foundation of China(No.91834303).
文摘This work investigated the pyrolysis reaction of waste resin in a fluidized bed reactor.It was found that the pyrolysis-generated ash would adhere to the surface of ceramic particles,causing particle agglomeration and defluidization.Adding kaolin could effectively inhibit the particle agglomeration during the fluidized pyrolysis reaction through physical isolation and chemical reaction.On the one hand,kaolin could form a coating layer on the surface of ceramic particles to prevent the adhesion of organic ash generated by the pyrolysis of resin.On the other hand,when a sufficient amount of kaolin(-0.2%(mass))was added,the activated kaolin could fully contact with the Na+ ions generated by the pyrolysis of resin and react to form a high-melting aluminosilicate mineral(nepheline),which could reduce the formation of low-melting-point sodium sulfate and thereby avoid the agglomeration of ceramic particles.
基金the National Natural Science Foundation of China(Grant No.42271436)the Shandong Provincial Natural Science Foundation,China(Grant Nos.ZR2021MD030,ZR2021QD148).
文摘The existing indoor fusion positioning methods based on Pedestrian Dead Reckoning(PDR)and geomagnetic technology have the problems of large initial position error,low sensor accuracy,and geomagnetic mismatch.In this study,a novel indoor fusion positioning approach based on the improved particle filter algorithm by geomagnetic iterative matching is proposed,where Wi-Fi,PDR,and geomagnetic signals are integrated to improve indoor positioning performances.One important contribution is that geomagnetic iterative matching is firstly proposed based on the particle filter algorithm.During the positioning process,an iterative window and a constraint window are introduced to limit the particle generation range and the geomagnetic matching range respectively.The position is corrected several times based on geomagnetic iterative matching in the location correction stage when the pedestrian movement is detected,which made up for the shortage of only one time of geomagnetic correction in the existing particle filter algorithm.In addition,this study also proposes a real-time step detection algorithm based on multi-threshold constraints to judge whether pedestrians are moving,which satisfies the real-time requirement of our fusion positioning approach.Through experimental verification,the average positioning accuracy of the proposed approach reaches 1.59 m,which improves 33.2%compared with the existing particle filter fusion positioning algorithms.
基金The work described in this paper was partially supported by grants from the Research Grant Council of the Hong Kong Special Administrative Region,China(Project Nos.HKU 17207518 and R5037-18).
文摘The traditional standard wet sieving method uses steel sieves with aperture?0.063 mm and can only determine the particle size distribution(PSD)of gravel and sand in general soil.This paper extends the traditional method and presents an extended wet sieving method.The extended method uses both the steel sieves and the nylon filter cloth sieves.The apertures of the cloth sieves are smaller than 0.063 mm and equal 0.048 mm,0.038 mm,0.014 mm,0.012 mm,0.0063 mm,0.004 mm,0.003 mm,0.002 mm,and 0.001 mm,respectively.The extended method uses five steps to separate the general soil into many material sub-groups of gravel,sand,silt and clay with known particle size ranges.The complete PSD of the general soil is then calculated from the dry masses of the individual material sub-groups.The extended method is demonstrated with a general soil of completely decomposed granite(CDG)in Hong Kong,China.The silt and clay materials with different particle size ranges are further examined,checked and verified using stereomicroscopic observation,physical and chemical property tests.The results further confirm the correctness of the extended wet sieving method.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金supported by the National Magnetic Confinement Fusion Energy R&D Program of China(No.2018YFE0310200)National Natural Science Foundation of China(Nos.11705242,11805138 and 12175156)the Fundamental Research Funds for the Central Universities in China(Nos.YJ201820 and YJ201954)。
文摘An E//B neutral particle analyzer(NPA)has been designed and is under development at Sichuan University and Southwestern Institute of Physics.The main purpose of the E//B NPA is to measure the distribution function of fast ions in the HL-2A/3 tokamak.The E//B NPA contains three main units,i.e.the stripping unit,the analyzing unit and the detection unit.A gas stripping chamber was adopted as the stripping unit.The results of the simulations and beam tests for the stripping chamber are presented.Parallel electric and magnetic fields provided by a NdFeB permanent magnet and two parallel electric plates were designed and constructed for the analyzing unit.The calibration of the magnetic and electric fields was performed using a 50 kV electron cyclotron resonance ion source(ECRIS)platform.The detection unit consists of 32lutetium-yttrium oxyorthosilicate(LYSO)detector modules arranged in two rows.The response functions ofα,hydrogen ions(H^(+),H_(2)^(+)and H_(3)^(+))andγfor a detector module were measured with^(241)Am,^(137)Cs and^(152)Eu sources together with the 50 kV ECRIS platform.The overall results indicate that the designed E//B NPA device is capable of measuring the intensity of neutral hydrogen and deuteron atoms with energy higher than 20 keV.
基金supported by the Major National Science and Technology Project(No.2016ZX05054011)。
文摘Deformable gel particles(DGPs) possess the capability of deep profile control and flooding. However, the deep migration behavior and plugging mechanism along their path remain unclear. Breakage, an inevitable phenomenon during particle migration, significantly impacts the deep plugging effect. Due to the complexity of the process, few studies have been conducted on this subject. In this paper, we conducted DGP flow experiments using a physical model of a multi-point sandpack under various injection rates and particle sizes. Particle size and concentration tests were performed at each measurement point to investigate the transportation behavior of particles in the deep part of the reservoir. The residual resistance coefficient and concentration changes along the porous media were combined to analyze the plugging performance of DGPs. Furthermore, the particle breakage along their path was revealed by analyzing the changes in particle size along the way. A mathematical model of breakage and concentration changes along the path was established. The results showed that the passage after breakage is a significant migration behavior of particles in porous media. The particles were reduced to less than half of their initial size at the front of the porous media. Breakage is an essential reason for the continuous decreases in particle concentration, size, and residual resistance coefficient. However, the particles can remain in porous media after breakage and play a significant role in deep plugging. Higher injection rates or larger particle sizes resulted in faster breakage along the injection direction, higher degrees of breakage, and faster decreases in residual resistance coefficient along the path. These conditions also led to a weaker deep plugging ability. Smaller particles were more evenly retained along the path, but more particles flowed out of the porous media, resulting in a poor deep plugging effect. The particle size is a function of particle size before injection, transport distance, and different injection parameters(injection rate or the diameter ratio of DGP to throat). Likewise, the particle concentration is a function of initial concentration, transport distance, and different injection parameters. These models can be utilized to optimize particle injection parameters, thereby achieving the goal of fine-tuning oil displacement.
基金the National Research Foundation of Korea(Nos.2018R1A5A7023490 and 2022R1A2C1003003)。
文摘A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.
基金supported by the Key Research and Development Plan of Shandong Province(the Major Scientific and Technological Innovation Projects,2021ZDSYS13)the Natural Science Foundation of Shandong Province(ZR2021MB135)Natural Science Foundation of Shandong Province(ZR2021ME224).
文摘It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on bubble motion and gas-liquid mass transfer process in a bubble column.Particle image velocimetry(PIV)technology was employed to analyze the flow field and bubble motion behavior,and dynamic oxygen absorption technology was used to measure the gas-liquid volumetric mass transfer coefficient(kLa).The effect of embedding wire mesh,adding solid particles,and wire mesh coupling solid particles on the flow characteristic and kLa were analyzed and compared.The results show that the gas-liquid interface area increases by 33%-72%when using the wire mesh coupling solid particles strategy compared to the gas-liquid two-phase flow,which is superior to the other two strengthening methods.Compared with the system without reinforcement,kLa in the bubble column increased by 0.5-1.8 times with wire mesh coupling solid particles method,which is higher than the sum of kLa increases with inserting wire mesh and adding particles,and the coupling reinforcement mechanism for affecting gas-liquid mass transfer process was discussed to provide a new idea for enhancing gas-liquid mass transfer.
基金supported by the European Research Council(ERC)under the European Union’s Horizon 2020 research and innovation program(Grant Agreement No.787539)funding from EPRSC(Grant Nos.EP/E035728,EP/C003586,and EP/P010059/1)supported by the National Sciences and Engineering Research Council of Canada(NSERC)and Compute Canada(Job:pve-323-ac,PA).
文摘Realizing the full potential of ultrahigh-intensity lasers for particle and radiation generation will require multi-beam arrangements due to technology limitations.Here,we investigate how to optimize their coupling with solid targets.Experimentally,we show that overlapping two intense lasers in a mirror-like configuration onto a solid with a large preplasma can greatly improve the generation of hot electrons at the target front and ion acceleration at the target backside.The underlying mechanisms are analyzed through multidimensional particle-in-cell simulations,revealing that the self-induced magnetic fields driven by the two laser beams at the target front are susceptible to reconnection,which is one possible mechanism to boost electron energization.In addition,the resistive magnetic field generated during the transport of the hot electrons in the target bulk tends to improve their collimation.Our simulations also indicate that such effects can be further enhanced by overlapping more than two laser beams.
基金Supported by National Key Research and Development Program of China(Grant No.2021YFB2500703)Science and Technology Department Program of Jilin Province of China(Grant No.20230101121JC).
文摘Accurate vehicle dynamic information plays an important role in vehicle driving safety.However,due to the characteristics of high mobility and multiple controllable degrees of freedom of drive-by-wire chassis vehicles,the current mature application of traditional vehicle state estimation algorithms can not meet the requirements of drive-by-wire chassis vehicle state estimation.This paper proposes a state estimation method for drive-by-wire chassis vehicle based on the dual unscented particle filter algorithm,which make full use of the known advantages of the four-wheel drive torque and steer angle parameters of the drive-by-wire chassis vehicle.In the dual unscented particle filter algorithm,two unscented particle filter transfer information to each other,observe the vehicle state information and the tire force parameter information of the four wheels respectively,which reduce the influence of parameter uncertainty and model parameter changes on the estimation accuracy during driving.The performance with the dual unscented particle filter algorithm,which is analyzed in terms of the time-average square error,is superior of the unscented Kalman filter algorithm.The effectiveness of the algorithm is further verified by driving simulator test.In this paper,a vehicle state estimator based on dual unscented particle filter algorithm was proposed for the first time to improve the estimation accuracy of vehicle parameters and states.
文摘Sediment accumulation on the bed of open sewers and drains reduces hydraulic efficiency and can cause localized flooding.Slotted invert traps installed underneath the bed of open sewers and drains can eliminate sediment build-up by catching sediment load.Previous three-dimensional(3D)computational studies have examined the particle trapping performance of invert traps of different shapes and depths under varied sediment and flow conditions,considering particles as spheres.For two-dimensional and 3D numerical modeling,researchers assumed the lid geometry to be a thin line and a plane,respectively.In this 3D numerical study,the particle trapping efficiency of a slotted irregular hexagonal invert trap fitted at the flume bottom was examined by incorporating the particle shape factor of non-spherical sewage solid particles and the thicknesses of upstream and downstream lids over the trap in the discrete phase model of the ANSYS Fluent 2020 R1 software.The volume of fluid(VOF)and the realizable k-turbulence models were used to predict the velocity field.The two-dimensional particle image velocimetry(PIV)was used to measure the velocity field inside the invert trap.The results showed that the thicknesses of upstream and downstream lids affected the velocity field and turbulent kinetic energy at all flow depths.The joint impact of the particle shape factor and lid thickness on the trap efficiency was significant.When both the lid thickness and particle shape factor were considered in the numerical modeling,trap efficiencies were underestimated,with relative errors of-8.66%to-0.65%in comparison to the experimental values of Mohsin and Kaushal(2017).They were also lower than the values predicted by Mohsin and Kaushal(2017),which showed an overall overestimation with errors of-2.3%to 17.4%.
基金Basic Research Funds for Colleges and Universities directly under the Inner Mongolia Autonomous Region:Desert Ecosystem Protection and Restoration Innovation Team(BR 22-13-03).
文摘Wind erosion is a geomorphic process in arid and semi-arid areas and has substantial implications for regional climate and desertification.In the Columbia Plateau of northwestern United States,the emissions from fine particles of loessial soils often contribute to the exceedance of inhalable particulate matter(PM)with an aerodynamic diameter of 10μm or less(PM10)according to the air quality standards.However,little is known about the threshold friction velocity(TFV)for particles of different sizes that comprise these soils.In this study,soil samples of two representative soil types(Warden sandy loam and Ritzville silt loam)collected from the Columbia Plateau were sieved to seven particle size fractions,and an experiment was then conducted to determine the relationship between TFV and particle size fraction.The results revealed that soil particle size significantly affected the initiation of soil movement and TFV;TFV ranged 0.304-0.844 and 0.249-0.739 m/s for different particle size fractions of Ritzville silt loam and Warden sandy loam,respectively.PM10 and total suspended particulates(TSP)emissions from a bed of 63-90μm soil particles were markedly higher for Warden sandy loam than for Ritzville silt loam.Together with the lower TFV of Warden sandy loam,dust emissions from fine particles(<100μm in diameter)of Warden sandy loam thus may be a main contributor to dust in the region's atmosphere,since the PM10 emissions from the soil erosion surfaces and its ensuing suspension within the atmosphere constitute an essential process of soil erosion in the Columbia Plateau.Developing and implementing strategic land management practices on sandy loam soils is therefore necessary to control dust emissions in the Columbia Plateau.
基金Project(51978585)supported by the National Natural Science Foundation,ChinaProject(2022YFB2603404)supported by the National Key Research and Development Program,China+1 种基金Project(U1734207)supported by the High-speed Rail Joint Fund Key Projects of Basic Research,ChinaProject(2023NSFSC1975)supported by the Sichuan Nature and Science Foundation Innovation Research Group Project,China。
文摘The problems associated with vibrations of viaducts and low-frequency structural noise radiation caused by train excitation continue to increase in importance.A new floating-slab track vibration isolator-non-obstructive particle damping-phononic crystal vibration isolator is proposed herein,which uses the particle damping vibration absorption technology and bandgap vibration control theory.The vibration reduction performance of the NOPD-PCVI was analyzed from the perspective of vibration control.The paper explores the structure-borne noise reduction performance of the NOPD-PCVIs installed on different bridge structures under varying service conditions encountered in practical engineering applications.The load transferred to the bridge is obtained from a coupled train-FST-bridge analytical model considering the different structural parameters of bridges.The vibration responses are obtained using the finite element method,while the structural noise radiation is simulated using the frequency-domain boundary element method.Using the particle swarm optimization algorithm,the parameters of the NOPD-PCVI are optimized so that its frequency bandgap matches the dominant bridge structural noise frequency range.The noise reduction performance of the NOPD-PCVIs is compared to the steel-spring isolation under different service conditions.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12175034 and 12005063)the National Key Research and Development Program of China (Grant No.2019YFE03030001)the Fundamental Research Funds for the Central Universities (Grant No.2232022G-10)。
文摘Effects of plasma equilibrium parameters on the alpha particle loss with the toroidal field ripple based on the CFETR steady-state scenario have been numerically investigated by the orbit-following code GYCAVA. It is found that alpha particle losses decrease and loss regions become narrower with the plasma current increasing or with the magnetic field decreasing. It is because the ripple stochastic transport and the ripple well loss of alpha particle are reduced with the safety factor decreasing. Decrease of the plasma density and temperature can reduce alpha particle losses due to enhancement of the slowing-down effect. The direction of the toroidal magnetic field can significantly affect heat loads induced by lost alpha particle. The vertical asymmetry of heat loads induced by the clockwise and counter-clockwise toroidal magnetic fields are due to the fact that the ripple distribution is asymmetric about the mid-plane, which can be explained by the typical orbits of alpha particle. The maximal heat load of alpha particle for the clockwise toroidal magnetic field is much smaller than that for the counter-clockwise one.
基金Project supported by the Qingdao National Laboratory for Marine Science and Technology(Grant No.2015ASKJ01)the National Natural Science Foundation of China(Grant Nos.11972212,12072200,and 12002213).
文摘Numerical simulation is employed to investigate the initial state of avalanche in polydisperse particle systems.Nucleation and propagation processes are illustrated for pentadisperse and triadisperse particle systems,respectively.In these processes,particles involved in the avalanche grow slowly in the early stage and explosively in the later stage,which is clearly different from the continuous and steady growth trend in the monodisperse system.By examining the avalanche propagation,the number growth of particles involved in the avalanche and the slope of the number growth,the initial state can be divided into three stages:T1(nucleation stage),T2(propagation stage),T3(overall avalanche stage).We focus on the characteristics of the avalanche in the T2 stage,and find that propagation distances increase almost linearly in both axial and radial directions in polydisperse systems.We also consider the distribution characteristics of the average coordination number and average velocity for the moving particles.The results support that the polydisperse particle systems are more stable in the T2 stage.
基金Funded by the National Natural Science Foundation of China(No.52172025)。
文摘The effects of high-volume slag-fly ash cement with different particle sizes on hydration degree,microstructure and mechanical properties were systematically studied,by means of laser particle size(DLS),X-ray diffraction (XRD),comprehensive thermal analysis (TG-DTA),scanning electron microscopy(SEM) and mechanical properties tests.The results show that suitable particle size distribution of cementitious material has significantly promoting effects on hydration reaction rate and mechanical properties.Compared with slag without further grinding,the slag after ball milling for 4 h has an obvious improvement in reactivity,which also provides a faster hydration rate and higher compressive strength for the cementitious material.When the slag milled for 1 and 4 h is mixed at a mass ratio of 2:1 (i e,slag with D_(50) of 7.4μm and average size of 9.9μm,and slag with D_(50) value of 2.6μm and average size of 5.3μm),and a certain amount of fly ash is added in,the most obvious improvement of compressive strength of cement is achieved.
基金financially supported by the National Natural Science Foundation of China(No.52104129)the Shandong Provincial Major Science and Technology Innovation Project,China(No.2019SDZY05)+2 种基金the key Laboratory of Mine Ecological Effects and Systematic Restoration,Ministry of Natural Resources(No.MEER-2022-09)the Double First-class Construction Project in Henan Province,China(No.AQ20230735)the Doctoral Fund of Henan Polytechnic University(No.B2021-59).
文摘The macroscopic flow behavior and rheological properties of cemented paste backfill(CPB)are highly impacted by the inherent structure of the paste matrix.In this study,the effects of shear-induced forces and proportioning parameters on the microstructure of fresh CPB were studied.The size evolution and distribution of floc/agglomerate/particles of paste were monitored by focused beam reflection measuring(FBRM)technique,and the influencing factors of aggregation and breakage kinetics of CPB were discussed.The results indicate that influenced by both internal and external factors,the paste kinetics evolution covers the dynamic phase and the stable phase.Increasing the mass content or the cement-tailings ratio can accelerate aggregation kinetics,which is advantageous for the rise of average floc size.Besides,the admixture and high shear can improve breaking kinetics,which is beneficial to reduce the average floc size.The chord length resembles a normal distribution somewhat,with a peak value of approximate 20μm.The particle disaggregation con-stant(k_(2))is positively correlated with the agitation rate,and k_(2) is five orders of magnitude greater than the particle aggregation constant(k1).The kinetics model depicts the evolution law of particles over time quantitatively and provides a theoretical foundation for the micromechanics of complicated rheological behavior of paste.
基金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.