创建另一个平面x-0,不选择From Surface列表中的symmetry,在New Surface Name文本框中输入x-0
创建XY图
Display/Plots/XY Plot 单击Set Up按钮,弹出Solution XY Plot对话框,从Surfaces列表中选择centre-line,设置Plot Direction为X=0,Y=0和Z=1,从Y Axis Function下拉列表中选择Velocity和Z Velocity,单击Axes按钮,弹出Axes-Solution XY Plot对话框,勾选Options选项组中的Major Rules复选框和Minor Rules复选框,同理对Y轴也做同样操作,最后单击Plot。
6 / 15
如图所示的速度图可以用来计算管道入口段的长度,速度达到99.9%时的距离即为入口段长度。由以上计算得到的图,入口段长度为15.19m(此处计算结果并不是这样,原因有待于进一步研究),这与本书参考文献的报告值相符。
计算出口的摩擦系数
摩擦系数被用以表示管道的压降,其定义为:
f?其中?p为沿管道长度的压降,
2?PD 2L?vD为管道直径,
v为管道横截面的平均速度,
?为流体密度
压力平衡后,剪切力使该式的形式变为:
f?8? ?v2式中,?为壁面的剪切应力。 创建壁面和出口交界面的线。 Surface/Iso-Surface
从From Surface列表中选择wall,从Surface of Constant下拉列表中选择Mesh和Z-Coordinate,在Iso-Values文本框中输入20,在New Surface Name文本框中输入wall-outlet。
报告面积分
Report/Result Reports/Surface Integrals
计算出口的平均剪切应力,从Report Type下拉列表中选择Area-Weighted Average,从Field Variable下拉列表中选择Wall Fluxes和Wall Shear Stress,从Surfaces列表中选择wall-outlet,单击Compute。
7 / 15
Area Weighted Average(pascal)中显示的剪切力值为0. 计算出口的平均速度
从Report Type下拉列表中选择Mass-Weighted Average,从Field Variable下拉列表中选择Velocity和Z Velocity,从Surfaces列表中选择outlet,单击Compute。
8 / 15
Mass-Weighted Average(m/s)中显示平均速度值。摩擦应力以下式计算
f?8? ?v2② 显示壁面Y+的值。
Display/Graphics and Animation/Contours,设置如图
9 / 15
结果如图
大多数计算域的5(这里似乎并不是这样),除了进口附近节点以外。这表明增强型壁面处理方式作为壁面函数是可以接受的。
Y+值大于
10 / 15