基于指数插值的浸没边界法在可压缩流模拟中的应用研究

由于浸没边界法在处理复杂外形和动态问题时具有极强的灵活性和便利性,因而受到广泛的关注。在当前浸没边界法的研究和应用中,主要工作集中在关于中/低雷诺数和不可压流动的模拟。本文引入 Choi 等提出的指数插值法对目标点速度进行插值,并根据多项式插值和 Crocce-Busemann 关系式等确定目标点压力和温度,并借鉴直角网格法中的虚拟网格法计算物面内的流动变量,建立了适用于可压缩流动的浸没边界法。通过亚、跨声速段翼型绕流和超声速圆柱绕流算例对方法进行验证,结果显示在 k 值取值合适时模拟结果与实验数据较为吻合。通过对相关算例网格收敛性分析,证明了本文发展的浸没边界法具有较高的精度。文中各算例的模拟结果表明,此方法能准确捕捉亚、跨、超声速绕流问题中的各类流动结构,适用于可压缩流动的模拟。     

可压缩算法在燃油冷却数值研究中的应用

为深入研究超燃冲压发动机的燃油主动冷却过程,在可压缩算法框架下,引入预处理技术、低耗散数值离散格式、一般流体状态方程、高压物性模型、燃油裂解模型等,发展了一套适用于求解一般流体低速可压缩问题的数值模拟方法。利用电加热管试验对数值方法进行了验证,并定量研究了PR-3燃油的冷却需求量及平板面板冷却槽道布置的影响。结果表明,数值结果与试验符合很好,证明当前数值方法在燃油冷却数值模拟中是准确有效的。针对350mm×90mm×6mm的平板面板冷却过程,需要约10g/s的燃油可确保1WM/m2热流下的防热安全。槽道布置对温度场影响明显,两侧进油方式比单侧进油的最大温差低约150K,均温性更好,更有利于降低结构热应力。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用.与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要.因此,基于两步四阶时间离散和高精度WE-NO重构,发展了具有四阶时间精度的气体动理学格式.该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力.目前...     

涡轮叶片表面气膜出流的数值模拟

深入理解气膜出流对涡轮叶栅气动性能和流场结构的影响对于冷却涡轮的设计具有重要意义。本文利用数值模拟手段对涡轮叶栅表面气膜出流与主流相互作用进行了研究,给出了冷却孔附近的三维流场,分析了叶片表面气膜冷却机理以及吹风比对气膜出流与主流相互作用的影响规律。研究结果表明,吹风比增大将导致叶栅内部流动损失增大,且压力面更为明显;当吹风比增大到一定程度,叶片两侧面上将出现掺混肾型涡流动结构,并且在压力面和吸力面诱导出旋向相反的诱导涡。     

用BGK格式对平面可压缩自由混合剪切流动的数值模拟

本文介绍了ＢＧＫ－Ｂｏｌｔｚｍａｎｎ格式的构造与用有限体积法计算网络内宏观量的方法。文中还介绍了可压缩平面混合剪切层计算域在不同出口马赫数下的边界条件。然后将上述方法应用于四种不同情况的可压缩平面混合剪切层的数值模拟。还详细介绍了计算结果,并与其他方法所得的的结果以及实验数据进行比较。表明本文的结果较好,方法是有效的。     

Fourier伪谱方法在不可压缩平板边界层研究中的应用

以三维不可压缩平板边界层为研究对象,扰动形式N-S方程为控制方程,从空间模式的角度,直接数值模拟了三维不稳定T-S波传播的过程.时间离散采用三阶精度混合显隐分裂格式,空间离散则结合Fourier伪谱方法及高精度紧致有限差分逼近,法向采用非等间距网格坐标变换,出口边界条件采用嵌边函数法,程序采用MPI(Message passing interface)并行方法编写.实例验证,该方法计算结果与流动稳定性分析的结果一致.      

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

两种雾化模型在横向流雾化数值模拟中的应用

通过对横向气流中燃料雾化的数值模拟,分析了两种常用的雾化模型TAB(Taylor analogy breakup)模型和Reitz波不稳定性模型中的经验参数变化对模拟结果的影响,为合理使用这两种雾化模型提供了参考.对雾化过程的模拟结果与实验结果的对比分析发现:虽然TAB和Reitz波模型都能对穿透深度进行很好的预估,但从液滴分布、SMD(Sauter mean diameter)和液滴速度等具体参数的计算结果来看,TAB模型的计算结果优于Reitz波模型,更适用于相关计算,并给出了适当的TAB模型的经验参数取值.最后,用TAB模型在相同经验参数下对不同燃料喷射速度的工况进行了计算,得到了满意的结果.这种先用实验数据校准经验参数,然后用校准后的经验参数计算相关未知结果工况的方法可以应用在今后的预估计算和研究中.      

用于可压缩壁面湍流LES及DNS入口边界条件的"回收-调节"方法综述

为了减小湍流流动问题的大涡模拟和直接数值模拟的入口湍流脉动边界条件对计算结果的影响,叙述了用于可压缩壁面湍流入口脉动生成的一种方法--"回收-调节"方法(recycling-rescaling method),叙述了这种方法的理论基础及简要推导过程,介绍了几种常用的"回收-调节"方法并指出了实现中的几个关键问题,分析了...     

基于数值模拟的流场附面层边缘识别方法

附面层边缘通常取在速度达到主流速度0.99倍的位置,而复杂流场中主流流动往往并不均匀,给附面层边缘的准确识别造成了困难。为解决此问题,提出了用"参考主流"代替实际主流识别附面层边缘的方法:通过零剪切力滑移壁面边界条件下数值模拟得到不受附面层干扰的参考主流,在根据附面层定义确定附面层边缘时以该参考主流中的速度代替实际的主流速度。通过斜楔压缩和弯曲压缩两个超声速压缩流场对该识别方法进行了验证,所得到的斜楔压缩出口截面上附面层厚度与采用实际主流速度判断得到的厚度相对误差仅4.1%。根据该方法的识别结果对弯曲压缩型面设计进行附面层修正后,弯曲激波高度与无黏设计值之间的误差从修正前的2.0%降低至0.3%,压缩面末端压力的相对误差从修正前的6.6%降低至2.3%。该方法避免了指定主流速度的主观性,识别结果较为准确。     

BL模式在超声速湍流边界层中应用的检验

湍流模型是目前解决工程中湍流问题不可缺少的工具.然而它终究是建立在简化与假设的基础上的,其准确性必须经过实际的检验.对超声速边界层,由于缺乏系统的实验数据,检验工作很难做.本文以BL模型为例,通过可压缩湍流边界层的三个典型算例,与直接数值模拟(DNS)的结果作对比.分析了其成功及不足之处,以及影响其精度的若干因素,并检验了用BL模型预测转捩位置的有效性.     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

进口畸变下非定常数值模拟进口边界数学模型

进口畸变下压气机内存在的不同空间和时间尺度的非定常流动之间的相互作用是进口畸变对压气机性能影响的重要流动机理.介绍了一种用于周向进口总压畸变下单转子通道内非定常数值模拟的进口边界简化模型,并基于该模型进行了三维非定常数值模拟.数值模拟结果显示了进口周向总压畸变的大尺度总压脉动和压气机转子通道非定常分离流场的相互作用过程.该模型为周向进口畸变的三维非定常数值模拟提出了一种简化方法.      

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动.该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递.分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法.流场...