涡轮叶片型面换热系数计算

本文用数值方法计算了无气膜冷却涡轮叶片上包括前驻点的整个型面上的换热系数。驻点区采用相似方程计算, 并用圆柱绕流近似表示驻点附近流场。较成功地模拟了边界层转捩和过渡区的换热情况。计算中考虑了来流湍流度的影响。      

翼型动失速问题研究

本文对俯仰振荡和快速上仰翼型的非定常动失速问题进行了计算研究。对绕流流场采用了高效率的分区精确数值模拟的研究方法,保证了计算结果的精度,大大节省了对非定常粘性流动至关重要的计算机时。文中单独划出前缘驻点区,对形成前缘动失速旋涡的前缘附近区域的流场进行精确计算。算例显示了动失速流场的形成、变化过程及其对气动特性的影响,计算结果和实验的比较是满意的。     

飞船外形高超音速流场N-S方程数值解

用二阶迎风TVD格式求解N-S访方程,对简化的飞船外形高超音速粘性绕流进行了数值模拟。给出了物面压力分布及流场结构,与实验及其他计算结果符合较好。计算表明,肩部圆弧半径对后体流场影响较小,但严重影响肩部附近热流分布。当肩部圆弧半径是头部半径的20％时,肩部热流低干驻点处热流;半径为5％时,肩部热流超过驻点处热流30％。     

TRIP2.0_SOLVER的开发与应用

本文介绍了“亚跨超CFD软件平台”2.0版本流场解算器部分(TRIP2.0_SOLVER)的开发工作和应用情况,重点是结构网格部分的开发工作,包括对接网格、拼接网格、重叠网格等多种网格拓扑结构。TRIP2.0_SOLVER采用有限体积法离散雷诺平均的N-S方程,数值模拟绕流飞行器的空间流场和飞行器的气动特性,本文从软件已经具备的基本功能,用户界面的开发、软件测试和和多种网格拓扑结构的应用等四个方面总结了亚跨超CFD软件平台流场解算器(TRIP2.0_SOLVER)的工作进展情况。     

M∞→1时圆球绕前体流场的数值模拟

本文根据来流马赫数Ｍ∞选取坐标变换函数，将Ｍ∞→１时的低超声速圆球绕流前体流场变换到矩形的计算区域，忽略粘性影响，采用时间相关法，用ＴＶＤ有限差分格式求Ｅｕｌｅｒ方程的定常解，得到了Ｍ∞＝１．０５、１．０１和１．００５的流场分布。结果与弹道靶的实验吻合较好。     

钝体驻点线附近粘性激波层的数值方法

一、引言 在研究高超音速钝头体粘性气体绕流时,驻点区域的气动力参数和热力学参数的预报也是十分重要的。当粘性层很薄时,一直采用古典边界层方程计算;但是,考虑粘性气体与无粘流的相互干扰时,边界层模型已不适用,可用粘性激波层方程来描述极其广泛的来流马赫数及雷诺数的流场。 近年来,不少人致力于驻点线附近流场的计算,有的跨激波,有的分离激波。但     

二维斜入射定常驻点流解析解——干扰剪切流动理论的一个应用

首先证实二维斜入射驻点流为干扰剪切流(ISF),导出该ISF控制方程组,并求得该斜入射驻点ISF的解析解,它与完全NS方程组得到的精确解完全一致,对此作了相关物理分析,为CFD可信度验证的精确解比较方法提供了一个典型的局部区域流动模型,并为该流动模型的数值求解提供了有效的途径.     

激光等离子体点源减阻技术中入射能量对气动阻力的影响

数值模拟了来流马赫数为5时的圆球绕流冷流场,以及来流与头部驻点附近注入激光能量而产生的等离子体相互耦合形成的流场的演化过程。计算结果表明:施加局部能量点源能够改变原有的弓形激波结构,使其变为斜激波,从而减小气动阻力;气动阻力随入射激光能量的增大而逐渐减小,当入射激光能量为1.1 J时,气动阻力减小的百分比高达40%。     

完全催化壁驻点高超声速流动加热地面模拟方法研究

在地面高焓风洞中准确模拟高超声速头部驻点区化学反应流动加热,从平衡边界层和冻结边界层驻点气动加热公式出发,分析了离解焓对不同催化特性驻点加热的影响、地面风洞非平衡来流条件下头激波后的流场参数变化规律和与天上飞行条件的差异,分析了壁面催化特性对流场参数的影响,建立了地面高焓风洞模拟完全催化壁驻点气动加热的模拟准则:只要在风洞条件下实现对驻点压力ps、驻点速度梯度βe、来流总焓hs的模拟,则能够复现天上状态的气动热载荷。针对不同半径的飞行器头部驻点开展化学非平衡流场数值模拟进行了验证,表明:地面风洞无法完全模拟飞行器头部绕流流场参数,材料的催化特性对壁面附近流场中的N、O组元和壁面热流影响较大;地面风洞采用半径1∶1的模型模拟完全催化壁驻点加热会偏低于天上状态,而根据三参数模拟准则确定的模型尺寸能够同时复现天上状态驻点区的焓、压力、热流,且随着飞行器头部半径增大头部驻点线近壁面附近的温度、组元N、O梯度与天上会趋于一致。     

入射激波与飞行器头部激波干扰流场的计算方法

本文论述了入射激波与飞行器头部激波相交后的二维或轴对称干扰流场的计算方法。它是文[3]建立的具有最小耗散系数、保持解的单调性、一阶与二阶自动调节的混合格式的推广。通过对绕楔、球头和锥三种干扰流场的计算表明,本文的方法是可靠的,具有较高的精度,可用于多重间断问题的流场计算。     

三维机翼亚临界定常流的局部线性化有限元法

本文使用了两种类型等参元,通过变分法,计算了三维机翼亚临界定常位势流中的压力分布,计算结果与试验结果符合较好。文中还使用了有限元法中的局部线性化理论,插值出元素中心点密度去处理迭代过程中元素系数阵的计算,其结果与通常的有限元计算相吻合,但计算时间却明显减少了。文中还对有限元法局部线性化理论和局部线性化有限元法做了阐述。     

惯性颗粒在非均匀加速流中的输运问题

基于新型碳基颗粒增强或碳基纤维材料的再入体热烧蚀问题对飞行器的气动力/热性能影响机制尚不明确。初步研究不考虑化学反应和气体引射效应,把热烧蚀问题中的颗粒剥离问题模化为驻点流中的离散颗粒启动问题。由于高速气流在驻点温度达上万度量级,局部马赫数接近于0,研究近似假设靠近壁面的流体动力学为不可压缩的,并使用不可压缩颗粒解析的直接数值模拟方法（PR-DNS）研究一个惯性颗粒在平面挤压流和壁面驻点流中的动力学,从而揭示在驻点流中的颗粒输运机制。研究发现颗粒在壁面驻点流中的输运机制是水平输运,几乎没有垂直输运,与平行剪切颗粒床输运非常不同,这是驻点流特有的垂直向下挤压作用导致的。     

To the Problem of Friction Minimization on Permeable Surfaces at Supersonic Flow Rate

In this paper, we consider the variational problem of minimizing the Newtonian drag on an airfoil in a supersonic gas flow. A new condition for Lagrange multipliers on the streamlined airfoil was found that can be used to construct an optimal law for injection control in the analytical form. It is shown that for the flow around wedges, this control is finite at the stagnation point that is important in terms of practical realization of the optimal control.     

一种设计亚音速叶型的解析方法

本文在近似于势、流函数的正交面上采用解析方法进行亚音速叶型设计。将改写后的无旋条件与连续方程联立求解,可在计算面上将叶型设计问题用标准拉普拉斯方程来表示。这一方法的优点是在流场计算中用解析解取代数值迭代计算过程,求解简便,节省计算机时,并为从理论上分析各种不同叶型的设计特性提供了便利条件。文中提供的算例表明,方法是有效和可行的。      

Influence of coupled boundary layer suction and bowed blade on flow field and performance of a diffusion cascade

Based on the investigation of mid-span local boundary layer suction and positive bowed cascade, a coupled local tailored boundary layer suction and positive bowed blade method is developed to improve the performance of a highly loaded diffusion cascade with less suction slot. The effectiveness of the coupled method under different inlet boundary layers is also investigated.Results show that mid-span local boundary layer suction can effectively remove trailing edge separation, but deteriorate the flow fields near the endwall. The positive bowed cascade is beneficial for reducing open corner separation, but is detrimental to mid-span flow fields. The coupled method can further improve the performance and flow field of the cascade. The mid-span trailing edge separation and open corner separation are eliminated. Compared with linear cascade with suction, the coupled method reduces overall loss of the cascade by 31.4% at most. The mid-span loss of the cascade decreases as the suction coefficient increases, but increases as bow angle increases. The endwall loss increases as the suction coefficient increases. By contrast, the endwall loss decreases significantly as the bow angle increases. The endwall loss of coupled controlled cascade is higher than that of bowed cascade with the same bow angle because of the spanwise inverse ‘‘C" shaped static pressure distribution. Under different inlet boundary layer conditions, the coupled method can also improve the cascade effectively.     

涡轮叶栅前缘气膜冷却对气动参数影响的数值研究

采用具有TVD性质的三阶精度有限差分格式、自由型曲面技术以及多区域网格算法,对前缘带有三排冷气孔的涡轮导叶进行了气膜冷却数值计算,分析了前缘冷气喷射对涡轮气动参数的影响,描述了叶型表面冷气射流的运动规律。结果表明,冷气喷射导致了流量、马赫数和温度较为明显的变化,前缘和吸力面获得的绝热效率高于压力面,压力面冷气射流的运动规律比较复杂。型面压力只在冷气孔区域有明显的波动。     

应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

高超声速化学反应源项雅可比矩阵的对角化

基于Von Neumann稳定性理论,对高超声速化学反应流中源项雅可比(Jacobian)矩阵应用不同对角化方法的耦合隐式格式进行稳定性分析。详细推导出对角化方法情况下数值方法的增长因子,并分析不同对角化方法条件下的增长因子与CFL (Courant-Friedrichs-Lewy)数及波数之间的关系。理论分析与数值验证表明,Kim方法较Eberhardt和Ju方法更适合源项雅可比矩阵的对角化,保证数值求解在较大的CFL数条件下也能够实现稳定快速计算。另外,结合计算得到的刚性参数,稳定性分析结果还指出在高超声速化学反应外部流场中,激波与驻点之间的亚声速高温高压状态是导致数值格式计算不稳定的物理因素。     

飞机水平尾翼水滴撞击特性及防冰热载荷计算

建立了利用欧拉法求解水滴撞击特性的方法,并基于Messinger质量和能量守恒方程建立了热载荷计算模型,以某型飞机三维水平尾翼为研究对象,展开了撞击特性和热载荷的计算.结果显示水平尾翼的局部水收集系数极值从翼根到翼尖逐渐增加,不同截面水滴收集系数的分布为设计防冰系统范围提供了依据.驻点附近及机翼上、下表面各存在一个热载荷较大的区域,进而确定该尾翼电加热防冰分为3个区域:中间连续加热区和上、下表面驻点附近的各一个间断加热区.