气泡堆积法生成曲线多边形区域非结构化网格及其应用

对气泡堆积法进行改进,发展了一种基于气泡堆积法生成曲线多边形区域非结构化网格的算法,将曲线边界映射为直线完成气泡添加和位置动态调整,再通过弧长参数化的方法将气泡位置映射回曲线边界,避免了移动调整过程中气泡偏离曲线的问题.提出了一种简便的判断点在复杂曲线多边形区域内的方法,简单且易于编程.通过添加不同大小的顶点气泡以及在内部设置人工点源,利用加权平均法实现了网格的局部加密.应用基于非结构化同位网格的SIMPLE算法对环扇形空腔顸盖驱动流进行了数值模拟,不同雷诺数下的计算结果与文献结果吻合较好.     

气泡堆积法生成局部加密非结构化网格

对生成非结构化网格的气泡堆积法(BPM)进行改进,提出了一种生成局部加密非结构化网格的算法.该算法通过向多边形顶点添加任意大小气泡以及在内部设置点源,并利用反距离加权平均法传递节点信息,简单地实现了局部加密/稀疏网格的生成.同时,还提出了气泡数量的控制方法,可对过于稠密或者疏松的气泡群进行气泡的删除和添加,使所有的气泡更好地与计算区域相匹配,消除了不规则三角形网格的生成,提高了生成网格的质量.用改进的BPM算法对矩形区域生成了边界加密、中心稀疏的非结构化网格,在雷诺数为1000的工况下对方腔项盖驱动流进行了数值模拟,与基准解对比非常吻合,证明生成的网格质量较好.     

基于气泡吸附的自由曲面三角形网格生成方法

为了满足自由曲面网格结构对网格均匀性的要求,提出了一种自由曲面的三角形网格自动生成方法——空间气泡法.首先,按一定规则在曲面上布置适量的网格点;然后,将网格点模拟为大小相同的弹性气泡,引入气泡间的相互作用力,得到各气泡的运动控制方程,并利用Velocity-Verlet算法求得气泡系统的平衡状态;接着,通过求气泡中心点的三维Voronoi图与曲面的交线得到曲面上的Voronoi图;最后,对曲面Voronoi图做对偶变化并删去冗余边,得到曲面上的Delaunay图,即所求的三角形网格.算例表明:该方法自动化程度高,适用范围广,生成的网格规整、均匀,为网格结构的设计提供了便利.     

基于Delaunay准则的三维网格自动插点算法

提出了一个新的三维网格自动插点算法 .该算法充分利用Delaunay的空圆特性来计算新点的位置并插入新点 ,维护三角化的Delaunay性质 ,使得生成的四面体网格和实体的几何边界完全保形 .该算法生成的四面体具有较好的性质和良好的密度分布 ,还采用了新的插点算子———线段插点、子面插点和四面体插点     

一种改进的Delaunay三角形化剖分方法

提出了一种基于Bowyer Watson算法的平面区域Delaunay三角化剖分的改进方法。它结合了前沿推进法的内部结点生成技术和Delaunay联点网格生成技术 ,使得每插入一点所破坏的单元尽可能地少。采用适当的数据结构 ,使Delaunay搜索过程限于局部 ,算法大为简化 ,易于编程 ,浮点计算量少 ,同时也避免了使用函数递归调用。采用在基网格上定义网格步长的办法控制网格的疏密 ,使网格疏密易于控制。几个算例表明 ,该算法是行之有效的。     

一种改进的Delaunay三角形化剖分方法

提出了一种基于Bowyer-Watson算法的平面区域Delaunay三角化剖分的改进方法。它结合了前沿推进法的内部结点生成技术和Delaunay联点网格生成技术,使得每播入一点所破坏的单元尽可能地少。采用适当的数据结构,使Delaunay搜索过程限于局部,算法大为简化,易于编程,浮点计算量少,同时也避免了使用函数递归调用。采用在基网格上定义网格步长的办法控制网格的疏密,使网格疏密易于控制。几个算例表明,该算法是行之有效的。     

激波诱导变形液滴外流场数值模拟

基于Delaunay非结构网格生成方法和局部阵面推进重新生成网格的方法,应用边界边剖分与合并技术,对二维含自由表面流动中运动边界网格的生成方法进行了研究,改善了局部网格重新生成过程中的网格质量;对激波诱导变形水液滴在高速气流作用下的受力情况做了数值模拟.结果表明:该文在处理二维自由表面的动边界问题中采用的方法是可行的.     

基于四叉树切削网格的N-S方程求解方法

以四叉树非结构化网格为基础,提出了背景区域采用正方形四叉树网格、边界区域采用切削网格的一种可以表达复杂几何形状的网格生成方法,该网格具有生成过程简单,正交性好等优点.在这种网格的基础上,采用非结构网格有限体积法进行离散得到了多种形状切削网格并存时Navier-Stokes（N-S）方程的求解算法,并以顶盖驱动斜方腔流和方腔内热圆柱自然对流为例,应用上述算法实现了网格生成和流动数值模拟,与基准解进行了比较,一致性较好.计算结果表明这种切削网格方法及其N-S方程求解方法具有可靠性和应用前景.     

基于Delaunay三角形的非结构化有限元GPR正演

介绍传统Bowyer-Watson三角网逐点插入法的原理与实现步骤,并将固定边界限制、Laplacian光顺、边压缩、边分裂、点插入等拓扑变换技术应用于网格剖分的优化;为了使数值解的误差在全域内接近于均匀分布,通过间隔函数法实现点源、线源等网格渐变控制,结合局部粗化或细化技术,建立高质量Delaunay三角形网格,实现自适应网格剖分。通过1个起伏地表与断层模型网格剖分实例验证非结构化网格对于物性参数分布复杂或几何特征不规则的地电模型的适应性。根据GPR有限元波动方程,应用三角形剖分、线性插值的Galerkin有限单元法进行求解。建立1个复杂GPR地电模型,利用Delaunay三角形对该GPR地电模型进行自适应网格剖分。研究结果表明:非结构化网格对于物性参数分布复杂或几何特征不规则的地电模型都具有良好的适应性;非结构化三角形网格剖分质量好,单元密度易控制,易于实现自适应有限元,能提高复杂模型正演精度。     

一种求解N-S方程的自适应直角网格方法

提出了一种用直角网格表达背景、切削网格表达边界的非结构化自适应直角网格方法.该方法采用四叉树保存网格数据,将切削简化成6种类型,用速度的旋度和散度作为自适应加密标准,从而可实现任意二维区域网格的自动生成和自适应加密.通过将极小网格边界化处理,利用SIM-PLE算法处理速度和压力的耦合,实现了该网格上N-S方程的离散和求解.算例表明,该方法网格生成简单,可以用于任意形状上的流动和传热模拟,相比非自适应方法,用一半的网格数目即可达到相同的计算精度.