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.

Simulation of the effect of hydrate adhesion properties on flow safety in solid fluidization exploitation

During the solid fluidization exploitation of marine natural gas hydrates,the hydrate particles and cuttings produced via excavation and crushing are transported by the drilling mud.The potential flow safety issues arising during the transport process,such as the blockage of pipelines and equipment,have attracted considerable attention.This study aims to investigate the impact of hydrate adhesion features,including agglomeration,cohesion,and deposition,on the flow transport processes in solid fluidization exploitation and to provide a reference for the design and application of multiphase hydrate slurry transport in solid fluidization exploitation.We established a numerical simulation model that considers the hydrate adhesion properties using the coupled computational fluid dynamics and discrete element method(CFD-DEM)for the multiphase mixed transport in solid fluidization exploitation.An appropriate model to simulate the adhesion force of the hydrate particles and the corresponding parameter values were obtained.The conclusions obtained are as follows.Under the same operating conditions,a stationary bed is more likely to form in the transport process due to the hydrate adhesion forces;adhesion forces can increase the critical deposition velocity of the mixture of hydrate particles and cuttings.Hydrate adhesion lowers the height of the solid-phase moving bed,while the agglomeration and cohesion of particles can intensify the aggregation and deposition of hydrate debris and cuttings at the bottom of the pipe.These particles tend to form a deposit bed rather than a moving bed,which reduces the effective flow area of the pipeline and increases the risk of blockage.

Fluidization characteristics of different sizes of quartz particles in the fluidized bed

Fluidization characteristics of quartz particles with different sizes are experimentally investigated in a fluidized bed with an inner diameter of 300 mm and height of 8250 mm.Results show that the average solid holdup increases with the increase in superficial gas velocity and the decrease in initial solid holdup in the dense zone of the fluidized bed.The average cross-sectional solid holdup decreases with increasing bed height and superficial gas velocity.The bed expansion coefficient increases with the increase in superficial gas velocity and the decrease in solid holdup.Correlations of average solid holdup,average cross-sectional solid holdup and bed expansion coefficient are also established and discussed.These correlations can provide guidelines for better understanding of the fluidization characteristics.

Oxide coating mechanism during fluidized bed reduction: solid-state reaction characteristics between iron ore particles and MgO

Experiments on the solid-state reaction between iron ore particles and MgO were performed to investigate the coating mechanism of MgO on the iron ore particles＇ surface during fluidized bed reduction. MgO powders and iron ore particles were mixed and compressed into briquettes and, subsequently, roasted at different temperatures and for different time periods. A Mg-containing layer was observed on the outer edge of the iron ore particles when the roasting temperature was greater than 1173 K. The concentration of Fe in the Mg-containing layer was evenly distributed and was approximately 10wt%, regardless of the temperature change. Boundary layers of Mg and Fe were observed outside of the iron ore particles. The change in concentration of Fe in the boundary layers was simulated using a gas–solid diffusion model, and the diffusion coefficients of Fe and Mg in these layers at different temperatures were calculated. The diffusion activation energies of Fe and Mg in the boundary layers in these experiments were evaluated to be approximately 176 and 172 k J/mol, respectively.

Predicting minimum fluidization velocities of multi-component solid mixtures

Employing well-established mixing rules for mean properties, appropriate expressions are derived for predicting minimum fluidization velocities of multi-component solid mixtures in terms of mono- component values for the velocity and the bed voidage at incipient fluidization. Based on flow regime and the mixing level of constituent species, it is found that these relationships differ significantly from each other, whether related to size-different or density-different mixtures. For mixed beds of size-different mixtures, the effect of volume contraction is accounted for by the mean voidage term, which is absent for segregated beds. Incorporating the volume-change of mixing leads to values of the mixture minimum fluidization velocities even lower than corresponding values for segregated bed, thus conforming to the trend reported in the literature. Size-different mixtures exhibit flow regime dependence irrespective of whether the bed is mixed or segregated. On the other hand, the mixing of constituent species does not affect the minimum fiuidization velocity of density-different mixtures, as the difference in the expres- sions for a segregated and a mixed system is rather inconsequential. Comparison with experimental data available in the literature is made to test the efficacy of the minimum fluidization velocity expressions derived here.

Effect of volume-contraction on incipient fluidization of binary-solid mixtures

Towards the development of a predictive model for computing the minimum fluidization velocity, the volume-contraction phenomenon arising from the mixing of unequal solid species is accounted for in the prediction of the bed void fraction of binary-solid mixtures at the incipient fluidization conditions. Com- parison with experimental data obtained from the literature clearly shows that significantly improved predictions are obtained except for cases where the stratification pattern whether arising from the slow defluidization or the difference in the densities of the two species affects the mixing of the constituent species.

Experimental investigation of ash deposits on convection heating surfaces of a circulating fluidized bed municipal solid waste incinerator

Incineration of municipal solid waste （MSW） is a waste treatment method which can be sustainable in terms of waste volume reduction, as well as a source of renewable energy. During MSW combustion, increased formation of deposits on convection heating exchanger surfaces can pose severe operational problems, such as fouling, slagging and corrosion. These problems can cause lower heat transfer efficiency from the hot flue gas to the working fluid inside the tubes. A study was performed where experiments were carried out to examine the ash deposition characteristics in a full-scale MSW circulating fluidized bed （CFB） incinerator, using a newly designed deposit probe that was fitted with six thermocouples and four removable half rings. The influence of probe exposure time and probe surface temperature （500, 560, and 700℃） on ash deposit formation rate was investigated. The results indicate that the deposition mass and collection efficiency achieve a minimum at the probe surface temperature of 560~C. Ash particles are deposited on both the windward and leeward sides of the probe by impacting and thermophoretic/ condensation behavior. The major inorganic elements present in the ash deposits are Ca, A1 and Si. Compared to ash deposits formed on the leeward side of the probe, windward-side ash deposits contain relatively higher Ca and S concentrations, but lower levels of A1 and Si. Among all cases at different surface temperatures, the differences in elemental composition of the ash deposits from the leeward side are insignificant. However, as the surface temperature increases, the concentrations of A1, Si, K and Na in the windward-side ash deposits increase, but the Ca concentration is reduced. Finally, governing mechanisms are proposed on the basis of the experimental data, such as deposit morphology, elemental composition and thermodynamic calculations.

Numerical simulation of flow hydrodynamics of struvite pellets in a liquid–solid fluidized bed

Phosphorus recovery in the form of struvite has been aroused in recent decades for its dual advantages in eutrophication control and resource protection.The usage of the struvite products is normally determined by the size which is largely depended on the hydrodynamics.In this study,flow behavior of struvite pellets was simulated by means of Eulerian–Eulerian two-fluid model combining with kinetic theory of granular flow in a liquid–solid fluidized bed reactor（FBR）.A parametric study including the mesh size,time step,discretization strategy,turbulent model and drag model was first developed,followed by the evaluations of crucial operational conditions,particle characteristics and reactor shapes.The results showed that a cold model with the mesh resolution of 16 × 240,default time step of 0.001 sec and first order discretization scheme was accurate enough to describe the fluidization.The struvite holdup profile using Syamlal–O＇Brien drag model was best fitted to the experimental data as compared with other drag models and the empirical Richardson–Zaki equation.Regarding the model evaluation,it showed that liquid velocity and particle size played important roles on both solid holdups and velocities.The reactor diameter only influenced the solid velocity while the static bed height almost took no effect.These results are direct and can be applied to guide the operation and process control of the struvite fluidization.Moreover,the model parameters can also be used as the basic settings in further crystallization simulations.

Experimental investigation on flow asymmetry in solid entrance region of a square circulating fluidized bed

To study the influence of back feeding particles on gas-solid flow in the riser, this paper investigated the flow asymmetry in the solid entrance region of a fluidized bed by particle concentration/velocity measurements in a cold square circulating fluidized beds （CFB）. The pressure drop distribution along the riser and the saturation carrying capacity of gas for Geldart-B type particles were first analyzed. Under the condition of u0 = 4 m/s and Gs = 21 kg/（m^2 s）, the back feeding particles were found to penetrate the lean gas-solid flow near the entrance （rear） wall before reaching the opposite （front） wall, thus leading to a relatively denser region near the front wall in the bottom bed. Higher solid circulation rate （u0 =4 m/s, Gs = 33 kg/（m^2 s）） resulted in a higher particle concentration in the riser. However the back feeding particles with higher momentum increased the asymmetry of the particle concentration/velocity profile in the solid entrance region. Lower air velocity （u0 =3.2 m/s） and Gs =21 kg/（m2 s）, beyond the saturation carrying capacity of gas, induced an S-shaped axial solid distribution with a denser bottom zone. This limited the penetration of the back feeding particles and forced the flnidizing air to flow in the central region, thus leading to a higher solid holdup near the rear wall. Under the conditions of uo = 4 m/s and Gs = 21 kg/（m^2 s）, addition of coarse particles （dp= 1145 μm） into the bed made the radial distribution of solids more symmetrical.

Fluidization stability and periodic fluctuations in gas–solid separation fluidized bed using Geldart A dense medium

Gas–solid separation fluidized bed is a typical method for coal separation without water utilization.Geldart A particles is also considered as the ideal dense medium to strengthen separation efficiency.Fluidization stability reflects the bed pressure fluctuations and the distribution of bubble and emulsion phases,affecting the separation performance.And the main frequency of pressure fluctuations can directly reflect the degree of pressure fluctuations.Therefore,the detailed fluidization stability is analyzed combined the method of standard deviation of pressure fluctuations,power spectral density,etc.,for Geldart A particles.The results showed that maintaining an appropriate gas velocity resulted in an average bed pressure of around 2000 Pa.The main frequency is mainly concentrated around 1–1.5 Hz.Finally,a prediction model of the main frequency of pressure fluctuations is established,and the error can be controlled within±0.15.The investigation further proved the stable fluidization of Geldart A particles and provides a method for predicting the main frequency of pressure fluctuations in the gas–solid separation fluidized bed.

Role of inter-particle collision on solids acceleration in riser

Collision among particles plays a significant role in governing the structure of gas-solids flow in a riser, especially in the dense and acceleration region. The inter-particle collision is the major cause not only for the kinetic energy dissipation （in terms of additional pressure drop beyond the solids hold-up） but also for the control of solids acceleration （in terms of a balancing force to prevent a free acceleration of solids）. A neglect of the balancing force of inter-particle collision against the hydrodynamic force in the solids momentum equation would simply overestimate the solids acceleration or concentration while underestimate the axial gradient of pressure along the riser by a large margin, typically by up to two orders of magnitude. This paper aims to illustrate the importance of the collision on the characteristic of the gas-solids riser flow. Our analysis shows that the collision force should be of the same order of magnitude as that of the drag force in the dense and acceleration region, which can be far beyond that of gravitational force on solids. A simple formulation of the collision force is therefore proposed to bear a similar format of drag force, with regard to the dependence upon local solids properties.With the inclusion of the proposed correlation of collision force in the solids momentum equation, our model would be able to yield reasonable phase distributions of gas-solid flows, which can be validated, in a bulk range, against available measurements of solids volume fraction and axial gradient of pressure.

Performance evaluation of circulating fluidized bed incineration of municipal solid waste by multivariate outlier detection in China

This first nationwide survey was conducted to evaluate the overall performance of the circulating fluidized bed （CFB） incineration of municipal solid waste （MSW） in 2014-2015 in China. Total 23 CFB incineration power plants were evaluated. The data for monthly average flue gas emission of particles, CO, NOx, SO2 and HCl were collected over 12 consecutive months. The data were analyzed to assess the overall performance of CFB incineration by applying the Mahalanobis distance as a multivariate outlier detection method. Although the flue gas emission parameters had met the Chinese national emission standards, there were 11 total outliers （abnormal behavior） detected in 6 out of 23 CFB incineration oower olants from the oersoective of the MSW incineration performance. The results demonstrate that it is more important for a better perlbrmance of CFBs to reduce the lrequencles ot the MSW load changes, rather than the magnitudes of the MSW load changes, particularly reducing the frequencies in the range of 10% and more of the load changes, under the same and stable conditions. Furthermore, the overloading occurs more often than the underloading during the operation of the CFB incineration power plants in China. The frequent overloading is 0% to 30% ot the designed capacity. To achieve the stable performance of CFBs in practice, an appropriately designed MSW storage capacity is suggested to build in a plant to buffer and reduce the frequency of the load changes.

Solids behavior in dilute zone of a CFB riser under turbulent conditions

The behavior of the solid phase in the upper zone of a circulating fluidized bed riser was studied using a phase Doppler anemometer. Glass particles of mean diameter 107μm and superficial gas velocities UE covering the turbulent and the beginning of the fast fluidization regime were investigated. Three static bed heights were tested. Ascending and descending particles were found co-existing under all oper ating conditions tested, and at all measurement locations. Superficial gas velocity proved/happened to have a larger effect on descending particles at the wall and on ascending particles in the central region. Transversal particle velocities in both directions （toward the center and toward the wall） behaved rela- tively equivalently, with only slight difference observed at the wall. However, observation of the number of particles moving in either transversal direction showed a change in bed structure when increasing Ug. Furthermore, a balance was constantly observed between the core zone and the annulus zone where the mutual mass transfer between these two zones occurred continuously. Transition from a slow to a fast particle motion was accompanied by a transition to high levels of velocity fluctuations, and was found corresponding to the appearance of significant solid particle flow rate.

An experimental investigation into modeling solids friction for fluidized dense-phase pneumatic transport of powders

Results are presented of an ongoing investigation into modeling friction in fiuidized dense-phase pneumatic transport of bulk solids. Many popular modeling methods of the solids friction use the dimen- sionless solids loading ratio and Froude number. When evaluated under proper scale-up conditions of pipe diameter and length, many of these models have resulted in significant inaccuracy. A technique for modeling solids friction has been developed using a new combination of dimensionless numbers, volu- metric loading ratio and the ratio of particle free settling velocity to superficial conveying air velocity, to replace the solids loading ratio and Froude number. The models developed using the new formalism were evaluated for accuracy and stability under significant scale-up conditions for four different prod- ucts conveyed through four different test rigs （subject to diameter and length scale-up conditions）. The new model considerably improves predictions compared with those obtained using the existing model, especially in the dense-phase region. Whereas the latter yields absolute average relative errors varying between 10% and 86%, the former yielded results with errors from 4% to 20% for a wide range of scale-up conditions. This represents a more reliable and narrower range of prediction that is suitable for industrial scale-up requirements.

Application of CFD-DEM to the study of solid exchange in a dual-leg fluidized bed

The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.

Syngas production via coal char-CO_2 fluidized bed gasification and the effect on the performance of LSCFN//LSGM//LSCFN solid oxide fuel cell

Fluidized bed reactor is widely used in coal char-CO2 gasification. In this work, the production of syngas by using a fluidized bed gasification technique was first investigated and then the effect of the produced syngas on the performance of the solid oxide fuel cell with a configuration of La0.4Sr0.6Co0.2 Fe0.7 Nb0.1O3-δ//La0.8Sr0.2Ga0.83Mg0.17O3-δ//La0.4Sr0.6Co0.2Fe0.7Nb0.1O3-δ（LSCFN//LSGM//LSCFN） was studied. During the syngas production, we found that the volume fraction of CO increased with the increment of gasification temperature, and it reached a maximum value of 88.8%, corresponding to a composition of 0.76% H2, 88.8% CO, and 10.44% CO2, when the ratio of oxygen mass flow rate to that of coal char （Mo2/Mchar） increased to 0.29. In the following utilization of the produced syngas in solid oxide fuel cells, it was found that the increasing CO volume fraction in the syngas results in a gradual increase of the peak power density of the LSCFN//LSGM//LSCFN cell. The maximum peak power density of 410 mW/cm^2 was achieved for the syngas produced at 0.29 of Mo2/Mchar. In the stability test, the cell voltage decreased by 4% at a constant current density of 0.475 A/cm^2 after 54 h when fueled with the syngas with the composition of 0.76% H2, 88.8% CO, and 10.44% CO2. It reveals that a carbon deposition with the content of 13.66% in the anode is attributed to the cell performance degradation.

Local heat transfer properties in co- and counter-current G-L-S magnetically stabilized fluidized beds

Heat transfer coefficients were measured by immersed probes in co- and counter-current G-L-S magnetically stabilized fluidized beds （MSFBs） using air, water and nickel-alloy particles as the gas, liquid and solid phases. Influences of major factors, including magnetic field intensity, superficial gas and liquid velocities, liquid viscosity and surface tension, on heat-transfer properties were studied experimentally, indicating that both co- and counter-current G-L-S MSFB can provide relatively uniform radial distribution of heat transfer coefficients under appropriate operation conditions, thus controlling operation temperature for highly exothermic multi-phase reaction systems. Two correlations were provided to estimate accurately heat transfer properties in both co- and counter-current G-L-S MSFB systems, with an average error of less than 10%.

Hydrodynamic behavior of an internally circulating fluidized bed with tubular gas distributors

To better understand the hydrodynamic behavior of an internally circulating fluidized bed, solids holdup in the down-comer （Eso）, solids circulation rate （Gs） and gas bypassing fraction （from down-comer to riser y^R, and from riser to down-comer yRD） were experimentally studied. The effects of gas velocities in the riser and in the down-comer （UR and UD）, orifice diameter in the draft tube （dor）, and draft tube height （HR） were investigated. Experimental results showed that increase of gas velocities led to increase in Gs and yDR, and slight decrease in yeD. Larger orifice diameter on the draft tube led to higher 8sD, Gs and yDR, but had insignificant influence on YRD. with increasing draft tube height, both Gs and YDR first increased and then decreased, while yRD first decreased and then increased. Proposed correlations for predicting the hydrodynamic parameters agreed reasonably well with experimental values.

Experimental investigation into transient pressure pulses during pneumatic conveying of fine powders using Shannon entropy

This paper presents the results of an ongoing investigation into transient pressure pulses using Shan- non entropy. Pressure fluctuations （produced by gas-solid two-phase flow during fluidized dense-phase conveying） are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals （Mittal, Mallick, ＆ Wypych, 2014）. Two different powders, namely fly ash （median particle diameter 45 μm, particle den- sity 1950 kg/m3. loosely poured bulk density 950 kg/m3） and cement （median particle diameter 15 p,m, particle density 3060 kg/m3, loosely poured bulk density 1070 kg/m3）, are conveyed through different pipelines （51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length）. The transient nature of pressure fluc- tuations （instead of steady-state behavior） is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities （i.e. 3-5 m/s） and very high velocities （i.e. 11-14 m/s） while low Shannon entropy corresponds to mid-range velocities （i.e. 6-8 m/s）.

Simulation of mass balance behavior in a large-scale circulating fluidized bed reactor

We determine using a compound model the influence of the mass of granular matter on the behavior of a supercritical circulating fluidized bed （CFB） reactor. Population balance enables a stationary-regime modeling of the mass flow of granular matter inside a CFB unit in a large-scale. The simulation includes some important dynamic processes of gas-particle flows in fluidized bed such as attrition, fragmenta- tion, elutriation, and fuel combustion. Numerical calculations with full boiler loading were performed of operational parameters such as furnace temperature, furnace pressure, feeding materials mass flows, and excess air ratio. Furthermore, three bed inventory masses were adopted as experimental variables in the simulation model of mass balance. This approach enables a sensitivity study of mass flows of granular matter inside a CFB facility. Some computational results from this population balance model obtained for a supercritical CFB reactor are presented that show consistency with the operational data for large-scale CFB units.