圆柱坐标系中非齐次热传导问题简单分解方法的精确解与数值解.pdfVIP

圆柱坐标系中非齐次热传导问题简单分解方法的精 确解与数值解 1 1 2 1,2 张钧波 ，张敏 ，John C. Chai ，许彬 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天学院，新加坡 (639798) E-mail ：zhangm@ 摘 要：在圆柱坐标系中，对非齐次热传导问题进行简单分解，并运用分离变量方法求得其 精确解。与此同时，采用有限容积方法，在结构化网格中进行数值计算。通过比较计算结果， 两者之间得到令人满意的一致性，这证明了简单分解方法的正确性和实用性。 关键词：非齐次热传导，数值计算，精确解 在求解由于偏微分方程和边界条件的非齐次造成的非齐次热传导问题时，一般采用简单 分解方法[1-4] 。在前面的文章中，我们已经在直角坐标系中，对非齐次热传导问题进行了简 单分解方法的分析和计算[5] 。现在，我们再来讨论一下在圆柱坐标系中，非齐次热传导问题 的简单分解求解。同时采用有限容积方法，在结构化网格中进行数值计算。 1. 理论分析 [5] 非齐次热传导问题简单分解方法的理论，在参考文献 中已做过具体介绍，在此不再累 述。这里我们直接对圆柱坐标系中非齐次热传导问题进行分析。 实心圆柱体（0 ≤r ≤b ），初始温度分布为F (r ) 。时间t 0 时，在该圆柱内以恒定的 速率g0 产生热量；r b 处的边界温度保持摄氏零度不变。本问题的数学描述为， 2 ∂ T 1 ∂T 1 1 ∂T 2 + + g0 0 =≤r ≤b t 0 (1.1a) ∂r r ∂r k α ∂t T 0 r b, t 0 (1.1b) T F r 0 =r b, t 0 (1.1c) ( ) 把该问题分解成一个T r 的稳态问题为， s ( ) 2 d Ts 1 dTs g0 2 + + 0 0 =≤r b (1.2a) dr r dr k T 0 r b (1.2b) s 与一个T x, t 的齐次问题为, h ( ) 2 ∂ T 1 ∂T 1 ∂T h + h h 0 =≤r ≤b t 0 (1.3a) ∂ 2 ∂ α ∂ r