三维非结构粘性网格生成方法

描述了一套适合粘性流动计算的三维非结构网格自动生成方法.在物面附近的粘性作用区域,采用推进层方法生成各向异性的"扁平"四面体网格,并通过一定的网格伸长控制参数,实现整个流场区域网格高度的平滑过渡.当粘性网格的推进高度达到预定要求时,推进层方法自动停止,转而采用阵面推进方法生成常规意义的尽量接近正四面体的各向同性网格.同时给出了利用该方法生成的M6机翼非结构粘性网格来求解机翼粘性绕流的简单算例.     

三维非结构网格自动生成

采用各向异性平面非结构网格生成技术对曲面的参数平面进行三角剖分,从而得到曲面的非结构网格,作为三维非结构网格的边界网格.应用推进面法生成网格内点,增量法将生成的内点逐点插入现有网格进行网格细化,得到三维计算域的Delaunay非结构网格.讨论了非结构网格质量优化方法.给出几个算例说明方法的应用.     

非结构网格的并行生成

研究了非结构网格的并行生成方法.改进了R.Lohner的"波阵面"区域分裂算法以使子网格及其边界更有益于网格并行生成,提出了边界并行优化策略,改善了边界的光滑性;完善了子区域内生成网格时接受新点及新单元的条件,在界面网格生成过程中提出只接受新单元而拒绝新点的策略,节省了机时.     

二维非结构网格生成及Euler方程计算的方法研究

用离翼型表面最小作为阵面推进法中的参数选择依据。生成二维问题的非结构网格。这种方法戏了传统的背景网格观念,直接提供网格生成过程中所需的背景信息。在求解Ｅｌｕｅｒ方程时,用格心格式的有限体积法作空间离散,用四步Ｒｕｎｇｅ－Ｋｕｔｔａ作时间推进,采用不同的加速收敛措施获得定常流动。提出了两这界条件的构造办法,并 不同边界条件对结果的影响。     

八叉树网格到非结构混合网格的转换

在数值模拟中, 非结构网格的优势是可以采用相同的数值格式统一处理任意复杂的计算区域, 但在网格生成过程中难度大, 并且不容易控制网格质量。树结构网格可以认为是介于结构网格和非结构网格之间的一种网格, 目前已经有相对成熟的方法快速在复杂区域内生成二维四叉树网格和三维八叉树网格。在实际应用中, 数值方法往往需要在连接协调的非结构网格上做离散, 树结构网格中不同尺寸的网格之间连接不是协调的, 在应用上会受到很多限制。文章实现了树结构网格到非结构混合网格的转换, 这种转换在二维情况下就是将四叉树网格转换为非结构三角形和四边形的混合网格, 三维情况下则将八叉树网格转换为非结构混合网格。这一转换过程的难点在于需要考虑数千种不同的八叉树单元, 并给出能实现连接协调的非结构混合网格划分。可以出现的网格单元包括六面体、三棱柱、金字塔和四面体这4种不同情况。通过特别的分类, 实现了程序的自动生成, 这种程序自动生成技术一方面可以避免人工编写大量程序时的失误, 另一方面也使得对数以千计的不同情况的处理成为可能。通过对几个简单网格的测试, 对网格数据转换方法做了初步的验证。      

三维复杂外形的非结构网格自动生成技术与应用

给出了一套三菜的非结构网络自动生成方法。任意外形由几种代表性的曲线和曲面来描述,将曲面变换为平面,在平面内由二维阵面推进法生成三角形网格,然后反变换为物形表面网格。空间的四面体网格由三维阵面推进法生成。四面体网格的优化采用了附加“关联质量”约束的节点松驰的局部的Ｄｅｌａｕｎａｙ变换技术。另外还发展了非结构风格的自适应技术。数值算例表明了方法的灵活性和通用性     

二维非结构网格生成及自动加密技术

介绍一种可对任意二维区域进行非结构网格划分的方法,提出关于网格调整光滑的新思想,可对Laplace型滤波算子赋予新的意义。在此基础上再有目的选择松弛因子的大小,可达到同样的光滑效果而迭代次数显著减少。利用文献[3]提供的自动加点算法,对该方法稍加改进,就能完成自适应过程。利用该方法编制的程序,只需提供边界点座标,就能快速地得到满意的网格。     

二维矢量边界推进非结构网格生成方法

提出二维矢量边界推进生成非结构三角形网格方法并证明其可行性.根据流场边界尺度布置边界节点并运用矢量边界推进方法生成背景网格,运用符号面积函数和概率筛选方法布置初始点阵,提出Spring-Laplace方法优化节点位置,同时利用边交换技术优化网格结构.该方法可包含任意点源、线源和内嵌边界,可自由进行局部自适应加密或稀疏,实现任意平面域内与尺度要求一致的高效光滑三角网格剖分.     

一种新的非结构四边形网格生成方法及软件

提出一种间接生成非结构四边形网格的方法.采用Delaunay三角形网格生成的算法生成三角形单元,在三角形单元中定义生成四边形网格的前沿,根据三角形的特点(边、角大小)以及生成四边形前沿的方向,将两个三角形合并为一个四边形单元.此方法的特点是不增加节点,仅改变点与点之间的拓扑连接,使单元基本变为四边形占优的单元.然后研究四边形占优的网格转换为100%的四边形网格的转换模版.并应用于实际物理问题.数值实验显示了该方法生成非结构四边形网格的特点.     

非结构子网格压力方法

纯拉氏流体力学计算失败的根源是砂漏运动和虚假涡旋,运用子网格拉氏质量与其相关的子网格密度得到子网格压力,子网格压力产生的扰动隅角力能消除这些虚假运动,以便消除困扰拉氏流体力学算法的网格畸变带来的数值计算困难.将包括子网格压力的相容流体算法推广到非结构网格,研制包括非结构子网格压力的相容流体程序,对Saltzman活塞问题进行数值模拟,从数值结果分析可知,非结构子网格压力方法能抑制虚假涡旋,消除砂漏畸变.     

三角B-B曲面在非结构网格生成中的应用

基于生成非结构化网格的Delaunay三角形化方法,应用三次三角Bernstein-Bézier曲面来控制网格点的分布情况。对于给定边界,利用三次非均匀B样条进行边界的拟合及重新离散。对初始化形成的Delaunay三角形,应用三次三角Bernstein-Bézier曲面来计算位于其中的点的长度标尺,通过对三角Bernstein-Bézier曲面边界点的定义来控制三角形内长度标尺的分布。对复杂通道网格剖分的实例表明,此方法可以很好地控制内部网格点的分布情况。     

含运动物体流场计算的自适应非结构二维网格生成方法

采用自动插点的Delaunay方法和局部网格重新生成方法,对含运动边界自适应非结构二维网格生成方法进行了系统研究,实现了带有非定常运动物体流场的数值计算,并进行了数值实验.     

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.     

三维非结构网格Euler方程的LU-SGS算法及其改进

将LU-SGS隐式时间推进格式运用到非结构网格Euler方程的求解中,并对传统LU-SGS格式进行改进,结合网格重排序,发展了一套效率更高的三维Euler方程求解器.以M6机翼及超临界LANN机翼的跨音速无粘流场为算例,将改进LU-SGS格式与四步龙格-库塔显式格式及传统LU-SGS格式进行了比较.计算结果表明:所有格式的计算结果与实验结果都符合很好;传统LU-SGS格式计算效率为显式格式的3倍多,而改进LU-SGS格式计算效率为显式格式的7倍多.     

非结构多边形网格解扭和重分方法

在多维流体动力学计算中，流体运动和计算网格的关系可以分为两种情况。一是Lagrangian方法，即网格跟随流体运动；二是Eulerian方法，即流体流过固定；下动的网格。一般计算网格的运动是任意的。这就对应于任意Lagrangian—Eulerian（ALE）方法。ALE方法的核心是通过调整网格运动，使得数值模拟的精度、效率有所提高。它的主要步骤是：显式Lagrangian步；网格重分，即得到新的计算网格；物理量重映，即将Lagrangian步的计算结果变换到新网格上。在这3步中，较少研究网格重分。数值模拟和网格重分的一个基本前提是网格是合理的，或者说网格不能发生翻转，网格应当是凸的。而Lagrangian步数值模拟会造成网格扭曲，因此在网格重分前进行网格解扭是十分必要的。文中描述了通用的网格解扭、重分算法，使得解扭、重分后的网格有较好的几何品质，同时尽可能接近Lagrangian网格。     

基于约束Delaunay三角化的二维非结构网格生成方法

给出基于局部重构和边交换技术的两种约束Delaunay三角剖分方法并证明其收敛性.采用边界指示法恢复流场形状;在预设尺度的指导下融合流场边界曲率、中轴线、梯度限制等信息修正流场尺度;运用Spring方法布置边界点,通过符号面积函数和概率筛选法布置计算区域节点;运用Spring-Laplace方法优化节点位置,伴同边交换和边吞噬技术优化网格结构.该方法可自由进行局部自适应加密或稀疏,并应用于映射曲面网格生成和移动网格技术.     

BGK方法在三维非结构网格上的初步应用

将BGK计算方法从二维拓展到三维并且应用于三维非结构网格,具有重要的理论价值和实用价值.采用旋转局部座标的方法,发展了一种针对三维非结构网格的BGK计算方法.在计算过程中,将最小二乘法应用于三维非结构网格的导数计算.对三维激波管和三维欠膨胀垂直冲击射流等两个算例进行了细致分析.这两个算例的计算结果表明,该方法在三维非结构网格上的初步应用是成功的 关键词： 气动BGK方法 三维 非结构网格     

非结构网格的生成及新型数据类型的应用

本文基于 Ｄｅｌａｕｎａｙ三角化方法及局部重组的阵面推进法，给出了任意形状平面区域的非结构化网格剖分方法。文中对加点方法及数据结构和存取模式进行了改进，同时应用链表结构、动态数组和派生类型的数据变量，对算法进行优化，节省了计算内存，提高了网格的生成效率。本文研究为特殊计算区域的三维非结构化网格的生成提供了基本方法。     

非结构网格自适应的新策略

对于用推进波阵面法生成二维非结构网格,提出了一种网格自适应的新方法。在网格需要加密处,布置一些小的三角形单元,并以这些单元的边作为初始推进波阵面的一部分。在推进过程中,小单元附近自然生长出尺寸逐渐变大的单元。对有激波的流动,在网格图形上,可形成十分清楚的激波图象。     

多重网格法在非结构网格中的应用

将多重网格法运用于非结构网格.网格是通过聚合法得到的,网格之间是相互关联的.方程的求解采用Jamson的有限体积法.给出了二维、三维情况的数值算例.