Large Eddy Simulation of Particle Wake Effect and RANS Simulation of Turbulence Modulation in Gas-Particle Flows

The turbulence enhancement by particle wake effect is studied by large eddy simulation （LES） of turbulent gas flows passing a single particle. The predicted time-averaged and root-mean-square fluctuation velocities behind the particle are in agreement with the Reynolds-averaged Navier-Stokes modeling results and experimental results. A semi-empirical turbulence enhancement model is proposed by the present-authors based on the LES resuits. This model is incorporated into the second-order moment two-phase turbulence model for simulating vertical gas-particle pipe flows and horizontal gas-particle channel flows. The simulation results show that compared with the model not accounting for the particle wake effect, the present model gives simulation results for the gas turbulence modulation in much better agreement with the experimental results.

SIMULATION OF SUDDEN-EXPANSION AND SWIRLING GAS-PARTICLE FLOWS USING A TWO-FLUID PARTICLE-WALL COLLISION MODEL WITH CONSIDERATION OF THE WALL ROUGHNESS

A two-fluid particle-wall collision model with consideration of wall roughness is pro- posed.It takes into account the effects of the friction,restitution and in particular the wall roughness, and hence the redistribution of Reynolds stress in different directions,the absorption of turbulent en- ergy from the mean motion and the attenuation of particle motion by the wall.The proposed model is used to simulate sudden-expansion and swirling gas-particle flows and is validated by comparing with experimental results.The results show that the proposed model gives better results than those obtained by the presently used zero-gradient condition.Hence,it is suggested that the proposed model should be used as the wall boundary condition for the particle phase in place of the presently used boundary condition.

Effect of Particle Size Distribution on Radiative Heat Transfer in High-Temperature Homogeneous Gas-Particle Mixtures

The weighted-sum-of-gray-gas(WSGG)model and Mie theory are applied to study the influents of particle size on the radiative transfer in high temperature homogeneous gas-particle mixtures,such as the flame in aero-engine combustor.The radiative transfer equation is solved by the finite volume method.The particle size is assumed to obey uniform distribution and logarithmic normal(L-N)distribution,respectively.Results reveal that when particle size obeys uniform distribution,increasing particle size with total particle volume fraction fvunchanged will result in the decreasing of the absolute value of radiative heat transfer properties,and the effect of ignoring particle scattering will also be weakened.Opposite conclusions can be obtained when total particle number concentration N0 is unchanged.Moreover,if particle size obeys L-N distribution,increasing the narrowness indexσor decreasing the characteristic diameter Dˉwith the total particle volume fraction fvunchanged will increase the absolute value of radiative heat transfer properties.With total particle number concentration N0 unchanged,opposite conclusions for radiative heat source and incident radiation terms can be obtained except for radiative heat flux term.As a whole,the effects of particle size on the radiative heat transfer in the high-temperature homogeneous gas-particle mixtures are complicated,and the particle scattering cannot be ignoring just according to the particle size.

Advances in LES of Two-phase Combustion(Ⅱ) LES of Complex Gas-Particle Flows and Coal Combustion

Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the sin- gle-scale second-order moment two-phase turbulence model superposed with a particle collision model （USM-θ model） in most regions.

Experimental study and numerical simulation of gas-particle flows with radial bias combustion and centrally fuel rich swirl burners

Numerical simulation is applied to gas-particle flows of the primary and the secondary air ducts and burner region, and of two kinds of swirl burners. The modeling results of Radial Bias Combustion (RBC) burner well agreed with the data from the three-dimensional Phase-Doppler anemometry (PDA) experiment by Li, et al. The modeling test conducted in a 1025 t/h boiler was to study the quality of aerodynamics for a Central Fuel Rich (CFR) burner, and the Internal Recirculation Zone (IRZ) was measured. In addition, gas-particle flows with a CFR burner were investigated by numerical simulation, whose results accorded with the test data fundamentally. By analyzing the distribution of gas velocity and trajectories of particles respectively, it is found that the primary air’s rigidity of CFR burner is stronger than that of RBC burner, and the primary air mixes with the secondary air later. Furthermore, high concentration region of pulverized coal exists in the burner’s central zone whose atmosphere is reduced, and trajectories of particles in IRZ of CFR burner are longer than that of RBC burner. They are favorable to coal’s ignition and the reduction of NOx emission.

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.

Distributed Nash Equilibrium Seeking Strategies Under Quantized Communication

This paper is concerned with distributed Nash equi librium seeking strategies under quantized communication. In the proposed seeking strategy, a projection operator is synthesized with a gradient search method to achieve the optimization o players' objective functions while restricting their actions within required non-empty, convex and compact domains. In addition, a leader-following consensus protocol, in which quantized informa tion flows are utilized, is employed for information sharing among players. More specifically, logarithmic quantizers and uniform quantizers are investigated under both undirected and connected communication graphs and strongly connected digraphs, respec tively. Through Lyapunov stability analysis, it is shown that play ers' actions can be steered to a neighborhood of the Nash equilib rium with logarithmic and uniform quantizers, and the quanti fied convergence error depends on the parameter of the quan tizer for both undirected and directed cases. A numerical exam ple is given to verify the theoretical results.

THE FLOW STRUCTURE OF DILUTE GAS-PARTICLE SUSPENSIONS BEHIND A SHOCK WAVE MOVING ALONG A FLAT SURFACE

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented. The Saffman lift force acting on a particle in a shear flow is taken into account. It is shown that particle migration across the boundary layer leads to intersections of particle trajectories. The corresponding modification of dusty gas model is proposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at a constant speed are obtained by using the method of matched asymptotic expansions. The method of the calculation of particle phase parameters in Lagrangian coordinates is described in detail. Some numerical results for the case of small particle concentration are given.

Equilibrium Strategy of the Pursuit-Evasion Game in Three-Dimensional Space

The pursuit-evasion game models the strategic interaction among players, attracting attention in many realistic scenarios, such as missile guidance, unmanned aerial vehicles, and target defense. Existing studies mainly concentrate on the cooperative pursuit of multiple players in two-dimensional pursuit-evasion games. However, these approaches can hardly be applied to practical situations where players usually move in three-dimensional space with a three-degree-of-freedom control. In this paper,we make the first attempt to investigate the equilibrium strategy of the realistic pursuit-evasion game, in which the pursuer follows a three-degree-of-freedom control, and the evader moves freely. First, we describe the pursuer's three-degree-of-freedom control and the evader's relative coordinate. We then rigorously derive the equilibrium strategy by solving the retrogressive path equation according to the Hamilton-Jacobi-Bellman-Isaacs(HJBI) method, which divides the pursuit-evasion process into the navigation and acceleration phases. Besides, we analyze the maximum allowable speed for the pursuer to capture the evader successfully and provide the strategy with which the evader can escape when the pursuer's speed exceeds the threshold. We further conduct comparison tests with various unilateral deviations to verify that the proposed strategy forms a Nash equilibrium.

Magnetopause properties at the dusk magnetospheric flank from global magnetohydrodynamic simulations,the kinetic Vlasov equilibrium,and in situ observations--Potential implications for SMILE

We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.

Numerical analysis for the free-boundary current reversal equilibrium in the AC plasma current operation in a tokamak

In recent decades, tokamak discharges with zero total toroidal current have been reported in tokamak experiments, and this is one of the key problems in alternating current(AC) operations.An efficient free-boundary equilibrium code is developed to investigate such advanced tokamak discharges with current reversal equilibrium configuration. The calculation results show that the reversal current equilibrium can maintain finite pressure and also has considerable effects on the position of the X-point and the magnetic separatrix shape, and hence also on the position of the strike point on the divertor plates, which is extremely useful for magnetic design, MHD stability analysis, and experimental data analysis etc. for the AC plasma current operation on tokamaks.

Kinetic equilibrium reconstruction with internal safety factor profile constraints on EAST tokamak

Reconstruction of plasma equilibrium plays an important role in the analysis and simulation of plasma experiments. The kinetic equilibrium reconstruction with pressure and edge current constraints has been employed on EAST tokamak. However, the internal safety factor(q) profile is not accurate. This paper proposes a new way of incorporating q profile constraints into kinetic equilibrium reconstruction. The q profile is yielded from the Polarimeter Interferometer(POINT)reconstruction. Virtual probes containing information on q profile constraints are added to inputs of the kinetic equilibrium reconstruction program to obtain the final equilibrium. The new equilibrium produces a more accurate internal q profile. This improved method would help analyze EAST experiments.

Effects of ionic liquids on the vapor-liquid equilibrium of 1,3,5-trioxane-water system at 101.3 kPa

Increasing the 1,3,5-trioxane(TOX) concentration in the equilibrated vapor phase of TOX-H_(2)O system has been recognized as a challenge for the azeotrope. Ionic liquids(ILs) were used to improve the relative volatility of TOX to H_(2)O and destroy the azeotrope in the TOX-H2O system. The vapor-liquid equilibrium of TOX-H2O system at 101.3 kPa was studied with the addition of 1-butyl-3-methylimidazolium hydrogen sulfate, 1-hexyl-3-methylimidazolium hydrogen sulfate and 1-butyl-3-methylimidazolium nitrate, respectively. The results showed that the volatility of TOX increased with the increase in IL dosage. And the volatility of water decreased with the increase in IL dosage. The relative volatility of TOX to H_(2)O was improved with the increase in ILs dosage. The azeotrope could be destroyed with an IL mole fraction of about 0.10. A non-random two-liquid(NRTL) model was successfully used to correlate the experimental data. The interaction parameters were obtained by fitting the experimental data with the model. The results indicated that a strong interaction existed between ILs and water. The strong interaction improved the volatility of TOX and inhibited the volatility of water, and then intensified the relative volatility of TOX to H_(2)O. The results showed that an ILs with strong polarity and hydrophilicity may be a potential additive to improve the TOX concentration in the equilibrated vapor phase.

Phase equilibrium data prediction and process optimizationin butadiene extraction process

In response to the lack of reliable physical parameters in the process simulation of the butadiene extraction,a large amount of phase equilibrium data were collected in the context of the actual process of butadiene production by acetonitrile.The accuracy of five prediction methods,UNIFAC(UNIQUAC Functional-group Activity Coefficients),UNIFAC-LL,UNIFAC-LBY,UNIFAC-DMD and COSMO-RS,applied to the butadiene extraction process was verified using partial phase equilibrium data.The results showed that the UNIFAC-DMD method had the highest accuracy in predicting phase equilibrium data for the missing system.COSMO-RS-predicted multiple systems showed good accuracy,and a large number of missing phase equilibrium data were estimated using the UNIFAC-DMD method and COSMO-RS method.The predicted phase equilibrium data were checked for consistency.The NRTL-RK(non-Random Two Liquid-Redlich-Kwong Equation of State)and UNIQUAC thermodynamic models were used to correlate the phase equilibrium data.Industrial device simulations were used to verify the accuracy of the thermodynamic model applied to the butadiene extraction process.The simulation results showed that the average deviations of the simulated results using the correlated thermodynamic model from the actual values were less than 2%compared to that using the commercial simulation software,Aspen Plus and its database.The average deviation was much smaller than that of the simulations using the Aspen Plus database(>10%),indicating that the obtained phase equilibrium data are highly accurate and reliable.The best phase equilibrium data and thermodynamic model parameters for butadiene extraction are provided.This improves the accuracy and reliability of the design,optimization and control of the process,and provides a basis and guarantee for developing a more environmentally friendly and economical butadiene extraction process.

Multi-Objective Equilibrium Optimizer for Feature Selection in High-Dimensional English Speech Emotion Recognition

Speech emotion recognition(SER)uses acoustic analysis to find features for emotion recognition and examines variations in voice that are caused by emotions.The number of features acquired with acoustic analysis is extremely high,so we introduce a hybrid filter-wrapper feature selection algorithm based on an improved equilibrium optimizer for constructing an emotion recognition system.The proposed algorithm implements multi-objective emotion recognition with the minimum number of selected features and maximum accuracy.First,we use the information gain and Fisher Score to sort the features extracted from signals.Then,we employ a multi-objective ranking method to evaluate these features and assign different importance to them.Features with high rankings have a large probability of being selected.Finally,we propose a repair strategy to address the problem of duplicate solutions in multi-objective feature selection,which can improve the diversity of solutions and avoid falling into local traps.Using random forest and K-nearest neighbor classifiers,four English speech emotion datasets are employed to test the proposed algorithm(MBEO)as well as other multi-objective emotion identification techniques.The results illustrate that it performs well in inverted generational distance,hypervolume,Pareto solutions,and execution time,and MBEO is appropriate for high-dimensional English SER.

Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

Intrinsic correlation between the generalized phase equilibrium condition and mechanical behaviors in hydrate-bearing sediments

The phase equilibrium and mechanical behaviors of natural gas hydrate-bearing sediment are essential for gas recovery from hydrate reservoirs.In heating closed systems,the temperature-pressure path of hydrate-bearing sediment deviates from that of pure bulk hydrate,reflecting the porous media effect in phase equilibrium.A generalized phase equilibrium equation was established for hydrate-bearing sediments,which indicates that both capillary and osmotic pressures cause the phase equilibrium curve to shift leftward on the temperature-pressure plane.In contrast to bulk hydrate,hydrate-bearing sediment always contains a certain amount of unhydrated water,which keeps phase equilibrium with the hydrate within the hydrate stability field.With changes in temperature and pressure,a portion of pore hydrate and unhydrated water may transform into each other,affecting the shear strength of hydrate-bearing sediment.A shear strength model is proposed to consider not only hydrate saturation but also the change in temperature and pressure of hydrate-bearing sediment.The model is validated by experimental data with various hydrate saturation,temperature and pressure conditions.The deformation induced by partial dissociation was studied through depressurization tests under constant effective stress.The reduction in gas pressure within the hydrate stability field indeed caused sediment deformation.The dissociation-induced deformation can be reasonably estimated as the difference in volume between hydrate-bearing and hydrate-free sediments from the compression curves.

Magnetic diagnostics layout design for CFETR plasma equilibrium reconstruction

Plasma equilibrium reconstruction provides essential information for tokamak operation and physical analysis.An extensive and reliable set of magnetic diagnostics is required to obtain accurate plasma equilibrium.This study designs and optimizes the magnetic diagnostics layout for the reconstruction of the equilibrium of the plasma according to the scientific objectives,engineering design parameters,and limitations of the Chinese Fusion Engineering Test Reactor(CFETR).Based on the CFETR discharge simulation,magnetic measurement data are employed to reconstruct consistent plasma equilibrium parameters,and magnetic diagnostics'number and position are optimized by truncated Singular value decomposition,verifying the redundancy reliability of the magnetic diagnostics layout design.This provides a design solution for the layout of the magnetic diagnostics system required to control the plasma equilibrium of CFETR,and the developed design and optimization method can provide effective support to design magnetic diagnostics systems for future magnetic confinement fusion devices.

Probabilistic Analysis of Slope Using Finite Element Approach and Limit Equilibrium Approach around Amalpata Landslide of West Central, Nepal

The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.