刚性方程的二阶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条件下导弹的层流场。计算结果表明与飞行速度相比，海拔高度和燃烧室总压对层流宽度有显著影响。通过对流场结构的分析与测量，建立了层流宽度与上述参数的拟合公式，为工程估算这一类型导弹层流影响边界提供参考。     

大跨屋盖悬挑结构的风荷载分析

针对悬挑结构的上下表面同时受到风荷载作用的特点，提出风洞试验中同步测量的实验方法，分析上下表面风压相关性及其物理意义，给出悬挑结构平均风压和脉动风压的计算方法，探讨了同步测量和异步测量情况下峰值因子的合理取值。