利用微气囊扰动进行流动控制的数值模拟研究

探讨一种可以用来模拟“三角翼上分布微气囊，从而控制流动获得滚转力矩”的数值方法。为微型飞机的气动设计提供一种工具。研究的内容包括：考虑微气囊的三角翼网格生成、流场NS方程计算、微气囊不同布局对流动的扰动等。     

基于给定压力分布的叶栅流场分析

介绍了采用时间相关的有限体积法，对不稳定的Euler方程求解；对已知叶型，在给定叶片表面压力分布的条件下计算其流场，对解法上的困难进行了仔细的分析，并在此基础上提出了一种处理叶片表面边界条件的新方法。     

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.     

超声速自由剪切尾流动的数值模拟和理论分析

本文从二维完全Ｎ－Ｓ方程出发，采用ＮＮＤ２Ｍ差分格式，计算了超声速自由剪切层流场，并且研究了上游扰动在剪切层内的传播规律，在数值模拟的基础上对流动涡结构进行了分析。     

超声速前机身及翼—身组合体的全速势方程数值计算

通过数值求解全速热方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某些一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

一类非线性演化方程族的可积耦合

寻求孤立于理论中的可积系统是可积系理论中的一个重要研究课题，其一般模式是现行的屠规彰格式，但寻求耦合系统的可积性问题，也仅是由马文秀博士研究的一个孤立子方程的可积耦合，对于一族孤立子方程的可积耦合目前尚未研究。本文通过构造一个新的Loop代数^～G和作一个恰当的Lax对变换，得到了一类非线性演化方程族的可积耦合。作为例子说明，本文求得了一类AKNS（4个人名字：Ablowitz Kaup,Newell Segur）族的可积耦合。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

尾支杆干扰非线性修正方法的初步研究

利用CFD方法 ,采用两套网格重叠的形式 ,研究了一种对尾支杆干扰进行修正的非线性方法。运用本方法 ,研究了尾支杆等直段长度、直径和锥角对模型气动特性的影响 ,研究了四个模型有攻角时的尾支杆干扰 ,并对GBM 0 4A模型的零阻进行了尾支杆干扰的修正。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

贫氧推进剂一次燃烧产物的热力计算

贫氧推进剂一次燃烧的物理和化学动力学因素决定了其一次燃烧为非平衡燃烧。说明了非平衡燃烧的特点，提出了计算贫氧推进剂一次燃烧的理论模型。以参与化学反应的推进剂组分为基础，把燃烧剩余物和燃烧产物分开，建立了燃烧产物的质量和化学平衡方程以及体系的能量平衡方程，采用布林克莱法解决了贫氧推进剂非平衡燃烧产物的热力计算。