《Water》：Elastic Water Column Modelling of an Undulating Pipeline Profile with Entrapped Air
编辑推荐：
中文翻译摘要：
排放具有起伏轮廓的压力管线时，若未设置进气阀，可能产生低于大气压的压力条件，从而危及结构完整性。迄今为止，针对该瞬变过程的一维模型均依赖于刚性水柱（RWC）假设，该假设忽略了液相弹性与管壁变形能力。本研究将弹性水柱模型（EWCM）应用于起伏轮
中文翻译摘要：
排放具有起伏轮廓的压力管线时，若未设置进气阀，可能产生低于大气压的压力条件，从而危及结构完整性。迄今为止，针对该瞬变过程的一维模型均依赖于刚性水柱（RWC）假设，该假设忽略了液相弹性与管壁变形能力。本研究将弹性水柱模型（EWCM）应用于起伏轮廓管线的排空过程，其中特征线法在两个独立移动网格上求解，并通过一个服从多方过程定律的单一截留空气团耦合。该模型经过了实验室设施中九种构型（每种记录两次）的验证，设施的两条支管以30°倾角上升至共同高点，形成倒V形。构型涵盖了6%开度与全开度下的对称空气团，以及全开度下的非对称空气团分布。EWCM预测最小压力的平均绝对误差为0.21%，而刚性公式为0.27%，两者均低于实验重复之间的离散度。主导参数取决于阀门开度水平：在全开度下，阀门操作曲线的形状起主导作用；而在限制性6%开度下，该形状不影响最小压力，多方指数成为影响最大的因素。波速在整个过程中可忽略不计。声学时间尺度与质量振荡时间尺度之比（可解析求得）在测试构型中介于0.014至0.028之间，两种公式的差异随该比值的增大而增大。
论文解读文章：
研究背景与问题：给水管网、灌溉系统和排水设施在调试、维护和修理过程中，充水与排空操作十分常见。排空过程会使截留空气团膨胀，绝对压力降至大气压以下，可能导致管道压瘪破坏。现场证据表明此类事件并非仅存在于实验室条件，运行系统中的突然减压可产生严重瞬变。起伏管线中的中间高点常是空气聚集之处，当管线排空且高点未设置进气装置时，截留空气团因两侧水柱后退而膨胀，产生的负压波谷决定了操作的安全性。以往针对截留空气瞬变的一维模型主要采用刚性水柱（RWC）公式，将液体视为不可压缩，忽略液弹性和管壁变形；弹性水柱（EWC）公式虽保留液体可压缩性，但大多针对充水过程。目前尚无针对“中间高点处单一截留空气团、两条支管独立排空、各自阀门操作”这一完整构型的弹性描述。为填补这一空缺，本论文将弹性水柱描述应用于具有中间高点的起伏管线排空问题，并通过实验数据验证模型的有效性。
研究内容与结论：研究人员构建了弹性水柱模型（EWCM），对每条支管采用特征线法（MOC）在两个独立变换网格上求解，两条水柱通过共享的截留空气团耦合，空气团遵循多方过程定律。模型在倒V形实验室设施中进行了九种构型（每种重复两次）的验证，涵盖6%与100%阀门开度、对称与非对称初始空气团分布。结果表明，EWCM预测最小压力的平均绝对误差为0.21%，RWC公式为0.27%，两者均小于实验重复之间的离散度（平均82.8 mm水柱）。模型与实验的均方根误差为69.8 mm，同样低于重复实验间82.8 mm的离散度。主导参数取决于阀门开度：全开度下阀门操作曲线形状影响最大，可使预测最小压力变化85.7至263.2 mm水柱；6%限制性开度下该影响消失，多方指数成为最大效应但影响仍小（7.0 mm水柱）；波速在所有构型中可忽略不计（最大影响0.97 mm水柱）。研究还导出了声学时间尺度与质量振荡时间尺度之比的闭式表达式，该比值在全部测试构型中介于0.014~0.028之间，约比两种公式可被测量区分的阈值（约0.11）低四倍。该比值可用于判断弹性效应的重要性，无需运行模型即可由几何参数、波速和初始空气团尺寸计算。
主要关键技术方法：研究人员在实验室设施中开展实验，该设施由两条63 mm内径透明PVC支管构成倒V形，两分支分别从两端水平段以30°倾角上升至共同高点，支管长度分别为3.70 m和3.50 m，总长7.38 m；高点处安装压力变送器（WIKA S-10，0~25 bar绝对量程）与PicoScope数据采集系统，采样间隔0.0062 s。数值方法上，采用特征线法（MOC）在Chaiko和Brinckman变换的归一化移动网格上求解水锤方程，每条支管独立网格，通过空气团多方方程耦合；阀门边界条件采用二分开方法求解含正负流量根的无量纲阀方程；时间步由两条支管中更严格的Courant条件决定，Courant数取0.80。多方指数经Run 4标定为1.1，其余八个构型使用同一值验证。
研究结果：弹性水柱模型成功再现了三种类型构型的压力演化。Type A（6%开度，对称空气团）呈现缓慢单调衰减、无振荡；Type B（100%开度，对称空气团）呈现持续振荡；Type C（100%开度，非对称空气团）呈现尖锐初始下降随后振荡，且随初始空气团非对称程度增大而更快衰减。模型在九种构型中最小压力预测的平均绝对误差为0.21%，其中四个构型优于RWC公式、四个略差、一个相同；最大差异出现在Run 9（最非对称构型），EWCM误差从0.59%降至0.24%。在Type A构型中，限制性开度使阀门开度曲线形状对最小压力无影响，预测值完全由热力学平衡态决定。波速敏感性分析显示，波速在250~1400 m/s范围内变化仅使最小压力变化1.0 mm水柱，证实系统处于准刚性状态，即压力波在单次质量振荡周期内多次穿越水柱，液体弹性对空气团动力学无显著贡献。多方指数对振荡相位高度敏感，1.1的标定值介于等温和绝热极限之间，符合长振荡周期内热交换显著的实际过程。网格收敛性测试表明，N=30与N=100网格间最大偏差仅0.41 mm水柱，收敛快速。时间尺度比在闭式线性化推导下，Run 4预测振荡周期0.600 s（实测0.609 s），全部构型比值介于0.014~0.028；当通过波速或初始空气团尺寸扫参使该比值增大时，两种公式预测的最小压力差异随之增大，在比值约0.11时达到压力记录分辨率3.9 mm水柱。数值扫参还表明该比值并非精确标度律，仅作为弹性效应重要性的体制指示器。
讨论与结论总结：讨论部分指出，EWCM与RWCM在本实验范围内保持紧密一致，符合Malekpour和Karney关于充水前沿可视为刚性柱的发现、Koppel和Laanearu将刚体解作为微扰首阶解的解析，以及Laanearu等人大型排空实验中水锤影响可忽略的观察。时间尺度比优势在于可预先由管道几何、波速和空气团尺寸闭式计算，无需运行模型。摩擦采用准稳态Darcy–Weisbach系数，未引入非稳态摩擦修正，但模型仍能再现振荡阻尼包络，残余偏差可能部分源于此简化。验证域受限于单座7.38 m倒V形设施、单一管材管径、单一几何布置，定量结果不应直接外推至其他轮廓、多空气团、管材或现场规模系统；多方指数仅基于Run 4标定，阀门操作曲线未直接测量，其影响在全开度下占主导，因此需要独立设施数据进行更广泛验证。研究结论可归纳为：1）瞬变幅值的主导参数取决于阀门开度水平——全开时阀门操作形状主导，6%限制开度时其影响消失，多方指数成为最大效应但所有效应均小；2）本设施中弹性与刚性描述等效，液体弹性对最小压力的影响仅0.46 mm水柱，低于测量分辨率，二者残余差异反映数值实现而非弹性描述本身；3）时间尺度比可由几何、波速和初始空气团尺寸闭式求解，两种公式差异随其增大，本实验比值范围0.014~0.028，约为可区分阈值的四分之一；4）空气团在两支管间的分布不控制瞬变幅值，但决定响应特征——限制开度下单调排空、对称构型下持续振荡、非对称构型下初始陡降后重分配；5）模型均方根误差（69.8 mm）低于实验重复间离散度（82.8 mm），剩余偏差处于实验可再现范围内；6）结果基于单一倒V形几何、实验室规模、单一截留空气团和经一个构型标定的多方指数，仍需独立设施验证，并实验确定两种描述开始可分时的比值。
