基于机器学习预测流场特征的网格生成技术研究进展

网格技术是影响数值模拟精度的一项重要技术。本文针对基于机器学习预测流场特征的网格生成框架,对流场特征指示器、机器学习预测流场、网格自动生成及自适应三项支撑技术进行了简要综述。现有的流场特征指示器与机器学习方法相结合有望成为提供网格生成先验参考的有效手段。在机器学习方面,适用于流体力学的物理嵌入方法是降低机器学习样本要求的可行方法。兼顾拓扑与密度分布的三维结构网格自动生成方法有待进一步研究。     

三维复杂几何外形的表面三角形网格自动生成

按点－线－面的层次数据结构,选取几种有代表性的曲线曲面,实现了三维复杂外形的表面描述。采用曲面到平面的保形变换,在平面内由阵面推进法自动生成三角形网格,实现了三维复杂外形的表面三角形网格自动生成。     

无限插值方法生成三维网格

 本文研究了生成三维网格的无限插值理论,应用它构筑了绕翼身组合体的三维C-O型网格,这是一种代数方法,适用于不同拓扑结构的网格生成,具有程序编写简单、计算时间短的优点。在这种方法生成的网格上进行了二维与三维的数值实验,结果表明,它是一种有效的网格生成方法。     

基于自由型曲面的贴体网格生成

本文建立了基于自由型曲面的曲面网格生成的数学模型,生成了曲面网格,把它作为三维Ｄｉｒｉｃｈｌｅｔ问题的边界条件,求解三维Ｐｏｉｓｓｏｎ方程,获得了满足数值模拟要求的三维贴体网格。采用上述方法,完成了基于自由型曲面的某直升机调座舱三维贴体网格的生成,结果令人满意。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

三维复杂域的Delaunay三角剖分

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

三维多块结构椭圆型贴体网格的生成

通过求解椭圆型方程生成三维贴体网格，采用分块技术生成多块结构网格体系。通过实验验证了方法对任意复杂网络区域生成三维多块结构网格能满足流场计算的要求，其质量是相当满意的。     

运输机机翼外形和吊舱位置一体化优化方法

针对运输机翼吊布局中机翼-吊舱之间有可能存在较大干扰阻力的问题,本文提出一种机翼扭转角和吊舱位置的同时优化的设计方法。气动分析模型采用高阶面元法程序PANAIR,应用参数化方法和CATIA二次开发技术实现飞机三维外形的自动生成;应用网格生成工具Gridgen的脚本语言实现气动网格的自动生成。为了减少调用气动分析模型的次数,在优化中应用了代理模型方法。算例的优化结果表明:优化后机翼的展向升力分布更加合理,显著地降低机翼和吊舱间的干扰阻力。     

复合固体推进剂颗粒的三维网格自动生成

为了分析圆柱颗粒对复合固体推进剂性能的影响,必须对复合固体推进剂中三维圆柱颗粒进行网格划分。根据"力平衡"网格划分原理设计了一种适合于三维圆柱颗粒的网格自动剖分方法。编程运行结果表明,方法能自动完成划分,网格质量较高,满足要求。     

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

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

RBFs-MSA Hybrid Method for Mesh Deformation

Simulating unsteady flow phenomena involving moving boundaries is a challenging task,one key requirement of which is a reliable and fast algorithm to deform the computational mesh.Radial basis functions(RBFs) interpolation is a very simple and robust method to deform the mesh.However,the number of operations and the requirement of memory storage will be increased rapidly as the number of grid nodes increases,which limits the application of RBFs to three-dimensional(3D) moving mesh.Moving submesh approach(MSA) is an efficient method,but its robustness depends on the method used to deform the background mesh.A hybrid method which combines the benefits of MSA and RBFs interpolation,which is called RBFs-MSA,has been presented.This hybrid method is proved to be robust and efficient via several numerical examples.From the aspect of the quality of deforming meshes,this hybrid method is comparable with the RBFs interpolation;from the aspect of computing efficiency,one test case shows that RBFs-MSA is about two orders of magnitude faster than RBFs interpolation.For these benefits of RBFs-MSA,the new method is suitable for unsteady flow simulation which refers to boundaries movement.     

飞机结构加筋壁板参数化有限元模型的实现

提出一种基于MSC.Nastran/Patran的加筋板结构有限元模型自动生成方法.方法对加筋条位置、加筋条几何尺寸、底板几何尺寸、材料物性以及网格密度完全参数化,给出标量参数、向量参数、字符参数及动态数组在PCL函数中的具体应用,载荷、约束以及材料特性全部施加于几何体以实现网格密度的自动疏密控制.实例表明,对于帽型加筋板,张角为60°时壁板位移最小.     

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

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

压力载荷下的结构拓扑-形状协同优化

 压力载荷作用下的结构轻量化设计是工程中的常见问题,由于压力加载面的可设计性,现有以固定载荷为基础的拓扑优化技术不能很好地处理这类问题。直接采用CAD参数化样条或B样条曲线描述压力加载面,通过拓扑和形状变量的联合优化满足了工程实际对结构轻量化与边界的功能性与光滑性设计要求。同时,为了避免结构边界形状变化时有限元网格刷新引起的定义拓扑伪密度变量的困难,用所提出的背景网格和密度点技术实现了每一步单元密度设计迭代结果的自动传递,并采用网格变形技术实现了形状设计变量灵敏度分析。采用4个数值算例验证了方法的有效性,其中发动机承力框架的设计结果充分说明该方法在航空结构设计中的重要应用价值。     

叶片参数化有限元建模与罩量优化设计

基于网格变形技术提出了一种叶片参数化有限元建模方法,该方法具有快速实现叶片参数化有限元建模的优点,同时能够在保证叶片有限元网格拓扑形式不变的前提下,有效解决叶片优化设计过程中的网格自动更新问题.在此基础上,结合叶片罩量优化设计的灵敏度分析技术,建立了叶片罩量自动优化一般设计流程,开展了叶片实例优化设计应用研究.优化结果表明,相对于初始设计,罩量优化设计可以有效地改善应力分布,降低叶片应力水平.      

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

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

结构自动离散的一种新方法

本文采用等参元素中的几何插值方法,以少量的输入数据,自动生成全部元素的节点坐标;采用“棋盘”网格,使元素和节点的编号简单易行;同时,采用“波前法”解方程组,使编号的程序简化。 本方法具有输入数据少,改变网格容易,边界节点坐标足够准确等优点,已用于轴对称体二维网格及透平叶片的三维网格,结果满意。