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




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




相关推荐
Martina_032116 小时前
AI生成3D场景导入 Unity/Unreal 后,NPC 总串层?用6步检查导航高度与跨层连接
人工智能·游戏·数学建模·3d·unity·自然语言处理·aigc
考研保研资料分享18 小时前
电子信息保研经验合辑:复旦直博与北邮集成电路,夏令营、预推免和专业面试
数学建模·面试·职场和发展·信息与通信
蒲公英eric3 天前
人口增长问题:从 Malthus 到 Logistic 的建模之旅
人工智能·python·数学建模·人口增长问题
Martina_03213 天前
同一片森林换个镜头就忽冷忽暖?用6步统一曝光、LUT与相机后处理
人工智能·游戏·数学建模·3d·自然语言处理·aigc·游戏策划
蒲公英eric4 天前
Python工具链与LaTeX排版:把兵器配齐
人工智能·python·算法·数学建模·数学建模工具
当时皓月4 天前
Spur Generation and Modulation(杂散的产生与调制)
数学建模·芯片
Zootopia6264 天前
多架 eVTOL集群排班调度方案思考
人工智能·算法·数学建模·matlab·动态规划·无人机·evtol
Doubbbbbbble云5 天前
算法优化中的数学建模与理论界限分析4
算法·数学建模
Zootopia6267 天前
美国载人飞船Crew-13明晚发射!
人工智能·算法·机器学习·数学建模·无人机·创业创新·信息与通信
Studying 开龙wu8 天前
2026年“华为杯”第二十三届中国研究生数学建模竞赛“D题:山区洪涝灾害下无人机运输与通信协同优化解答总结
数学建模·无人机