数学物理方程练习题

σf 4dSdt.
根据热量平衡有 故所求边界条件为
−k
∂u ∂n
dSdt
=
σu4dSdt
−
σf
4dSdt.
−k
∂u ∂n
=
σ(u4
− f 4).
齐海涛 (SDU)
数学物理方程
2012-10-3 12 / 49
1. 热传导方程及其定解问题的导出 2. 初边值问题的分离变量法 3. 柯西问题 4. 极值原理、定解问题解的唯一性和稳定性 5. 解的渐近性态
dQ = −βQ, dt Q(0) = Q0,
⇒ Q(t) = Q0e−βt.
易知 t1 到 t2 时刻, 砼内任一区域 Ω 中的热量的增加等于从 Ω 外部流入 Ω 的热量及砼中的水化热之和, 即
齐海涛 (SDU)
数学物理方程
2012-10-3 7 / 49
热传导方程及其定解问题的导出
∫ t2 cρ ∂u dtdxdydz =
.
热传导方程
.
Heat Equations
齐海涛
山东大学（威海）数学与统计学院
htqisdu@
齐海涛 (SDU)
数学物理方程
2012-10-3 1 / 49
目录
1. 热传导方程及其定解问题的导出 2. 初边值问题的分离变量法 3. 柯西问题 4. 极值原理、定解问题解的唯一性和稳定性 5. 解的渐近性态
热传导方程及其定解问题的导出
.E.xample 1.2
.试直接推导扩散过程所满足的微分方程.
解: 设 N(x, y, z, t) 表示在时刻 t, (x, y, z) 点处扩散物质的浓度, D(x, y, z) 为 扩散系数, 在无穷小时间段 dt 内, 通过无穷小曲面块 dS 的质量为
dm
=
−D(x,
解: ∑∞
u(x, t) = Cke−k2π2t sin kπx.
k=1
Ck = 2 ∫
1
∫
2
ξ sin kπξdξ +
1(1 − ξ) sin kπξdξ
=
0
4 kπ k2π2 sin 2
=
1 2
0,
4(−1)n (2n + 1)2π2
,
k = 2n, k = 2n + 1,
n = 0, 1, 2, . . .
齐海涛 (SDU)
数学物理方程
2012-10-3 9 / 49
热传导方程及其定解问题的导出
解: 与第1题类似, 取导线轴为 x 轴, 在时刻 t1 到 t2 介于 [x1, x2] 的导线段 的热量增加为: 从导线的其它部分流入的热量, 从侧面流入的热量以及电流通
过 [x1, x2] 这段产生的热量之和, 即
dQ1 = −k1(u − u1)πl∆x;
单位时间从 x 处, x + ∆x 处流入的热量为
dQ2
=
−k(x)
∂u ∂x
(x,
t)
·
πl2 4
,
dQ3
=
k(x
+
∂u ∆x) ∂x (x
+
∆x,
t)
·
πl2 4
,
齐海涛 (SDU)
数学物理方程
2012-10-3 3 / 49
热传导方程及其定解问题的导出
另外, 从时刻 t1 到 t2, Ω 中该物质的增加为
[N(x, y, z, t2) − N(x, y, z, t1)]dxdydz =
Ω
∫ t2 ∂N dtdxdydz. Ω t1 ∂t
根据质量守恒, 并注意到 Ω, t1, t2 的任意性, 得所求方程为
∂N ∂t
=
∂ ∂x
( ∂N ) D ∂x
y,
z)
∂N ∂n
dSdt.
因此从时刻 t1 到 t2 流入区域 Ω (Γ 为 Ω 的表面) 的质量为
∫ t2
D(x, y, z) ∂N dSdt = ∫ t2
t1
Γ
∂n
t1
div(DgradN)dxdydzdt.
Ω
齐海涛 (SDU)
数学物理方程
2012-10-3 5 / 49
热传导方程及其定解问题的导出
用分离变量法求解热传导方程的初边值问题:
.
ut = uxx u(x, 0) =
({t
> 0, x, 1−
0< x,
x
0
1 2
< 1), <x≤ <x<
1 2
,
1,
u(0, t) = u(1, t) = 0 (t > 0).
齐海涛 (SDU)
数学物理方程
2012-10-3 15 / 49
初边值问题的分离变量法
齐海涛 (SDU)
数学物理方程
2012-10-3 2 / 49
1. 热传导方程及其定解问题的导出 2. 初边值问题的分离变量法 3. 柯西问题 4. 极值原理、定解问题解的唯一性和稳定性 5. 解的渐近性态
齐海涛 (SDU)
数学物理方程
2012-10-3 3 / 49
热传导方程及其定解问题的导出
故单位时间流入 (x, x + ∆x) 的热量为
( ∂
) ∂u
πl2
dQ = dQ1 + dQ2 + dQ3 = ∂x
k(x) ∂x
·
x∗
4 ∆x − k1(u − u1)πl∆x.
综上, 从时刻 t1 到 t2 流入位于 [x1, x2] 杆段的热量为
∫ t2
t1
∫ x2
x1
[ ∂ ∂x
(
)
∂u
k(x) ∂x
πl2 4
−
k1(u
−
] u1)πl dxdt.
而在这段时间内 [x1, x2] 杆段内各点温度从 u(x, t1) 变到 u(x, t2), 其吸收热量
为
∫ x2
x1
cρ(u(x, t2)
−
u(x,
t1
))
πl2 4
dx
=
∫ t2
t1
∫ x2
x1
πl2 4
cρ
∂u ∂t
dxdt.
根据热量守恒, 并注意到 x1, x2, t1, t2 的任意性, 得所求方程为
齐海涛 (SDU)
数学物理方程
2012-10-3 11 / 49
热传导方程及其定解问题的导出
.E.xample 1.5
设物体表面的绝对温度为 u, 此时它向外界辐射出去的热量依斯特藩-玻耳兹曼 (Stefan-Boltzmann) 定律正比于 u4, 即
dQ = σu4dSdt.
假设物体和周围介质之间只有热辐射而没有热传导, 又假设物体周围介质的绝 .对温度为已知函数 f(x, y, z, t), 求此时该物体热传导问题的边界条件.
初始时刻所储的热量, 则 dQ .c, 密度为 ρ, 热传导系数为dt
= −βQ, 其中 β 为正常数. 又假设砼的比热为 k, 求它在浇筑后温度 u 满足的方程.
齐海涛 (SDU)
数学物理方程
2012-10-3 7 / 49
热传导方程及其定解问题的导出
.E.xample 1.3
砼(混泥土)内部储藏着热量, 称为水化热, 在它浇筑后逐渐放出, 放热速度和
⇒
u(x,
t)
=
∑∞
4(−1)n (2n + 1)2π2
+
∂ ∂y
( ∂N ) D ∂y
+
∂ ∂z
( D
∂N ∂z
)
.
齐海涛 (SDU)
数学物理方程
2012-10-3 6 / 49
热传导方程及其定解问题的导出
.E.xample 1.3
砼(混泥土)内部储藏着热量, 称为水化热, 在它浇筑后逐渐放出, 放热速度和
它所储藏的水化热成正比. 以 Q(t) 表示它在单位体积中所储的热量, Q0 为
解: 设 u(x, t) = X(x)T(t), 则
T ′ + λa2T = 0,
X ′′ + λX = 0, X(0) = X ′(π) = 0.
⇒ λk
=
( k
+
1 )2 2
,
k = 0, 1, 2, . . .
齐海涛 (SDU)
数学物理方程
2012-10-3 13 / 49
初边值问题的分离变量法
齐海涛 (SDU)
数学物理方程
2012-10-3 13 / 49
初边值问题的分离变量法
.E.xample 2.1
用分离变量法求下列定解问题的解:
.
ut = a2uxx (t > 0, 0 < x < π), u(0, t) = ux(π, t) = 0 (t > 0), u(x, 0) = f(x) (0 < x < π).
t1
∂t
(Q(t1) − Q(t2))dxdydz+
Ω∫ t2
[( )
Ω( )
( )]
∂ ∂u ∂ ∂u ∂ ∂u
k + k + k dxdydzdt
t1
∂x ∂x ∂y ∂y ∂z ∂z
Ω ∫ t2 dQ
=−
dtdxdydz+
t1 dt
∫ t2 Ω
[∂
( ∂u ) k+
∂
( ∂u ) k+
∂
( k
解: 考察边界上的面积微元 dS. 在 dt 时间内, 经边界微元流出的热量为 (k