DLR-F6翼身组合体数值计算

使用自行研制的混合网格亚跨超声速流场解算器程序MFlow计算了AIAA第三届阻力会议提供的DLR-F6及其带整流装置FX2B情况下的阻力。重点分析了两者的网格收敛特性、阻力极曲线以及压力分布等,并与阻力会议的各个软件的计算结果进行比较。详细分析了DLR-F6翼身结合处后缘附近网格精度对分离气泡计算的影响。计算结果表明,MFlow程序的计算精度与国外软件相近,能够为阻力计算提供比较精确的结果。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。     

尖拱弹身大攻角N-S方程数值计算

本文采用有限体积的NND格式,数值模拟了弹身超音速有攻角绕流,分析了周向网格和法向网格对弹身背风区分离流态及气动性能的影响,计算与实验结果比较,两者相当一致.     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

发动机进排气效应对民机构型气动特性影响

基于非结构混合网格技术,通过数值求解Navier-Stokes方程,分析了航空发动机进排气效应对民机构型气动特性的影响。通过设置合适的进排气边界条件来模拟发动机进排气效应的影响,采用单独发动机短舱风洞试验模型,验证了计算方法的可靠性。在此基础上,分析了发动机进排气效应对翼吊式和尾吊式两种典型民机构型气动特性的影响,结果表明:翼吊式民机构型发动机进排气效应对升力的干扰主要由发动机尾喷流的引射效应对机翼的干扰引起;尾吊式民机构型随着发动机进气流量的增大,机翼上表面压力逐渐减小,激波位置逐渐后移;不同发动机进气流量会对飞机的阻力特性产生较大影响,在飞机详细设计阶段需充分考虑机体与动力装置之间的干扰影响。      

一个非常规前体机身的流动显示研究

描述了有类似Ｅｒｉｃｋｓｏｎ前体的非常规机身的初步流动显示研究结果,包括测力、表面油流和激光片光流动显示结果和用层流Ｎ－Ｓ方程进行数值模拟的结果。计算和试验的参数范围：α＝０°～５０°,β＝０°～２０°,虽然计算与试验所用的外形在后部有一些不同,但是两者在涡的位置方面显示了良好的一致性。同时研究也表明,大攻角的流动特性可以通过改变机身前体形状进行控制。通过研究还表明,这类前体在改善大攻角横侧方向安定性方面具有很大的潜力。     

非结构网格体心梯度求解方法的精度分析

对于空间二阶精度的非结构网格求解器,无黏通量和黏性通量的计算都需要流场变量在单元体心处的梯度值。体心梯度求解的精度很大程度上决定了求解器整体的精度。本文推导了几种常用的梯度求解方法(包括格林高斯法,节点格林高斯法和最小二乘法)求解体心梯度的理论精度,并从数值方面对上述理论精度进行了验证。在数值测试过程中,通过以当前单元体心为基准进行坐标局部缩放的做法,保证了在非结构网格上做精度测试时,网格拓扑能够严格保持不变。理论推导和数值测试一致发现:使用保线性权插值到节点的节点格林高斯法、使用共面单元的最小二乘法和使用共点单元的最小二乘法,不管网格拓扑关系如何,都能保证梯度求解为一阶精度。而格林-高斯法和使用距离倒数权插值到节点的节点格林-高斯法在一般的网格中只有零阶精度,只有在特定的网格和权函数关系下,才能有一阶精度。最后用NACA 0012翼型亚声速绕流模拟进一步验证了梯度求解方法的收敛精度。