明渠层流近床面颗粒拖曳力研究

泥沙颗粒受到的拖曳力是泥沙运动的主要驱动力，而当前应用于计算流体力学-离散颗粒法(CFD-DPM)耦合模型进行水沙运动模拟的泥沙颗粒拖曳力公式均没有考虑明渠流底床边壁作用的影响。求解不可压缩Navier-Stokes方程，对明渠层流不同雷诺数条件下床面附近不同高度处颗粒所受拖曳力进行了模拟，根据模拟结果变化规律，提出了综合考虑床面和水流惯性对标准拖曳力影响的修正拖曳力计算公式。与常用的单颗粒标准拖曳力公式和考虑遮蔽效应的多颗粒拖曳力公式相比，采用本文修正公式得到的水沙作用力更接近高精度数值解，应用于CFD-DPM输沙模拟获得的输沙结果与输沙率公式结果一致，应用分析表明输沙模拟应当采用粗糙底床边界。     

近床面水平圆柱局部冲刷二维数值模拟

为研究近床面水平圆柱的局部冲刷问题,基于N-S方程和有限体积法,在FLUENT中通过二次开发建立了局部冲刷二维数值模型。模型采用标准κ-ε紊流模型来计算水平圆柱周围的流场,同时借助FLUENT软件中的自定义函数功能提取床面剪应力参数来计算该时刻的推移质输沙率及床面节点位移变化值,然后运用动网格技术来模拟床面地形的变化,通过物理模型试验来进行模型验证。结果表明,计算结果与试验结果基本一致,从而证实了冲刷模型的准确性。     

透水床面明渠湍流时空平均特征的大涡模拟与分析

基于大涡模拟数据,研究了理想粗糙透水床面明渠湍流的时空平均特性. 考虑到空间异构性,对比分析了不同位置的时空平均流速、雷诺剪应力、构造剪应力、脉动幅度的垂线分布. 结果表明:第一,顶层床面之上,空间异构性的影响较小,不同位置的双平均流速符合类似的对数分布,但由于透水床面影响,卡门常数较不透水床面小;在床面附近,空间异构性影响较大,不同位置的双平均流速分别符合线性分布与多项式分布;在透水河床内部,靠近底层球孔的双平均流速为上部球孔双平均流速的1.55 倍. 第二,床面之上,雷诺剪应力占总剪应力的95% 以上,占有主体地位;床面附近,紊动较大,构造剪应力不能忽略,其值大约占总剪应力的15%.由于流场的各向异性,纵向与垂向的脉动幅度有所差异.      

压敏漆测压在低速射流撞击壁面的应用研究

压敏漆(pressure sensitive paint,PSP)测压是一种基于发光分子氧猝灭效应的测压技术,与传统测压方式相比具有成本小、空间分辨率高等特点.实验搭建了压敏漆测压硬件平台,开发了数据处理软件系统.利用驻点处的压强与灰度关系拟合了标定曲线,研究了低速射流撞击平板表面的压强分布特性,分析了喷管距离壁面高度和雷诺数的影响.结果表明,当喷管高度h/d<6时,压强分布曲线的峰值压强和半宽度随喷管高度变化基本不变;当喷管高度h/d>6时,随着高度的增大,峰值压强线性减小,半宽度线性增加.采用峰值压强和半宽度对不同喷管高度的压强分布曲线进行无量纲化,压强分布曲线呈现自相似特性.进一步研究表明,在实验雷诺数范围内,压强分布受雷诺数影响较小.      

床面上直立圆柱的三维湍流数值模拟

从数值预报桥墩等结构物床面局部冲刷的角度发展绕直立圆柱的三维湍流的数值模拟技术. 基于Wilcox的k-ω两方程湍流模式，采用基于有限体积法的压力修正SIMPLE算法, 计算了绕床面直立圆柱的三维湍流流场，分析了光滑和粗糙床面两种情况下的流动情况. 通过系列的验证计算，表明该计算模型能够比较准确地反映不同外来流条件下绕直立圆柱的流场. 计算结果揭示了床面粗糙度对绕圆柱的湍流流动的影响.     

沙粒跃移运动的数值模拟

沙粒在风场作用下形成的风沙运动,直接导致一系列环境问题.现有的数值模拟并没有揭示风沙流的二维性质.本文提出一种合理有效的模拟方法来揭示风沙流时空的变化.气相运动采用大涡模型,而固相(沙粒)处理为高散体系(DEM),并利用已有的沙粒与床面碰撞的激溅函数,以获得整个风沙系统运动特性.结果表明:其发展过程与以往的试验结果相比较符合:风沙流达到稳定后,不同截面输沙率沿高度分布与平均输沙率相比差别很大,并大约在9cm高处,平均输沙率出现极值,单宽输沙率沿风向分布同样具有非均匀性.     

明渠湍流边界层中颗粒的运动与分布

利用直接数值模拟、点球浸入边界法和颗粒离散元法相结合的方法, 模拟了颗粒在明渠湍流边界层中的运动, 并对颗粒的瞬时位置进行了Voronoi 分析, 定量研究了颗粒在湍流边界层中的运动和分布规律. 研究发现:颗粒的输运对湍流的统计特征有影响, 其运动与近壁区湍流拟序结构密切相关, 在"喷发"结构作用下被带离壁面, 在"扫掠" 结构和自身重力作用下回到壁面; 在湍流边界层中, 颗粒倾向于聚集在低流速带, 呈条带状分布;颗粒在大部分时间处于"簇"状态, 偶尔跳跃到"空" 状态, 但能够很快回到邻近低速区域.      

不同高度柔性植被对坡面流水动力特性的影响

基于系列坡面流试验,研究了柔性植被的高度(4 cm,8 cm和15 cm)对坡面流水动力特性的影响.研究表明:覆盖率相同时,4 cm植被阻水效果最好,15cm最差,而该差异随着坡度增大而减小.4 cm和8 cm植被覆盖时,水流阻力系数随雷诺数增加先减小后增加,这是由于雷诺数较小时水流阻力主要以床面阻力为主,雷诺增大时主要以植被阻力为主.另外,通过染色剂法测得光滑床面平均流速修正系数为0.2～0.4,有植被覆盖床面为0.4～0.8,且修正系数随平均流速增大而增大.     

基于CFD流态分析技术的座便器虹吸管道优化设计

在Fluent软件中应用RNG $k$-$\varepsilon$湍流模型及流体体积函数(VOF)对座便 器内部三维湍流流动进行了雷诺平均N-S方程的数值模拟，得到 了座便器内流场三维流动形态，研究了座便器虹吸管内流动规律，分析了虹吸管形状对流速 分布、压力分布及虹吸性能的影响，以三维湍流场的分析结果为依据，实现了座便器虹吸管 道的优化设计. 通过PIV测试数据验证了三维湍流数值模拟结果的准确性.     

稳定分层流密度界面处湍流混合与分形结构

采用室内混合箱研究稳定分层流(上层淡水、下层盐水) 无剪切密度界面处的湍流混合与分形结构. 湍流通过浸没在盐水层中的振动格栅产生, 密度界面结构通过在盐水层中添加荧光剂或染料可视化, 共进行12 组实验. 实验观测并记录了:(1) 淡、盐水密度界面距混合箱底部的平均高程(h);(2) 淡、盐水层的密度(ρ₀,ρ), (3) 淡、盐水密度界面. 其中, 淡、盐水密度界面通过照片、录像进行记录. 观测结果用于计算:(1) 盐水层密度;(2) 卷挟速度, (3) 整体理查孙数(Ri_o), (4) 二维、三维密度界面, (5) 二维、三维密度界面的分形维度. 结果分析发现:(1) 湍流卷挟率随Rio 增大而减小, 并且满足Ri_o的-3=2 或-7=4 幂律;表明随着湍流强度的减弱, 混合的速度也越来越缓慢;(2) 二维密度界面分形维度大于1, 三维密度界面分形维度大于2;表明二维、三维密度界面存在分形结构;(3) 分形维度随Ri_o的增大而减小;表明随着湍流强度的减弱, 密度界面也越来越趋于光滑.      

Natural Low-Frequency Riblets

A turbulent flow at modest Reynolds numbers (26,000–32,000, as defined in the text) over flat erodible sandy beds was investigated. The beds were composed of sand grains of uniform size with a diameter smaller than the viscous sublayer of the flow. When a modest sand transport (1.2 × 10-3–8.5 × 10-3 g/s·cm) was present, the near-wall coherent structures of the flow produced a bed form with a streaky longitudinal pattern. This pattern had a lateral spacing of about 100 viscous units and had less drag than a flat wall.The mechanism causing the drag reduction is due to a low-frequency riblet pattern which allows the flow to sustain higher wall shear stresses without additional erosion. That means during the riblet regime the sediment transport remains quasi-constant.The formation of the sand streaks as well as the drag reducing mechanism are discussed.     

基于滑移理论的超疏水表面减阻性能分析

采用滑移理论计算和数值模拟方法研究超疏水表面的减阻性能.层流流动状态的数值模拟和理论计算结果的一致性较好,并可证明滑移速度与通道内流体速度之比与无量纲压降比的数量级相当;湍流状态的数值仿真结果表明,在目前可实现的滑移速度范围内,超疏水表面对水下航行器的流体阻力影响不显著.     

The loss of stability of laminar flow in open channel and the mechanism of sand ripple formation

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｒｉｎｇｔｈｅｔｉｍｅｏｆｂｅｄｓｕｒｆａｃｅｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔａｔｉｏｎ,ｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｆｂｅｄｌｏａｄｍａｋｅｔｈｅｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｖａｒｉｏｕｓｄｉｆｆｅｒｅｎｔｆｏｒｍｓｏｎｔｈｅｒｉｖｅｒｂｅｄｓｕｒｆａｃｅ,ｔｈｕｓｔｈｅｆｏｒｍｏｆｂｅｄｓｕｒｆａｃｅｃｈａｎｇｅｓｃｏｎｓｔａｎｔｌｙ .Ｔｈｉｓｋｉｎｄｏｆｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｎｔｈｅｂｅｄｓｕｒｆ…     

粘土悬浮液明渠流的减阻现象

水和粘土悬浮液在水槽中的流动试验证明，在粗糙边界上的粘土悬浮液明渠流阻力明显小于清水流．在同样水力条件下粘土悬液的时均流速比清水流速大得多．当悬液浓度较高时，流动阻力系数比清水流阻力系数的一半还小．实验表明粘土悬液的高粘性和屈服应力抑制了紊动的发展，减小了紊动剪力，导致阻力降低．对于光滑边界湍流，粘土悬液的阻力系数与清水流的相近，不发生减阻现象     

Numerical and experimental investigation of turbulent characteristics in a drag-reducing flow with surfactant additives

Cetyltrimethyl ammonium chloride (CTAC) surfactant additives, because of their long-life characteristics, can be used as promising drag-reducers in district heating and cooling systems. In the present study we performed both numerical and experimental tests for a 75 ppm CTAC surfactant drag-reducing channel flow. A two-component PIV system was used to measure the instantaneous streamwise and wall-normal velocity components. A Giesekus constitutive equation was adopted to model the extra stress due to the surfactant additives, with the constitutive parameters being determined by well-fitting apparent shear viscosities, as measured by an Advanced Rheometric Expansion System (ARES) rheometer. In the numerical study, we connected the realistic rheological properties with the drag-reduction rate. This is different from previous numerical studies in which the model parameters were set artificially. By performing consistent comparisons between numerical and experimental results, we have obtained an insight into the mechanism of the additive-induced drag-reduction phenomena.

Our simulation showed that the addition of surfactant additives introduces several changes in turbulent flow characteristics: (1) In the viscous sublayer, the mean velocity gradient becomes gentler due to the viscoelastic forces introduced by the additives. The buffer layer becomes expanded and the slope of the velocity profile in the logarithmic layer increases. (2) The locations where the streamwise velocity fluctuation and Reynolds shear stress attain their maximum value shifted from the wall region to the bulk flow region. (3) The root-mean-square velocity fluctuations in the wall-normal direction decrease for the drag-reducing flow. (4) The Reynolds shear stress decreases dramatically and the deficit of the Reynolds shear stress is mainly compensated by the viscoelastic shear stress. (5) The turbulent production becomes much smaller and its peak-value position moves toward the bulk flow region. All of these findings agree qualitatively with experimental measurements.

Regarding flow visualization, the violent streamwise vortices in the near wall region become dramatically suppressed, indicating that the additives weaken the ejection and sweeping motion, and thereby inhibit the generation of turbulence. The reduction in turbulence is accomplished by additive-introduced viscoelastic stress. Surfactant additives have dual effects on frictional drag: (1) introduce viscoelastic shear stress, which increases frictional drag; and (2) dampen the turbulent vortical structures, decrease the turbulent shear stress, and then decrease the frictional drag. Since the second effect is greater than the first one, drag-reduction occurs.     

圆柱绕流局部冲刷机制的实验研究

采用粒子图像测速技术和粒子图像识别技术研究了冲刷坑的演化过程及局部冲刷的流体动力机制。首先研究了动床条件下圆柱周围的床面演化过程，给出垂直对称面输沙率随时间变化以及圆柱下游垂直对称面的沙颗粒浓度和泥沙分布规律；其次针对典型冲刷坑形状固定底床，测量了圆柱下游的瞬时流场，并给出平均流场。结果表明尾流区平均剪应力较低，但垂向涡量和湍流强度较高。分析指出平均剪应力较低的尾流区域，局部冲刷的主要动力机制是尾涡和湍流的综合作用。     

Influence of curvature on drag reduction by opposition control in turbulent flow along a thin cylinder

Opposition controlled fully developed turbulent flow along a thin cylinder is analyzed by means of direct numerical simulations. The influence of cylinder curvature on the skin-friction drag reduction effect by the classical opposition control (i.e., the radial velocity control) is investigated. The curvature of the cylinder affects the uncontrolled flow statistics; for instance, skin-friction coefficient increases while Reynolds shear stress (RSS) and turbulent intensity decrease. However, the control effect in the case of a small curvature is similar to that in channel flow. When the curvature is large, the maximum drag reduction rate decreased. However, the optimal location of the detection plane is the same as that in a flat plate. Further, the drag reduction effect is achieved even on a high detection plane where the drag increases in the flat plate. Although a difference in the drag reduction effect can be observed with a change in the curvature, its mechanism considered in this analysis based on the transport of the Reynolds stress is similar to that of the flat plate.     

Turbulent mass transport and attenuation in Stokes waves

The motion of turbulent Stokes waves on a finite constant depth fluid with a rough bed is considered. First and second order turbulent boundary layer equations are solved numerically for a range of roughness parameters, and from the solutions are calculated the mass transport velocity profiles and attenuation coefficients. A new mechanism of turbulent mass transport is found which predicts a reduction and reversal of drift velocity in shallow water in agreement with experimental observations under turbulent conditions. This transpires because the second order Stokes wave motion, in a turbulent boundary layer, can directly influence the mass transport velocity by mode coupling interactions between different second order Fourier modes of oscillation. It is also found that the Euler contribution due to the radiation stress of the first order motion is reduced to half of it's corresponding laminar value as a consequence of the velocity squared stress law. The attenuation is found to be of inverse algebraic type with the reciprocal wave height varying linearly with either distance or time. The severe wave height restriction applicable to the Longuet-Higgins [4] solution is shown not to apply to progressive waves on a finite constant depth of fluid. The existence of sand bars on sloping beaches exposed to turbulent waves is predicted.     

NUMERICAL STUDY ON THE GROWTH PROCESS OF SECONDARY AEROSOL IN THE FOG

雾环境气溶胶之间的碰撞和凝并是气溶胶迁移、生长的动力学基础. 基于颗粒群平衡方程和多重蒙特卡洛方法,针对雾环境中二次液态气溶胶生长过程开展数值研究,着重分析了湍流和布朗作用机制对单分散和指数分布的二次气溶胶碰撞凝并的影响. 结果表明:雾环境气溶胶之间的碰撞凝并使得颗粒的总数目减少,颗粒平均体积逐渐增大. 对于初始尺度为自由分子区、单分散的气溶胶,布朗运动时间为600 s时,二次气溶胶平均尺度增大至初始的202.7倍左右,数目降至初始的0.006. 对于湍流作用下有布朗运动的二次气溶胶,在较短的时间内(如100 s)气溶胶尺度增至初始的163倍,颗粒数目降至初始的0.025倍, 说明雾环境流场的湍流运动引起二次气溶胶较强的输运和聚集,导致颗粒碰撞凝并概率增加,颗粒尺度增大.