《热传导与热辐射》3、圆柱坐标系中的分离变量法.pptVIP

二、含热源项非稳态问题 控制方程： 边界条件： 初始条件： 该问题可以化为： 稳态非齐次问题 课堂练习（书中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 三、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 含热源问题 非齐次问题 一、含热源项多维稳态问题 控制方程: 利用第二章介绍过的方法，方程的解 为齐次方程 的通解 和 3-30 的特解 之和: 非齐次方程的特解 p r, ?, z 可以由 表3-4 查到。 利用Bessel函数的积分求模 见梁昆淼：数学物理方法 ， 在 r b 处对三类边界条件求出 特征函数: ， 模: 及特征值: 讨论在 r b 处第三类边界条件的情况: 3-10a 的解为 3-16 要使上述解满足 r b 处的第三类边界条件，即 3-10b ， 特征值 应该是下列方程的正根: 求模 由模的定义 3-13b ，利用Bessel函数的积分 3-15 及解 3-16 ，得到 对 r