均匀流中近壁面垂直流向振荡圆柱水动力特性研究

对均匀来流下靠近壁面处在垂直流向做强迫振荡运动的光滑圆柱的水动力特性进行了试验研究. 试验在拖曳水池中进行, 雷诺数为2× 10⁵, 通过采集顺流向和垂直流向的力, 得到了阻力系数、升力系数、相位角等与间隙比、振荡频率和振幅之间的关系. 通过研究得到如下结论: 1)振荡圆柱的平均阻力系数在近壁面处随间隙比的减小而骤降; 2)振荡圆柱泄涡受到完全抑制的临界间隙比要小于静止圆柱; 3)近壁面的存在对振荡圆柱的能量传递有着重要的影响, 自由边界圆柱强迫振荡所得到的水动力系数不能用来预报海底管道的涡激振动; 4)对于振荡圆柱, 附加质量系数只有在一定的频率范围内才是定值, 且在低频率区域其绝对值随间隙比减小而增大; 5)圆柱在进行强迫振荡时, 其平均阻力系数、振荡阻力系数和振荡升力系数均随无因次振幅的增加而增大. 关键词： 海底管道 强迫振荡 水动力特性 涡激振动     

Viscous Interaction Between a Vortex Tube and a Rotating Spherical Particle

The transient advection of a cylindrical vortex tube in a viscous incompressible flow field and its interaction with a rotating/spinning spherical particle has been investigated numerically at Reynolds numbers in the range of 20≤ Re≤200 for angular velocities of 0≤Ω≤0.5. The effects of vortex parameters such as size, circulation strength and initial position relative to the particle, on the temporal behavior of the lift and drag forces are studied. Vortex‐sphere interactions bring about major changes in the flow field particularly when coupled with particle rotation. It is observed that the forces acting on the particle are significantly influenced during the time that the vortex core is in the vicinity of the particle. The extent of these local changes are about ±30% in the drag coefficient and about ±200% in the lift coefficient as compared to flow over a rotating solid sphere with no vortex interaction. It is also found that a vortex with core radius between one and two particle diameters creates the strongest temporal variations in the lift and drag coefficients. Furthermore, maximum lift variations occur for the vortex‐particle head on collision, while a vortex with an offset distance of about one diameter from the principal flow axis generates the maximum drag variations.     

煤粉颗粒所受Magnus力的数值模拟

本文采用数值模拟方法计算了煤粉颗粒所受的Ｍａｇｎｕｓ力,考察了颗粒的旋转速度、流动Ｒｅ数对Ｍａｇｎｕｓ力的影响,并且给出了相应的关系公式．通过比较煤粉颗粒在运动过程中所受到的Ｍａｇｎｕｓ力与气动力的大小,得出：当煤粉颗粒的旋转速度为１８００转／分时,其大小约为气动力的１％,因此在实际的煤粉颗粒受力分析中是完全可以忽略该力的。     

电磁力控制湍流边界层分离圆柱绕流场特性数值分析

在亚临界区高雷诺数Re=1.4×105下,采用脱体涡模拟结合湍流分离的方法对弱电解质中电磁力作用下湍流边界层分离圆柱绕流场及其升(阻)力特性进行了数值模拟和分析.结果表明,电磁力可以提高圆柱体湍流边界层内的流体动能,延缓圆柱体湍流边界层的流动分离,减弱圆柱体湍流绕流场中在流向和展向上大尺度漩涡的强度,减小圆柱体阻力时均值及其升力脉动幅值.当电磁力作用参数大于某个临界值后,湍流边界层流动分离消失,在圆柱体尾部产生射流现象,从而电磁力对圆柱体产生净推力作用,出现负阻力现象,而且升力脉动幅值接近于零,出现圆柱体升力消失现象.     

疏水表面减阻的格子Boltzmann方法数值模拟

采用格子Boltzmann方法的多松弛模型和Shan-Chen多相流模型对雷诺数为100的疏水表面方柱绕流进行了数值模拟, 分析了疏水表面接触角和来流含气率对方柱绕流流场的影响. 研究结果表明: 疏水表面接触角一定时, 来流含气率在一定范围内, 疏水表面具有减阻的能力, 超出这一范围时会出现阻力系数、升力系数升高的现象, 同时在方柱近壁面处伴随涡的形成产生了气团脱落; 当来流含气率处于适当水平时, 接触角越大, 绕流物体近壁面处含气率越稳定, 减阻效果越明显. 分析发现疏水表面减阻的关键在于保证近壁面处气层的稳定性, 此时接触角越大, 减阻效果越明显. 本文从含气率角度出发分析疏水表面的减阻现象, 为进一步探索疏水表面减阻机理提出了新的思路.     

Drag and Lift Forces Acting on a Sphere in Shear Flow of Power-Law Fluid

Laminar flow of a power-law fluid over a sphere is considered for unbounded shear flow. The Navier–Stokes equations with power-law viscosity are solved numerically using an in-house developed CFD package. Vorticities structures downstream of particle are suppressed for powerlaw fluid. The shear rate influence on drag force is negligible for power index close to unit, and the drag force appreciably decreases with falling power index. For small Reynolds numbers, the lift force coefficient monotonically decreases against the power index and exhibits an opposite behavior for moderate values of Reynolds numbers. The results of the parametric studies are used to derive correlations for the drag force and to detect the hydrodynamic differences from uniform flow. The investigation parameters varied within the following ranges: power-law index 0.3 ≤ n ≤ 1, Reynolds number 0 < Re ≤ 150, and dimensionless shear rate 0.05 ≤ s ≤ 0.4.     

Generation of vortices by a streamwise oscillating cylinder

The wake of a streamwise oscillating cylinder is presently investigated. The Reynolds number investigated is 300, based on the cylinder diameterd. The cylinder oscillates at an amplitude of 0.5d and a frequencyf _e/f_s=1.8, wheref _e is the cylinder oscillating frequency andf _s is the natural vortex shedding frequency of a stationary cylinder. Under these conditions the flow is essentially two dimensional. A two-dimensional direct numerical simulation (DNS) scheme has been developed to calculate the flow. The DNS results display a street of binary vortices, each containing two counter-rotating vortical structures, symmetrical about the centerline, which is in excellent agreement with measurements. The drag and lift on the cylinder have been examined. The time averaged drag and lift are 1.4 and 0, respectively, which are the same as those on a stationary cylinder at the sameRe. However, the fluctuating drag was high, about 2.68. It has been found that, being symmetrically formed about the centerline, the binary vortices induce an essentially zero fluctuating lift, which may have a profound implication in flow control and engineering.     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.     

Mechanism of controlling turbulent channel flow with the effect of spanwise Lorentz force distribution

A direct numerical simulation(DNS) is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction(or even drag increase) is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.     

滴形颗粒在通道中沉降的格子Boltzmann-虚拟区域方法模拟

采用格子Boltzmann-虚拟区域方法对滴形颗粒在垂直通道中的沉降过程进行直接数值模拟.通过格子Boltzmann方法求解N-S方程,流体与固体之间的相互作用通过虚拟区域方法描述.研究雷诺数在10^-2到100范围内颗粒形状因子对其摩擦系数和阻力系数的影响.为便于比较,给出了圆形颗粒沉降的结果.结果发现,当雷诺数小于1的时候,颗粒的摩擦系数始终保持常数,而当雷诺数大于1时摩擦系数随雷诺数的增大而增大.此外,当雷诺数小于约30时圆形颗粒的摩擦系数和阻力系数均小于滴形颗粒,而当雷诺数大于30时情况正好相反.颗粒周围的压力分布证明了这一结论.     

Numerical simulation of flow around a circular cylinder with curved sectional grooves

In a circular cylinder with uniform flow, a sudden decrease in the drag force occurs at a high Reynolds numbers; however, it is known that the same phenomenon occurs at a lower Reynolds number in the case where there exist grooves or roughness on the circular cylinder surface. To clarify the flow characteristics around a circular cylinder in the case of changing the shape of grooves we analyzed the drag coefficient, lift coefficient, turbulent kinetic energy, vorticity and pressure by applying the RNGk-ε turbulent model. The shapes of the grooves were arced, triangulated and curved. The results showed that the separation point for a circular cylinder with curved sectional grooves shifts to the most downstream side and the drag coefficient becomes the smallest among circular cylinders with grooves.     

EXPERIMENTAL INVESTIGATION ON THE HEAT TRANSFER OF TURBULENT FLOW IN A PIPE OSCILLATING AT LOW FREQUENCIES

An experimental study was performed focusing on heat transfer and friction coefficient associated with turbulent oscillating tube flow. For this goal an oscillating mechanism was designed. Experiments were conducted for the low oscillating frequency in the range of 0.008–1.988 Hz and dimensionless amplitude was chosen as X₀ = 0.3, 0.6, and 0.9. Reynolds number was changed from 0.5 × 10⁴ to 2.5 × 10⁴. The bulk temperature of the fluid at the exit of the oscillating section was fond to be increasing with oscillating frequency and amplitude. For the oscillating cases, heat transfer enhancement is obtained 52% for f = 1.988 s^?1, 40% for f = 1.320 s^?1, and 28% for f = 0.008 s^?1, in comparison with the smooth pipe at the highest Reynolds number. The results also showed that Nusselt number and friction coefficient also increased with increasing frequency and amplitude.     

Forces acting on a stationary sphere in power-law fluid flow near the wall

The analysis and evaluation of the forces acting on the particle in a linear shear flow of power-law fluid (PLF) in the presence of the wall were performed. Using the results of a series of computations for a model problem with a spherical particle near a flat wall in the Reynolds number range of 0?200 and the distance to the wall from 0 to 20 particle diameters, the correlation formulas for calculating the coefficients of drag force and lift force were obtained. Special attention was paid to the behavior of the forces acting on the particle approaching the wall.     

通道中串列旋转圆柱的格子Boltzmann-虚拟区域方法的模拟

通过结合格子Boltzmann方法(LBM)和虚拟区域(Fictitiou sDomain)思想,建立格子Boltzmann-虚拟区域(LB-DF/FD)方法.采用两套网格系统,欧拉网格用于流体,拉格朗日网格用于固体.原有的LBM在计算运动固体的受力方面存在数据振荡,LB-DF/FD方法改进了此缺陷.为验证该方法,模拟圆柱绕流、圆形颗粒在无限长通道中平动及在无限大流场中转动三种情况,结果与其他数值解及理论解符合得很好.利用该方法模拟低雷诺数下通道中串列旋转圆柱周围的流场,分析圆柱间距(g)及雷诺数(Re)对流场结构的影响.给出Re=0.001,0.1和10下,0.2≤g≤8.0的流线结构、圆柱升力、阻力以及力矩等数值结果.结果表明,g对流场的结构及圆柱的受力有显著影响,Re对圆柱阻力及Stokes单元数目的影响较大.     

Frenkel-Kontorova模型中基底势振动的影响

基于一维Frenkel-Kontorova模型, 研究了振动的基底势对系统纳米摩擦现象的影响. 分别在相邻原子间的距离与周期势场的周期比为不公度(incommensurate)、可公度(commensurate)两种情形下, 探讨了基底势振动的振幅和频率对滞回现象(hysteresis)、最大静摩擦力以及超滑现象的作用机理. 两种情形下, 固定频率, 随着振幅的增大, 滞回区域的面积以及最大静摩擦力都将减小, 对于不同的频率, 减小的趋势不同. 系统甚至产生了超滑现象. 但当频率过大时, 振幅的改变不会影响滞回区域的面积以及最大静摩擦力的大小, 此时与基底不加振动时的情形一致; 当振幅固定, 随着频率的增大, 滞回区域的面积将增大, 对于不同振幅, 增大的趋势不同. 特别地, 对于某些固定的振幅, 最大静摩擦力随着振动频率的增大先逐步减小直至出现超滑现象, 再进一步增大频率, 最大静摩擦力又转而逐步增大. 这一现象类似于共振, 表明存在最佳的振动频率促进系统内所有原子的共同运动, 使得整个系统的最大静摩擦力几乎消失. 另外, 两种情形的区别是, 对于某些固定的频率(如ω= 0.5)和不同的小振幅, 不可公度情形往往具有相同的平均终止速度, 而可公度情形则不同, 表明相同前提下后者具有更复杂的动力学行为. 关键词： Frenkel-Kontorova模型 滞回 最大静摩擦力 超润滑     

振荡射流改善翼型气动性能的实验研究

本文针对采用振荡射流控制流动分离改善大攻角下翼型气动特性的问题,在NACA0015翼型上进行了多种工况的风洞实验。结果表明:在失速攻角附近,振荡射流抑制流动分离提高翼型升力系数的作用十分明显,可将翼型失速攻角推迟2°左右。存在最佳的振荡射流频率段、射流动量范围和射流位置,使得翼型性能的改善最大。实验还得到了振荡射流的频率、动量和施加位置等参数对翼型气动性能的影响规律。     

具有多孔介质覆面层的柱体绕流流场分析

建立多孔结构覆面柱体绕流模型, 采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对覆盖多孔介质层的方柱绕流进行数值模拟, 研究多孔介质对钝体绕流流场特性的影响。结果表明: 相比于不可渗透壁的柱体, 引入合适参数的多孔介质覆面层后可以有效降低其升力脉动幅值, 但阻力有所增加。同时, 较高雷诺数下多孔方柱的数值模拟表明: 多孔介质壁面使得尾迹区域的剪切层相距更远, 降低了尾流处湍动能, 并将雷诺应力的峰值移动到尾迹区域, 抑制了方柱两侧的动量交换, 使动量交换的位置发生在尾迹区域, 继而使得尾迹的涡街更加规则化。     

Flow pattern and lift evolution of hydrofoil with control of electro-magnetic forces

The initial responses and evolutions of the flow pattern and lift coefficient of a hydrofoil under the action of electro-magnetic (Lorentz) force have been studied experimentally and numerically, and trace particle methods are employed for them. With the introduction of BVF (boundary vortex flux), the quantitative relation among Lorentz forces, BVF and lifts is deduced. The influences of flow patterns on the hydrofoil lift coefficient have been discussed based on the BVF distribution, and the flow control mechanism of Lorentz force for a hydrofoil has been elucidated. Our results show that the flow pattern and lift of the hydrofoil vary periodically without any force. However, with the action of streamwise Lorentz forces, the separation point on the hydrofoil surface moves backward with a certain velocity, which makes the flow field steady finally. The streamwise Lorentz force raises the foil lift due to the increase of BVF intensity. On the other hand, Lorentz force also increases the hydrofoil surface pressure, which makes the lift decrease. However, the factor leading to the lift enhancement is determinant, therefore, the Lorentz force on the suction side can increase the lift, and the stronger the Lorentz force, the larger the lift enhancement. Our results also show that the localized Lorentz force can also both suppress the flow separation and increase the hydrofoil lift coefficient, furthermore, the Lorentz force located on the tail acts better than that located on the front.     

翼型绕流的洛伦兹力控制机理

在翼型上翼面壁面附近流场中形成的流向洛伦兹力,可提升翼型的升力减小阻力,然而制约其推广应用的主要瓶颈是极为低下的控制效率,为提高洛伦兹力的控制效率,需研究其控制机理.以翼型绕流的洛伦兹力控制为例,利用双时间步Roe格式及水槽对其进行数值及实验研究.结果表明:洛伦兹力的控制效果随着来流速度的增加而下降,升力增幅和阻力减幅与来流速度大小呈反比关系,但升力增加和阻力减小的规律不变,都是升力先急剧增加随后缓慢增加,而阻力先急剧减小然后再缓慢增加,基本原因为升力和阻力先受洛伦兹力推力的影响而分别增加和减小,随后洛伦兹力作用增加翼面壁面摩擦力,导致升力减小和阻力增加,流向洛伦兹力还导致翼型壁面压力下降,增加翼型升力和压差阻力;壁面摩擦力导致的升力降幅比壁面压力变化导致的升力增幅小,壁面压力变化起主导作用;洛伦兹力推力对阻力的降幅比压差阻力的增幅大,洛伦兹力推力起主导作用,因此阻力减小.     

纵向受迫振荡圆柱绕流问题的数值模拟

用有限体积法对平行于均匀来流方向受迫振荡的圆柱绕流问题进行了二维数值模拟.雷诺数选取Re=200、855、4000等几种亚临界雷诺数情况.通过研究不同振幅和振动频率下的流场结构和一些重要流动参数如升阻力系数、Strouhal数等随Re数、KC数、Stokes数的变化关系,验证了实验中观察到的一定条件下发生的"频率锁定"现象,并将涡脱落方式划分为三种主要模态.文中引入网格速度,对常用的处理速度与压力耦合的SIMPLE算法作了适当的补充和修改,以适应随时间变化的网格坐标.