边界层流动的非平行稳定性研究

采用极为有效的抛物化稳定性方程（PSE）方法研究边界民支的非平行稳定性。并利用抛物线坐标变换、有限远延拓边界条件、及流向和法向全差分的数值方法，对不同频率的二维Tollmien-Schichting波的非平行稳定尾进行了计算和分析，计算结果与Orr-Sommerfeld方程（OES）的解以及谱配置方法的PSE结果作了详细的比较，得到了满意的结果。     

跨音速小扰动方程的多重解

本文利用Engquist-Osher格式,离散跨音速小扰动方程(TSD)和边界条件,选择不同初场计算NACA0012翼型二维定常位势流场,得到了与试验值较符合的结果,同时用Murman-Cole非守恒和守恒两格式计算并进行了比较。着重用数值实验进一步探讨多重解问题,发现并总结出了几点很有参考价值的规律和现象。     

跨音多级轴流压气机对突变型扰动的动态响应

基于多级轴流压气机的逐级特性详细研究了跨音多级轴流压气机对突变型温度扰动、突变型压力－温度联合扰动的动态响应和气动稳定性。同时还给出了动态扰动下多级轴流压气机失稳识别的新的预估方法。     

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

弹道导弹制导计算中扰动引力的快速赋值

主动段扰动引力是引起弹道导弹制导方法误差的主要因素。因此,要提高导弹的制导精度,就必须能够在弹上实时计算扰动引力。但现有方法在计算快速性和存储量之间无法得到有效协调。为此,把广义延拓逼近思想引入有限元逼近方法中,将插值单元周围节点的信息也包含到单元内一点扰动引力的计算当中,建立了一种新的数学模型。对所选发射空域,在发射坐标系中进行了直角坐标划分。计算结果表明,这种方法能够更加精确地逼近弹道导弹主动段的扰动引力,在600 km×250 km×6 km的主动段飞行区域内,只需要保存60个节点数据,就能使由逼近误差导致的落点偏差小于10 m,是一般有限元逼近方法精度的4倍以上。     

大攻角非对称流动的非定常弱扰动控制

 研制了一种新的大攻角细长旋成体非对称涡的主动控制技术, 即在细长旋成体头部施加非定常弱扰动来控制头部非对称背涡。应用七孔探针测量的空间截面流场揭示了非定常控制下非对称涡变成对称涡的流态特征。测力试验研究结果表明该方法不仅能完全消除背涡的非对称性及其产生的侧向力, 并且有效控制攻角范围从30b直到80b。     

混合移频静刚度法和直接扰动迭代法

本文对混合移频动柔度法的失代式进行等价变换，从而得到一个“变种”，该变种不包含混合移频系统的动柔度式F^△（λ*），这便消除了F^△（λ*）带来的近似性，另外，对直接扰动法也进行等价变换，引出可以省时地用于许多特征向量导数计算的直接扰动迭代法。最后，为了供读者参数，本文还导出了三种混合移频动柔度法的非迭代式，其中有两种恰好是另外两种动柔度法经混合移频后的非迭代式，这实际上已将另两种动柔度法也扩展到密集根情况了。     

低纬地磁Z分量的半年变化

利用武汉、广州、泉州和琼中等4个低纬地磁站连续多年的地磁资料，计算了各月5个磁静日Z分量日均值与中午1100---1300时段平均值之差(Dz)，对每年12个Dz采用多元回归分析方法，得到各年的半年变化幅度和相位．结果表明：4个站的Dz每年都有半年变化现象；半年变化幅度与太阳活动有关，一般来说，太阳活动高年Dz半年变化幅度明显大于太阳活动低年；太阳活动本身的半年变化，对Dz半年变化幅度有显著的调制作用；Dz半年变化的相位在3—4月(或9—10月)，即极大值出现在分季；低纬地区地磁Z分量存在显著的半年变化，能够反映赤道电急流也有明显的半年变化，这再一次证明，赤道电急流幅度的半年变化，通过“喷泉效应”使得电离层，f0F2产生半年变化，其是产生，f0F2半年变化的一个主要因素．     

神舟2号大气密度探测器的探测结果(Ⅰ)日照和阴影区域热层大气密度变化

选用了神舟2号（SZ-2）大气密度探测器在2001年2—4月间的探测数据,进行日照和阴影区域热层大气密度变化的探讨.结果表明:在高度410km附近,日照和阴影区域大气密度变幅为2—3倍,变幅的大小与地磁活动程度呈负相关关系.日照面大气密度峰区位于星下点地方时1400—1500LT的纬度处,峰值大小与太阳活动程度呈正相关关系.阴影面大气密度谷区位于星下点地方时0400-0500的纬度处,同时在±10°纬度区域中还出现了阴影面峰区.