复杂外形的动态混合网格生成方法

本文将层推进方法、四分树方法和阵面推进法相结合,建立了二维复杂外形的动态混合网格生成方法.首先用层推进方法在运动物体附近生成贴体的四边形网格,然后用四分树方法在外场生成自适应的矩形网格,而中间由三角形网格过渡.三角形网格由阵面推进法生成,并由节点松弛和"对角线交换"方法进行优化.当物体运动时,贴体的四边形网格随物体运动而运动,而外场的矩形网格保持静止.中间过渡层的三角形网格随物体运动而变形.当物体运动位移较大导致出现质量较差的单元甚至导致网格相交时,则在局部重新生成网格.利用上述动态混合网格生成方法,生成了多种含相对运动的复杂外形的动态混合网格,网格质量令人满意.     

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

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

生成非结构网络的改进阵面推进法

通过对非结构化网格生成的阵面推进方法进行分析,对该算法进行了优化设计,采用笛卡尔直角坐标网格来代替背景非结构化网格,其上网格尺度参数通过控制体有限元法求解Possion方程的Dirchlet问题得到,提高了网格生成的效率和质量。同时利用Laplacian滤波算子和对角交换方法来改善三角形网格的质量,尤其是对角交换技术对消除畸形三角形单元效果明显。对固壁边界采用三次样条拟合后重新离散的方式形成初始阵面,改进了网格的贴体性,提高了网格的光滑性。在编制通用程序时利用链表结构和派生类型的数据变量,节省了计算内存,提高了网格生成效率。     

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

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

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

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

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.     

Delaunay非结构网格生成之布点技术

Ｄｅｌａｕｎａｙ 非结构网格生成技术是 Ｃ Ｆ Ｄ 领域广为应用的有效技术之一。如何布置流场中的计算点是 Ｄｅｌａｕｎａｙ 三角剖分的一个关键问题。本文以二维非结构网格生成为例， 介绍三种布点策略， 即重心布点、外接圆心布点和 Ｖｏｒｏｎｏｉ 边布点， 算例表明三种方法均可生成满意的非结构网格， 其中 Ｖｏｒｏｎｏｉ 边布点方法生成的网格质量更好一些。三种方法均可方便地推广到三维问题。     

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

二维欧拉方程在非结构网格上的自适应有限元算法

发展了一种易于和自适应网格加密技术相结合的非结构三角形网格生成技术。这种方法的特点是域内节点可任意规定;根据欧拉方程解的误差估计,在流场某些区域增加节点后,能重新快速生成新的具有高质量三角形单元的非结构网格。欧拉方程采用Ｇａｌｅｒｋｉｎ空间离散和Ｒｉｃｈｔｍｙｅｒ两步显式时间推进相结合的有限元法求解,并以沿流线熵的变化量作为物理判据进行网格自适应加密,提高解的精度。     

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

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

基于离散曲率计算的三角网格模型优化调整

采用面积夹角加权的三角网格模型顶点法矢及三角片质心权值对Taubin的三角网格模型离散曲率计算方法进行了改进,在此基础上提出了一种新的三角网格模型优化调整方法。用该方法调整三角网格模型,在模型上曲率变化较平缓的平坦区域及曲率变化较剧烈的特征区域都能取得较好的调整效果。     

复合材料自动铺带技术研究——曲面铺带轨迹算法

为满足自动铺带成型工艺要求,在自由曲面上规划预浸带轨迹。采用三角形网格划分逼近自由曲面方法,将商业CAD软件导出曲面的表面三角化数据格式（STL）文件用于曲面分析,对网格信息重构,建立点、边、面之间的关联信息。算法中应用“自然路径”原理在曲面上规划预浸带轨迹,研究了曲面边界处轨迹计算问题,提出了过三角形顶点出现“奇异性”问题的解决方法。通过在曲面上求解预浸带轨迹,证明了该算法的可行性。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

飞机蒙皮机器人化制造曲面分片算法研究

为实现飞机曲面机器人加工的自动轨迹规划,研究了基于模型几何特征的飞机曲面分片算法。通过分析飞机构件加工的工艺参数信息,以表面加工允许的最大曲率或最大偏角为阈值,提出基于近似曲率的三角网格模型分片方法。以三角网格中任意面片为起点,沿网格邻接关系向外拓展联接,形成若干近似平面化的模型分片结构,当超过阈值时则停止拓展,从而得到一整片曲面。通过对飞机CAD模型进行分片,发现采用本文提出方法可以保证分片结果的均匀性,为后期加工轨迹自动规划规划提供保障。     

基于激波拟合与柔性壁的反设计数值方法研究

将激波和待求激波压缩面作为运动间断边界拟合,研究了超声速来流下的型面反设计数值计算方法.算法秉承拟合法思想,采用二阶格式时间推进欧拉方程组,在运动边界和动网格条件下,按预估的波系结构划分网格块,反求可产生指定压力分布的三维型面,并获得精确的激波型面.通过包括凸包进气道在内的三个反设计算例验证了本文算法的合理性及其对网格质量的适应性,证明了在型面反设计中,采用激波拟合法可大大提高型面精度.本算法可用于超声速进气道压缩面和唇口完全三维反设计.     

动网格技术在增升装置优化中的应用

将基于Denaulay背景图的动网格方法应用于增升装置二维缝道参数优化设计过程中。针对二维增升装置的构型特点和翼型运动的特点,采用了添加辅助点的双层Denaulay背景图动网格方法,并详细论述了辅助点添加、运动技术。为了进一步提高动态网格质量,提出了针对前后缘运动位置中比较极端的两种位置生成两个初始网格,并采用两种初始网格进行动态网格的生成。结果表明新动态网格方法可以为二维增升装置提供高质量的动态网格。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。     

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

用隐式多重网格法计算实用外形无粘绕流

网格展弦比的好坏严重地影响多重网格的效率。文中细致地处理了多重网格算法的细节,了网格展弦比的影响,发展了一种计算实用外形绕流的隐式多重网格解法,计算效率提高了５倍。