基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

基于非结构网格的涡轮气热耦合数值方法研究

利用课题组自主开发的三维非结构隐式N-S计算软件CU_Turbo,采用气热耦合计算方法,对MarkⅡ内冷径向涡轮导向叶片、带气膜冷却涡轮导叶MT1的流场和温度场进行了数值模拟。计算过程中,隐式时间推进中Jacobi ans矩阵采用对Roe通量的一种近似方法求解。结果表明,计算值与试验值吻合良好,验证了气热耦合计算方法的实用性和有效性,为涡轮工程设计提供了一种新的计算分析方法;涡轮叶片通道内附面层的不同流动状态及气膜冷却,对当地换热都有很大的影响。     

扑翼非定常运动对平尾影响及俯仰力矩特性研究

为了研究微型扑翼飞行器尾流对其平尾设计、飞行器稳定性以及飞行控制的影响,选取微型扑翼飞行器ASN-211为原模型,采用将其简化为二维的前后串列翼模型进行具体计算和分析。首先以Fluent动网格技术为背景,在用户自定义函数控制扑翼的非定常运动条件下进行不可压、非定常二维流动的计算,并研究在扑翼非定常运动条件下的模型的俯仰力矩特性。然后通过计算不同来流攻角、扑翼扑动频率、扑翼与平尾间距以及不同力矩中心下扑翼、平尾及总的力矩系数,讨论各个参数对力矩特性的影响。最后得出不同的扑翼扑动频率以及扑翼与平尾间距将会对平尾与扑翼的俯仰力矩间的相位差产生影响,所得结论为扑翼飞机的重心布置设计提供参考。     

结构网格方法对高升力构型的应用研究

采用"亚跨超CFD软件平台"(TRIP)数值模拟了梯形翼全展长与半展长高升力构型的复杂流场,并与试验做了对比分析。对应风洞试验是在NASA Langley 14×22英尺亚声速风洞(FST)和NASA Ames 12英尺增压风洞(PWT)中完成的。计算中采用一方程SA湍流模型和MUSCL-ROE格式,并综合运用对接/重叠/拼接网格技术,数值模拟了两种高升力构型的气动特性,给出了典型站位的压力分布,并对比研究了不同结构网格技术对此类高升力构型的计算差异。研究表明,与修正后的试验数据相比,数值模拟得到的气动力系数和典型剖面的压力分布均与试验结果吻合良好。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

复杂装备维修方式决策方法研究*

针对复杂装备维修方式难以决策的问题,提出了一种三角模糊数、层次分析和熵权相结合的模糊综合决策法,克服了传统模糊层次分析法准确程度不高以及在权值确定方面处理重要程度不确定情况的问题。首先在对不同装备选择维修方式考虑因素分析的基础上,建立了复杂装备通用维修决策准则体系和评价集。然后根据准则间以及准则与评价集之间对应的模糊判断矩阵求出层次分析与熵权权值,并由两权值计算得出综合权值。最后根据综合权值与模糊综合决策矩阵进行运算（综合决策集）,确定装备的维修方式。     

基于动网格法的翼型启动过程数值模拟

采用动网格法对NACA4421翼型以15°攻角启动过程进行了数值模拟。计算给出了启动过程中尾缘启动涡的生成、脱落与绕翼型环流充分发展的瞬态流场及气动力特性变化曲线,并对升力数据进行了拟合。计算结果表明,启动瞬间,上翼面最大负压和上下翼面最大压强差均出现在翼型后半段,随后逐渐向前缘移动,最终稳定在前缘点附近。下翼面最大正压点和上下翼面压差随弦向位置的最大变化率则始终维持在前缘点附近。加速过程中,整个翼型受到的升力近似于瞬时速度的二次幂指数的规律变化。加速段结束后,翼型转入匀速运动的瞬间出现升力小幅下降的现象,之后逐渐回升至稳定升力。     

复合固体推进剂颗粒的三维网格自动生成

为了分析圆柱颗粒对复合固体推进剂性能的影响,必须对复合固体推进剂中三维圆柱颗粒进行网格划分。根据"力平衡"网格划分原理设计了一种适合于三维圆柱颗粒的网格自动剖分方法。编程运行结果表明,方法能自动完成划分,网格质量较高,满足要求。     

构造特征函数拟合波浪曲面

自然景物的模拟是近年计算机图形学中最具挑战性的研究方向之一。作为其中一个重要组成部分,波浪的模拟受到了人们越来越多的关注。目前常用的一种方法是从求解水力学方程组出发来模拟波浪;另一种方法是通过构造特征函数来模拟波浪曲面。通过构造一个波浪函数,利用纹理映射技术可以快速高效地模拟波浪。     

高超声速流动分区对接网格算法研究

本文发展了一种能适应强间断流场的满足通量守恒的分区对接网格算法，并给出了保证通量守恒的通量插值系数应满足的条件。采用有限差分NND格式，求解层流NS方程，数值模拟了马赫数6—20绕圆柱的高超声速粘性流场。数值试验表明：方法具有良好的稳定性和收敛性，强激波能自由地穿越内边界。