粘性超声速绕体流动的有限谱法计算

首先推导出粘性可压缩流体动力学基本方程组在广义坐标下矢通量分裂的一般形式，并将有限谱ENO格式用于求解超声速流动问题。通过典型算例，验证了这一格式的精度和稳定性。然后进行了圆柱和机翼超声速绕流的数值模拟，得到了与实验和其他计算相符的结果。这表明有限谱法不单具有精度高，应用灵活等多方面的优点，而且突破了传统谱方法全域性的局限，可以成功地应用到一般形状流场可压缩流动的数值计算中。     

基于当地流动结构的RANS/LES混合模型

针对传统雷诺平均Navier-Stokes/大涡模拟(RANS/LES)混合模型固有的RANS和LES区域衔接难题,借鉴Vreman提出的亚格子(SGS)模型,提取其中描述流动状态的判断函数,利用不可压缩槽道湍流标定模型系数,构造了一个新的RANS/LES混合模型。评估算例包括不可压缩槽道湍流、稳态超声速平板边界层、NACA4412翼型绕流和亚临界雷诺数圆柱绕流。该模型作为壁面模化LES,解决了传统RANS/LES方法中速度型偏离对数律(LLM)的缺陷,且拓宽了原Vreman模型在粗糙网格上的适用性,具有较高精度。该模型同时解决了模化应力折损问题,且对大尺度分离流动的模拟精度相较传统的RANS/LES混合模型如脱体涡模拟(DES)有进一步提升。     

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。     

基于高阶间断伽辽金数值方法的浸没边界法

为了降低复杂边界条件下贴体网格生成难度，近年来基于笛卡尔网格的浸没边界法逐渐成为研究此类问题的主要数值方法之一，然而计算精度和计算效率仍然是此类方法目前面临的挑战。与传统的空间2阶精度的有限体积格式相比，空间精度为3阶或3阶以上的高精度方法具有空间精度高、数值分辨率高、数值耗散小的优点，而作为高精度数值方法之一的间断伽辽金方法在浸没边界方面的应用仍较少。本文将高阶间断伽辽金方法的高精度优势与浸没边界法无需贴体网格的优势结合起来，提出了适用于可压缩流动的高精度浸没边界法。其中边界条件采用体积惩罚方式实现，同时采用牛顿法迭代以及MPI并行提高计算效率，物面的数据重建采用插值点处的逆距离权重（IDW-IP）方法代替高阶格式下多项式插值方法。本文基于笛卡尔网格测试验证了二维定常和非定常情况下的数值模拟效果，并与传统贴体网格的计算结果进行了对比。     

有限体积格子Boltzmann方法用于近空间连续流区绕流模拟

针对航天器再入解体形成残骸碎片的近空间绕流计算问题,拓展了格子Boltzmann方法在可压缩流动模拟的能力。引入有限体积隐式格式求解耦合双分布函数格子Boltzmann模型方程及采用圆函数为基础构造的D2Q13离散速度模型;引入IMEX-RK格式进行了时间项离散解决模型方程的源项刚性问题;对Riemann问题、平板双马赫反射问题、RAE2822翼型跨声速绕流等近空间连续流区、可压缩典型案例进行了数值模拟。通过比较分析,初步验证了耦合双分布函数有限体积格子Boltzmann方法对连续流区、可压缩流动的模拟能力;进一步开展了方柱形解体残骸的超声速绕流模拟,得到与N-S方程计算结果吻合一致的绕流结果,证实经改进的耦合双分布函数有限体积格子Boltzmann方法对解体残骸碎片绕流问题具有较好的模拟能力。     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

基于四叉树网格下的浸入边界-格子Boltzmann方法

提出了一种基于非均匀四叉树网格下的浸入边界-格子Boltzmann方法.在不同层网格界面上,为了保证物理量的连续性,需要在时空方向上进行插值.由于四叉树同层网格在空间上步长相等、且相邻时间层上时间步长相等,所以在时空方向上采用平均值插值,这样做的优点是时空方向插值既可以达到二阶精度,又可以简化计算过程,节省资源,对任意边界加密下的网格,格子Boltzmann方法的实现比较容易.为了充分利用均匀笛卡尔网格的优势,物面边界的处理采用了速度修正法,与传统浸入边界-格子Boltzmann方法中的直接力法、动量交换法相比,无滑移边界条件得到了较好的保证.同时,把大涡模型加入到浸入边界-格子Boltzmann方法中,实现了在四叉树网格数据结构、边界处理技术、大涡模拟几种模型相结合下,绕障碍物的较大雷诺数流动的模拟.通过不可压缩粘性流中圆柱绕流算例验证,结果与其它方法结果吻合良好.     

SST湍流模型的激波与可压缩效应改进

由于缺乏对某些重要流动特征的考虑,针对不可压流发展的标准SST湍流模型在描述超声速流场时存在明显的局限性。为改善SST模型在吸气式高超声速推进系统内流等复杂超声速流场中的预测精度,基于流动特征结构定向开展了激波和可压缩效应改进。通过激波/湍流边界层判别函数和可压缩效应判别函数定位标准SST模型参数或建模假设失效的区域,针对性地改进湍流模型。采用超声速平板边界层流动、超声速压缩拐角分离流动、超声速隔离段复杂激波串流动以及HIFiRE-2超声速内流等算例进行了测试,结果表明改进模型具有与标准SST模型一致的边界层预测能力,但显著提高了对激波干扰流动及逆压分离流的预测精度。     

求解可压缩流动的同位网格SIMPLE方法研究

在Rhie-Chou动量插值的基础上,推导了同位网格可压缩SIMPLE算法.经过无粘流超音速凸包算例和激波/湍流边界层干扰算例计算发现,如果对流项采用高阶有界HLPA格式,密度插值采用一阶迎风和中心差分的混合格式,这种算法能够很好地模拟凸包超音流的流动现象,在采用了新型GAO-YONG湍流模型后也能够较好地模拟激波/湍流边界层干扰.      

一种求解低速跨流域流固耦合问题的离散统一气体动理学格式

运动边界及流固耦合问题是低速及高速、连续流域及稀薄流域流动中常面临的多场耦合问题。本文在原始离散统一气体动理学格式基础上，采用任意拉格朗日-欧拉方法及动网格技术，并耦合计算结构动力学，实现适用于低速流动运动边界及流固耦合问题的计算方法。在方法验证方面，进行连续流域圆柱绕流强迫振动和自由振动数值模拟，计算结果与参考文献结果吻合良好，准确捕捉圆柱在锁频区和非锁频区的流动特性。在稀薄流域范围，将方法中的连续流D2Q9离散速度模型替换成跨流域离散速度模型，对比低速圆柱绕流强迫振动算例在两种离散速度模型下的数值结果。虽主要关注低速流动问题，但发展的算法仍具备模拟高超声速运动边界问题的能力，对此进行了Ma=5.0稀薄流振动圆柱研究。低速和高速运动边界模拟能力显示了方法在考虑稀薄气体效应流固耦合问题方面的潜在应用价值。     

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

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

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

基于本征正交分解的流场快速预测方法研究

对于飞行器气动设计,通过实验或CFD计算获得所有设计状态处流场信息成本高、耗时长,难以满足工程应用要求。利用基于本征正交分解(POD)降阶模型结合相应的插值方法,只需要通过实验或CFD计算获取有限个数的采样流场,可以快速预测出具有满足精度要求的设计状态处流场信息。首先,证明对POD基系数插值获取预测流场与对采样流场直接插值获取预测流场的等价性,并分析采样流场分布的要求。然后,以三段翼型流场作为研究对象,POD降阶模型分别与三次样条插值和径向基插值结合,实现设计状态处流场信息快速预测。最后,通过将流场预测结果与相同条件下的流场CFD计算结果对比,表明了POD降阶模型结合两种插值方法的各自特性。     

台阶流模型流场的局部DQ-Lagrange求解方法

以Lagrange多项式作为Differential Quadrature(DQ)方法的基函数,应用局部DQ-Lagrange方法对Navier-Stokes方程进行了数值求解,建立了基于局部迎风DQ-Lagrange方法的复杂台阶流流场求解模型.并在此研究基础上,提出了本质和自然边界条件的处理方法.分析了支持域的选择、边界条件处理方法、迎风机制等对求解精度的影响.结果表明:支持域较大时,Lagrange高阶插值易出现振荡现象,可以采用局部替代全局的插值方法;迎风支持域的加入可以提高求解精度;建模时速度支持域的选择应该大于压力支持域.      

燃烧室内二维两相湍流流场数值分析

利用控制容积法求解固体火箭发动机燃烧室内的二相不可压湍流流动，采用多重区域网格法，将复杂的物理域分成两个简单的物理域，分别生成贴体网格。     

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

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

用于复杂流场数值模拟的重叠／嵌套网格方法

介绍了一类用于复杂流场数值模拟的重叠／嵌套网格方法。它是一类多块结构化网格,可用于有或无多重网格加速收敛技术的有限体积、有限差分等计算。其基本思想是把复杂流场的计算区域化分成相对简单的子区域,每个子区域上的网格可以独立生成;网格界面之间的流场信息交换通过三线性插值来实现,通过修改流场解算方法来排除洞中的点。     

基于响应面插值的非线性气动弹性计算

飞翼、大展弦比低雷诺数气动布局容易在小迎角的条件下出现气流分离,会带来明显的非线性气动力问题,同时气动弹性带来的影响亦不可忽略。针对此类布局提出一种建立基于径向基函数插值的非线性压力系数分布模型的方法,利用径向基函数插值建立面元上压力系数对迎角导数的响应面,将压力系数积分并通过无限板样条(IPS)方法进行气动结构多次迭代插值实现非线性气弹分析。结果验证了该方法对于静气动弹性分析的有效性,同时能准确地反映弹性带来气动效率的降低和变形对升力及阻力的影响。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。