flow3d官方培训教程中的实例中文说明.docxVIP

Flow3D 学习——3 算例 1 Aerospace Tutorial Aerospace Tutorial 新建一个项目，Model Setup Tab-Meshing Geometry Tab-Subcomponent Tab-Geometry Files-c:\Flow3D\gui\stl_lib\tank.stl，Type and Potential 使用缺省选项，因为将引入其它形状作为固体，Subcomponent 1 中坐标范围 (Min/Max)为： X: 5.0～15.0, Y: 5.0～15.0, Z: 0.0～15.0 tank.stl 的单位对 FLOW-3D 来说是未知的，可能是英寸、英尺、毫米等，现在 假设模型是 SI（国际单位），那么流体或固体的属性都应该是 SI 的。（这里有 些糊涂，FLOW-3D 会使用 STL 文件中的单位么？） 模拟的情况为从圆柱形底部入口向球形水箱内充水，计算域应该和此形状范围 相近，略大一点但不能紧贴着形状边界。 底边界的位置和边界条件类型有关，如果入口处流速已知那么模拟多少入口长 度没有关系，因为断面形状是固定的，但是如果特定位置的压力是已知的，那 么要把边界放在该位置处因为压力会受入口长度的重力和粘性效应影响而变化。 建议计算域要大于最大几何尺寸的 5%，底边界除外，可以小于 5%，这样计算域 底部和入口交叉，不会挡住水流，因此计算域定义为 X: 4.95～15.05 Y: 4.95～15.05 Z: 0.05～15.05 在 Mesh-Cartesian 的 Block 1 中按上面参数修改计算域尺寸，然后在 Block 1 上右键选择 Update Mesh 更新显示。 Re = Reynold 数 = Inertial Force/Viscous Force = UL/ν Bo = Bond 数 = Gravitational Force/Surface Tension Force = gΔρL^2/σ We = Weber 数 = Inertial Force/Surface Tension Force = LU^2ρ/σ U 是特征流速，L 是特征长度，g 是重力加速度，ρ 是密度，σ 是表面张力系数。 这个问题中大约用 100s 充满水，冲水体积 540 立米，入口直径 2m，入口流速为 540/(100*π*1^2)=1.7m/s，Reynolds 数为 Re=1.7*2/(1.0e-6)=3.4e6 由于 Reynolds 数远大于 1，因此惯性力比粘性力更重要，即不需要很密的网格 考虑粘性力，也不需要指定粘性特性。考虑表面张力影响的系数为 Bo = gΔρL2/σ = 9.8m/s^2 * 1000 kg/m^3 * (2m)^2/(0.073kg/s^2) = 5.4e5 We = LU2ρ/σ = 2m * (1.7 m/s)^2 * 1000 kg/m^3 / (0.073kg/s^2) = 7.9e4 可以看出表面张力也不需要考虑。 缺省网格在 X、Z 方向为 10，Y 方向为 1。 在 Block 1 上右键选择 Automatic Mesh，通过网格大小或总网格数修改网格。 由于入口直径 2m，最少应该有 5 个网格，因此网格大小为 0.4m，选择所有方向 均为 0.4 生成网格。 当前水箱为固体，在 Geometry 的 Component 1 下面展开 Subcomponent 1，把 Subcomponent Type 由 Solid 变为 Complement（什么意思？），网格定义变为 正常情况。 选择工具栏中的 Favor 图标，在 FAVORize 下选择 using computational mesh， 在 view component volume 中选择 Open，单击 Render，保存项目显示。用 cutting planes 或 mesh slice 可以看到网格内部。（这里是什么作用？为什么 第二个例子里没有这步？） 选择 Model Setup 中的 Boundaries tab 定义边界条件，选择 Z-min 边界，选择 Specified Velocity 按钮然后输入 W-velocity 为 1.7，其它边界条件不变。 General tab：设置计算结束时间为 100s，并在 Finish condition 中选择使用 Finish time。低速流入时不必考虑气体影响，因此在 Number of fluids 中仅 选择一种流体。Interface Tracking 选择 Free surface or sharp interface，如果是气体注入则用 No sharp interf