化工问题的建模与数学分析方法 第四章习题及答案.doc

第四章习题 1.（√）判别以下方程的类型，并指出变系数中自变量取值范围 （1） （2） （3） （4） 解：（1）a＝1，b＝2，c＝3，b2－ac0，是双曲型方程 （2）b＝0.5，b2－ac0，是双曲型方程 （3）b2－ac＝x2y2，当x＝0或y＝0时，是抛物线型方程，否则是双曲型方程 （4）b2－ac＝xy，当x＝0或y＝0时，是抛物线型方程，当x和y同号时是双曲型方程，否则是椭圆型方程 2. （√）证明： （1）圆形区域上Laplace方程在圆对称情况下的通解为 式中r为径向极坐标，A、B为任意常数 （2）球形区域上Laplace方程在球对称情况下的通解为 式中r为径向球坐标，A、B为任意常数 证明：（1）在极坐标下，圆型区域内，laplace方程的表达式为 在圆对称情况下 ， 原方程可化为 ， 解得 （2）同理可证 3. （√） 用分离变量法求解以下一维热传导方程的定解问题 (1) 解：设 由（4） （5）代入（2） 其中An由下式确定： (2) 解：首先将方程化齐，为此，令 代入方程，得 解得 于是 用分离变量法解问题（12） 得特征值 特征函数 及 因此，问题的解为 4. 解下列矩形域的拉普拉斯方程定解问题 解： 首先假设问题的解具有以下变量分离的形式 （1） 代入原方程，得到关于X(x)的特征值问题和Y(y)的方程 （2） （3） 类似与书P199中对于方程（4.2.5）、（4.2.6）的讨论，这里的参数λ只能取正值，否则只能得到零解。因此方程（2）的通解是 （4） 由边界条件X(0)=0知c2=0，再由知，若要有非零解, c1≠0，必须＝0，由此λ应取以下值 （n=1,2,3,…） （5） 由此得到 （n=1,2,3,…） （6） 由方程（3）可得 （7） 将λ值代入方程（7），得 （n=1,2,3,…）（8） 从而由方程（1）得到， （9） 式（9）中An=A1nCn，Bn= B1nCn，均为任意常数。 根据叠加原理构造以下级数形式的解 （10） 令上式满足y的边界条件，得到确定系数An， Bn的方程 将以上两式的右端展开为的傅立叶级数，然后逐项比较系数，得到 由此解得 代入方程（10），于是问题的解最终为 （n=1,2,3,…） 5．对于第一章习题2所述的池塘结冰问题，如果空气温度TW每天呈现周期性变化，其规律用以下方程描述 求冰层中的温度分布及厚度l的时间变化趋势，若给定冰的导温系数，TW＝－10℃，(T＝5℃，( ＝86400t的Fourier级数以反映天气的逐日变化和中长期变化的影响。 解： a）可用一维热传导方程来描述该过程，假定厚度l不变 假定定解为下列形式，T=F(t)G(x)。由于温度变化是周期的，则F(t)中解的指数为虚数，原一维热传导方程可化为 即 解得 由此将产生4个特解 当时，T有界，故将T4和T3舍去，将T1和T2相加得到另一个特解 用三角函数进行表示， 依照边界条件确定常数A，B和λ 比较得 在x＝0处，其平均温度是零瞬间时初始瞬变的结果。这一变化的阶梯解是： 将此式子附加至解，得: 上式并不满足初始条件，所以不能用于很小的时间范围。但是对于很大的t值，趋于零，结果形成的表示式可以很好地表示在初始瞬变已消失的情况下的周期温度分布： b）厚度l的时间变化趋势 由一维热传导公式可得： 其中为潜热 两端分别积分可得 将具体数值代入上式，用matlab解得t＝6882779.6s＝79.6天。 随着时间的增加，负积温逐渐累积，冰层不断加厚，但是速度越来越缓慢。 6．（√）求图示的半环形区域内的稳态温度分布() 边界条件为：当r＝c时，u＝u0，其余边界保持0度。 解：极坐标系下的热传导方程为 （提示：本题作法与第四章第2节情况3介绍的圆上Laplace方程的第一