外弹道数据处理中航向角精确计算方法

针对国内靶场在外测事后数据处理中航向角计算存在的问题，深入分析了传统计算方法存在的缺陷，提出了两种新的航向角计算方法——坐标法和速度法，并提供了两种方法的算法和模型。在理论分析的基础上，采用实测数据对新旧三种方法比较分析，最后验证速度法是目标瞬时航向角最佳计算方法。     

模拟真空羽流场的特征线法

介绍了特征线法（MOC法）对真空羽流场进行数值模拟的过程，用MOC法进行了喷管出口流动不均匀情况真空羽流场算例的求解，所得结果与理论分析相符合，与DSMC法计算结果吻合，提高了计算效率，节省了计算时间。结果表明：用MOC法可代替DSMC法对连续介质区进行模拟，并可在用N-S和DSMC耦合数值模拟中确定耦合边界。     

振动信号的幅值和相位计算方法分析

给出了基于Kaiser窗FFT频谱法、正弦逼近法和全相位时移相位差法计算振动信号幅值和相位的计算方法,并使用三种方法分别对单一的正弦信号以及混有噪声和其它频率的正弦信号进行了测试。通过测试数据对三种方法的信号处理效果进行了分析比较,阐述了不同计算方法的优缺点。比较结果显示具有抑制频谱泄漏和＂相位不变性＂的全相位谱分析结合高精度相位差法而得的全相位时移相位差法精度更高,应用范围更广。     

有限体积法的分析与改进

本文对以Jameson为代表提出的常用的一种有限体积法进行了较深入的分析，指出了它的适用范围，网格对计算精度的影响并提出了改进措施。算例表明本文提出的分析与改进措施是正确、有效的。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

基于相对流面理论的液体火箭发动机泵内部流动计算

应用相对流面理论对某液体火箭发动机离心泵内叶轮机械部件的内部流动进行计算和分析.流面的求解采用流线曲率法.根据不同部件内的流动数值模拟编制了模块化计算程序,使各部件的计算相对独立.不同部件之间采用混合平面法进行处理,将非定常问题转化为定常问题,简化计算.黏性的影响通过几种损失系数进行修正.数值计算的结果显示诱导轮和离心叶轮增加了液体的压强,而导流通道主要改变了液体的流动方向,这与理论分析一致      

燃烧室温度计算方法分析及应用

介绍了航空发动机燃烧室出口温度常用的两种计算方法，结合批产某型发动机试车数据进行计算，对两种计算结果进行对比分析。用涡轮流量法对某燃气发生器的燃烧室出口温度进行计算，其结果与外方给的数值基本一致。提出用涡轮流量法计算燃烧室出口温度较切合实际，为发动机寿命评估及其工作可靠性提供依据。     

用Simons法计算真空羽流

姿态控制发动机真空环境下工作产生的羽流对航天器造成的影响非常明显，用Simons法对真空羽流流场进行了解析计算，包括等熵核心区和边界层膨胀区的计算，得到了符合流动规律的结果，这是一种快速计算真空羽流场物理参数分布的方法，工作对卫星、飞船和空间站上的姿态控制发动机的工程设计和科学实验具有参考价值。     

化学非平衡高超声速电离流动数值研究

本文对非平衡化学区应高超声速绕流零攻角轴对称钝头体电离空气粘性流动进行了研究，通过带有化学源项的Ｎａｖｉｃｒ－Ｓｔｏｋｅｓ方程对流动作数值模拟，总的连续方程由单个组元守恒方程代替。高温电离空气中存在７个主要的组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ，和ｅ￣－。化学动力学模型由７个组元和这些组元之间６个可能的反应来描述。发展了一种数值方法以求解这组完全耦合的与时间有关的偏微分方程组。利用一类新的有限差分方法来求解这组以守恒形式表示的控制方程。这是一种以ＴＶＤ格式来离散粘性流通量的方法。在数值实验中采用了ＴＶＤＹｅｅ－Ｈａｒｔｅｎ格式，对全部算例都是用激波捕捉法。关于化学组元，我们试图处理完全催化壁、非催化壁和有限催化壁等三种不同的表面边界层条件。计算是对ＲＡＭＣ－Ⅱ试验条件进行的，得出了７１ｋｍ高空情形的典型流场数值解。对表面压力，Ｓｔａｎｔｏｎ数和热流率给出了化学反应气体结果与完全气体结果的比较。计算的流场电子密度与ＲＡＭＣ－Ⅱ实验数据进行了比较，计算结果与实验数据符合是相当好的。     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

模态跟踪技术在颤振计算中的应用

以p-k法为基础，在特征值选取过程中，采用了多点预测方法及优化选择方法进行模态跟踪，有效解决了模态之间的窜支问题。介绍了基于模态跟踪技术的防窜支颤振计算方法；编制了防窜支p-k法颤振计算程序；对某工程实例，分别使用该计算程序与MSC．Nastran软件进行计算，经结果对比分析，证明该方法能有效解决窜支问题。     

燃烧室流量分配计算方法适用性研究

为考察有效面积法、流阻法在燃烧室流量分配设计中的适用性,分别采用2种方法完成了某环形燃烧室流量分配计算,并采用数值模拟的方法进行了对比分析,以试验手段验证了其可靠性。结果表明:有效面积法及流阻法理论计算简化条件基本可以得到满足,计算精度相当。采用有效面积法得到主燃孔进气量为32.54%,低于实际值3.12%,得到掺混孔进气量为20.43%,高于实际值0.3%;采用流阻法得到主燃孔进气量为35.11%,低于实际值0.55%,得到掺混孔进气量为17.32%,低于实际值2.81%。2种方法大孔总进气量均低于实际值约3%。数值模拟结果与试验结果吻合较好,可用于燃烧室流量分配的精确计算。     

“Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析

针对高超声速进气道前缘"Ⅳ型"激波干扰产生的气动加热与结构传热多物理场耦合计算问题，发展了一种基于有限体积法的流-热-固一体化计算方法。该方法采用一体化控制方程组统一离散求解外部高速流场与内部结构温度场，规避了传统分区耦合算法在时间域内交替迭代的繁琐数据交换策略。另外，提出一种新的双温阻模型计算流-固交界面的物性参数以保证计算准确性，采用LU-SGS隐式时间迭代和自适应时间步长以提高计算效率。采用经典高超声速二维圆管流-热-固耦合算例对该一体化方法进行验证，计算结果与试验值和参考文献数据吻合较好，证明了该方法的可靠性和正确性。利用一体化方法对高超声速前缘"Ⅳ型"激波干扰流-热-固耦合问题进行定常/非定常计算与分析，给出了温度与热流的时变特性，计算结果表明，激波干扰作用产生的超声速"喷流"不断冲击壁面，使得壁面最大压力系数增大约9倍，壁面最大热流增大约4.7倍，给高速飞行器的热防护设计与选材带来严峻挑战。同时，也表明了一体化计算方法可以较好地用于长航时飞行条件下与复杂飞行环境下的高超声速热防护系统的热环境特性分析与综合性能评估。     

用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

序列响应面法及其在飞机结构可靠性分析中的应用

响应面法是当前求解大型结构可靠性问题最有希望的方法之一。本文阐述了响应面法的基本思想，提出了一种与重要取样法相结合的序列响应面法，并介绍了基于这种方法的飞机结构可靠性分析程序RSM。文中给出了用RSM对一系列结构问题（包括著名的AGARD机翼模型）进行可靠性分析的结果，并通过与精确解及随机有限元计算结果的对比，指出了本文所述方法与程序的实用性。     

超声速机翼—尾翼对突风动态响应的数值方法

本文以Ｇｒｅｅｎ函数为基础的边界元法，结合时间历程法，计算机翼－尾翼组合体对突风动态响应的气动性能。本文介绍的方法，特别适于研究、分析各类飞行器非定常动态响应，为飞行器的设计及其非定常气动特性分析提供了一个行之有效的计算方法。编制的计算程序快捷、精度好。本文提供了单个机翼对突风动态响应的几个算例，有的与Ｌｏｍａｘ解进行了比较，吻合得很好，说明本文所使用的方法是可行的，程序是成功的。对于机翼－尾翼组     

液体火箭发动机实验台液路系统工作过程仿真

从一维可压缩液体的瞬变管流方程出发,采用有限体积法,对液体火箭发动机实验台液路系统工作过程进行了一维数值仿真,并将计算结果与实验结果进行了比较与分析。结果表明,该方法较好的描述了发动机实验台液路系统的工作过程,能为实验台系统的设计、故障诊断、优化以及实验发动机的性能分析提供指导。      

D φ—ψ流线坐标系中求解跨音速绕流问题

本文提出了在一种新的流线坐标系，位函数一流函数坐系系中求解跨音速流问题的数值方法。由于跨音速流动对边界条件的变化比较敏感，现在数值计算跨音速流动已大量采用边界拟合坐标系，流线会标系也是一种边界拟合坐标系，且由于流线随着计算的进而而自动调整。因此，本方法与近来发展较快的自适应网格技术有许多相似之处。数值计算表明，本方法在稳定性和收敛性方面都比较好，与实验结果及其它数值计算结果相比，相互吻合程度令人满     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。