Sample size adaptive strategy for time-dependent Monte Carlo particle transport simulation

When multiphysics coupling calculations contain time-dependent Monte Carlo particle transport simulations, these simulations often account for the largest part of the calculation time, which is insufferable in certain important cases. This study proposes an adaptive strategy for automatically adjusting the sample size to fulfil more reasonable simulations. This is realized based on an extension of the Shannon entropy concept and is essentially different from the popular methods in timeindependent Monte Carlo particle transport simulations, such as controlling the sample size according to the relative error of a target tally or by experience. The results of the two models show that this strategy can yield almost similar results while significantly reducing the calculation time. Considering the efficiency, the sample size should not be increased blindly if the efficiency cannot be enhanced further. The strategy proposed herein satisfies this requirement.

Numerical simulation of predicting and reducing solid particle erosion of solid-liquid two-phase flow in a choke

Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.

Validation of the current and pressure coupling schemes with nonlinear simulations of TAE and analysis on the linear stability of tearing mode in the presence of energetic particles

Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under different approximations.First,when considering only the perturbed distribution function of energetic particles(EPs),the equivalence can be proved analytically.Second,when both the variations of the magnetic field and the EP distribution function are included,the current and pressure coupling schemes numerically produce the same result in the nonlinear simulations.On this basis,the influences of co-/counter-passing and trapped EPs on the linear stabilities of the m/n=2/1 tearing mode(TM)have been investigated(where m and n represent the poloidal and toroidal mode numbers,respectively).The results of scanningβh of EPs show that the co-passing and trapped EPs are found to stabilize the TM,while the counter-passing EPs tend to destabilize the TM.The behind(de)stabilization mechanisms of the TM by EPs are carefully analyzed.Furthermore,after exceeding critical EP betas,the same branch of the high-frequency mode is excited by co-/counterpassing and trapped EPs,which is identified as the m/n=2/1 energetic particle mode.

Dynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device

The shear-induced migration of neutrally-buoyant non-colloidal circular particles in a two-dimensional circular Couette flow is investigated numerically with a distributed Lagrange multiplier based fictitious domain method.The effects of inertia and volume fraction on the particle migration are examined.The results indicate that inertia has a negative effect on the particle migration.In consistence with the experimental observations,the rapid migration of particles near the inner cylinder at the early stage is observed in the simulation,which is believed to be related to the chain-like clustering of particles.The migration of circular particles in a plane Poiseuille flow is also examined in order to further confirm the effect of such clustering on the particle migration at early stage.There is tendency for the particles in the vicinity of outer cylinder in the Couette device to pack into concentric rings at late stage in case of high particle concentration.

Application of DSAPSO Algorithm in Distribution Network Reconfiguration with Distributed Generation

With the current integration of distributed energy resources into the grid,the structure of distribution networks is becoming more complex.This complexity significantly expands the solution space in the optimization process for network reconstruction using intelligent algorithms.Consequently,traditional intelligent algorithms frequently encounter insufficient search accuracy and become trapped in local optima.To tackle this issue,a more advanced particle swarm optimization algorithm is proposed.To address the varying emphases at different stages of the optimization process,a dynamic strategy is implemented to regulate the social and self-learning factors.The Metropolis criterion is introduced into the simulated annealing algorithm to occasionally accept suboptimal solutions,thereby mitigating premature convergence in the population optimization process.The inertia weight is adjusted using the logistic mapping technique to maintain a balance between the algorithm’s global and local search abilities.The incorporation of the Pareto principle involves the consideration of network losses and voltage deviations as objective functions.A fuzzy membership function is employed for selecting the results.Simulation analysis is carried out on the restructuring of the distribution network,using the IEEE-33 node system and the IEEE-69 node system as examples,in conjunction with the integration of distributed energy resources.The findings demonstrate that,in comparison to other intelligent optimization algorithms,the proposed enhanced algorithm demonstrates a shorter convergence time and effectively reduces active power losses within the network.Furthermore,it enhances the amplitude of node voltages,thereby improving the stability of distribution network operations and power supply quality.Additionally,the algorithm exhibits a high level of generality and applicability.

Electromagnetic Particle-in-Cell Simulations of Electron Holes Formed During the Electron Two-Stream Instability

Previous electrostatic particle-in-cell （PIC） simulations have pointed out that elec- tron phase-space holes （electron holes） can be formed during the nonlinear evolution of the electron two-stream instability. The parallel cuts of the parallel and perpendicular electric field have bipolar and unipolar structures in these electron holes, respectively. In this study, two-dimensional （2D） electromagnetic PIC simulations are performed in the x - y plane to investigate the evolution of the electron two-stream instability, with the emphasis on the magnetic structures associated with these electron holes in different plasma conditions. In the simulations, the background magnetic field （Bo = Boer） is along the x direction. In weakly magnetized plasma （Ωe 〈ωpe, where Ωe and ωpe are the electron gyrofrequency and electron plasma frequency, respectively）, several 2D electron holes are formed. In these 2D electron holes, the parallel cut of the fluctuating magnetic field δBx and δBz has unipolar structures, while the fluctuating magnetic field δBy has bipolar structures. In strongly magnetized plasma （Ωe 〉 ωpe）, several quasi-lD electron holes are formed. The electrostatic whistler waves with streaked structures of Ey are excited. The fluctuating mag- netic field δBx and δBz also have streaked structures. The fluctuating magnetic field δBx and δBy are produced by the current in the z direction due to the electric field drift of the trapped elec- trons, while the fluctuating magnetic field δBz can be explained by the Lorentz transformation of a moving quasielectrostatic structure. The influences of the initial temperature anisotropy on the magnetic structures of the electron holes are also analyzed. The electromagnetic whistler waves are found to be excited in weakly magnetized plasma. However, they do not have any significant effects on the electrostatic structures of the electron holes.

Measurement and simulation of the two-phase velocity correlation in sudden-expansion gas-particle flows

In this paper the present authors measured the gas-particle two-phase velocity correlation in sudden expansion gas-particle flows with a phase Doppler particle anemometer （PDPA） and simulated the system behavior by using both a Reynolds-averaged Navier-Stokes （RANS） model and a large-eddy simulation （LES）. The results of the measurements yield the axial and radial time-averaged velocities as well as the fluctuation velocities of gas and three particle-size groups （30μm, 50μm, and 95μm） and the gasparticle velocity correlation for 30μm and 50μm particles. From the measurements, theoretical analysis, and simulation, it is found that the two-phase velocity correlation of sudden-expansion flows, like that of jet flows, is less than the gas and particle Reynolds stresses. What distinguishes the two-phase velocity correlations of sudden-expansion flow from those of jet and channel flows is the absence of a clear relationship between the two-phase velocity correlation and particle size in sudden-expansion flows. The measurements, theoretical analysis, and numerical simulation all lead to the above-stated conclusions. Quantitatively, the results of the LES are better than those of the RANS model.

Refinement of Adaptive Dynamical Simulation of Quantum Mechanical Double Slit Interference Phenomenon

We applied adaptive dynamics to double slit interference phenomenon using particle model and obtained partial successful results in our previous report. The patterns qualitatively corresponded well with experiments. Several properties such as concave single slit pattern and large influence of slight displacement of the emission position were different from the experimental results. In this study we tried other slit conditions and obtained consistent patterns with experiments. We do not claim that the adaptive dynamics is the principle of quantum mechanics, but the present results support the probability of adaptive dynamics as the candidate of the basis of quantum mechanics. We discuss the advantages of the adaptive dynamical view for foundations of quantum mechanics.

Numerical Simulation and Analysis of Solid-liquid Two-phase Flow in Centrifugal Pump

The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.

Numerical simulation of the sedimentation of cylindrical pollutant particles in fluid

The sedimentation of cylindrical pollutant particles which fall through a fluid is investigated. Differing from previous research work, particle oscillation and effect of particle on the fluid are considered, and the torque exerted on a particle when viscous fluid flow around a particle is got through experiment and included in the numerical simulation. The computational results showed that the sedimentation velocities of particle increase slowly with the increase of particle aspect ratio . For disk like particle, when the motion direction of particle is parallel to axis of particle, particle falls more slowly than the case of perpendicular to axis of particle; while for rod like particle, it is inverse. For sedimentation of a crowd of high frequency oscillating cylindrical particles with arbitrary initial orientation, both vertical velocity and horizontal velocity oscillate dramatically, the degree of oscillation of the former is stronger than the later. A crowd of particles fall more quickly than an isolated particle. Particles tend to strongly align in the direction of gravity. The computational results agreed well with the experimental ones and helpful for controlling of pollutant particles.

Numerical Simulation of Particle Concentration in a Gas Cyclone Separator

The particle concentration inside a cyclone separator at different operation parameters was simulated with the FLUENT software. The Advanced Reynolds Stress Model （ARSM） was used in gas phase turbulence modeling. Stochastic Particle Tracking Model （SPTM） and the Particle-Source-In-Cell （PSIC） method were adopted for particles computing. The interaction between particles and the gas phase was also taken into account. The numerical simulation results were in agreement with the experimental data. The simulation revealed that an unsteady spiral dust strand appeared near the cyclone wall and a non-axi-symmetrical dust ring appeared in the annular space and under the cover plate of the cyclone. There were two regions in the radial particle concentration distribution, in which particle concentration was low in the inner region （r/R≤0.75） and increased greatly in the outer region （r/R〉0.75）. Large particles generally had higher concentration in the near-wall region and small particles had higher concentration in the inner swirling flow region. The axial distribution of particle concentration in the inner swirling flow （r/R≤0.3） region showed that there existed serious fine particle entrainment within the height of 0.SD above the dust discharge port and a short-cut flow at a distance of about 0.25D below the entrance of the vortex finder. The dimensionless concentration in the high-concentration region increased obviously in the upper part of the cyclone separation space when inlet particle loading was large. With increasing gas temperature, the particle separation ability of the cyclone was obviously weakened.

Physical Simulation and Experimental Examination of ε-Cu Particles Dissolution Evolution During Welding of Copper Precipitation Strengthening Steel

The kinetics of ε-Cu particles dissolution in the matrix during welding of a copper-precipitation strengthening steel was determined by a combination of GleebleTM physical simulation, TEM examination and hardness meas urement. The ε-Cu particles underwent a coarsening and part dissolution and then complete dissolution reaction as the peak temperature increased from 750 to 1 000℃, which resulted in the decrease in the number density of ε-Cu particles and hardness in the heat affected zone （HAZ）. The results can be used to understand the evolution of this transformation and a softening behavior of the HAZ during welding of this type of steel.

Numerical simulation of profile control by clay particles after polymer flooding

A three-dimensional,two-phase,five-component mathematical model has been developed to describe flow characteristics of clay particles and flocs in the profile control process,in which the clay particle suspension is injected into the formation to react with residual polymer.This model considers the reaction of clay particles with residual polymer,apparent viscosity of the mixture,retention of clay particles and flocs,as well as the decline in porosity and permeability caused by the retention of clay particles and flocs.A finite difference method is used to discretize the equation for each component in the model.The Runge-Kutta method is used to solve the polymer flow equation,and operator splitting algorithms are used to split the flow equation for clay particles into a hyperbolic equation for convection and a parabolic equation for diffusion,which effectively ensures excellent precision,high speed and good stability.The numerical simulation had been applied successfully in the 4-P1920 unit of the Lamadian Oilfield to forecast the blocking capacity of clay particle suspension and to optimize the injection parameters.

Particle dispersion modeling in ventilated room using artificial neural network

Due to insufficiency of a platform based on experimental results for numerical simulation validation using computational fluid dynamic method(CFD) for different geometries and conditions,in this paper we propose a modeling approach based on the artificial neural network(ANN) to describe spatial distribution of the particles concentration in an indoor environment.This study was performed for a stationary flow regime.The database used to build the ANN model was deducted from bibliography literature and composed by 261 points of experimental measurement.Multilayer perceptron-type neural network(MLP-ANN) model was developed to map the relation between the input variables and the outputs.Several training algorithms were tested to give a choice of the Fletcher conjugate gradient algorithm(TrainCgf).The predictive ability of the results determined by simulation of the ANN model was compared with the results simulated by the CFD approach.The developed neural network was beneficial and easy to predict the particle dispersion curves compared to CFD model.The average absolute error given by the ANN model does not reach 5%against 18%by the Lagrangian model and 28% by the Euler drift-flux model of the CFD approach.

DEM simulation of particle flow on a single deck banana screen

A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.

NUMERICAL SIMULATION OF THE TUBE EROSION RESULTED FROM PARTICLES IMPACTS

In this paper,a tube erosion caused by the turbulent flow of a dilute particle-laden gas isstudied numerically.Eulerian equations are used to describe the gas-phase motion with the turbulenceviscosity evaluated from a k-ε model of turbulence.The effect of the turbulence with a stochasticdispersion model has been taken into account for the prediction of impact particle velocity and itstrajectory,The particle impact and rebound model and the erosion model of ductile alloys obtainedby Tabakoff et al.are used to predict the particle rebound phenomena and the erosion suffered bythe tubes.The results obtained in this study include the distributions of particle collision frequencyand erosion of tube surface.

Numerical simulation of the blocking process of gelled particles in porous media with remaining polymers

Gelled particles can be transferred deeply inside oil reservoirs to block water channels due to their physicochemical characteristics, including swelling, deformation, and synergetic effect （reacting with polymers）, and then the injection profiles are significantly modified. At present, research on gelled particles is mainly focused on laboratory studies of drive mechanisms, and rarely on mathematical models describing the blocking process of gelled particles. In this paper, the blocking process of gelled particles is divided into two sub-processes： deposition and desorption due to particle deformation. A mathematical model based on filtration theory is proposed considering the effect of characteristics of gelled particles on the blocking process. Blocking laws were simulated and researched using the mathematical model. Results of the simulation of the blocking of gelled particles are quite consistent with the experimental results, which confirms the reliability of the mathematical model developed.

NUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL CYCLONE

The three-dimension gas-particle flow in a spiral cyclone is simulated nu- merically in this paper. The gas flow field was obtained by solving the three-dimension Navier-Stokes equations with Reynolds Stress Model （RSM）. It is shown that there are two regions in the cyclone, the steadily tangential flow in the spiral channel and the combined vortex flow in the centre. Numerical results for particles trajectories show that the initial position of the particle at the inlet plane substantially affects its trajectory in the cyclone. The particle collection efficiency curves at different inlet velocities were obtained and the effects of inlet flow rate On the performance of the spiral cyclone were presented. Numerical results also show that the increase of flow rate leads to the increase of particles collection efficiency, but the pressure drop increases sharply.

Pyrolysis of single large biomass particle: Simulation and experiments

Pyrolysis and heat transfer characteristics of single large biomass particle were investigated using threedimensional unsteady heat transfer model coupled with chemical reactions.The consumption of biomass and the production of products were simulated.Some experiments were designed to provide model parameters for simulation calculations.The simulation was verified by pyrolysis experiments of large biomass particle in a vertical tube furnace.The simulation results show the internal heat and mass transfer law during the pyrolysis of large biomass particle.When the biomass particle diameter is between 10 and 30 mm,for every 5 mm increase in particle diameter,the time required for complete pyrolysis will increase on average by about 50 s.When the pyrolysis temperature is between 673 K and 873 K,a slight decrease in the pyrolysis temperature will cause the time required for the biomass to fully pyrolyze to rise significantly.And the phenomenon is more obvious in the low temperature range.The results indicate that the numerical simulation agrees well with the experimental results.

Particle path tracking method in two- and three-dimensional continuously rotating detonation engines

The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.