Three-dimensional Simulation of Gas/Solid Flow in Spout-fluid Beds with Kinetic Theory of Granular Flow

A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of lm and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Besides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macroscopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.

A transient production prediction method for tight condensate gas wells with multiphase flow

Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and pressure in the full-path of tight condensate gas well is proposed,and a model for predicting the transient production from tight condensate gas wells with multiphase flow is established.The research indicates that the relationship curve between condensate oil saturation and pressure is crucial for calculating the pseudo-pressure.In the early stage of production or in areas far from the wellbore with high reservoir pressure,the condensate oil saturation can be calculated using early-stage production dynamic data through material balance models.In the late stage of production or in areas close to the wellbore with low reservoir pressure,the condensate oil saturation can be calculated using the data of constant composition expansion test.In the middle stages of production or when reservoir pressure is at an intermediate level,the data obtained from the previous two stages can be interpolated to form a complete full-path relationship curve between oil saturation and pressure.Through simulation and field application,the new method is verified to be reliable and practical.It can be applied for prediction of middle-stage and late-stage production of tight condensate gas wells and assessment of single-well recoverable reserves.

CFD modeling of gas−liquid flow phenomenon in lead smelting oxygen-enriched side-blown furnace

A validated numerical model was established to simulate gas−liquid flow behaviors in the oxygen-enriched side-blown bath furnace.This model included the slip velocity between phases and the gas thermal expansion effect.Its modeling results were verified with theoretical correlations and experiments,and the nozzle-eroded states in practice were also involved in the analysis.Through comparison,it is confirmed that the thermal expansion effect influences the flow pattern significantly,which may lead to the backward motion of airflow and create a potential risk to production safety.Consequently,the influences of air injection velocity and furnace width on airflow behavior were investigated to provide operating and design guidance.It is found that the thin layer melt,which avoids high-rate oxygen airflow eroding nozzles,shrinks as the injection velocity increases,but safety can be guaranteed when the velocity ranges from 175 to 275 m/s.Moreover,the isoline patterns and heights of thin layers change slightly when the furnace width increases from 2.2 to 2.8 m,indicating that the furnace width shows a limited influence on production safety.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Model Prediction and Optimal Control of Gas Oxygen Content for A Municipal Solid Waste Incineration Process

In the municipal solid waste incineration process,it is difficult to effectively control the gas oxygen content by setting the air flow according to artificial experience.To address this problem,this paper proposes an optimization control method of gas oxygen content based on model predictive control.First,a stochastic configuration network is utilized to establish a prediction model of gas oxygen content.Second,an improved differential evolution algorithm that is based on parameter adaptive and t-distribution strategy is employed to address the set value of air flow.Finally,model predictive control is combined with the event triggering strategy to reduce the amount of computation and the controller's frequent actions.The experimental results show that the optimization control method proposed in this paper obtains a smaller degree of fluctuation in the air flow set value,which can ensure the tracking control performance of the gas oxygen content while reducing the amount of calculation.

A developed transient gas-liquid-solid flow model with hydrate phase transition for solid fluidization exploitation of marine natural gas hydrate reservoirs

The multiphase flow characteristic is one of the most concerning problems during solid fluidization exploitation of marine natural gas hydrate reservoirs.In this research,a new transient gas-liquid-solid multiphase flow model with hydrate phase transition was developed.Meanwhile,this model considered the coupling relationship among convective heat transfer,hydrate dynamic decomposition,and multi-phase flow.The model can simulate the change of flow pattern from solid-liquid to gas-liquid-solid flow,and describe the distribution character of volume fraction of phase,wellbore temperature and pressure,and hydrate decomposition rate during transportation.The simulation results indicate that the hydrate decomposition region in the wellbore gradually expands,but the hydrate decomposition rate gradually decreases during the solid fluidization exploitation of hydrate.When mining time lasts for 4 h,and the bottom hole pressure decreases by about 0.4 MPa.Increasing NaCl concentration in seawater helps expand hydrate decomposition regions and improves the wellbore hydrate decomposition rate.When the Nacl mass fraction in seawater reaches 15%,it will raise the hydrate decomposition regions to the whole wellbore.In addition,the higher the wellhead backpressure,the lower the decomposition area and decomposition rate of hydrate in the wellbore.When wellhead backpressure reaches 2 MPa,the volume fraction of gas near the wellhead will reduce to about 12%.This work is expected to provide a theoretical basis for the development of marine hydrate reservoirs.

Criterion of gas and solid dual-phase flow atomization crash in molten metal

A self-invented atomization process, in which molten metal is atomized into powder by a high-velocity gas stream carrying solid particles as the atomization medium, was introduced. The characteristics of powders prepared by common gas atomization and dual-phase flow atomization under similar conditions were compared. The experimental results show that the dual-phase flow-atomized powders have average particle sizes that are one-half that of the common gas-atomized particles;additionally, they possess a finer microstructure and higher cooling rate under the same atomization gas pressure and the same gas flow. The Weber number in the crash criteria of liquid atomization is adopted to measure the crash ability of the atomization media. The Weber number of the dual-phase flow atomization medium is the sum of that of the gas and the solid particles. Furthermore, the critical equation of the crash model in dual-phase flow atomization is established, and the main regularities associated with this process were analyzed.

Effect of gas blowing nozzle angle on multiphase flow and mass transfer during RH refining process

A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.

Experimental study on solid particle migration and production behaviors during marine natural gas hydrate dissociation by depressurization

Sand production is one of the main obstacles restricting gas extraction efficiency and safety from marine natural gas hydrate(NGH)reservoirs.Particle migration within the NGH reservoir dominates sand production behaviors,while their relationships were rarely reported,severely constrains quantitative evaluation of sand production risks.This paper reports the optical observations of solid particle migration and production from micrometer to mesoscopic scales conditioned to gravel packing during depressurization-induced NGH dissociation for the first time.Theoretical evolutionary modes of sand migration are established based on experimental observations,and its implications on field NGH are comprehensively discussed.Five particle migration regimes of local borehole failure,continuous collapse,wormhole expansion,extensive slow deformation,and pore-wall fluidization are proved to occur during depressurization.The types of particle migration regimes and their transmission modes during depressurization are predominantly determined by initial hydrate saturation.In contrast,the depressurization mainly dominates the transmission rate of the particle migration regimes.Furthermore,both the cumulative mass and the medium grain size of the produced sand decrease linearly with increasing initial methane hydrate(MH)saturation.Discontinuous gas bubble emission,expansion,and explosion during MH dissociation delay sand migration into the wellbore.At the same time,continuous water flow is a requirement for sand production during hydrate dissociation by depressurization.The experiments enlighten us that a constitutive model that can illustrate visible particle migration regimes and their transmission modes is urgently needed to bridge numerical simulation and field applications.Optimizing wellbore layout positions or special reservoir treatment shall be important for mitigating sand production tendency during NGH exploitation.

Simultaneous measurement of velocity profile and liquid film thickness in horizontal gas–liquid slug flow by using ultrasonic Doppler method

Horizontal gas-liquid two-phase flows widely exist in chemical engineering,oil/gas production and other important industrial processes.Slug flow pattern is the main form of horizontal gas-liquid flows and characterized by intermittent motion of film region and slug region.This work aims to develop the ultrasonic Doppler method to realize the simultaneous measurement of the velocity profile and liquid film thickness of slug flow.A single-frequency single-channel transducer is adopted in the design of the field-programmable gate array based ultrasonic Doppler system.A multiple echo repetition technology is used to improve the temporal-spatial resolution for the velocity profile.An experiment of horizontal gas-liquid two-phase flow is implemented in an acrylic pipe with an inner diameter of 20 mm.Considering the aerated characteristics of the liquid slug,slug flow is divided into low-aerated slug flow,high-aerated slug flow and pseudo slug flow.The temporal-spatial velocity distributions of the three kinds of slug flows are reconstructed by using the ultrasonic velocity profile measurement.The evolution characteristics of the average velocity profile in slug flows are investigated.A novel method is proposed to derive the liquid film thickness based on the instantaneous velocity profile.The liquid film thickness can be effectively measured by detecting the position and the size of the bubbles nearly below the elongated gas bubble.Compared with the time of flight method,the film thickness measured by the Doppler system shows a higher accuracy as a bubble layer occurs in the film region.The effect of the gas distribution on the film thickness is uncovered in three kinds of slug flows.

Numerical simulation of gas–liquid flow in the bubble column using Wray–Agarwal turbulence model coupled with population balance model

In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..

Inter-layer interference for multi-layered tight gas reservoir in the absence and presence of movable water

Due to the dissimilarity among different producing layers,the influences of inter-layer interference on the production performance of a multi-layer gas reservoir are possible.However,systematic studies of inter-layer interference for tight gas reservoirs are really limited,especially for those reservoirs in the presence of water.In this work,five types of possible inter-layer interferences,including both absence and presence of water,are identified for commingled production of tight gas reservoirs.Subsequently,a series of reservoir-scale and pore-scale numerical simulations are conducted to quantify the degree of influence of each type of interference.Consistent field evidence from the Yan'an tight gas reservoir(Ordos Basin,China)is found to support the simulation results.Additionally,suggestions are proposed to mitigate the potential inter-layer interferences.The results indicate that,in the absence of water,commingled production is favorable in two situations:when there is a difference in physical properties and when there is a difference in the pressure system of each layer.For reservoirs with a multi-pressure system,the backflow phenomenon,which significantly influences the production performance,only occurs under extreme conditions(such as very low production rates or well shut-in periods).When water is introduced into the multi-layer system,inter-layer interference becomes nearly inevitable.Perforating both the gas-rich layer and water-rich layer for commingled production is not desirable,as it can trigger water invasion from the water-rich layer into the gas-rich layer.The gas-rich layer might also be interfered with by water from the neighboring unperforated water-rich layer,where the water might break the barrier(eg weak joint surface,cement in fractures)between the two layers and migrate into the gas-rich layer.Additionally,the gas-rich layer could possibly be interfered with by water that accumulates at the bottom of the wellbore due to gravitational differentiation during shut-in operations.

A new numerical approach of coupled modeling for solid deformation and gas leak flow in multi-coal-seams

From the viewpoint of interaction mechanics for solid and gas, a coupled mathematical model was presented for solid coal/rock deformation and gas leak flow in parallel deformable coal seams. Numerical solutions using the SIP (Strong Implicit Proce- dure) method to the coupled mathematical model for double parallel coal seams were also developed in detail. Numerical simulations for the prediction of the safety range using protection layer mining were performed with experimental data from a mine with potential danger of coal/gas outbursts. Analyses show that the numerical simulation results are consistent with the measured data in situ.

Non-symmetric distributions of solids deposition for solid-water stratified flow in deviated tubing strings

Dynamic performance on solids flow with water in deviated tubing is essential for the reliability of pump and normal operation of horizontal and directional wells.Compared with coal-water flow in vertical tubing and sand-oil flow with high production in deviated tubing,solids deposition with water shows obvious non-symmetric distributions in deviated tubing from simulations and experiments.The mathematical model of two phase flow was developed under coupling conditions of deviated tubing,low flow rate and viscosity based on the kinetic theory of granular flow and first-order discrete scheme.The results show that solid-water stratified flow in deviated tubing can be divided into two zones of suspension bed and the moving bed throughout the flow field.The solid flowing velocity with water is negative and particles slide down at the bottom of moving bed zone.The process of solids flow with water in deviated tubing will produce pressure loss and consume the kinetic energy.The thickness of deposited layer and the flowing velocity of solids flow downward with water at the moving bed zone enhance with the decreased inlet flow rate and the increased particle size,tubing inside diameter(ID)and inclination angle.Solids are easier into suspension from the upper part of moving bed zone to suspension bed zone and more solid particles flow with water towards the tubing outlet with the increase of inlet flowing velocity.The decision is made to reduce the screen width,tubing ID and inclination angle to maintain to be greater than critical deposition velocity in order to prevent solids settling.And it provides the theoretical basis and technical reserves for solid control and offers an effective approach to enhance tubing cleaning in deviated strings.

Numerical Simulation of Oil and Gas Two-Phase Flow in Deep Condensate Gas Reservoirs in Bohai Buried Hills

The BZ19-6 gas field is characterized by high temperature and high pressure (HTHP), high condensate content, little difference between the formation pressure and dew point pressure, and large amount of reverse condensate liquid. During the early stage of depletion development, the production gas-oil ratio (GOR) and production capacity remain relatively stable, which is inconsistent with the conventional reverse condensate seepage law. In view of the static and dynamic conflict in development and production, indoor high-temperature and high-pressure PVT experiment was carried out to reveal the mist-like condensation phenomenon of fluids in the BZ19-6 formation. And the seepage characteristics of condensate gas reservoirs with various degrees of depletion under the condition of HTHP were analyzed based on production performance. The change rule of fluid phase state was analyzed in response to the characterization difficulties of the seepage mechanism. The fluid state was described using the miscible mechanism. And the interphase permeability interpolation coefficient was introduced based on interfacial tension. By doing so, the accurate characterization of the “single-phase flow of condensate gas-near-miscible mist-like quasi single-phase flow-oil-gas two-phase flow” during the development process was achieved. Then the accurate fitting of key indicators for oilfield development was completed, and the distribution law of formation pressure and the law of condensate oil precipitation under different reservoir conditions are obtained. Based on research results, the regulation strategy of variable flow rate production was developed. Currently, the work system has been optimized for 11 wells, achieving a “zero increase” in the GOS of the gas field and an annual oil increase of 22,000 cubic meters.

Numerical Simulation of Gas-Solid Flow in Square Cyclone Separators with Downward Exit

Cyclone separators are widely used in industrial applications. The separation efficiency and pressure drop are the most important parameters to evaluate the performance of processing system. In the simulations,the flow behavior of gas and particles within a square cyclone separator is simulated by means of computational fluid dynamics. The RNG k- ε model and the Reynolds stress model( RSM) are used to model gas turbulence. The flow behavior is examined in the term of tangential velocity components,static pressure and pressure drop contour plots for flow field and solid volume fraction. The effects of the turbulence model and solid volume fraction on the square cyclone are discussed. The results indicate that the pressure drop increases with the increase of solid volume fraction,and increase with the increase of inlet velocities for two turbulence models, moreover,the simulations results are compared with pressure field. For all runs,the RSM model gives a higher pressure drop compared to the RNG k- ε model. The RSM model provides well the forced vortex and free vortex,and captures better the phenomena occurring during intense vortex flow in the presence of walls within cyclone separators.

Permeability evolution and gas flow in wet coal under non-equilibrium state:Considering both water swelling and process-based gas swelling

Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.

Simulation of Gas-Water Two-Phase Flow in Tight Gas Reservoirs Considering the Gas Slip Effect

A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.

Effect of gas flow on the nanoparticles transport in dusty acetylene plasmas

This article presents simulation results on the effects of neutral gas flow for nanoparticle transport in atmospheric-pressure,radio-frequency,capacitively-coupled,and acetylene discharge.The acetylene gas is set to flow into the chamber from the upper showerhead electrode.The internal energy of the gas medium therein is transferred into kinetic energy so the gas advection can be triggered.This is represented by the pressure volume work term of the gas energy converse equation.The gas advection leads to the gas temperature sink at the gas inlet,hence a large gas temperature gradient is formed.The thermophoresis relies on the gas temperature gradient,and causes the profile of nanoparticle density to vary from a double-peak structure to a single-peak one.The gas advection influences the properties of electron density and temperature as well and causes the drift-ambipolar mode profile of electron density asymmetric.In the bulk region,i.e.away from the inlet,the gas advection is more like one isovolumetric compression,which slightly increases the temperature of the gas medium at consuming its kinetic energy.

Gas liquid cylindrical cyclone flow regime identification using machine learning combined with experimental mechanism explanation

The flow regimes of GLCC with horizon inlet and a vertical pipe are investigated in experiments,and the velocities and pressure drops data labeled by the corresponding flow regimes are collected.Combined with the flow regimes data of other GLCC positions from other literatures in existence,the gas and liquid superficial velocities and pressure drops are used as the input of the machine learning algorithms respectively which are applied to identify the flow regimes.The choosing of input data types takes the availability of data for practical industry fields into consideration,and the twelve machine learning algorithms are chosen from the classical and popular algorithms in the area of classification,including the typical ensemble models,SVM,KNN,Bayesian Model and MLP.The results of flow regimes identification show that gas and liquid superficial velocities are the ideal type of input data for the flow regimes identification by machine learning.Most of the ensemble models can identify the flow regimes of GLCC by gas and liquid velocities with the accuracy of 0.99 and more.For the pressure drops as the input of each algorithm,it is not the suitable as gas and liquid velocities,and only XGBoost and Bagging Tree can identify the GLCC flow regimes accurately.The success and confusion of each algorithm are analyzed and explained based on the experimental phenomena of flow regimes evolution processes,the flow regimes map,and the principles of algorithms.The applicability and feasibility of each algorithm according to different types of data for GLCC flow regimes identification are proposed.