《热传导与热辐射》2、直角坐标系中的分离变量法.ppt

小 结 1 方程为齐次的，只有一个非齐次 边界条件，其解为 u . 2 方程为齐次的，有多个非齐次边 界条件：化为只有一个非齐次边界 问题的解的和 3 方程为非齐次的 含内热源 稳态 问题 一个齐次方程的通解 θ + 一个非齐次方程特解 p 4 非齐次非稳态问题 一个非齐次稳态问题的解 Ts + 一个齐次非稳态问题的解 Th 注意：上述问题中的非齐次部分都不 随时间变化。 完整数学描述 求解方程组 完 全 解 分离变量法求解热传导问题的步骤： 判断和分解 控制方程 边界条件 初始条件 （非）稳态、（非）齐次 分离变量法 特征函数 范 数 特征值 系 数 示 意 图 例2-6： 某一立方体， 0 ≤ x ≤ a, 0 ≤ y ≤ b ,0 ≤ z ≤ c 。立方体内有恒定的热源 q0 。在 x 0 , y 0 及 z 0 处绝热，在 x a , y b 及 z c 处以对流的方式向温度为零的周围环境散热。 初始温度为 F x, y, z 。求τ 0 时，立方体内 的温度分布 T x, y, z, τ 。 控制方程： 边界条件： 初始条件： 将方程分解为稳态、非齐次解 Ts 和非稳态、齐次解 Th 。 首先研究稳态、非齐次解： s s s 方程 2 的解 稳态齐次通解 θ + 非齐次问题特解 p 查表 2-4 ,得到 得到 将上式代入方程（2），得到 将 3 分解成仅有一个非齐次边界条件的两个方程，即： 将方程 4 分离变量 X1,Y1,Z1，查 表 2-2 ,求解得到 代入方程 4 的非齐次边界条件，得 用算子作用在方程的两边，利用特征函数的正交性，可以确定系数 Cmn 积分 7 式，得到 现在求解非稳态、齐次方程： 方程 10 的解为 查表 2-2 最后，我们得到方程 1 的解: The End 换一换？？？ ? 在实际工程中经常遇到两个或两个以上的非齐次边界问题，该类问题可以化为几个具有一个非齐次边界问题解的叠加。 分解为一组较为简单的解： 其中： 最终的解： 例2-4: 平板具有二个非齐次边界,见图 2-7 。 图 2-7 u, v 分别满足问题 2 , 3 。 二 稳态、齐次边界、非齐次方程（含有内热源）问题的求解 非齐次方程的解 T x,y,z 非齐次方程的特解 p x,y,z + 齐次方程的通解 θ x,y,z 即 式中 为非齐次方程 的特解，对各种 的特解由表 2－4 给出 p.72 。 三 稳态、非齐次边界、非齐次方程问题的求解 问题的解 T x,y,z 非齐次方程的特解 p x,y,z + 齐次方程非齐次边界条件的 通解 θ x,y,z 四 非稳态、非齐次边界、非齐次方程问题的求解 该问题的解可分为: 稳态、非齐次方程、非齐次边界的解 加上非稳态、齐次方程、齐次边界条件的解 非稳态、非齐次边界、 非齐次方程问题 将该问题分解成两个较为简单的问题： 一组按温度 定义 的稳态非齐次问题； 2.一个按温度 定义的齐次非稳态问题。 稳态非齐次问题 其中： 齐次非稳态问题 最终的解： 所以 即满足 2-27b 当 x L , 将上式及 2-27a 代入 2-27c ，经 整理，得到 特征值 βm 是上述超越方程的根。 根据模的定义，确定出 见第二章的注释2 p.86 对不同的边界条件，计算得到的 , 及 已经列于表 2-2 p.37 ，可直接查表。 半无限大物体的解: 半无限大物体的特点是只有一个 边界条件。 此时特征值 为从零到无限大连续域内的任意值，则解为 解法同上。 对于 x 0 的各种边界条件对应的特征函数 X β, x 和模 N β 列于表 2-3 p.43 。 无限大物体 –∞ x +∞ ,β值应从 – ∞ ~ + ∞积分，因为 –∞ x 0 和 0 x + ∞是对称的，所以 β 值从 0 ~ +∞积分即可。 （二）多维齐次 边条 问题 可以分解为几个一维齐次问题求解。 例 2-3: 考虑一矩形区域，初始和边界条件如下图中所示，求矩形区域的温度分布。见图2-3。 图2-3 矩形区域的温度分布 解: 控制方程： 边界条件： 初始条件： 上述问题可以分解为两个一维问题： 及 完全解的形式为: 用算子 及算子 依次作用在上式的两边，利用特 征函数的正交性求出 模 及特征值 可以从表 2-2 中查到。 半无限长的带子，见图 2-4 , 图 2－4 对 y 方向为一维有限大物体，对 x 方向是半无限大物体，可利用 表 2-2 及 2-3 。