圆柱坐标系中的分离变量法.ppt

二、含热源项非稳态问题 控制方程： 边界条件： 初始条件： 该问题可以化为： 稳态非齐次问题 课堂练习（书中P146） ? 应用练习——肋片热分析 应用练习——肋片热分析 应用练习——肋片热分析 对于非等截面的三角肋、环肋，沿热流方向肋的横截面是变的，对流换热面积从肋基与肋端的关系也不是线性的。 图3-5 肋的形状 下面以环肋为例，进行分析，建立方程及求解。 环肋轴对称，且肋很薄，所以温度只随 r 变化。 取微元，建立能量平衡方程，见 图(3-5). 导入微元体的热量 = 导出微元体的热量 + 微元体对流损失热流 其中 其中，A为导热面积， S为对流换热 面积, 代入能量平衡方程： 得到： 方程化为 这是零阶修正Bessel方程(书中P91 式(3-51c))。 ? 边界条件： 可求：(1)肋的温度分布；(2) 通过肋 的传热量。 本 章 小 结 （一） 齐次问题 (1) 变量为 (r, τ) 的齐次问题 0 ≤ r ≤ b，见表3-1(P97)； a ≤ r ≤ b，见表3-2(P102)； a ≤ r ∞，见表3-3(P105)； (2) 变量为 (r, z, τ) 的齐次问题； (3) 变量为 (r, ? ,τ)的齐次问题 0 ≤ ? ≤ 2π 的整圆柱; 0 ≤ ? ?0 2π 的部分圆柱; (4) 变量为 (r, ? , z, τ) 的齐次问题 实心圆柱；空心圆柱；部分实 心圆柱；部分空心圆柱。 （二）乘积解：多维齐次问题的解法 （三）多维非齐次问题 同第2章 The End Homework II 习题6—8，10月28日交 重要时间 考试：11月18日上午10：00~12：00 计算大作业：11月11日 三、0 ≤ r ∞ 区域问题 特征方程： 若 F(r) 在区域内有界， 是绝对收敛的，那么，任意函数 F(r) 可以用Bessel函数表示为： 令 , 则上式可写成: (3-26) 例3-2: 0 ≤ r ∞，τ= 0，T(r,τ)= F(r)， 试求：τ 0 时该区域的 T(r,τ) 的表达式。 控制方程： 初始条件： 进行分离变量： (3)，(4)的解分别为 求系数 C(β) 由初始条件式，得到 (6)式与(3-26)式对照比较，得到: 将(8)代入(5)，得到 具体积分见附录III(p714)。 四、a ≤ r ∞ 区域问题 图3-3 a ≤ r ∞的区域 一维轴对称方程： 用(3-27a)问题的解 表示任意函数 F(r)的结果为： 对于在 r = a 处为第一、二、三类 边界条件时的 列于 表(3-3)。 例3-3: 在圆柱坐标系中，a ≤ r ∞ ， τ = 0 ， T(r, τ) = F(r)， 在 r = a 处 T = 0，试求：τ 0 时该区域 的T(r,τ)的表达式。 控制方程： 边界条件： 初始条件： 进行分离变量 时间变量函数的解为 空间变量函数应为满足以下问题的解 解为 由初始条件(3)式，得到： 将(8)式与(3-28)式对照比较，得到系数 C(β)： 将(9)代入(7)，得到： 式中, 和 从表(3-3)第三条查到。 五、部分圆柱问题 图3-4 部分圆柱区域 定义的部分圆柱为0 ≤ ? ≤ ?0 2π, 该问题的基本解为 和 其中，特征函数 模 和特征值 ，可以由表(3-1)和 表(3-2)得到. 上述方程和直角坐标系中的特征值问题 是相同的。 对于在? = 0 及? = ?0 处的各种边界条 件的组合，其对应的特征函数Φ( ν, ?)， 模 N(ν) 及特征值 ν 可以由表(2-2)查 到。 §3.3 乘积解 在第二章直角坐标系中的乘积解的原理和条件，同样适合于圆柱坐标系中的多维齐次热传导问题。 若初始条件可以表示为各个空间变 量的函数的乘积，则圆柱多维齐次 热传导问题的解可以化为一维圆柱 热传导问题的解的乘积。 注意：方程和边界条件都是齐次的。 §3.4 含热源问题(非齐次问题) 一、含热源项多维稳态问