脉动作用下波壁管内非牛顿流体流动特性的研究

在热量和质量传递过程中,流体的流动状态对整个传递过程有着重要的影响,为此,通过实验方法对以聚丙烯酰胺（PAM）溶液为代表的非牛顿流体在波壁管内的流动特性进行压力差和可视化研究。研究发现,PAM溶液存在一个最佳减阻浓度,在该浓度下与牛顿流体相比PAM溶液的转捩点明显提前。此外在脉动流场下,在一个脉动周期内,波壁管内流体存在明显的稳定流动和不稳定流动两种结构;相同净雷诺数条件下,脉动流场下波壁管内流体的流动混合情况比定常流动状态下强烈很多,意味着脉动流场具有更加优越的质量传递特性。     

基于气体射流的气液两相流动减阻特性

为了明确射流表面减阻特性,基于气液两相流动原理,建立具有射流功能的气液两相流动减阻计算模型,采用数值计算方法,通过RNG k-ε湍流模型对射流表面不同射流方向角、以及不同气体射流速度情况下的流场进行数值模拟及实验验证,分析气体射流对壁面黏性阻力、减阻率及流场的影响,研究气体射流的气液两相的减阻特性。结果表明:气体射流表面能够减小壁面黏性阻力,具有较好的减阻效果;四种射流方向角下,射流表面减阻率均在气体射流速度为5m/s时达到极大值,并且当射流方向角在xoz平面为-30°时,减阻率最大,为7.18%;数值模拟的减阻率效果略高于实验结果;气体射流的气液两相流动对壁面边界层进行了有效控制,减小了壁面的黏性剪应力和雷诺应力,导致壁面的黏性阻力减小,使得射流表面具有减阻性能。     

沟槽面湍流减阻数值评估方法

为解决雷诺平均Navier-Stokes(RANS)方法模拟沟槽面减阻的实际应用问题,受Wilcox对粗糙度模化方法的启发,通过深入分析k-ω剪切应力输运(SST)湍流模型近壁区ω值的作用效果,发现在黏性底层内增加壁面ω值可使对应近壁区的湍动能、湍流黏度、雷诺应力均减小,这种对近壁区流动特性的作用效果与真实沟槽面一致。以经典的对称V型沟槽面(高度=间距)减阻实验数据为基础,通过引入减阻效果影响因子——压力梯度与偏航角,建立复杂流动条件下沟槽面几何尺寸到壁面ω值的模化函数方程,将数值模拟结果与实验结果对比,验证了模化沟槽面可以达到与真实沟槽面相同的减阻效果,并依此给出沟槽面减阻在工程实践中的应用方法,以C919翼身组合体巡航构型为例,完成其沟槽面减阻的总体设计与评估。     

基于被动控制的D型钝体减阻大涡模拟研究

钝体减阻在航空航天等多个领域内具有重要的应用潜力,也是重要的基础研究热点问题。为了精确捕捉流动控制的流动细节,并发展效果优良的流动控制减阻方案,论文围绕抽象出的D型钝体,采用大涡模拟方法,开展了被动控制减阻高精度数值模拟研究。首先基于前期研究成果,对D型钝体尾迹区剪切层附近放置一个光滑小圆柱的被动减阻方法开展了数值模拟,发现总阻力减小17.7%,与试验结果吻合很好,同时数值预测的速度场分布也与试验结果吻合良好。在此基础上,进一步提出了采用齿槽型表面结构的小圆柱对D型钝体尾迹区进行扰动,并开展了数值验证,发现总阻力减小21.4%,优于前期的减阻方法。最后研究了在三种雷诺数工况下两种小圆柱扰动情况下的减阻效果,均表现出良好的减阻效果,两种小圆柱扰动下总阻力最大降幅分别为19.6%和23.1%,同时基于大涡模拟计算结果对减阻流动机理进行了探讨。上述研究结果表明,通过进一步优化流场结构,可以得到更优的流动减阻方案。     

小间距梯形扰流柱通道内的流动换热数值计算

为了深入了解涡轮叶片尾缘扰流柱区域的流动结构和换热情况,通过数值计算研究了小间距梯形扰流柱通道内的流动和换热.数值计算结果表明:①数值计算得到的端壁静压分布、弦向出流量分布、总压系数与实验结果符合良好,端壁换热分布规律总体上与实验结果一致.②扰流柱区域的流动方向基本为弦向.沿弦向湍流度是逐渐升高的,靠近弦向出口扰流柱后的湍流度比较低.数值计算结果对涡轮叶片内部冷却计算具有重要的参考价值.      

移动壁对压气机叶栅间隙流动的影响

采用分区网格处理方法并发展了基于分块结构的三维N-S方程流体求解器,用以模拟压气机间隙区域复杂的流动现象。以一压气机叶栅为研究对象,计算分析了端壁静止和端壁移动情况下叶栅内部的三维流场,并考虑了间隙大小的影响。通过详细分析叶栅通道内部及尾迹区域的流动并与已有的实验结果对比分析表明,程序较好地模拟出了移动壁效应对间隙区域二次流动形成和发展的影响。数值和实验结果均表明,端壁移动情况下间隙区域的二次流动特征是刮削泄漏涡的形成和发展。      

不同直径光滑圆管中黄原胶溶液流动减阻特性的实验研究

采用实验的方法,测试了黄原胶溶液在不同直径的光滑管道流动中的减阻特性和管径效应。其中变化参数为：（1）管道直径,共有3种管径,分别为5、10和20mm;（2）黄原胶溶液的浓度,变化范围为50~550ppm;（3）流速,用广义雷诺数（ReM）来表征,变化范围为0~50000。实验测量了不同参数下的管道压降Δp和体积流量Q,得到了黄原胶溶液浓度与减阻率的关系,以及ReM数对沿程阻力系数λ的影响,观测到黄原胶溶液减阻具有很明显的浓度效应,即随着溶液浓度的升高减阻率不断增大,直至达到最佳饱和减阻浓度,减阻率基本保持恒定。还观测到黄原胶溶液在不同管径中高低流速下表现出减阻类型的差异。     

大型飞机后体流动控制及减阻机理研究

通过数值模拟方法研究了大型飞机机身后体的流动特征和机理,并采用涡流发生器对机身后体流动分离进行控制,分析了其控制减阻的机理,讨论了不同参数对减阻效率的影响。计算结果表明:大型飞机机身后体流动分离是导致巡航阻力增加的一个重要因素;在后体底部附近安装涡流发生器,当来流迎角为负时减阻效果较明显,随着迎角的逐渐增大,减阻效果降低;涡流发生器的尺寸、周向距离、安装在机身后体上的前后位置、安装角度等参数对减阻效果具有显著影响。     

超声速内流道摩擦阻力分析及减阻技术研究

针对带有后向台阶的等截面受限空间,通过三维数值模拟开展了超声速内流道摩擦阻力分析及减阻技术研究。分析对比了飞行马赫数为5、6、6.5及7对应的燃烧室入口条件下相同质量氢气喷注、燃烧对壁面摩擦阻力的影响机制以及不同喷注压力对喷孔下游壁面剪应力的影响。研究结果表明,同等质量的氢气,低速喷注优于高速喷注(507、50.7 kPa喷注压力分别得到10%、5%左右的减阻效果)。近壁区燃烧得到接近70%的减阻效果;气流经过突扩结构之后,壁面剪应力呈现规律地不均匀变化,最大差异达100%;剪应力与密度变化趋势基本吻合。因此,发动机内流道减阻的关键在于营造近壁区低密度场;稳定、有效的减阻区域发生在靠后方的位置,但由于流动掺混、燃料的燃烧消耗,减阻效果沿流向逐渐减弱。      

对转压气机轮缘端壁抽吸流场特性数值分析及实验研究

以某双排对转轴流压气机为研究对象,通过区域缩放法开展非定常数值研究。数值结果显示近喘点附近转子2(R2)叶尖泄漏流动造成了高叶展范围内的流动阻塞以及高熵流动区域的产生,下游出口导流叶片承受着非常不利的来流条件。通过在转子叶尖近高熵流动核心区轮缘端壁处设置抽吸孔进行低能流体的移除能够明显改善叶尖区域的流场品质,提高压气机工作特性。根据数值研究结果,重新设计并加工了用于端壁抽吸的轮缘机匣,搭建了附面层抽吸系统及测控系统。在80%设计转速下,针对转子2叶尖进行端壁附面层抽吸,相对抽吸流量为1.0%。实验结果显示,端壁附面层抽吸应工作在低于原始特性最高效率点所在流量的小流量范围内。在此流量范围内端壁抽吸会大幅提高工作效率,改善压气机负荷水平。而在高于此流量的范围内,端壁附面层抽吸会带来负面影响。     

二维翼型混合层流控制减阻技术试验研究

选择NACA0006系列层流翼型作为物理模型,使用FLUENT商用软件计算分析翼型表面压力梯度,结合对翼型后缘做局部优化修形增大顺压梯度范围以及在翼型前缘布置吸气控制单元并配套吸气装置形成混合层流控制减阻技术。风洞试验中应用红外成像技术测量翼型表面层流区域,探索研究了混合层流控制减阻技术的实用效果。试验结果表明：对翼型实施混合层流控制减阻技术后,明显增大了翼型表面的层流面积。     

典型激波针减阻降热特性及流动机理

为研究激波针对超/高超声速钝头飞行器进行减阻降热的相关特性,采用数值模拟方法对6种典型激波针构型的流动特征开展系统研究。给出了激波针长度、来流马赫数对流动特征的影响规律,并对其形成机理进行了探讨。结果显示:马赫数较低时,头部有扰流物的5种激波针在回流区即将分裂时减阻率最大;马赫数较高时,减阻率在回流区分裂前后出现局部峰值,但最大减阻率将出现在回流区分裂后更长的激波针长度下。马赫数为3时,6种构型减阻率达最大时的相对长度在0.8~1.2之间,相比而言,球型、半球型和双锥型的减阻效果最好,最大减阻率为45%~50%;圆锥型最差,为20%~25%,明显低于头部有扰流物的构型。相同的激波针长度下,头部有扰流物构型的减阻率随马赫数增大而增大,圆锥型则相反。流场回流区、分离激波、弓形激波、局部膨胀流动等导致的压力分布变化是构型整体阻力变化的主要成因。      

槽道壁湍流的展向振荡电磁力减阻控制的PIV研究

壁湍流主要表现为条带和流向涡为主要特征的拟序结构和间歇性的湍流猝发事件,这些都会导致壁面阻力的增加,因此,为了实现减阻,需控制或消除壁面附近的流向涡进而抑制湍流的猝发。在该研究中,利用展向振荡电磁力对槽道湍流的近壁流向涡进行控制,以达到减阻的目的;并利用PIV测试系统对其进行了实验研究,讨论了这种电磁力的减阻效果及其减阻机理。结果表明：展向振荡电磁力具有减小壁面阻力的功能,其减阻机理为展向振荡的电磁力可以使条带倾斜,在流场中产生附加的负展向涡,导致近壁区域平均流向速度梯度的减小,因此,壁面阻力减小。     

Study of riblet drag reduction for an infinite span wing with different sweep angles

This paper studies the riblet drag reduction effect for an infinite swept wing under a low Reynolds number using a large-eddy simulation. The results show that the drag reduction ratio is not linear under different sweep angles. The maximum drag reduction ratio in this study is 9.5% for a wing with a 45° sweep angle. The local surface streamline angle and turbulence quantities are calculated to analyze the drag reduction mechanism. The results demonstrate that the riblets considerably suppress the Reynolds stresses above the wing upper surface, while the turbulence kinetic energy in the near wake is increased. A possible relaminarization phenomenon is observed at the middle part of the wing. Quasi-two-dimensional flow structures are observed near the wall, and a peak frequency is considered as the dominant frequency of the region.     

反旋进气混合式减涡结构流动特性数值计算

提出了一种将反旋进气孔与减涡管相结合的混合式减涡器,数值研究了其减阻引气效果,分析了旋转雷诺数、无量纲入口质量流量对内部流场结构和压力损失的影响。研究发现:在混合式减涡器引气结构中,静压沿径向平缓降低,在周向分布均匀;随着无量纲入口质量流量或旋转雷诺数的增加,引气结构总压降呈现单调上升的趋势,其中在高旋转雷诺数、低无量纲质量流量工况下具有突出的减阻性能,其对应的湍流参数为0.106 4~0.324 5。相比于简单盘腔,反旋进气孔式及管式减涡器的压力损失分别降低62.5%、60.5%,混合式减涡器可降低80.4%,体现出良好的减阻引气效果。      

Experimental and numerical investigation on opposing plasma synthetic jet for drag reduction

Experimental and numerical studies are carried out to validate the potential of opposing Plasma Synthetic Jet (PSJ) for drag reduction for a hemisphere. Firstly, flow field changes of opposing PSJ are analyzed by comparing the experimental schlieren images and simulation results in a supersonic free stream of Mach number 3. As PSJ is a kind of unsteady pulsed jet, the shock standoff distance increases initially and then decreases under the control of PSJ, which corresponds to the change of the strength of PSJ. Accordingly, the amount of drag reduction of the hemisphere increases initially and then decreases. It is found that there is a short period of “drag rise” during the formation of PSJ before the drag reduction, which is induced by the generation of normal shock waves and the area difference of the cavity wall of PSJ Actuator (PSJA). Secondly, the effects of five parameters, including exit diameter, discharge energy of PSJA, Mach number, static pressure of incoming flow and angle of attack, on drag reduction of opposing PSJ were studied in detail by using numerical method. It is found that the Maximum Pressure Ratio (MPR) has a significant impact on the average drag reduction for a configuration-determined PSJA. For the configuration selected in this study, the flow field of opposing PSJ shows typical Short Penetration Mode (SPM) in a control cycle of PSJ when the MPR is less than 0.89. However, the flow field shows typical Long Penetration Mode (LPM) at some time when the MPR is bigger than 0.89. Relatively better drag reduction is achieved in this case.