一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

双流混合边界的正交贴体网格生成

通过数值实验表明，用泊松方程变换形式生成的贴体坐标网格，其边界点的β值与源项呈连续单调的函数关系。根据这种关系，可以用方程迭代求根方法中的增值寻根法和对分法，使网格边界正交化。由薄壁结构相隔的两种流体混合边界正交网格生成的实例表明，这种方法对于生成适应于红外抑制器全流场计算的网格是很有效的。     

用能量优化法选取细分网格顶点

分析了混合曲面初始细分网格的特点，将网格线进行分类。在这一基础上．利用曲面光顺的网格能量法，给出了计算网格顶点的优化模型。针对优化模型的特征，优先计算网格中的关键点，并把优化模型转化为线性方程组求解。这个计算网格顶点的方法实现了初始细分网格中顶点选取的自动化。实验表明，用文中方法选出的顶点能使混合曲面具有良好的光顺性，且有较高的计算效率。     

三维Euler方程的自适应八叉树结构直角网格算法

给出了三维自适应八叉树结构直角网格的生成方法，采用了以几何外形为基础的网格加细方法，并解决了网格生成过程中的若干难点，使网格生成更简便迅速，产生的网格更适合于Ｅｕｌｅｒ方程计算。在计算中采用以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法和以面为基础的通量计算方法，减少了工作量。在物面边界处采用通量分解的方法实现边界条件，简单易行。本文对前机身、单个及多个外挂等问题进行了数值实验。结果表明，计算与实验符合较好。     

粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

自动生成有限元四边形网格的铺筑法

对一种新的网格生成技术Paving(铺筑 )法 ,进行了分析和研究 ,铺筑法用全部都是四边形的网格生成任意的二维几何图形 ,文中研究了算法的理论依据和技术要点 ,并评估了算法的使用价值。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。     

三角网格模型的数据分块算法

对测量获得的三角网格模型进行分块。首先根据三角片中心点的高斯曲率和平均曲率的符号将三角片标示成8种类型。然后任取一个三角片作为生长点，沿着边界向外扩展“生长”，将具有相同类型的三角片组成一个连通的块。最后采用两种方法对块实施优化：根据平均曲率差值最小的原则，将小块融合到邻近块中．用两点间的最短路径来光滑边界。本方法中的三角网格和块的拓扑均可以是任意的，分块结果符合零件的特征构造。     

用任意空间曲线控制的复杂组合曲面数控加工复法研究

提出了一种基于参数曲面三角化的复杂曲面刀具轨迹生成方法。即利用给定的任意空间曲线作为控制边界,构造新的参数曲面,在此简单参数曲面上生成刀具轨迹。     

计入桨叶运动的旋翼CFD网格设计技术

提出了一个考虑桨叶实际的挥舞和变距运动的旋翼CFD网格生成方法。该方法首先采用改进的几何法在桨叶表面生成初始网格，通过求解Poisson方程得出贴体正交网格，然后结合嵌套网格和动态网格的思想，给出了一种运动嵌套网格来计入桨叶的挥舞和变距运动．并针对嵌套网格中关键的贡献单元搜寻问题，提出了一种新的“伪贡献单元搜寻法(PSDE)”。该搜寻方法不仅可以适用于在结构网格中搜索贡献单元，而且不作修改即可应用于非结构网格，还可克服“InverseMap”法在结构网格中需要单独处理的奇异线(面)问题。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

Body-Fitted Momentum Source Method for Predicting Rotor Aerodynamics Characteristics in Hover

To gain high efficiency for the simulation of the aerodynamic characteristics of the rotor in hover,body-fitted momentum source(BFMS)method is proposed. In this method,the actual blade geometry is represented by the single layer of volume grid surrounding the blade. Aiming at correctly simulating the aerodynamic characteristics of the discrete cells along the chordwise of blade airfoil section,a new distributed force model is proposed. For comparison,the RANS method with S-A turbulence model and the steady rotor momentum source(SRMS)method based on embedded grid systems are established,respectively. And the grid connecting methodology is improved to embed the blade into the background grids for the three methods. Then,simulations are performed for the hovering Caradonna-Tung rotor by these methods,and the calculated results are compared with the available experimental data. Moreover,the pressure distributions along the blade are compared with the conventional momentum source methods. It is demonstrated that the BFMS method can be employed as an effective approach to predict rotor aerodynamic characteristics with a low computational resource and reasonable accuracy.     

壁面冷却的三股混合气流的数值计算

对双层壁扩压器与波瓣喷管组合的三股气流的引射混合进行了流场温度场的数值计算。计算域划分为各子域 ,通过求解拉普拉斯方程 ,分别生成三维贴体曲线坐标网格 ,边界网格加密且波瓣处正交 ,各子域组合成整个网格。为避免因波瓣造成的网格强烈的非正交而引起解的发散 ,本文采用了同位网格 SIMPLEC计算方法和Chen-Kim修正的 k-ε湍流模型。固体区域采用大粘性解决流固耦合 ,还结合了线性欠松弛和拟瞬态欠松弛。结果证明 :双层壁间有外界冷气流被引射进入 ,形成壁面的冷却气流 ,相对单层壁扩压器 ,双层壁扩压器的壁面温度明显降低。对于小扩压角 ( 0～ 5°) ,随着扩压角的增大 ,壁面平均温度降低 ;但当扩压角较大时 ( 5～ 1 0°) ,壁面平均温度则增加     

GENERATION AND APPLICATION OF UNSTRUCTURED ADAPTIVE MESHES WITH MOVING BOUNDARIES

INTRODUCTIONMany engineering applications require thenumerical prediction of strongly unsteady flowsinteracting with moving bodies.A case in pointis that the detailed calculation of sphere motionsin an unsteady field induced by blast waves is akey issue to determine the drag force exerting up-on the sphere by the surrounding gas and,inturn,to reveal the mechanism of unsteady effectson the drag coefficient.To cope with such pro-blems,one should simulate unsteady flows withmoving boundary co…     

Numerical Simulation of Aerodynamic Characteristics of Helicopter Rotor with Tip Slots

Effects of tip slots on the aerodynamic characteristics of helicopter rotor were investigated numerically by solving three-dimensional Navier-Stokes equations based on unstructured overset grids algorithm.Improved delayed detached eddy simulation (IDDES) based on the Spalart-Allmaras turbulence model and adaptive grid refinement technique were employed.Several slots in the rotor blade tip were designed on the base of Caradonna-Tung rotor to study the effect of tip slots.Numerical results show that tip slots are able to introduce the airflow from the leading edge and turn it in the spanwise direction to be ejected out of the face at the rotor blade tip,which can reduce the strength of the rotor blade tip vortex and accelerate the dissipation process.Although tip slots may lead to the decrease of airfoils' lift coefficient at the root of the rotor blade,it can increase the lift coefficient of airfoils at the rotor blade tip,so the lift of the rotor with tip slots is almost the same as that of the rotor without tip slots.In addition,tip slots can also reduce the intensity of the tip shock wave,which is beneficial to reduce the wave drag of the rotor.     

非平行直纹面的电解磨削

为解决整体叶轮扭曲叶片型面的精加工难题,本文进行了五轴联动数控展成电解磨削的基础研究.在分析平行直纹展成电解磨削整体叶轮扭曲叶片的过切误差、得到过切误差计算公式的基础上,提出用圆锥磨轮和组合式五轴联动数控方案来消除过切误差的方法.文中还介绍了数控展成电解磨床的结构与运动、经济型多轴数控系统及其联动控制方法,并建立了电解磨削非平行直纹展成曲面的数学模型,开发了五轴联动数控展成电解磨削自动数控编程系统,对航空发动机整体叶轮的叶片型面进行电解磨削.结果表明,加工效率比手工修磨、抛光提高了12倍以上.     

基于嵌套网格和计入尾迹影响的旋翼流场数值模拟

给出了三维Euler方程数值模拟悬停状态旋翼流场的方法和模型。为充分考虑旋翼尾迹对流场的影响,采用了由动量理论导出的远场边界条件和由尾迹诱导速度代入滑移条件而构造的新翼面边界条件。为减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。针对嵌套网格中关键的贡献单元搜寻问题,本文采用一种伪贡献单元搜寻法。应用上述方法进行了算例计算,给出了有无尾迹修正时的旋翼桨叶表面压力分布和沿展向升力分布的计算结果,并与可得到的实验数据进行了对比。此外,从计算的等涡量线图上分析了旋翼尾迹流场畸变的发展趋势,表明了Euler方程法在旋翼涡尾迹捕捉方面的能力。     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格，该方法的特点是在生成单元的同时引入结点，与其他预布结点的方法相比，其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构，有效地提高了程序运行效果，最后给出了Ｍ６机翼的数值算例。     

一种简易的曲面隐藏线消去法

本文介绍一种简易的曲面隐藏线消去法。先将曲面上各网点分成可见与不可见两类,然后根据线与线相关性,判断各网格线段的可见性,最后只显示或绘制所有的可见网格线段,得到消隐后的曲面立体图形,本法比较简易,调用方便,通用性较强,不管曲面是根据什么方法构作的,只要顺序输入两个方向的网点数及所有网点的坐标值,调用本子程序,一般均能进行曲面消隐。