减阻用表面活性剂溶液分子动力学模拟研究进展

减阻用表面活性剂在能源动力及化工领域有着广泛应用,在管道流体中加入少量表面活性剂可以使流动阻力大大降低从而节约能源,对于表面活性剂减阻机理的讨论也是近些年学者关注的热点之一.本文不仅对课题组前些年在表面活性剂溶液流变性、湍流减阻、减阻与传热的相关性、布朗动力学模拟方面的工作进行了概述,而且详细介绍了近三年来在表面活性剂粗粒化分子动力学模拟方面的研究成果.粗粒化模拟是近年来发展起来的方法,目前已广泛应用于化学、生物等诸多领域.在粗粒化分子动力学模拟方面的工作包括:表面活性剂溶液的流变性能与微观结构、表面活性剂溶液湍流减阻机理研究、湍流减阻失效分析三个部分.通过对表面活性剂溶液分子动力学模拟研究进展的回顾,作者认为,利用粗粒化分子动力学模拟方法可以合理揭示表面活性剂胶束的结构与流变性的对应关系,对胶束的断裂与再连接能力进行多维度的评价,如胶束的拉伸能、断裂能、最大拉伸长度、结合能、ζ电势、疏水基驱动作用等方面.并对"黏弹说"减阻机理进行分子模拟层面的验证,对实际应用中的湍流减阻失效原理进行初步分析.最后,根据对近几年分子动力学模拟工作的总结,展望了未来粗粒化分子动力学模拟在表面活性剂方面的研究方向.     

添加剂湍流减阻流动与换热研究综述

添加剂湍流减阻是指在液体的管道湍流中添加少量的高分子聚合物或某种表面活性剂从而使湍流阻力大大降低的现象.从其被发现至今,经过近半个世纪的研究(实验研究、理论分析、数值模拟和实际系统的应用研究),尽管对这一现象及其实际应用价值已有了较为深入的认识,但仍有许多方面尚有欠缺,例如对湍流减阻的机理仍然在探索中.本文归纳评述了高分子聚合物或表面活性剂添加剂湍流减阻流动与换热现象的研究现状,从湍流减阻剂的特性、减阻剂的湍流减阻机理、湍流减阻发生时的换热机理、减阻流动速度场分布和换热控制等几个方面综述了添加剂湍流减阻流动与换热特性,并综述了湍流减阻剂在实际工业系统中的应用情况,在对添加剂湍流减阻机理、有湍流减阻发生时的对流换热机理等的理解方面进行了新的总结.      

表面活性剂减阻溶液湍流流动研究进展

与聚合物添加剂相比,表面活性剂具有寿命长,不受机械力和高温影响发生降解的特点,目前被认为是最具有适用价值的减阻添加剂．针对表面活性剂湍流减阻机理,尽管已经开展了大量的实验研究、理论分析和数值模拟,但仍处于探索阶段,尚未定论．文中对目前国内外有关表面活性剂减阻溶液湍流流动的研究如平均速度、湍流强度、雷诺应力和相关系数等湍流统计量以及取得的成果进行了归纳和总结．目前普遍认为阻力减小是由流动垂直方向的湍流强度受抑和脉动速度分量解耦导致雷诺应力极大降低引起的．分析了目前研究存在的问题,对今后继续开展研究提出了自己的观点．      

水下减阻技术研究综述

在简要回顾早期减阻研究的基础上,对现有典型水下湍流减阻技术进行了较深入地分析.重点介绍了脊状表面减阻、微气泡减阻和疏水/超疏水表面减阻的研究现状.分别从实验研究和理论研究两方面对其进行了阐述,并着重强调了各自的减阻机理. 此外,还简要介绍了柔顺壁面减阻、壁面振动减阻等其它减阻技术.展望了水下减阻技术今后的研究重点及其应用前景.      

凹坑形表面在空气介质中的减阻性能研究

采用标准k-ε湍流模型对一种新型的非光滑表面———凹坑形表面在空气介质中不同条件下的流动进行了数值模拟,可以得到光滑表面和半径0.4mm的半球状凹坑分布的非光滑表面在4~48m/s范围内的来流速度下的阻力系数.计算结果表明该凹坑形非光滑表面在这个速度区间内均能产生一定的减阻效果.在速度为24m/s时,它的减阻率达到最大的7.2%.最后本文对该非光滑表面的减阻机理进行了研究.     

矩形槽道内表面活性剂减阻流体流场特性

采用粒子图像测速仪对矩形槽道内表面活性减阻流体在流动方向(x方向)与壁面垂直方向(y方向)所在平面的流场进行了测量,分析了速度、涡量、速度脉动相关量在流场内的瞬态分布,以及对500幅相同工况的流场进行了统计平均. 结果显示: 与牛顿流体相比, 表面活性剂减阻流体接近于层流流动,横向速度脉动被大幅减弱,导致湍流输运减弱,雷诺应力远远小于水. 减阻流体流向速度脉动呈条带特征,沿流动方向发展,反映了减阻流体不同于水的湍流输运特征.      

脊状表面减阻机理研究

针对脊状表面流场的特点,通过实验测量和数值模拟的方法对脊状表面微观流场进行了深入研究,获得了脊状表面湍流边界层的时均速度分布曲线、湍流度分布曲线和微观流场结构.为了得到脊状结构对壁面物性的影响,对脊状表面进行了疏水性测试,获得了液滴在脊状表面上的表观接触角,并通过水洞试验验证了脊状表面的减阻效果.研究表明,与光滑表面相比,脊状表面微观流场结构中存在"二次涡",近壁区的黏性底层厚度比平板的要厚得多,湍流度显著降低,且脊状表面表现出明显的疏水性.由此提出了基于壁面隔离效应、增大湍流阻尼效应和改变壁面物性效应的减阻机理.     

二维平板湍流边界层减阻研究

利用一套以位移传感器为主的弹性平衡装置,测量几种不同表面形状设计的平板在充分发展的二维湍流边界层中的减阻效果,并对其减阻机理进行研究。结果表明,大涡破碎器 LEBU(Large Eddy Break-up)和其它减阻装置的形状和布置对表面摩擦阻力有较大影响,在有些设计状态下的平板得到了净减阻。     

矩形槽道内表面活性剂减阻流体场特性

采用粒子图像测速仪对矩形槽道内表面活性减阻流体在流动方向(x方向)与壁 面垂直方向(y方向)所在平面的流场进行了测量,分析了速度、涡量、速度脉 动相关量在流场内的瞬态分布,以及对500幅相同工况的流场进行了统计平均. 结 果显示: 与牛顿流体相比, 表面活性剂减阻流体接近于层流流动,横向速度脉动被大幅 减弱,导致湍流输运减弱,雷诺应力远远小于水. 减阻流体流向速度脉动呈条带 特征,沿流动方向发展,反映了减阻流体不同于水的湍流输运特征.     

边界层对壁面湍流流动高分子减阻的影响

以“粘性”机制为理论基础,近年来在壁面湍流高分子减阻研究中提出了一种拉伸的高分子会产生自洽的等效粘度模型,这种等效粘度随离开壁面的距离而改变.通过等效粘度模型与Navier-Stokes方程的结合,运用雷诺应力模型计算壁面湍流减阻,并与基于高分子有限拉伸的非线性弹性哑铃模型的直接数值模拟结果进行比较,进一步校验了此等效粘度理论.通过肋条破坏槽道流中的边界层,显示了边界层对高分子减阻的影响,结果表明只有形成稳定的边界层,高分子才能有减阻作用.边界层是高分子减阻的首要条件,边界层中的粘性底层和对数率分布区之间的缓冲层可能是减阻的主要影响区域.     

Turbulent drag reduction by spanwise wall oscillations

In the present work a technique is numerically investigated, which is aimed at reducing the friction drag in turbulent boundary layers and channel flows. A cyclic spanwise oscillation of the wall with a proper frequency and amplitude is imposed, allowing a reduction of the turbulent drag of up to 40%. The present work is based on the numerical simulation of the Navier-Stokes equations in the simple geometry of a plane channel flow. The frequency of the oscillations is kept fixed at the most efficient value determined in previous studies, while the choice of the best value for the amplitude of the oscillations is evaluated not only in terms of friction reduction, but also by taking into consideration the overall energy balance and the power spent for the motion of the wall. The analysis of turbulence statistics allows to shed some light on the way oscillations interact with wall turbulence, as illustrated by visual inspection of some instantaneous flow fields. Finally, a simple explanation is proposed for this interaction, which leads to a rough estimate of the most efficient value for the frequency of the oscillations.     

Mechanism of drag reduction by spanwise oscillating Lorentz force in turbulent channel flow

Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.     

The role of elongational viscosity in the mechanism of drag reduction by polymer additives

The role of elongational viscosity in the mechanism of drag reduction by polymer additives is investigated qualitatively by means of direct numerical simulations of a turbulent pipe flow. For the polymer solution, a generalised Newtonian constitutive model is utilised in which the viscosity depends on the second and third invariant of the rate-of-strain tensor via an elongation parameter. This elongation parameter is capable of identifying elongational type of regions within the flow. The simulations show that complementary to stretching of the polymers, also compression must be incorporated to have drag reduction, contrary to many suggestions done in the literature on the mechanism which assume that stretching of the polymers is most important.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Surfactant systems for drag reduction: Physico-chemical properties and rheological behaviour

The first part of the work presents an overview of the physical chemistry of surfactants which in aqueous solutions reduce the frictional loss in turbulent pipe flow. It is shown that these surfactants form rodlike micelles above a characteristic concentraionc _t. The experimental evidence for rodlike micelles are reviewed and the prerequisites that the surfactant system must fulfill in order to form rodlike micelles are given. It is demonstrated by electrical conductivity measurements that the critical concentration for the formation of spherical micelles shows little temperature dependence, whereasc _t increases very rapidly with temperature. The length of the rodlike micelles, as determined by electric birefringence, decreases with rising temperature and increases with rising surfactant concentration. The dynamic processes in these micellar systems at rest and the influence of additives such as electrolytes and short chain alcohols are discussed.In the second part, the rheological behaviour of these surfactant solutions under laminar and turbulent flow conditions are investigated. Viscosity measurements in laminar pipe and Couette flow show the build-up of a shear induced viscoelastic state, SIS, from normal Newtonian fluid flow. A complete alignment of the rodlike micelles in the flow direction in the SIS was verified by flow birefringence. In turbulent pipe flow, drag reduction occurs in these surfactant systems as soon as rodlike micelles are present in the solution. The extent and type of drag reduction, i.e. the shape of the friction factor versus Reynolds number curve, depends directly on the size, number and surface charge of the rodlike micelles. The friction factor curve of each surfactant investigated changes in the same characteristic way as a function of temperature. For each surfactant, independent of concentration, an upper absolute temperature limit,T _L, for drag reduction exists which is caused by the micellar dynamics.T _L is influenced by the hydrophobic chain length and the counter-ion of the surfactant system. A first attempt is made to explain the drag reduction of surfactants by combining the results of these rheological measurements with the physico-chemical properties of the micellar systems.     

材料表面润湿性调控及减阻性能研究

设计合成不同结构的自组装分子,使其可以在不改变表面粗糙度的情况下改变表面的润湿性能;利用低表面能涂层修饰粗糙表面得到超疏水表面.采用流变仪和水洞试验分别在层流和湍流流动状态下测试了具有不同润湿行为的亲、疏水材料的减阻性能.结果表明:在层流流动状态,随着不同表面的接触角从13°增加到45°、113°和161°,减阻率随之从1.8%增大到7.2%、7.9%和14.9%;在湍流流动状态下,自组装涂层接触角为13°、45°和113°的三组模型的平均减阻率为0.8%、1.9%和6.8%,最大减阻率分别可达3.6%、9.2%和18.0%.两种流体流动中均存在材料表面水接触角增加减阻效率增大的行为.     

柔管自激振动特性与减阻效果的实验研究

为弄清柔管自激振动的湍流减阻效果的初步机理,在通过实验确认柔管确有湍流减阻效果的基础上,采用双重管结构和激光测位仪,对柔管的自激振动特性及湍流减阻效果及其两者的关系进行了实验研究．结果表明:管的壁厚越小,管外壁的自激振动的脉动位移有效值越大,湍流减阻效果越好;管壁外为压力平衡空气且雷诺数约为17500时,壁厚为2mm,3mm及4mm柔管的自激振动减阻率依次约为12％,10％,9％．这将为开发有效的管道流体输送技术提供参考．     

On the mechanism of drag reduction in dilute polymer solutions

Experimental evidence is given that drag reducing polymer molecules are preferentially collected by strained vortices. This can explain why extremely small amounts of additives can be so effective. They become concentrated in areas of a turbulent flow where they are most efficient.     

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一, 对缓解日益严峻的能源危机极为重要. 在重力式管道实验系统中, 测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态, 对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力, 沟槽结构采用机械方法加工, 并在表面喷涂超疏水涂层. 结果表明, 持续通气能解决超疏水沟槽表面气膜层流失问题, 实现气膜层长时间稳定维持; 恒定雷诺数下, 随通气速率增大, 超疏水沟槽表面气膜铺展更趋均匀, 减阻率上升; 由于通气速率影响气膜横向扩展能力, 致使恒定通气速率下, 减阻率随雷诺数的变化呈现两种模式; 在固定雷诺数及通气速率时, 减阻率随沟槽尺寸的扩大先增后减, $S^{+}\approx 76$时减阻率最大. 分析其原因在于, 沟槽结构增大沾湿面积的同时, 显著提升了通气状态下超疏水表面气膜层的稳定性, 因而展示出与超疏水表面和沟槽表面均不相同的减阻规律, 且效果更佳.