一种有限区域运动目标分布模型及计算

潜艇的运动模型对于反潜效能的评估具有重要意义。针对应召搜索时概知潜艇的初始位置、但潜艇运动不确定性的特点,建立了一种有限区域内潜艇的运动和分布模型,并利用数值积分的方法从理论上进行了计算分析,给出了计算的结果。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

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

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

迟滞非线性二元机翼颤振特性分析

 采用多项式迟滞非线性模型建立二元机翼气动弹性运动方程，并用数值积分法进行求解。通过系统响应振幅随来流速度变化的分叉图和频谱分析发现，俯仰方向由于含有非线性因素，振动中的高阶分量随速度提高不断增加，并引起高次分叉。重点研究“机翼/空气”质量比以及“沉浮/俯仰”两个自由度的自然振动频率比对非线性颤振速度边界的影响，并提出可以通过提高自然振动频率比来减小迟滞非线性因素的不利影响。     

用多重网格方法计算旋翼跨声速无粘流场

发展了一种加快悬停旋翼无粘流场计算收敛速度的多重网格方法。由于悬停旋翼流场中存在不可压区域，同时旋翼尾涡系统的发展需要较长的时间，使得旋翼流场的收敛速度远低于固定翼流场，因此研究旋翼流场的多重网格算法具有重要意义。空间离散格式采用了中心有限体积方法，时间推进应用了五步龙格-库塔法。采用3层网格的V循环，对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明：尽管多重网格方法对旋翼流场的加速收敛作用不如对固定翼流场的加速收敛效果，但是多重网格方法仍然可以显著地加快旋翼流场收敛。     

Adaptive Gaussian sum squared-root cubature Kalman filter with split-merge scheme for state estimation

The paper deals with state estimation problem of nonlinear non-Gaussian discrete dynamic systems for improvement of accuracy and consistency. An efficient new algorithm called the adaptive Gaussian-sum square-root cubature Kalman filter（AGSSCKF） with a split-merge scheme is proposed. It is developed based on the squared-root extension of newly introduced cubature Kalman filter（SCKF） and is built within a Gaussian-sum framework. Based on the condition that the probability density functions of process noises and initial state are denoted by a Gaussian sum using optimization method, a bank of SCKF are used as the sub-filters to estimate state of system with the corresponding weights respectively, which is adaptively updated. The new algorithm consists of an adaptive splitting and merging procedure according to a proposed split-decision model based on the nonlinearity degree of measurement. The results of two simulation scenarios（one-dimensional state estimation and bearings-only tracking） show that the proposed filter demonstrates comparable performance to the particle filter with significantly reduced computational cost.     

北斗-罗兰C组合导航的ASF修正计算与研究

介绍讨论了修正ASF的米林顿法和CLC数值积分法,并用MATLAB编程计算;用米林顿法计算的(SF+ASF),与军标"长波地波传输信道计算方法"计算结果非常吻合;CLC数值积分法计算的结果与英国PaulWilliams的结果也很吻合;而CLC数值积分法计算的(SF+ASF)与米林顿法有较大差别,但在用以推求陆地路径ASF时,二者差别就很小。试验表明,CLC数值积分法有助于复杂地形路径的ASF修正计算。     

J2项和大气阻力摄动作用下卫星编队构型的演化分析

研究卫星编队构型在J2项摄动和大气阻力摄动作用下的演化进程。从相对运动动力学模型出发,分别加入J2项摄动和大气阻力摄动因素,采用数值积分对所获得的运动方程求解,得到参考星和环绕星在惯性坐标系中的绝对运动规律;将两星的绝对运动规律作差后,利用坐标转换矩阵,将惯性坐标系中两星的相对运动规律转换到Hill坐标系中。由于该方法没有经过任何简化,故研究摄动的影响时不存在误差。最后通过仿真分析,给出了J2项摄动和大气阻力摄动作用下卫星编队构型演化的相关结论,并提出进一步研究的方向。     

圆柱形霍尔推力器内等离子体数值模拟

利用一维磁流体动力学模型,对圆柱形霍尔推力器的放电等离子体进行了数值模拟。考虑了等离子体的电离、中和、碰撞、玻姆扩散及阳极鞘层的影响,由龙格-库塔方法得到离子速度、离子数密度、电子的温度等分布,其与实验结果有很好的一致性。经分析可知,此分布与电磁场的分布、粒子碰撞及电子阻抗等因素有关,并分析了电子温度分布、电子速度与电磁场的关系。结果表明,离子数密度沿通道方向增加,但在出口附近略有下降;而中性粒子数密度逐渐降低;离子速度在出口达到最大值,电子速度在下游有较大的梯度分布。     

时域有限体积法的格式稳定性分析

以模型方程为例,采用MUSCL插值的S-W通量分裂格式,分别结合二步、四步龙格一库塔时间推进格式建立数值求解方法.用Von Neumann傅里叶级数分析方法,通过对不同格式的放大因子及色散因子的计算,分析了时域有限体积法(FVTD)格式的稳定性.