突然关阀，水为什么会冲一下？
管里的水原本正在流。阀门突然关上时，附近的水先慢下来，远处的水却还在往前走。这个变化需要沿着管道一段段传开，不能让整管水同时停住。
水会微微被压缩，管壁也会微微鼓开，压力变化就这样向远处传播。这种快速的压力变化叫水锤，它还会在管道两端反射回来。
跟住颜色变化和箭头：颜色表示压力的变化，箭头表示水向哪里流。画面把形变和时间放大，帮助你看清；手机上竖着排的是同一根水平管。
在画面里比较快关与慢关
- 先快关。 看靠近阀门与远离阀门的水，哪边先改变？
- 再慢关。 其余条件相同，比较两次最大的压力变化。
- 跟住反射。 暂停、倒回，看变化怎样传到水库，又怎样回来。
用动画做比较即可，不需要在家反复猛关水龙头。
两次反射和一个时间尺度
展开：更多对照与操作说明
- 关闭时间。 比较 0.02 s 与 0.90 s，两次保持同一管长、管径、初速、压力和开度曲线形状。看压力峰值怎样改变，不把手柄位移直接当作流量。
- 管壁刚度。 从 200 GPa 降到 20 GPa，波速和第一次快速关闭的峰值都变小，往返时间变长。这里改变的是理想线弹性模量，不模拟真实塑料的黏弹性。
- 初始流速。 首次快速停止产生的压力增量与初速有关；零初速时，关闭没有产生水锤波。
- 初始绝对压力。 降低它，返回的降压波更容易碰到蒸气压边界。模型在这里停止，不继续画出不可信的负绝对压力。
- 时间与测点。 在同一次试验中前后拖动时间，或沿管移动测点。查看同一位置的流速、绝对压力与整个管路状态如何对应。
每次更改装置条件都会从同一规定初态重新计算。没有弹性支架、管网支路、气罐或真实阀门曲线被暗中加入。
水库和阀门为什么反射出不同结果？
恒压水库固定上游压力，因此到达它的压力扰动必须与反射波相消。压力扰动反号返回，流速却可以反向。关闭的末端阀门约束的是零流速，所以那里反射的压力扰动保持同号，用来消除到达边界的流速变化。
“同号”是相对于到达阀门的那次扰动，不代表阀前压力永远升高。快速关闭的初始升压波在水库处变成降压波，回到已关阀门后可使阀前压力低于初始值。
波速为 a，管长为 L。单程时间为 L/a，回到阀门需要 2L/a。默认参数下分别约为 0.089 s 和 0.179 s。演示把这些亚秒过程放慢，显示的物理时间来自求解器。
独立模型具体求什么？
以 p′ = p − p₀ 表示压力增量，u 表示轴向流速，模型求解：
∂p′/∂t + ρa² ∂u/∂x = 0
∂u/∂t + (1/ρ) ∂p′/∂x = 0
合成顺应性为 C = 1/K + cD/(Ee)，波速 a = 1/√(ρC)。水的体积模量 K 描述压缩，管壁弹性模量 E、管径 D、壁厚 e 和约束系数 c 描述壁面顺应性。作者取 c = 1，不求解轴向管壁运动或泊松耦合。USACE HEC 的有压管流技术说明给出这些传播方程及管壁弹性修正。
作者参数：L = 120 m，D = 0.08 m，e = 0.004 m，ρ = 998.2 kg/m³，K = 2.2 GPa。默认 E = 200 GPa，初速 1 m/s，上游绝对压力 1.8 MPa，下游绝对压力 0.101325 MPa。它们是教学条件，不是现场管道的设计或耐压评级。
忽略管摩阻、对流项、气体、黏弹性、支架振动和非线性材料效应。由 u ± p′/(ρa) 沿 x ± at 不变，采用 96 段特征线网格，满足 aΔt/Δx = 1。管内改变只能经过有限时间传到别处，没有全管同步停流的捷径。
阀门并不按指定流量硬关
有效通流系数比例 τ 从 1 平滑降到 0：关闭期间 τ = 1 − 3r² + 2r³，r = t/T꜀。图中的矩形阀口表示这一指定系数比例，不代表任意阀门的真实转角或升程关系。孔口规律同时依赖压差：
边界把它与到达阀门的特征关系联立求解，允许仍开启的阀门在压差反向时出现反向流动。完全关闭后直接施加 u = 0。这里没有把阀门开度百分比误标成流量百分比。
焦科夫斯基关系什么时候成立？
简单管路中，在返回波影响前快速改变流速，首个压力增量满足 Δp = −ρaΔu。完全停止时，升压幅度为 ρaU₀。缓慢关闭使反射在阀门尚未关完时就返回；本场景的平滑孔口规律因此给出更低的峰值。
关闭时间和 2L/a 的比较很有用，但不是一条涵盖所有阀型、气体、管网和空化过程的安全公式。尤其在液柱分离与汽腔塌陷后，压力峰值可能超过简单的首波值。Simpson 与 Wylie 的实验和分析专门研究了这种情况。本模型遇到蒸气压即停止，不声称预测了这些后续冲击。
能量、体积与汽化边界怎样核对？
线性系统的动能密度为 ρu²/2，增量弹性能密度为 Cp′²/2；其中水和管壁分别占 (1/K)p′²/2 与 [D/(Ee)]p′²/2。瞬时关闭、零摩阻、恒压水库和闭阀边界下，离散总能量守恒。有限时间关闭时，边界可以交换增量能量，不能把所有减少都叫作管壁发热。
压缩与扩张的库存增量为 A∫Cp′ dx，随时间的变化等于入口体积流量减去出口体积流量。测试独立核对该关系，并用 32、64、128 段解与 1024 段参考解比较，检查渐进关闭后的波形收敛。
20 °C 水的蒸气压取约 2.338 kPa。每一步先检查首次触及这一绝对压力的位置，插值到边界并冻结最后有效状态，之后不继续传播单相解。Bergant、Simpson 与 Tijsseling 的液柱分离研究综述解释了为什么汽腔与重新撞合需要额外的物理模型。
主图管壁形变放大 500 倍，局部形变放大 800 倍，数据仍显示真实线性模型的微小应变。时间插值只用于连续呈现；数值验证在求解网格上完成。