高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

一种跨音速翼型和机翼设计方法的新进展

 讨论一种基于N-S方程、欧拉和速势方程加混合附面层修正的跨音速流动准确控制方程逆解法，及其在跨音速翼型和机翼设计中的应用。分析了两种紊流模型对分离的适应性，探讨了N-S和欧拉方法的阻力估算精度问题，列举了若干超临界和层流翼型及机翼的设计实例以及软件验证。     

跨声速翼身组合体欧拉方程计算

本文利用Ｊａｍｅｓｏｎ等提出的有限体积法，迭代求解三维定常欧拉方程，探讨了加快收敛的途径，并用这种计算了大展弦比翼身组合体跨声速定常无粘绕流，计算结果与有关实验及文献进行了比较，获得了满意结论。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

对《离散数学》中三个公式、定义及定理的研究

本文对《离散数学》中三个公式、定义及定理进行了充分的研究及透彻的分析，分别给出了这三个公式、定义及定理的简化形式及其规律，使得它们变的浅显易懂。而这些无论对老师还是对学生都是十分重要的。     

线性系统采样Z变换的一个递推公式

本文给出了计算1/(s a)~n的采样z变换的递推公式,证明了其分子多项式系数的对称性质,外对α为复数的情况进行了推广。     

双马赫反射的数值模拟与光学干涉定量测量的比较(M_s=4.62)

本文用欧拉方程计算了楔形角为40°,M_s=4.62双马赫反射的流场,计算采用了矢通量分裂和 TVD 两种格式,并考虑了γ=1.4和高温平衡气体两种情况。流场的密度场测量采用双曝光全息、M-Z、差分三种干涉技术。壁面密度分布的计算与测量结果比较表明:三种干涉方法测量结果差别不大,都处在γ=1.4和平衡气体的计算结果范围内。     

基于弹载捷联惯性导航系统精确导航的双欧拉全姿态方法

垂直发射的地空导弹在大机动飞行时,其四元数法姿态角输出存在奇异性。为解决该问题,提出了正、反欧拉角微分方程切换的双欧拉全姿态解算方法。对正、反欧拉角微分方程进行了详细推导,论述了正、反欧拉角微分方程精华区与奇异区的倒置关系。在此基础上,建立了正、反欧拉角的转换公式。根据正欧拉角的俯仰角对正、反欧拉角无缝切换,有效解决了俯仰角在±90°附近时姿态奇异性问题。仿真和半物理试验的结果表明:双欧拉全姿态算法解算的实际姿态与理论姿态保持一致,在奇异点附近的姿态角仍能正常输出,未出现突变现象,该方法得出的导航精度比四元数法得出的导航精度提高了9.81%。     

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

CE／SE并行计算研究及在Euler方程中的应用

时空守恒元和解元方法(CE/SE)是一种高精度的数值求解方法.采用改进的CE/SE方法和并行分区技术,针对非结构网格,发展了一套求解二维Euler方程的并行程序.对NACA0012翼型的亚跨声速流动和多段翼型复杂绕流进行了分区并行求解,多区并行计算的结果与实验的结果吻合较好,表明了本文研究方法的正确性和有效性.同时文中的分区方法能实现各处理器之间的负载均衡,有效地节省计算时间.