流场自适应叉树网格数值模拟三维复杂超声速流场

建立了基于分区结构网格的三维贴体叉树形网格的数据结构,并阐述了在此基础上进行自适应分裂/合并判别方法.为节省网格量和保证流场结构捕捉质量,提出对自适应程度进行区域性控制,以及对流场结构进行"保护"性预加密的优化方式.通过应用该网格对三维复杂超声速流场算例的计算,证明该方法对网格加密控制方便,对流场结构分辨率高.     

跨音压气机中非结构网格自适应的研究

跨音压气机中存在激波/边界层相互干涉、激波/叶尖泄漏涡相互干涉等复杂流动结构,为了实现高精度的数值模拟,精细的捕捉到激波以及叶尖泄漏涡附近的流动结构,从而探究其复杂的流动机理,本文发展了非结构网格自适应技术,其可以自动的加密关注流动特征附近的网格,从而实现对激波、叶尖泄漏涡及激波与其他流动现象相互作用的流动结构的精细捕捉。本文将网格自适应技术应用于RAE-2822翼型二维激波算例以及真实跨音速压气机Rotor-37算例的数值模拟中,计算结果表明网格自适应技术能自动的捕捉到激波、叶尖泄漏涡及其与其他流动现象相互干涉的精细流动结构。     

水斗曲面三维非正交贴体网格数值生成

为实现冲击式水轮机水斗内非定常流动的数值解析,需要高密度匀布非正交贴体坐标网格来拟合水斗曲面.本文开发了一种用于水斗曲面三维非正交贴体坐标匀布网格的数值生成程序,其主要算法是基于曲线与曲面的微分几何性质提出的四点投影法.程序可根据需要计算出不同密度的匀布贴体网格点.通过对两不同水斗的算例结果显示数值生成的网格很好地拟合了水斗不规则曲面,同时计算出各网格点的切向矢量和单位法向矢量.     

基于能带结构的布里渊区非均匀四面体网格生成技术

基于非局域赝势能带计算及四面体网格单元中能量满足线性关系,提出一种布里渊区非均匀四面体网格产生方法,在满足精度条件下能自动得出数目最少的四面体网格,使布里渊区积分计算的精度和效率大为提高．通过对硅、锗两种金刚石结构半导体简约布里渊区所产生网格的比较,表明该方法可以根据能带结构的特点自动生成优化的非均匀网格．对现有的态密度四面体计算公式进行了补充完善,并根据生成的网格和完善后的公式计算了硅和锗导带第一、二能带的态密度．     

一种基于预估-修正策略的非定常流网格自适应加密/稀疏方法

给出一种非定常流动数值模拟的网格自适应处理方法.在"求解流动方程-自适应调整网格"的流程中,引入预估-修正步.根据自适应周期内每个时间步上的流场预估解,计算单元上的事后误差估算值.建立考虑解演变的网格自适应指示器,并进行多层次单元加密-稀疏的动态网格自适应处理.在自适应网格上重新计算流场.每个自适应周期中,流动演变区域的网格获得加密;而前一个周期中的特征现象已离开区域的网格被稀疏.应用边界非协调的当地DFD(Domain-Free Discretization)方法求解流动方程.为验证网格自适应处理方法,针对静止圆柱和自推进游鱼的流动进行了数值实验.     

二维非结构网格无粘流场数值计算

1前言非结构化网格流场数值模拟是近十几年来计算流体力学发展的重要标志之一。非结构网格与结构网格相比，在网格划分的灵活性上远远超出，这正好是复杂区域划分所必需的。非结构网格可在某一局部加密，修改而不影响其余部分，用于自适应网格时其方便处也是结构网格所无法比拟的。非结构网格划分的灵活性和加密、修改的方便，使它得以在复杂几何结构流动的计算中得以充分应用，十几年来成为计算流体力学界的潮流，得到了大发展［‘－‘］本文在三角划分的计算域上，进行了二维无粘流场数值计算。采用Roe格式离散二维Euler方程组，三步Run…     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

模型加力燃烧室燃烧流场的并行计算

本文采用基于MPI的并行算法,采用动态内存分配、分区算法和多点重合交错网格系统,在贴体网格下对带V形槽稳定器模型加力燃烧室素流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。计算结果表明并行计算对复杂形状的化学反应流动计算效率很高,是模拟大规模的燃烧问题的有力工具。     

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

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

基于自适应非结构化网格的VOF方法

本文提出了一种基于自适应非结构化网格的VOF算法,根据相界面网格的相函数值对相界面网格进行自适应细化与合并,通过基于非结构化网格的界面构造方法构造相界面,在适量增加网格单元数量的情况下提高了计算的精度。该方法随着时间及相界面的变化无需重新整体生成网格,算法效率较高。经典算例的验证结果表明,本文自适应网格方法计算得到的结...     

强爆炸火球辐射流体自适应网格高精度数值模拟

提出一种自适应网格方法,应用于基于Euler方法的强爆炸辐射流体高精度数值求解.通过与Zinn数值结果对比,验证该方法的正确性.研究自适应网格对冲击波和光辐射输出模拟精度的影响,对比不同网格尺度下的计算耗时.在相同的条件下,使用自适应网格与均匀网格加密3倍得到的冲击波超压分布、光辐射输出演化接近,计算效率提升约8.5倍...     

气固两相各向同性湍流的自适应网格模拟

气固两相湍流是一个复杂的多尺度,多物理场的难题.为了捕捉流场的小尺度涡以及其对颗粒的作用,传统上采用了静态的高分辨率网格,消耗极大的计算资源.本文基于小波算法对网格进行加密或疏化,从而实现网格的自适应.基于自适应网格的框架上,建立颗粒与流场分块的映射,搭建了欧拉-拉格朗日的气固两相流模拟的并行平台.通过对均匀各向同性湍...     

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

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

先进燃烧室涡流杯流场计算

本文采用贴体坐标系（Ｂｏｄｙ－ＦｉｔｔｅｄＣｏｏｒｄｉｎａｔｅＳｙｓｔｅｍ）数值研究先进燃烧室涡流杯内速度场和紊流特性,并利用ＴＴＭ法生成贴体网格。由于双级旋流器形状复杂,本文提出了型线定点法来确定边界网格、ｕｖ线法生成交错网格以及整体不分区进行流场计算。计算结果表明;计算方法合理,计算机程序是可靠的。     

贴体坐标系下模型加力室的大涡模拟

本文利用贴体网格对带V形槽稳定器模型加力燃烧室素流化学反应流流动进行大涡模拟的研究。采用区域法生成模型加力燃烧室的二维贴体网格,并采用多区域耦合法进行区域之间的数据传递,求解加力室整体流场。采用k方程亚网格尺度模型和亚网格EBU燃烧模型分别估算其亚网格紊流粘性和化学反应速率,用热通量辐射模型估算辐射通量,并用交错网格下SIMPLE算法和混合差分格式求解离散方程,壁面函数处理固壁边界条件。计算结果显示了稳定器后面的回流区气流结构,所得的热态流场模拟结果与实验比较吻合,表明采用贴体网格对模型加力燃烧室进行大涡模拟能真实反映流体流动及燃烧过程。     

一种高精度保界的守恒重映方法

Lagrange方法中，当流场发生大变形时，跟踪流体运动的Lagrange网格发生扭曲，使计算无法进行下去，此时必须重分网格，把网格修复成较好的形状。另外，网格自适应技术中的重构、合并与加密，以及同一问题不同程序相继计算的连接，并行计算中相邻块边界区域的数据传递等，这些情况都需要利用旧网格上的物理量来确定新网格上的物理量，是一个物理量重映过程。质点重映方法是基于物理上守恒规律的一种离散的物理量守恒映射方法，既可实现分片常数分布的一阶精度重映计算，又可实现分片线性分布的二阶精度重映计算。这种方法可严格保证守恒量的守恒性，且可以实现任意多边形网格以及节点上物理量的守恒重映。但是，基于分片线性分布的二阶精度重映方法，如果新网格的守恒量没有进行保界调整，那么相应的强度量有可能在其局部的限制范围之外，破坏了原网格物理量的单调性。因而，对二阶精度的质点重映方法进行了进一步研究。在分片线性分布的基础上，将基于结构网格的保界算法扩展到非结构网格上，给出了二阶保界的质点守恒重映方法。     

结构网格与非结构网格联合技术及在内爆问题中的应用

对纯拉氏方法，格式的健壮性、精度和效率强烈地依赖于初始网格的生成。结构网格与非结构网格的联合技术，不仅可适应于复杂区域的网格生成，而且可根据实际问题的特殊性，在初始自适应的生成网格。本文在研究结构网格和非结构网格生成的基础上，研究了结构网格与非结构网格联合使用的自适应网格生成技术。首先是根据实际问题的特征，把整个计算区域划分成若干子区域，其次是把每个子区域分成若干子块，最后根据每个子块的特征生成网格，再把子块拼接起来。     

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

粘性不可压流体流动问题用直角坐标网格的贴体解法

研究一种新的全贴体的求解粘性不可压流体流动问题的非结构化直角坐标网格方法.该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的.这有别于目前流行的各种非结构化直角坐标网格方法.通过对两个典型流动问题的计算对该数值方法进行验证.对比结果表明,本方法计算的结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法是可行的.还对二维钝头体周围的流场进行了计算,计算的流场与STAR-CD的结果相当吻和,说明该算法还可计算较复杂的流动现象.     

平面斜激波固壁反射的镶嵌网格解法

镶嵌网格法是自适应网格法的一种。本文提出的镶嵌网格法以Ni的二阶有限体积法为基础,用密度一阶差分作为察觉激波的特征参数,给出网格局部加密(一分为四)的原则,构造了较节省存贮的数据结构。用于平面斜激波固壁反射的计算结果表明:镶嵌网格对激波的分辨能力不仅比粗网格强得多,而且也稍好于细网格;计算时间仅为细网格的一半稍多些.可见,此方法对解决计算效率和精度的矛盾是有利的。