贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

移动网格变形法的应用及其修正方法

运用Liao提出的动态型移动网格变形法求解二维双曲型守恒律方程组.由于计算过程中误差的影响,我们无法严格控制网格节点分布,导致实际获得的网格与理论所要求的网格存在偏差,因此为了衡量变形法移动网格的质量,给出质量参数的定义,结合稳态型移动网格变形法对已得网格进行修正,并将该修正方法运用到具体算例,从数值结果可以看出,本方法是可行而且有效的.     

多重网格方法及其在高超声速进气道中的应用

随着计算流体力学的发展,多重网格在求解欧拉方程、N—S方程过程中得到了广泛应用。将这种方法应用于高超声速进气道的数值模拟中,通过具体的算例,验证了该方法的有效性。     

欧拉方程多段翼型绕流计算的嵌套网格方法

本文采用嵌合体埋入技术,将主段翼型和襟翼分成两子域,对多段翼型的复杂构型进行了分区迭代求解,在有限体体积离散的基础上,求解欧拉方程计算各子域流场的解,通过各流场的耦合迭代到多段翼型绕流的解,典型的带３０％襟翼ＧＡ（Ｗ）－１翼型算例表明,该计算方法稳定,能较地模拟多段翼型中主段翼型与襟翼交接处的复杂流动现象计算结果与实验值有较好的一致性。     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

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

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

欧拉方程多重网格方法及三维跨声速大迎角流动数值计算

本文给出一种粗、细网格时间步长可调整匹配至最优组合的多重网格方法以获得在最细网格上最大的加速收敛效益,二重网格中节省机时可达３５％。文中的数值方法可用于飞行器跨声速大迎角无粘有旋的涡干扰流场的模拟研究。     

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程.此方法在一系列粗网格和细网格之间进行值的传递,并且粗网格和细网格之间结构上无关.在每一重网格上用有限体积法进行计算.采用本文的方法数值求解了NACA0012翼型的绕流.     

多重网格法在混合网格中的应用

设计了一种高效的混合网格聚合方法,在此基础上应用多重网格技术来加速Navier-Stokes方程的收敛过程.方程求解采用Jamson的中心格式有限体积法.对NACA0012翼型和ONERA M6机翼的粘性绕流进行了数值模拟,算例结果与试验结果符合良好,加速收敛效果显著.     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

基于TLBM-FVM耦合的飞行器舱内热环境跨尺度预测方法

舱内热环境的有效预测是优化飞行器热控与防热设计、减小系统冗余并保障飞行器热安全的重要基础。在国家数值风洞（NNW）工程支持下,针对目前舱内热环境多尺度、多机制复合传热特点及其数值预测面临的精度与效率提升难题,发展了多区域协同推进的时空耦合模型及流/固界面的自适应分辨率识别算法,建立了基于热格子Boltzmann方法（TLBM）与有限体积法（FVM）相互耦合的舱内复合传热跨尺度预测方法,开展了典型飞行器仪器舱的综合热分析,验证了耦合方法的计算精度及效率。研究表明,相关方法可实现舱内热环境的局部精细化与整体大规模的协同模拟,用于开展整体自然对流与设备局部热量传递的多尺度数值模拟,掌握不同环境参数对舱内热质传递过程的影响规律,从而为飞行器热防护/热管理一体化设计提供重要技术支撑。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

基于迎风格式的三角翼涡流场数值模拟研究

采用基于有限体积方法的高精度迎风格式对层流N-S方程进行了数值模拟,研究了普通三角翼和双三角翼亚音速绕流的流场特性,并对双三角翼的几何外形、网格分布及远场边界对数值模拟的影响等问题进行了初步探讨.     

用高精度紧致差分格式分块耦合求解二维粘性不可压缩复杂流场

本文提出了张量形式的粘性不可压缩流场方程的高精度分块耦合求解方法.在分块算法中用重叠一层的方法来实现子块间的数据交换,同时克服子块剖分所造成的数值奇性.利用本文所提出的算法成功地模拟了椭圆柱绕流问题的时间发展过程,如尾流,旋涡的产生,脱落过程以及Karman涡街的形成等.     

弱可压缩流体与可压缩流体模型的管路水击研究

为提升液体推进系统模型对可压缩流体的仿真能力,针对供应管路的瞬变流动,采用一维有限体积法,使用Modelica语言基于MWorks平台分别开发了适用于补燃循环发动机仿真模型库的弱可压缩流体与可压缩流体的管路模型,不同于控制方程中显含声速的弱可压缩流体模型,可压缩流体模型通过将守恒方程与真实流体的SRK型状态方程耦合求解来计算压力波的传播。在管路供应系统中,分别以液氢、液氧及液甲烷为工质进行关阀水击的仿真研究,结果表明:使用SRK状态方程避免了可压缩流体模型对内能拟合公式的依赖,减轻了建模难度;可压缩流体模型对液氧及液甲烷的仿真结果较为准确,最大误差不超过2%,对液氢的最大计算误差约9%;相对于可压缩流体模型,弱可压缩流体模型计算得到的水击压力在振幅及频率上均偏大;可压缩流体模型计算得到的流体密度波动不超过3%,因此工程实践中处理水击问题时将液体假设为弱可压具有一定的合理性。     

模拟飞机尾喷焰辐射传输的非结构有限体积法

有限体积法因其物理意义明确,在流动与传热的数值计算中有广泛的应用。在利用有限体积法求解辐射传输问题时,不仅能够保证在每一控制体内辐射能量守恒,还可以保证在每一个控制立体角内辐射能量守恒,因而使得模拟结果具有较高的精度;同时,非结构化网格技术是处理复杂几何形状物理问题的一种有效方法。本文编写了求解辐射传输的非结构化有限体积法程序,将计算结果与其他文献的算例进行了对比,证明了模型和程序的正确性,分析了各种因素对程序计算精度和计算时间的影响。在此基础上,将求解辐射传输的非结构化网格有限体积法程序集成到Fluent软件中,计算了飞机尾喷焰的温度场分布。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一.对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%.本文还分析了计算结果与串行计算时的差异,以及利用区域重叠     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

基于响应面法的机翼气动/结构一体化优化设计研究

基于响应面法进行了机翼气动/结构一体化优化设计研究,流动控制方程为三维欧拉方程,采用有限体积法进行数值求解,应力和结构变形采用有限元方法计算,静气动弹性分析采用强耦合迭代方式,响应面模型采用二次多项式来构造.以跨音速M6机翼为初始机翼,进行了多目标、多约束情形下的气动/结构综合优化设计,优化后所得到的新机翼具有更佳的气动/结构综合性能,升阻比增加了9.25%,而重量减轻了4.84%;响应面模型精度满足设计要求,拟合误差均不超过1%;这说明本文所发展气动/结构综合优化设计方法是成功且有效的,具有广泛的应用前