飞行器动态稳定性参数的数值计算方法

本文提出了一种确定飞行器动态稳定性参数的数值方法。本文首先研究了气动力与飞行器运动状态参数之间的关系，在此基础上给出了稳定性参数的定义。根据这一定义并结合计算流力学和参数辨识方法，构造了计算飞行器稳定性参数的完整数值方法，数值实验表明，该方法具有通用性强和准确度较高的优点。     

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性     

航天器贮箱气液自由界面追踪数值模拟

主要讨论了航天器贮箱在轨道航行时的微重力状态下其液体推进剂在贮箱内的形态分布及控制.文中采用VOF方法,加入了表面张力的效应,追踪气液两相流的自由界面,对液面在微重力条件下的位形变化进行了数值模拟.通过比较不同重力加速度及接触角下的两相流的相图,分析了影响贮箱中液体推进剂位形变化的主要因素及对其有效的控制方法.      

三相线钉扎的大液滴蒸发实验

利用实践十号返回式科学实验卫星蒸发对流箱,开展了三相线处于钉扎状态且接触半径大于毛细长度的无水乙醇大滴在加热PTFE表面蒸发的地基科学实验.实验发现,液滴体积随时间线性递减,但钉扎大液滴蒸发过程中没有出现恒定接触角(CCA)阶段.与小液滴蒸发的恒定接触半径(CCR)阶段相同,大液滴的平均蒸发速率也与初始体积无关,表明受...     

球锥俯仰阻尼导数的数值计算

本文从非定常Ｅｕｌｅｒ方程出发，应用谐振摄动法计算了球锥体超声速、高超声速下的俯仰阻尼导数和。Ｅｕｌｅｒ方程的求解利用ＳＣＭ方法和ＭａｃＣｏｒｍａｃｋ的两步算法。文中用任意攻角下不同俯仰轴的气动导数转换公式来简化计算，用较粗网格可以得到较精确的结果，程序编制相对简单。一次计算可以得到定常绕流解及一系列钝头率、不同俯仰轴下的多种动导数结果，在ＥＬＸＳＩ６４００机上约需２０～４０分钟。所得结果与试验结果符合程度令人满意。     

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

求解双曲型守恒律的高精度,无波动样条逼近有限体积方法

本文从三次样条逼近出发，提出了求解流体力学双曲型守恒律的一种高精度、无波动的数值方法。针对样条函数的特性，本文提出了一种新的通量限制技术，使该方法在光滑区可以达到四阶精度，在流动参数的空间分布出现拐点或极值点时分别退化为二阶或一阶精度的格式。数值实验表明，该方法对流场中的激波和接触间断有很高的分辨率，优于二阶精度的ＴＶＤ格式。     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果