Simulation of gas-solid flow characteristics of the circulating fluidized bed boiler under pure-oxygen combustion conditions

Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the conventional OFC technology usually depends on the flue gas recirculation system,which faces significant investment,high energy consumption,and potential low-temperature corrosion problem.Considering these deficiencies,the direct utilization of pure oxygen to achieve particle fluidization and fuel combustion may reduce the overall energy consumption and CO_(2)-capture costs.In this paper,the fundamental structure of a self-designed 130 t·h^(-1) pure-oxygen combustion circulating fluidized bed(CFB)boiler was provided,and the computational particle fluid dynamics method was used to analyze the gas-solid flow characteristics of this new-concept boiler under different working conditions.The results indicate that through the careful selection of design or operational parameters,such as average bed-material size and fluidization velocity,the pure-oxygen combustion CFB system can maintain the ideal fluidization state,namely significant internal and external particle circulation.Besides,the contraction section of the boiler leads to the particle backflow in the lower furnace,resulting in the particle suspension concentration near the wall region being higher than that in the center region.Conversely,the upper furnace still retains the classic core-annulus flow structure.In addition to increasing solid circulation rate by reducing the average bed-material size,altering primary gas ratio and bed inventory can also exert varying degrees of influence on the gas-solid flow characteristics of the pure-oxygen combustion CFB boiler.

THE NONLINEAR STABILITY OF PLANE PARALLEL SHEAR FLOWS WITH RESPECT TO TILTED PERTURBATIONS

The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.

Identifying the enhancement mechanism of Al/MoO_(3) reactive multilayered films on the ignition ability of semiconductor bridge using a one-dimensional gas-solid two-phase flow model

Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.

Effect of shear-induced contact area and aperture variations on nonlinear flow behaviors in fractal rock fractures

This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid flow and visualization tests is performed on four transparent fracture specimens with various shear displacements of 1 mm,3 mm,5 mm,7 mm and 10 mm under a normal stress of 0.5 MPa.Four granite fractures with different roughnesses are selected and quantified using variogram fractal dimensions.The obtained results show that the critical Reynolds number tends to increase with increasing shear displacement but decrease with increasing roughness of fracture surface.The flow paths are more tortuous at the beginning of shear because of the wide distribution of small contact spots.As the shear displacement continues to increase,preferential flow paths are more distinctly observed due to the decrease in the number of contact spots caused by shear dilation;yet the area of single contacts in-creases.Based on the experimental results,an empirical mathematical equation is proposed to quantify the critical Reynolds number using the contact area ratio and fractal dimension.

Global Well-Posedness for Aggregation Equation with Time-Space Nonlocal Operator and Shear Flow

In this paper,we consider the two-dimensional aggregation equation with the shear flow and time-space nonlocal attractive operator.Without the advection,the solution of the aggregation equation may blow up in finite time.We show that the shear flow can suppress the blow-up.

V Flow技术测量颈动脉壁面剪应力的一致性研究

目的评价血流向量成像(V Flow)技术在测量健康成年人颈动脉壁面剪应力(WSS)中的一致性。方法于2021年2月至2021年3月招募健康成年志愿者20人,由2名不同年资的超声医师使用配备V Flow功能的Mindray Resona 7超声仪和3~9 MHz线阵探头进行双侧颈动脉扫查,分别采集双侧颈总动脉远段、颈总动脉分叉至颈内动脉起始部两段动脉的动态V Flow图像,测量两侧颈总动脉远段的近心端、远心端、分叉处及颈内动脉起始部的前、后壁的WSS,使用组内相关系数(ICC)和Bland-Altman图评估组内一致性及组间一致性。结果双侧颈动脉前、后壁的4个不同节段之间WSS值差异均有统计学意义(P<0.05)。高年资超声医师2次测量结果的一致性结果显示,左侧颈总动脉远段的远心端一致性极好(ICC 0.779),右侧颈总动脉远段的近心端(ICC 0.605)、远心端(ICC 0.585)、颈内动脉起始部(ICC 0.457)、左侧颈总动脉分叉处(ICC 0.606)及颈内动脉起始部(ICC 0.702)一致性均较好;不同年资超声医师的测量结果显示,仅右侧颈总动脉分叉处(ICC 0.486)及左侧颈总动脉远段的远心端(ICC 0.576)一致性较好。结论V Flow技术可显示不同位点间颈动脉WSS之间的差异,其组内一致性较好,但在不同年资超声医师间存在一定的差异。

Investigation of high rate mechanical flow followed by ignition for high-energy propellant under dynamic extrusion loading

Investigating the ignition response of nitrate ester plasticized polyether(NEPE) propellant under dynamic extrusion loading is of great significant at least for two cases. Firstly, it helps to understand the mechanism and conditions of unwanted ignition inside charged propellant under accident stimulus.Secondly, evaluates the risk of a shell crevice in a solid rocket motor(SRM) under a falling or overturning scene. In the present study, an innovative visual crevice extrusion experiment is designed using a dropweight apparatus. The dynamic responses of NEPE propellant during extrusion loading, including compaction and compression, rapid shear flow into the crevice, stress concentration, and ignition reaction, have been firstly observed using a high-performance high-speed camera. The ignition reaction is observed in the triangular region of the NEPE propellant sample above the crevice when the drop weight velocity was 1.90 m/s. Based on the user material subroutine interface UMAT provided by finite element software LS-DYNA, a viscoelastic-plastic model and dual ignition criterion related to plastic shear dissipation are developed and applied to the local ignition response analysis under crevice extrusion conditions. The stress concentration occurs in the crevice location of the propellant sample, the shear stress is relatively large, the effective plastic work is relatively large, and the ignition reaction is easy to occur. When the sample thickness decreases from 5 mm to 2.5 mm, the shear stress increases from 22.3 MPa to 28.6 MPa, the critical value of effective plastic work required for ignition is shortened from 1280 μs to 730 μs, and the triangular area is easily triggering an ignition reaction. The propellant sample with a small thickness is more likely to stress concentration, resulting in large shear stress and effective work, triggering an ignition reaction.

Boulder-induced form roughness and skin shear stresses in a gravel-bed stream

Boulder spacing in mountain rivers and near-wake flow zones within the boulder array is very useful for fish habitat and growth of aquatic organisms.The present study aims to investigate how the boulder array and spacing influence the near-bed flow structures in a gravel-bed stream.Boulders are staggered over a gravel-bed stream with three different inter-boulder spacing namely(a)large(b)medium and(c)small spacing.An acoustic Doppler velocimeter was used for flow measurements in a rectangular channel and the results were compared with those acquired from numerical simulation.The time-averaged velocity profiles at the near-wake flow zones of boulders experience maximum flow retardation which is an outcome of the boulder-induced form roughness.The ratio of velocity differences associated to form and skin roughness and its positive magnitude reveals the dominance of form roughness closest to the boulders.Form roughness computed is 1.75 to 2 times higher than the skin roughness at the near-wake flow region.In particular,the collective immobile boulders placed at different inter-boulder spacings developed high and low bed shear stresses closest to the boulders.The low bed shear stresses characterised by a secondary peak developed at the trough location of the boulders is attributed to the skin shear stress.Further,the spatial averaging of time-averaged flow quantities gives additional impetus to present an improved illustration of fluid shear stresses.The formation of form-induced shear stress is estimated to be 17%to 23%of doubleaveraged Reynolds shear stress and partially compensates for the damping of time-averaged Reynolds shear stress in the interfacial sub-layer.The quadrant analysis of spatial velocity fluctuations depicts that the form-induced shear stresses are dominant in the interfacial sub-layer and have no significance above the gravel-bed surface.

Seismic anisotropy and upper mantle dynamics in Alaska:A review of shear wave splitting analyses

Shear wave splitting(SWS)is regarded as the most effective geophysical method to delineate mantle flow fields by detecting seismic azimuthal anisotropy in the earth's upper mantle,especially in tectonically active regions such as subduction zones.The Aleutian-Alaska subduction zone has a convergence rate of approximately 50 mm/yr,with a trench length reaching nearly 2800 km.Such a long subduction zone has led to intensive continental deformation and numerous strong earthquakes in southern and central Alaska,while northern Alaska is relatively inactive.The sharp contrast makes Alaska a favorable locale to investigate the impact of subduction on mantle dynamics.Moreover,the uniqueness of this subduction zone,including the unusual subducting type,varying slab geometry,and atypical magmatic activity and composition,has intrigued the curiosity of many geoscientists.To identify different sources of seismic anisotropy beneath the Alaska region and probe the influence of a geometrically varying subducting slab on mantle dynamics,extensive SWS analyses have been conducted in the past decades.However,the insufficient station and azimuthal coverage,especially in early studies,not only led to some conflicting results but also strongly limited the in-depth investigation of layered anisotropy and the estimation of anisotropy depth.With the completion of the Transportable Array project in Alaska,recent studies have revealed more detailed mantle structures and characteristics based on the dense station coverage and newly collected massive seismic data.In this study,we review significant regional-and continental-scale SWS studies in the Alaska region and conclude the mantle flow fields therein,to understand how a geometrically varying subducting slab alters the regional mantle dynamics.The summarized mantle flow mechanisms are believed to be conducive to the understanding of seismic anisotropy patterns in other subduction zones with a complicated tectonic setting.

Theory of Flexural Shear, Bending and Torsion for a Thin-Walled Beam of Open Section

Aspects of the general Vlasov theory are examined separately as applied to a thin-walled channel section cantilever beam under free-end end loading. In particular, the flexural bending and shear that arise under transverse shear and axial torsional loading are each considered theoretically. These analyses involve the location of the shear centre at which transverse shear forces when applied do not produce torsion. This centre, when taken to be coincident with the centre of twist implies an equivalent reciprocal behaviour. That is, an axial torsion applied concentric with the shear centre will twist but not bend the beam. The respective bending and shear stress conversions are derived for each action applied to three aluminium alloy extruded channel sections mounted as cantilevers with a horizontal principal axis of symmetry. Bending and shear are considered more generally for other thin-walled sections when the transverse loading axes at the shear centre are not parallel to the section = s centroidal axes of principal second moments of area. The fixing at one end of the cantilever modifies the St Venant free angular twist and the free warping displacement. It is shown from the Wagner-Kappus torsion theory how the end constrained warping generates an axial stress distribution that varies with the length and across the cross-section for an axial torsion applied to the shear centre. It should be mentioned here for wider applications and validation of the Vlasov theory that attendant papers are to consider in detail bending and torsional loadings applied to other axes through each of the centroid and the web centre. Therein, both bending and twisting arise from transverse shear and axial torsion applied to each position being displaced from the shear centre. Here, the influence of the axis position upon the net axial and shear stress distributions is to be established. That is, the net axial stress from axial torsional loading is identified with the sum of axial stress due to bending and axial stress arising from constrained warping displacements at the fixing. The net shear stress distribution overlays the distributions from axial torsion and that from flexural shear under transverse loading. Both arise when transverse forces are displaced from the shear centre.

Review of gas-solid two phase flow rate-concentration detection technology

The indirect detection method basic principle of rate and concentration,application range and research results on gassolid two phase flow were discussed.The present development situation and the existing problems of rate and concentration detection technology were analyzed and summarized.Emphatically analyzed the existing problems in the industrial application and research status of electrostatic method in measuring phase concentration.Design criterion of electrostatic phase concentration sensor is given,the superiority and wide industrial application prospect of the sensor used for phase concentration measurement are clarified.

Vortex-Induced Vibrations of A Long Flexible Cylinder in Linear and Exponential Shear Flows

A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.

An application of interacting shear flows theory: exact solution for unsteady oblique stagnation point flow

An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.

Particle dynamics in dense shear granular flow

The particle dynamics in an annular shear granular flow is studied using the discrete element method, and the influences of packing fraction, shear rate and friction coefficient are analyzed. We demonstrate the existence of a critical packing fraction exists in the shear granular flow. When the packing fraction is lower than this critical value, the mean tangential velocity profile exhibits a rate-independent feature. However, when the packing fraction exceeds this critical value, the tangential velocity profile becomes rate-dependent and varies gradually from linear to nonlinear with increasing shear rate. Furthermore, we find a continuous transition from the unjammed state to the jammed state in a shear granular flow as the packing fraction increases. In this transforming process, the force distribution varies distinctly and the contact force network also exhibits different features.

Numerical Simulation of the Vortex-Induced Vibration of A Curved Flexible Riser in Shear Flow

A series of fully three-dimensional(3 D) numerical simulations of flow past a free-to-oscillate curved flexible riser in shear flow were conducted at Reynolds number of 185–1015. The numerical results obtained by the two-way fluid–structure interaction(FSI) simulations are in good agreement with the experimental results reported in the earlier study. It is further found that the frequency transition is out of phase not only in the inline(IL) and crossflow(CF) directions but also along the span direction. The mode competition leads to the non-zero nodes of the rootmean-square(RMS) amplitude and the relatively chaotic trajectories. The fluid–structure interaction is to some extent reflected by the transverse velocity of the ambient fluid, which reaches the maximum value when the riser reaches the equilibrium position. Moreover, the local maximum transverse velocities occur at the peak CF amplitudes, and the values are relatively large when the vibration is in the resonance regions. The 3 D vortex columns are shed nearly parallel to the axis of the curved flexible riser. As the local Reynolds number increases from 0 at the bottom of the riser to the maximum value at the top, the wake undergoes a transition from a two-dimensional structure to a 3 D one. More irregular small-scale vortices appeared at the wake region of the riser, undergoing large amplitude responses.

Shear flow induced vibrations of long slender cylinders with a wake oscillator model

A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations（VIV） when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.

Discrete Particle Simulation of Gas-Solid Flow in Air-Blowing Seed Metering Device

In this paper,the gas and seed flow characters in the air-blowing seed metering device are investigated by using the coupled computational fluid dynamics and discrete element method(CFD-DEM)in three dimensions(3D).The method of establishing boundary model based on the computer-aided design(CAD)drawing,has been used to build the boundary model of seed metering device.The 3D laser scanning technique and multi-element method are adopted to establish the particle model.Through a combined numerical and experimental effort,using 3D CFD-DEM software,which is based on the in-house codes,the mechanisms governing the gas and solid dynamic behaviors in the seed metering device have been studied.The gas velocity field and the effect of different rotational speeds and air pressures on the seeding performance and particle velocity have been studied,similar agreements between numerical and experimental results are gained.This reveals that the 3D CFD-DEM model established is able to predict the performance of the air-blowing seed metering device.It can be used to design and optimize the air-blowing seed metering device and other similar agriculture devices.

Flow Behaviors of Gas-Solid Injector by 3D Simulation with Kinetic Theory of Granular Flow

A computational study on the flow behavior of a gas-solid injector by Eulerian approach was carried out. The gas phase was modeled with k-ε turbulent model and the particle phase was modeled with kinetic theory of granular flow. The simulations by Eulerian two-fluid model （TFM） were compared with the corresponding results by discrete element method （DEM） and experiments. It was showed that TFM simulated results were in reasonable agreement with the experimental and DEM simulated results. Based on TFM simulations, gas-solid flow pattern, gas velocity, particle velocity and the static pressure under different driving jet velocity, backpressure and convergent section angle were obtained. The results showed that the time average axial gas velocity sharply decreased and then slightly increased to a constant value in the horizontal conveying pipe. The time average axial particle velocity increased initially and then decreased, but in the outlet region of the convergent section the particle velocity remarkably increased once more to the maximal value. As a whole, the static pressure distribution change trends were found to be independent on driving gas velocity, backpressure and convergent section angle. However, the static pressure increased with increase of convergent section angle and gas jet velocities. The difference of static pressure to backpressure increased with increasing backpressure.

Modeling of gas-solid flow in a CFB riser based on computational particle fluid dynamics

A three-dimensional model for gas-solid flow in a circulating fluidized bed（CFB） riser was developed based on computational particle fluid dynamics（CPFD）.The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport（DPT） regime and a fast fluidization（FF） regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.

Numerical Simulation on Gas-Solid Two-Phase Turbulent Flow in FCC Riser Reactors(Ⅰ) Turbulent Gas-Solid Flow-Reaction Model

Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.