飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

Numerical simulation of compression corner flows at Mach number 9

A hypersonic vehicle encounters a wide range of conditions during its complete flight regime. These flight conditions may vary from low to high Mach numbers with varying angles of attack. The near-wall viscous dissipation associated with flows at combined high Mach and Reynolds numbers leads to significant wall heat transfer rates and shear stresses. The shock wave/boundary-layer interaction results in a flow separation region, which commonly augments total pressure losses in the flow and lowers the efficiency of aerodynamic control surfaces such as fins installed on a vehicle. The standard turbulence models, when used to resolve such flows, result in incorrect separation bubble size for large separated flows. Therefore, it results in an inaccurate aerodynamic load, such as the wall pressures, skin friction distribution, and heat transfer rate. In previous studies, the application of the shock-unsteadiness correction to the standard two-equation k-ω turbulence model improved the separation bubble size leading to an accurate pressure prediction and shock definition with the assumption of constant Prandtl number. In the present work, the new shock-unsteadiness modification to the k-ω turbulence model is applied to the hypersonic compression corner flows. This new model with variable Prandtl number is based on the model parameter, which depends upon the local density ratio. The computed wall pressures, heat flux and flow field are compared to the experimental data. A parametric study is carried out by varying compression deflection angles, free stream Reynolds number and wall temperatures to compute the flow field and wall data accurately, particularly in the shock boundary layer interaction region. The new shock-unsteadiness modified k-ω model with variable Prandtl number shows an accurate prediction of initial pressure rise location, pressure distribution in the plateau region and heat flux in comparison to the standard k-ω model.     

时间矢量推进方法及其在叶轮机非定常流分析中的应用

为提高周期性非定常流动的求解效率,将非定常计算的初值问题转换为边值问题,提出了时间矢量推进方法,并将该方法应用于叶轮机多排的非定常流动分析中。通过对两排对转风扇进行非定常仿真,并对比双时间步计算方法的计算结果,讨论了新方法的计算效率,研究了该方法对排间干扰捕捉的准确性和可靠性。得到了以下结论:在物理时间步长相等的情况下,新方法与双时间步方法的计算结果基本一致,且计算时间约为双时间步方法的1/8;时间矢量推进方法能够很好捕捉排间的势扰动、熵扰动和涡扰动以及主流和扰动之间的非线性作用;时刻样本数较少会使时间矢量推进方法捕捉到的非定常变化幅值变小,且无法解析时间尺度较小的非定常流动现象。     

超高速高焓流动研究进展

本文介绍了中国科学院高温气体动力学重点实验室在超高速高焓流动模拟技术和试验方法方面取得的研究进展.文章主要包括三部分研究内容:第一部分是关于发展先进的超高速试验模拟技术,包括爆轰驱动高焓激波风洞和爆轰驱动高焓膨胀管.高焓激波风洞产生的超高速气流速度的范围是3.5km/s～6.0km/s,高焓膨胀管能够模拟速度为6.5km/s～10km/s的超高速气流.第二部分介绍高焓激波风洞喷管流场诊断结果,用来检验喷管产生的超高速流场的流场品质及其与飞行条件的差异.第三部分是关于超高速流动的试验方法和数值技术研究,包括高焓流动中真实气体效应对飞行器俯仰力矩变化的影响;热化学反应流动中表面催化效应诱导的气动热变化规律;喷管流场的气流非平衡效应对试验结果可能产生的影响.     

倾斜单孔射流轨迹与速度分布的研究

本文采用实验和数值模拟的方法,在吹风比分别为0．5、1．0、1．5、2．0的情况下,研究了圆柱单孔二次射流的贴壁性及速度分布。通过对比不同截面位置的流动轨迹和速度分布,发现圆柱孔射流流场的速度分布与圆柱扰流流场的类似。射流出口的最大速度并不位于射流孔中心位置,u值也不是在射流孔中心位置最高。     

微重力下容器圆形倒角处的毛细驱动流

空间探索计划的需求对空间流体管理技术提出了新的挑战.流体管理设备的内角是主要的流体传输"管道",因此研究微重力下容器内角处的毛细驱动流具有重要的意义.由于设计和制造等原因,完美的尖角并不常见,而是带有一定圆角过渡的内角.笔者设计不同的带有圆形过渡内角的容器,选用10cSt硅油为实验流体介质,通过一系列的落塔实验,研究了容器内角处不同的圆角半径对毛细驱动流的影响,并分析得到了有价值的实验规律.     

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

二元内流通道内湍流附面层分离的数值与实验研究

本文应用激光多普勒测速仪测量了二维扩压器中不可压湍流附面层分离流动,得到了时均速度和表面压强分布,表面摩擦系数和时均速度分布经验关系式做了比较。预测中采用准同时粘性/无粘流匹配模型,并比较了两种形状因子H_1经验关系式所得预测结果。在最大反流以前,所预测的附面层位移厚度及形状因子与实验吻合。     

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。