Mathematical modeling and simulations of stress mitigation by coating polycrystalline particles in lithium-ion batteries

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.

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.

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.

Study of deep transportation and plugging performance of deformable gel particles in porous media

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.

Modeling analysis of cobalt-based Fischer-Tropsch catalyst particles

The influences of particle size,shape,and catalyst distribution on the reactivity and hydrocarbon product selectivity of a cobalt-based catalyst for Fischer-Tropsch synthesis were investigated in the present work.A self-consistent kinetic model for Fischer-Tropsch reaction proposed here was found to correlate experimental data well and hence was used to describe the consumption rates of reactants and formation rates of hydrocarbon products.The perturbed-chain statistical associating fluid theory equation of state was used to describe vapor-liquid equilibrium behavior associated with Fischer-Tropsch reaction.Local interaction between intraparticle diffusion and Fischer-Tropsch reaction was investigated in detail.Results showed that in order to avoid the adverse influence of intraparticle diffusional limitations on catalyst reactivity and product selectivity,the use of small particles is necessary.Large eggshell spherical particles are shown to keep the original catalyst reactivity and enhance the selectivity of heavy hydrocarbon products.The suitable layer thickness for a spherical particle with a diameter of 2 mm is nearly 0.15 mm.With the same outer diameter of 2 mm,the catalyst reactivity and heavy product selectivity of hollow cylindrical particles with a layer thickness of 0.25 mm are found to be larger than eggshell spherical particles.From the viewpoint of catalytic performance,hollow cylindrical particles are a better choice for industrial applications.

Influence of syngas components and ash particles on the radiative heat transfer in a radiant syngas cooler

Radiant syngas cooler(RSC)is widely used as a waste heat recovery equipment in industrial gasification.In this work,an RSC with radiation screens is established and the impact of gaseous radiative property models,gas components,and ash particles on heat transfer is investigated by the numerical simulation method.Considering the syngas components and the pressure environment of the RSC,a modified weighted-sum-of-gray-gases model was developed.The modified model shows high accuracy in validation.In computational fluid dynamics simulation,the calculated steam production is only 0.63%in error with the industrial data.Compared with Smith's model,the temperature decay along the axial direction calculated by the modified model is faster.Syngas components are of great significance to heat recovery capacity,especially when the absorbing gas fraction is less than 10%.After considering the influence of particles,the outlet temperature and the proportion of radiative heat transfer are less affected,but the difference in steam output reaches 2.7 t·h^(-1).The particle deposition on the wall greatly reduces the heat recovery performance of an RSC.

Particle size spatial distribution in landslide dams

The particle composition and spatial distribution of landslide-induced dam bodies are critical geotechnical parameters for studying the hazards of dam-break floods.However,current research often neglects the influence of the initial particle composition and spatial distribution of the landslide on the particle composition and spatial distribution of the landslide dam.This study investigated the impact of initial particle size distribution,volume,and sliding length on the energy and velocity changes of characteristic particles during the sliding process and the spatial distribution of particle sizes in the landslide dam body.Numerical simulations and physical models were employed to examine the effects of sequential gradient arrangements(where particle sizes decrease from top to bottom)and four other different initial particle arrangements on the energy and velocity changes of particles and the spatial distribution of particle sizes in the dam body.The study reveals the characteristics of translational and rotational energy of different particles and the laws of mechanical energy conversion,obtaining the spatial distribution patterns of particle sizes in landslide-induced dams.The results show that under the sequential gradient arrangement,the energy dissipation of the landslide movement is lower,with larger particles mainly distributed at the distal end and smaller particles at the proximal end of the landslide dam.In contrast,under the reverse gradient arrangement,the energy dissipation of the landslide movement is higher,and the distribution pattern of the dam particles is opposite to that of the sequential gradient arrangement.For the other arrangement modes,the spatial distribution of dam particles falls between the aforementioned two.There is a positive correlation between particle size and translational kinetic energy within the particle flow during the landslide process,and rotational motion increases energy dissipation.Under constant slope conditions,sliding length does not affect the movement pattern of the particle flow or the spatial distribution of particles in the dam body.The findings of this study provide a scientific basis for the accurate simulation and prediction of dam-break flood processes.

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.

Cascade refrigeration system synthesis based on hybrid simulated annealing and particle swarm optimization algorithm

Cascade refrigeration system(CRS)can meet a wider range of refrigeration temperature requirements and is more energy efficient than single-refrigerant refrigeration system,making it more widely used in low-temperature industry processes.The synthesis of a CRS with simultaneous consideration of heat integration between refrigerant and process streams is challenging but promising for significant cost saving and reduction of carbon emission.This study presented a stochastic optimization method for the synthesis of CRS.An MINLP model was formulated based on the superstructure developed for the CRS,and an optimization framework was proposed,where simulated annealing algorithm was used to evolve the numbers of pressure/temperature levels for all sub-refrigeration systems,and particle swarm optimization algorithm was employed to optimize the continuous variables.The effectiveness of the proposed methodology was verified by a case study of CRS optimization in an ethylene plant with 21.89%the total annual cost saving.

Location and Capacity Determination Method of Electric Vehicle Charging Station Based on Simulated Annealing Immune Particle Swarm Optimization

As the number of electric vehicles(EVs)continues to grow and the demand for charging infrastructure is also increasing,how to improve the charging infrastructure has become a bottleneck restricting the development of EVs.In other words,reasonably planning the location and capacity of charging stations is important for development of the EV industry and the safe and stable operation of the power system.Considering the construction and maintenance of the charging station,the distribution network loss of the charging station,and the economic loss on the user side of the EV,this paper takes the node and capacity of charging station planning as control variables and the minimum cost of system comprehensive planning as objective function,and thus proposes a location and capacity planning model for the EV charging station.Based on the problems of low efficiency and insufficient global optimization ability of the current algorithm,the simulated annealing immune particle swarm optimization algorithm(SA-IPSO)is adopted in this paper.The simulated annealing algorithm is used in the global update of the particle swarm optimization(PSO),and the immune mechanism is introduced to participate in the iterative update of the particles,so as to improve the speed and efficiency of PSO.Voronoi diagram is used to divide service area of the charging station,and a joint solution process of Voronoi diagram and SA-IPSO is proposed.By example analysis,the results show that the optimal solution corresponding to the optimisation method proposed in this paper has a low overall cost,while the average charging waiting time is only 1.8 min and the charging pile utilisation rate is 75.5%.The simulation comparison verifies that the improved algorithm improves the operational efficiency by 18.1%and basically does not fall into local convergence.

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 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）.

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.

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.