结合部网格系统在机身-边条-机翼组合体绕流计算中的应用

对应用于嵌套网格技术中的Collar网格思想进行了拓宽，提出了建立结合部网格系统的思想，使得对于复杂的组合体，Collar网格的生成都非常容易，并将其与虚拟网格和其他网格相结合，在保证计算网格的生成方便快捷而且网格质量高的前提下，成功地解决了嵌套网格技术中的结合部问题。机身－边条－机翼组合体绕流计算的实例表明，将结合部网格系统应用在嵌套网格技术中能保证几何外形不发生变化，能有效地处理各种复杂组合体外形的结合部问题，对于嵌套网格的自动化生成具有一定的实际意义。     

网格生成软件GRGES的设计与实现

首先,本文说明了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。其次,本文论述了数值网格生成的基本概念和基本步骤,提出了网格生成软件GRGES的设计思想,简单介绍了GRGES的总体结构和五种网格生成方法。最后给出了三个网格生成算例,它们均可用于计算实际流场。     

无限插值方法生成三维网格

 本文研究了生成三维网格的无限插值理论,应用它构筑了绕翼身组合体的三维C-O型网格,这是一种代数方法,适用于不同拓扑结构的网格生成,具有程序编写简单、计算时间短的优点。在这种方法生成的网格上进行了二维与三维的数值实验,结果表明,它是一种有效的网格生成方法。     

非定常数值模拟方法的发展及其在动态绕流中的应用

基于混合通量分裂的思想,通过应用Gauss-Seidel迭代求解差分方程,构造了一种时空二阶精度、无条件稳定的隐式迭代NND算法,并讨论了时间精度与稳定性,亚迭代收敛判则,几何守恒律的应用以及动壁边界条件等相关问题.通过引入一种简便易行的加权函数来综合刚性动网格生成技术和超限插值生成动网格这两种方法的优点,发展了一种计算量小、比较实用的加权动网格生成技术.作为应用实例,本文给出了多个动态物体绕流的数值模拟算例,计算结果表明了本文数值方法的成功.     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

Delaunay非结构网格生成之布点技术

Ｄｅｌａｕｎａｙ 非结构网格生成技术是 Ｃ Ｆ Ｄ 领域广为应用的有效技术之一。如何布置流场中的计算点是 Ｄｅｌａｕｎａｙ 三角剖分的一个关键问题。本文以二维非结构网格生成为例， 介绍三种布点策略， 即重心布点、外接圆心布点和 Ｖｏｒｏｎｏｉ 边布点， 算例表明三种方法均可生成满意的非结构网格， 其中 Ｖｏｒｏｎｏｉ 边布点方法生成的网格质量更好一些。三种方法均可方便地推广到三维问题。     

叶片机叶片平截面型线坐标的快速插值法

建立了一种快速而精确的代数算法,用于插值计算柱坐标系下叶片机叶片平截面空间点坐标,并给出了相应算例。适用于计算流体动力学（ＣＦＤ）工作者构造各种复杂的计算网络,并可以推广到各种复杂的空间曲线点的插值计算。     

发动机进口流场等值线极坐标绘图法

整个绘图系统分三层。第一层包括(MA-INLA, MAINLB等子程序)的圆上或环上插值曲面计算;第二层包括(PLOTWA, MAPLIN子程序)的极坐标下矩形网格(即扇形网格)棱边与等值线的交点计算;第三层把等值线与棱边的交点绘成平面上折线和光滑曲线。一、二层是平行关系,第三层是由第二层调用。具体结构见图。     

PVM平台下通量守恒的嵌入自适应网格解激波和底部流

运用嵌入式自适应网格技术来提高结构网格求解复杂流场的分辨率。该方法根据流场的形成,将密集网格嵌入到物理量剧烈变化的区域,使求解精度得以提高的同时,避免较大的计算开销,从而达到求解效率和求解精度之间的平衡。 Ｐ Ｖ Ｍ 平台使得被嵌入的网格可以在网络中不同的机器上完成,大幅度降低了墙上时间。原始网格与被嵌入网格之间数值通量守恒的实现增强了该方法在捕捉间断时的有效性和可靠性。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

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

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

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

A fast and automatic full-potential finite volume solver on Cartesian grids for unconventional configurations

To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algo rithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual (GMRES) algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution con verging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body (BWB).The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

含孤岛型腔铣削加工的螺旋刀轨生成算法

当型腔含有孤岛时,使用传统环切铣削方法加工往往难以避免拐角方向存在突变的问题,故提出一种螺旋刀轨生成算法,实现含凸边界孤岛型腔的快速铣削加工。首先在孤岛和型腔的加工轮廓上设置放样点,将放样点连接成放样线,线性插值放样线生成螺旋折线;接着插入以指数函数规律分布的控制顶点,并以其所定义的任意阶B样条曲线规划最终刀轨。所生成的刀轨能适应孤岛偏心的情况。实验结果显示,与传统环切相比,高阶连续的B样条曲线刀轨可直接用于具有NURBS插补功能的加工中心,孤岛轮廓的加工精度高,可有效避免刀轨方向突变,降低切削力的变化幅度,提高加工效率。