类前缘防热层流场和热响应耦合计算研究.pdfVIP

中国工程热物理学会 传热传质学 第十届 年会 编号：013164 类前缘防热层流场与热响应耦合计算研究+ 桂业伟袁湘江 (中国空气动力研究与茇展中心．四川绵阳621000) (联系电话：0816—2466025 e—mail：gywy神etang．com) 摘 要本文在以往对类前缘防热层的热响应计算、分析的基础上，进一步研究实现了外流 场高超音速Ns方程数值计算表面气动加热与防热层结构热响应的耦合计算．这对于研究和工 程应用中常用的非耦合计算方法来说是～大进步，也为进一步开展外流场／结构热响应／热 应力全耦台一体化计算研究和防热层表面畎气强化传热问题的场协同研究打下了基础。 关键词热结构NS方程热响应数值模拟 1引言 目前，世界航天发达国家都在致力发展完全可重复使用运载器，其中关键是可重 复使用热结构。由于相似性和成本等原因，在先期研究中常采用计算机模拟，以弄清结 构在强烈气动加热下的热响应情况“1】。这是典型的高超音速气动加热与结构热响应耦 合问题，国际上目前在沿再入轨道研究此类问题中常采用松耦合方法，即在几个给定轨 道点上进行迭代，以确定气动热与表面温度的变化关系；在工程计算中．通常沿轨道给 定表面冷壁热流，实际计算中按表面温度再作热壁修正。事实上，由于再入过程中强烈 的气动加热，使再入体表面温度发生迅速的变化，这会对表面气动加热带来影响，采用 前述的松耦合计算就存在较大偏差。本文正是针对前缘防热层进行了Ns方程与结构热 响应全耦合计算，其中对防热层热响应计算引入了计算流体力学常用的任意坐标变换和 代数网格“]，用差分法成功地进行了计算：外流场则求解完全NS方程得到。通过耦合 计算，本文探索了前缘防热层外流场气动热和结构热响应在耦合条件下的特征，为下一 步开展外流场／结构热响应／结构热应力的全耦台一体化计算研究打下了基础，同时也 为进一步分析研究在表面带有局部吹气时的流场和温度场的协同特征打下了基础。 2模型与方法 根据前缘类热结构特征和模型示意图‘”，由于前缘类结构在展向的变化很小，因此 可以按二维问题进行处理嘲。本文所研究的对结构热防护起关键作用的防热层的形状如 ’国家自然科学基金项目 国家重大基础研究(973)项目：G2000026301 图1所示，防热层总厚度取4cm。 按文献[4]引入任意坐标交换，把复杂物理平面变换成矩形的计算平面。变换后的 热传导方程为 譬警=昙kL曼阮习+量妻恢砷+毒晏旧砷+与妻旧叫 +善l毒垤砷+靠妻恢砷+旬妻旧7j+。毒％叫m 式中下标表示偏导数，具体表达式详见文献[4]，，为坐标变换3acobian行列式，p、 c、^分别为常物性的密度、比熟和导热系数。在文献[4]基础上，生成贴体代数网格， 考虑到加热面附近温度变化大，因此对加热面附近法向网格进行加密。本文按照 Frankel—Duffoft格式对方程进行离散化，离散后方程形式见文献[4]，计算格式具有 时间和空问二阶精度．适于非定常计算。 对于外流场计算，引入无困次变换口】，则可得二维守恒型无因次NS方程： 掣+娑+竽+孕+譬：0 (3) at 瓠 av 。]c av 一 式中U，F、G等各物理量的表达式详见文献【5】。 外流场计算的边界条件如下：在网格外缘边界，为来流条件：在类前缘物形的后部 出口边界，采用零梯度条件(即物理量沿x方向的梯度．沿Y方向的梯度分别为零)； 在固壁边界，“=v=0，壁面温度与结构热响应计算耦合给出。 将方程(3)变换到一般曲线坐标系(具体见文献[5】)．再利用通量分裂技术 F=F++F一，G=G++G一(上标“+、一”分别表示正负通量)．方程可写成： au+翌+笠+笪+堑+盟+亟：o (4) af a善 刁善 a可 口叩 a善 a，7 将方程(4)分裂为如下两下“一维”方程