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.     

Torsional static and vibration analysis of functionally graded nanotube with bi-Helmholtz kernel based stress-driven nonlocal integral model

A torsional static and free vibration analysis of the functionally graded nanotube(FGNT)composed of two materials varying continuously according to the power-law along the radial direction is performed using the bi-Helmholtz kernel based stress-driven nonlocal integral model.The differential governing equation and boundary conditions are deduced on the basis of Hamilton’s principle,and the constitutive relationship is expressed as an integral equation with the bi-Helmholtz kernel.Several nominal variables are introduced to simplify the differential governing equation,integral constitutive equation,and boundary conditions.Rather than transforming the constitutive equation from integral to differential forms,the Laplace transformation is used directly to solve the integro-differential equations.The explicit expression for nominal torsional rotation and torque contains four unknown constants,which can be determined with the help of two boundary conditions and two extra constraints from the integral constitutive relation.A few benchmarked examples are solved to illustrate the nonlocal influence on the static torsion of a clamped-clamped(CC)FGNT under torsional constraints and a clamped-free(CF)FGNT under concentrated and uniformly distributed torques as well as the torsional free vibration of an FGNT under different boundary conditions.     

Analysis on nonlinear effect of unsteady percolation in the inhomogeneous shale gas reservoir

The nonlinear effects of unsteady multi-scale shale gas percolation,such as desorption,slippage,diffusion,pressure-dependent viscosity,and compressibility,are investigated by numerical simulation.A new general mathematical model of the problem is built,in which the Gaussian distribution is used to describe the inhomogeneous intrinsic permeability.Based on the Boltzmann transformation,an efficient semi-analytical method is proposed.The problem is then converted into a nonlinear equation in an integral form for the pressure field,and a related explicit iteration scheme is constructed by numerical discretization.The validation examples show that the proposed method has good convergence,and the simulation results also agree well with the results obtained from both numerical and actual data of two vertical fractured test wells in the literature.Desorption,slippage,and diffusion have significant influence on shale gas flows.The accuracy of the usual technique that the product of viscosity and compressibility is approximated as its value at the average formation pressure is examined.     

Active control of multiscale features in wall-bounded turbulence

This study experimentally investigates the impact of a single piezoelectric(PZT)actuator on a turbulent boundary layer from a statistical viewpoint.The working conditions of the actuator include a range of frequencies and amplitudes.The streamwise velocity signals in the turbulent boundary layer flow are measured downstream of the actuator using a hot-wire anemometer.The mean velocity profiles and other basic parameters are reported.Spectra results obtained by discrete wavelet decomposition indicate that the PZT vibration primarily influences the near-wall region.The turbulent intensities at different scales suggest that the actuator redistributes the near-wall turbulent energy.The skewness and flatness distributions show that the actuator effectively alters the sweep events and reduces intermittency at smaller scales.Moreover,under the impact of the PZT actuator,the symmetry of vibration scales’velocity signals is promoted and the structural composition appears in an orderly manner.Probability distribution function results indicate that perturbation causes the fluctuations in vibration scales and smaller scales with high intensity and low intermittency.Based on the flatness factor,the bursting process is also detected.The vibrations reduce the relative intensities of the burst events,indicating that the streamwise vortices in the buffer layer experience direct interference due to the PZT control.     

Asymptotics for filtration of polydisperse suspension with small impurities

A model for deep bed filtration of a polydisperse suspension with small impurities in a porous medium is considered.Different suspended particles move with the same velocity as the carrier water and get blocked in the pore throats due to the size-exclusion mechanism of particle retention.A solution of the model in the form of a traveling wave is obtained.The global exact solution for a multiparticle filtration with one high concentration and several low concentrations of suspended particles is obtained in an explicit form.The analytic solutions for a bidisperse suspension with large and small particles are constructed.The profiles of the retained small particles change monotony with time.The global asymptotics for the filtration of a polydisperse suspension with small kinetic rates is constructed in the whole filtration zone.     

一种Cartesian气动网格的自适应划分算法

提出了针对空间Cartesian结构化网格的自适应划分算法,以满足飞机总体设计过程中进行多学科设计优化时根据方案变化生成CFD计算网格的需要,使气动特性能够得以自动化分析,并且保证优化迭代过程的顺利进行.分析了Cartesian网格划分过程中的两个关键问题,并给出了解决方法.在此基础之上,从理论上推导了能够自动适应飞机布局形式以及几何尺寸任意变化的Cartesian网格生成算法和实际划分步骤,并给出了其中部分主要的计算公式,同时论述了自适应划分过程中的若干重要问题.阐述了自适应划分所需的飞机及其部件特征尺寸的获取途径.最后展示了几个由自适应划分算法生成的Cartesian网格实例以说明其实际应用.     

低温超导核磁共振陀螺仪模型

建立了单工作物质的三自由度3He低温超导核磁共振陀螺仪结构;利用了量子力学和经典动力学,经过严密的力学分析和数学演算,给出了三轴陀螺仪的工作原理和结构示意图;针对陀螺仪的交叉轴角速率耦合问题,给出了附加磁场线圈解耦法和冗余设计结构,最后建立了基于超导量子干涉仪探测磁矩,并采用最小二乘估计法来推导陀螺进动频率.陀螺仪测速范围可以达到10-9～103 rad/s,漂移为10-4 (°)/h.该陀螺仪结合低温超导技术具有高精度的前景.需要进一步对超导量子干涉磁矩检测仪的精度与陀螺性能进行研究.     

DEM investigation of macro-and micro-mechanical properties of rigid-grain and soft-chip mixtures

We investigated the macro-and micro-mechanical properties of rigid-grain and soft-chip mixtures(GCMs)through numerical simulations using the discrete element method.We present a novel framework for the discrete modeling of soft chips and rigid grains in conjunction with calibration processes.Several numerical triaxial tests were also performed on GCMs with 0%,10%,20%,and 30%volumetric chip contents,P.The simulation results demonstrate that increasing P leads to higher GCM toughness,higher deviatoric peak stress,and higher corresponding shear strain.Higher P also contributes to more volume contraction and less dilation.The friction angles at both the peak and residual state significantly increase with increasing P.In view of the micro-mechanical features,strong contact force chains develop along the loading direction,which results in considerable anisotropy in the peak and residual states.Both the formation of strong force chains and rotation of grains decrease with increasing P,whereas the grain sliding percentage increases.The tensile force is mobilized with shearing and higher P leads to less mobilization of the tensile force.These findings are useful for better understanding the internal structure of GCMs with different soft-chip contents,especially in granular mixture mechanics and geomechanics.     

Interaction effects on heat conduction and thermal stress in an infinite elastic plane with two circular holes

Numerous researches have focused on the physical behavior of an elastic material in the vicinity of a single hole under the assumption that the interaction effects arising from the introduction of multiple holes remain negligible if the holes are placed sufficiently far from each other.In an effort to understand hole interaction effects on heat conduction and thermal stress,we consider the case when two circular holes are embedded in an infinite elastic material and use complex variable methods together with numerical analysis to obtain solutions describing temperature and elastic fields in the vicinity of the two circular holes.The results indicate that the interaction effects on temperature distribution and stress strongly depend on the relative size of the two holes and the distance placed between them but not on the actual size of the holes.     

?????????

编辑部的同志要我写一点回忆,以纪念<力学与实践>创刊30年.往事悠悠,使我又回到20世纪70年代那动荡不安,却又是激情四射的年代.     

Using the discrete element method to assess the mixing of polydisperse solid particles in a rotary drum

Despite the wide applications of powder and solid mixing in industry, knowledge on the mixing of polydisperse solid particles in rotary drum blenders is lacking. This study investigates the mixing of monodisperse, bidisperse, tridisperse, and polydisperse solid particles in a rotary drum using the discrete element method. To validate the model developed in this study, experimental and simulation results were compared. The validated model was then employed to investigate the effects of the drum rotational speed, particle size, and initial loading method on the mixing quality. The degree of mixing of polydisperse particles was smaller than that for monodisperse particles owing to the segregation phenomenon. The mixing index increased from an initial value to a maximum and decreased slightly before reaching a plateau for bidisperse, tridisperse, and polydisperse particles as a direct result of the segregation of particles of different sizes. Final mixing indices were higher for polydisperse particles than for tridisperse and bidisperse particles. Additionally, segregation was weakened by introducing additional particles of intermediate size. The best mixing of bidisperse and tridisperse particles was achieved for top–bottom smaller-to-larger initial loading, while that of polydisperse systems was achieved using top–bottom smaller-to-larger and top–bottom larger-to-smaller initial loading methods.     

Experimental and computational studies on flow behavior of gas-solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally.Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695--0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

A solids mixing rate correlation for small scale fluidized beds

A new first degree solids mixing rate is proposed to evaluate the mixing of solids in small scale fluidized beds. Particle mixing experiments were carried out in a 2D fluidized bed with a cross-section of 0.02 m × 0.2 m and a height of 1 m. White and black particles with average diameters of 850 and 450 μm were used in our experiments. Image processing was used to measure the concentration of the tracers at different times. The effects of four representative operating parameters (superficial gas velocity, ratio of tracer particles to bed particles, tracer particle position, and particle size) on mixing are discussed with reference to the mixing index. We found that the Lacey index depends on the concentration of the tracers. The position of the tracers affects the initial mixing rate but not the final degree of mixing. However, the new mixing rate equation does not depend on the initial configuration of the particles because this situation is considered to be the initial condition. Using the data obtained in this work and that found in literature, an empirical correlation is proposed to evaluate the mixing rate constant as a function of dimensionless numbers (Archimedes, Reynolds, and Froude) in small scale fluidized beds. This correlation allows for an estimation of the mixing rate under different operating conditions and for the detection of the end point and/or the time of mixing.     

Experimental and computational studies on flow behavior of gas–solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally. Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695-0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

Numerical simulation of tetrahedral particle mixing and motion in rotating drums

A regular tetrahedron is the simplest three-dimensional structure and has the largest non-sphericity. Mixing of tetrahedral particles in a thin drum mixer was studied by the soft-sphere-imbedded pseudo-hard particle model and compared with that of spherical particles. The two particle types were simulated with different rotation speeds and drum filling levels. The Lacey mixing index and Shannon information entropy were used to explore the effects of sphericity on the mixing and motion of particles. Moreover, the probability density functions and mean values and variances of motion velocities, including translational and rotational, were computed to quantify the differences between the motion features of tetrahedra and spheres. We found that the flow regime depended on the particle shape in addition to the rotation speed and filling level of the drum. The mixing of tetrahedral particles was better than that of spherical particles in the rolling and cascading regimes at a high filling level, whereas it may be poorer when the filling level was low. The Shannon information entropy is better than the Lacey mixing index to evaluate mixing because it can reflect the real change of flow regime from the cataracting to the centrifugal regime, whereas the mixing index cannot.     

Investigation of dense slurry suspensions with coaxial mixers: Influences of design variables through tomography and mathematical modelling

The coaxial mixers enhance the suspension of concentrated slurries in an agitated reactor. In this research work, the complex slurry suspension and dissemination behavior in a coaxial slurry mixing system (comprised of a close clearance anchor rotating with a low speed and an inner axial impeller rotating with a high speed) was analyzed employing ERT (electrical resistance tomography, a non-intrusive flow visualization technique), and computational fluid dynamics (CFD). The numerical models were validated by comparing the axial solid concentration profiles generated using the ERT data and the CFD simulation results. The influences of various important parameters such as the diameter of the inner axial impeller, the inner impeller type, and the inner impeller spacing on the hydrodynamic characteristics of the slurry suspensions in a coaxial mixing vessel were thoroughly analyzed. The radial and axial velocity profiles of solid particles were generated using the validated mathematical models. The assessment of energy loss due to the solid–solid collisions, the particle–fluid frictions, and the particle–vessel wall collisions was conducted. The evaluation of optimum inner impeller clearance and inner impeller diameter is essential to attain a high degree of solids suspension and dissemination in a coaxial slurry mixing system.     

Comparison of near-infrared spectroscopy and DEM simulations for tracing spatial dispersion of water during mixing

The aim of this work was to validate a wet mixing process, in which a liquid spray is used to impregnate particles during mixing. The experimental results obtained using a bladed-mixer with a near-infrared sensor were compared with the results obtained using a 1:1 discrete element method simulation. The porous particles used in both cases absorbed the sprayed liquid for a process time of about 18 min. Multiple sensors attached to the mixer wall continuously monitored the liquid contents of passing particles. The sensors were modeled in the simulation and the resulting signals were analyzed and compared with the experimental results. We show that the algorithms used for spray and liquid absorption can be used to predict the moisture distribution inside granular materials in chemical and pharmaceutical processes. Such simulations can help to save money, e.g., in resource-intensive experimental plans and equipment design studies, and by varying material parameters.     

Mixing and segregation assessment of bi-disperse solid particles in a double paddle mixer

A double paddle blender's flow patterns and mixing mechanisms were analyzed using discrete element method (DEM) and experiments. The mixing performance of this type of the blender containing bi-disperse particles has been rarely studied in the literature. Plackett-Burman design of experiments (DoE) methodology was used to calibrate the DEM input parameters. Subsequently, the impact of the particle number ratio, vessel fill level, and paddle rotational speed on mixing performance was investigated using the calibrated DEM model. The mixing performance was assessed using relative standard deviation and segregation intensity. Mixing performance was significantly affected by the paddle rotational speed and particle number ratio. Moreover, the Peclet number and diffusivity coefficient were used to evaluate the mixing mechanism in the blender. Results revealed that the diffusion was the predominant mixing mechanism, and the best mixing performance was observed when the diffusivity coefficients of 3 mm and 5 mm particles were almost equal.     

Effect of mixing structure on the hygroscopic behavior of ultrafine ammonium sulfate particles mixed with succinic acid and levoglucosan

Understanding the interactions between water and atmospheric aerosols is critical for estimating their impact on the radiation budget and cloud formation. The hygroscopic behavior of ultrafine (<100 nm) ammonium sulfate particles internally mixed with either succinic acid (slightly soluble) or levoglucosan (soluble) in different mixing structures (core-shell vs. well-mixed) were measured using a hygroscopicity tandem differential mobility analyzer (HTDMA). During the hydration process (6–92% relative humidity (RH)), the size of core-shell particles (ammonium sulfate and succinic acid) remained unchanged until a slow increase in particle size occurred at 79% RH; however, an abrupt increase in size (i.e., a clear deliquescence) was observed at ∼72% RH for well-mixed particles with a similar volume fraction to the core-shell particles (80:20 by volume). This increase might occur because the shell hindered the complete dissolution of the core-shell particles below 92% RH. The onset RH value was lower for the ammonium sulfate/levoglucosan core-shell particles than the ammonium sulfate/succinic acid core-shell particles due to levoglucosan's higher solubility relative to succinic acid. The growth factor (GF) of the core-shell particles was lower than that of the well-mixed particles, while the GF of the ammonium sulfate/levoglucosan particles was higher than that of ammonium sulfate/succinic acid particles with the same volume fractions. As the volume fraction of the organic species increased, the GF decreased. The data suggest that the mixing structure is also important when determining hygroscopic behavior of the mixed particles.