平板表面薄圆柱绕流摩擦力矢量场全局测量

针对壁面摩擦力矢量场测量问题,基于剪切敏感液晶(SSLC)涂层技术建立了一种测量平板表面摩擦力矢量场的方法。该方法基于多视角测量原理,采用六台同步相机从不同方向同时采集SSLC涂层在摩擦力作用下的颜色变化,与采用单台相机相比能够降低测量噪声,并且具有测量非定常流动的摩擦力场的潜力。应用该方法测量了平板表面薄圆柱绕流的摩擦力矢量场,结果表明:(1)SSLC涂层能够以彩色方式定性显示壁面摩擦力信息;(2)通过对不同方向观测的SSLC涂层颜色进行分析处理,该方法能够高分辨率测量薄圆柱绕流的摩擦力矢量场,详细地捕获了流动特征;(3)同一份SSLC涂层可以重复使用并且可用于测量不同的摩擦力矢量场。     

三维荧光油流技术的试验研究

荧光油流摩擦力场测量技术可以得到模型表面的全局摩擦力信息。该技术基于传统荧光油流技术发展而来,运用荧光油流技术作为原始数据获取手段,对试验图像运用HS光学流动算法进行图像处理得到模型表面油膜厚度随时间轴的变化量,依据动量定理对油膜厚度变化量进行计算得到当地的相对摩擦力信息实现表面摩擦力的可视化测量。本文运用该技术在低速风洞中对平板模型、75°平板三角翼、平板-翼型角区三维模型的表面摩擦力场分布情况进行测量,得到各模型表面的相对摩擦力分布和摩擦力线,并与平板模型Blasius层流解和三角翼模型经典流场结构进行对比。试验结果表明:在低速环境下该测量技术可以应用于模型表面摩擦力的可视化测量,所得摩擦力分布及表面流动情况基本可靠。     

绕方柱类凸起物分离流研究

本文发表了一组物面上按装方形凸起物引起流动分离的试验研究结果.全部试验在高超声速来流(M=5)条件下完成,流场中的方形凸起物在改变几何高度时,引起激波与边界层干扰流场的形态与物面压力分布规律的改变,说明方形凸起物的高度参数H/D在不同的范围内,对流场的影响不同.上述结果给出方形凸起物对应分离流场的特点,界定了竖直柱干扰区的分离边界与压力峰值.     

圆柱-平板角区湍流流动控制的风洞试验研究

角区流动中马蹄涡系的存在通常会造成不良影响。对圆柱-平板角区流动,在圆柱上游放置一倾斜的小圆棒能够改变角区流动结构。利用油流法和平板表面压力测量方法探讨了湍流流态下不同的小圆棒对平板表面的摩擦力线和压力分布的影响。油流实验揭示了倾斜棒能够改变角区的三维分离,新的分离线由倾斜棒和圆柱共同作用引起;倾斜棒对角区的作用可归类为两种不同的物理现象;倾斜棒能够引起圆柱侧面的分离线向下游发生极大的迁移,导致圆柱底部区域尾迹变窄。平板表面压力测量实验揭示附加的倾斜棒能够极大地改变压力分布情况,角区的逆压梯度相应减小;由此,逆压梯度引起的三维分离必然被削弱。     

局部凸起对可压缩平板边界层稳定性的影响

在采用直接数值模拟(DNS)方法计算得到来流马赫数为4.5的带局部凸起的可压缩平板边界层基本流的基础上,基于线性稳定性理论(LST) 进行稳定性分析,通过求解扰动方程研究了不同频率的扰动波与不同高度的局部凸起的相互作用,定义了穿透系数来研究局部凸起对流动稳定性的影响.研究结果表明:当局部凸起高度小于0.2个边界层厚度时局部凸起对来流的影响是局部的;局部凸起使频率小于最不稳定波频率的扰动波更不稳定,对频率大于最不稳定波频率的扰动波有稳定作用,其中最不稳定波频率是指局部凸起中心位置处平板边界层中第2模态最不稳定波的频率;穿透系数定量刻画了局部凸起对稳定性的影响,高度为0.2个边界层厚度的局部凸起引起e-N方法中N值的修正量达到0.8,接近转捩位置处N值的10%,在转捩预测中必须考虑.      

GRECO中棱边绕射场计算的改进

图形电磁计算 (GRECO)方法是计算复杂目标高频区雷达散射截面 (RCS)的有效方法之一。分析了原始GRECO方法在判定目标图象棱边象素的不足之处 ,给出了相应的改进措施。改进后的软件能够更准确、充分地判定目标的棱边象素及获得棱边参数。在边缘绕射场的计算方面 ,指出了相关文献中存在的错误 ,给出了基于等效电磁流法 (MEC)和物理绕射理论 (PTD)的边缘绕射场计算式 ,及与物理光学 (PO)场叠加求取RCS的完整表达式。计算实例表明 ,新的方法具有更高的准确度 ,与实验测量值吻合     

基于微型压力传感器阵列的边界层分离点检测方法

为了检测边界层的分离点,探讨了根据表面压力分布变化判定分离点的检测方法,提出了基于微机电系统(MEMS)技术和柔性衬底的微型压力传感器阵列的结构方案,可实现在线实时测量,并满足非平面表面的流体测量要求;同时提出了1种分离点检测的判定方法,并通过二维圆柱绕流仿真验证了该方法的正确性和有效性。微型压力传感器阵列的引入,拓宽了边界层分离点检测的解决途径。     

Improvements of Edges Diffraction Computing in GRECO

图形电磁计算(GRECO)方法是计算复杂目标高频区雷达散射截面(RCS)的有效方法之一。分析了原始GRECO方法在判定目标图象棱边象素的不足之处，给出了相应的改进措施。改进后的软件能够更准确、充分地判定目标的棱边象素及获得棱边参数。在边缘绕射场的计算方面，指出了相关文献中存在的错误，给出了基于等效电磁流法(MEC)和物理绕射理论(PTD)的边缘绕射场计算式，及与物理光学(PO)场叠加求取RCS的完整表达式。计算实例表明，新的方法具有更高的准确度，与实验测量值吻合。     

微型凸起物高超声速气动热特性研究

利用N-S数值解方法对高超声速气流中微型凸起物的气动热环境进行了计算研究,并对前斜坡峰值压力和热流跃升进行了分析,结果表明对于高雷诺数流动,由于其边界层很薄,气流对凸起物的冲击作用比低雷诺数流动大很多,给凸起物前缘和上表面带来较严重的气动热问题。研究还表明利用斜面侧扩张角可以降低凸起物后斜坡的气动热环境。最后对球锥与平板物型的差异给微型凸起物热环境特性带来的影响进行了研究,结果表明在本文计算状态下,风洞平板实验模型得到的凸起物前缘峰值压力和热流跃升比实际球锥外形要高,因此地面实验结果外推到天上实际情况还是比较保守的。     

运动激波绕射过程中的粘性效应

应用 TVD格式求解二维非定常全 N- S方程 ,以研究运动激波绕射现象。考虑粘性后 ,绕凹角流动时 ,压缩拐角附近表面压力分布与实验结果更为符合 ;绕凸起物时 ,出现反射激波因边界层干扰根部加宽 ,前、后角点附近有局部分离与再附 ,并引起局部压力上升 (与无粘流场模拟相比 )等现象。结果表明 ,全 N- S方程可较好地模拟实际粘性流中复杂激波障碍物干扰流场特性     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

矢量推力喷流对飞行器气动载荷影响的研究

矢量推力喷流对飞行器绕流具有干扰作用,引起飞行器气动载荷变化,对飞行器操纵性也有潜在影响。这种气动干扰和载荷变化是矢量推力技术研究的一个重要方面。与普通喷流相比,矢量推力喷流具有喷流偏转的特点,喷流偏转对气动干扰有重要贡献。因此本文采用多块搭接网格模拟了矢量推力飞行器无粘绕流流场,通过力系数和力矩系数对比,探讨了飞行器气动载荷随喷流偏转角的变化规律,这些规律具有理论和实际意义。     

圆柱绕流涡系演变的DPIV测试

利用数字式PIV技术在阵风环境风洞中对Re＝500时的圆柱绕流进行测量，给出了圆柱后2．2倍直径、圆柱两侧各一倍直径的区域内一系列瞬态的速度矢量场、涡量场和流线图。根据可视化图像和量化结果，着重分析了该区域的涡系结构和一个拟周期内旋涡的发展演化过程。     

探月返回器稀薄气体热化学非平衡特性数值模拟

通过发展直角/非结构混合网格下,适于复杂返回器含内能激发的五组元热化学非平衡稀薄气体直接模拟蒙特卡洛(DSMC)方法,数值模拟了钝体返回器稀薄流区7.5 km/s和10.6km/s再入速度下的非平衡流场特性和对物面气动力/热分布的影响。计算结果表明在月球返回速度下稀薄过渡流区存在着强烈的热化学非平衡,近连续滑移区仅在激波层区域存在着较大的平动、转动和振动非平衡度,而在高稀薄流区,热力学非平衡遍布绕流物体四周,包括整个压缩区和尾迹区。高温真实气体效应对表面热环境影响显著,对表面压力和摩擦力的影响相对较弱。计算结果显示探月返回器再入绕流过程100 km飞行高度仍需要考虑热化学非平衡的影响。     

侧向喷流干扰工程估算模型研究

喷流的直接力控制适用于全速域和全空域，具有响应快、效率高等特点，是控制飞行器飞行的有效手段之一。喷流与飞行器流场的相互干扰十分复杂，准确地获得喷流干扰气动力非常困难。笔者调研了国内外大量的喷流实验和局部凸起物干扰的分离流动实验，利用二维平板及轴对称凸起物绕流的分离特性，初步建立了计算喷流干扰的工程模型。该模型可计算喷流穿透高度、分离区长度、二维分离区边界、分离区平台压力系数，并可对分离特性进行三维修正。利用这一模型，笔者进行了某飞行器的超声速喷流的气动干扰特性计算，并与数值计算结果进行了比较。     

PIV技术在降落伞绕流实验中的应用

本文利用PIV技术在风洞中对五种不同类型降落伞模型的绕流空间流场进行了测量，得到了绕降落伞的空间流场结构。PIV测量结果可以为数值模拟降落伞绕流流场提供可靠的数据，为降落伞的改进优化提供指导和依据。     

破裂涡流中非定常现象与频率特性实验研究

通过流动显示、表面动态压力测量及热线测量等实验手段,对三角翼破裂涡流中的多种频率成分进行了分析。频谱分析确定了破裂点脉动和螺旋波的频率特征。实验结果表明,螺旋波主频随着弦向位置的增大先是迅速而后平缓减小。前缘涡破裂点振动具有准周期性,在不同的弦向位置上主频大小几乎没有改变,在靠近破裂点的位置有较大的振动能量。实验分析还表明,在破裂涡的流动状态下,虽然没有形成完全分离流,三角翼绕流流场中已经存在涡脱落的现象。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

热辐射计内稀薄气体流动的蒙特卡罗模拟

本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.