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.     

亚临界圆柱绕流的DES方法比较

针对地面结冰试验设施喷雾系统尾迹湍流模拟需求,考察了3种脱体涡模拟(DES)方法对三维亚临界圆柱绕流预测的准确性,比较和分析了瞬时流动特征和流场统计量。研究发现:1)从瞬时流动特征来看,k-ω雷诺平均NavierStokes(RANS)控制的区域对剪切应力输运k-ωDES(SSTk-ωDES)求解的准确性有显著影响;2)从流动参数统计量的误差范围来看,回流区长度和流向最小速度的预测质量能在一定程度上反应圆柱绕流数值模拟的准确性;3)综合比较,SSTk-ωDES预测值与试验值和大涡模拟(LES)数据最为吻合,具有应用于地面结冰试验设施喷雾系统尾迹湍流计算的潜力。     

合成射流微扰动对后台阶湍流分离流动控制的实验研究

后台阶流动是流体力学中一个经典的研究课题,代表着工程中一类横截面突扩的钝体绕流问题。后台阶流动分离会导致一些不利的影响,如高速旋涡的形成、流动损失、压力脉动以及气动噪声等。基于阵列式合成射流激励器对二维矩形后台阶湍流分离再附流动控制进行了研究,综合应用表面测压、七孔探针、粒子图像测速仪(PIV)和热线等多种实验手段,获取了后台阶的表面压力分布和非定常流场结构。结果表明:利用在台阶前缘形成的合成射流微扰动可使无量纲再附点长度降低25%,合成射流控制使得沿台阶下游的湍动能和雷诺应力增强,提高了台阶下游流场的混合效率。热线结果表明,频率是后台阶分离流动控制的重要参数,当频率为260 Hz,扰动频率与剪切层涡脱落频率之比为1.32时,合成射流控制可使位于1/2倍频的剪切层能量增强,仅需消耗较小的能量即可实现流动控制的目的。     

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟,评估了Spalart-Allmaras, 标准k-ε, RNG (re-normalization group) k-ε,realizable k-ε,标准k-ω,SST(shear stress transport) k-ω,v2-f,Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确,其中SST k-ω模型比v2-f更能准确地预测壁面压力,然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

喷嘴内湍流运动对火焰浮起长度影响的数值模拟

为考虑喷嘴内部湍流运动对燃油雾化和火焰浮起长度的影响,将喷嘴内部的湍流流动以权重的形式加入初次破碎模型中,并对二次破碎模型进行了修正。建立了完整的燃油雾化和燃烧的数学模型。通过与实验数据对比来验证燃油雾化模型的准确性,并讨论了喷嘴内湍流运动对燃油雾化过程的影响。结果表明,湍流运动会加快液滴破碎和蒸发的速率,从而减小燃油蒸气贯穿距。火焰浮起长度的计算采用本文建立的燃油雾化模型,成功计算了火焰浮起长度随氧气体积分数、气体密度、气体温度和入射压力变化的规律。同时发现在不同气体密度和氧气体积分数的工况下,喷嘴内部湍流运动对火焰浮起长度的影响基本保持不变,分别为9%和13%;入射压力和气体温度的升高会导致喷嘴内部湍流对火焰浮起长度的影响逐渐变大。      

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

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

Large eddy simulation of turbulent diffusion flame combustion using a conserved scalar methodology

The present paper describes an LES prediction of turbulent diffusion flame combustion in a simplified axi-symmetric combustor geometry.The calculations are carried out using a well-tested finite volume incompressible LES code which has been modified to handle variable density and reacting flows.The basic mixture fraction conserved scalar method is used with the chemical state relationships described by fast chemistry.The turbulence-chemistry interaction is modelled by a sub-grid PDF method and the PDF is assumed to follow a Beta-function shape.The LES predictions have been time-averaged over 3.5 flow-through times to generate the mean radial profiles of mixture fraction,product mass fraction,temperature,axial velocity and axial rms.The agreement of the LES predictions with the experimental data is good for all the above quantities at four different axial positions with largest differences at the first measurement plane.The LES method also provides information on the unsteady nature of turbulent diffusion combustion. For turbulent reacting flows with large density ratio,it was found necessary to use a relaxation method in order to remove unphysical high-frequency fluctuations and to maintain numerical stability.      

近壁区Reynolds应力输运方程的理论模型

在湍流边界层的近壁区,采用对称共振三波的理论模型描述相干结构,根据理论模型对Reynolds应力输运方程逐项进行计算和分析,结果与直接数值模拟符合很好.这不仅在理论上有益于对湍流物理机制的了解,而且展现了一种可能性,即根据相干结构的理论知识来改造湍流模型,使之具有更清晰的物理内涵以提高近壁区的预估精度.     

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

超高速高焓流动研究进展

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