FLUENT水锤模拟和验证

10月前浏览863

正文共：1404字 8图预计阅读时间：4分钟

1 前言

水锤和汽锤都是管路系统中常见的现象，他们均会对管路的长久安全运行造成威胁。“所谓水锤是在突然停电或者在阀门关闭太快时,由于压力水流的惯性，产生水流冲击波，就像锤子敲打一样，所以叫水锤。水流冲击波来回产生的力，有时会很大，从而破坏阀门和水泵”^[1]。之前我们做了几个关于汽锤的案例，今天我们继续做一个水锤的案例。

2 建模与网格

建立如下的二维轴对称管道模型，管道内径20mm，长度15.22m。划分四边形网格，单元数16500（15×1100），控制y+＜5。

3 求解设置

启用SST k-Omega湍流模型。

复制材料库中的液态水材质，并将密度设置为可压缩液体模型。可压缩液体模型方法解释详见帮助文档（以下截图），本案例就不做翻译了，详细读者应该都能看得懂。

入口设置为压力入口，全压为449988.56Pa（水头46m）。

出口设置为速度入口，参考之前汽锤的模拟方法，采用流量的突变来模拟阀门关闭，因为阀门内部的流动细节对于水锤模拟式可以忽略的。采用如下UDF定义边界，速度初始值为0.423m/s，阀门关闭耗时0.018s，阀门关闭特性为线性曲线。实际上对于可压缩流动，采用速度边界并不是很合适，采用质量边界更好。但是有研究表明，这两种边界对于水锤的模拟结果不敏感^[2]。

#include "udf.h"DEFINE_PROFILE(out_velocity,t,i){real velocity;real the_current_time;face_t f;the_current_time=CURRENT_TIME;if(the_current_time>=0&&the_current_time<0.018){velocity=-0.423*(1-pow(the_current_time/0.018,1));}else{velocity=0;}begin_f_loop(f,t){F_PROFILE(f,t,i)=velocity;}end_f_loop(f,t)}