拉普拉斯变换法在求解非稳态热传导问题中的应用-论文完整版.pdf

拉普拉斯变换法在求解非稳态热传导问题中的应用 1 2 1 1 2 2 许彬 ，John C. Chai ，谭俊杰 ，张敏 ，刘晶 ，王莎莎 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天学院，新加坡 (639798) E-mail ：norkistar@126.com 摘 要：在结构化网格中，用基元中心有限容积法和全隐时间格式求解非稳态热传导问题， 并通过拉普拉斯变化法求解该问题,得到带有误差函数的精确解，将数值解和精确解进行比 较，并得到令人满意的结果。 关键词：拉普拉斯变换，有限容积法，误差函数 拉普拉斯变换法是求解偏微分方程行之有效的方法之一。在应用这个方法时，先利用拉 普拉斯变换把对时间变量的偏导数消掉，由此得到的方程求解温度的变换量，再将温度的变 换量进行反变换，以找到温度场的分布解。这种方法在原理上是简单易行的，但是，除非变 换在拉普拉斯变换表中可以查到，一般来说，反变换很困难[1-6] 。 本文以非稳态热传导偏微分方程为例，在结构化网格中，采用基元中心有限容积法和全 隐时间格式，求解非稳态热传导问题，并将数值计算得到的结果同拉普拉斯变换求得的精确 解进行比较，以此展示两种求解方法的精彩之处。 1 基本概念和方程 拉普拉斯变化常用于初始值问题，即已知某个物理量在初始时刻t 0 的值f (0) ，而求 解它在初始时刻之后的变化情况f (t) ，至于它在初始时刻之前的值，我们并不感兴趣。不 妨置f (t) 0 (t 0) ，为了获得宽松的变换条件，把f (t) 加工为g (t) [2]， −σt g (t) e f (t) (1.1) 这里e−σt 是收敛因子，就是说，正的实数σ 的值选得如此之大，以保证g (t) 在区间(−∞,∞) 上绝对可积。于是，可以对g (t) 施行傅立叶变换 1 ∞ 1 ∞ iωt (σ iω) ( ) ( ) − − + (1.2) G ω ∫ g t e dt ∫ f (t)e dt 2π −∞ 2π 0 将σ +iω记作s ，并将G(ω) 改记为f (s) / 2π ，则 ∞ f (s) ∫ f (t)e−st dt (1.3) 0 ∞ 其中，积分∫ f (t)e−st dt 称为拉普拉斯积分，f (s) 称为f (t) 的拉普拉斯变换函数，式（1.3） 0 代表着从f (t) 到f (s) 的一种积分变换，称为拉普拉斯变换（简称拉氏变换），e−st 称为拉 普拉斯变换的核。 拉普拉斯变换，即式(1.3)中的积分存在的条件是（1）在0 ≤t ∞的任一区间上，除了 有限个第一类间断点外，函数f (t) 及其导数是处处连续的。(2)存在常数M 0 和σ ≥0 ， 使得对于任何t 值（0 ≤t ∞），有 σt