颗粒物质在黏性流体中运动过程的数值分析

为了研究颗粒物质在黏性流体中整个下落过程的运动规律，对不同半径的颗粒物体运动过程进行了数值模拟分析。在给出颗粒物质下落过程模型的基础上，考虑黏性流体对颗粒物质的黏滞阻力，并通过受力分析建立了下落运动微分方程；利用计算机数值分析了不同半径的颗粒物质在同一种黏性流体中下落速度、入潜深度随时间变化的情况，并分析了不同半径的颗粒物质下落达到的稳态速度与所消耗时间的变化关系。计算结果表明：当颗粒物质半径 r<1mm 时，其在黏性流体中的下落距离与时间近似呈线性关系；随着颗粒物质尺寸的增加，其下落距离会呈现非线性增加；当颗粒物质半径 r>2mm 时，其达到稳态速度后的入潜深度与颗粒物质的半径呈非线性变化，且颗粒物质越大，达到稳态速度后的入潜深度越深；颗粒物质在黏性流体中下落后达到的稳态速度与颗粒物质的半径并非呈线性关系，且颗粒物质越小，达到的稳态速度越小，与黏性流体的深度无关，其相应入潜深度与所需时间呈非线性关系。     

NUMERICAL SIMULATION OF DRAG REDUCTION CAUSED BY PATTERNED WALL

用加入近邻相互作用势的格子玻尔兹曼方法模拟了表面修饰了去润湿性颗粒的二维管道内流体的流动. 结果显示,管壁的粗糙度、修饰颗粒的润湿性、管内压强以及管内流体黏滞性等都是影响管内流体流速的因素. 在一定压强下,管壁修饰疏水性的颗粒可以明显减小阻力,模拟结果也解释了开采地下油藏时存在启动压力梯度的原因.     

环形激波聚焦流场特性的数值研究

针对环形激波聚焦过程产生的高温、高压特性,采用间断有限元方法模拟了环形激波在同轴圆柱 形激波管内的聚焦流场特性。计算结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的二次 激波、涡环、三波交点和球面双马赫反射等主要流动特征。此外,通过改变环形管道内外半径对聚焦流场进行 模拟发现,环形管道外径对中心轴线上聚焦峰值压力的大小和位置影响较小,环形管道内径对中心轴线上聚 焦峰值压力的大小和位置影响较大。计算结果可以为工程应用提供一定的理论指导。     

柔性圆柱涡激振动流体力系数识别及其特性

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前,海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据,存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用,建立了柔性圆柱涡激振动流体力模型,运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数,设计并开展了拖曳水池模型实验,实验用柔性圆柱模型的质量比为1.82,长径比为195.5.通过与刚性圆柱流体力系数对比,深入分析了柔性圆柱流体力系数的特性.结果表明:柔性圆柱在一阶模态控制区,流体力系数随约化速度变化趋势与刚性圆柱大致相似;二阶模态控制区,升力系数和脉动阻力系数显著增大;附加质量系数在响应频率较低时与振动位移的相关性增强;当响应频率较低时,振动位移较大区域为能量耗散区,当响应频率较高时,振动位移较大区域为能量输入区.     

微型管中液滴流动的三维数值模拟

对于微型设备中的低雷诺数流动,毛细力和黏性力起主导作用. 应用相场方法,引 入自由能泛函,研究了二相流体在微型管中流动问题及表面浸润现象,并给出了微型管中二 相流体的无量纲输运方程. 针对方形微管道,利用差分法给出了输运方程的数值求解方法. 最后,模拟了方形直管中的液滴流动和变形的过程,并给出了液滴前后压力差与其它主要物 理参数之间的变化关系. 结果表明,压力差随液滴半径增大而增加,而随毛细管系数的增大 而减小.     

圆柱薄壳稳定性的一个修正理论

著名的唐乃尔(Donnell)——穆什塔利的简化壳体理论只能较精确地适用于较短圆柱壳稳定性计算.其近似性误差随长度与半径之比的增加而增大.本文考虑了横向切力的影响,对非完善型圆柱壳体推导了几何非线性理论的基本方程,建立了对各种长度半径比的圆柱壳体稳定性均适用的修正理论.     

具有黏滞阻力时的最速降线

针对实际应用中存在黏滞阻力的最速降线的问题, 首先推导出适于此类问题的解除约束的 广义变分原理, 它适用于具有摩擦阻尼和多自由度系统优化的问题. 得到描述有黏滞阻力情 况下最速降线相关函数的微分方程, 它在黏滞阻力为零时即退化为滚轮线. 利用MATLAB数值 计算给出了最速降线受黏滞阻力的影响: 在黏滞阻力系数较小时最速降线趋于变凹, 当阻力 系数增大到一定值之后最速降线趋于平缓, 当阻力系数很大时最速降线趋于直线.     

一种滞弹簧耗能的新型离散元滚动阻力模型研究

颗粒间滚动阻力对颗粒体系的稳定性起着重要作用. 在传统的离散元法中, 滚动阻力模型通常由转动弹簧、转动黏壶和摩擦元件表达, 颗粒滚动动能由黏滞力(矩)和摩擦力做功耗散. 由于黏滞力(矩)与滚动速度相关, 临近静止状态的颗粒滚动速度变小, 动能耗散减弱, 传统的离散元模拟得到颗粒由滚动到静止耗费的时间比试验观测的结果要长. 为解决这一问题, 基于摩擦学理论分析了滚动阻力产生的材料滞弹性机理, 将其引入离散元滚动阻力模型, 提出了一种速度无关型动能耗散的滞弹簧, 给出了滞弹簧的弹性恢复力计算公式, 建立了一种新型的离散元滞弹性滚动阻力模型(HDEM). 为验证新型滚动阻力模型的正确性, 通过一个光学物理试验对单个圆形颗粒试件的自由滚动过程进行了测量, 将测量数据与新型的滞弹型离散元模型和传统离散元模型计算结果进行了对比. 结果显示, 基于滞弹性滚动阻力模型HDEM计算结果与试验数据吻合程度更高, 而且模拟得到的颗粒摆动频率更符合试验现象.      

一种新的高阶弹簧-阻尼-质量边界——-无限域圆柱对称波动问题

提出一种描述力-位移时间卷积关系的高阶弹簧-阻尼-质量模型,并将其作为人工边界条件直接应用于弹性动力学无限域圆柱对称运动问题的时域数值求解. 该人工边界条件不存在旁轴近似、多次透射等位移型外推人工边界条件普遍存在的高、低频失稳问题;与黏性、黏弹性边界等应力型人工边界条件相比,它具有高阶精度,且是严格高、低频双渐近的,也可以退化到黏性、黏弹性边界;该边界可以像黏性、黏弹性边界一样利用商用有限元软件中内置的并联弹簧-阻尼器、质量单元和时间积分求解器在商用软件中方便地实现,便于研究人员和工程师应用. 分析的几个简单数值算例也验证了该边界条件的上述优点.      

双圆柱绕流特性的模拟研究

采用格子Boltzmann方法对低雷诺数下气体绕流圆柱的规律进行了研究. 对比计算了双圆柱在不同圆心距、不同Re数、不同来流速度与双圆柱圆心连线角度的情况下,各个圆柱的受力大小和曳力系数. 结果表明,若Re数为20, 改变圆柱间距,圆柱间距在1.2d和1.4d之间时,下游圆柱所受曳力有极小值;双圆柱间距为1.6d时,双圆柱受到总曳力最小;圆柱间距大于2d时,上游颗粒受到的曳力不再受到下游颗粒的影响. 若圆柱间距为1.2d, 改变雷诺数,Re数在30和40之间,下游圆柱所受曳力有极小值. 另外,来流速度角度对圆柱的受力影响很大. 上述规律为低Re数下圆柱绕流的深入研究与应用打下基础.      

Isothermal flow past a blowing sphere

Steady, axisymmetric, isothermal, incompressible flow past a sphere with uniform blowing out of the surface is investigated for Reynolds numbers in the range 1 to 100 and surface velocities up to 10 times the free stream value. A stream-function-velocity formulation of the flow equations in spherical polar co-ordinates is used and the equations are solved by a Galerkin finite-element method. Reductions in the drag coefficients arising from blowing are computed and the effects on the viscous and pressure contributions to the drag considered. Changes in the surface pressure, surface vorticity and flow patterns for two values of the Reynolds number (1 and 40) are examined in greater detail. Particular attention is paid to the perturbation to the flow field far from the sphere.     

Cell models for micropolar flow past a viscous fluid sphere

The Stokes axisymmetrical flow of an incompressible micropolar fluid past a viscous fluid sphere and the flow of a viscous fluid past a micropolar fluid sphere are investigated. The appropriate boundary conditions are taken on the surface of the sphere, while the proper conditions applied on the fictitious boundary of the fluid envelope vary depending on the kind of cell-model. These problems are solved separately in an analytical fashion, and the velocity profile and the pressure distribution inside and outside of the droplet are shown in several graphs for different values of the parameters. Numerical results for the normalized hydrodynamic drag force acting, in each case, on the spherical droplet-in-cell are obtained for various values of the parameters representing volume fraction, the classical relative viscosity, the micropolarity and spin parameters are presented both in tabular and graphical forms. Results of the drag force are compared with the previous particular cases.     

A general formula for the drag on a sphere placed in a creeping unsteady micropolar fluid flow

In the present work, we investigate the creeping unsteady motion of an infinite micropolar fluid flow past a fixed sphere. The technique of Laplace transform is used. The drag formula is obtained in the physical domain analytically by using the complex inversion formula of the Laplace transform. The well known formula of Basset for the drag on a sphere placed in an unsteady viscous fluid flow and that of Ramkissoon and Majumdar for steady motion in the case of micropolar fluids are recovered as special cases. The obtained formula is employed to calculate the drag force for some micropolar fluid flows. Numerical results are obtained and represented graphically.     

The inlet effect on the drag factor of a sphere in a tube

The creeping motion Ground a sphere situated axisymmetrically near the entrance of asemi-infinite circular cylindrical tube is analyzed using infinite series solutions for thevelocity components. pressure and the stream function. Truncating the infinite series. thecorresponding coefficients in the series are determined by a collocation technique. The dragfactor and the stress distribution on the surface of the sphere are calculated for the sphere inmotion in quiescent fluid and for the flow with uniform velocity at the entrance past a rigidlyheld sphere. The results indicate that a sphere near the entrance which has a uniformentrance velocity profile will suffer larger drag than that in infinite tube.Theconvergence of the collocation technique is tested by numerical calculation. It is shown thatthe technique has good convergence properties.     

Drag of a rotating sphere

We consider the problem of steady incompressible viscous fluid flow about a rotating sphere, with the flow specified on a sphere of finite radius, which reduces to the solution of the complete Navier-Stokes equations.The dimensionless stream functions and circulai velocity are sought in the form of series in powers of the Reynolds numbers, which converge for small values of this number. Recurrence formulas are derived for determining the coefficients of these series. The pressure, rotational resistance torque, and drag are determined. It is established that the rotating sphere has higher drag than a stationary sphere. The leading term of the series in powers of the Reynolds number for the drag and resistive torque is calculated.     

Slow motion of a permeable sphere in a viscous fluid

Flow of a viscous fluid past a permeable sphere is investigated in the Stokes approximation. An example of such a flow is flow past a perforated or meshed spherical surface. The elements of the sphere contain rigid impermeable sections and openings through which the fluid can flow. The interaction of the sphere with the flow is described by two drag coefficients, which established the connection between the flow velocity of the fluid at the sphere and the stress tensor on it. The dependence of the flow pattern and also the drag and flow rate of the fluid on these coefficients is investigated. In special cases, the obtained solution describes flow past solid and liquid spheres.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 165–167, September–October, 1982.     

Hyperbolic Approximation of the Navier-Stokes Equations for Viscous Mixed Flows

Simplified two-dimensional Navier-Stokes equations of the hyperbolic type are derived for viscous mixed (with transition through the sonic velocity) internal and external flows as a result of a special splitting of the pressure gradient in the predominant flow direction into hyperbolic and elliptic components. The application of these equations is illustrated with reference to the calculation of Laval nozzle flows and the problem of supersonic flow past blunt bodies. The hyperbolic approximation obtained adequately describes the interaction between the stream and surfaces for internal and external flows and can be used over a wide Mach number range at moderate and high Reynolds numbers. Examples of the calculation of viscous mixed flows in a Laval nozzle with large longitudinal throat curvature and in a shock layer in the neighborhood of a sphere and a large-aspect-ratio hemisphere-cylinder are given. The problem of determining the drag coefficient of cold and hot spheres is solved in a new formulation for supersonic air flow over a wide range of Reynolds numbers. In the case of low and moderate Reynolds numbers a drag reduction effect is detected when the surface of the sphere is cooled.     

Lift force on rotating sphere at low Reynolds numbers and high rotational speeds

The lift force on an isolated rotating sphere in a uniform flow was investigated by means of a three-dimensional numerical simulation for low Reynolds numbers (based on the sphere diameter) (Re&lt;68.4) and high dimensionless rotational speeds (Г5). The Navier-Stokes equations in Cartesian coordinate system were solved using a finite volume formulation based on SIMPLE procedure. The accuracy of the numerical simulation was tested through a comparison with available theoretical, numerical and experimental results at low Reynolds numbers, and it was found that they were in close agreement under the above mentioned ranges of the Reynolds number and rotational speed. From a detailed computation of the flow field around a rotational sphere in extended ranges of the Reynolds number and rotational speed, the results show that, with increasing the rotational speed or decreasing the Reynolds number, the lift coefficient increases. An empirical equation more accurate than those obtained by previous studies was obtained to describe both effects of the rotational speed and Reynolds number on the lift force on a sphere. It was found in calcttlations that the drag coefficient is not significantly affected by the rotation of the sphere. The ratio of the lift force to the drag force, both of which act on a sphere in a uniform flow at the same time, was investigated. For a small spherical particle such as one of about 100μm in diameter, even if the rotational speed reaches about 10^6 revolutions per minute, the lift force can be neglected as compared with the drag force.     

Drag and Separation Flow Past Solid Sphere with Porous Shell at Moderate Reynolds Number

Axisymmetric viscous, two-dimensional steady and incompressible fluid flow past a solid sphere with porous shell at moderate Reynolds numbers is investigated numerically. There are two dimensionless parameters that govern the flow in this study: the Reynolds number based on the free stream fluid velocity and the diameter of the solid core, and the ratio of the porous shell thickness to the square root of its permeability. The flow in the free fluid region outside the shell is governed by the Navier–Stokes equation. The flow within the porous annulus region of the shell is governed by a Darcy model. Using a commercially available computational fluid dynamics (CFD) package, drag coefficient and separation angle have been computed for flow past a solid sphere with a porous shell for Reynolds numbers of 50, 100, and 200, and for porous parameter in the range of 0.025–2.5. In all simulation cases, the ratio of b/a was fixed at 1.5; i.e., the ratio of outer shell radius to the inner core radius. A parametric equation relating the drag coefficient and separation point with the Reynolds number and porosity parameter were obtained by multiple linear regression. In the limit of very high permeability, the computed drag coefficient as well as the separation angle approaches that for a solid sphere of radius a, as expected. In the limit of very low permeability, the computed total drag coefficient approaches that for a solid sphere of radius b, as expected. The simulation results are presented in terms of viscous drag coefficient, separation angles and total drag coefficient. It was found that the total drag coefficient around the solid sphere as well as the separation angle are strongly governed by the porous shell permeability as well as the Reynolds number. The separation point shifts toward the rear stagnation point as the shell permeability is increased. Separation angle and drag coefficient for the special case of a solid sphere of radius r = a was found to be in good agreement with previous experimental results and with the standard drag curve.     

Motion of a permeable shell in a spherical container filled with non-Newtonian fluid

This paper presents an analytical study of creeping motion of a permeable sphere in a spherical container filled with a micro-polar fluid. The drag experienced by the permeable sphere when it passes through the center of the spherical container is studied.Stream function solutions for the flow fields are obtained in terms of modified Bessel functions and Gegenbauer functions. The pressure fields, the micro-rotation components,the drag experienced by a permeable sphere, the wall correction factor, and the flow rate through the permeable surface are obtained for the frictionless impermeable spherical container and the zero shear stress at the impermeable spherical container. Variations of the drag force and the wall correction factor with respect to different fluid parameters are studied. It is observed that the drag force, the wall correction factor, and the flow rate are greater for the frictionless impermeable spherical container than the zero shear stress at the impermeable spherical container. Several cases of interest are deduced from the present analysis.