入射涡与圆柱相互作用的数值模拟

本文用快速涡方法对入射涡与圆柱的相互作用进行了数值模拟,观察到了入射涡在圆柱表面上诱导的二次分离和三次分离现象。二次涡的产生,与入射涡配对,改变了它们的运动轨迹。二次涡是入射涡“回跳”现象的主要原因。本文还对不同入射涡强度及相互位置作了计算,并分析了不同参数对涡运动轨迹的影响。这些现象与涡的无粘圆柱绕流有着本质的差异。     

Gas Flow in Porous Media With Klinkenberg Effects

Gas flow in porous media differs from liquid flow because of the large gas compressibility and pressure-dependent effective permeability. The latter effect, named after Klinkenberg, may have significant impact on gas flow behavior, especially in low permeability media, but it has been ignored in most of the previous studies because of the mathematical difficulty in handling the additional nonlinear term in the gas flow governing equation. This paper presents a set of new analytical solutions developed for analyzing steady-state and transient gas flow through porous media including Klinkenberg effects. The analytical solutions are obtained using a new form of gas flow governing equation that incorporates the Klinkenberg effect. Additional analytical solutions for one-, two- and three-dimensional gas flow in porous media could be readily derived by the following solution procedures in this paper. Furthermore, the validity of the conventional assumption used for linearizing the gas flow equation has been examined. A generally applicable procedure has been developed for accurate evaluation of the analytical solutions which use a linearized diffusivity for transient gas flow. As application examples, the new analytical solutions have been used to verify numerical solutions, and to design new laboratory and field testing techniques to determine the Klinkenberg parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs from The Geysers geothermal field. We show that this new approach and the traditional method of Klinkenberg yield similar results of Klinkenberg constants for the laboratory tests; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

粘弹性厚壁筒问题的辛本征解方法

将哈密顿体系引进到粘弹性力学厚壁筒问题中,在辛体系下重新描述了基本问题,即建立了正则方程组。借助于积分变换,得到了拉伸、扭转和弯曲等问题的解以及有边界局部效应的解。将原问题归结为辛几何空间中的零本征值本征解和非零本征值本征解问题,从而建立了一种有效的分析问题方法和数值方法。为解决同类问题提供了一条可行的路径。     

Modeling the effect of diffusion of the ambient medium on the long-term strength of a hollow cylinder under uniaxial tension

The problem of the long-term strength of an extended thick-walled tube containing a corrosive medium in the internal cavity is solved. The diffusion of this medium into the tube material is analyzed. The diffusion equation is solved approximately by introducing the diffusion front, and the error of the solution is estimated. The dependence of the time of fracture of the tube on the variable tensile stress and the concentration of the medium filling the cavity is obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 4, pp. 88–93, July–August, 2007.     

Modeling fluid and heat transport in the reactive, porous bed of downdraft (biomass) gasifier

A fluid flow and heat transfer model has been developed for the reactive, porous bed of the biomass gasifier to simulate pressure drop, temperature profile in the bed and flow rates. The conservation equations, momentum equation and energy equation are used to describe fluid and heat transport in porous gasifier bed. The model accounted for drag at wall, and the effect of radial as well as axial variation in bed porosity to predict pressure drop in bed. Heat transfer has been modeled using effective thermal conductivity approach. Model predictions are validated against the experiments, while effective thermal conductivity values are tested qualitatively using models available in literature. Parametric analysis has been carried out to investigate the effect of various parameters on bed temperature profile and pressure drop through the gasifier. The temperature profile is found to be very sensitive to gas flow rate, and heat generation in oxidation zone, while high bed temperature, gas flow rate and the reduction in feedstock particle size are found to cause a marked increase in pressure drop through the gasifier. The temperatures of the down stream zones are more sensitive to any change in heat generation in the bed as compared to upstream zone. Author recommends that the size of preheating zone may be extended up to pyrolysis zone in order to enhance preheating of input air, while thermal insulation should not be less than 15 cm.     

CH4-O2 混合气中爆燃爆震转捩的数值模拟

运用化学流体力学基本理论和两步燃烧反应模型原理，建立了一维封闭体系可燃气爆燃爆震转变现象的数学模型，利用拉格朗日质量坐标变换下的Ｌａｘ－Ｗｅｎｄｒｏｆ和Ｍｅｃｏｒｍｉｃ气动差分与Ａｄａｍｓ化学差分格式，求解基本方程，成功地完成了过程的数值模拟，清楚地说明了可燃气中ＤＤＴ现象由压缩波到激波达到稳态爆震的发生机制和火焰带引生爆震波的过程行为。     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

Further study on vortex shedding flow by a circular cylinder

In the present paper the mechanism involved in vortex shedding flows is investigated in detail. In the early stage of the unsteady separated flow the interaction between secondary vortices and primary vortices is quite strong. In the later stage of the flow, corresponding to the vortex shedding the recirculating flow region on each side of the aft body goes through such a cycle: growth-contraction-growth, the secondary separation occurs periodically rather than continuously. The reduction of circulation is taken into account in three cases with different decay factors to study its influence on the prediction of main flow characteristics. Results show that to simulate vortex shedding flow it is necessary to include the reduction of circulation to bring the calculated results into good agreement with experiments. An improved discrete vortex model is suggested in which both the secondary separation and the reduction are incorporated. The processes of vortex shedding, the forces prediction and other flow characteristics are given and some discussions are made. Porject is supported by National Natural Science Foundation of China.     

静压气体球轴承支承球形转子的干扰力矩分析

中心小孔供气单向受载球面气体轴承可以用于球形转子的静平衡测量。对作用于转子上的干扰力矩进行估算是平衡装置设计的重要部分。本在一定假设条件下推导了由粘性剪切应力和气膜支撑力引起的作用于转子上的干扰力矩。干扰力矩以轴承包角、中心气室张角、气膜压力、转子转速、转子旋转轴位置、失中度、转子非球表示。以干扰力矩最小为准则分析了这些参数的影响，结论有助于静平衡装置的优化设计以及对精度的进一步分析。     

Elastodynamic solution of a non-homogeneous orthotropic hollow cylinder

A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions and the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to homogeneous ones. By virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally obtained, which can be degenerated in a rather straightforward way into the solution for a homogeneous orthotropic hollow cylinder and isotropic solid cylinder as well as that for a non-homogeneous isotropic hollow cylinder. Using the present method, integral transform can be avoided and it can be used for hollow cylinders with arbitrary thickness and subjected to arbitrary dynamic loads. Numerical results are presented for a non-homogeneous orthotropic hollow cylinder subjected to dynamic internal pressure. The project supported by the National Natural Science Foundation of China (10172075 and 10002016)