三维多块结构椭圆型贴体网格的生成

通过求解椭圆型方程生成三维贴体网格，采用分块技术生成多块结构网格体系。通过实验验证了方法对任意复杂网络区域生成三维多块结构网格能满足流场计算的要求，其质量是相当满意的。     

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法，求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

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

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

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...     

二维间断有限元混合网格混合算法研究

针对二维流场进行结构／非结构混合网格划分,在靠近物面处生成结构化网格,在结构化网格以外和远场之间以非结构化网格填充。使用间断有限元方法（DGM）离散Euler方程。在混合网格混合算法中,结构化网格区域使用结构DGM,非结构化网格区域使用非结构DGM,保证结构／非结构对接面处的通量守恒,并加入moment限制器消除激波前后伪震荡。对NACA0012和RAE2822翼型跨音速无粘流动的数值模拟结果说明了混合算法的有效性以及在捕捉间断方面的能力。     

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

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

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

二维有限元网格全自动生成方法

主要研究二维有限元网格全自动生成方法。提出一种新的节点生成方法 ,通过实例验证该方法可靠。采用三角形形状因子评定三角形网格质量 ,可减少三角化的时间 ,优化三角网格质量。     

空间曲面的非结构网格生成

介绍一种空间曲面的 Delaunay三角形网格生成技术。采用参数表示法将空间曲面变换为平面参数域,引入相关的度量矩阵后将平面参数域进行三角剖分,最后得到空间曲面的三角剖分。空间曲面的边界点作为已知点被给定,内点采用一类自动生成的方法计算,例如按边插点法,Voronoi边插点法等。该方法可以很好地控制曲面网格的质量,使用方便,可以直接用于物理问题数值模拟中的曲面网格生成,也可以用于三维非结构网格生成中提供边界面三角剖分。文中给出的算例说明了本方法的应用情况。