Carry HE

handouts(3) - heat transfer

· Carry HE

导热数值解法基础

  1. 导热数值解法的思想(简答题常考)

    1. 将时间上和空间上连续的物理量用离散的点的集合来代替

    2. 将微分和微商用差分和差商代替

    3. 按照一定的方法(泰勒展开,控制容积)去建立节点(离散方程)

    4. 用一定的方法去解结点方程

  2. 导热数值解法的优势(为什么用数值解不用解析解)

    1. 物体的几何条件复杂(分析法解不出来)

    2. 物体的传热过程复杂,比如说含相变过程(分析法可能解不出来)

    3. 物体的物性参数不是常数,比如说是温度的函数(分析法可能解不出来)

    4. 物体边界条件复杂,比如说边界条件是对流传热和辐射的复合形式(分析法可能解不出来)

  3. 建立离散方程的两大方法(不列方程)

    1. 泰勒级数法(高等数学和考研数学的内容,具体见泰勒展开)

    2. 控制容积法(列方程都是用控制容积法,因为控制容积法适用范围更加广,物理意义更加清晰)

  4. 非稳态导热的数值解中比稳态导热中多了一个非稳态项,扩散项的离散写法两者是一致的

    1. 时间项(非稳态项)向前差分$\frac{\partial t}{\partial \tau}|_{i,k}=\frac{t^{k+1}-t^{k}}{\triangle \tau}$

    2. 时间项(非稳态项)向后差分$\frac{\partial t}{\partial \tau}|_{i,k}=\frac{t^k - t^{k-1}}{\triangle \tau}$

    3. 白话:向前差分就是想求$k$ 时刻的温度,我用未来时刻的温度来差分;向后差分就是用过去时刻的温度来差分

    4. 时间项不可以用中心差分(这个在泰勒展开中有详细解释,具体可以参考建工传热学对应教材),因为会导致不稳定的结果

    5. (简答题)向前差分又称显式格式(pros计算简单,不需要知道当前$k$ 时刻的所有节点温度,但是也可以知道目标节点$k+1$ 时刻的温度),cons步长的选择会存在限制(后面会进一步解释)

    6. (简答题)向后差分又称隐式格式(pros步长选择没有限制,但是计算会比较复杂,需要知道$k$ 时刻所有节点的温度,才能求出$k+1$ 时刻目标节点的温度)

    7. 步长可以分为距离步长和时间步长,时间步长仅仅出现在非稳态的数值数值计算中

    8. 向前差分的步长限制(一维内部节点):$Fo\leqslant \frac{1}{2}$

    9. 向前差分的步长限制(二维内部节点):$Fo\leqslant \frac{1}{4}$

    10. 根据上面你可以推到出来三维内部节点的步长限制

    11. 向前差分的步长限制(一维边界节点):$Fo\leqslant \frac{1}{2Bi+2}$

    12. 向前差分的步长限制(二维边界节点):$Fo\leqslant \frac{1}{2Bi+4}$

    13. 相信你能看出规律,推出三维边界节点的步长限制条件

    14. 步长限制的主要原因是,如果在数值解中步长不满足限制条件,计算出来的结果可能会离谱,违背热力学第二定律

brief summary - diagram

graph LR 导热 --> 基础理论 基础理论 --> 区分热传导热对流和热辐射 基础理论 --> 金属和气体导热机理 基础理论 --> 傅里叶定律和导热微分方程 基础理论 --> 单值性条件 导热 --> 稳态导热 稳态导热 --> 一维无限大平壁 一维无限大平壁 --> 复合平壁导热 一维无限大平壁 --> 导热系数随温度变化 稳态导热 --> 圆柱体导热 圆柱体导热 --> 临界热绝缘直径 稳态导热 --> 肋片导热 肋片导热 --> 肋片效率 导热 --> 非稳态导热 非稳态导热 --> 周期性导热 非稳态导热 --> 瞬态导热 瞬态导热 --> 无限大平壁导热 瞬态导热 --> 半无限大平壁导热 瞬态导热 --> 集总参数法 导热 --> 数值解法 数值解法 --> 思想和优势 数值解法 --> 建立离散节点方程 建立离散节点方程 --> 泰勒展开 建立离散节点方程 --> 控制容积 数值解法 --> 稳态和非稳态方程建立 数值解法 --> 时间项差分 时间项差分 --> 显式格式 显式格式 --> 限制条件 显式格式 --> 原因 时间项差分 --> 隐式格式