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

拉普拉斯变换法在求解非稳态热传导问题中的应用 许彬 1，John C. Chai2，谭俊杰 1，张敏 1，刘晶 2，王莎莎 2 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 ∞ G(ω) = g(t)e ?iωt dt = f (t)e ?(σ +iω) dt (1.2) 2π ∫?∞ 2π ∫0 将σ + iω 记作 s ，并将G(ω) 改记为 f (s) / 2π ，则 ∞ ?st f (s) = f (t)e 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 f (