N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

空腔噪声的马赫数敏感性研究

敏感度分析在评估参数的重要程度以及计算不确定度方面具有重要作用。通过风洞试验开展了亚跨超声速下的空腔噪声马赫数敏感性研究。亚跨声速下，通过调节调压阀的开度改变马赫数，马赫数名义增量为0.010。超声速下，通过改变模型迎角实现马赫数的连续变化，迎角增量为1°。利用总压耙测量空腔入口马赫数。结果表明:空腔后部测点脉动压力系数在亚跨声速下随着马赫数的增加而增加，而在超声速下随着马赫数的增加而减小。跨声速下，脉动压力系数对马赫数的敏感性导数最大。不考虑模态切换的情况下，不同速域的主导声模态St对马赫数的敏感性导数均为负数。主导声模态谱峰在亚声速下随着马赫数的增加而增加，而在超声速下随着马赫数的增加而降低。敏感度研究结果不仅可用于内埋武器舱气动噪声载荷的不确定度评估，也有助于更好地认识空腔噪声特性。     

基于隐式嵌套重叠网格技术的阻力预测

 采用一种多层多块隐式嵌套重叠网格技术,对美国国家航空航天局通用化研究模型(NASA-CRM)翼身平尾(WBT)组合体进行了数值模拟与分析。多层多块隐式嵌套重叠网格技术是结合多层多块嵌套重叠网格处理策略和隐式切割方法,在建立重叠网格之间的流场信息传递关系时,基于网格单元切割准则选择"最优"重叠单元而无需人工设定插值边界。对美国AIAA委员会召开的第4届阻力预测研讨会(DPW-4)提供的CRM WBT组合体生成4种不同密度的结构化多层多块嵌套重叠网格,并采用计算流体力学(CFD)方法进行数值计算和阻力预测,计算结果与CFL3D和OVERFLOW的结果进行了对比。数值模拟结果表明:计算得到的压力分布和极曲线与CFL3D和OVERFLOW的结果几乎相同,说明了隐式嵌套重叠网格技术的有效性,同时也验证了流场求解方法与程序的可靠性。当迎角增大到3°左右时,在机身与机翼、尾翼连接处出现明显的分离涡,影响CRM WBT组合体的气动特性。在阻力预测方面,增加网格密度能够提高阻力预测的精度。采用不同的湍流模型会导致升、阻力系数的计算结果存在一定的差异,因此,湍流模型的选择也是阻力预测需要考虑的因素。     

新型等离子体激励器的数值模拟

为了对新型多电极对等离子体激励器进行可靠的数值模拟,通过求解带有体积力源项的Navier—Stokes方程,模拟单电极对等离子体激励器在二维不可压流中对周围流场的影响,计算结果表明所使用唯象模型的参数基本正确。在此基础上,改进原有模型,对新型双极性等离子体激励器诱导流动进行计算。通过与实验结果的对比,发现改进后的模型能较好地模拟新型双极性等离子体激励器的控制作用,可为今后新型等离子体激励器的应用提供可靠的数值模拟方法。     

基于改进CST参数化方法和转捩模型的翼型优化设计

为提高翼型优化设计效率,增大设计空间,采用B样条基函数替代传统的形状类别函数(CST)方法中的Bezier多项式,增强了对翼型参数化表达的局部控制能力并提高了翼型局部表达精度。为了确保翼型在优化设计过程中的几何光顺特性和代理模型的准确性,采用小波分解技术提出了多分辨率翼型的局部光顺方法。采用基于本征正交分解(POD)的流场数值代理模型,并结合γ-Reθt转捩模型实现了快速准确的气动力与流动转捩预测。采用小波技术光顺的CST翼型参数化建模、POD流场数值计算代理模型以及γ-Reθt转捩模型,结合遗传算法建立了完整的翼型气动优化设计系统。针对低速层流翼型与超临界翼型进行优化设计,优化设计后的翼型升阻比分别提高了47.42%和45.85%,且对改进前后参数化建模方法的优化性能进行了对比,结果表明本文构建的翼型气动优化设计系统具备很高的优化效率。     

新概念轮翼气动特性实验研究

研究一种关于小型飞行器的新型升力和推力系统——轮翼的气动特性。通过自行设计并搭建的基于虚拟仪器的实验平台,进行地面试验,试验主要研究了轮翼转速对系统气动特性的影响。结果表明：随着轮翼转速的增加,系统提供的升力及推力增大。由于轮翼系统消耗的能量绝大部分用来产生升力及推力,因此轮翼系统具有较高的工作效率,为进一步研制小型无人直升机探索新型动力模式。     

蔡晋生 流体力学、飞行器气动布局设计专家

请您简要介绍一下计算流体力学(Computational Fluid Dynamics,CFD)以及飞行器气动布局设计技术在航空航天领域的应用情况.     

战斗机后体流场数值模拟与减阻优化设计

采用数值计算方法模拟零迎角跨声速来流下战斗机后体流场,并进行后体外形的减阻优化设计。通过数值求解二维轴对称Navier-Stokes方程、k-ωSST湍流模型和气体组分方程,研究战斗机后体绕流与尾喷流相互耦合的流场特性,对三种欠膨胀喷管压比下的喷流进行数值模拟,同时对比多组分气体喷流和理想气体喷流对后体阻力的影响。采用梯度法对轴对称后体外形进行减阻优化设计,提出一种优化设计加速算法,其基本思路是通过逐步增加设计控制点个数并根据外形曲率合理分布设计控制点的位置,从而改善梯度法的优化效率。计算结果发现,超声速喷流会在后体尾部附近形成复杂的波系结构;与组分气体喷流相比,采用理想气体作为喷流介质时的后体阻力系数略高;在跨声速来流状态下,后体阻力系数值随喷压比的增大而减小。优化结果显示,优化后的的后体阻力系数可以降低13%左右;与一次性均匀分布优化控制点的梯度法优化方法相比,采用提出的优化加速算法可以缩短优化计算时间40%左右,并且可以提高优化设计的精度。     

CFD数学模型的线性化方法及其应用

计算流体力学(CFD)方法不仅仅起到数值模拟的作用,它本身是一个复杂的非线性系统。在流动稳定性分析、气动弹性分析、优化设计以及流动控制等领域,从系统的角度出发,对CFD数学模型线性化后,可以对模型的系统矩阵进行定量分析获得更多的系统特性。但是CFD数学模型往往非常复杂且阶数很高,因此其线性化系统矩阵的获得比较困难。鉴于此,采用人工编程和自动微分相结合,构造有限体积法并行CFD模型的线性化系统矩阵。其中自动微分只被用来得到每个界面通量的局部雅可比矩阵,而采用人工编程方法来实现并行环境下的稀疏雅可比矩阵的组装。线性化系统的并行求解采用了块雅可比预处理的广义最小残量法,每个并行进程内部则采用零填充不完全LU分解预处理。为了验证这种线性化方法,上述方法被用于:① NACA 0012翼型的非定常绕流线性系统构造与求解;② NACA 0012翼型稳态流动的伴随方程构造与求解;③ AGARD wing 445.6机翼颤振问题降阶建模。上述三个算例的结果与CFD模拟的吻合一致。     

新概念轮翼气动特性实验研究

研究一种关于小型飞行器的新型升力和推力系统--轮翼的气动特性.通过自行设计并搭建的基于虚拟仪器的实验平台,进行地面试验,试验主要研究了轮翼转速对系统气动特性的影响.结果表明:随着轮翼转速的增加,系统提供的升力及推力增大.由于轮翼系统消耗的能量绝大部分用来产生升力及推力,因此轮翼系统具有较高的工作效率,为进一步研制小型无人直升机探索新型动力模式.