三维分区对接网格Euler方程的并行数值模拟

本文根据CFD中分区对接网格的特点，应用基于工作站集群系统的MPI网络并行环境，实现了典型双三角翼模型亚声速绕流Euler方程的并行计算。文中对计算结果的收敛情况、耗时、加速比及效率等并行性能指标进行了讨论。结果表明，多域多机与单域单机的计算结果基本吻合；2节点并行效率达到88．21％，4节点并行效率达到78．99％。     

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

基于Hamilton-Jacobi方程的编队飞行控制

 针对椭圆参考轨道的编队飞行控制系统,基于哈密尔顿 雅克比（Hamilton-Jacobi）方程,给出一种新的生成函数半解析近似解的计算方法,同时推导了编队飞行控制系统相对运动的线性状态转移矩阵,并且验证了模型的精度。在此基础上,结合一个具体的编队飞行任务,采用线性规划算法得到了队形保持的闭环控制器,该算法得到的是脉冲形式的控制序列,因此更具有实际意义。同时还提出了一种简单的避免生成函数奇异点的方法,从而减少了生成函数间相互转换所带来的计算量。最后的非线性仿真结果验证了控制方法的有效性。     

三维舵面绕流的N-S方程计算研究

基于粘性非结构混合网格,采用格心格式的有限体积法,应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟。操纵面铰链力矩的计算结果与风洞试验结果的对比分析表明,计算方法对粘性流动具有良好的适应性。对舵面铰链力矩随迎角、马赫数、雷诺数变化的规律也进行了分析研究。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

考虑飞行器动力系统进排气效应的设计参数灵敏度分析研究

面向飞行器内外流一体化设计,基于自主研发的大规模并行化结构化网格RANS求解器以及离散伴随方程求解器,开展了考虑推进系统动力状态下进排气边界条件的变分研究,通过链式求导法则避免直接对守恒变量变分,进一步引入中间变量大幅度简化了进排气边界条件变分的难度,建立了考虑进排气效应的设计变量灵敏度高效分析方法,并通过TPS标准模型计算验证了进排气数值模拟精度,与有限差分对比验证了灵敏度计算精度,以翼上发动机气动布局进排气影响数值模拟为例,系统分析了低速、高速、定攻角、定升力状态,推进系统有无动力工况灵敏度的变化以及影响机理。     

超燃冲压发动机燃烧室的准一维计算与分析

在考虑有限速率化学反应的准一维Euler方程基础上，通过增加截面面积变化、壁面摩擦和添质的源项，发展了适用于超燃燃烧室性能分析的准一维计算方法。依次以中国空气动力研究与发展中心（CARDC）和日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为验证算例，分别采用传统的一步反应模型和发展的有限速率反应模型，模拟了燃烧室流场，并基于NAL燃烧室，计算分析了不同当量比和进口压强对燃烧室流动特性的影响。结果表明：两种方法都能得到与实验数据吻合良好的结果；和一步反应模型相比，有限速率反应模型不仅可以更细致地捕捉流场细节，而且能够初步分析化学非平衡效应的影响；对于NAL燃烧室，当量比≥0.6时，压强随当量比的升高而增大，当达到1.0时，反压已推进隔离段，且推进速度随当量比增大而增加；进口压强不大于110.444kPa时，反压随进口压强增大而升高，且当反压不小于82.833kPa时，反压被隔离在等直段燃烧室入口处；过小的当量比和过大的进口压强均会导致燃烧室出口马赫数严重下降，甚至出现亚声速出流状态。     

加筋板破坏载荷计算中的蒙皮有效宽度计算方法适应性验证

飞机机身和机翼结构多采用中长加筋板,中长加筋板在轴压载荷作用下的破坏载荷一直是工程中研究的热点;目前中长加筋板破坏载荷计算多采用Johnson-Euler方程,涉及到蒙皮有效宽度计算。对国内外常用的有效宽度计算方法进行分析,说明其产生的根源,并通过机身加筋板试验结果对其适应性进行验证。结果表明：对于薄蒙皮铆接Z型机身加筋板破坏载荷计算,本文介绍的有效宽度计算方法都能使计算误差控制在 10%,可为飞机设计人员对轴压加筋板设计提供帮助。