三维管内激波/湍流附面层干扰流场的数值模拟

给出了可压缩性及曲率修正两方程湍流模型，用它对跨音速三维喷管管内激波／湍流附面层干扰流场进行了数值模拟；给出了不同纵向截面上的激波结构、等马赫线图和四周固壁上的摩擦力线谱。将计算所得的摩擦力线谱、激波结构与实验结果进行了比较。分析了横截面上的流动结构，计算与实验吻合较好，进一步证实可压缩性及曲率修正两方程湍流模型能够较好地模拟激波／湍流附面层干扰流场，满足工程设计和分析要求。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

低速高湍流度90°弯管流动数值模拟

分析了曲率对弯曲管道流动的影响 ,给出了曲率修正的双层 k-ε湍流模型 ,并数值模拟了非均匀来流低速高湍流度 90°弯管内流动。文中比较了曲率修正双层 k- ε湍流模型的计算结果与实验结果 ,以及双层 k- ε湍流模型和逆压力梯度修正双层 k- ε湍流模型的计算结果。通过比较 ,发现经过曲率修正后的双层 k- ε湍流模型更好地模拟了非均匀来流低速高湍流度 90°弯管内流动。同时比较了壁面参数的插值方法 ,发现在弯曲管道流动模拟中对壁面参数进行一阶插值比零阶插值能更好地模拟真实流场     

跨声速风洞扩散段数值模拟

运用时间相关的全隐ＡＦ方法，计算二维及轴对称亚声速高速扩散段内流场。流动控制方程使用全Ｎ－Ｓ方程。湍流模型选择了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，和两层κ－ε二方程模型。经过同国外实验结果的比较，证明本文计算结果合理。在此基础上，计算了中国气动力研究与发展中心２．４ｍ×２．４ｍ引射式跨声速风洞第二扩散段在不考虑喷流干扰及常规运行状态下的流场。     

用隐-显法解超音速紊流流过凹角的流动

本文介绍一种用于对可压缩的Navier-Stokes方程求数值解的隐—显结合的数值方法,并用此法对二元的激波与紊流附面层干扰问题作了研究。计算表明此方法是可行的。 所采用的隐式格式为“delta”形式并保持了守恒,此格式只用于近壁面的法线方向上,而其余地区及沿物面方向均采用显式。为了减小非线性振荡,还采用了开关函数,此函数在整个流场均匀变化,使得格式在激波附近保持一阶精度,而其余地区仍为二阶精度。 紊流模型为简单代数旋涡模型,并对压力梯度及紊流记忆的影响作了修正。用此模型较好地预估了分离点的位置及压力升高。     

汽车气动力计算

用Ｆａｖｒｅ平均ＮＳ方程和Ｂ／Ｌ代数湍流模型对Ｓａａｂ９００汽车绕流流场进行了数值模拟,给出了汽车绕流流场的基本特征,探讨了车速对绕流流场,如车尾流动分离点,尾涡及汽车的阻力特性的影响,并建立了尾涡长度与汽车阻力系数之间的关系。     

竖直上升通道内微气泡湍流流动的数值模拟

为了提高气泡浓度分布预测精度,考虑流体和气泡之间的双向耦合作用,对竖直上升通道内含有微小气泡的湍流流动进行了数值模拟研究。采用直接数值模拟方法模拟流体流动,并用拉格朗日法跟踪气泡受流体和重力作用所产生的运动轨迹。气泡受力包括相间阻力、虚质量力、升力、壁面升力、压力梯度和重力等。数值模拟结果表明,通过加入壁面升力和采用双向耦合方法,可以较准确地预测气泡浓度分布和揭示气泡流动对流体湍流结构的影响。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型，采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式，对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟，探讨了内外喷流及喷管内的流动特征，并与实验数据进行了比较，计算结果与实验结果吻合较好，验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性，从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

激波穿过孔板系统的传播及其诱导的流场

基于Ｅｕｌｅｒ方程，采用高精度、高分辨率的ＧＲＰ有限差分法及算子分裂技术，对激波通过孔板系统的传播、绕射和反射，以及诱导的流场进行了数值研究。计算结果清楚地给出了激波传播、绕射、反射、及相互干扰的复杂流场结构。与实验结果比较，早期两者吻合较好，后期由于湍流发展的影响，必需考虑粘性修正。     

LEBU湍流减阻的展向扩展实验研究

利用大涡破碎器（ＬＥＢＵ）测量和考察了跨越平板湍流表面边界屋的速度型和动量厚度沿展向的分布，结果表明，ＬＥＢＵ的存在有可能扩大了层内流动参数沿展向的非均匀性；此外，还介绍了一个两对反向旋转的纵向涡模型，真实存在于风洞的横截面上，影响了展向非均匀性。使用本文介绍的一种沿展向非均匀性处理方法，对表面湍流摩阻估算的结果与使用摩擦阻力天平直接测量的结果基本一致。     

逆压力梯度下的自由流湍流度对湍流边界层特性的影响

本文对有逆压力梯度下的自由流湍流度对湍流边界层特性的影响做了较为仔细的测量,做了边界层内平均速度分布、壁面摩擦系数、压力梯度、边界层内三个方向脉动速度以及脉动总压的测量。试验结果指出:自由流湍流度的增加可大大削弱逆压力梯度对湍流边界层的不利影响,可使边界层速度分布饱满,尾迹强度减弱,从而延迟边界层分离。试验还首次得出自由流湍流度的增加可使边界层内的速度脉动和压力脉动分布均匀以及流动趋于稳定,因而可同时改善扩压器出口的静、动态特性。     

微型涡流发生器对超临界翼型减阻机理实验与数值分析

针对安装在超临界翼型后部的微型涡流发生器减阻问题，先用风洞实验测出微型涡流发生器对超临界翼型升阻特性的影响，然后采用RANS方程和κ-ε湍流模型进行数值模拟，分析安装在超临界翼型后部的微型涡流发生器减阻原因。研究发现:微型涡流发生器使下游近壁面处低能气体向上卷起与外层高能气体掺混，近壁面平均湍动能增加、翼型后部脉动压强增大，压差阻力减小；湍流应力由速度梯度、湍流粘性系数和脉动压强共同决定，虽然气流掺混，弦向速度法向梯度减小、湍流粘性系数减小，但展向速度法向梯度和脉动压强增大，湍流应力增大，摩擦阻力增大；微型涡流发生器尺寸很小，完全浸没于附面层内，仅掺混与它高度相当的附面层内流体，对附面层厚度影响小，对翼型升力影响小。     

新型高超声速进气道边界层人工转捩方法研究

为确保高超声速进气道的安全工作,其压缩面边界层在进入其内流道前必须完成转捩。针对高超声速进气道边界层转捩需要,依据二维高超声速边界层转捩机理,尝试了一种新型低阻高效的边界层人工转捩方法,在FD-07风洞中开展了试验验证。试验中首先通过进气道对称面压力分布和激波纹影获得进气道的自起动情况,进而推断进气道入口前的边界层转捩情况。试验包括进气道前体边界层自然转捩和人工转捩,试验结果表明在Ma=5、6,迎角α=0°来流条件下,使用同一波长的人工转捩带可以成功实现进气道边界层转捩,验证了基于线性稳定性理论设计的人工转捩带在宽马赫数范围的适用性。     

沟槽面对湍流边界层流动特征影响的实验研究

利用二维激光多普勒测速仪对零压力梯度下光滑面和小尺度沟槽面的沟槽及峰上部区域中湍流边界层流场进行了对比测量。笔者着重考察了在相同流动条件下两位置的平均速度、流向速度脉动强度、高阶矩以及雷诺剪应力的分布特性。实验中发现:沟槽能够增加粘性底层的厚度,减小壁面剪切力,有减阻效应。而此处的流向速度脉动的强度、高阶矩以及雷诺剪应力在距离壁面的不同区域中均有减小或降低,说明沟槽具有削弱湍流湍动的功能;而峰上部的湍流统计平均量则表现出与槽相反的结果或趋势。还讨论了该沟槽面减阻的原因。     

湍流边界层近壁区马蹄涡拟序结构的数值模拟与分析

采用大涡模拟对低雷诺数槽道湍流进行了数值模拟,进口边界条件采用了给定扰动波方法,经过短过渡段得到完全发展湍流.计算中采用了动态亚格子模型和预处理技术,进口采用的动量厚度雷诺数是670.大涡模拟对计算得到的平均速度剖面以及脉动分量与DNS解进行对比,验证了计算结果的可靠性.计算结果显示了湍流场中马蹄涡的形成及演化过程,其中包括形成单腿马蹄涡和形成对称涡腿的马蹄涡,同时发现流场中存在由亚谐波引起的拟序结构的交错现象,并在此基础上分析了湍流边界层近壁区马蹄涡结构的演化.     

平板边界层参数水槽测量与仿真分析研究

利用LDV测试技术，在小型水槽中对零压力梯度的光滑平板边界层进行了平均速度剖面测量。利用测得的速度数据进行从壁面到对数律层尾区全壁面律的拟合求解获得壁面摩擦速度和其他边界层流动参数。在实验测量之外还开展了平板绕流CFD仿真分析，并将实验结果、仿真计算结果和平板经验公式计算结果进行比较。结果表明，基于LDV的全壁面律拟合求解平板表面边界层流动参数具有较高精度，结合仿真分析，可为利用平板开展水中MEMS壁面剪应力传感器标定提供理想的输入。     

低湍流度风洞中的斜置翼三维边界层测量试验研究

在西工大低湍流度风洞中，于不同后掠角和湍流度下，以４５°单丝热线探头，在斜置翼上测量三维边界层速度型并计算特征参数的情况表明：横流影响的规律性明显，三维速度型扭曲相当复杂；后掠角和湍流度增大时，速度型和厚度的发展、不同层次的流线偏移和特征参数的变化受横流影响更为明显。     

高超声速圆锥边界层转捩纹影显示

为发展高超声速边界层转捩的试验研究手段，在中国空气动力研究与发展中心（CARDC）超高速空气动力研究所的FD-14激波风洞上开展了基于纹影显示技术的边界层转捩特性研究。试验模型为半锥角为7°的钝锥，头部钝度R_n有2.0和0.5mm两种。试验的名义马赫数为8和10，单位雷诺数变化范围为1.6×107~4.4×107m-1，高速相机的拍摄帧频20kHz。纹影显示捕捉到了第二模态不稳定波和湍流斑的空间结构。对纹影图像的灰度分布进行了功率谱密度（PSD）分析，结果表明第二模态波长约为边界层厚度的2倍。对纹影图像序列的分析表明，湍流斑波前传播速度大于波尾，并且略大于边界层外缘速度。     

湍流边界层的湍流参数测量和研究

利用DANTEC三维激光多普勒测速仪对水槽壁面边界层的湍流参数进行了测量和分析，并与经典实验曲线做了比较。同时就摩擦速度的计算与涡粘系数的分布曲线展开了讨论。