隐式格式求解拟压缩性非定常不可压Navier-Stokes方程

采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟.     

基于可观测性分析的SINS/CNS降维设计

针对长航时无人机的长航时和组合导航系统固有的非线性特性,建立了SINS/CNS组合导航系统的非线性模型,该模型能够更加趋近于真实模型。模型建立之后对其进行可观测性分析,并根据可观测性分析的结果进行降维设计,只对可观测性好的状态进行状态反馈。并选择了精度更高的SCKF算法进行滤波,仿真结果表明,选择的滤波算法精度更高,降维设计之后组合导航系统即能够保证导航精度又能够大大提高实时性。     

两种湍流模型在高速旋转翼身组合弹箭中的对比研究1）

弹箭设计、弹道计算和稳定性研究都需要准确预测旋转弹箭的马格努斯力和力矩,国内针对旋转弹箭气动特性的数值模拟工作集中在旋成体上,对带翼外形进行完全时间相关的非定常研究鲜有见到;国外虽然有对带翼外形开展研究,但以验证方法为主,对湍流模型在复杂外形弹箭旋转中的研究未曾见到.采用完全时间相关的非定常N-S方程,对带翼弹箭开展计算,对比了一方程SA（Spalart-Allmaras）湍流模型和两方程k-!SST（shear-stress-transport）湍流模型对马格努斯效应产生的影响,并分析了旋转导致的边界层和涡非对称畸变,以及周向压力分布和剪切应力分布非对称畸变.结果表明：旋转引起的物面流场参数变化主要体现在弹体中后部,SA和SST湍流模型预测的全弹马格努斯特性与阿诺德工程发展中心（Arnold Engineering Development Center,AEDC）实验及陆军研究实验室（Army Research Laboratory,ARL）的计算结果一致性很好,对动导数而言两湍流模型计算精度相当.两湍流模型计算的弹体左侧流场参数差异比右侧大,分析认为正向旋转使左侧壁面速度方向与来流速度相反,相互阻碍使气流脉动效应更强.壁面附近湍流黏性系数SA结果大于SST结果,y=0截面物面压力SA结果小于SST结果、最大相差6%,摩阻系数SA结果大于SST结果、最大相差35%.SA对旋转产生的分离抑制作用强于SST.     

S形扩张管道三维流场的数值模拟

采用改进的Beam-Warirling格式离散完全Navier-Stokes方程数值模拟了矩形变截面S形扩张管道内部亚跨声速层流流动,给出了流场的细致结构,并与文献[1]的计算结果进行了比较.结果表明:流场基本特性一致,二次流则较为清晰合理.     

进气道入口形状对冲压发动机性能影响数值研究

进气道是冲压发动机设计中非常重要的环节,其入口设计会直接影响冲压发动机的总压恢复系数σ和流量系数φ等重要性能。本文采用二阶迎风隐式TVD格式,内外流分区耦合求解可压缩N-S方程,数值模拟了掺混段出口压力与来流压力比为P/P∞=4.2超临界状态下,二元方形截面“X”型布局进气道、弹体和掺混段一体化通气模型复杂流场。计算比较了两种不同二元进气道入口形状流场和冲压发动机总压恢复系数σ、流量系数φ的结果,对其产生的原因进行了分析。