高速旋转再入飞行器防热层热滞后数值计算.pdfVIP

第十届垒国计算流体力学会议论文集 高速旋转再入飞行器防热层热滞后数值计算 陈兰 胡汉东 中国空气动力研究与发展中心，匹川绵阳621000 摘要 本文利用数值方法研究了高速飞行器有攻角旋转再^时出现防热层热滞后现象。该方法在字间 离散上采用分块的有限体积法，时间推进采用四阶龙椿一库塔积分．实践表明该数值计算稳定，时间精度高． 町吐报好地模拟复杂防热结构瞬态温度场非定常热响应过程。 O引 言 当高速飞行器以一定速率旋转作有攻角再八，其防热层表面每一点所受到的气动加热各 不相同．在迎风面上气动加热量大，而转到背风面后，其所受气动加热又降低，因此，表面 热流边界是随时间作连续周期性变化的。这一气动加热的特点给飞行器表面防热层的温度分 布规律带来新的特征：防热层表面温度的最大值与表面气动加热量的最大值并不“同步”，从 E行器周向来看，随着飞行器的转动．表面温度的变化相对于表面热流状态的变化有一个相 位滞后，这就是由于飞行器旋转而出现的热滞后现象。滞后角的出现会对飞行器的烧蚀防热 户生影响，从而影响飞行器的气动力特性，因此，弄清热滞后现象在飞行器设计过程中很有 2：。 意义lIf 本文提出了一种求解复杂防热结构瞬态温度场的数值方法，用有限体积法解二维热传导 方程，时间推进采用四阶龙格一库塔积分。通过坐标变换和网格分块技术，使边界的处理更加 精确．可以方便地应用到各种复杂形状的热传导问题中。典型算例的计算表明该方法数值计 算稳定，时问精度高，可以很好地模拟非定常热响应过程。同时，本文也给出了再入龟行器 防热层滞后角沿再入轨道的变化曲线。 1数值方法 1．1控制方程 设有温度场r(u，f)=r(x，y，=，r)，场中某一点M，任取包含点M内部的一闭曲面S， 其所围区域为r，i为s的外法线方向上的单位矢量。在闭曲面S上任取面积微元dS，由 热量守恒定律．单位时间内通过曲面S流入区域V的全部热量等于区域y内各点温度从 ，kY，‘，)变化到，0，Y，而，+At)升温所需的热量： dv=§砜T·hdS 聃lpCPi01 一392— 公式中的p、CP、k分别表示密度、比热和热传导系数。 蚂争2炒阶R喇+e蒯q Ev=J毽；Ev+；yFv+{zG。、 ‘=J(仉巩+巩鼻+仉q) G。=J《；E。+‘vFv+‘zG，、 于：刀 这里，J为雅可比行列式，=制 E=t豢，u=t筹Gv=‘iOT 该控制方程可采用有限体积方法离散。 1．2时间推进的四阶龙格一库塔法 将控制方程记为 iOT：Rh4r) 利用标准四阶龙格一库塔递推公式有 at(：q+2止+2兀+^) Tn-I=r”+-。1 T”、 ft=Rhs(t f2=冗bk +％出，，”+％护^) f3=兄bk +％△f’r+％出·，2) +At，T”+ f4=Rbb △f+^J 高阶龙格一库塔算法精确度高，绝对稳定性区域较大，不过每一时间步需多次计算右端项。 1．3初始和边界条件 初始时刻t=0时，给出求解域内各点温度分布情况孔一。=“M) 若空间区域V的边界记为F，阻P表示边界r上的动点，考虑了第一、二类两类边界，即