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.     

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

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

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

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

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年代那动荡不安,却又是激情四射的年代.     

Two non-probabilistic set-theoretical models for dynamic response and buckling failure measures of bars with unknown-but-bounded initial imperfections

This paper is concerned with the problem of comparison of two non-probabilistic set-theoretical models for dynamic response and buckling failure measures of bars with unknown-but-bounded initial imperfections. Two kinds of non-probabilistic set-theoretical models are convex models and interval analysis models. In convex models and interval analysis models, the uncertain quantities are considered to be unknown except that they belong to given sets in an appropriate vector space. In this case, all information about the dynamic response and buckling failure measures of bars is provided by the set of dynamic responses and buckling failure measures consistent with the constraints on the uncertain quantities. The dynamic response estimate is actually a set in appropriate response space rather than a single vector. The set estimate is the smallest calculable set which contains the uncertain dynamic response, but it is usually impractical to calculate this set. Two set estimate methods are developed which can calculate the time varying box or hyperrectangle, i.e. interval vector in the response space that contains the true dynamic response. Comparison between convex models and interval analysis models in mathematical proofs and numerical calculations shows that, under the condition of the outer enclosed ellipsoid from a hyperrectangle or an interval vector, the set dynamic response predicted by interval analysis models is smaller than that yielded by convex models; under the condition of the outer enclosed hyperrectangle or an interval vector from an ellipsoid, the dynamic response set calculated by convex models is smaller than that obtained by interval analysis models.     

Extension of convex models and its improvement on the approximate solution

In this paper, by means of combining non probabilistic convex modeling with perturbation theory, an improvement is made on the first order approximate solution in convex models of uncertainties. Convex modeling is extended to largely uncertain and non-convex sets of uncertainties and the combinational convex modeling is developed. The presented method not only extends applications of convex modeling, but also improves its accuracy in uncertain problems and computational efficiency. The numerical example illustrates the efficiency of the proposed method. The project supported by the National Outstanding Youth Science Foundation of China and the National Post Doctor Science Foundation of China     

基于Taylor展式的不确定结构复特征值问题两种非概率方法比较研究

代替传统的处理不确定问题的概率统计方法,将利用区问数学和凸模型理论研究具有有界不确定参数的非比例阻尼结构复特征值所在区域问题．区间数学将有界不确定结构参数用超长方体即区问向量进行定量化,而凸模型理论则用椭球对有界不确定参数进行定量化．在不用知道不确定变量的概率统计特性的条件下,区间分析方法和凸模型理论都可以确定出有界不确定结构参数的非比例阻尼结构复特征值所在区域．通过数学证明和数值算例来说明,在凸模型理论中的椭球在由区间分析中的超长方体—区间向量来确定的条件下,由区间数学所确定出不确定结构复特征值实部和虚部的宽度要比凸模型所确定出的范围的宽度要小,而这正是工程技术人员所要求的结果。     

Vertex solution theorem for the upper and lower bounds on the dynamic response of structures with uncertain-but-bounded parameters

The aim of this paper is to evaluate the effects of uncertain-but-bounded parameters on the dynamic response of structures. By combining the interval mathematics and the finite element analysis, the mass matrix, damping matrix, stiffness matrix and the external loads are represented as interval matrices and vector. With the help of the optimization theory, we present the vertex solution theorem for determining both the exact upper bounds or maximum values and the exact lower bounds or minimum values of the dynamic response of structures, in which these parameters reach their extreme values on the boundary of the interval mass, damping, stiffness matrices and the interval extemal loads vector. Three examples are used to illustrate the computational aspects of the presented vertex solution theorem.     

椭球凸模型非概率可靠性度量和区间安全系数的关系

研究了椭球凸模型非概率可靠性度量和区间模型安全系数的关系。根据基于椭球凸模型的非概率可靠性指标和非概率可靠度的定义,建立了两者之间的函数关系;按照区间模型安全系数的定义,给出了由椭球参数确定的3种区间模型安全系数,分析了它们的意义;建立了非概率可靠性指标和区间模型安全系数之间的解析关系,讨论了它们在评估结构可靠性或安全程度上的意义;通过数值算例验证了分析结果。     

THE UNCERTAIN EIGENVALUES ANALYSIS OF STRUCTURES WITH UNCERTAIN PARAMETERS

This paper presents a method for estimating the upper and lower bounds of eigenvalues of structures with uncertainties. The uncertain parameters are described by the convex model. A numerical example of the frame structure is given to illustrate the efficiency of the method. This work is supported by the National Natural Science Foundation of China(No. 19872028).     

Combination of structural reliability and interval analysis

In engineering applications, probabilistic reliability theory appears to be presently the most important method, however, in many cases precise probabilistic reliability theory cannot be considered as adequate and credible model of the real state of actual affairs. In this paper, we developed a hybrid of probabilistic and non-probabilistic reliability theory, which describes the structural uncertain parameters as interval variables when statistical data are found insufficient. By using the interval analysis, a new method for calculating the interval of the structural reliability as well as the reliability index is introduced in this paper, and the traditional probabilistic theory is incorporated with the interval analysis. Moreover, the new method preserves the useful part of the traditional probabilistic reliability theory, but removes the restriction of its strict requirement on data acquisition. Example is presented to demonstrate the feasibility and validity of the proposed theory.     

Stability and strict positive realness of convex polytopes of interval polynomials

For an uncertain system described by convex combination of interval polynomials, its Hurwitz-stability can be guaranteed by certain subset composed of vertices and edges. Furthermore, the testing set does not increase when the order of given system increases. Foundation item: the National Natural Science Foundation of China (69925307) Biography: WANG Zhi-zhen (1974-)     

不确定初始几何缺陷壳的动态屈曲分析

研究了动载作用下含有不确定初始几何缺陷壳的动态屈曲问题.给出了不确定初始几何缺陷在椭球描述和区间描述下,壳的动态响应和安全因子的界限估计.从数学证明和数值算例两方面,讨论了基于不确定性初始几何缺陷两种描述的凸模型方法和区间分析方法所得壳动态响应及安全因子的界限的关系,为判断壳的动态屈曲失效提供了依据.