[2018年最新整理]flow3dHydraulics教程.docVIP

Flow3d 9.3.2 Hydraulics Tutorial水力教程 本练习的目的是模拟水从水库通过薄壁堰流进下游水池。 图1 水流模拟Re (Reynold’s Number) 雷诺数 =Inertial Force/Viscous Force =贯性力/黏性 =UL/ν Bo (Bond Number) 邦德数 =Gravitational Force/Surface Tension Force =重力/表面张力 =gΔρL2/σ We (Weber Number) 韦伯数 =Inertial Force/Surface Tension Force =惯性力/表面张力 =LU2ρ/σ 这里U是特征速度，L是特征长度，g是重力加速度，ρ是密度，σ 是表面张力系数。 对本问题，水从18cm高堰流过，水流在堰底的速度可近似按自由落体运动分析得出： Velocity = sqrt(2*980*18) = 187.8 cm/s 流体的雷诺数为： Re = 30cm x 187.8cm/s / 10-2cm^2/s = 5.6 x 105 雷诺数大，意味着与贯性力相比，黏性力不可忽略。因此，我们不需要精细的网格求解壁黏性剪切层。当然，由于流态的紊乱，液体内部有很多黏性剪切力，因此，需要在模型中指定黏性参数。 邦德数按下式求得： Bo = 980cm/s^2 * 1 gm/cc * (30cm)^2/(73gm/s^2) = 1.2 x 104 韦伯数按下式求得： We = 30cm * (187.8 cm/s)^2 * 1gm/cc / (73gm/s^2) = 1.45 x 104 再者，大的邦德数和大的韦伯数表明，与重力和惯性力相比，表面张力可忽略。模型是这种情况时，不考虑表面张力。 问题的大小（模型运行的时间）可以利用堰中心顺水流平面的对称特性进行简化。因此，我们仅仅需要模拟整个范围的一部分（即堰的后半部分），就可也得到堰的全部信息。我们已经对问题进行了简化，下面是如何建立这些条件，如何确定几何条件，利用flow3d求解问题。 建模 总体参数 点击“Model Setup”表的“General”表，“General”是确定整个问题的参数，如结束时间、结束条件、界面追踪Simulation units”标题菜单中，选CGS单位（厘米。克。秒），其它设置采用缺省设置。 在“总信息表Global tab”的底部注释中，你可以在第一行为问题指定一个名字。名字会出现在所有输出文件和图形上。本例名称为“Flow over a Weir”（过堰流体）。 建立几何体Geometry Setup 我们将添加元件定义堰体。首先，我们输入一个已有的STL文件，weri1.stl，该文件放在目录“c:\flow3d\gui\stl_lib”。切换到“几何与分网Meshing Geometry ”表，单击工具条STL图标，会打开标题为“几何Geometry”的对话框，点击添加，打开对话框，找到并选择weri1.stl。在“Geometry File”点击ok，接受缺省设置。在之后出现的添加部件对话框中接受缺省设置。现在STL文件已经输入，并且出现在工作空间中。输入文件也被列在树形结构表中。 下面，我们将通过“FLOW-3D”简单建模创建另一个组件，来添加上游水库河床。在工具栏点击box（盒子）图标，盒子对话框显示如图2。 图2 盒子组件选项FLOW-3D”组件中的所有子部件共享相同的参数。因此，在盒子子部件对话框的下拉菜单中，选择“New Component 2”，输入盒子尺寸，如图2。其它按缺省设置。 在盒子对话框中，单击ok，接受缺省参数。下一个对话框出现（添加组件）。盒子会作为新子部件在左侧树形结构及工作空间中出现。 注：对本教程，我们提供了河床的范围。为了你能够自己确定河床尺寸，你可以通过树形结构查看STL文件的尺寸。这时，通过单击“+”号，打开部件1数据，打开子部件1分支，查看部件的x、y、z方向的最大、最小尺寸。 分网 在进行任何模拟时，最重要的工作之一是考虑如何定义计算网格。网格单元的数量，取决于定义边界的尺寸。并且，网格单元的数量极大地影响计算结果、运行时间、计算精度。因此，问题的范围必须仔细选择。总之，计算人在问题未做好之前，应当谨慎勾画问题的简图。 对本问题，我们需要定义的问题是2个：堰后流体及堰前面流体流过的范围，当然还有堰本身。注意，不要将范围定得太小，如图3、图4。如果堰上游范围太小，见图3，计算结果可能不稳定，因为可能会出现突然的加速度。如果堰下游定得太小，见图4，边界条件会影响流态。同样，边界范围不应定得太大，因为没有必要增加计算规模（计算时间）。 图3 计算结果不稳定