水箱排水模拟器 返回
流体力学模拟器

水箱排水模拟器 — 托里拆利定理

用托里拆利定理实时可视化通过底部孔口排水的水箱水位下降。改变横截面积、孔口面积和流量系数,学习非定常流出的特性。

参数设置
水箱横截面积 A
孔口面积 a
cm²
初始水位 h₀
m
流量系数 Cd
经过时间 t
%

「经过时间 t」是相对于排空时间的百分比(%)。

计算结果
当前水位 h
流出速度 v
当前流量 Q
排空时间
水箱与水位

蓝色=当前水位 / 底部黄色箭头=从孔口的流出

水位随时间变化 h(t)

横轴=经过时间 t / 纵轴=水位 h(黄色点=当前时刻)

理论与主要公式

处理横截面积为 A 的水箱通过面积为 a 的孔口排水的非定常流出。

托里拆利定理给出的流出速度与流量(Cd 为流量系数):

$$v = C_d \sqrt{2gh}, \qquad Q = a\,v = C_d\,a\,\sqrt{2gh}$$

水位的控制方程(质量守恒):

$$A\,\frac{dh}{dt} = -Q = -C_d\,a\,\sqrt{2gh}$$

求解后,水位以时间的二次函数减小:

$$h(t) = \left(\sqrt{h_0} - \frac{C_d\,a}{A}\sqrt{\frac{g}{2}}\;t\right)^{2}$$

排空时间通过在上式中令 h(t) = 0 求得:

$$t_{empty} = \frac{A}{C_d\,a}\sqrt{\frac{2h_0}{g}}$$

什么是水箱排水模拟器

🙋
拔掉浴缸的塞子后,水流起初很急、然后逐渐变慢,对吧?这是为什么呢?
🎓
注意得很好。简单来说,是因为流出的速度由「水位高度」决定。根据托里拆利定理,流出速度为 $v = \sqrt{2gh}$。水位 h 越高越快,越低越慢。试着把上方模拟器的「初始水位 h₀」调高再按「排水动画」,起初的势头完全不同。
🙋
真的!可是,为什么会变成和下落物体速度相同的公式呢?$\sqrt{2gh}$ 是自由落体的公式吧。
🎓
很敏锐。其实把伯努利方程应用于「水面」和「孔口出口」两点,就会发现位能 $gh$ 直接转化为动能 $v^2/2$。所以会变成和自由落体相同的公式。从能量守恒的角度看,它等于下落了高度 h 的水的动能。
🙋
原来如此!可是把「流量系数 Cd」从 0.62 调到 1.0 后排空时间明显变短了。Cd 是什么?
🎓
Cd 是实际流量与理论流量之比。对锐缘孔口而言,流动在出孔后会猛地收缩,形成「缩流」,使有效截面减小,再加上粘性摩擦。所以现实中的 Cd 约为 0.6。反之换成光滑的喷嘴可升到 0.95 以上。在实务中设计排水管路时,正确估算这个 Cd 决定了排空时间预测的精度。

常见问题

控制方程 A dh/dt = -Cd·a·sqrt(2gh) 是可分离变量的微分方程。两边除以 sqrt(h) 后积分,sqrt(h) 成为时间的一次函数,因此 h 以时间的二次函数减小。图形之所以是下凸的抛物线就是这个原因:满水时斜率(流量)最大,越接近排空斜率越平缓。
排空时间与孔口面积 a 成反比。直径增大一倍,面积变为四倍,因此排空时间约缩短为四分之一。在本模拟器中把孔口面积从 10cm² 改为 40cm²,可以确认排空时间约缩短为四分之一。
不是。本模拟器以横截面积 A 不随水位变化的圆柱(棱柱)水箱为前提。对圆锥水箱或球形水箱而言,横截面积成为水位的函数 A(h),因此需要重新对控制方程 A(h)·dh/dt = -Cd·a·sqrt(2gh) 积分,水位随时间的变化不再是二次函数。
托里拆利定理中的 h 是「从自由表面到孔口中心的高度」。当孔口位于侧面中部时,水位下降到孔口高度后流出就会停止。本模拟器处理孔口位于水箱底部这一最基本的情形,因此 h 与从水箱底部测得的水位一致。

实际应用

给排水设备设计:直接用于估算蓄水箱、游泳池、净水厂沉淀池等的排空时间。预测因维护或清洗而排空水箱所需的时间,并据此确定排水阀和排水管的尺寸。

水位计量与流量计:通过孔口或喷嘴的流动与水位(差压)之间的关系,正是孔板流量计、文丘里流量计的原理本身。托里拆利定理为这些流量计量仪器提供理论基础。

化工厂的批次操作:用于计算将液体从反应槽或储罐输送到下一工序所需的时间。理解非定常流出的特性,对批次工艺的排程和管路设计很重要。

水力模型实验与教学:托里拆利定理是流体力学入门中最先学习的非定常问题之一。作为能够一次性体验伯努利方程、连续性方程和微分方程求解的教材,在许多大学的水力实验中都有涉及。

常见误解与注意事项

最常见的误解是误以为流出速度依赖于水箱横截面积 A。托里拆利定理 $v = C_d\sqrt{2gh}$ 中并不出现 A。流出速度仅由水位 h 决定,无论是粗水箱还是细水箱,在相同水位下都以相同速度流出。A 影响的是「排空时间」,横截面积越大排空越久。请在模拟器中只改变 A,确认流出速度 v 不变。

其次是试图用「除以平均流量」的简单计算来敷衍的错误。由于流量随水位时刻变化,是非定常问题,用满水时的流量去除总水量会大大低估排空时间。正确做法是使用求解微分方程得到的 $t_{empty} = (A/C_d a)\sqrt{2h_0/g}$。事实上,可以从公式确认正确的排空时间恰好是「用满水流量去除」所得值的两倍。

最后是把流量系数 Cd 当作 1.0 处理的陷阱。直接使用理论公式 $v = \sqrt{2gh}$ 会比实际高估 30~40% 的流量。对锐缘孔口而言,缩流和粘性使 Cd ≈ 0.62,而现实管路中还要加上接头和阀门的损失。设计中务必使用实测或可靠文献值的 Cd,并预留保守的余量。