复杂外形的非结构四面体网格生成算法

本文改进一类经典的三维Delaunay网格生成算法,给出鲁棒的Bowyer-Watson增量插点内核;针对保形边界恢复和约束边界恢复难题提出完备的解决方案。结合光滑化和拓扑变换,对初始Delaunay网格进行后处理,提升了最终网格的几何质量。实验表明,本文算法鲁棒可靠,可用于复杂外形飞行器的高质量计算网格生成。     

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

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

一个解析三维非结构网格生成方法

给出了一个三维非结构网格生成方法。提出弱区域指示函数描述三维封闭流场形状,并证明了该函数的存在不唯一性。利用该函数基于Delaunay方法依据流场尺度生成网格结构。其中,边界引力有效解决了流场边界恢复和内嵌边界问题,弹簧振子和正四面体趋进技术有效调整网格结构并提高网格质量。最后,弹簧振子和边界引力随着弱区域指示函数的演变与监测函数的指导成功应用到自适应网格生成。     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

空间曲面的非结构网格生成

介绍一种空间曲面的 Delaunay三角形网格生成技术。采用参数表示法将空间曲面变换为平面参数域,引入相关的度量矩阵后将平面参数域进行三角剖分,最后得到空间曲面的三角剖分。空间曲面的边界点作为已知点被给定,内点采用一类自动生成的方法计算,例如按边插点法,Voronoi边插点法等。该方法可以很好地控制曲面网格的质量,使用方便,可以直接用于物理问题数值模拟中的曲面网格生成,也可以用于三维非结构网格生成中提供边界面三角剖分。文中给出的算例说明了本方法的应用情况。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

用代数法生成飞机翼型进气道等复杂形状物体的网格

本文论述了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。着重介绍了网格生成的多面方法,编制了多面法网格生成程序,生成了飞机翼型、进气道等复杂形状物体的网格.该方法可快速生成二维或三维边界正交网格,同时还可有效地控制区域内网线的分布和边界上网格点的分布.     

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

三维复杂域的Delaunay三角剖分

介绍一类自动生成内点的Delaunay三角剖分算法 ,其中包括重心插点法 ,外心插点法和Voronoi边插点法。该类方法适用于二维和三维复杂计算域的非结构网格生成 ,所需人工干预少 ,使用方便 ,计算效率高 ,生成的网格质量也很好。文中并给出几个算例说明本方法的应用。     

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

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

二维Delaunay网格的一个约束边恢复算法

针对二维约束Delaunay网格生成中约束边恢复问题,提出了一个基于平面扫描策略的约束边恢复算法,证明其收敛性。给出了算法的计算复杂度,说明了算法能有效减少＂相交测试＂的次数。实现算法并用复杂算例验证了所得结论。     

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

基于“各向异性”四面体网格聚合的复杂外形混合网格生成方法

复杂外形粘性流动的计算网格生成一直是阻碍CFD实际工程应用的瓶颈问题之一,虽然三棱柱/四面体混合网格的出现极大地缓解了这一问题,但是在飞行器几何外形非常复杂的情况下,在粘性边界层内生成高质量的贴体三棱柱网格仍然是非常困难的事情。本文从"各向异性"四面体网格的生成算法中得到启示,发展了一种基于"各向异性"四面体网格聚合的混合网格生成方法。通过对DLR-F6翼身组合体和F16战斗机外形的混合网格生成及数值模拟表明,该方法生成的混合网格具有网格质量好、自动化程度高等特点,对提高复杂外形飞行器关键气动参数的模拟精度具有重要意义。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

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

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

大幅运动复杂构形扑翼动态网格生成的一种新方法

 基于Delaunay图映射的动态网格生成方法无需迭代计算,效率高,稳定性好。但对大幅运动复杂构形的动边界问题,背景图容易交叉,重新生成背景图和重新定位网格点信息不仅费时而且会导致网格质量的严重下降。提出一种双重Delaunay图映射的动网格生成方法,通过在初始背景图中添加辅助点,生成虚实两种背景图和虚实两种映射关系。分别根据虚映射关系和虚背景图、实映射关系和实背景图,移动辅助点和网格点。几个复杂构形的扑翼算例表明,双重图映射方法多付出极少的内存代价即可避免背景图交叉引起的问题,提高了动网格生成的效率和质量,增强了处理大变形复杂动边界问题的强健性。