基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

网格自适应策略在高升力构型计算中的应用

提出了一种基于多尺度网格误差分析的新型混合网格自适应局部加密策略。在初始网格基础上,通过多重网格计算,比较次细网格和细网格的计算结果,揭示出网格敏感区域,然后对这些网格敏感区域加密,计算自适应加密后的网格得到最终计算结果。相比于常规网格自适应策略,提出的新策略特别适合多种流动现象并存的高升力构型绕流流场计算中的敏感区域判断。NASA高升力全展构型的数值试验充分表明这种新策略的有效性。     

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

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

三维Euler方程非结构自适应网格投影和光顺技术研究

研究了三维Euler方程非结构自适应网格的投影和光顺技术.以ONERA M6机翼为计算模型,非结构网格为计算网格,守恒变量的范数为自适应加密准则,对采用一般自适应计算策略得到的网格和数值结果进行了分析,并发展了一种考虑网格投影和光顺的自适应计算策略.数值试验结果表明:通过对自适应所生成网格的投影与光顺,实现了自适应计算过程中的几何保真性,获得了相对光滑的压力分布曲线.     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

用混合网格模拟绕X-38粘性流动

采用混合网格，利用DLR所发展的计算软件TAU对X-38高超声速粘性绕流进行了数值计算，网格自适应技术被用来生成有效的计算网格及提高对流场细节的分辨率。计算结果与现有的实验结果和用结构网格得到的结果吻合很好。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

CFD计算网格误差分析的一个算例

为了给一般工程问题的CFD仿真的网格划分提供参考,研究了光滑直圆管中粘性稳态等温流动在不同尺寸、不同纵横比网格情况下的CFD计算数值解的误差.采用FLUENT软件进行数值计算,应用CFD仿真验证和确认的基本方法,与精确解比较,计算了各种网格的数值计算结果的误差,得到了网格独立解的网格尺寸,给出了误差随网格尺寸变化的规律.     

国家数值风洞(NNW)工程中的黏性自适应笛卡尔网格方法研究进展

相比传统的结构和非结构网格,笛卡尔网格具有自动化和高质量兼顾的优势,是未来网格技术发展的重要方向。面向笛卡尔网格生成软件开发需求,依托国家数值风洞（NNW）工程,介绍了黏性自适应笛卡尔网格技术相关课题的研究工作进展,重点从网格生成技术、自适应方法、黏性物面边界处理方法等方面展开。网格生成技术方面,基于全线程树数据结构的思想,开展了叉树数据结构的优化;针对物面离散单元信息的快速检索,构建了稳定、平衡、高效的叉树式物面网格数据结构;对于网格单元类型的判断,发展了染色方法和包围盒方法以提高效率;针对大规模网格量,发展了高性能并行计算技术。自适应技术方面,发展了基于几何特征和流场解特征的自适应方法,建立了流场结构捕捉判据,并开展了三维复杂构型的应用;为了降低全局网格规模,发展了各向异性自适应方法,构建了各向异性特征识别的综合判据;发展了边界层法向射线加密技术,以减少边界层内的网格数目。物面边界处理方面,发展了重叠笛卡尔网格方法,构建了重叠网格尺度匹配技术和数据交互方法;开展了壁面函数方法在笛卡尔网格框架下的适应性应用,以放宽壁面网格尺度限制、降低边界层网格规模,并通过典型算例验证了方法的可靠性。     

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

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

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

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

三维复杂无粘流场的自适应八叉树结构直角网格算法模拟

运用自适应八叉树结构直角网格对三维复杂流场进行了划分 ,并采用点面判别法进行外法向量的判断 ,使与之关联的与物面相交网格的体积计算更为简便准确。在此基础上 ,采用以中心差分为基础的 Jameson的有限体积法进行 Euler方程求解 ,在通量计算过程中 ,采用以面为基础的算法 ,减少了计算工作量。对多个外挂、翼身组合体及全机等复杂外形绕流流场进行了模拟 ,均取得了较好的效果     

Adaptive mesh refinement computation of acoustic radiation from an engine intake

A block-structured adaptive mesh refinement (AMR) method was applied to the computational problem of acoustic radiation from an aeroengine intake. The aim is to improve the computational and storage efficiency in aeroengine noise prediction through reduction of computational cells. A parallel implementation of the adaptive mesh refinement algorithm was achieved using message passing interface. It combined a range of 2nd- and 4th-order spatial stencils, a 4th-order low-dissipation and low-dispersion Runge–Kutta scheme for time integration and several different interpolation methods. Both the parallel AMR algorithms and numerical issues were introduced briefly in this work. To solve the problem of acoustic radiation from an aeroengine intake, the code was extended to support body-fitted grid structures. The problem of acoustic radiation was solved with linearised Euler equations. The AMR results were compared with the previous results computed on a uniformly fine mesh to demonstrate the accuracy and the efficiency of the current AMR strategy. As the computational load of the whole adaptively refined mesh has to be balanced between nodes on-line, the parallel performance of the existing code deteriorates along with the increase of processors due to the expensive inter-nodes memory communication costs. The potential solution was suggested in the end.     

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

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