数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

带有双层壁扩压器的波瓣喷管混合流场的数值计算

采用三维贴体曲线坐标网格,边界网格加密且正交,在整个计算区域进行全场计算。为避免因波瓣造成的网格强烈的非正交引起解的发散,采用了Chen和Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法,对带有双层壁扩压器的波瓣喷管混合流动进行了数值计算和分析。对波瓣及双层壁,采用大粘性的方法解决流固耦合。计算结果表明:在双层壁间有外界冷气流被引射进入,形成了壁面的冷却气流;自波瓣出口截面沿流向产生的环流速度场,强化了主次流的掺混,速度分布渐趋均匀。计算结果与实验数据符合良好,二者在离开波瓣52 mm的混合管内相对主流速度的最大误差为21.65 %。     

计算矢量喷管流场的有限元法

本文通过采用扫描有限元法及多步龙格－库塔时间推进格式求解非定常平均Ｎ－Ｓ方程来模拟二元矢量喷管的内部跨速流场,成功地计算了二维有激波的跨声速喷管内部流场及仰 仰角３０°的三维矢量喷管内部流场。     

利用微气囊扰动进行流动控制的数值模拟研究

探讨一种可以用来模拟“三角翼上分布微气囊，从而控制流动获得滚转力矩”的数值方法。为微型飞机的气动设计提供一种工具。研究的内容包括：考虑微气囊的三角翼网格生成、流场NS方程计算、微气囊不同布局对流动的扰动等。     

民用飞机颤振试飞技术研究

为了能够顺利完成ARJ21-700飞机颤振试飞,研究民用运输类飞机颤振试飞的依据与要求,提出民用运输类飞机颤振试飞的测试方法、试飞方法、激励方法和数据处理方法;结合ARJ21-700飞机合格审定试飞进行试飞验证,给出该型飞机颤振特性试飞的符合性验证结果。结果表明:该套方法适合民用运输类飞机合格审定颤振试飞。     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

非定常风沙流中风速脉动对沙粒相瞬时浓度影响的实验研究

在大气边界层风洞内采用高速CMOS相机记录了时均风速线性上升的来流条件下运动沙粒的图像。根据二进制图像标记像素的连通性,从数字图像中提取相机拍摄区域内沙粒的数量N,由此计算得到风沙流的瞬时含沙量C。获取图像的同时,使用恒温式热线风速仪MiniCTA配合55R49热膜探头测量观测域上方主流位置和跃移层顶部附近的风速变化。实验结果表明,风沙流中的瞬时含沙量总是非稳定的,其变化频率至少在100Hz以上。风沙之间的相互作用存在一定的滞后效应,大气湍流脉动对含沙量的作用与沙粒粒径相关,粒径越小,风速脉动作用越大。当粒径≤160μm时,风速的微小变化即可引起含沙量的明显波动,二者高度相关。当粒径较大的时候,如300～500μm,只有较长时间较大幅度的风速脉动才能使含沙量有较大幅度的变化,研究时均风速变化对含沙量的影响更具有意义。     

对流层折射修正中水汽压公式对比研究

针对当前对流层大气折射效应评估计算的工程应用中,水汽压的计算方法不尽统一、精度参差不齐的现状,分别在-50℃～0℃气温范围内利用各平冰面饱和水汽压公式、在0℃～50℃气温范围内利用各平液面饱和水汽压公式对大气折射率湿项进行了仿真计算,通过对仿真结果的比较可知：在平冰面和平液面情况下,Buck（巴克）1981公式和Buck1996公式与Goff-Gratch（戈夫-格雷奇）公式的偏差均较小,计算精度较高,建议作为测控系统外测数据处理中大气折射修正的首选经验公式;该分析结果也可为其他需要考虑水汽影响的工程领域提供参考.