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

Simulation of two-phase debris flow scouring bridge pier

Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Numerical simulation of three-dimensional combustion flows

A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.

Numerical simulation of slurry jets using mixture model

Slurry jets in a static uniform environment were simulated with a two-phase mixture model in which flow-particle interactions were considered. A standard k-e turbulence model was chosen to close the governing equations. The computational results were in agreement with previous laboratory measurements. The characteristics of the two-phase flow field and the influences of hydraulic and geometric parameters on the distribution of the slurry jets were analyzed on the basis of the computational results. The calculated results reveal that if the initial velocity of the slurry jet is high, the jet spreads less in the radial direction. When the slurry jet is less influenced by the ambient fluid （when the Stokes number St is relatively large）, the turbulent kinetic energy k and turbulent dissipation rate e, which are relatively concentrated around the jet axis, decrease more rapidly after the slurry jet passes through the nozzle. For different values of St, the radial distributions of streamwise velocity and particle volume fraction are both self-similar and fit a Gaussian profile after the slurry jet fully develops. The decay rate of the particle velocity is lower than that of water velocity along the jet axis, and the axial distributions of the centerline particle streamwise velocity are self-similar along the jet axis. The pattern of particle dispersion depends on the Stokes number St. When St = 0.39, the panicle dispersion along the radial direction is considerable, and the relative velocity is very low due to the low dynamic response time. When St = 3.08, the dispersion of particles along the radial direction is very little, and most of the particles have high relative velocities along the streamwise direction.

NUMERICAL SIMULATION FOR AIR AND AIR-PM FLOW IN WALL FLOW DIESEL PARTICULATE FILTERS

Numerical simulations are performed both for the single airflow and air-PMtwo-phase flow in wall flow diesel participate filters (DPF) for the first time. The calculationdomain is divided into two regions. In. the inlet and outlet flow channels, the simulations areperformed for the steady and laminar flow; In the porous filtration walls, the calculation model forflow in porous media is used. The Lagrange two-phase flow model is used to calculate the air-PMflow in DPF, for the dispersed phase (PM), its flow tracks are obtained by the integrating of theLagrange kinetic equation. The calculated velocity, pressure distribution and PM flow tracks in DPFare obtained, which exhibits the main flow characteristics in wall flow DPF and will be help for theoptimal design and performance prediction of wall flow DPF.

Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.

NUMERICAL SIMULATION ON 2-D WATER-AIR TWO-PHASE FLOW OVER TOP OUTLET

Flood discharge over top outlet of dam is simulated by 2 dimension water air two phase mathematical model. Distribution of dynamic pressure, turbulent kinetic energy (k), turbulent dissipation rate (ε), free water surface and velocity field have been obtained. The simulated results were tested by physical model, which shows that the computed results are identical with that of the physical model.

Numerical simulation of liquid desiccant evaporator driven by heat pump

The standard k-ε turbulence model and discrete phase model (DPM) were used to simulate the heat and mass transfer in a liquid-desiccant evaporator driven by a heat pump using FLUENT software, and the temperature field and velocity field in the device were obtained. The performance of the liquid-desiccant evaporator was studied as the concentration of the inlet solution varied between 21% and 30% and the pipe wall temperature between 30 and 50 ℃. Results show that the humidification rate and the humidification efficiency increased with the inlet air temperature, the solution flow rate, the solution temperature, and the pipe wall temperature. The humidification rate and humidification efficiency decreased with increasing moisture content in inlet air and the concentration of inlet solution. The humidification rate increased substantially but the humidification efficiency decreased as the inlet air flow rate increased. The error between the simulations and experimental results is acceptable, meaning that our model can provide a theoretical basis for optimizing the performance of a humidifying evaporator.

Experimental and numerical investigation of liquid-solid binary fluidized beds： Radioactive particle tracking technique and dense discrete phase model simulations

Liquid-solid binary fluidized beds are widely used in many industries. However, the flow behavior of such beds is not well understood due to the lack of accurate experimental and numerical data. In the current study, the behavior of monodisperse and binary liquid-solid fluidized beds of the same density but dif- ferent sizes is investigated using radioactive particle tracking （RPT） technique and a dense discrete phase model （DDPM）. Experiments and simulations are performed in monodisperse fluidized beds containing two different sizes of glass beads （0.6 and I mm） and a binary fluidized bed of the same particles for vari- ous bed compositions. The results show that both RPT and DDPM can predict the mixing and segregation pattern in liquid-solid binary fluidized beds. The mean velocity predictions of DDPM are in good agree- ment with the experimental findings for both monodisperse and binary fluidized beds. However, the axial root mean square velocity predictions are only reasonable for bigger particles. Particle-particle interac- tions are found to be critical for predicting the flow behavior of solids in liquid-solid binary fluidized beds.

Modeling and simulation of chemically reacting flows in gas-solid catalytic and non-catalytic processes

This paper gives an overview of the recent development of modeling and simulation of chemically react- ing flows in gas-solid catalytic and non-catalytic processes. General methodology has been focused on the Eulerian-Lagrangian description of particulate flows, where the particles behave as the catalysts or the reactant materials. For the strong interaction between the transport phenomena （i.e., momentum, heat and mass transfer） and the chemical reactions at the particle scale, a cross-scale modeling approach, i.e., CFD-DEM or CFD-DPM, is established for describing a wide variety of complex reacting flows in multiphase reactors, Representative processes, including fluid catalytic cracking （FCC）, catalytic conversion of syngas to methane, and coal pyrolysis to acetylene in thermal plasma, are chosen as case studies to demonstrate the unique advantages of the theoretical scheme based on the integrated particle-scale information with clear physical meanings, This type of modeling approach provides a solid basis for understanding the multiphase reacting flow problems in general.

Multi-Scale Numerical Simulation of Flow,Heat and Mass Transfer Behaviors in Dense Gas-Solid Flows:A Brief Review

Dense gas-solid flows are very common in actual production and industrial fields,so it is significant to understand their hydrodynamic characteristics and heat and mass transfer behaviors.This article provides a brief review of multi-scale numerical simulation of flow,heat and mass transfer behaviors in dense gas-solid flows.It describes multiscale models(direct numerical simulation,discrete particle model,and two-fluid model)and the results of related research.Finally,it discusses possible future developments in research on the flow,heat and mass transfer characteristics of dense gas-solid two-phase flows.

Fluid Flow Behavior of Liquid in Cylindrical Vessels Stirred by One or Two Air Jets

Based on the two-phase model (Eulerian-Eulerian model), the three dimensional fluid flow in water and that liquid steel systems stirred by one or two multiple gas jets are simulated. In the Eulerian-Eulerian two-phase model, the gas and the liquid phase are considered to be two different continuous fluids interacting with each other through the finite inter-phase areas. The exchange between the phases is represented by source terms in conversation equations. Turbulence is assumed to be a property of the liquid phase. A new turbulence modification - model is introduced to consider the bubbles movement contribution to and . The dispersion of phases due to turbulence is represented by introducing a diffusion term in mass conservation equation. The mathematical simulation agrees well with the experiment results. The study results indicate that the distance of two nozzles has big effect on fluid flow behavior in the vessel. Using two gas injection nozzles at the half radii of one diameter of the bottom generates a much better mixing than with one nozzle under the condition of the same total gas flow rate.

Scheduling of high-speed rail traffic based on discrete-time movement model

In this paper, a new simulation approach for solving the mixed train scheduling problem on the high-speed double-track rail line is presented. Based on the discrete-time movement model, we propose control strategies for mixed train movement with different speeds on a high-speed double-track rail line, including braking strategy, priority rule, travelling strategy, and departing rule. A new detailed algorithm is also presented based on the proposed control strategies for mixed train movement. Moreover, we analyze the dynamic properties of rail traffic flow on a high-speed rail line. Using our proposed method, we can effectively simulate the mixed train schedule on a rail line. The numerical results demonstrate that an appropriate decrease of the departure interval can enhance the capacity, and a suitable increase of the distance between two adjacent stations can enhance the average speed. Meanwhile, the capacity and the average speed will be increased by appropriately enhancing the ratio of faster train number to slower train number from 1.

Numerical simulation of fluid flow and alloy melting in RH process for electrical steels

Based on the Eulerian-Lagrangian approach,a mathematical model was established to describe the gas-liquid flow behavior in the Ruhrstahl-Heraeus(RH)degasser.The momentum source and the turbulent kinetic energy source due to the motion of gas bubbles were considered for the liquid flow.The effect of the expansion of gas bubbles on the liquid velocity,recirculation rate,and mixing time was quantitatively evaluated.After the fluid flow reached the steady state,the melting and mixing processes of aluminum alloys in the RH degasser were also investigated.The results indicate that the expansion of gas bubbles has a significant infuence on the recirculation rate and the mixing time in the RH process.Increasing the superheat of liquid steel and decreasing the initial diameter of alloy particles are beneficial to promote the melting and mixing of alloy particles.Due to the existence of solidified steel shells,the maximum diameter of the alloy particle is about 1.5 times its initial diameter.

A Review of the Dynamic Modeling Approaches for Characterizing Fluid Flow in Naturally Fractured Reservoirs

Fluid flow in fractured media has been studied for decades and received considerable attention in the oil and gas industry because of the high productivity of naturally fractured reservoirs.Due to formation complexity and reservoir heterogeneity,characterizing fluid flow with an appropriate reservoir model presents a challenging task that differs relatively from homogeneous conventional reservoirs in many aspects of view,including geological,petrophysical,production,and economics.In most fractured reservoirs,fracture networks create complex pathways that affect hydrocarbon flow,well performance,hence reservoir characterization.A better and comprehensive understanding of the available reservoir modeling approaches is much needed to accurately characterize fluid flow behavior in NFRs.Therefore,in this paper,a perspective review of the available modeling approaches was presented for fluid flow characterization in naturally fractured medium.Modeling methods were evaluated in terms of their description,application,advantages,and disadvantages.This study has also included the applications of these reservoir models in fluid flow characterizing studies and governing equations for fluid flow.Dual continuum models were proved to be better than single continuum models in the presence of large scale fractures.In comparison,discrete models were more appropriate for reservoirs that contain a smaller number of fractures.However,hybrid modeling was the best method to provide accurate and scalable fluid flow modeling.It is our understanding that this paper will bridge the gap between the fundamental understanding and application of NFRs modeling approaches and serve as a useful reference for engineers and researchers for present and future applications.

Numerical simulation of spray performance based on the Euler-Lagrange Approach

Numerical simulations have been carried out to investigate the liquid atomization and spray process using the Discrete Phase Model of the commercial CFD code combined with the Wall-Film boundary conditions. The effects of spray parameters on droplets Santer mean diameter （SMD）, droplet collision speed, the thickness of liquid-film, the surface temperature and its uniformity were analyzed in the present study. The simulation results and the experimental data obtained in the available literature agree within 13.8%, The computational results show that the spray pressure is the main factor to realize the atomization. Increasing the mass flux and the spray pressure, the droplet collision speed increases while the corresponding maximum film thickness on the heated surface declines. The surface temperature changes indistinctively with the increase of the spray distance, but the temperature distribution tends to be uniform.

Numerical investigation of fluid phase momentum transfer in carbonate acidizing

This work mainly studies the effect of fluid phase momentum transfer mechanisms on the acidizing results,including the retardation effect of the porous structure and the interaction between the fluid phase,such as viscous dissipation and inertial effect.The results show that the acid fluid momentum transfer is influenced by the complex porous structure and fluid viscous dissipation.Eventually,the Stokes-Darcy equation is recommended to be adopted to describe the fluid phase momentum transfer in the following numerical simulation studies of the carbonate acidizing process.Based on this model,a parametric research is carried out to investigate the impact of acid on rock physical characteristics in the stimulation process.Increasing the acid concentration appears to minimize the quantity of acid consumed for the breakthrough.The acid surface reaction rate has a considerable impact on the pore volume to breakthrough and the optimum acid injection rate.The influence of permeability on the acidizing results basically shows a negative correlation with the injection rate.The difference between the acidizing curves of different permeability gradually becomes insignificant with the decrease in injection rate.The existence of isolated fracture and vug significantly reduces acid consumption for the breakthrough.

MODELING TWO-PHASE FLOW IN PULSED FLUIDIZED BED

Mathematical models for pulsed fluidization are systematically discussed. Several undetermined constitutive relationships are included in the General Two-Fluid Model (GTFM), the adjustable parameters of which are always chosen at will to some extent. Although there are no adjustable parameters in the Basic Two-Fluid Model (BTFM), its eigenvalues are complex numbers and it is ill-posed for initial-value problems. The Local Equilibrium Model (LEM), a further simplification of BTFM, is discussed at length. Although the model is very simple, it is highly capable of simulating complex processes in pulsed fluidization over a broad range of operating parameters, and its numerical results well fit experimental results in both the variation of bed height and the distribution of particle concentration as fluidizing velocity varies.