拉普拉斯变换法在求解非稳态热传导问题中应用(II).pdfVIP

拉普拉斯变换法在求解非稳态热传导问题中的应用（II ） 1 2 1 1 2 2 许彬 ，John C. Chai ，谭俊杰 ，张敏 ，刘晶 ，王莎莎 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天学院，新加坡 (639798) E-mail ：norkistar@126.com 摘 要：在非结构化网格中，用基元中心有限容积法和全隐时间格式求解非稳态热传导问题， 并通过拉普拉斯变化法求解该问题,得到带有误差函数的精确解，将数值解和精确解进行比 较，并得到令人满意的结果。 关键词：非结构化网格，拉普拉斯变换，有限容积法，误差函数 拉普拉斯变换法是求解偏微分方程行之有效的方法之一。在前篇文章中，以非稳态热传 导偏微分方程为例，在结构化网格中，采用基元中心有限容积法和全隐时间格式，求解非稳 态热传导问题，并将数值计算得到的结果，同拉普拉斯变换求得的精确解进行比较，以此展 示拉普拉斯变换法求解方法的精彩之处。本文将在非结构化网格中，采用与前篇同样的算例 [1-6] 和方法，来进一步验证拉普拉斯变换法的可行性和实用性 。 1 算例分析和数值计算 一个半无限大物体（x ≥0 ），初始温度是零度。当时间t 0 时，在x 0 的边界表面 处，具有随时间变化的温度 [1] T f (t) 。该问题的数学描述为 ， 2 ∂ T ( x, t) 1 ∂T ( x, t) 2 0 =x ∞, t 0 (1a) x α t ∂ ∂ T (x ,t) f (t) x 0, t 0 (1b) T (x ,t) 0 x →∞, t 0 (1c) T (x ,t) 0 x ≥0, t 0 (1d) 对于方程式(1)时间t ，进行拉普拉斯变换，可以得到， 2 d T (x, s) s T ( x, s) 0 0 x (2a) − = ∞ dx2 α T (x ,s) f (s) x 0 (2b) T ( x, s) 0