气动加热下飞行器油箱温度场耦合分析

气动加热与结构热传递耦合问题在航空航天领域非常重要.分析了飞行器持续受热过程中的耦合性,指出了耦合性在分析持续气动加热问题中的必要性.针对存在的耦合性,提出了一种处理方法,即在“飞机结构三维温度场分析”程序中增加气动热流工程算法,同时采用对流传热来模拟油箱内燃油与油箱内壁面传热.对流换热系数以及出口油温通过与油箱壁面结构温度耦合计算来确定,即在每个时间点上,将当前结构内壁面温度咒代入对流换热系数α以及油温L的计算公式,求得该时间点的α和t,以此实现气动/温度场/流场三个物理场的实时耦合.计算结果表明,耦合计算结果与试验结果吻合得很好.     

类乘波体防热壁板气动加热-温度场耦合特性分析

针对类乘波体防热壁板,采用松耦合算法推进气动热和传热迭代计算,研究气动加热-温度场的耦合特性.结果表明:在壁板达到辐射平衡前,冷壁热流和辐射平衡热流与真实气动加热的误差分别达+55.1%和-15.4%以上,必须将气动加热和温度场进行耦合分析;当地绝热壁面温度不随时间变化,表面传热系数是耦合效应的关键参数;采用平均表面传热系数进行瞬态气动加热-温度场耦合计算只进行2次气动加热计算,壁面温度预测误差在2.5%以内,可有效提高气动加热-温度场耦合计算的效率.      

钝体外形气动加热与结构传热一体化数值模拟

为了提高气动加热与热防护结构传热多物理场数值模拟的稳态计算效率与计算精度,发展了一种基于有限体积法的气动加热与结构传热一体化数值计算方法。该方法将高速流场与结构温度场统一到同一物理场,基于统一的控制方程组,采用基于LU-SGS隐式时间迭代和自适应时间步长的有限体积方法进行求解,避开了传统气动加热与结构传热耦合求解方法在时间域内的所需繁琐数据交替迭代策略。对二维/三维钝体进行一体化数值计算分析,计算结果表明:二维钝体非稳态下,得到2s时圆管驻点温度最高达到390.2K,驻点热流密度和结构温度与参考文献和实验值吻合较好,证明了方法的可靠性和可行性。同时分析了三维钝体应用算例的流-固-热稳态计算特征,计算得到稳态时钝头体结构外壁表面最高温度达到535.6K,表明一体化计算方法可用于长航时飞行条件下的气动加热-结构传热多物理场耦合计算分析,为高速飞行器热防护结构设计与选材提供一定的理论与技术支持。     

流-固耦合计算的应用研究

为了解决用常规算法求出的流体与固体对流换热系数的不真实性问题,采用了流固耦合计算方法计算固体壁面的温度.其中采用模型模拟湍流流动,用壁面函数法修正处于流场内部的固壁,通过固体和流体双向耦合换热计算,得出了整个流场、温度场包括(固体部分和流体部分)的分布.     

不同热边界对超音速机翼气动热数值分析的影响

飞行器在超音速飞行时受到的气动加热效应给结构强度及热防护设计带来极大影响，且真实状态下的气动热环境需要考虑外流场与结构的耦合及内壁面边界条件的影响。采用S-A湍流模型求解Navier-Stokes方程，通过流场与固体壁面交界处的信息传递，实现外流场与结构场的耦合数值分析。针对三种不同翼型的超音速绕流气动加热进行耦合数值研究，对比翼型内壁面在不同热边界条件下的气动热效应，结果表明：不同翼型具有与气动力相似的的气动热效应；内壁面考虑对流换热的边界条件最接近真实；考虑机翼燃油箱满油时，三维机翼前缘驻点处热流密度最高可达4200w/m2。     

有限催化对返回舱气动热环境影响

基于给定催化效率的方法,发展了包含离子组分的有限催化边界条件,并结合多组分化学非平衡N-S方程数值求解,建立了有限催化条件下的高超声速飞行器气动热环境计算方法。采用不同催化效率对返回舱外形典型再入工况气动热环境开展了计算分析,研究了壁面有限催化对该返回舱气动热环境的影响规律。结果表明:壁面催化效率对返回舱气动加热影响显著,采用低催化效率壁面材料可有效缓和返回舱气动热环境;扩散热流相对传导热流对壁面催化效率更加敏感,是影响气动热的主要机制,但热流并不随催化效率增加而线性增大;壁面有限催化对气动热的影响不仅与壁面材料催化效率有关,也与流场离解电离程度、壁面密度、温度等当地流动参数相关。     

双工况氢氧发动机燃烧与传热数值分析

应用三维湍流N-S方程以及颗粒轨道模型描述双工况氢氧发动机内部喷雾两相燃烧流动过程。两相之间的质量、能量交换由液滴蒸发模型计算，气相化学反应速率由Arrhenius公式计算。通过耦合求解气液两相的模型方程，对发动机转工况前后的三维流场进行了数值计算，并耦合计算了燃气与壁面之间的传热以确定壁面的温度和热流分布。另外还对分别采用同轴离心式喷嘴和直流式喷嘴得到的燃烧流场与燃烧效率进行了比较。计算结果表明转工况前的壁面温度与热流都比转工况后大。离心式喷嘴的雾化混合效果与燃烧效率都比直流式喷嘴好。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

高超声速流场与结构温度场一体化计算方法

在对国内外流场与结构温度场一体化计算方法的不足进行细致讨论的基础上,提出了一种高超声速流场与结构温度场一体化计算方法.采用统一的积分方程组作为气动加热和结构传热物理过程的控制方程,对整个物理场进行统一的迎风格式有限体积方法离散,给出了流场与结构交界面上温度、温度梯度及导热系数等参数的计算方法.在时间推进方面,定常状态采用多步龙格库塔迭代格式,非定常状态则采用双时间步长方法.采用发展的一体化算法对二维圆管模型的气动加热和结构传热问题进行了数值模拟.结果表明:2s时驻点物面温度为390K,与其他文献的误差在3.1K范围内;稳定时驻点物面温度为647K.      

考虑导热对流和辐射作用的轴对称收扩喷管壁温计算

基于N-S方程求解了包括引射流、加力燃烧室在内的轴对称收扩喷管内外流一体化的流场,建立了考虑导热、对流换热和辐射换热作用的轴对称收扩喷管各层壁温分布的计算模型,包括隔热屏、喷管简体、外调节片三层结构.对某航空发动机进行了多工况计算,燃气的物性随压力、温度和油气比的变化采用了一系列精度较高的计算公式和计算方法.计算结果表明:基于密度求解器求解包含引射流在内的轴对称收扩喷管内外流一体化的跨声速流场是成功的.隔热屏和收扩喷管筒体沿流向温度逐渐升高,喷管筒体壁温在喉部达到最大,在扩张段逐渐降低;收扩喷管外调节片壁温与收扩喷管简体的壁温变化规律相同,但是壁温最大值则位于喉部前某一位置.计算结果与经过试验验证程序的结果符合良好.      

旋转固体火箭发动机喷管热结构分析

以某固体火箭发动机喷管为例,研究其旋转情况下的流动、传热及热结构响应性能。利用雷诺应力湍流模型,结合增强型壁面函数,求解N-S方程。以喷管内壁面温度分布为边界条件,求解二维轴对称瞬态热传导方程。将瞬态温度场按时间步长加载,进行瞬态静力学分析,得到不同转速下的流场、温度场及应力场。结果表明,喷管在旋转情况下,喷管内壁面尤其是喉部及扩张段温度分布在转速为100 r/min及600 r/min左右时变化剧烈;旋转内流场与喷管结构的耦合作用加剧了喷管的传热,尤其是喉衬的烧蚀,特征点温度值随转速增大而升高;最大热应力位于喷管最外层尾端,整体热应力在转速低于100 r/min时得到释放,随着转速的不断增大,喷管整体最大热应力及扩张段特征点应力随之增大,而喉衬特征点由于旋转导致了温度梯度降低,其应力值随之减小。     

准一维可压缩瞬变管流的有限体积模型(Ⅱ)管壁温度场的有限体积模型

在流场模型的基础上,通过对圆柱坐标系下轴对称管道壁面划分的二维有限体积网格,建立了一种计算管壁瞬变传热的有限体积模型,可处理对流换热和辐射换热两种边界情况,可处理具有包覆层或真空夹层结构的变物性管壁传热.温度场模型是流场建模思想的自然延伸,对二者的组合运用即为所发展的可仿真准一维可压缩流管内瞬变流动的有限体积模型,一方面,在流场仿真的体系内发展了传热计算的部分并最终扩展成为统一的流动/传热仿真体系,另一方面,结合阀芯节流模型,从此模型出发可推导出管路系统常见元件的流场和温度场模型.对某发动机试验台液氧贮箱增压系统的建模与仿真表明,提出的模型体系具有很好的适用范围和良好的仿真精度.      

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

TBCC燃烧室/喷管一体化壁面温度计算

基于Navier Stokes(N-S)方程组对包括隔热屏、隔热屏内外流、大气外流在内的涡轮基组合动力(TBCC)发动机燃烧室/喷管进行了一体化的气/热耦合数值模拟,考虑了燃气组分输运、辐射换热等影响,研究了其在某典型飞行状态下TBCC冲压发动机燃烧室/喷管筒体及隔热屏内外壁壁面温度、辐射换热热流及对流换热热流分布.结果表明:燃烧室/喷管筒体与对称面上下交线的壁面温度在轴向距离为0.5~2.6m内变化较小,在轴向距离为2.6~3.1m内急剧增加,在轴向距离为3.1~3.5m内急剧下降.之后,上交线筒体壁面温度沿流向减小,下交线筒体壁面温度先升高后降低.筒体壁面温度最高点在喷管下调节板收缩段,为1577K.隔热屏内壁面辐射热流在370~500kW/m2变化,上下交线处的辐射热流较外壁面的辐射热流约高300kW/m2,辐射热流沿流向先减小后增加.隔热屏外壁面辐射热流在50~200kW/m2范围内分布.      

客机APU舱温度场的数值计算

为研究大型客机APU舱传热过程和进行结构热分析,采用数值方法计算了APU舱温度场。采用商业CFD软件Fluent,选用Realizable k-ε湍流模型和S2S热辐射模型建立了APU舱内流动和传热数值计算模型,得到了流场、温度场和热流信息。APU舱中主要传热过程是辐射传热;APU防火罩和排气管隔热罩起到了热辐射遮热罩的作用;防火罩内的APU冷却空气对于降低防火罩温度有重要作用;排气管内的引射冷却空气有效降低了壁面温度。     

低温推进剂贮箱增压系统分布参数数值仿真(Ⅰ)贮箱的有限体积模型

从准一维可压缩瞬变管流的有限体积模型出发,通过数值拓展和集成建立了一种考虑流体与箱壁传热的低温推进剂贮箱分布参数模型.流场采用一维交错网格的有限体积模型,可考虑变物性、轴向热传导、重力场、壁面传热影响;壁面温度场采用圆柱坐标系下的轴对称二维有限体积模型,可处理对流换热和辐射换热两种边界情况,可处理具有包覆层或真空夹层结构的变物性管壁传热.模型全面考虑了贮箱内的分布参数效应,适用于贮箱增压计算领域.      

复合喷管热结构耦合计算的一种新的策略

基于流固耦合完成复合喷管流动换热、传热过程温度场的数值模拟,编写了计算流体力学(CFD)与有限单元法(FEM)的接口程序,快速实现温度载荷和机械载荷数据的传递,提出了复合喷管热结构耦合计算一种新的策略。利用直接约束法考虑了复合喷管部件界面间高温碳化带来的结构边界非线性因素。计算结果表明耦合方法得到的复合喷管温度场是合理的,高温下界面间结构功能退化问题对于变形特征以及应力场影响很显著。     

结构形式对双层微通道热沉传热性能的影响

设计一种双层微通道热沉,即热沉上下层通道结构相异B（上层为矩形微通道,下层为复杂结构微通道）,旨在降低压降的同时保持良好的换热性能,并与上下层通道均为复杂结构微通道的热沉A作对比。用数值方法研究了顺流和逆流时、不同通道结构（PA、PB、CA、CB）对其传热性能的影响,并用场协同原理分析流场与温度场的协同关系对传热的影响。结果表明:当0相似文献   

摆盘流动与换热性能数值仿真

针对受燃烧室出口高温燃气冲刷的摆盘装置的冷却水回路,建立二维、三维物理仿真模型,模拟了冷却水在内部流道的 流动与换热过程,考察了冷却水进口压力、冷却水流量对流阻和换热性能的影响.结果表明:①摆盘冷却水进口压力由2.3×105Pa提高到8.3×105Pa,摆盘壁温变化微小,结构1壁温升高5K,结构2壁温升高7K;②冷却水进口速度由0.5m/s提高到5.3m/s,结构1壁温降低约120K,结构2壁温降低约100K,冷却效果明显;③结构2通过缩小流道的流通面积,能在更小的冷却水流量的工况下得到更好的换热效果.