双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

强制偏流喷管的三维网络生成

运用泊松方程实强了制偏流喷管内的三维网络生成，方程中的源项采用了与边界形状和边界网格点分布无关的数学形式，源项内的自由参数由横向网格线与边界的交角及横向网格线在边界处的局部曲率两个约束条件来确定，实现了由边界上的网格分布直接控制内部的网格点，使生成的网格具有良好的正交性和合理性。     

用非结构网格和欧拉方程计算多段翼型绕流

用Ｄｅｌａｕｎａｙ法划分三角形构造非结构网格时，如果布点不当，会产生边界穿透，使网格生成失败，本文分析了几种布点情形，提出不容易产生边界穿透的布点方式。这一原理已用于生成绕多段翼型的非结构网格。     

固定网格中基于边界描述的CAD模型快速分析方法

固定网格技术避免了沿结构边界繁琐的网格划分过程，在复杂几何结构分析时具有很大优势。建模精度和分析精度是采用固定网格进行结构分析的两大难题。文中针对固定网格分析中的建模精度问题，提出了一种由CAD模型至分析模型的快速转换方法，首先将CAD模型采用边界描述，然后通过射线交点法来区分固定网格中单元的属性，并借助四叉树/八叉树技术实现对边界网格的细分。最后采用加权B样条有限胞元法来保证CAD模型在固定网格中的分析精度。本文建立了一种由结构CAD模型至固定网格仿真分析的设计框架，数值算例表明了该方法的有效性。     

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

固体火箭发动机燃烧室三维流场数值模拟

利用ＳＩＭＰＬＥ方法数值求解层流Ｎ－Ｓ方程，模拟了具有移动边界的固体火箭发动机燃烧室内的三维流场。网格划分采用交错网格技术，流场计算时不进行网格划分，而是直接读取网格点的值，将网格划分与流场计算分开，解决了现有计算机内存不够的问题，计算获得了满意结果。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

一个解析三维非结构网格生成方法

给出了一个三维非结构网格生成方法。提出弱区域指示函数描述三维封闭流场形状,并证明了该函数的存在不唯一性。利用该函数基于Delaunay方法依据流场尺度生成网格结构。其中,边界引力有效解决了流场边界恢复和内嵌边界问题,弹簧振子和正四面体趋进技术有效调整网格结构并提高网格质量。最后,弹簧振子和边界引力随着弱区域指示函数的演变与监测函数的指导成功应用到自适应网格生成。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

用随机选取法（RCM）计算高温平衡气体的激波管流动

本文用RCM方法计算了高温平衡气体的激波管流动。在求解Riemann问题时,激波关系式采用精确解,而对于稀疏波关系式,本文建立了等效比热比的局部理想气体近似关系式,从而使RCM方法大为简化,并具有足够精度。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

采用二阶迎风格式的 PISO 算法研究

在ＰＩＳＯ算法基础上,发展了一种新的二阶迎风格式以提高计算精度又保证算法的稳定性。新格式应用于液体火箭发动机燃烧流动过程计算,其数值稳定性很好。     

复杂几何域中不可压流动的多重网格计算

本文首次将多重网格中的全近似格式FAS用到一般曲线坐标系下交错网格布局的SIM-PLEC算法中,采用三个流动结构:渐扩平面通道、收扩管流和轴对称弯曲管流作为算例,通过与单网格在不同Re数、不同网格点数和不同几何结构上的迭代性能的对比,证明多重网格克服了曲线坐标系下单网格SIMPLEC算法的不足,大大提高收敛速度,节省CPU时间。      

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

COMPUTATION OF COMPRESSIBLE FLOW PAST SLENDER WING－BODY USING FULL－POTENTIAL EQUATION

ＣＯＭＰＵＴＡＴＩＯＮＯＦＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＰＡＳＴＳＬＥＮＤＥＲＷＩＮＧ－ＢＯＤＹＵＳＩＮＧＦＵＬＬ－ＰＯＴＥＮＴＩＡＬＥＱＵＡＴＩＯＮＨｕａｎｇＭｉｎｇｋｅ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮａｎｊｉｎｇＵｎｉｖｅ...     

用全位势方程计算细长翼-身组合体的可压缩绕流

用Ｈ－Ｏ网格、守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕大后掠细长翼－身组合体的可压缩流动。对ＡＧＡＲＤ－Ｂ翼－身组合体的计算表明,计算的自由流Ｍａｃｈ数可从亚音速直到低超音速。算出的升力与力矩特性与实验数据符合较好     

三维翼跨音速守恒型全位势流高效差分算法的应用

本文对跨音速守恒型全位势流,用Holst的高效有限差分算法,研制出分析三维翼跨音速流的实用程序。与Holst的TWING程序不同,本文用快速保角转绘生成二维正交O型贴体坐标网格,并以此为基础构造三维贴体坐标网格,因而节省网格生成的机时。本文不仅将方法用到有中等后掠的厚翼,还提供大后掠大根梢比薄翼的计算经验。对根梢比很大,前后缘后掠角相差非常大的机翼,指出本算法的困难,并分析了困难原因。     

翼-身组合体跨声速流场计算

本文把Holst的AF2因式分解格式推广到求解机翼-机身跨声速绕流的全速势方程。利用求解椭圆型偏微分方程组来数值生成贴体的计算网格。算例表明本方法运算可靠,费时较少,可以向工程实用部门推广使用。