基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

基于直接配点法的月球软着陆轨道快速优化

介绍了直接配点法在月球探测器软着陆轨道快速优化问题中的应用.首先给出了软着陆最优化控制问题模型,其中状态方程为量纲为1的三自由度模型,性能指标选为燃料消耗最小,控制变量则为推力攻角和推力.终端状态受到速度和高度的约束.然后,应用直接配点法将最优控制问题离散化为非线性规划问题,即将动态优化问题转化为静态参数最优化问题.选取各节点和配点上的状态量和控制量作为优化参数.最后应用基于Matlab语言的SNOPT软件包对参数最优化问题进行求解,该软件包对于求解大型非线性规划问题具有很好的收敛性.仿真结果表明直接配点法对于月球探测器软着陆轨道初始状态量和控制量的取值不敏感,且求解过程具有一定的实时性.因此,直接配点法对于再入轨迹优化问题的求解是可行的.     

小卫星临近作业轨道和姿态联合控制

建立了小卫星相对轨道和姿态的误差动力学模型,根据作业任务和姿态指向要求确定了小卫星相对轨道和姿态的期望运动;采用相对轨道和姿态联合控制策略,并考虑小卫星作业过程中质量和转动惯量的不确定性,设计了自适应全状态反馈控制律;数值仿真结果表明该控制律具有较强的自适应性。     

基于修iE Rodrigues参数和UKF的姿态估计算法

以修正Rodrigues参数作为姿态参数,利用无味Kalman滤波(uKF)设计了一种全姿态估计算法.修正Rodrigues参数在姿态估计中具有简洁高效的优点,然而,因其奇异性不能用于全姿态运动情况.利用修正Rodrigues参数与其影子参数相互切换的方法可以避免奇异现象的发生.UKF直接利用非线性模型而不需要线性化,适用于高精度的姿态估计.本文针对"陀螺 矢量观测"这种典型的姿态确定模式,以修正Rodrigues参数与其影子参数相互切换的方法描述姿态,利用UKF设计了姿态估计器.状态协方差阵在姿态参数切换时发生突变,引起姿态估计器的滤波波动.为此,提出了在姿态参数切换时对姿态估计器进行修正的方法.最后通过仿真验证了算法的有效性.     

一种基于Markov模型的星载计算机系统可靠性评测的新方式

可靠性在星载计算机系统中具有非常重要的作用。现以Markov链为基础,提出了一种新的估计星载计算机系统可靠性的统计使用检测方法。首先通过对故障的传播进行分析,对检测过程中Markov链上由故障产生的状态转移描述方法进行了改进,将“故障不一定产生失效”的现象用Markov链描述出来。其次,提出了程序执行路径的概率计算方法,并且用该方法对星载计算机可靠性进行了估计。最后通过仿真实验表明了该计算方法的有效性。通过对故障状态转移的分析和路径概率计算方法的提出,使星载计算机系统可靠性评测更能反映实际情况。     

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

拉普拉斯反变换的算法和程序

拉氏反变换可直接用查变换表法或传递函数的部分分式展开法求得.在用部分分式展开法时,往往很难确定部分分式展开式的系数值.本文根据C.F.Olsen和C.G.Brockus提出的一种含有多重实极点或复极点函数的拉氏反变换的新算法,结合劈因子求根法,用PASCAL和BASIC两种语言编制了程序(分别在DUAL 68000和PDP-11/23计算机上调试通过),这样就能很好地解决这个困难,完成拉氏反变换. 本文只介绍用PASCAL语言编制程序所用的框图,并对程序中用到的两个公式作简单证明.在程序运行中发现了一些问题,提出了解决的工程方法.     

逐次逼近摄动法及其在力学中的应用

提出逐次逼近摄动法，将预估的零级近似解逐渐加以修正，使之逐渐逼近于其真实解。有非线性问题的算例。算例表明，解法的精度和其它奇异摄动解法（如L—P法）的基本相同，但方法允许进行一次选代，使其精度更加提高。按照方法，可求得任一级近似解，其精度逐级提高，计算过程也很简便。     

典型激波针减阻降热特性及流动机理

为研究激波针对超/高超声速钝头飞行器进行减阻降热的相关特性,采用数值模拟方法对6种典型激波针构型的流动特征开展系统研究。给出了激波针长度、来流马赫数对流动特征的影响规律,并对其形成机理进行了探讨。结果显示:马赫数较低时,头部有扰流物的5种激波针在回流区即将分裂时减阻率最大;马赫数较高时,减阻率在回流区分裂前后出现局部峰值,但最大减阻率将出现在回流区分裂后更长的激波针长度下。马赫数为3时,6种构型减阻率达最大时的相对长度在0.8~1.2之间,相比而言,球型、半球型和双锥型的减阻效果最好,最大减阻率为45%~50%;圆锥型最差,为20%~25%,明显低于头部有扰流物的构型。相同的激波针长度下,头部有扰流物构型的减阻率随马赫数增大而增大,圆锥型则相反。流场回流区、分离激波、弓形激波、局部膨胀流动等导致的压力分布变化是构型整体阻力变化的主要成因。