本文介绍了如何在双曲面肖伯格漏斗中产生三种不同的水涡流态,阐述了这些流态最重要的特征,并说明了如何计算相关参数,例如氧气传递速率。
本文介绍了如何在双曲面肖伯格漏斗中产生三种不同的水涡流态,阐述了这些流态最重要的特征,并说明了如何计算相关参数,例如氧气传递速率。
自由液面涡流在工业中广泛应用于流量调节、能量耗散和能量生成。尽管已有大量研究,但关于自由液面涡流的详细实验数据仍然缺乏,尤其是涉及界面湍流方面的数据。本文报道了一种由沃尔特·绍伯格(Walter Schauberger)在20世纪60年代首次提出的特殊类型自由液面涡流,其氧体积传质系数超过类似系统的数值。这种特殊类型的涡流在双曲漏斗中形成。通过不同的水力特性可稳定多种不同的流动状态。该技术的其他优势包括高能效、结构简单以及可扩展性。双曲漏斗内的流动具有强烈湍流特征,并显著增加了气-水界面的表面积。表面局部压力变化显著,导致明显的波浪状气-水边界层。由于螺旋流动的存在,这些扰动向内移动,并带动边界层一同运动。由此产生的压力梯度将一定体积的空气吸入水涡中。本文展示了基本双曲漏斗装置的构建方法及运行实例,包括针对三种不同稳定状态的高速可视化结果。
我们的生活与螺旋结构密切相关。螺旋几乎存在于万物之中,无处不在,例如贝壳和菊石的结构,以及飓风、龙卷风和漩涡的形成1,2。在宇宙学尺度上,星系的形成与演化遵循对数螺旋的原理3。其中最著名的螺旋是黄金螺旋和斐波那契螺旋4,其应用范围广泛,从描述植物生长和某些固体的晶体结构,到开发计算机数据库搜索算法。斐波那契数列是一个以0和1开始的数字序列,其后续每一项均为前两项之和。在计算兔子繁殖速率时也能发现这一数列。螺旋是智人(Homo sapiens)绘制的最古老的几何图形之一,例如在哥伦比亚和澳大利亚发现的同心圆(公元前40,000–20,000年1)。列奥纳多·da Vinci5曾尝试利用螺旋叶片(源自希腊语 ἕλιξ πτερόν,即 helix pteron,意为螺旋翼)设计一种直升机形状的飞行器。遵循同样的原理,飞机设计师 Igor Sikorsky 在450年后制造出了第一架投入系列生产的直升机6。
许多其他实例表明,螺旋流动结构可能非常高效且节约成本,因为这种流动形式在自然界中普遍存在。20世纪初,奥地利林务员兼哲学家维克多·绍博格(Viktor Schauberger)认识到了这一点。他提出,人类应当研究自然并向其学习,而非试图改造自然。基于这一理念,他建造了相当独特的原木滑道用于木材运输;这些滑道并未采用两点之间最短的直线路径,而是沿山谷和溪流的蜿蜒路线延伸。这种设计使水流沿轴线呈螺旋状旋转,从而形成涡流,不仅减少了用水量,而且运输速率显著超过了当时认为的正常水平7。
受其父亲影响,维克托的儿子瓦尔特(Walter)利用水涡流8开发了多种新技术,应用于饮用水处理、工业流程、池塘与水道的修复、池塘及小型湖泊的增氧,以及河流的调控与生态恢复。其中一项技术近年来受到广泛关注,即利用双曲漏斗8进行水处理;该装置仅依靠水流本身形成涡流,无需任何搅拌装置。该方法已被证实对地下水中铁的氧化具有显著效果9,10。然而,该技术的局限性在于对低pH值的水处理效率较低11。
荷兰的大量饮用水来自地下水源12,其中铁浓度可高达每升数十毫克13,而根据标准,可接受的限值为 0.2 mg/L14。大多数饮用水处理厂在水处理过程的初期阶段采用曝气法以降低水中铁的浓度。在大多数情况下,曝气的目的是增加溶解氧含量,去除水中的气体及其他相关物质,或同时实现这两项目标15。曝气将氧气引入液体介质的方法有多种,包括使用搅拌器或涡轮搅拌液体表面,以及通过宏观孔口或多孔材料释放空气16。
van de Griend17证明了铁氧化的化学过程,其中氧分子从亚铁离子夺取一个电子,并与一个游离质子反应生成水,同时铁离子被氧化(方程[1]):
, (1)
铁离子随后沉淀为 铁(OH)3 因其与水反应会释放质子(方程式[2]):
(2)
总反应由方程(3)给出:
. (3)
在曝气过程中,最常采用的技术包括跌水、塔式、喷淋和板式曝气系统18,19。这些技术的缺点是其能耗占污水处理设施总能耗的50%至90%20,并占设施运行与维护预算的最高达40%21。
使用双曲面漏斗进行曝气可显著降低该过程的成本并提高效率。由于双曲面漏斗的几何结构以及无运动部件,其对堵塞的敏感性较低,这意味着能量仅消耗于水泵运行。此类系统可通过多个参数进行表征,例如漏斗每小时的水流量(φ)、平均停留时间(MRT)、水力停留时间(HRT)、氧容积传质系数(KLa20)(校正至标准温度20°C)、标准氧转移速率(SORT)以及标准曝气效率(SAE)。漏斗的流量用于计算在特定时间内可处理的水体积。MRT通过在特定工况下漏斗中水流量与其体积之比,利用公式(4)计算得出:
(4)
其中 V 表示反应器中的液体体积。
HRT 可通过示踪技术进行实验测定22 通过 其停留时间分布函数。水力停留时间(HRT)为混合过程、滞留现象及分离行为提供了基本的洞察。23由 Donepudi 证明24 水射流距离入口越远,向出口移动的速度越快。在初始时刻,水沿漏斗上部圆柱段的切线方向被泵入。随后,在重力和系统几何结构的共同作用下,切向速度减小,轴向速度增大。氧传质体积系数 KLa20 (单位倒数时间),表示系统促进氧气向液相传递的能力10可以计算得出25,26 根据公式(5):
(5)
其中 C出 是主体液体中的溶解氧(DO)浓度, C在 是进料中的溶解氧浓度, Cs 是饱和状态下的 DO 浓度,且 T 是水温。
SORT值是指系统向液相传递氧气的标准速率,由公式(6)27确定:
(6)
其中
表示在 20 °C 温度下溶解氧的饱和值。SOTR 值可针对某一特定过程进行定义,此时方程(6)中使用的体积需归一化为 1 小时的处理时间(过程特异性 SOTR),从而实现小试规模曝气方法与实际规模系统的比较。对于漏斗中某一特定工况的能力评估,必须计算系统特异性 SOTR,该值使用漏斗内水体体积及(工况特异性)水力停留时间。该参数在计算给定漏斗中某一工况的实际曝气能力时具有重要意义。
SAE 是 SOTR 与曝气所消耗功率的比值。由于能量仅消耗在将水抽至漏斗顶部并提供形成涡流所需的水流上,因此可通过公式(7)计算,即每小时抽水体积在对应于漏斗长度高度处的势能与水形成涡流所需的动能之和27:
(7)
其中,Pp 是将水泵送至漏斗高度所需克服的势能功率(单位为 kW),Pk 是在漏斗顶部的水流获得足够流速以形成涡流所需的动能功率(单位为 kW)。通常情况下,公式(7)应采用系统特定的 SOTR。若改用过程特定的 SOTR,则得到的是水力停留时间为 1 小时的(理论)系统的能耗。
这些参数足以评估该技术应用的有效性和可行性,但不足以描述该过程本身。需要指出的是,涡旋是流体力学中最不被充分理解的现象之一,因此该领域投入了大量研究工作。在流体力学中寻找涡旋普遍规律的主要挑战之一在于,几何边界条件始终存在变化,这些变化会影响涡旋的发展,并显著影响其形成与动力学行为。因此,可以合理地认为,自由表面涡旋(FSV)不能类比于实验室中受限的涡旋。然而,Mulligan 等人28在泰勒-库埃特流(TCF)的研究中表明,若将自由表面涡旋的气核视为一个以气核相同转速旋转的虚拟内圆柱,则两者可被类似处理。通过这种方法,描述自由表面涡旋流场的方程可被虚拟圆柱的角速度条件所替代,从而转化为适用于TCF系统的方程。此外还证明,当假想圆柱的旋转速度增加到一定程度时,会出现类似泰勒的二次流涡旋结构28,并在接近壁面时消失。
Niemeijr29研究表明,在Schauberger漏斗中可产生三种不同类型的水流涡旋(扭曲形、直线形和受限形)(图1 和 图2),这些涡旋具有不同的水力参数特征。在此基础上,Donepudi24采用了与Mulligan等人28相同的方法,利用计算流体动力学(CFD)模拟涡旋流态,进而分析其流场结构,以揭示其背后的物理机制。该系统具有高度湍流特性,二次流场极不稳定,并表现出大量类泰勒(Taylor-like)涡旋的形成。气体从气相向液相的传输过程受扩散、对流和反应的共同作用。因此,为提高该过程的效率,必须增加气相浓度梯度或增强液体的体积流动。后者直接取决于系统中类泰勒涡旋所体现的湍流强度,这些涡旋有助于将界面处饱和的流体单元输运至主体液体中。在另一项相关研究9中,比较了不同涡旋流态下的主要参数,如水流速率、KLa20 和 SOTR。该研究表明,相较于其他用于水体曝气的方法,该技术在气体传递速率方面展现出显著优势,具有广阔的应用前景。
本文的目的是提供并演示一种在双曲面绍伯格漏斗(小型:高26 cm,顶部直径15 cm;中型:高94 cm,顶部直径30 cm;大型:高153 cm,顶部直径59 cm)中产生不同水涡流态的方法,以实现高效的水体增氧。
1. 一般建议
2. 实验装置
3. 操作(中号漏斗)
在施陶格伯格双曲面漏斗中,水漩涡在不同工况下(扭曲态、直线态和受限态)形成(图1)。由此,水体被大气中的氧气富集,促进了水中化学物质的氧化。该系统除将水抽送至双曲面漏斗上部所需能量外,无需额外能耗。
扭曲流态具有双螺旋形状,且水与空气之间的界面最大。形成该流态需要施加中等水流速度(小漏斗为75-78 L/h,中漏斗为1,194-1,218 L/h,大漏斗为4,834-5,032 L/h)。其在漏斗上部圆柱段的高度不应超过小漏斗2 cm、中漏斗7 cm、大漏斗16 cm。
直线流态呈平滑的直线状,水与空气之间的界面较小。该流态需要最大的水流速(小型漏斗为93-100 L/h,中型漏斗为1,314-1,338 L/h,大型漏斗为5,102-5,289 L/h)。对于所有漏斗,其水流高度均可达到顶部盖板。
根据水位的不同,受限流态可呈现为扭曲状或直线状的涡流。然而,该流态的特点在于其长度会随着背压的施加而发生变化,这与之前不施加压力的各种流态不同。该流态在漏斗顶部形成;但随着背压的增加,其尾部开始缩短,涡流逐渐从底部消失。其水流量极小(小漏斗为58–70 L/h,中漏斗为882–936 L/h,大漏斗为2,351–2,634 L/h),其高度可根据漏斗几何形状达到最小或最大值。
不同的流态可通过调节水流速率、背压和系统几何结构实现稳定并相互转化。表征曝气效率的参数包括水流速率、氧的体积传质系数和标准氧转移速率。可以看出,在低水流速率的扭曲涡流条件下,KLa20 最高(图4),比K高出数倍La20 在直线和受限流态下,该指标比用于湖泊和河流曝气的常规系统(如空气喷射、叶轮、桨板)高出数十倍,而这些常规系统能耗更高。随着水流进一步增加,KLa20 逐渐减少,但水位,即系统中的水量,却增加了。当超过某个阈值后,螺旋流态转变为直线流态。对于每种流态,均存在体积和水力参数保持不变的稳态条件。
然而,在比较小、中、大三种漏斗的相似条件时,各系统在水流速率和体积上的差异显著。但与此同时,K 的比值La20 数值变化不大。83 h 的最大值-1 对于小型漏斗,60 小时-1 用于中号漏斗,以及79 h-1 在扭曲状态下,大漏斗实现了上述结果。
同时,当 KLa20 随水流增加而降低时,MRT 增加,表明水通过漏斗所需时间更长,Donepudi24对此有详细描述。然而,与 KLa20 相似,在不同漏斗中,扭曲流态与直流通态下的 MRT 值大致相同。小漏斗的 MRT 范围为 10 s 至 43 s,中等漏斗为 14 s 至 30 s,大漏斗为 24 s 至 43 s(表 1)。

图1:26 cm 高的玻璃双曲面绍伯格器漏斗中的水涡流态。 (A)扭曲涡流(75 L/h),(B)直行涡流(100 L/h),(C)受限涡流(70 L/h)。 请点击此处查看此图的放大版本。

图2:在94 cm高的玻璃双曲面绍伯格漏斗中形成的水流涡流状态。(A)扭曲状(1,194 L/h),(B)直线状(1,314 L/h),(C)受限状(882 L/h)。请点击此处查看此图的放大版本。

图3:用于协议步骤3.1–3.3中所述实验的装置示意图。(1) 地下水泵;(2) 控制阀;(3) 水流流量计;(4, 5) 用于溶解氧检测的聚合物光纤;(6, 7) pH探针;(8) 温度传感器;(9) 含示踪剂的注射器;(10) 阀门;(11) 绍伯格双曲面漏斗;(12) 光纤氧信号发射器;(13) 笔记本电脑;(14) 数据记录仪;(15) 夹具;(16) 排水口。 请点击此处查看该图的放大版本。

图 4:大型漏斗装置的照片。 (1)地下水储槽;(2)水泵;(3)水流计;(4)含示踪剂的注射器;(5, 6)带有氧传感器斑点的玻璃接头;(7)、(8)pH 探针;(9)绍伯格双曲漏斗;(10)排水口。 请点击此处查看此图的放大版本。
| 漏斗 | 流型 | φ (L/h) | HRT (s) | MRT (s) | KLa20 (h-1) | V (L) | 液位 (cm) | Cin (mg/L) | Cout (mg/L) | SOTR (g O2/h) | SAE (g O2/kWh) |
| 小号 | 扭曲型 | 75 | 10 | 15 | 83 | 0.2 | 0.5 | 0.0 | 1.8 | 0.2 | 2801 |
| 78 | 20 | 24 | 41 | 0.4 | 2 | 0.0 | 1.9 | 0.2 | 2932 | ||
| 直通型 | 93 | 31 | 24 | 25 | 0.8 | 4 | 0.0 | 1.8 | 0.2 | 2688 | |
| 100 | 43 | 32 | 18 | 1.2 | 6 | 0.0 | 1.7 | 0.2 | 2635 | ||
| 限流型 | 58 | 18 | 23 | 14 | 0.3 | 1 | 0.0 | 0.6 | 0.0 | 872 | |
| 70 | 53 | 31 | 2 | 1.0 | 5 | 0.0 | 0.3 | 0.0 | 459 | ||
| 中号 | 扭曲型 | 1194 | 14 | 13 | 60 | 4.5 | 2 | 0.0 | 1.8 | 2.4 | 784 |
| 1218 | 19 | 19 | 37 | 6.3 | 5 | 0.0 | 1.6 | 2.1 | 667 | ||
| 直通型 | 1314 | 29 | 29 | 18 | 10.7 | 11 | 0.0 | 1.2 | 1.8 | 509 | |
| 1338 | 30 | 31 | 18 | 11.0 | 11.7 | 0.0 | 1.2 | 1.8 | 500 | ||
| 限流型 | 882 | 21 | 24 | 17 | 5.1 | 3 | 0.0 | 0.9 | 0.8 | 348 | |
| 936 | 37 | 36 | 5 | 9.7 | 9 | 0.0 | 0.5 | 0.4 | 180 | ||
| 大号 | 扭曲型 | 4834 | 24 | 23 | 79 | 32 | 11 | 0.7 | 4.1 | 22.9 | 1113 |
| 5032 | 34 | 26 | 52 | 48 | 16 | 0.6 | 4.0 | 22.6 | 1054 | ||
| 直通型 | 5102 | 38 | 29 | 31 | 54 | 19.5 | 0.7 | 3.0 | 15.0 | 690 | |
| 5289 | 43 | 30 | 19 | 64 | 22.5 | 0.6 | 2.3 | 10.8 | 479 | ||
| 限流型 | 2351 | 58 | 43 | 16 | 38 | 7 | 0.6 | 2.6 | 5.5 | 557 | |
| 2634 | 95 | 50 | 7 | 70 | 19 | 0.6 | 2.0 | 4.2 | 380 |
表1:小型(图1)、中型(图2)和大型漏斗的基本水力特性与曝气效率参数。
如果地下水泵功率过大而系统无法承受压力,可在控制阀前增加一个排水口以降低压力。为了获得可靠的结果并确保示踪实验使用快速探头,校准传感器非常重要。如果探头响应过慢,则会扭曲水力停留时间(HRT)的测量结果。此外,若在直线流态下,HRT 显著小于平均停留时间(MRT),这可能表明漏斗的切向入口明显低于水位,示踪液进入漏斗后有一部分直接流向排水口,从而导致 HRT 降低。
双曲沙伯格漏斗中的水涡旋对水流速率非常敏感。系统越小,对流量变化的依赖性就越大。如果流动状态稳定,漏斗中的水位应不随时间变化。若情况并非如此,水位将会上升或下降。因此,应注意观察水位,以避免发生水溢出、因漏斗内部压力升高而导致的裂纹,或出现非预期的流动状态转变。
为了 确定涡流的状态(方案步骤 3.1.3.1–3.1.3.3)及其稳定性,使用透明的漏斗较为有利。因此,本研究中采用了玻璃漏斗 。在运输、操作和安装过程中必须格外小心,应注意不要过度拧紧盖子上的螺丝,以免造成损坏(方案步骤 2.1.2)。
为了确定水力停留时间(HRT),应尽可能多次重复步骤3.2.2–3.2.3(至少10次),因为系统存在 高度湍流及次级流动(类似泰勒涡流),示踪剂射流可能分离并沿不同路径通过漏斗。例如,Donepudi 等人24 和 Mulligan 等人28 的研究表明,水层越接近玻璃壁,其流向排水口的速度越快。探头每次使用后均应用去离子水清洗并擦拭,以避免样品与保存液混合,否则可能破坏数据并降低电极保存质量。
进行 DO 实验时,需确保系统出口处的氧气浓度值稳定(实验步骤 3.3.2.2)。若系统尚未达到稳定状态但波动不显著,则应对测得的数值进行平均处理。此外,盖子上必须留有通气孔,以使空气流入系统,实现进一步曝气。
尽管该系统的KLa20 值和能量效率较高,但由于现有漏斗的水流速率较低,其SOTR值相较于其他方法仍较低26,这目前限制了双曲面漏斗在工业水体曝气中的应用。然而,已有研究表明,通过使用大、中、小型不同尺寸的漏斗,可在不同规模下实现系统的高效率。由此可知,通过改变漏斗的几何结构(如尺寸、进出口直径、壁面曲率),可在不降低曝气效率的前提下,显著提高水处理的速度和处理量。此外,从表1中可以看出,当漏斗长度增加1.1 m时,SOTR提升了100倍以上。考虑到某些水处理厂中水位差可达数米,因此(部分)曝气过程可能以远低于当前的成本实现。因此,明确漏斗不同几何参数对涡流状态下水流速率及KLa20 的影响,将有助于开发一种低成本且具有竞争力的地下水曝气技术。此外,如Schauberger31所示,曝气还可用于改善水库、湖泊和河流的水质。
作者声明,他们不存在任何已知的相互竞争的经济利益或可能被认为影响本文所报告工作的个人关系。
本工作是在Wetsus欧洲可持续水技术卓越中心(www.wetsus.eu)的应用水物理研究主题合作框架下完成的。Wetsus由荷兰经济事务部、基础设施与环境部、弗里斯兰省以及荷兰北部各省共同资助。本研究获得了欧盟“地平线2020”研究与创新计划在玛丽·斯克沃多夫斯卡-居里基金协议No. 665874以及Gilbert-Armstrong实验室的资助。我们非常感谢Maarten V. van de Griend对本工作的支持。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 单通道/双通道变送器 | Endress+Hauser | CM442 | 数据记录仪 |
| 控制阀 | +GF+ | 625DN20 | 型号514 |
| 数据记录仪 | Endress+Hauser | CM442 | Liquiline |
| 光纤氧变送器 | PreSens | SACN0002000005 | Fibox 3 |
| 玻璃弯头连接件 | 定制 | - | 管道适配器 |
| 地下水泵 | SAER | 3637899 | H/150 |
| 笔记本电脑 | 任意品牌 | 任意型号 | Windows 10 或更高版本 |
| 大型玻璃漏斗 | 定制 | - | 高94 cm |
| 氧气计算器 | PreSens | v. 3.1.1 | 软件 |
| 氧传感器斑点 | PreSens | NAU-D5-YOP | SP-PSt3 |
| pH连接件 | 定制 | - | pH探头适配器 |
| pH传感器 | Endress+Hauser | CPS11 | Orbisint CPS11 |
| 聚合物光纤 | PreSens | POF-L2.5-2SMA | OXY-1 SMA |
| 橡胶垫片 | ERIKS | 11535207 | 141×197×2 mm |
| 橡胶垫片 | ERIKS | 12252766 | 273×340×3 mm |
| 小型玻璃漏斗 | 定制 | - | 高26 cm |
| 水流量计 | Endress+Hauser | P7066819000 | Picomag |
| 水流量计 | Kobolt | 5NA15AC34P | MIK |
| 水温连接器 | PreSens | - | Pt100 |