多点择优推进阵面法生成曲面三角形网格

在现有曲面非结构网格生成法的基础上,提出了一种新的曲面网格生成法——多点择优推进阵面法.它可在曲面上直接进行三角形网格划分,克服了映射法的网格变形问题,并且可以在网格生成结束后,对曲面网格直接进行Laplace格点松弛光顺.该方法使用简单,不受曲面块类型的限制,且网格质量高,可以为三维非结构网格生成提供高质量的初始阵面,并给出了若干个算例.     

三维非结构网格自动生成

采用各向异性平面非结构网格生成技术对曲面的参数平面进行三角剖分,从而得到曲面的非结构网格,作为三维非结构网格的边界网格.应用推进面法生成网格内点,增量法将生成的内点逐点插入现有网格进行网格细化,得到三维计算域的Delaunay非结构网格.讨论了非结构网格质量优化方法.给出几个算例说明方法的应用.     

三维非结构粘性网格生成方法

描述了一套适合粘性流动计算的三维非结构网格自动生成方法.在物面附近的粘性作用区域,采用推进层方法生成各向异性的"扁平"四面体网格,并通过一定的网格伸长控制参数,实现整个流场区域网格高度的平滑过渡.当粘性网格的推进高度达到预定要求时,推进层方法自动停止,转而采用阵面推进方法生成常规意义的尽量接近正四面体的各向同性网格.同时给出了利用该方法生成的M6机翼非结构粘性网格来求解机翼粘性绕流的简单算例.     

基于Steiner点扰动和矢量边界推进技术的三维约束非结构网格生成方法

提出一种结合Steiner点扰动和矢量边界推进三角化技术的三维约束非结构四面体网格生成方法.在由物面节点生成的Delaunay网格基础上,利用Conforming方法恢复计算域形状,然后利用点扰动技术迫使所有Steiner点全部从约束边、面上转移,并利用矢量边界推进三角化方法重构约束面三角剖分,进而得到完全恢复所有约束条件的Constrained网格,并证明方法的收敛性和稳定性.     

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.     

非结构网格的生成及新型数据类型的应用

本文基于 Ｄｅｌａｕｎａｙ三角化方法及局部重组的阵面推进法，给出了任意形状平面区域的非结构化网格剖分方法。文中对加点方法及数据结构和存取模式进行了改进，同时应用链表结构、动态数组和派生类型的数据变量，对算法进行优化，节省了计算内存，提高了网格的生成效率。本文研究为特殊计算区域的三维非结构化网格的生成提供了基本方法。     

二维非结构网格生成及Euler方程计算的方法研究

用离翼型表面最小作为阵面推进法中的参数选择依据。生成二维问题的非结构网格。这种方法戏了传统的背景网格观念,直接提供网格生成过程中所需的背景信息。在求解Ｅｌｕｅｒ方程时,用格心格式的有限体积法作空间离散,用四步Ｒｕｎｇｅ－Ｋｕｔｔａ作时间推进,采用不同的加速收敛措施获得定常流动。提出了两这界条件的构造办法,并 不同边界条件对结果的影响。     

基于能带结构的布里渊区非均匀四面体网格生成技术

基于非局域赝势能带计算及四面体网格单元中能量满足线性关系,提出一种布里渊区非均匀四面体网格产生方法,在满足精度条件下能自动得出数目最少的四面体网格,使布里渊区积分计算的精度和效率大为提高．通过对硅、锗两种金刚石结构半导体简约布里渊区所产生网格的比较,表明该方法可以根据能带结构的特点自动生成优化的非均匀网格．对现有的态密度四面体计算公式进行了补充完善,并根据生成的网格和完善后的公式计算了硅和锗导带第一、二能带的态密度．     

八叉树网格到非结构混合网格的转换

在数值模拟中, 非结构网格的优势是可以采用相同的数值格式统一处理任意复杂的计算区域, 但在网格生成过程中难度大, 并且不容易控制网格质量。树结构网格可以认为是介于结构网格和非结构网格之间的一种网格, 目前已经有相对成熟的方法快速在复杂区域内生成二维四叉树网格和三维八叉树网格。在实际应用中, 数值方法往往需要在连接协调的非结构网格上做离散, 树结构网格中不同尺寸的网格之间连接不是协调的, 在应用上会受到很多限制。文章实现了树结构网格到非结构混合网格的转换, 这种转换在二维情况下就是将四叉树网格转换为非结构三角形和四边形的混合网格, 三维情况下则将八叉树网格转换为非结构混合网格。这一转换过程的难点在于需要考虑数千种不同的八叉树单元, 并给出能实现连接协调的非结构混合网格划分。可以出现的网格单元包括六面体、三棱柱、金字塔和四面体这4种不同情况。通过特别的分类, 实现了程序的自动生成, 这种程序自动生成技术一方面可以避免人工编写大量程序时的失误, 另一方面也使得对数以千计的不同情况的处理成为可能。通过对几个简单网格的测试, 对网格数据转换方法做了初步的验证。      

非结构网格自适应的新策略

对于用推进波阵面法生成二维非结构网格,提出了一种网格自适应的新方法。在网格需要加密处,布置一些小的三角形单元,并以这些单元的边作为初始推进波阵面的一部分。在推进过程中,小单元附近自然生长出尺寸逐渐变大的单元。对有激波的流动,在网格图形上,可形成十分清楚的激波图象。     

二维矢量边界推进非结构网格生成方法

提出二维矢量边界推进生成非结构三角形网格方法并证明其可行性.根据流场边界尺度布置边界节点并运用矢量边界推进方法生成背景网格,运用符号面积函数和概率筛选方法布置初始点阵,提出Spring-Laplace方法优化节点位置,同时利用边交换技术优化网格结构.该方法可包含任意点源、线源和内嵌边界,可自由进行局部自适应加密或稀疏,实现任意平面域内与尺度要求一致的高效光滑三角网格剖分.     

嵌套网格技术中的Collar网格和虚拟网格方法

针对嵌套网格技术中的结合部问题,提出了Collar网格和虚拟网格方法.用双曲型微分方程生成的处于结合部的有两个边界分别处于不同固定曲面上的Collar网格,在保证计算网格的生成方便快捷而且网格质量高的前提下,解决了为物体结合部的内外边界点提供有效插值单元的问题.虚拟网格为紧贴物面的面网格,它的作用是将物面边界条件传递给其它网格的边界面,而其本身不作流场计算.计算实践表明,将Collar网格和虚拟网格结合起来应用在嵌套网格技术中能保证几何外形不发生变化,有效地处理各种复杂组合体外形的结合部问题.     

三维多介质欧拉方法混合网格的有限体积计算格式

在采用Youngs界面重构技术的基础上，对三维欧拉方法混合网格的计算格式进行了研究。运用Youngs技术确定界面后，混合网格内每一部分物质一般不再是正规的六面体结构，可能是非规则的四面体、五面体、六面体或七面体。本文采用对非规则网格适应性很强的有限体积法对每一部分分别进行计算，给出了混合网格内每种物质的压力、速度、能量等的计算公式，比较有效地解决了混合网格的计算问题。     

结构网格与非结构网格联合技术及在内爆问题中的应用

对纯拉氏方法，格式的健壮性、精度和效率强烈地依赖于初始网格的生成。结构网格与非结构网格的联合技术，不仅可适应于复杂区域的网格生成，而且可根据实际问题的特殊性，在初始自适应的生成网格。本文在研究结构网格和非结构网格生成的基础上，研究了结构网格与非结构网格联合使用的自适应网格生成技术。首先是根据实际问题的特征，把整个计算区域划分成若干子区域，其次是把每个子区域分成若干子块，最后根据每个子块的特征生成网格，再把子块拼接起来。     

一种变分方式的平面网格生成与网格重分方法

在对流体力学方程进行数值求解时，需要应用到网格，当在物理区域上生成网格时，为了克服不同物理区域所带来的求解问题的复杂性，可以把物理区域转换为一个简单的逻辑区域，此时物理区域上的一组网格可以看作由逻辑区域内的一组点通过一个变换得到，这一变换通常要求是保边界的，同时变换的Jacobian也要求是非零的。如何选取这样的变换，使得通过变换在物理区域中生成均匀、光滑并尽可能正交的网格，是网格生成的主要目标之一。     

水斗曲面三维非正交贴体网格数值生成

为实现冲击式水轮机水斗内非定常流动的数值解析,需要高密度匀布非正交贴体坐标网格来拟合水斗曲面.本文开发了一种用于水斗曲面三维非正交贴体坐标匀布网格的数值生成程序,其主要算法是基于曲线与曲面的微分几何性质提出的四点投影法.程序可根据需要计算出不同密度的匀布贴体网格点.通过对两不同水斗的算例结果显示数值生成的网格很好地拟合了水斗不规则曲面,同时计算出各网格点的切向矢量和单位法向矢量.     

二维拉氏自适应流体动力学（LAD2D）软件研制

二维拉氏自适应流体动力学软件LAD2D，是采用建立在拉氏自适应结构和非结构网格上的有限体积格式，可以计算平面二维和柱对称二维多物质大变形弹塑性流体动力学问题。LAD2D软件系统主要由5部分组成：主控程序、数据模块、前处理模块、主体计算模块、网格模块和后处理模块。其中主体计算采用了结构网格与非结构网格联合使用的拉氏网格体系，计算格式采用了有限体积格式。网格模块包括网格生成、自适应网格加密（AMR技术）和网格重分技术，以及网格改变后物理量守恒重映技术。LAD2D软件系统由主体程序、二维网格生成程序（GRID2D）、二维自适应网格加密程序（AMR2D）、二维自适应程序（ADAPT2D）f[I-维物理量重映程序（REMAP2D）组成。     

适应强扭曲网格的扩散方程时空二阶精度计算

以全局支撑算子方法为基础,通过引入面通量,构造了具有局部模板点的时空二阶精度格式。对于大变形扭曲网格,格式采用法向修正技术和合理的单元角体积计算方法,可以保持通量的精确性。算例表明该格式在非凸网格上能够精确获得线性解; 在非光滑网格上可以达到时空二阶精度; 能够较好地保持对称性; 并适合于三维非结构网格上的求解。     

复杂外形静动态混合网格生成技术研究新进展

网格生成技术是CFD复杂工程应用的重要环节, 网格质量的好坏直接影响计算结果的精准度, 因此其已成为CFD的重要研究领域.在张涵信院士的指导下, 作者自20世纪90年代初开始开展非结构网格、混合网格技术和相应的计算方法研究, 并逐步发展至动态混合网格技术及非定常计算方法.在张涵信院士80华诞之际, 对近年来作者及团队在网格生成技术及应用方面所取得的进展进行了简要的综述, 分别介绍了静动态混合网格生成、定常/非定常计算方法、网格技术的应用等方面的进展情况.最后, 就网格生成技术目前还存在的问题, 展望了未来的发展方向.作者谨以此文表达对张涵信院士25年来的培养、关怀和帮助的崇高敬意.      

一种非结构/结构多层混合网格方法及其应用

对于粘性绕流的数值模拟,在自适应直角网格基础上,结合三角形非结构网格和结构化网格,利用其各自的优势和特点,提出一种生成混合杂交网格的思路和方法.在物面附近生成适合粘性流计算的大长宽比结构化网格,在远场分布自适应直角网格,快速离散计算空间.对于复杂的多体问题,采用三角形网格来连接各体网格,并运用网格合并的方法,保证各网格之间的光滑过渡与连接,提高网格质量.针对一些二维、三维外形的绕流问题,在上述网格基础上,采用B-L代数湍流模型和中心有限体积法,完成Navier-Stokes和Euler方程数值模拟的对比计算,结果表明网格生成和流场计算是正确的.