热传导问题
物理问题:空间某个物体或静止流体内温度分布不均匀,引起热量流动及温度的变化。
理想化假设:
- 物体由同一介质构成,且介质均匀分布、各向同性
- 介质的密度、比热和热传导系数均为常数。
物理定律:
- 能量守恒定律
- 傅里叶热传导定律
数学建模(1):
-
在空间取定直角坐标系
-
取各点在t时刻的温度 u = u ( t , x , y , z ) u=u(t,x,y,z) u=u(t,x,y,z)为热运动的表征量
-
**微元分析:**在介质内任取微元 d V = [ ( x , y , z ) , ( x + d x , y + d y , z + d z ) ] dV=[(x,y,z),(x+dx,y+dy,z+dz)] dV=[(x,y,z),(x+dx,y+dy,z+dz)],考察微元 d V dV dV在时间间隔 [ t , t + d t ] [t,t+dt] [t,t+dt]内的温度变化
-
微元满足能量守恒定律, [ t , t + d t ] [t, t+dt] [t,t+dt]
外界流入热量 + 内部热源产热 = 温度升高所需热量
Q 流 入 + Q 热 源 = Q 温 度 升 高 Q_{流入}+Q_{热源}=Q_{温度升高} Q流入+Q热源=Q温度升高
傅里叶热传导定律:热量从高温处向低温处流动,沿某方向流动热量的多少与温度在该方向的减少率成比例。
q → = − k ∇ u = { q x = − k ∂ u ∂ x q y = − k ∂ u ∂ y q z = − k ∂ u ∂ z \overrightarrow q=-k\nabla u = \begin{cases} q_x = -k\frac{\partial u}{\partial x} \\ q_y = -k\frac{\partial u}{\partial y} \\ q_z = -k\frac{\partial u}{\partial z} \end{cases} q=−k∇u=⎩⎪⎨⎪⎧qx=−k∂x∂uqy=−k∂y∂uqz=−k∂z∂u
其中 q → \overrightarrow q q是热流密度矢量,表示单位时间沿单位面积的法向流出的热量。∴ \therefore ∴
Q 左 右 = q ∣ x ⋅ d t ⋅ d y d z − q ∣ x + d x ⋅ d t ⋅ d y d z = ( q ∣ x − q ∣ x + d x ) ⋅ d t d y d z = − ∂ q ∂ x d x ⋅ d t d y d z = − ∂ ∂ x ( − k ∂ u ∂ x ) ⋅ d t d V = k ∂ 2 u ∂ x 2 d t d V Q_{左右}=q|_x·dt·dydz-q|_{x+dx}·dt·dydz \\ =(q|_x-q|_{x+dx})·dtdydz \\ =-\frac{\partial q}{\partial x}dx·dtdydz =-\frac{\partial}{\partial x}(-k\frac{\partial u}{\partial x})·dtdV \\ =k\frac{\partial^2u}{\partial x^2}dtdV Q左右=q∣x⋅dt⋅dydz−q∣x+dx⋅dt⋅dydz=(q∣x−q∣x+dx)⋅dtdydz=−∂x∂qdx⋅dtdydz=−∂x∂(−k∂x∂u)⋅dtdV=k∂x2∂2udtdVQ 前 后 = q ∣ y ⋅ d t ⋅ d x d z − q ∣ y + d y ⋅ d t ⋅ d x d z = ( q ∣ y − q ∣ y + d y ) ⋅ d t d x d z = − ∂ q ∂ y d y ⋅ d t d x d z = − ∂ ∂ y ( − k ∂ u ∂ y ) ⋅ d t d V = k ∂ 2 u ∂ y 2 ⋅ d t d V Q_{前后}=q|_y·dt·dxdz-q|_{y+dy}·dt·dxdz \\ =(q|_y-q|_{y+dy})·dtdxdz \\ =-\frac{\partial q}{\partial y}dy·dtdxdz = -\frac{\partial}{\partial y}(-k\frac{\partial u}{\partial y})·dtdV \\ =k\frac{\partial^2 u}{\partial y^2}·dtdV Q前后=q∣y⋅dt⋅dxdz−q∣y+dy⋅dt⋅dxdz=(q∣y−q∣y+dy)⋅dtdxdz=−∂y∂qdy⋅dtdxdz=−∂y∂(−k∂y∂u)⋅dtdV=k∂y2∂2u⋅dtdV
Q 上 下 = k ∂ 2 u ∂ z 2 ⋅ d t d V Q_{上下}=k\frac{\partial^2 u }{\partial z^2}·dtdV Q上下=k∂z2∂2u⋅dtdV
Q 流 入 = k ( ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 ) ⋅ d t d V = k Δ u ⋅ d t d V Q_{流入}=k(\frac{\partial^2u}{\partial x^2}+\frac{\partial^2u}{\partial y^2}+\frac{\partial^2u}{\partial z^2})·dtdV=k\Delta u·dtdV Q流入=k(∂x2∂2u+∂y2∂2u+∂z2∂2u)⋅dtdV=kΔu⋅dtdV
∵ \because ∵
Q 热 源 = g ( x , y , z , t ) d t d V Q_{热源}=g(x,y,z,t)dtdV Q热源=g(x,y,z,t)dtdV
g ( x , y , z , t ) g(x,y,z,t) g(x,y,z,t)表示单位体积内部热源的产热率(单位时间单位面积产热量)
Q 升 温 = c ⋅ ρ d V ⋅ [ u ( t + d t , x , y , z ) − u ( t , x , y , z ) ] = c ⋅ ρ d V ⋅ ∂ u ∂ t d t = c ρ ∂ u ∂ t d t d V Q_{升温}=c·\rho dV·[u(t+dt,x,y,z)-u(t,x,y,z)] \\=c·\rho dV·\frac{\partial u}{\partial t}dt \\=c\rho \frac{\partial u}{\partial t}dtdV Q升温=c⋅ρdV⋅[u(t+dt,x,y,z)−u(t,x,y,z)]=c⋅ρdV⋅∂t∂udt=cρ∂t∂udtdV
由能量守恒定律得
k Δ u ⋅ d t d V + g ( x , y , z , t ) d t d V = c ρ ∂ u ∂ t d t d V k c ρ Δ u + g ( x , y , z , t ) c ρ = ∂ u ∂ t k\Delta u·dtdV + g(x,y,z,t)dtdV = c\rho \frac{\partial u}{\partial t}dtdV \\\frac{k}{c\rho}\Delta u+\frac{g(x,y,z,t)}{c\rho}=\frac{\partial u}{\partial t} kΔu⋅dtdV+g(x,y,z,t)dtdV=cρ∂t∂udtdVcρkΔu+cρg(x,y,z,t)=∂t∂u热传导方程(扩散方程):
∂ u ∂ t = a 2 Δ u + f ( t , x → ) , a = κ c ρ , f ( t , x → ) = g ( t , x → ) c ρ \frac{\partial u}{\partial t}=a^2\Delta u+f(t,\overrightarrow x),\space a=\sqrt{\frac{\kappa}{c\rho}},f(t,\overrightarrow x)=\frac{g(t,\overrightarrow x)}{c\rho} ∂t∂u=a2Δu+f(t,x), a=cρκ,f(t,x)=cρg(t,x)
其中, κ \kappa κ为热扩散系数。量纲分析:
[ k ] [ c ] ⋅ [ ρ ] = J / ( s ⋅ m ⋅ K ) J / ( k g ⋅ K ) ⋅ k g / m 3 = m / s [ u ] [ t ] = [ a 2 ] ⋅ [ u ] [ x 2 ] ⟹ [ a 2 ] = m 2 / s \sqrt{\frac{[k]}{[c]·[\rho]}} = \sqrt{\frac{J/(s·m·K)}{J/(kg·K)·kg/m^3}}=m/\sqrt{s} \\ \frac{[u]}{[t]}=[a^2]·\frac{[u]}{[x^2]} \implies [a^2]=m^2/s [c]⋅[ρ][k]=J/(kg⋅K)⋅kg/m3J/(s⋅m⋅K)=m/s[t][u]=[a2]⋅[x2][u]⟹[a2]=m2/s
根据量纲可知,扩散传播距离与时间之间的关系:
x 2 ∝ a 2 t x^2 \propto a^2t x2∝a2t
数学建模(2):
-
在空间取定直角坐标系
-
取各点在t时刻的温度 u = u ( t , x , y , z ) u=u(t,x,y,z) u=u(t,x,y,z)为热运动的表征量。
-
在介质内任取微元 d V = [ x , x + d x ] × [ y , y + d y ] × [ z , z + d z ] dV=[x,x+dx]\times [y,y+dy]\times [z,z+dz] dV=[x,x+dx]×[y,y+dy]×[z,z+dz],考察微元 d V dV dV在时间间隔 [ t , t + d t ] [t,t+dt] [t,t+dt]内的温度变化。
-
根据能量守恒定律,物体温度升高所需热量等于外部流入热量和内部热源产生热量之和。
热量的流动遵循傅里叶热传导定律:热量从温度高处向低处,沿某方向流动热量的多少与温度在该方向的减少率成比例,其数学表示式为
Q n = − k ( x , y , z ; n ) ∂ u ∂ n n \bold Q_n=-k(x,y,z;n)\frac{\partial u}{\partial n}\bold n Qn=−k(x,y,z;n)∂n∂un
其中, Q n \bold Q_n Qn为 n \bold n n方向的热流密度矢量,即单位时间沿 n \bold n n方向通过单位面积的热量; k ( x , y , z ; n ) k(x,y,z;n) k(x,y,z;n)为介质的热传导系数,在介质均匀,各项同性假设下是常数,记为k。在 [ t , t + d t ] [t,t+dt] [t,t+dt]时间间隔内通过微元的左右面传入的热量为
− k ∂ u ∂ x ∣ ( t , x , y , z ) d t d y d z + k ∂ u ∂ x ∣ ( t , x + d x , y , z ) ≈ k ∂ 2 u ∂ x 2 ∣ ( t , x , y , z ) d t d x d y d z -k\frac{\partial u}{\partial x}|_{(t,x,y,z)}dtdydz+k\frac{\partial u}{\partial x}|_{(t,x+dx,y,z)} \approx k\frac{\partial^2 u}{\partial x^2}|_{(t,x,y,z)}dtdxdydz −k∂x∂u∣(t,x,y,z)dtdydz+k∂x∂u∣(t,x+dx,y,z)≈k∂x2∂2u∣(t,x,y,z)dtdxdydz
同样可以求出通过前后和上下面流入的热量分别为
k ∂ 2 u ∂ y 2 d t d x d y d z 和 k ∂ 2 u ∂ z 2 d t d x d y d z k\frac{\partial^2u}{\partial y^2}dtdxdydz \quad 和 \quad k\frac{\partial^2u}{\partial z^2}dtdxdydz k∂y2∂2udtdxdydz和k∂z2∂2udtdxdydz
如果介质内部有热源,其热源密度,即单位时间单位体积热源流出的热量为 g ( t , x , y , z ) g(t,x,y,z) g(t,x,y,z),则在 [ t , t + d t ] [t,t+dt] [t,t+dt]时间间隔内,微元内部热源流出热量为
g ( t , x , y , z ) d t d x d y d z g(t,x,y,z)dtdxdydz g(t,x,y,z)dtdxdydz
而微元温度升高所需的热量为
c ρ [ u ( t + d t , x , y , z ) − u ( t , x , y , z ) ] d x d y d z ≈ c ρ ∂ u ∂ t d t d x d y d z c\rho[u(t+dt,x,y,z)-u(t,x,y,z)]dxdydz \approx c\rho \frac{\partial u}{\partial t}dtdxdydz cρ[u(t+dt,x,y,z)−u(t,x,y,z)]dxdydz≈cρ∂t∂udtdxdydz
这些等式中都忽略了高阶无穷小量讲这些量代入能量守恒定律,便得方程
c ρ ∂ u ∂ t = k ( ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 ) + g ( t , x , y , z ) c\rho \frac{\partial u}{\partial t}=k(\frac{\partial^2u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2}+\frac{\partial^2 u}{\partial z^2})+g(t,x,y,z) cρ∂t∂u=k(∂x2∂2u+∂y2∂2u+∂z2∂2u)+g(t,x,y,z)
即热传导方程
∂ u ∂ t = a 2 Δ u + f ( t , x , y , z ) \frac{\partial u}{\partial t}=a^2\Delta u +f(t,x,y,z) ∂t∂u=a2Δu+f(t,x,y,z)
其中, Δ = ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 \Delta=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2} Δ=∂x2∂2+∂y2∂2+∂z2∂2为三维拉普拉斯算子, a = k c ρ , f ( t , x , y , z ) = g ( t , x , y , z ) c ρ a=\sqrt{\frac{k}{c\rho}}, f(t,x,y,z)=\frac{g(t,x,y,z)}{c\rho} a=cρk,f(t,x,y,z)=cρg(t,x,y,z)如果考虑侧面绝热杆的温度,或柱上与高度无关的温度变化,同样可导出热传导方程,只是拉普拉斯算子相应地取为一维或二维。
总结:热传导方程的建立基于能量守恒和热传导两条基本物理定律。像气体扩散、杂志在固体或液体中扩散这些物理过程,其机理与热传导相似,都是由浓度的不均匀引起不同物质分子的位置变换,变换过程中每种物质的总量保持不变。选取适当的未知函数,导出的方程与热传导方程有相同形式,因此也称热传导方程为扩散方程。
波动方程和热传导方程分别描述了双向传播和单向传播两种完全不同的物理过程。它们都与时间t有关,称为发展方程。如果考虑热传导方程的稳恒状态,即
∂
u
∂
t
≡
0
\frac{\partial u}{\partial t}\equiv 0
∂t∂u≡0,它就成为泊松(Poisson)方程
Δ
u
=
−
1
a
2
f
(
x
,
y
,
z
)
\Delta u=-\frac{1}{a^2}f(x,y,z)
Δu=−a21f(x,y,z)
当
f
(
x
,
y
,
z
)
≡
0
f(x,y,z)\equiv 0
f(x,y,z)≡0时,就是Laplace方程(也称调和方程)
Δ
u
=
0
\Delta u=0
Δu=0