气冷涡轮气热弹耦合有限差分算法地研究.pdf

Classified Index: TK471 U.D.C.: 621.438 Dissertation for the Doctoral Degree in Engineering RESEACH ON METHODOLOGY OF THERMAL-FLOW-ELASTIC COUPING NUMERIAL SIMULATION IN AIR-COOLED TURBINE WITH THE FINITE DIFFERENCE Candidate ： Guo Zhaoyuan Supervisor ： Prof. Feng Guotai Associate Supervisor ： Prof. Wang Songtao Academic Degree Applied for ： Doctor of Engineering Speciality： Power Machine and Engineering Affiliation ： Energy Science and Engineering Date of Defence ： Oct., 2009 Degree-Conferring-Institution ： Harbin Institute of Technology 摘 要 摘 要 航空工业的迅速发展对航空发动机提出了更高的要求，即高的推重比和热效 率，导致涡轮的进口温度不断提高，目前该温度已远远高于涡轮材料的耐热极限， 这就需要采用复杂的冷却技术来保证涡轮叶片安全的工作和较长的工作寿命。准 确预测叶片内热传导以及热应力的分布是提高涡轮冷却效率、保证叶片工作寿命 的关键性问题，气热弹数值仿真技术是解决该问题的有效方法和工具之一。本文 主要对应用有限差分方法进行气冷涡轮气热弹数值仿真的关键性问题进行研究， 并对气冷涡轮内气热弹耦合现象进行初步探讨。 首先通过深入分析气冷涡轮内复杂的多场耦合现象，研究了气热弹多场耦合 的物理模型，给出了以柱坐标系下参数为不变量的任意曲线坐标系下的无量纲时 间平均雷诺N-S 方程，给出了直角坐标系中考虑热传导、热辐射以及热应变对能 量耗散影响的温度场方程并将其在任意曲线坐标系中展开；推导了直角坐标系和 任意曲线系中考虑热变形影响以及不考虑热变形影响的以位移为求解量的弹性固 体应力场平衡微分方程的表达式；讨论了气热、气弹、热弹、气热弹等四类多场 耦合边界条件的给定方式，从而建立了多场耦合的计算模型。由于采用了多块结 构化网格离散气冷涡轮气热弹耦合的计算域，又给出了四类不同耦合计算网格块 之间的数值传递方法，实现了多场耦合网格之间的数据传递。 然后研究了流场、温度场、弹性固体应力场求解的差分格式及其数值方法， 这是应用有限差分法进行涡轮气热弹多场耦合计算中的重要问题。三维粘性流场 控制方程采用具有 Godunov 性质的三阶 TVD 格式离散，并采用 AF 方法求解。温 度场采用了稳态和非稳态的两种求解方式：稳态温度场采用隐式 ADI 交替方向方 法求解；对于非稳态温度场，构造紧致格式离散一次和二次交叉偏微分项，并应 用高阶格式的 ADI 法求解离散后的方程，该方法具有时间方向一阶精度、空间方 向四阶精度，经验证表明，该方法较显格式和 C-N 格式具有更高的求解精度。引 入位移法求解二维弹性应力方程，为了验证该方法的可行性和可靠性，对具有解 析解的悬臂梁进行受力分析，结果表明：第二类边界条件的求解精度对整个计算 的结果有重大的影响，平衡方程二阶精度的计算值与解析解吻合程度好，而一阶 精度的计算值偏离解析解较大。然后对三维弹性固体应力场控制方程进行了无量 纲化。在 Possion 模型方程的基础上推导了交替方向求解方式，分别构造了空间 方向三点（二阶精度）及五点（四阶精度）的差分格式，在求解中，有变量分离 式以及耦合式两种方法供选择。计算表明：高阶格式具有较高的求解精度；变量