刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

空间飞行体远尾区内场和密度的非稳态非线性计算

对空间飞行体远尾区内等离子体与场之间非稳态非线性相互作用的耦合方程，在考虑了时间导数项后，进行了数值模拟计算．计算中考虑到场量具有轴对称的特点，利用二维三分量的初始条件，采用FTCS有限差分格式方法，得到了场和密度扰动的分布．结果表明，这种非线性相互作用也会引起场的塌缩并出现密度空洞，但塌缩的速度较快．通过对密度空洞变化的探测，可以对隐身飞行体进行跟踪探测．     

粗差拟准检定法在星载GPS低轨卫星定轨中的应用

利用自适应卡尔曼滤波进行星载GPS低轨卫星定轨时,必须解决量测方程中经常存在的粗差问题．在分析以往方法的优缺点后,用拟准检定法来探测和修正量测方程中存在的粗差．该法的优点是辨识粗差准确率高,能同时定位多个粗差．另外,为了克服星载GPS低轨卫星定轨的滤波器可能出现的数值不稳定性及发散现象,还采用了UD分解算法及Sage自适应滤波器．最后用一个CHAMP卫星的模拟算例验证本方法的可行性和有效性．     

太阳风质子束流的演化

太阳神飞船观测表明，太阳风高速流中质子束分量相对于核分量的密度随日心距离增加而增加．提出解释这一观测现象的机理并给出二维数值模拟结果．由于阿尔芬波速随日心距离增加而减少，第二支左旋波将与更多的质子共振，把部分原来属于核分布的质子拉到束分布中来．用数值模拟方程方法求解回旋波共振导致的准线性扩散方程，数值结果与观测结果相符合．     

异型封闭腔辐射几何学问题

讨论了解决异型封闭腔辐射几何学问题的几种方法:异型封闭腔内,表面间可见性的判断方法;优化异型封闭腔角系数计算的方法;确定异型封闭腔内定向辐射的壁面辐射源的方法.计算了两种结构喷管封闭腔内的角系数,并从角系数的互换性和完整性定理验证了方法的正确性和有效性,为异型封闭腔的辐射计算解决了几何问题.     

短喷管型面设计

用特征线法计算出超音速短喷管的型并做了附面层位移厚度修正，对该型面的喷管流场求解了ＮＳ方程，证明流场指标符合要求，流场校测结果得到一致结论。     

一类高阶有理差分动力系统

讨论了高阶有理差分方程xn 1=(αxn-k 1)/(-1 xnxn-1…xn-k 1),for n=0,1,…,其中α<0,k∈{1,2,…},x-k 1,…,x0都是实数,并且使得x-k 1…x0≠1。给出了这类方程解的明确表达式,并对其解的稳定性和收敛性也做了进一步讨论。     

压电材料壳的广义变分原理和状态向量方程

将三维弹性材料的H-R变分原理引入到具有机-电耦合效应的三维压电弹性材料圆柱壳问题中,建立了对应于正交各向异性压电材料圆柱壳的Hamilton型广义变分原理。通过变分运算得到了柱坐标系下的状态向量方程,该状态向量方程也是Hamilton正则方程。最后指出纯正交各向异性弹性材料壳Hamilton型广义变分原理是变分公式的特例。变分原理的建立将有利于压电材料圆柱壳静力学和动力学问题的有限元方法或半解析法的推导。     

某型导弹尾喷流形状的数值模拟

从三维薄层近似N—S方程出发，采用EN0差分格式，对某型导弹层流场进行了数值模拟。计算模拟了飞行马赫数Mα0．8—1．6，海拔高度H＝2km—15km，燃烧室总压为20MPa和40MPa条件下导弹的层流场。计算结果表明与飞行速度相比，海拔高度和燃烧室总压对层流宽度有显著影响。通过对流场结构的分析与测量，建立了层流宽度与上述参数的拟合公式，为工程估算这一类型导弹层流影响边界提供参考。     

叶轮机械叶栅流道几何处理的计算方法

叙述了叶轮机械叶栅流道网格生成前的几何处理的计算方法,通过将处理后的数据保存为PLOT3D文件格式,可以在其他CFD软件中直接使用,或在商业网格生成软件中快速获得叶栅流道的实体造型,与传统的在CAD类软件中进行交互式几何处理相比,可以大大缩短CFD计算周期。