管路·流体条件
— 判定中 —
阀门位置处的压力-时间波形（简化矩形波近似）。红虚线 = 许可压力（PN16=1.6MPa）参考值。
阀门关闭时间 tc 对最大水锤压 ΔP 的影响。竖虚线 = 当前 tc、横虚线 = 临界关闭时间 Tr。
管路长 L 对临界关闭时间 Tr 和最大水锤压的影响（相同 V₀）。
理论·主要公式
急速关闭情况下（tc ≤ 2L/a）的最大水锤压力升高（Joukowsky公式）：
\(\Delta P = \rho \cdot a \cdot \Delta V\)
弹性管内的压力波速度（Halliwell公式）：
\(a = \sqrt{\frac{K/\rho}{1 + (K \cdot D)/(E \cdot e)}}\)
其中 K = 液体的体积弹性系数，E = 管材的杨氏模量，D = 内径，e = 管壁厚。
临界（往复）时间和缓慢关闭时的水锤压（Allievi近似）：
\(T_r = \frac{2L}{a}, \quad \Delta P_{slow} \approx \Delta P_{Joukowsky} \cdot \frac{T_r}{t_c} \quad (t_c > T_r)\)
水锤的原理 — 对话理解
🙋「水锤」是指水龙头关闭时发出的"砰！"的声音吧？但为什么在工厂配管中会成为问题呢？
🎓是的，那就是水锤。家用水龙头虽然只是发出声音，但工厂和发电厂的大口径配管流速很大，若阀门急速关闭就会产生数MPa的压力冲击波。这会导致配管破裂、接头损伤和泵反转等严重事故。历史上有多个水力发电站和化工厂因此发生过重大事故。
🙋Joukowsky公式 ΔP = ρaΔV 中的"a"是什么？就是音速吗？
🎓接近水中的音速，但准确说是"管内压力波速度"，它同时受液体的体积弹性系数K和管壁弹性E的影响。纯水中的a约为1480 m/s，但钢管内由于管壁弹性会减慢，大约1200-1400 m/s。橡胶管或HDPE管中更慢，约200-500 m/s，因此水锤压力也会降低。
🙋什么是"临界关闭时间"？只要超过这个时间慢慢关闭就没事了吗？
🎓Tr = 2L/a 是压力波在管路往返的时间。当阀门关闭时间 tc > Tr 时，为"缓慢关闭"，反射波返回前阀门已部分打开，所以水锤压力会降低。当 tc < Tr 时为"急速关闭"，会产生Joukowsky的最大值 ΔP = ρaΔV。例如L=100m、a=1400m/s时，Tr=0.14s。若超过这个时间缓慢关闭就能避免完全的水锤冲击。
🙋防止水锤只需要阀门慢慢关闭就可以了吗？还有其他对策吗？
🎓缓慢关闭是最基本的对策。其他还有膨胀罐（在管路上设置开放式罐体来吸收压力）、空气室（密闭空气罐来缓冲压力冲击）、安全阀/泄压阀（超压时释放）、带飞轮的泵（防止急速停止）等。通常需要综合使用多种对策。
🎓在正压波之后会出现负压波，管内压力可能降到大气压以下。这会引发"液柱分离（column separation）"，液体汽化和空穴化，之后再次碰撞时产生"二次水锤"，这是极其危险的现象。长管路和起伏配管特别容易出现这种情况。
🙋经常听说MOC（特征线法）分析，和这个模拟器有什么区别呢？
🎓这个工具采用Joukowsky公式做"简单的一维分析"。实际设计中使用MOC（特征线法）这种数值解法，逐步计算管路各区间的压力和流速随时间的变化。可以包含分支管、泵和阀门的特性曲线，对整个复杂管网系统进行精确模拟。AFT Impulse和HAMMER（Bentley）等专业软件就是这样的工具。
常见问题
临界关闭时间 Tr = 2L/a。当 tc ≤ Tr 时为急速关闭（产生最大 ΔP = ρaΔV），tc > Tr 时为缓慢关闭（ΔP 降低为 1/Tr×tc 倍）。如果L=200m、a=1400m/s，则Tr=0.29s。超过这个时间缓慢关闭可以抑制最大水锤压力。实际设计中通常以Tr的2-3倍以上作为关闭时间目标。
钢管（E=210GPa）：1200-1400 m/s；铸铁管：1000-1200 m/s；混凝土管：1000-1200 m/s；HDPE管（E=1GPa）：200-400 m/s；橡胶管：100-300 m/s。弹性系数越低的材料，波速越低，根据Joukowsky公式水锤压力也会降低。采用HDPE或可挠管的水锤降低设计也很常见。
当负压波使管内压力低于大气压（绝对压力低于液体蒸汽压）时，液体汽化形成空穴（vapor cavity）。空穴崩溃时会产生"二次水锤"，这是比初始水锤更大的冲击压力。长管路、起伏管路和高速泵停止时特别危险。需要用MOC分析进行包含空穴模型的精确模拟。
泵急速停止时，逆止阀关闭，产生类似阀门关闭的水锤现象。ΔV与泵额定流速相近。防止措施包括飞轮（通过惯性减缓停止）、缓闭逆止阀、旁通管路等。泵保护需要在MOC分析中将止回阀的动特性也纳入计算。
空气室是在管路上安装的密闭容器，内部密封空气。压力波来临时空气被压缩，吸收能量，降低压力峰值。负压时空气膨胀补充液体，防止液柱分离。设计需要优化空气容积、初始压力和安装位置，计算时需考虑压缩气体的状态方程（多变过程）。
在压力容器设计规范（如GB150）和工业配管设计规范中，水锤被视为过渡工况。城市供水规范（CJJ/T 111）对给水管道系统的阀门操作时间有限制。石油行业参考API 651，电力行业参考ASME B31.1等国际规范。设计时需以MOC分析结果的最大压力乘以适当安全系数，选择相应的法兰定级压力（ANSI/GB等级）。
水锤（ウォーターハンマー）分析模拟器简介
本模拟器的物理模型基于非稳定流中水锤现象的Joukowsky基础理论进行分析。当管内流速发生急变时，会产生压力波，最大水锤压 \(\Delta p\) 由 \(\Delta p = \rho c \Delta v\) 给出。其中 \(\rho\) 是流体密度，\(c\) 是压力波速度，\(\Delta v\) 是流速变化量。压力波速度 \(c\) 考虑管的弹性，表示为 \(c = \sqrt{\frac{K}{\rho}} \left(1 + \frac{K D}{E e}\right)^{-1/2}\)。\(K\) 是流体的体积弹性率，\(D\) 是管内径，\(E\) 是管材的杨氏模量，\(e\) 是管厚。当阀门关闭时间 \(t_c\) 短于压力波往返时间 \(2L/c\)（\(L\) 为管路长）时，判定为急速关闭，最大水锤压直接作用。本模型基于这些关系式，从输入参数实时计算最大水锤压、压力波速度、临界关闭时间，为安全设计提供指标。
实际应用案例
工业中的实际应用
在化工厂和石油精炼厂中，高压泵出口阀门急速关闭时产生的水锤压力评估采用本模拟器，基于Joukowsky公式进行实时计算。例如，中国石化和中国国家石油公司的配管设计中，利用此工具评估不同管材（碳钢、不锈钢等）对许可应力的安全裕度，通过调整阀门关闭时间来抑制最大水锤压力，防止法兰接头和膨胀接头的破损风险。
教育研究中的应用
在大学流体工程实验和配管设计课程中作为可视化教学工具被采用。清华大学、浙江大学等机械工程系学生通过改变流速和管径，观察压力波形变化，理解临界关闭时间的概念。相比仅用理论公式，通过交互式操作更容易理解瞬态现象的特性，有效促进主动学习。
CAE分析的联动和实务定位
本模拟器作为详细的一维和三维CAE（如Ansys Fluent和STAR-CCM+的瞬态分析）前期筛选工具。实务中，在设计初期通过本工具快速获取水锤压估算值，缩小危险工况范围，之后进行本格CAE分析，可降低计算成本。通过高效筛选危险条件，从而有效验证工业配管耐压设计中安全系数的合理性。
常见误解和注意事项
常见误解是"阀门关闭时间短于临界关闭时间就一定危险"，但Joukowsky公式以瞬时关闭为前提，即使关闭时间短于临界时间，根据关闭方式（线性或指数）、管路弹性和液体压缩性等因素，最大水锤压也可能低于理论值。另外，"流速越快水锤压越高"的理解是对的，但压力波速（水锤波传播速度）取决于管内径、管材杨氏模量、液体体积弹性率等，即使相同流速，由于管道材料和支承方式不同，水锤压也可能差异很大。需要特别注意的是，"模拟结果安全就意味着实机也安全"的过度信任是危险的。本工具简化模型为单一直管和简单阀门关闭，实际工业配管中分支、弯头、泵的非线性特性和气穴化影响不容忽视，必须在设计中适当考虑安全系数。