三维凸域中候选卫星点集的确定算法

局部网格生成是无缝有限元并行计算的关键问题,然而,不恰当的局部网格生成算法会产生不一致网格.借助Delaunay三角剖分的唯一性,指出了三维凸域中Delaunay网格局部生成的不一致性,分析其产生的原因,建立避免局部网格不一致性的关键算法--搜索候选卫星点集的优化探索球算法,为网格并行生成做好数据准备.并从理论上证明了任意节点的优化探索球域包含其所有卫星点,从而,保证全域网格没有不一致性.通过算例分析了算法的正确性和局部性.     

基于图的DTAD的简易稳健自动生成算法

二维任意域(含任意多个空腔)的Delaunay三角剖分(简记为DTAD)广泛应用于有限元网格生成等方面。本文提出了一种基于图的DTAD的简易稳健自动生成算法。该算法首先构造任意域中所有边界点集的约束最小生成树,然后将同时满足本文所给三个简易约束算法的边逐次引入初始三角形网格表,形成边界点集的三角剖分,再通过稳健局部优化算法和贴体生成核插入算法自动生成DTAD。本文同时分析了退化及数值计算误差对剖分结果的影响,在提高局部优化的稳健性方面对该算法做了进一步完善,使之能更好地满足任意域网格剖分的要求。文末给出具体应用实例以说明本文所提算法的有效性。     

任意平面域的三角形网格和混合网格生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网格和混合网格自动生成的方法。该方法包含了易于生成与调控的矩形背景网格技术,简便的外形输入和初始阵面划分技术,加速查寻的数据结构,以及节点松弛和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量三角形网格和适于粘性计算的三角形／四边形混合网格,具有较好的通用性与实用性。     

任意平面域的三角形网络和混合网络生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网络和混合网络自动生成的方法。该方法包括了易于生成与调控的矩形背景网络技术,简便的外形输入和初始陈面划分技术,加速查寻的数据结构,以及节点松驰和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量用角形网格和适于粘性计算的三角形／四边形混合网络,具有较好的通用性与实用性。     

一种基于约束框架的棱柱网格生成方法

提出了一种基于约束框架的棱柱网格生成方法,即利用物面网格节点与背景结构框架网格之间的单一映射关系线性插值生成棱柱网格。背景结构框架网格基于TFI方法生成。物面网格节点与背景框架网格之间的映射关系通过一种线性投影方法建立。然后通过V型槽、圆球、M6机翼和F6翼身组合体四种不同外形对该棱柱网格生成方法进行了测试,测试结果表明该方法能够避免棱柱网格生成中的网格交叉现象,保证棱柱网格的成功生成。该方法在棱柱网格质量控制、棱柱网格局部修改等方面具有较好的效率和优势,但背景框架网格的生成依然是该方法的一个瓶颈,需要进一步的研究。     

复合材料板壳任意位置处载荷识别

提出了一种受载单元判别法,根据遗传算法生成的染色体所对应的载荷作用点的横、纵坐标和每个单元四条边间的几何关系,判断出集中载荷作用的单元,再根据力的等效分配原则,把载荷等效分配到所在单元的节点上,进而基于有限元分析和遗传算法自身优越的全局搜索能力实现作用于复合材料结构网格节点或单元内部任意位置处的载荷识别.与现有的神经网络识别法和有限元反分析法相比,该方法具有能够识别非网格节点处的载荷的优点,且仿真算例表明其识别精度更高.      

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

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

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

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

三维非结构网格的生成及优化

叙述了一种新的生成三维非结构网格的Delaunay方法,该方法能够高效地生成优质四面体网格。由于采用了局部剖分的方法来保证边界的完整性,因而不需要添加新的节点到原始表面三角形中。该方法在四面体外接球的球心处插入新的节点,在插入全部内点后,采用删除边、边/面交换、网格光顺等方法对生成的网格进行了优化,基本能删除质量差的四面体。     

面向高性能数值模拟的并行网格生成软件研制

网格是开展大多数数值模拟计算的必要条件。面向高性能数值模拟应用,几何结构日趋复杂,局部特征日趋精细,网格规模将达到十亿规模以上,这就要求网格生成既要高效处理复杂模型,又要自适应刻画局部特征,更要考虑大规模网格向解算器的低开销传输,因此研制并行网格生成软件对推动高性能数值模拟应用具有重要价值。首先从算法研究、软件产品化两个层次论述相关领域的国内外情况,接着介绍了自主可控软件SuperMesh的架构设计和关键技术进展,最后探讨了并行网格生成技术的未来发展方向。     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

基于局部无网格的混合算法研究

研究了一种求解任意外形绕流的局部无网格混合算法.该算法计算区域整体采用了直角网格,只在物面附近嵌入局部无网格处理.整体直角网格的规则性与正交性,不仅使得网格生成快,而且计算格式实施简单,避免了与非正交性相关的计算项,从而引入误差相对较小.局部无网格处理只要求物面附近布点离散,具有灵活性,使得发展的算法适合处理任意外形.先以理想不可压流动问题为例,给出了其控制方程Laplace方程的具体求解实施过程,并给出了直槽道内柱体绕流的计算结果;接着把方法推广用于求解可压缩流动的Euler方程,对喷管等典型流动问题进行了数值模拟,并与传统有限体积法结果进行了比较,在捕捉的波系结构、表面压力系数分布等方面具有良好的一致性,表明发展的新方法是可行的.     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

A SELF－ADAPTIVE FINITE ELEMENT METHOD FOR SOLVING 2－D EULER EQUATIONS ON THE UNSTRUCTURED MESHES

ＡＳＥＬＦ－ＡＤＡＰＴＩＶＥＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧ２－ＤＥＵＬＥＲＥＱＵＡＴＩＯＮＳＯＮＴＨＥＵＮＳＴＲＵＣＴＵＲＥＤＭＥＳＨＥＳＺｈｏｕＣｈｕｎｈｕａ；ＹａｎｇＺｕｏｓｈｅｎｇ（ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙ...     

网络连接机群上CFD计算的一种负载平衡方法

讨论在网络连接的分布式存储并行机群上进行CFD计算的一种负载平衡方法和相关问题。当数值计算的方法是区域分解法时,各节点机上的负载与其处理的子区域的网格节点数成正比,所以同样大小的子区域可保证负载平衡。采用负载再分配即区域再剖分的方法来达到动态负载平衡,构造了一种对多块网格(每块内为结构网格但各块间的连接关系可以是非结构的)进行近似平均分划的方法,并用一种数据结构和算法保证了剖分出的多个子区之间正确的互联和边界信息传递关系。给出了用这种剖分算法对常见多块网格进行分划的例子,并对三维流场进行了实际计算。结果证明本文方法是正确的,能有效地均衡计算负载,达到较理想的并行效率。     

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。