问题:半球体的容器中盛满水,容器底部有一个小孔,水从小孔流出。给出水体积的变化量 V 随水面高度 h 变化的微分关系式。
在微小的时间间隔
[
t
,
t
+
d
t
]
[t, t+\mathrm{d}t]
[t,t+dt] 内,水面高度由
h
h
h 降至
h
+
d
h
,
(
d
h
<
0
)
h+\mathrm{d}h, (\mathrm{d}h<0)
h+dh,(dh<0),水的体积变化量近似一个扁平的圆柱,所以可以利用圆柱的体积公式(
V
=
π
r
2
h
V=\pi r^2h
V=πr2h)写出
d
V
\mathrm{d}V
dV 和
d
h
\mathrm{d}h
dh 的关系:
d
V
=
−
π
r
2
d
h
\mathrm{d}V=-\pi r^2\mathrm{d}h
dV=−πr2dh
体积变化量
d
V
>
0
\mathrm{d}V>0
dV>0,高度变化量
d
h
<
0
\mathrm{d}h<0
dh<0,所以前面加负号。
假设半球的高度为 1 m 1m 1m,底部小孔横截面的面积为 1 c m 2 1cm^2 1cm2,则可以推出水面高度 h h h 随时间 t t t 变化的规律,并计算水流完所需的时间:
由流体力学水从孔口流出的流量(水的体积 V V V 对时间 t t t 的变化率) Q = d V d t = k S 2 g h Q=\dfrac{\mathrm{d}V}{\mathrm{d}t}=kS\sqrt{2gh} Q=dtdV=kS2gh. ( k k k 为流量系数, k = 0.62 k=0.62 k=0.62, S S S 为孔口横截面积, g g g 为重力加速度。)
t
t
t 时刻水面半径:
r
=
1
2
−
(
1
−
h
2
)
=
2
h
−
h
2
r=\sqrt{1^2-(1-h^2)}=\sqrt{2h-h^2}
r=12−(1−h2)=2h−h2,代入以上微分关系式:
d
V
=
−
π
(
2
h
−
h
2
)
d
h
\mathrm{d}V=-\pi(2h-h^2)\mathrm{d}h
dV=−π(2h−h2)dh
整合流量公式得:
k
S
2
g
h
d
t
=
−
π
(
2
h
−
h
2
)
d
h
kS\sqrt{2gh}\mathrm{d}t=-\pi(2h-h^2)\mathrm{d}h
kS2ghdt=−π(2h−h2)dh
上式即未知函数
h
=
h
(
t
)
h=h(t)
h=h(t) 满足的微分方程,分离变量得:
d
t
=
−
π
k
S
2
g
(
2
h
1
2
−
h
3
2
)
d
h
\mathrm{d}t=-\dfrac{\pi}{kS\sqrt{2g}}(2h^{\frac{1}{2}}-h^{\frac{3}{2}})\mathrm{d}h
dt=−kS2gπ(2h21−h23)dh
等式两端同时求积分得:
t
=
−
π
k
S
2
g
(
4
3
h
3
2
−
2
5
h
5
2
+
C
)
t=-\dfrac{\pi}{kS\sqrt{2g}}\Big(\dfrac{4}{3}h^\frac{3}{2}-\dfrac{2}{5}h^\frac{5}{2}+\mathrm{C}\Big)
t=−kS2gπ(34h23−52h25+C)
代入初值条件:
h
∣
t
=
0
=
1
h|_{t=0}=1
h∣t=0=1 得:
C
=
−
14
15
\mathrm{C}=-\dfrac{14}{15}
C=−1514.
于是有:
t
=
14
π
15
k
S
2
g
(
1
−
10
7
h
3
2
+
3
7
h
5
2
)
t=\dfrac{14\pi}{15kS\sqrt{2g}}\Big(1-\dfrac{10}{7}h^\frac{3}{2}+\dfrac{3}{7}h^\frac{5}{2}\Big)
t=15kS2g14π(1−710h23+73h25)
代入
k
=
0.62
,
S
=
1
0
−
4
m
2
,
g
=
9.8
m
/
s
2
k=0.62,\ \ S=10^{-4}m^2,\ \ \ g=9.8m/s^2
k=0.62, S=10−4m2, g=9.8m/s2 得:
t
=
1.068
×
1
0
4
(
1
−
10
7
h
3
2
+
3
7
h
5
2
)
t=1.068\times10^4\Big(1-\dfrac{10}{7}h^\frac{3}{2}+\dfrac{3}{7}h^\frac{5}{2}\Big)
t=1.068×104(1−710h23+73h25)
当水全部流出时,
h
=
0
h=0
h=0,代入上式便可得水流完所需得时间为:
t
=
1.068
×
1
0
4
s
=
2
h
58
m
i
n
t=1.068\times10^4\mathrm{s}=2\mathrm{h}\ 58\mathrm{min}
t=1.068×104s=2h 58min.