二维非结构网格有限体积热响应计算的研究.pdfVIP

中国工程热物理学会 传热传质学 1 2003年学术会议 编号：03301 二维非结构网格有限体积热响应计算研究 贺立新桂业伟王安龄杨永键 (中国空气动力研究与发展中心 四川绵阳621000) wal rE—mail tel lh990601@)aotmail．com08162463267) 摘要 本文采用非结构网格和有限体积方法计算了二维热结构的热响应问题，如控制 舵和弧形翼截面(二维问题)在气动加热条件下内部热响应，并与采用有限差分方法计算 结果进行比较。计算结果表明该方法具有较高的精度和网格收敛性，易于处理复杂形状问 题。本文主要给出了有关的模型、计算方法、计算结果和相应的讨论分析。 关键词 非结构网格．有限体积方法，热结构，热响应 0引言 在热防护问题中，结构的热响应是一个最基本的物理现象，热响应特征对结构的防 热性能、抗破坏性能、流场气动加热情况等都有着耦合的作用和影响。因此，弄清结构 的热响应特征，对于飞行器设计来说是十分重要的。由于试验成本、耦合的物理现象的 相似性等原因，在飞行器型号先期研究中，采用计算机数值模拟的方法。弄清结构在飞 行器再入过程中受到强烈气动加热条件下的热响应和结构应力情况，弄清结构的热响应 特征、与流场气动熟的耦合影响规律等，是一项投入较少，同时又很有意义的工作，国 际上已经开展了很多工作“。1。 对于弧形翼外形类问题。由于该计算模型需要多种坐标系，在半园平板内，需用柱 坐标系，在矩形平板内，需用直角坐标系。计算该类模型温度场时，如果直接采用有限 差分法进行数值模拟，存在困难。首先，在计算矩形平板与半园平板的交界处节点温度 时，由于涉及到柱坐标系和直角坐标系下的节点，使计算复杂；其次，在矩形平板和半 园平板区内及它们的边界，需要采用好几套有限差分计算公式。针对以上问题，我们采 用非结构网格有限体积思想来构造计算格式，使得上述问题得以很好地解决。该方法不 仅适用于本模型，而且对于类似的问题具有通用性。 1 二维温度场有限体积计算方法 1．1控制方程 孵詈=尝(t罢)+昙(七号] ㈨ 1．2有限体积计算格式 (1)扩散项的离散 对(1)式在控制单元内进行积分可以得到： ∥，警s=炯彻=扣∥和‘ (2) 其中S为单元面积，n为单元边界的外法向，N为单元边总数，，^为单元k边的长 度，磊为单元k边的热流。 (2)时间项的离散 时间方向采用二阶TVD-Runge—Kutta方法“1进行离散可以得到如下形式 RHS,=∑(一嚣)·碗 (3) 正1；王一脚JRHS?fLpO ∥=j1而ha(互”+互1+删?) 这里i为单元序号，n表示n时刻。 (3)单元边的温度梯度计算方法 单元边上的温度梯度由在包围该单元边的面积内应用高斯定理得到，即对任意函数 f(x．y)I昏 兽：丢弘，凼 (4) OX O。 S为由所计算的单元边相邻的两个单元的中心分别连接该单元边的顶点所围成的四 边形的面积，m是J方向的外法向。同样可以得蛩]OWOy,由此计算得到单元边中心点 的a，／Ox．a∥a^计算单元边上的温度梯度用到了单元顶点值，该值由共享该顶