分离流中旋涡运动的热行为及其控制的数值模拟

以压缩拐角流动为例，对分离流中旋涡运动的热行为及其控制进行了数值模拟研究。研究结果表明，光滑表面上流动的分离特性及其热行为受到壁面温度的极大影响，提高壁面温度可以使分离点前移，再附点后移，分离范围扩大，涡心提高，而再附点附近及峰值热流则经历了一个先上升后下降的非单调变化过程；分离区和现附点附近的温度型与附着流差别明显，具有很高的峰值和很靠近物面的峰值点，从而造成高热流。而且温度梯度随壁温的不敏感性     

尖楔结构低速高/中温双路气流组合热试验方法

在典型高超声速飞行工况下,数值模拟分析了高温合金尖楔前缘结构沿气流方向大温度分布梯度将带来严重的热强度问题,产生大温度分布梯度的根本原因是尖楔结构头部区域平均热流密度与后段平板区域平均热流密度之差,而受头部区域热流密度具体分布的影响不大;进而提出了一种低速高/中温双路气流组合热试验方案,并通过数值模拟方法证明了该方案具有两股气流参数可以独立调节分别满足尖楔结构头部驻点区域及后段平板区域大、小两种热流密度的优点,进而解决单喷口低速高温燃气流热模拟试验难题,满足尖楔结构高超声速飞行工况下大温度梯度模拟要求.同时,该方案通过高/中温气源的合理组合搭配可以大大降低尖楔结构热试验所需高温气源发生功率,推广应用于电弧风洞可拓展其热试验范围.      

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

复杂外形飞行器热流的NND有限元数值计算方法

有限元方法是一种加权积分的方法,其特点和权函数的特性使得有限元方法对于复杂外形的边界处理和计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文采用张涵信的NND思想发展了复杂流场的NND有限元计算方法,算例表明,在相同的网格雷诺数下,采用积分的方法提高了物体表面热流的计算精度,对高超声速复杂外形飞行器的绕流获得了良好的数值模拟结果.     

跨音速压气机在非设计工况下动静叶相互作用的数值模拟

本文计算了某跨音速压气机的一级在非设计工况下动静叶干扰产生的非定常流动。计算以准三维的Navier-Stokes方程为模型,湍流模式采用Baldwin-Lomax代数模式。主控方程用高分辨率的MUSCL型Roe格式离散,离散方程用高效的LU-SSOR方法求解。本文提出了一种完全守恒的动、静叶交界面处理方法。周期性边界处理采用Rai的方案。根据计算结果,讨论了非设计工况下动静叶干扰的机理。      

带副翼三维机翼绕流的N-S方程解

提出了一种求解带副翼偏转三维机翼绕流的N-S方程计算方法。采用对接分区网格与分区求解算法的结合,有效地求解绕此外形的复杂流动。提出了一种满足通量守恒的内边界耦合条件。数值方法中选用vanLeer分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD性质,采用中心差分格式来离散粘性通量项。数值算例表明该方法是求解带操纵面偏转的机翼绕流的有效方法。     

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。      

防热试验用高温超声速燃气流场热环境分析

对用于防热系统试验考核的高温超声速燃气流场的流场结构和热环境特性进行了计算和分析.将CFD模拟方法与工程计算方法相结合,对两类锥形喷管形成的燃气流场进行了数值模拟计算和流场结构分析,并与试验过程中高速拍摄的图片进行了比对.利用模拟得到的流场气动参数,计算了沿流场轴向位置球头驻点模型的热流和压力,并与试验测量数据进行了对比,结果吻合较好.锥形喷管形成压缩波和膨胀波交替的流场结构,导致流场气动参数与模型热流和压力形成正弦变化的特点.     

低速高温燃气流热模拟试验方法和设备

对比分析了两种气流状态参数和两种加热情况下典型前缘部件表面热流密度的相似性,论证了利用亚声速高温燃气流加热方式进行近地空间高超声速飞行工况气动热模拟试验的可行性.针对高超声速飞行器典型钝头锥结构提出“小喷口低速高温燃气流+石英灯”组合热试验方案.通过采用新型高效双腔蒸发管型燃气发生器、新型带保温夹层和耐高温陶瓷内衬的水冷不锈钢高温管道结构,同时引入电加热器预热及燃烧室两路供油方案,使所建低速高温燃气流热试验设备产生燃气流温度达到2100K,φ250mm喷口处平均径向温度分布梯度约3K/mm,具有线性温度控制功能且稳态控制温差约46K,满足24km、马赫数为6典型高超声速飞行器工况驻点区域高温/大热流密度气动热试验要求.      

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

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

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

对流传热场协同原理在高速可压缩边界层流动中的推广

应用理论分析方法对适用于不可压缩层流与湍流流动的对流传热场协同原理进行了可压缩层流与湍流流动上的推广。分析结果表明,与不可压缩流动的对流传热场协同原理不同,可压缩层流与湍流的对流传热取决于流动当地单位体积的动量与总焓梯度的协同。用当地单位体积的动量与总焓梯度的协同研究可压缩流动的壁面传热问题,对层流热流,不但计及了高速流动的密度变化对热流的作用,而且包括了静焓梯度、压力梯度、壁面分子黏性剪切效应对热流的影响;对湍流问题,除了高速流动的密度变化、压力梯度、壁面分子黏性剪切效应对热流的影响外,还计及了雷诺剪切应力对热流的作用。另外,对黏性影响不能忽略的不可压缩流动的对流传热问题,用速度向量与总温（总焓）梯度协同更精确。     

不同淬熄结构下RQL燃烧室流场特性的数值研究

采用雷诺平均Navier-Stokes(RANS)方法对不同淬熄结构的富油/淬熄/贫油(RQL)燃烧室内定常掺混流动过程进行了数值模拟,分析了不同截面的流场结构,研究了动量通量比及淬熄孔排布方式的影响。结果表明:射流进入RQL燃烧室后发生偏转且偏转方向和旋流方向有关,回流区长度、穿透深度和高湍动能区随着动量通量比的增大而增加;当J≥80时,淬熄射流下游出现高湍动能区,逐渐与上游融为一体;当J≥120时,回流区的长度基本不再发生明显变化。此外,非中心对冲结构C3具有更大的高湍动能区,掺混更均匀。      

新型热流计测量技术的探索研究

目前的热流计在测量热量时需要测量供热管道的流量以及进出口温度,然后再对其进行积分计算.其技术非常复杂,而且精度低,仪器成本高,难以普及.为了解决这一系列问题,本文主要依据能量守恒定律提出了一种新型的能量测量装置,并对该装置的测量原理进行了理论分析.通过数值模拟,在理论上对其进行了验证.数值模拟的结果证实该型热流计测量结果与系统流体的流量无关,这也是本测量方法不同于其它常规热流计的优点之一,摆脱了流量对能量测量的影响,同时也简化了计算过程,大大降低了成本.     

轴对称矢量喷管内流特性的数值模拟研究

利用矢通量分裂有限体积格式和B-L湍流模型,采用TTM网格生成和局部网格加密技术,求解三维N-S方程,对轴对称矢量喷管的内流流场进行数值模拟,并将计算结果与模型试验结果进行了分析比较。研究结果表明,在矢量状态下,计算结果和试验结果相比,壁面静压分布的相对误差不大于13%,气动矢量角的相对误差不大于10%,推力系数的相对误差不大于2%,所开发的计算方法和计算程序可用于轴对称矢量喷管的工程设计。      

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

超临界压力正癸烷在竖直圆管内的流动阻力实验

研究了加热条件下超临界压力下正癸烷在竖直细圆管内的流动阻力特性,分析了热流密度、系统压力以及进口温度对摩擦压降的影响规律。实验结果表明:当流体温度较低时,流体摩擦压降随着热流密度的升高而逐渐降低;当流体温度逐渐接近拟临界温度时,由于物性的剧烈变化,流体摩擦阻力压降随着热流密度的升高而增加。并且系统压力越低,拟临界点附近的变化就越剧烈。当流体温度高于拟临界温度时,物性变化趋于平稳,流体摩擦压降随热流密度升高而增加。基于实验结果,总结出了一种针对不同温度、不同热流密度与质量流速比值范围的摩擦阻力系数经验关系式。      

稀薄流航天器鼻锥迎风凹腔气动力和气动热性能研究

本文采用直接模拟蒙特卡罗(DSMC)方法,对高超声速稀薄流中航天器鼻锥迎风凹腔气动力与气动热性能进行了数值研究。得到了鼻锥外壁面、凹腔侧壁面以及凹腔底面的热流密度分布,分析了不同凹腔深宽比对鼻锥冷却效率以及凹腔腔体内气体参数的影响;以深宽比为1的凹腔为基准,研究了凹腔唇口钝化半径对航天器气动热与气动力的影响。数值结果表明,高超声速稀薄流中迎风凹腔能够降低鼻锥外壁面的热流密度;当凹腔深宽比达到1之后,凹腔侧壁面热流变化趋于一致,热流密度最低点的轴向位置不随深宽比改变,且凹腔底部热流很小;凹腔近底部气体均由稀薄流转化为连续流,腔内气体压力不断振荡;唇口钝化没有明显优势,虽然可以降低鼻锥峰值热流,但是会带来严重的气动力性能下降。     

锯齿扁管内沸腾换热试验

建立以等热流密度方式进行试验件加热的沸腾换热试验系统,分别对当量直径为1.28mm和1.59mm锯齿扁管内R134a工质的沸腾换热特性进行研究,试验参数范围:制冷剂质量流率为68.5~305.5kg/(m2·s),工作饱和压力为0.27~0.46 MPa,加热热流密度为9~42kW/m2。试验结果表明:相同结构的通道,当量直径小换热能力更强;热流密度和饱和压力对沸腾换热的影响与一个干度值有关。当干度小于此值时,沸腾换热系数会随着热流密度及饱和压力增大而增大;而当干度大于此值时,沸腾换热系数随着干度增大而急剧下降,热流密度和饱和压力对换热的影响较小;该干度值会随着热流密度或饱和压力增大而逐渐变小。质量流率对沸腾换热的影响与热流密度有关,随着热流密度增大,质量流率的影响趋向大干度区域。通过分析各参数对沸腾换热的影响,建立了一个预测试验工况下微小尺寸锯齿扁管的沸腾换热系数计算经验公式。