COMPUTATION OF TURBULENCE FLOW ON HYBRID GRIDS USING k-ω TURBULENCE MODEL AND OSHER SCHEME

An unstructured Reynolds-averaged NavierStokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equatlon turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.     

关于薄膜反射镜成形机理的研究

聚酰亚胺薄膜反射镜是近几年发展起来的新概念。以Karman方程为基础讨论了薄膜反射镜的成型原理，分析了外载荷对薄膜反射镜面型的影响；给出了非线性耦合Karman方程的求解方法；以一种聚酰亚胺圆薄膜为例，用有限元法对Karman方程进行了数值求解，用zemike多项式拟合方法得到了圆薄膜在外载荷作用下的近似解析解，证明了把Karman方程解析解表示为多项式的合理性；用有限元法得到的不同外载荷作用下的薄膜的变形分析证明了通过控制外载荷的加载方式和载荷强度可以控制薄膜反射镜的面形，改善反射镜像差。     

短喷管型面设计

用特征线法计算出超音速短喷管的型并做了附面层位移厚度修正，对该型面的喷管流场求解了ＮＳ方程，证明流场指标符合要求，流场校测结果得到一致结论。     

压电材料壳的广义变分原理和状态向量方程

将三维弹性材料的H-R变分原理引入到具有机-电耦合效应的三维压电弹性材料圆柱壳问题中,建立了对应于正交各向异性压电材料圆柱壳的Hamilton型广义变分原理。通过变分运算得到了柱坐标系下的状态向量方程,该状态向量方程也是Hamilton正则方程。最后指出纯正交各向异性弹性材料壳Hamilton型广义变分原理是变分公式的特例。变分原理的建立将有利于压电材料圆柱壳静力学和动力学问题的有限元方法或半解析法的推导。     

基于N-S方程的翼型实验侧壁边界层控制

使用三维可压、雷诺平均N-S方程数值计算了NACA0012翼型在二元翼型风洞中不同迎角和马赫数状态下的压力分布,改进了划分计算网格的Hilgenstock方法,发展了一种利用N S方程计算结果,来确定翼型实验中为消除侧壁边界层干扰的抽气压力的方法。通过比较利用本文介绍的方法计算的NACA0012翼型实验侧壁抽气压力和相同状态下,在西北工业大学0.3m×0.1m二维跨音速风洞中经油流显示证明为合适的抽气压力,结果表明该方法是可行的。     

基于变域变分有限元的翼型反设计

反问题作为一种自由(未知)边界问题近年来得到了广泛的研究。变域变分通过把自由边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而自由边界求解和流场分析可以完全耦合进行。本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型，计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果。文中还采用了两种计算网格(H型和O型)来设计翼型，比较了它们的精度，计算表明本方法完全能用于翼型的反设计。     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

高速风洞支架干扰数值修正研究

提出了一种计算高速风洞支架系统对飞行器模型纵向气动力干扰量的数值计算方法 ,从跨声速全位势积分方程出发 ,编制了适用于飞行器全机模型及其带支架情况下的跨声速绕流计算程序。通过对双垂尾模型和GBM 0 3模型两个算例的计算 ,讨论了尾支撑位置及其几何外形参数对模型气动力的影响 ,并对GBM 0 3模型带短支杆情况下的纵向实验结果进行了修正。表明该方法对于分析研究风洞模型支架干扰问题并进行支架干扰修正是可行的、有效的 ,可以作为选择尾支撑位置及其几何外形参数和对跨声速风洞纵向实验结果进行支架干扰修正的工具。     

考虑测量约束的非合作目标近距离快速导引控制方法

研究了考虑测量约束的非合作目标近距离快速导引控制方法。首先,基于C-W方程建立了解析形式的径切联合近距离快速导引控制模型,并采用微分修正方法将解析结果进行修正;然后推导了任意初始条件下,径切联合控制双脉冲与视场角的解析关系;最后基于相对偏心率矢量建立了径切联合控制的控后效果评估模型。仿真结果表明:提出的控制方法和控后效果评估模型有效。