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

盐岩储库区地面沉降预测与控制研究现状与展望

利用盐岩地层建设地下油气储库群已成为各国发展能源储备的重要方向,然而关于储库区地面沉降的研究至今仍处于探索阶段。由于地面沉降所造成的重大事故屡见不鲜,因此有必要集中力量研究盐岩库区地面沉降问题,为我国盐岩储库群的建设提供可靠的理论依据。本文主要从以下3个方面总结了近几十年来国内外学者在盐岩储库区地面沉降预测方面所做出的努力与成绩： （1）现场监测——法国、德国、美国先后建立了现场监测网络,获得了大量的长期监测数据;  （2）理论计算——我国学者于20世纪90年代末将随机介质理论引入到盐矿水溶开采地面沉降预测中,从而为储库区地面沉降预测提供了新的理论基础;  （3）数值模拟——欧洲学者利用FLAC2D 程序对盐岩储库区的地面稳定性进行了数值模拟,而我国学者则利用在新概率积分法的基础上建立起的预测模型成功进行了沉降预测模拟。此外,为了有效地控制库区地面变形,预留矿柱式水溶开采法在国外众多盐岩地下储库群建设中得以广泛应用。本文最后提出在后续研究中将重点着眼于地表变形随时间的发展过程,并对现有理论模型进行适当的修正以及利用现场实测数据进一步完善预测所需参数的选取方法,从而使预测结果更为准确可靠。     

Spontaneous Switching between Permeability Enhancement and Degradation in Fractures in Carbonate: Lumped Parameter Representation of Mechanically- and Chemically-Mediated Dissolution

Principal mechanical and chemical processes contributing to the observed spontaneous switching from net decrease in permeability to net increase in a fracture in carbonate are examined. The evolution of permeability, and related fracture aperture, is represented through a lumped parameter model. The significant processes of pressure solution beneath bridging asperities, transport of dissolved mass to the fracture void, and subsequent precipitation or dissolution within the fracture void enable the principal characteristics of observed behavior to be followed. The evolution of dissolved mass concentration in the pore fluid is followed for arbitrary applied stress, temperature, and pH conditions, with appropriate feedback to the evolution of fracture permeability. Comparisons with experimental measurements in limestone (Polak et al., 2004, Water Resour. Res. Vol. 40, W03502, doi: 10.1029/2003GL017575) show satisfactory agreement for the evolution of fracture aperture and to a lesser degree in calcium concentrations in the effluent pore fluid. Importantly, the spontaneous switching in permeability change, from aperture reducing to aperture increasing, with no change in environmental conditions, is replicated without the need for an ad hoc trigger. Although this switch is accurately replicated, the lumped parameter model is incapable of replicating the rapid observed growth in permeability that directly follows. This inability results from the assumed form of the lumped asperity model, that is incapable of representing the spatially distributed change in aperture that is seen to occur within the fracture. Despite this inconsistency, the model is shown capable of representing the principal behaviors evident in the response.     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

Investigation of power-type variational principles in liquid-filled system

Starting from the basic equations of hydrodynamics, the maximum powertype variational principle of the hydrodynamics of viscous fluids was established by Weizang CHIEN in 1984. Through long-term research, it is clarified that the maximum power-type variational principle coincides with the Jourdian principle, which is one of the common principles for analytical mechanics. In the paper, the power-type variational principle is extended to rigid-body dynamics, elasto-dynamics, and rigid-elastic-liquid coupling dynamics. The governing equations of the rigid-elastic-liquid coupling dynamics in the liquid-filled system are obtained by deriving the stationary value conditions. The results show that, with the power-type variational principles studied directly in the state space, some transformations in the time domain space may be omitted in the establishing process, and the rigid-elastic-liquid coupling dynamics can be easily numerically modeled. Moreover, the analysis of the coupling dynamics in the liquid-filled system in this paper agrees well with the numerical analyses of the coupling dynamics in the liquid-filled system offered in the literatures.     

On the spatial localization of the laminar boundary layer in a non-Newtonian dilatant fluid

In dilatant fluids the shear perturbation propagation rate is finite, in contrast to Newtonian and pseudoplastic fluids in which it is infinite [1]. Therefore, in certain dilatant fluid flows, frontal surfaces separating regions with zero and nonzero shear perturbations may be formed. Since, in a sense, the boundary layer is a “time scan” of the nonstationary shear perturbation propagation process, in dilatant fluids the boundary layer should definitely be spatially localized. This was first mentioned in [2] where, however, it was mistakenly asserted that boundary layer spatial localization does not take place in all dilatant fluids and is absent in so-called “hardening” dilatant fluids. In [3], the solutions of the laminar boundary-layer equations for speudoplastic and “hardening” dilatant fluids were investigated qualitatively. The formation of frontal surfaces in dilatant fluid flows is usually mathematically related with the existence of singular solutions of the corresponding differential equations [4]. However, since the analysis performed in [3] was inaccurate, in that study singular solutions were not found and it was incorrectly concluded that in “hardening” dilatant fluids there is no spatial boundary layer localization. The investigation performed in [5] showed that in fact in “hardening” dilatant fluids boundary layers are spatially localized, since there exist singular solutions of the corresponding differential equations. Subsequently, this result was reproduced in [6], where an attempt was also made to carry out a qualitative investigation of the solutions of the laminar boundary-layer equations for other types of dilatant fluids. The author did not find singular solutions in this case and mistakenly concluded that in these fluids there is no spatial boundary layer localization. This misunderstanding was due to the fact that in [6] it was not understood that in dilatant fluid flows the formation of frontal surfaces can be mathematically described not only in relation to the existence of singular solutions.     

具有弱夹层的封闭油藏中的瞬态渗流

多层油藏分为层间无窜流的油藏和层间有窜流的油藏。本文研究了层间有弱渗透夹层的两层油藏中两种情形的瞬态渗流,包括单层打开的瞬态渗流和一层注入、一层采出的两层油藏的瞬态渗流。用面积平均方法求出了弱渗透夹层封闭油藏中瞬态渗流的平均压降。分析表明:在渗流前期,层间窜流量逐渐增大,在后期层间窜流量趋于稳定;在单层开采的情况下:初期两层压降不同,非打开层的压降有滞后现象,并且隔层渗透率越小,油井工作强度越弱,滞后时间越长;后期两层压降趋于一致,可以作为一层油藏处理。在一注一采的油藏中,当注采量相同时,两层压降在后期趋于不同的稳定值;当注采量不同时,两层压降在后期趋于一致,但不能达到稳定。     

岩石隧道围岩变形时空效应分析

岩石隧道围岩变形具有时空效应特征。根据围岩变形速率,岩石隧道围岩变形一般可划分为3个阶段,即急剧变形阶段、稳定变形阶段和流变阶段。通过总结分析围岩变形3阶段的特点,结合中梁山隧道D-5H量测剖面的实测数据,对围岩变形的空间效应和时间效应进行了分析。空间效应集中发生在急剧变形段,空间效应段主要靠围岩自身以及初次支护克服围岩发生破坏变形,时间效应则主要体现在流变段。以华蓥山隧道等76个隧道实例为统计样本,分别对围岩变形时空效应与围岩类别和塌方事故的关系进行了相关性分析。结果表明:80％以上的塌方发生于急剧变形段,13％发生在稳定变形段,只有7％左右的塌方发生在流变段。其中Ⅳ类和Ⅴ类围岩在3个阶段都可能发生塌方,Ⅲ类围岩则很少在流变段发生塌方。Ⅰ类和Ⅱ类围岩则基本不会发生较大规模的塌方。对深入了解隧道围岩的变形规律,为隧道灾害防治、选择恰当的支护时机和支护方式很有意义。     

Flow and heat transfer over a backward-facing step with a cylinder mounted near its top corner

A study was made to see if it is possible to enhance the heat transfer in the downstream region of a backward-facing step, where heat transfer is normally deteriorated, by the insertion of a cylinder near the top corner of the step. Cylinder size and streamwise position of the cylinder were kept constant but the cross-stream position of the cylinder was changed in three steps. Results of the heat transfer experiment, flow visualization, and measurement of the averaged and fluctuating flow fields were reported. When the cylinder was mounted at a position, a little higher than the top surface of the step, a jet-like flow pattern emerged in the averaged velocity profile beneath the cylinder and the recirculating flow was intensified. Therefore, the velocity of recirculating flow near the wall is increased at some streamwise positions. Additionally, the velocity fluctuation was intensified not only in the shear layer between the jet-like flow and the recirculating flow regions but also in the near wall region, resulting in the effective augmentation of heat transfer in this case. Therefore, it is concluded that the mounting of a cylinder is effective in the enhancement of deteriorated heat transfer in the recirculating flow region, if its is mounted in a proper position.     

延性动态损伤的细观描述

对动载荷下材料的延性动态损伤的细观力学研究现状进行了评述，详细讨论了现有的几类典型的动态延性损伤的微孔生长模型。针对目前该领域工作中存在的问题，指出了今后应开展的工作方向。     

人体生物摩擦学的研究现状与展望

介绍了人体生物摩擦学的概念,分析了头发、眼睛、皮肤、牙齿和关节等系统的生物摩擦学研究现状,就血液、人工心脏泵、人工心脏瓣膜和食品味觉的生物摩擦学研究方向进行了展望.指出当前人体生物摩擦学研究的关键科学问题主要包括人体摩擦副的特殊摩擦学行为及机理、模拟人体生物环境条件的人工器官摩擦磨损试验方法研究、人工材料在体内的生物摩擦学行为反应预测及磨粒在体内的生物反应和临床病症的评定.