桥墩绕流问题计算-北京大学教务部.doc

桥墩绕流问题计算 Computation on Flow Around the Pier 北京大学 力学与工程科学系 工程结构分析专业 00级 王嘉 Tel：86-10Email：kingadd@ 指导老师：陈耀松 李植 [摘要] 本文研究现实桥墩绕流问题。该问题的难点在于必须考虑自由面。本文突破处理自由面的传统方法（VOF），将问题分解为大小两个范围的层次，并在边界衔接的方法来处理：对大范围潮流采用长波理论，流动参量用沿水深的平均值代替，这就是Boussinesq方程，对此已不必考虑自由面；在桥墩或其他障碍物附近的局部流动，特别靠近底部，海底和桥墩必须同时考虑，此处的水涡运动是三维的，它对工程基础的破环作用甚大，必须按三维如实计算。只是在这第二个层次中我们可以规定一个计算区，以第一层次算得的压力和水平速度为边界值，用三维的粘性动力学方程进行内部计算，避开了处理“自由边界”的困难。课题经如此分解，即可利用商品程序单独解决，两个层次之间的衔接亦可利用软件提供“自设功能(UDF)”设法解决。最后形成的解题“包”，分层和衔接对最终用户来说完全是“透明”的。本文采用的软件是Fluent/6.1，所设计的算法在海洋以及其他水工设计中有较为光广泛的实用截止价值。 关键词：自由面，Boussinesq方程，N-S方程，水汽比拟方程 [引言] 三维粘性流动研究的典型课题是围绕平板加圆柱开展的，它有相当普遍的代表性，在飞行物附近有它，在海洋工程中有它，但对后者来说尚有自由面的作用要考虑。本文的目的在于针对后者提出一种可行的数值模拟计算方法。 1．物理数学模型 1.1计算区域 设有一长为120m，宽为60m,平均高度为15m的平直河道。无穷远处的来流速度为2m/s。在河道的中央，有一直立椭圆桥墩，桥墩长轴为6.0m，短轴为2.0m。其水平截面和垂直截面分别如图（Fig1）、（Fig2）所示。由于流场长宽尺寸相对于桥墩长短轴尺寸分别为20倍和30倍，所以计算区域可以近似的看成无穷大流场。 水平截面（Fig1） 垂直截面（Fig2） 1.2 水气比拟方程 其中，对于一般情况（如图Fig3） （1） (2) 这里的和h分别表示上下物质面函数。其中，物质面不一定是水气交界面或者是壁面。 对上面两个方程求物质导数得， （3） （4） （Fig3） （5） 沿着z方向从-h到积分可得， （6） 将方程（3，4）带入上式，并且利用变限积分得， （7） （8） 垂直向上动力学方程， （9） 假设一：w方向流体质点产生的加速度很小，则可得， （10） 即， （11） 水平方向上的动力学方程为， （12） （13） 假设二：u、v和z无关，并且由方程（11）可得， （14） （15） 方程（12，13）化为， （16）