二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

高速磁浮列车会车压力波数值计算研究

根据可压缩粘性流体的N-S方程和k-ε双方程湍流模型,运用移动网格技术,采用有限容积法对磁浮列车以500 km/h速度运行时的会车压力波进行了数值计算研究,计算结果与现车试验结果相差10%以内,表明了本文计算模型和方法的正确性;文中还研究了交会过程中列车周围的流场分布,结果表明:当列车交会时,两车之间的流场相互叠加,而两车之外的流场分布在整个交会过程中变化不大。     

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

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于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°攻角启动过程进行了数值模拟。计算给出了启动过程中尾缘启动涡的生成、脱落与绕翼型环流充分发展的瞬态流场及气动力特性变化曲线,并对升力数据进行了拟合。计算结果表明,启动瞬间,上翼面最大负压和上下翼面最大压强差均出现在翼型后半段,随后逐渐向前缘移动,最终稳定在前缘点附近。下翼面最大正压点和上下翼面压差随弦向位置的最大变化率则始终维持在前缘点附近。加速过程中,整个翼型受到的升力近似于瞬时速度的二次幂指数的规律变化。加速段结束后,翼型转入匀速运动的瞬间出现升力小幅下降的现象,之后逐渐回升至稳定升力。     

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

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