用逐点插入法自动生成全四边形的自适应有限元网格

本文给出一种有限元网格全自动划分的方法。方法包括三部分：１）用单元尺寸场控制生成网格的疏密分布；２）用基于Ｄｅｌａｕｎａｙ三角化过程的逐点插入法生成纯三角形网格；３）将三角形网格转化为全四边形网格。相应的网格生成器具有良好的用户界面。只须最小限度的边界描述即可自由、快速地生成单元形状良好、疏密分布任意的网格。     

用逐点插入法自动生成全四边形的自适应有限元网络

本文给出一种有限元网格全自动划分的方法。方法包括三部分；１．用单元尺寸控制生成网格的疏密分布１；２．用基于Ｄｅｌａｕｎａｙ三角化过程的逐点插入法生成纯三角形网格；３．将三角形网格转化为全四边形网格。相应的网格生成器具成良好的用户界面。只须最小限度的边界描述即可自由，快速地生成单元形状良好，疏密分布任意的网格。     

二维自适应网格生成的改进AFT与背景网格法

在基于重生成的自适应有限元网格生成算法研究中,将推进波前法（AFT）与背景网格法结合并提出改进方法,有效地解决了网格生成和单元尺寸计算这两个关键问题。改进的AFT方法,将前沿区分为活跃前沿和非活跃前沿两类,在选取目标前沿时既考虑前沿尺寸又考虑前沿分类。改进的背景网格法,利用结构化栅格对背景网格进行管理,在栅格中直接存放背景网格中的单元,既提高了新单元尺寸的计算速度,又从数值上保证了新生成网格中单元之间尺寸合理过渡。     

一种基于双前沿推进的边界层网格生成算法

高雷诺数粘性流动模拟对边界层内的网格正交性有特殊要求.对于复杂外形,这类问题的网格自动化生成十分困难.面向该问题,提出一种双前沿推进思想,并形成一种面向复杂几何外形的边界层网格全自动生成算法.结合多种网格技术处理局部几何特征以保证边界层网格的质量.双前沿推进思想同时适用于多块结构网格和混合网格的边界层网格生成.多个模型...     

前沿推进曲面四边形网格生成算法

将一类前沿推进三角形曲面网格生成算法拓展到曲面四边形网格生成。新算法逐个单元推进前沿,避免了铺路法（Paving）逐排推进方式的鲁棒性问题,针对四边形算法在理想点计算、候选点列表构建和新单元生成等环节存在的特殊技术问题给出了系统的解决方案。数值实验表明,本文算法能针对复杂的组合参数曲面自动生成全四边形网格,网格质量优于...     

适应复杂外形的三维粘性混合网格生成算法

提出一类适应复杂外形的粘性混合网格生成算法。表面网格由前沿推进三角形曲面网格程序获得,边界层布置各向异性的三棱柱体网格,远物面区域采用Delaunay方法生成四面体网格。针对模型的复杂几何特征,综合采用了各种网格处理技术,以保证边界层网格的质量,并避免算法失效问题。网格实例及计算结果表明了本文算法的实用性及和效性。     

基于人工神经网络的非结构网格尺度控制方法

网格自动化生成和自适应是制约计算流体力学发展的瓶颈问题之一, 网格生成质量、效率、灵活性、自动化程度和鲁棒性是非结构网格生成的关键问题. 在非结构网格生成中, 网格空间尺度分布控制至关重要, 直接影响网格生成质量、效率和求解精度. 采用传统的背景网格法进行空间尺度分布控制需要在背景网格上求解微分方程得到背景网格上的尺度分布, 再将网格尺度从背景网格插值到真实空间点, 过程十分繁琐且耗时. 本文从效率和自动化角度提出两种网格尺度控制方法, 首先发展了基于径向基函数(RBF)插值的网格尺度控制方法, 通过贪婪算法实现边界参考点序列的精简, 提高了RBF插值的效率. 同时, 还采用人工神经网络进行网格尺度控制, 初步引入相对壁面距离和相对网格尺度作为神经网络输入输出参数, 建立人工神经网络训练模型, 采用商业软件生成二维圆柱和二维翼型非结构三角形网格作为训练样本, 通过训练和学习建立起相对壁面距离和相对网格尺度的神经网络关系. 进一步实现了二维圆柱、不同的二维翼型的尺度预测, RBF方法和神经网络方法的效率与传统背景网格法相比提高了5~10倍, 有助于提高网格生成的效率. 最后, 将方法推广应用于各向异性混合网格尺度预测, 得到的网格质量满足要求.      

黏性边界层网格自动生成

高雷诺数黏性流动在壁面附近存在边界层,在计算模拟中自动生成可靠且有效的计算网格仍然是计算流体力学存在的瓶颈.三棱柱/四面体混合网格技术在一定程度上缓解了这个困难.然而,对于复杂外形的情况,在边界层内自动高效生成高质量的三棱柱单元仍然十分困难.常用的层推进法在凹凸区域及角点处生成的边界层网格单元质量较差,边界层网格最外层尺寸不均匀.针对这些问题,发展了一种黏性边界层网格自动生成方法,通过顶点周围边的二面角识别物面网格特征确定多生长方向,预估并调整生长高度处理相交情况.同时提出一种三维前沿尺寸调节方式,提高了边界层网格单元的正交性,保证了边界层网格与远场网格尺寸的光滑过渡.通过ONERA M6翼型以及带发动机短舱的DLR-F6翼身组合体等外形的网格生成实例及绕流数值模拟,将计算值与标准实验值进行对比,结果表明:该方法能够自动高效地生成满足数值计算需求的混合网格.     

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

AutoCAD平台上的有限元前处理软件AutoMESH

概述了ＡｕｔｏＣＡＤ平台上的有限元前处理软件ＡｕｔｏＭＥＳＨ的主要功能,着重介绍了该软件在二维和三维有限元网格自动生成自适应方面的技术特点,给出了几个网格实例,ＡｕｔｏＭＥＳＨ是工程技术仍员进行计算机辅助有限元分析的理想工具。     

用逐点插入法生成Delaunay四面体自适应网格

介绍一种基于Delaunay算法的四面体自适应网格的自动划分方法。该方法用单元尺度场控制生成网格的疏密分布,在不满足尺度场要求的单元面形心处插入新节点,同时计算新节点单元尺寸参数,实现三维实体的Delaunay四面体自动划分。此方法具有几个特点:一是表面网格与体内网格同步划分,无需区分两者;二是结点与单元同时生成;三是生成网格自适应性好,疏密分布任意。另外,还介绍了三维网格划分中两个相关算法:一个是约束面恢复算法,该算法基于约束面不允许有单元边与之相交的性质而提出的;另一个是将二维射线法推广至三维空间,判断一个点是否在一多面体内,实现了凹多面体的划分。最后通过算例对单元质量进行了评价。本文所述方法是一种有效的四面体自适应单元生成算法。     

A Voronoi‐based ALE solver for the calculation of incompressible flow on deforming unstructured meshes

An Arbitrary Lagrangian–Eulerian method for the calculation of incompressible Navier–Stokes equations in deforming geometries is described. The mesh node connectivity is defined by a Delaunay triangulation of the nodes, whereas the discretized equations are solved using finite volumes defined by the Voronoi dual of the triangulation. For prescribed boundary motion, an automatic node motion algorithm provides smooth motion of the interior nodes. Changes in the connectivity of the nodes are made through the use of local transformations to maintain the mesh as Delaunay. This allows the nodes and their associated Voronoi finite volumes to migrate through the domain in a free manner, without compromising the quality of the mesh. An MAC finite volume solver is applied on the Voronoi dual using a cell‐centred non‐staggered formulation, with cell‐face velocities being calculated by the Rhie–Chow momentum interpolation. Advective fluxes are approximated with the third‐order QUICK differencing scheme. The solver is demonstrated via its application to a driven cavity flow, and the flow about flapping aerofoil geometries. Copyright © 2010 John Wiley & Sons, Ltd.     

脆性材料破坏模拟的网格生成算法

脆性材料的破坏过程具有随机性,当前的网格生成算法没有充分考虑脆性材料破坏时裂纹扩展和碎块生成的随机性。在Persson网格生成算法与Delaunay随机网格剖分理论基础上,提出了一种可根据模拟需要动态控制网格品质的网格生成算法。通过对随机分布点的Delauna三角化,生成初始网格,然后将网格体系比拟为桁架结构,网格节点即为桁架节点。桁架节点在虚拟力作用下可动态调整位置,并最终达到整个体系受力平衡。对Persson 算法中的尺寸分布函数和收敛条件进行了修正,从而提高了收敛速度,并适用于任意形状对象的网格剖分。 基于VC++平台开发了算法程序。通过实例对算法进行了验证,表明算法能够满足脆性材料破碎模拟的需要。     

基于AFT-Delaunay的二维解耦并行网格生成算法

面向平面任意几何区域网格生成,提出了一种将波前法AFT(Advancing Front Technique)与Delaunay法相结合的解耦并行网格生成算法。算法主要思想是沿着求解几何区域惯性轴,采用扩展的AFT-Delaunay算法生成高质量三角形网格墙,递归地将几何区域动态划分成多个彼此解耦的子区域;采用OpenMP多线程并行技术,将子区域分配给多个CPU并行生成子区域网格;子区域内部的网格生成复用AFT-Delaunay算法,保证了生成网格的质量、效率和一致性要求。本算法优先生成几何边界与交界面网格,有利于提高有限元计算精度;各个子区域的网格生成彼此完全解耦,因此并行网格生成过程无需通信。该方法克服了并行交界面网格质量恶化难题,且具有良好的并行加速比,能够全自动、高效率地并行生成高质量的三角网格。     

AUTOMATIC MESH GENERATION OF 3—D GEOMETRIC MODELS

In this paper the presentation of the ball-packing method is reviewed, and a scheme to generate mesh for complex 3-D geometric models is given, which consists of 4 steps: (1) create nodes in 3-D models by ball-packing method, (2) connect nodes to generate mesh by 3-D Delaunay triangulation, (3) retrieve the boundary of the model after Delaunay triangulation, (4) improve the mesh.     

高质量点集的快速局部网格生成算法

高效及高质量的局部网格生成算法是基于节点有限元并行方法设计的关键。泡泡布点算法能够在复杂区域上不经过人工干预生成高质量的节点集,本文提出了基于该方法所生成的节点集的快速局部网格生成算法。该算法充分利用泡泡布点方法提供的节点集及节点邻接链表信息,避免了桶数据结构的建立以及节点的局部搜索过程,只需应用Delaunay三角剖分的外接圆准则从中心节点的邻接链表中去除极少数的非卫星点,可快速地生成局部网格,比现有的局部网格生成算法更为快捷。算例结果表明,该算法高效可靠,生成网格与Delaunay三角剖分网格一致。     

面向复杂几何模型的并行四面体网格生成方法

面向大规模工程计算等数值模拟领域,提出了一种支持复杂几何模型的大规模四面体网格并行生成方法。该方法以复杂几何模型作为输入,首先采用串行网格生成方法生成初始四面体网格,然后通过两级区域分解方法将初始网格分解为多个子网格并分配到相应的进程中,进程间并行地提取出子网格的表面网格,并基于几何模型对面网格进行贴体加密,最后对加密后的面网格采用Delaunay方法重新生成四面体网格,该方法可以更好地适应高性能计算机体系结构,较好地克服了并行方法中并行性能和网格质量不能兼顾的问题。对三峡大坝模型进行测试和验证,证明该方法具有良好的并行效率和可扩展性,可以在数万处理器核上并行生成数十亿高质量四面体网格。     

三维约束Delaunay三角化的边界恢复和薄元消除方法

提出一种有效的三维约束Delaunay三角剖分的边界恢复算法，该算法综合了P．L．George算法和N．P．Weatherill算法的优点，通过将约束边和约束面加以恢复，保持了实体边界的完整性，解决了经典Delaunay算法不能剖分凹域的问题，从而实现了复杂三维实体的网格剖分。提出了一种简易而有效的消除薄元方法——薄元分解法，彻底解决了三维Delaunay三角剖分过程中所产生的薄元问题。实践证明，本文提出的边界恢复算法和薄元消除算法健壮有效，生成网格的质量高，并且易于实现。     

Two‐dimensional implicit time‐dependent calculations on adaptive unstructured meshes with time evolving boundaries

An implicit multigrid‐driven algorithm for two‐dimensional incompressible laminar viscous flows has been coupled with a solution adaptation method and a mesh movement method for boundary movement. Time‐dependent calculations are performed implicitly by regarding each time step as a steady‐state problem in pseudo‐time. The method of artificial compressibility is used to solve the flow equations. The solution mesh adaptation method performs local mesh refinement using an incremental Delaunay algorithm and mesh coarsening by means of edge collapse. Mesh movement is achieved by modeling the computational domain as an elastic solid and solving the equilibrium equations for the stress field. The solution adaptation method has been validated by comparison with experimental results and other computational results for low Reynolds number flow over a shedding circular cylinder. Preliminary validation of the mesh movement method has been demonstrated by a comparison with experimental results of an oscillating airfoil and with computational results for an oscillating cylinder. Copyright © 2005 John Wiley & Sons, Ltd.