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

1、Flow3D 学习3 算例 1 Aerospace TutorialAerospace Tutorial新建一个项目,Model Setup Tab-Meshing & Geometry Tab-Subcomponent Tab-Geometry Files-c:Flow3Dguistl_libtank.stl, Type and Potential使用缺省选项，因为将引入其它形状作为固体，Subcomponent 1中坐标范围（Min/Max）为：X: 5.0 15.0, Y: 5.015.0, Z: 0.015.0tank.stl的单位对FLOW-3DB说是未知的，可能是英寸、英尺

2、、毫米等，现在 假设模型是SI （国际单位），那么流体或固体的届性都应该是SI的。（这里有些糊涂，FLOW-3氐使用STL文件中的单位么？）模拟的情况为从圆柱形底部入口向球形水箱内充水，计算域应该和此形状范围相 近，略大一点但不能紧贴着形状边界。底边界的位置和边界条件类型有关，如果入口处流速已知那么模拟多少入口长度 没有关系，因为断面形状是固定的，但是如果特定位置的压力是已知的，那么要把边界放在该位置处因为压力会受入口长度的重力和粘性效应影响而变化。建议计算域要大丁最大几何尺寸的5%底边界除外，可以小丁 5%这样计算域底部和入口交义，不会挡住水流，因此计算域定义为X: 4.95 15.05 Y

3、: 4.95 15.05 Z: 0.05 15.05在Mesh-Cartesian的Block 1中按上面参数修改计算域尺寸，然后在 Block 1 上右键选择Update Mesh更新显示。Re = Reynold 数=Inertial Force/Viscous Force = UL/ vBo = Bond 数=Gravitational Force/Surface Tension Force = g p LA2/ oWe = Weber数=Inertial Force/Surface Tension Force = L#2p / gU是特征流速，L是特征长度，g是重力加速度，p是密度，b

4、是表面张力系数。 这个问题中大约用100s充满水，冲水体积540立米，入口直径2m入口流速为540/(100* 攵 *1A2)=1.7m/s , Reynolds 数为Re=1.7*2/(1.0e-6)=3.4e6由于Reynolds数远大于1,因此惯性力比粘性力更重要，即不需要很密的网格 考虑粘性力，也不需要指定粘性特性。考虑表面张力影响的系数为Bo = g p L2/ b = 9.8m/sA2 * 1000 kg/mA3 * (2m)A2/(0.073kg/sA2) = 5.4e5 We = LU2p / b = 2m * (1.7 m/s)A2 * 1000 kg/mA3 / (0.07

5、3kg/sA2) = 7.9e4 可以看出表面张力也不需要考虑。缺省网格在X、Z方向为10, Y方向为1。在Block 1上右键选择Automatic Mesh ,通过网格大小或总网格数修改网格。 由于入口直径2m最少应该有5个网格，因此网格大小为0.4m,选择所有方向 均为0.4生成网格。当前水箱为固体，在 Geometry的Component 1下面展开Subcomponent 1,把 SubcomponentType由Solid变为Complement什么意思？),网格定义变为正 常情况。选择工具栏中的Favor图标，在FAVORize下选择using computational mes

6、h , 在 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。低速流入时不必考虑气体影响，因此在 Nu

7、mber of fluids 中仅选 择一种流体。Interface Tracking 选择 Free surface or sharp interface , 如 果是气体注入则用No sharp interface 选项并用2种流体让气体扩散。Flow Mode 选择不可压流体。Initial tab :这个模拟中不需要设置初始条件。Physics tab :设置使用的模块，因为是高 Reynolds数，不需要考虑粘性，所以 仅使用 Gravity , Z-direction 中输入-9.8 。Viscosity 和 Surface tension 模 块不必考虑。Fluids tab :

8、仅需要考虑流体密度，在左侧 Density中的Fluid 1 的Density 中输入1000 (或者找到 Water at 20 C si用Load Fluid 1 载入)Output tab : Flow-3D始终输出速度、温度、流体fraction 等，缺省输出间隔 为1/10结束时间，暂时取缺省值。Simulate tab : Run Simulation 。然后可以观看结果。下面进行热启动计算，关闭入口边界，定义x方向的加速度使流体摇晃，在计算 域中嵌套一个块体并计算流体摇晃时流体对块体的作用力。转到Navigator下，File>Add Restart Simulation

9、,创建热启动模拟，模拟名 称取缺省值，选择OK在General tab中选择Restart按钮，数据源(flsgrf.* )应该正确，热启动开 始时间为100s,选择reset time to zero 保证非惯性输入从0开始，其它取缺 省选项。修改结束时间为12s。热启动既可以从输入文件中读入新的数据也可以从上次计算结果中读取数据，这些都在刚才的对话框(Restart Options )中设置。热启动时间输出在 hd3msg中，可以在Restart Options 的Times下拉歹U表中选择，还可以在 Diagnostics>Solver Messages 中查看。在Boundari

10、es tab 中设置Z-min边界为symmetry边界条件。在Physics Tab中选择Non-Inertial Reference Frame按钮，激活非惯性作用， 并使用Tabular Input ,输入数据0 5 0.0 9.8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.03-5 0.0 9.8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.06 5 0.0 9.8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.09 -5 0.0 9.8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.010 0 0.0 9.8 0.0 0.

11、0 0.0 0.0 0.0 0.0 0.0 0.012 0 0.0 9.8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0在 Meshing and Geometry 中增力口 Component 2, 选择 Subcomponent-Cylinder 并选择NewComponent2，圆柱半径设置为 0.25 , Z-low为7, Z-high为12,单 击transform 按钮在X和Y方向均平移10。设置完后Geometry歹U表中多出 Component 2项。(每个 Component表示具有不同性质的材料，同一 Component 下的每个Subcomponent都

12、是相同材料的)要计算障碍物上的作用力，创建嵌套的加密网格包含障碍物，这样计算结果更准确，嵌套计算网格范围为：X: 9.0 to 11.0 Y: 9.0 to 11.0 Z: 6.5 to 12.5在Mesh-Cartesian上右键选择 Add Mesh Block ,输入上面的计算范围。嵌套网格要保持比例在2:1以内，因此在Block 2上右键选择Auto mesh,设置所有方 向网格大小为0.2。为保证插值准确，在 Block 1中增加Block 2的边界点，在 Mesh Cartesian的 Block 1的X Direction 上右键 Add,增加9和11两个点，在 Y方向相同，Z方

13、 向增加6.5和12.5。(这是只增加点，不必重新划分网格？)在 Numerics tab 下选择 Pressure solver options 下的 Convergence controls , 缺省 inter-block boundary type coefficient 为 0.25，修改为 0。(如果不为 0是什么含义？)为计算障碍物上的作用力，包括压力和粘性剪切应力(不是没有考虑粘性？)，需要手动编辑输入文件，File-Edit simulation 在$grafic namelist下增加障碍物上的力(这是什么意思？可能需要参考Input变量的说明) nwinf=1,xf1(1

14、)=9.4, xf2(1)=10.6, yf1(1)=9.4, yf2(1)=10.6, zf1(1)=6.9,zf2(1)=12.1,修改 Model Setup-Output 中的 Restart Data-time interval 为 0.24 ,这样每 计算0.24s输出一次结果，12s会有50次输出。Simulate-Run simulation开始模拟在Analyze tab下选择space_restart打开结果文件，作用力输出在Probe Tab下，在Data source中选择General History ,可以检查障碍物上 各方向的作用力，注意这些

15、作用力包括压力和粘性剪应力，在列表中选择 x-force window 1, Render查看。Flow3D学习4 算例 2 Casting Tutorial创建新的 workspace，名称为 Casting Example，仓U建新的 simulation 名称为 casting。导入几个stl文件，每个文件都是Component 1的Subcomponent即材料特性将应用丁里面的每I注入口（sprue）先作为固体导入以便观看，然后修改为complement:,这样网格区域内除几何形状在的区域都变为固体，即对那部分进行取反操作。另外 cad生成的stl文件是毫米单位，这里用 单位，需要缩

16、放1/10。导入 c:Flow3Dguistl_libsprue2.stl 和 part.stl ，都作为 Component 1 的 Subcomponent,: part.stl 的类型改为 Hole。把每个 Subcomponent的 Global Magnifications 变为 0.1。Subcompoi 2需要平移，X方向Translations 为3, Z方向0.5。每次输入后按Enter键可以更新显示。增加一个 Box作为 Hole,参数为 X 4.5-8.5, Y -0.4-0.4, Z -0.4-0.6使用3块网格计算，修改计算网格参数# Cells Low HighX:

17、 6 -1.5 1.5Y: 1 -0.1 0.1Z: 30 -3.4 11.6X: 21 1.5 12.0Y: 1 -0.1 0.1Z: 6 -2.4 0.4X: 16 2.5 10.5Y: 1 -0.1 0.1Z: 20 0.4 10.4此时可以把Subcomponent 1从Solid变为complement:,这样会填充包围它的网格。记住，几何形状和网格块无关，几何形状会自动占据网格块内的空间。为检查网格，进行前处理 Simulate>Preprocess Simulation，然后在Analyze tab 下选择prpgrf.casting ，在2-D Tab中的右下角Mesh

18、Block按钮弹出的对话框中选择所有块，Render示可以用于检查的图形。（这里给出了检查的方式）Block 1顶部注入口边界流速是变化的，不易指定，可以使用定常的 stagnation pressure表示上 注入溶液高度固定（大溶液池）或者 stagnation pressure 随时间变化，即上面液面高度逐渐降低 最简单的是线性变化。先假设压力为常数，在边界中使用压力18375.0并勾选Stagnation pressur 注入温度需要指定，在 Physics Tab中选择Heat Transfer , 1st Order , OK然后重新在刚才边 条件对话框中单击 Thermal In

19、formation 按钮，输入900.0 , OKFlow-3D会自动设置网格块连接的边界条件为Mesh Block ,也可以自己设置。（哪里设置？）在General Tab中设置结束时间为0.5s , Flow-3D始终在结束时间到达后停止计算，也可以在Finicondition中指定填充比例作为结束条件，如果不设置Fill fraction 而在入口处使用流速边界条时在注满溶液后求解器可能因为要继续注入溶液而出错退出，这里设置为1.0即全部填充并把Ur改为CG卯位。因为空间内开始时没有液体，因此这个例子中也就不存在初始条件。Physics Tab中需要使用的模块包括：模块名称理由Visco

20、sity湍流引起的墙面剪应力 Wall shear force and turbulent mixing ,选择牛顿粘性Renormalized Group （RNG兢型，此模型健壮性、易用性和精度均较好。（紊流模型都推荐这个吗 其它取缺省值Energy equation（流体中热传导）计算金届中的温度 To compute the temperature in the m（Heat Transfer and Conduction in obstacles砂模和其中流体的热交换 Heat exchange betw（the sand mold and the fluid and within

21、mold能量对流有2个选项，二阶计算更精确但是更费 通常在浇铸中使用一阶即可。流体和固体热传导中的无热传导选项用于隔热模拟，后面的2个选顼 示温度从流体传递给温度恒定的固体，因为浇铸时间很长，砂模温度变化很大，因此这里需要使用 Full Energy Equation选项求解完全的能量方程。Gravity重力加速度 Acceleration due to gravity CGS 单位中 Z 向加速度是-980。Solidification模拟注入过程中的金届凝固 Simulate metal freezing during the filling能没有充满就凝固了，这里激活此部分并取默认值（模

22、拟没有shrinkage的情况）。Defect Tracking记录氧化物相对浓度和缺陷的可能位置Records the relative concentratioxides and thus the probable location of oxide related defects此模块计算液态金届暴露空气中的时间,能表示溶液前锋氧化物含量水平,在氧化物缺陷位置预测时有用。由于氧化物生成度没有试验确定，因此这里用1表示氧化物相对浓度。在Fluids tab 中载入alloy AL 356 in SI届性并在对话框中选择转化为 CG驴位。推荐在设置模后立刻设置材料届性，因为有些模块如Soli

23、dification和特定的流体特性有关，这样可以保证需的特性从数据库中载入了。（材料特性是否必须立刻载入，能否设置完所有模块后再载入？）下面在Meshing & Geometry Tab中的Tools菜单中选择 Solids Database，从歹0表中选择 Gener Silica Sand at 1000K ，载入，此时 thermal conductivity 和 density*specific Heat 等都应该 数值了，然后在左侧 Component 1的Initial Conditions中修改温度为293。运行，然后载入flsgrf.casting 看结果。先在Ana

24、lyze Tab右下角的MeshBlock中选择所有网格 否则只会显示一块。Contour变量用Temperature , Render后可以看到结果。由于速度很快，温度 于熔点，因此溶液没有凝固现象。如果结果中有空洞可能由两个原因，一是Fill Fraction Finish T 条件只检查空网格，而不是部分空网格引起，因此可能在没有空网格但是有部分空网格时计算结束 二是温度图把流体体积比小于 0.5 （液体比例小于50%的网格屏蔽了，这个比例可以在Analyze 左下角的Advanced对话框中设置。把等值线变量修改为Surface Defect Concentration , Rend（

25、后可以看出氧化物在金届表面累积情况，可以看出顶部较高，可以设置一个冒口把这部分溶液排出把计算网格Y方向尺度变为-2.02.0，网格大小为0.5，把模型变为三维情况。计算时间修改为1.0 边界条件变为随时间变化的情况，线性变化到 0.7s时2450,此后压力值保持不变。运行后观看， Iso-Surface 变量设置为fraction of fluid，表示溶液表面，Color Variable用选择的变量覆盖体的iso-surface ，选择Temperature , Render。(对这些设置糊涂，哪里有说明？)创建热启动模拟凝固,在 General tab 中选择从0.81s开始，并选中re

26、start time to zero保讯始的非惯性力从0开始(上个例子也是这样，这是什么意思？)，其它用缺省选项，结束时间为12 凝固模拟时假设上部没有流体注入， 在Boundaries Tab中Block 1的Z-max边界设置Fluid Fracti 为 0 (应该是 F Fraction )。在 General tab 中修改 Finish Condition 为 Solidified Fluid Fract 允许FLOW-3训模拟结束时间前停止。模拟120s要很长时间，凝固与热传导有很大关系，在 Numerics Tab中缺省是Explicit 方式求解 需要很小的时间步，用Impli

27、cit方式可以避免这个问题(Explicit/implicit solver options叶Heat transfer 选项)。假设在凝固过程中 z向流速很小，在 Fluid flow solver options中选择I Zero Velocity Field可以大大减少计算时间。并在输出选项中选择Solid fraction 。开始模拟。计算完成后在Analyze Tab中载入flsgrf.dat ，选择2-D Tab,选择3个网格块和XZ平面，在D Source中选择Selected，并用Solid Fraction做等值线变量，可以看出凝固过程是从右侧开始的然后向上凝固。还有其它变量

28、可以检查浇铸过程，在Data Source中选择Solidification ,可用等值线变量包括Solidification Time网格中温度降低到凝固温度下Solidus Velocity(VEL) 表示solidus front 运动的速度,-1表示没有数据，如果两个solidus fro 碰撞时FLOW-3DC法确定速度或者开始凝固时发生。Cooling Rate(SCR)在凝固温度时单位时间内能量的变化Local Solidification Time(LST)网格在液态和固态之间变化使用的时间还有一些其它变量，因不会用到，略过。Flow3D学习5 算例 3 Hydraulics

29、TutorialHydraulics Tutorial这个问题中流体从18cm高的坝上流下，根据自由落体运动，水流流速为Velocity = sqrt(2*980*18) = 187.8 cm/sReynold 数为 Re = 30cm x 187.8cm/sA2 / 10-2cmA2/s = 5.6 x 10A5因此粘性应力可以忽略，不需要用很密的网格描述墙面粘性剪切应力，但是由丁 存在紊流，流体内会存在粘性剪切应力，因此需要在模拟中使用粘性模型。Bond数和Weber数为Bo = 980cm/sA2 * 1 gm/cc * (30cm)A2/(73gm/sA2) = 1.2 x 10八4W

30、e = 30cm * (187.8 cm/s)-2 * 1gm/cc / (73gm/s2) = 1.45 x 10八4即不需要考虑表面张力。根据问题的性质，采用对称方式模拟。创建workspace和simulation ，在General中设置结束时间为1.0s，单位CGS 在下面Notes中的第一行是标题，将出现在所有的输出文件和图形中， 这里写为“Flow over a Weir ” 。在 Meshing&Geometry tab 中导入 c:flow3dguistl_libweir1.stl,按缺省设置。下面增加上游水库的底面，单击工具栏Box图标，为定义不同的届性(如 糙率参

31、数),选择 New Component X -10 1, Y -20 20, Z 0.0 8。选取计算域时上游不能太小，否则会因为突然出现的加速度引起计算不稳定，下 游如果计算范围小会使边界条件影响流体流态，但是计算范围太大会使计算时间 增加。这个例子先用一块网格计算，然后嵌套加密的网格，网格定义可以用手工或图形 方式。设置计算网格范围X -1020, Y 010, Z 018，各方向网格分别为30、 10、18,为在对称处加密，在 Y方向的Pt(1)选项中设置网格大小为0.5。Boundaries Tab :这个问题中X两端边界被设置为 hydrostatic pressure 边界 条件，

32、如果边界的“fluid height ”比计算域内高则会向计算域内流入，否则会 流出。在X Min边界上选择Specified Pressure 边界条件同时选择 Stagnation Pressure，并设置 Fluid Height 为 15.5 ;在X Max中同样使用 Specified Pressure 和择 Stagnation Pressure ，并设置 Fluid Height 为 1.7 ,设置 F fraction 为0.0 ,这样可以防止流体从这个边界进入计算域，但是允许流体自由流出。（开 边界？）Initial tab :由于在X两端使用了静水压力边界，必须激活静水压力

33、选项才能 保证边界条件正确设置，在Initial pressure field 中选择Hydrostatic pressure in Z-direction ,这将把网格中所有流体初始化为静水压力，同时指 定垂向压力边界条件为静水压力边界条件。然后创建网格中的初始流体，单击 Add fluid ，在对话框右下方的 Fluid options 中选择Add fluid , X High输入 0, Z High输入15 （其它不必输入？默认为计算网格极值？）（上游，边界条件 水位15.5），单击OK然后再次增加一个流体区域 X Low为1, Z High为1.7 （下游，边界条件水位1.7），这2

34、个流体区域在网格中增加0时刻的初始流体。Physics tab :选择 Viscosity and turbulence 模块，使用 Newtonian viscosity , 在 Wall shear boundary conditions 中使用 No-slip or partial slip 条件（墙 面剪切应力模型，不允许流体沿墙面滑移），单击 OK在Gravity中Z向为 -980（cm/s2）。Fluids tab :需要流体的密度和粘性，载入 water at 20 C。如果不在选择模型 后立刻设置流体届性，在关闭simulation时没有指定届性的模块会关闭。Output ta

35、b :在 restart data 下的 Time interval 中输入 0.05，即每隔 0.05s 输出一次。Numerics tab :不需要修改。Simulate->Preprocess Simulationion，前处理，在 Analyze tab 中打开prpgrf.dat ，在 3D tab 中的 Component Iso-surface Overlay 中选择 Solid volume, Render,检查。模拟后打开flsgrf.dat 观看。增加初始条件：下面在出口前端增加一块流体并赋予 x方向的初始速度好让流态 尽快稳定。Add simulation copy

36、，在初始条件中增加一块流体，X 0.02.5 , Y-5.0 5.0 , Z 13.0 15.0 , OK 在 Fluid initial state 中设置 U-velocity 为20。检查，模拟，观看。增加嵌套网格：复制 simulation ，在Meshing&Geometry tab中选择工具栏上的 Create Mesh（先要使用2D投影方式），用鼠标创建嵌套网格区域，然后在左侧 修改网格范围-3<x<3, 0<y<5, and 8<z<18 ，网格数分别为12、10、20。前处理 后可以在Diagnostics->Preproce

37、ssor Report中找到每块网格的边界条件。excessive convection failures - calculation terminated双流体模拟：复制simulation , General tab 中选择Two Fluids ,尽管空气是 可压缩的，但是由于 Mach数很低，可以看作是不可压缩流体。Fluids tab中载 入fluid 2的届性为Air at 15 degrees ，单击Load Fluid 2 载入，展开左侧参数检查。运行，观察。（用上次的结果进行双流体模拟失败？）需要指出，计算结果并不是物理上真实的，主要问题在于每个网格内只计算了一 个速度。当两种流体密度差别较大，比如大于 10倍时，因为同一网格中不同流 体的流速差别较大，即同一网格内空气和水的速度有较大差别，这种混合速度计 算方法就不精确，后果就是两种流体间的滑移条件没有精确模拟，可能导致流体界面运动模拟不正确。当两种流体的密度差别到一定程度时，双流体模型就不精确（但是对很多CFD软件来说，这种双流体模型是计算自由表面的唯一方法）。为精确模拟自由表面， FLOW-3映用了 true VOF方法，这