半球体容器漏水体积微分问题




问题:半球体的容器中盛满水,容器底部有一个小孔,水从小孔流出。给出水体积的变化量 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 kS2gh dt=−π(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=15kS2g 14π(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.




相关推荐
一只码代码的章鱼1 小时前
粒子群算法 笔记 数学建模
笔记·算法·数学建模·逻辑回归
艾思科蓝 AiScholar5 小时前
【连续多届EI稳定收录&出版级别高&高录用快检索】第五届机械设计与仿真国际学术会议(MDS 2025)
人工智能·数学建模·自然语言处理·系统架构·机器人·软件工程·拓扑学
胡萝卜不甜9 小时前
数学建模论文通用模板(细节方法二)
数学建模
金融OG1 天前
99.8 金融难点通俗解释:净资产收益率(ROE)
大数据·python·线性代数·机器学习·数学建模·金融·矩阵
spssau1 天前
2025美赛倒计时,数学建模五类模型40+常用算法及算法手册汇总
算法·数学建模·数据分析·spssau
C灿灿数模2 天前
2025美赛数学建模B题思路+模型+代码+论文
数学建模
金融OG2 天前
99.12 金融难点通俗解释:毛利率
python·算法·机器学习·数学建模·金融
小笼包数模3 天前
2025年美国大学生数学建模竞赛赛前准备计划
数学建模
Better Rose3 天前
2025美赛Latex模板可直接运行!O奖自用版
数学建模·latex·template method·美赛
Better Rose3 天前
【数学建模美赛速成系列】O奖论文绘图复现代码
数学建模·matlab