基于非结构重叠网格的二维N-S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

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

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

二维Euler方程的自适应非结构网格生成方法

发展了一种用于数值模拟二维Euler方程的自适应非结构网格的生成方法。这一方法根据背景网格提供的信息直接确定内部结点的位置，然后依照Delaunay准则，将新结点加入到已有网格划分中。用于平观流及管道流的实际计算中，结果表明，该方法不受流场形状及流动特征的限制，具有较好的通用性。     

AECSC-J ASMIN湍流燃烧仿真软件研发和检验

航空发动机燃烧室几何结构复杂,湍流和化学反应存在强烈非线性相互作用,需要对流动和燃烧及其相互作用进行高精度高时空分辨率的刻画,目前燃烧室湍流燃烧数值模拟仍然是高难度的瓶颈问题之一。介绍了由北京航空航天大学航空发动机数值仿真研究中心、北京应用物理与计算数学研究所和中国工程物理研究院高性能数值模拟软件中心联合研发的AECSC-JASMIN软件主要框架、算法以及针对该软件的算例检验。在Sandia射流火焰、支板火焰和单头部燃烧室检验算例中,对比实验数据,射流和支板火焰预测结果与实验值一致;支板算例的平均相对误差在15%之内;单头部燃烧室模拟结果符合物理实际,总压损失与实验值基本吻合。说明AECSC-JASMIN软件可用于复杂结构高分辨率高精度湍流燃烧数值模拟。     

Scramjet整机流场的非结构网格并行数值模拟

为了得到某支板喷氢超燃冲压发动机内流场的细节,采用非结构网格CFD隐式解算器,耦合求解Na-vire-Stokes方程、组分连续方程和M enter的k-ω应力剪切传输(SST)湍流模型,对整机模型的冷流和三个当量油气比的湍流反应流进行了三维并行模拟。给出了流线和水的分布云图,分析了支板附近旋涡和推进流道一维质量加权平均的马赫数、静温、组分质量分数及燃烧效率等参数随当量油气比的变化,并将壁面压力和实验值进行了比较,吻合较好。除=0.8在支板附近出现了小范围亚声速外,其余状态整个内流道超声速占主流。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

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

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

三维动态非结构重叠网格Navier-Stokes方程并行算法

 提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。     

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

非结构网格生成技术方法研究

利用点源或线源所影响的势函数分布规律,通过建立势函数与网格背景信息之间的关系,提供非结构网格生成过程中所需的网格背景参数;在物体边界节点已知的情况下,适当布置点源和线源个数,用最小二乘法确定各点源和线源的强度;通过对部分点源或线源的影响势函数施加方向控制开关,使所生成的网格在期望的区域得到加密;算例表明本文方法是一种有效、快速的处理方法,与传统用背景网格提供背景信息的方法相比,在计算量、程序量、光滑性和灵活性方面都能表现出良好的性质。     

基于非结构网格的跨声速静气动弹性计算

基于非结构网格开发了静气动弹性计算程序,通过气动力与弹性变形的交替迭代计算实现流固耦合。气动力求解采用RANS方程,结构变形计算采用柔度系数法,通过面样条(IPS)方法进行柔度系数插值,使用改进弹簧近似法实现网格变形。通过对AGARD 445.6机翼的静气弹计算验证了方法的正确性,同时对机翼的静气动弹性特性进行了分析,计算表明,弹性变形使其升力线斜率下降、焦点前移。     

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

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

基于动网格的涡轮流场仿真方法及特性分析

为提高涡轮流场仿真计算的收敛性和结果的精度,基于动网格的流场仿真技术,探讨了涡轮瞬态流场仿真的 2步方法——先稳态后瞬态计算。计算结果表明:在转子—静子间滑动界面上出现了计算方程守恒引起的流场参数阶跃,流场参数变化趋势与一维等熵流动理论分析结果一致。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法,前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率,而后者保证了网格间边界附近流场求解精度;本文非定常Euler方程的求解采用有限体积法,通量计算采用Osher格式,时间积分采用双时间步长推进。为验证本文提出的方法,分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率,保证了网格边界处的精度,具有较高的工程应用价值。     

基于DTF和REDIM技术构建的亚网格燃烧模型及其验证

将反应扩散流形(REDIM)技术与动态火焰增厚(DTF)模型结合起来,构建了一种新的湍流预混燃烧模型。用该模型对PRECCINSTA燃烧器全局当量比分别为0.75和0.83的两种贫燃预混旋转火焰开展大涡模拟计算,将计算结果与原DTF模型计算结果以及实验结果进行对比研究,着重分析了速度、温度以及CO2质量分数的径向分布情况。研究表明:新模型的计算结果与实验结果吻合良好,由于新模型求解的组分输运方程数目减少,计算效率比原DTF模型约提高15%;理论上,新模型可推广应用于部分预混火焰的计算,只要采用合适的REDIM化学反应表格,如二维的REDIM表格。     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

基于贝叶斯因子和二阶概率方法的圣地亚热传导模型确认挑战问题

针对美国圣地亚国家实验室为促进模型确认方法的发展所提出的热传导挑战问题,根据现代模型确认思想总结了该挑战问题解答中应包含的核心内容。在确认度量环节采用贝叶斯因子考察了实验数据是否支持所给模型,在此基础上通过二阶概率方法得到了模型预测的分布,以此计算出模型预测结论的置信度。该过程中考虑了模型参数的随机不确定性和认知不确定性,最后通过灵敏度分析辨识了模型参数的不确定性对模型预测的影响。研究表明该挑战问题中的实验数据支持所给模型,模型预测受导热系数的不确定性影响最大,模型预测材料在调控条件下的失效概率不满足调控要求,该结论的置信度为99.97%。     

仿真模型验证实用方法研究

仿真模型VV&A是仿真模型开发中的一项重要工作 ,贯穿于模型开发的全生命周期。本文重点讲述了实际数据序列和仿真数据序列比较法和理论比较法 ,并给出了相应的例子来说明这两种方法的使用情况。