低速长航时无人机垂直突风过载分析

低速中空长航时无人机的升力非线性较强,利用标准离散突风模型和连续紊流模型计算得到的突风过载差别较大.针对该问题,本文采用改进的离散突风模型求解突风载荷,即直接求解突风条件下原始的运动学方程,并在运动历程中考虑非线性升力的影响.结果表明:改进的离散突风模型计算得到的突风载荷明显低于标准离散突风模型所得的载荷,更贴近无人机的实际情况,且改进的离散突风模型与连续紊流模型的计算结果也相近.     

脉冲爆震发动机旋流式气动阀工作机理和扰流器阻力特性的实验研究

在脉冲爆振发动机（Pulse detonation engine,PDE)上采用旋流式气动阀,并获得成功应用。在此基础上,通过测量旋流式气动阀正反向流动阻力系数,研究了旋流式气动阀的工作机理。同时,还研究扰流片形状、阻塞比、扰流片片数及间距对扰流器阻力系数的影响。研究表明：PDE旋流式气动阀反向流动阻力系数仅比正向流动阻力系数大5％～36％,并通过试验结果建立了PDE旋流式气动阀工作机理。试验中发现：扰流器的阻力系数和扰流片型式、堵塞比、扰流片数目和间距有关。这些研究结果对PDE旋流式气动阀和扰流器设计有较大的参考价值。     

具有地面效应理想车体三维绕流的研究

考虑理想车体与地面相对运动的基础上,建立地面效应模型。采用标准k—ε两方程湍流模型对理想车体外流场进行数值模拟。并对数值模拟结果进行了研究分析。     

Numerical simulation of airfoil flows with a turbulence model

Turbulent flows over AS240 and NACA4412 airfoil were simulated numerically using a two-equation turbulence model named k-ξ model. The predictions of velocity profiles and the pressure coefficient of airfoil AS240 at 8°/19° attack angle and NACA4412 at 13.87° attack angle were calculated. The results were compared with those using k-ε and k-ω models,as well as experimental data. It indicates that the new k-ξ model offers more realistic prediction than the other two models. The main finding shows that the new k-ξ model is good at predicting separated flows around airfoils, and it captures the flow feature of pressure-induced separation adequately. All calculations are implemented as per openFOAM 1.7.1(open source field operation and manipulation).      

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

主流湍流度对涡轮导向叶片气膜冷却特性影响的实验

采用基于窄带热色液晶的瞬态全表面传热测量技术,研究了主流湍流度对涡轮导向叶片吸力面圆柱形孔排气膜冷却特性的影响规律.结果表明:在实验工况范围内,主流湍流度从0.59%提高至6.85%,可以在气膜出流的上游区域促进气膜贴向壁面并扩大展向覆盖面积,从而改善气膜覆盖效果,但是在主流湍流度较大的工况下,气膜覆盖效果迅速变差;在气膜出流的下游区域,主流湍流度的提高使得气膜冷却效率逐渐降低;主流湍流度的增大,增强了无气膜冷却光滑叶片表面的对流换热;在气膜冷却条件下,气膜出流对叶片表面对流换热的增强效果随着主流湍流度的增大呈现出明显的区域性特点:表面传热系数比在上游区域是先增强后减弱;中游区域是逐渐减弱;下游区域则是逐渐增强.     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

基于SA和SSTγ-■模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-■转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-■转捩模型预测的气动力特性和流场分布。结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-■转捩模型在流动细节上依然存在着一些差异。