Interaction between two spheres placed in tandem arrangement in steady and pulsating flow

The interaction among two spheres in tandem formation are studied for a Reynolds number of 300 using both steady and pulsating inflow conditions. The purpose is to further investigate the force characteristics as well as the shedding patterns of the two spheres as the separation distance is changed from 1.5 to 12 sphere diameters. The method used for the simulations is the volume of solid (VOS) method, an approach based on the volume of fluid (VOF) method. Comparisons with other computational methods have shown VOS to accurately resolve the flow field around solid spheres. The results show that the separation distance plays a significant role in changing the flow patterns and shedding frequencies at moderate separation distances, whereas effect on drag is observed even at a separation distance of 12 diameters.     

Numerical study of the axially symmetric motion of an incompressible viscous fluid in an annulus between two concentric rotating spheres

The research reported herein involved the study of the transient motion of a system consisting of an incompressible Newtonian fluid in an annulus between two concentric, rotating, rigid spheres. The primary purpose of the research was to study the use of a numerical method for analysing the transient motion that results from the interaction between the fluid in the annulus and the spheres which are started suddenly by the action of prescribed torques. The problems considered in this research included cases where: (a) one or both spheres rotate with prescribed constant angular velocities and (b) one sphere rotates due to the action of an applied constant or impulsive t?orque. In this research the coupled solid and fluid equations were solved numerically by employing the finite difference technique. With the approach adopted in this research, only the derivatives with respect to spatial variables were approximated with the use of the finite difference formulae. The steady state problem was also solved as a separate problem (for verification purposes), and the results were compared with those obtained from the solution of the transient problem. Newton's algorithm was employed to solve the algebraic equations which resulted from the steady state problem, and the Adams fourth-order predictor–corrector method was employed to solve the ordinary differential equations for the transient problem. Results were obtained for the streamfunction, circumferential function, angular velocity of the spheres and viscous torques acting on the spheres as a function of time for various values of the system dimensionless parameters.     

旋转圆球入水空泡特性与流场结构的大涡模拟研究

圆球旋转入水过程对于基于先导物投放的新型入水降载方式具有重要研究价值. 采用大涡模拟方法结合均质多相流模型和VOF界面捕捉算法, 对低弗劳德数条件下疏水圆球高速旋转入水的自由运动过程进行了数值模拟, 研究了转速对入水空泡演化、流场结构和水动力特性的影响. 采用动网格与滑移网格技术实现圆球的自由运动, 并基于试验结果对比验证了数值模拟的可靠性与正确性. 旋转运动的升力效应导致圆球入水弹道发生偏转并从水面携带横向楔形射流侵入空泡内部. 采用入水砰击速度与转速进行归一化分析, 结果表明入水转速的增加显著改变了圆球的动力特性: 水平方向的速度和位移以及升力峰值都随入水转速的增加而变大, 但升力峰值受到入水速度的限制; 而垂直方向的速度和加速度以及空泡断裂深度几乎不受转速增加的影响, 并且空泡深闭合发生前圆球转速变化不大. 入水转速的增加也使液面飞溅环和空泡断裂的非对称性增强, 在较低转速时发生空泡表面闭合, 而在较高转速时则发生空泡深闭合. 对于空泡深闭合模式, 入水转速的增加带来更强的横向楔形射流, 并且抑制了空泡断裂产生的高压以及相应涡结构的生成, 致使圆球在入水砰击阶段承受更低的侧向压力.      

涡激诱导并列双圆柱碰撞数值模拟研究

圆柱类结构物的涡激振动是工程中较为常见的一种现象,如果圆柱结构物之间的距离较小, 就会产生涡激诱导碰撞现象,而涡激碰撞会比涡激振动对结构物疲劳破坏产生更严重的威胁.采用浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时出现数值求解不稳定问题,采用有限元方法对圆柱的运动和碰撞进行求解,通过数据回归方法建立了流体流动条件下的润滑模型,对不同间隙比下涡激诱导并列双圆柱振动及碰撞过程进行了数值模拟, 数值结果表明,如果两圆柱产生了碰撞将会有连续的碰撞发生, 碰撞时出现了多阶频率,振动主频率要比无碰撞时大, 两圆柱碰撞时的相对速度比自由来流速度小;当两圆柱相互接近时, 随着涡环分离角度的逐渐倾斜, 横向流体力先逐渐减小,当两圆柱间涡环开始相互影响发生挤压时, 横向流体力开始逐渐增大;当两圆柱开始反弹时, 两圆柱间形成了低压区, 改变了横向流体阻力的方向,使两圆柱又产生了接近运动,如此反复从而产生了碰撞后横向流体力和圆柱速度的振荡现象.      

三峡水库岸边排污的特性及数值模拟研究

数值模拟是研究和解决三峡水库污水岸边排污的重要手段之一．针对三峡水库岸边污水排放的特点,分析并阐述了排污口附近的水流和污染物分布特性及其所带来的数值困难．基于这些分析,讨论了深度平均二维模型和不同分层三维模型在模拟三峡库区岸边排污问题中的优势和不足,为库区岸边侧向排污数值模型的进一步发展奠定了基础．     

下倾管-立管水气严重段塞流数值模拟

针对海洋油气传输中常见的下倾管-立管系统, 采用Brackbill模型模拟气液相界面间表面张力, VOF方法追踪气液两相运动界面, 提出了管内气液两相流数值模拟方法. 在低气液相进口折算速度下, 数值模拟了该种管型下的严重段塞流动现象, 分析了相关物理参数的变化特性. 结果表明, 在严重段塞流下, 管内流型流态、压力、液塞运动速度、立管出口气液相平均速度、下倾管及立管内含气率等均具有明显周期性特征, 而且一个周期内严重段塞流可分为4个阶段, 进而给出了各阶段中相关参数的变化特性. 数值模拟结果与相关文献中的实验结果吻合良好,表明了该数值模拟方法的有效性.      

A capillary force model for interactions between two spheres

Based on the method of energy principle, an analytical approach for computing the capillary force for sphere/sphere geometry is presented in this paper. In modeling the capillary force, we consider spheres with both equal and non-equal radii, for both symmetric and asymmetric configurations at liquid/solid interfaces. We use numerical analysis to investigate the validity and efficiency of the derived model. The effect of various parameters including humidity, distance between two spheres, radii of spheres and contact angles on the meniscus force are investigated. Finally the results obtained from the model are compared with experimental measurements, and the accuracy and precision of the presented approach is verified.     

Squeeze flow of a second-order fluid between two parallel disks or two spheres

The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds‘ lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neelected.     

近海水域海洋流体动力学方程的数值模拟

近海水域海洋流体动力学方程数值模拟的研究,最早可追溯到1919年Defant的工作,50年代末开始得到蓬勃发展。60年代和70年代主要发展二维全流模式。现在对于二维全流模式数值模拟的研究应该说是相当完善的,不仅有系统的理论,提出了许多有效的方法,而且广泛地应用于诸如潮汐环流、风暴潮等的工程和预报的实际计算。三维模式是70年代开始发展的,近10年来发展很快。三维模式比较起二维模式来,其数值模拟自然要困难许多。这引起了研究者广泛的兴趣。现在三维模式的数值模拟,虽然提出的方法很多,也已有了不少成果,但研究工作仍是方兴未艾。      

Motion of two pulsating spheres in an ideal incompressible fluid

The problem of the interaction of two pulsating spheres in an ideal incompressible fluid was first investigated in detail by Bjerknes [1]. However, his and subsequent studies on this subject [2–5] were restricted to the interaction forces between the spheres, whereas the law of their motion was not considered because of the much greater complexity of the corresponding problem. The aim of the present paper is to find an approximate analytic solution to the problem of the motion of two pulsating spheres in an ideal incompressible fluid filling the entire space exterior to the spheres under the assumption that the flow of the fluid is irrotational.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 159–162, May–June, 1983.     

CFD-DEM simulation of the hole cleaning process in a deviated well drilling: The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid–particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle−particle, particle−wall, and particle−drill pipe are taken into account with the Hertz–Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the multi-sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid–particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid–solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.     

复合材料水翼水动力与结构强度特性数值研究

针对复合材料水翼存在的流固耦合求解问题,结合其自身特有属性,对复合材料水翼结构变形特性进行了数值仿真计算研究.研究建立了复合材料水翼流固耦合数值计算模型,并将数值计算结果与Zarruk等的实验结果进行对比,验证模型的正确性,得出复合材料水翼尖端扭转角随雷诺数的增加而增加的研究结论.基于数值计算模型,系统地研究了不同铺层角对复合材料水翼水动力特性及强度特性的影响,结果表明:不同铺层角复合材料水翼的尖端扭转角,随铺层角的增大而减小,而其尖端位移量随铺层角的增大先减小后增大.为了削弱工程常数的影响对复合材料水翼变形的影响,研究提出了无量纲扭转角和无量纲位移量,进一步探究复合材料水翼结构的弯扭耦合作用对其变形特性的影响.最后利用蔡$\!$-$\!$-$\!$吴失效准则进行复合材料水翼强度特性的判断和分析,结果表明:不同铺层角复合材料水翼的蔡$\!$-$\!$-$\!$吴系数,随铺层角的增大呈现先减小后增大的趋势,其中0$^\circ$铺层时的复合材料水翼蔡$\!$-$\!$-$\!$吴系数最小,50$^\circ$铺层时的复合材料水翼的蔡$\!$-$\!$-$\!$吴系数最大.      

Freely rising light solid spheres

This paper examines the behavior of spheres rising freely in a Newtonian fluid when the ratio between the density of the spheres and that of the surrounding fluid is about 0.02. High-speed imaging is used to reconstruct three-dimensional trajectories of the rising spheres. From the analysis of the trajectories the magnitudes of the drag and lift forces exerted by the surrounding fluid are deduced. It is argued that the two main contributions to the drag force are (i) a viscous drag that may be estimated from the standard drag curve by evaluating the Reynolds number using the actual value of the velocity, and (ii) an inertial drag that arises essentially by the same mechanisms that cause the lift-induced drag familiar from wing theory. Estimates of both contributions, the latter using visualizations of the wakes of the spheres, give a favorable agreement with the measured drag forces. These findings are closely related to recent numerical results of in the literature on the forces experienced by oblate ellipsoidal bubbles rising in quiescent water.     

电磁力诱导二次流对微流体混合效果影响的研究

设计了一种微流道中液体混合的方法, 基于电磁力的作用, 诱导微流道中二次流动 发生, 从而产生Dean涡, 使流道中流体出现往复运动及流体界面的弯曲延伸, 使不 同流体的接触面积大大地增加, 从而提高混合效率. 对3种结构的工况进行了比较. 建立了微流道混合系统的理论模型并进行了数值模拟,对流体混合过程进行了分析并对混合效率进行了评价. 进行了验证性的试验研究, 对模拟结果与试验结果进行了比较, 结果相吻合.     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

基于虚拟区域法的黏弹性流体中微生物游动的数值模拟和应用

微生物是自然生态系统的重要组成部分,掌握微生物在复杂流体中的运动特性可以为微型器件的设计制造提供理论指导.壁面效应是微生物游动研究中的重要问题之一,已有研究表明微生物在壁面附近存在复杂的行为特征.然而已有研究大多集中于微生物在牛顿流体中的游动模拟,仅有少数涉及黏弹性流体等非牛顿流体.本文采用直接力虚拟区域法与乔列斯基分解相结合的数值方法,引入Squirmer微生物游动模型,研究了微生物在黏弹性流体中的游动问题.首先给出求解黏弹性流体本构方程的数值格式;并将该方法应用于研究微生物游动中的壁面效应.研究结果表明,游动方向是影响微生物颗粒壁面效应的重要因素.流体弹性应力会对微生物产生一个反向转矩,影响微生物的游动方向,从而阻碍微生物逃离壁面.微生物颗粒在黏弹性流体中与壁面作用时间较长,几乎达到牛顿流体的两倍以上.     

混合式CRP面元法计算对比

目前,对于混合式CRP(contra-rotating propeller)的数值研究大部分是基于黏流方法,这主要是由于混合式CRP组成构件较多,运用势流方法求解存在一定的困难.但势流方法求解效率较高,在混合式CRP的设计方面,具有独特的优势.为建立混合式CRP的高效数值分析方法,从而为混合式CRP的设计工作奠定基础,研究首先将混合式CRP分成前桨和吊舱推进器两部分,以单个螺旋桨面元法和吊舱推进器整体面元法为基础,建立了混合式CRP的迭代面元法.随后对混合式CRP的几何结构以及桨叶面元奇点强度的变化特点进行分析,提出了将混合式CRP各构件进行整体计算的面元法,推导了整体面元法的计算公式.编译完成迭代面元法以及整体面元法的计算程序,并对混合式CRP的敞水性能进行了计算分析,两种面元法的计算结果与试验值的对比表明;迭代面元法与整体面元法计算得到的推力、扭矩系数相对误差均能控制在5%以内,但整体面元法的计算时间比迭代面元法节约1/3左右,同时整体面元法省去了迭代步骤,适宜于混合式CRP的设计.最后对研究中所建立的整体面元法的理论误差进行了分析.     

Squeeze flow of interstitial Herschel–Bulkley fluid between two rigid spheres

Interaction between two spheres with an interstitial fluid is crucial in discrete element modeling for simulating the behaviors of ‘wet’ particulate materials. The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial Herschel–Bulkley fluid was studied on the basis of Reynolds’ lubrication theory, resulting in analytical integral expressions of pressure distribution and the viscous force between the two spheres. According to the variation of shear stress, the fluid was divided into yielding and unyielding regions, followed by a discussion on the thickness of the two regions. The result of this paper could be reduced to either the power-law fluid or the Bingham fluid case.     

FFP在黏土中贯入过程的CFD模拟

圆锥静力触探法(Cone penetration test, CPT)广泛应用于土的抗剪强度测试中。在此基础上发展的自由落体式贯入仪(Free fall penetrometer, FFP),依靠自由下落获得的动能和自身重力势能贯入土中,不需要借助外部加载装置,提高其使用的便捷性。但在动力贯入过程中FFP与土的相互作用更为复杂,涉及到土体的率效应和拖曳阻力等的影响。因此,对FFP各项受力和相关参数的准确分析有助于提高其实用性以及测量的准确性。作者采用基于计算流体动力学(computational fluid dynamics, CFD)的分析软件ANSYS CFX 17.0模拟FFP在均质黏土中的贯入过程,借助动网格的大变形分析方法来模拟FFP贯入过程中的运动边界问题。提出了薄层单元法模拟FFP与土体的界面摩擦接触行为。在CFD模拟中,土体材料采用非牛顿流体来模拟,其剪切强度受土体切应变率的影响(即土体的率效应)。通过模拟贯入仪在黏土中以不同的速度贯入的过程,研究FFP的端部阻力和侧壁阻力与贯入速率、土体强度和密度、界面摩擦系数以及率效应参数之间的关系,建立了端部承载力系数、端部和侧壁率效应参数及拖曳系数的表达式,并提出了土体不排水抗剪强度的预测方法,为FFP测试数据的解析提供依据。      

THE TRANSVERSE RESPONSE OF SANDWICH PANELS TO AN UNDERWATER SHOCK WAVE

The response of submerged structures to an underwater shock wave involves both structural and fluid behaviour. For sandwich structures the response to an initial shock wave in the transverse direction is significantly different from that of a homogeneous structure. This is due to the elastic properties of the core in the sandwich. For a homogeneous structure one cavitation zone is initially developed and the position of this zone is dependent on the assumed cavitation pressure. At a sandwich structure two cavitation zones initially appear, one adjacent to the structure and another away from the structure, depending again on the assumed cavitation pressure. The response of a sandwich panel in the transverse direction is also investigated, using a combination of a numerical method and finite elements, developed for the fluid–structure interaction problem. The method includes the appearance of cavitation in the fluid, it is found that, as the sandwich section is moving, the faces in the sandwich oscillate about the core. The nonlinear properties of the faces give rise to considerable transverse strain which may very well be high enough to cause delamination in the faces, debonding between the faces and the core, or core failure.