大气密度分层计算时的再入弹道解析解

给出了在考虑阻力系数为速度的函数的情况下,并将大气密度进行分层计算时,估算再入飞行器弹道主要参数的一组解析式。按上述方法进行了仿真计算,并将计算结果,Ｂａｒｂｅｒａ解,有关文献解同数值解进行了比较。结果表明,由于大气密度采用分层计算,更真实地反映了飞行器的再入环境、解析估算值较其它两种解更接近于数值解,完全可用于再入飞行器的再入弹道初步估计和再入性能分析。     

高超声速再入飞行器气动轴研究

根据根据力学中的力螺旋概念，论述了再入飞行器的气动轴存在的一般性，分析了再入飞行器气动轴的三种类型，提供了确定气动轴的方法，建立了计算非对称再入飞行器飞动轴的近似方法和数值方法，近似方法基于修正的牛顿流理论和推广的内伏牛顿流理论，数值方法基于Ｊ．Ｅ．Ｄａｙｗｉｔｔ教授发展的分离系数矩阵欧拉方程时间相关方法和张鲁民研究员发展的三维定常分离系数矩阵欧拉方程的空间推进方法。     

火星大气非平衡条件下驻点热流计算公式的探讨

针对目前绝大多数火星再入飞行器驻点热流公式未考虑热力学非平衡效应故而适用性及准确性存疑的问题,分析了典型的零攻角驻点热流公式Fay-Riddell公式在火星再入飞行器热力学非平衡驻点热流计算中的局限性;然后采用数值计算方法,对涵盖火星再入飞行器飞行状态的大范围计算状况开展了热力学非平衡模拟,并对模拟得到的驻点热流公式进行拟合,得到了适用于火星再入飞行器热力学非平衡条件下的零攻角驻点热流计算公式;利用Mars Pathfinder飞行数据对公式进行了验证,计算表明,对于验证的计算状态,该公式的计算误差小于10%,符合工程预估的可接受误差范围。     

非对称体超音速无粘绕流流场数值计算

 本文给出了非对称再入飞行器的超音速三维无粘流场的数值计算方法。它归结为求解三个自变量—阶定常拟线性双曲型偏微分方程组的初边值问题。采用MacCormack二步显式差分格式。 为了说明计算方法的实际能力,对弯头钝锥和带控制翼的再入体作了计算。其精确度是令人满意的。     

地球轨道再人飞行器最小气动加热外形优化设计

本文研究了地球轨道再入飞行器气动外形优化设计问题。文中选用飞行器再入总加热量为目标函数,考虑再入飞行器纵向气动稳定性能的设计要求,在给定再入重量和约定长细比的条件下,用内点罚函数方法求解满足规定纵向气动稳定性能要求的最小再入加热气动外形。计算结果表明,平头双锥体是一种很有吸引力的地球轨道再入飞行器气动外形,它既可以最大限度地降低再入总加热量和减轻再入防热层重量,又可以有效地提高纵向气动稳定性能。     

航天飞机三维最优再入轨道及气动加热过程

本文应用现代控制理论研究了航天飞行器三维最优再入轨道和与轨道参数密切相关的气动加热过程。文中选择飞行器迎角和倾斜角作为控制变量,以飞行器气动加热率和飞行过载沿轨道积分最小作为优化性能指标,按极大原理导出最优再入轨道有约束控制的非线性两点边值问题。采用了数值优化方法——共轭梯度法求解有升力飞行器的最优再入轨道及其热过程。文中以允许误差法讨论了权系数和罚函数的选取方法;对不同速度范围研究了不同的加热模型;按热平衡方程与优化轨道同步迭代的方法求得了算例数值结果。算例的数值结果与文献[13]的量值是一致的。     

高超声速再入飞行器气动特性的快速预示—局部方法的推广应用

本文将计算高超声速稀薄气流过渡领域中气动特性的局部方法,推广应用到连续介质中弹头型高超声速再入飞行器气动力特性的快速估算。由激波风洞中M_∞=9.9时,一个8°钝锥的气动力测量结果,导出这一实验条件下的领域系数,并以此来估算不同锥角、不同钝度比及不同外形弹头型再入飞行器的气动力和力矩系数,其结果与无粘数值解及实验结果作了比较,在攻角2°～14°范围内吻合得很好。局部方法可用于弹头型高超声速再入飞行器气动特性的快速预示。     

关于非对称外形再入飞行器的压力中心系数的新定义

本文给出了非对称外形再入飞行器的压力中心系数的一种新的定义方法。用它计算出的无侧滑的非对称外形再入飞行器的压力中心系数,随攻角的变化是连续的,即不再出现∞型间断。对于对称外形再入飞行器的原来的压力中心系数的定义方法,只是新定义对于对称外形再入飞行器应用的特例。     

Weissman准则在升力式再入飞行器设计中的应用

针对升力式再入飞行器的气动特性和控制方式,分析了将传统Weissman判据推广应用到升力式再入飞行器设计中面临的一些主要问题。由赫尔维茨稳定性定理推导了采用反馈控制的滚转反逆准则的计算形式,并以HL-20飞行器为对象,采用非线性仿真和Kriging建模方法,实现了对Weissman图不同稳定性敏感区域的划分。仿真分析表明,升力式再入飞行器对应的Weissman图与传统航空飞行器存在明显不同,可以将其推广应用到升力式再入飞行器设计中。     

有动力空间再入打击飞行器航迹优化

提出了有动力再入打击飞行器的航迹优化方法。基于标准的大气模型、数值的气动模型,考虑地球自转效应,建立了有动力再入打击飞行器的航迹优化模型。将再入航迹分为助推加速段和无动力滑翔段分别优化,从而得到优化航迹。算例表明,该方法得到的航程比最大升阻比航迹增加7.3%。因此,该方法用于有动力再入航迹的优化是有效的、可行的。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

跨音速机翼结构动力响应的计算研究

通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Ｅｕｌｅｒ方程,并用Ｊａｍｅｓｏｎ提出的有限体积方案,四步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进求解;结构模型采用二自由度ＰＡＰ模型。经分析,计算结果合理。     

不等熵一维定常管流数值解法

不等熵一维定常管流控制方程的数值解法,已有许多讨论,但大多局限于某些特定情况,特别是控制方程在奇点(M=1)附近的计算,还没有一个简便的方法.本文引进新变量,把方程变为非奇异的,提出了对各种不等熵流都适用的数值解法.由于引进了喉部条件,保证了数值解的唯一性.经过对模型方程的计算精度分析,证明该数值解除了在奇点附近稍有误差外,其它地方与精确解一致.本文最后给出了该方法在无喷管发动机内弹道计算中的应用.     

多重网格方法及其在高超声速进气道中的应用

随着计算流体力学的发展,多重网格在求解欧拉方程、N—S方程过程中得到了广泛应用。将这种方法应用于高超声速进气道的数值模拟中,通过具体的算例,验证了该方法的有效性。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。