非结构网格生成技术

描述了非结构网格生成的过程,网格的生成方法基于Delaunay 三角形方法,并对物面附近网格进行拉伸以便获得在边界层内高拉伸比的三角形单元。网格的生成过程分两部分:(1) 按Delaunay 三角形方法生成适于无粘情况下的非结构网格。(2) 对近物面的部分三角形单元进行拉伸,使生成的网格能适用于粘流。使用上述方法本文构造了二维弯管内绕叶栅区域的非结构网格     

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

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

一种基于块结构的非结构网格快速生成方法

以阵面推进法为基础,提出了一种用预制模块的方法实现在平面域内生成非结构网格的新方法,即模块方法。模块是由一组三角形构成,文中以三角形模块为例,介绍了网格模块的预制方法。网格的生成是根据空间步长的背景信息,从预存的网格模块中选取当地合适的模块,直接复制或变形到计算区域实现的。由于生成过程中每一步引入的是预制的网格模块而不是传统的单个三角形单元,因此可期望新方法要比传统方法快得多。算例显示了新方法的快速网格生成过程,该方法具有推广应用前景。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

基于切割单元直角坐标网格的DSMC方法及优化技术

对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。     

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

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

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

二维非结构网格DSMC方法的实现及其应用

研究了非结构网格DSMC方法实现的过程，将Bird位置元方案中的子网格思想引入到非结构网格上来，只存储子网格的总体标识号，做到了利用较少的计算网格提高分子的分辨率和计算精度；提出了将面积元坐标搜索算法与交替数字二叉树搜索算法(ADT)相结合来跟踪模拟分子在网格之间的迁移；使用ADT方法判别分子与物面是否相交，节约了计算时间；最后对过渡流域高超声速绕流进行了数值试验．结果表明这种方法的可行性、高效性及算法的通用性。     

一种新的非结构动网格的生成方法

弹簧方法是一种常被采用的非结构动网格生成方法,文中采用一种新的非结构动网格生成方法即弹性体方法(Elastic solid method,ESM)。和弹簧方法相比,虽然弹性体方法的计算量大,但是可以处理大位移的情况,尤其善于处理网格扭转很厉害的情况。另外,采用弹性体方法可以减小网格重新生成的频率,从而可以减小映射误差,同时也可以减小网格重新生成带来的计算量。将两种方法的数值模拟结果进行对比,进一步表明该弹性体方法在处理动网格问题中的优势。     

二维非结构网格生成及其数值模拟

提出一种二维非结构网格生成算法以及一种新的绕翼型的数值计算格式。首先,在Delaunay生成算法的框架下结合了阵面推进算法的优点,从覆盖整个计算区域的三角形着手,给出判断准则,由这些三角形的边直接生成新点,让新点与符合要求的三角形的边的两点相连构成三角形,这样生成的三角形网格具有较好的几何性质;其次,针对二维N0012翼型跨声速绕流问题提出了一种新的计算格式。通过建立插值模板的概念,在每一个模板上运用最小二乘和MUSCL思想相结合的办法针对初始变量构造线性插值多项式,改进了原有的权系数的计算准则,将各插值多项式通过加权平均得到在整个模板集上通用的插值多项式,从而构造出高阶的有限体积格式;最后,计算结果表明本文工作具有一定实际应用价值。     

三维复杂外形非结构网格的生成及应用

以 Delaunay方法为基础生成了三维复杂外形的四面体非结构网格。这种方法的特点是先布结点 ,然后按 Delaunay准则将其连接起来。与其他方法相比 ,其好处是能较好地根据计算需要控制网格点的疏密 ,提高计算效率 ,其运行过程也较稳定。在布点过程中采用了背景网格 ,使点的生成自动而合理。文末给出了 M6机翼和简化航天飞机的数值算例。     

基因方法在网格结点位置优化中的应用

将基因方法应用于网格结点位置的优化中。文中首先简单介绍了基因优化方法中基于达尔文进化论和Ｍｅｎｄｅｌ基因理论的基本原理,其中包括插索空间表达、三个基因作用器（选择、交配和变异）等要点;然后着重阐述了相关偏微分方程的离散误差和三角形网格几何形状的适应度函数的定义、结点位置的二进制基因表达及基因方法的优化进程。离散误差是在二次非连续彭鼓包（ｂｕｍｐ）函数的空间中近似定义的,并且在点移动过程中相关解的二     

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

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

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

等边三角形网格加筋壳轴压承载研究

随着新一代中型运载火箭研制的开展,贮箱焊接筒段的精细化设计成为结构减重、火箭提高运载能力的关键。对比分析国内外主要液体火箭贮箱筒段的不同网格加筋形式,发现三角形网格对轴压承载有着显著优势。随即对三角形网格加筋金属壳体的轴压承载能力展开算法分析、有限元计算和全尺寸级破坏试验的研究,最终验证等刚度铺层计算方法正确可行。理论分析和试验结果均表明,相比传统的斜置正交网格形式,三角形网格拥有更高的结构效率。按照结构最优原则设计的等重量三角形网格圆筒的轴压承载能力高出斜置正交网格15.3%左右。采用三角形网格加筋可有效减轻贮箱结构的重量。     

复合材料弹翼－弹身组合结构强度计算

本文利用自行开发的复合材料结构应力分析软件（CMSSA）对某型导弹的复合材料弹冀－弹身组合结构进行了静力分析。本文的计算软件采用动态链表数据结构存贮总体刚度矩阵，具有很高的存贮和计算效率。文中给出了任意四边形壳单元和三角形壳单元的形函数及几何矩阵，这两种单元能适应结构中的复杂曲线边界及单元过渡区域。     

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

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

用于弹性旋翼流场模拟的网格变形方法研究

建立了一个基于改进"Ball-Vertex"的具有拓扑无关性的网格变形方法,用来精确模拟桨叶发生弹性变形后贴体网格的运动以及旋翼近场气动力环境,为更好地分析弹性变形对桨叶贴体网格的影响提供了一套解决方案.该网格变形方法从基本的弹簧模拟思路出发,用一网状弹簧系统来模拟桨叶贴体网格的变形和运动.在网格结点周围引入"Ball-Vertex"概念,通过添加相关的冗余约束,该方法保证网格结点在变形过程中不会越过与该点所对应的对角面,从而避免了畸形网格单元的产生,可有效地用于较大变形情况下的流场数值计算.使用所建立的方法,首先对绕二维翼型的结构网格和非结构网格进行了网格变形模拟,以表明方法的有效性;然后,将该方法应用于振荡翼型的流场数值求解,进一步验证了所建立的方法;最后,结合该网格变形方法和旋翼流场求解方法对弹性旋翼流场进行了数值计算.在此基础上,得出了一些有意义的结论.     

基于扫掠法的六面体网格生成算法及实现

提出了一种基于扫掠法的六面体网格生成算法,该算法通过源和目标曲面参数域间的近似最小二乘仿射映射将源曲面的网格投影到目标曲面,并利用背景网格插值的方法计算内结点,从而将2.5维实体生成六面体网格。该算法避免了空间网格投影过程中耗时的“寻根”过程,并且不需要导致效率降低的网格优化,能够处理各种2.5维实体,生成的网格质量较好。