flow3d官方培训教程中的实例中文说明文档格式.docx

1、Bo = Bond数 = Gravitational Force/Surface Tension Force = gL2/We = Weber数 = Inertial Force/Surface Tension Force = LU2/U是特征流速，L是特征长度，g是重力加速度，是密度，是表面张力系数。这个问题中大约用100s充满水，冲水体积540立米，入口直径2m，入口流速为540/（100*12）=1.7m/s，Reynolds数为Re=1.7*2/（1.0e-6）=3.4e6由于Reynolds数远大于1，因此惯性力比粘性力更重要，即不需要很密的网格考虑粘性力，也不需要指定粘性特性。考虑

2、表面张力影响的系数为Bo = gL2/ = 9.8m/s2 * 1000 kg/m3 * （2m）2/（0.073kg/s2） = 5.4e5We = LU2/ = 2m * （1.7 m/s）2 * 1000 kg/m3 / （0.073kg/s2） = 7.9e4可以看出表面张力也不需要考虑。缺省网格在X、Z方向为10，Y方向为1。在Block 1上右键选择Automatic Mesh，通过网格大小或总网格数修改网格。由于入口直径2m，最少应该有5个网格，因此网格大小为0.4m，选择所有方向均为0.4生成网格。当前水箱为固体，在Geometry的Component 1下面展开Subcomp

3、onent 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-velocit

4、y为1.7，其它边界条件不变。General tab：设置计算结束时间为100s，并在Finish condition中选择使用Finish time。低速流入时不必考虑气体影响，因此在Number of fluids中仅选择一种流体。Interface Tracking选择Free surface or sharp interface，如果是气体注入则用No sharp interface选项并用2种流体让气体扩散。Flow Mode选择不可压流体。Initial tab：这个模拟中不需要设置初始条件。Physics tab：设置使用的模块，因为是高Reynolds数，不需要考虑粘性，所以仅

5、使用Gravity，Z-direction中输入-9.8。Viscosity和Surface tension模块不必考虑。Fluids tab：仅需要考虑流体密度，在左侧Density中的Fluid 1的Density中输入1000（或者找到Water at 20 Csi用Load Fluid 1载入）Output tab：Flow-3D始终输出速度、温度、流体fraction等，缺省输出间隔为1/10结束时间，暂时取缺省值。Simulate tab：Run Simulation。然后可以观看结果。下面进行热启动计算，关闭入口边界，定义x方向的加速度使流体摇晃，在计算域中嵌套一个块体并计算流体

6、摇晃时流体对块体的作用力。转到Navigator下， Restart Simulation，创建热启动模拟，模拟名称取缺省值，选择OK。在General tab中选择Restart按钮，数据源（flsgrf.*）应该正确，热启动开始时间为100s，选择reset time to zero保证非惯性输入从0开始，其它取缺省选项。修改结束时间为12s。热启动既可以从输入文件中读入新的数据也可以从上次计算结果中读取数据，这些都在刚才的对话框（Restart Options）中设置。热启动时间输出在hd3msg中，可以在Restart Options的Times下拉列表中选择，还可以在Diagnost

7、icsSolver Messages中查看。在Boundaries 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

8、.0 0.0 0.010 0 0.0 9.8 0.0 0.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并选择New Component 2，圆柱半径设置为0.25，Z-low为7，Z-high为12，单击transform按钮在X和Y方向均平移10。设置完后Geometry列表中多出Component 2项。（每个Component表示具有不同性质的材料，同一Component下的每个

9、Subcomponent都是相同材料的）要计算障碍物上的作用力，创建嵌套的加密网格包含障碍物，这样计算结果更准确，嵌套计算网格范围为： 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方向相同，

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

11、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下选择flsgrf.aerospace_restart打开结果文件，作用力输出在Probe Tab下，在Data source中选择General History，可以检查障碍物上各方向的作用力，注意这些作用力包括压力和粘性剪应力，在列表中选择x-force wind

12、ow 1，Render查看。Flow3D学习4 算例2 Casting Tutorial创建新的workspace，名称为Casting Example，创建新的simulation名称为casting。导入几个stl文件，每个文件都是Component 1的Subcomponent，即材料特性将应用于里面的每个固体。注入口（sprue）先作为固体导入以便观看，然后修改为complement，这样网格区域内除几何形状存在的区域都变为固体，即对那部分进行取反操作。另外cad生成的stl文件是毫米单位，这里用CGS单位，需要缩放1/10。导入c:Flow3Dguistl_libsprue2.stl

13、和part.stl，都作为Component 1的Subcomponent，把part.stl的类型改为Hole。把每个Subcomponent的Global Magnifications变为0.1。Subcomponent 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 High 6 -1.5 1.5Y: 1 -0.1 0.1Z: 30 -3.4 11.6 21 1.5 12.0 -

14、2.4 0.4 16 2.5 10.5 20 0.4 10.4此时可以把Subcomponent 1从Solid变为complement，这样会填充包围它的网格。记住，几何形状和网格块无关，几何形状会自动占据网格块内的空间。为检查网格，进行前处理SimulatePreprocess Simulation，然后在Analyze tab下选择prpgrf.casting，在2-D Tab中的右下角Mesh Block按钮弹出的对话框中选择所有块，Render后显示可以用于检查的图形。（这里给出了检查的方式）Block 1顶部注入口边界流速是变化的，不易指定，可以使用定常的stagnation pr

15、essure表示上面注入溶液高度固定（大溶液池）或者stagnation pressure随时间变化，即上面液面高度逐渐降低，最简单的是线性变化。先假设压力为常数，在边界中使用压力18375.0并勾选Stagnation pressure。注入温度需要指定，在Physics Tab中选择Heat Transfer，1st Order，OK，然后重新在刚才边界条件对话框中单击Thermal Information按钮，输入900.0，OK。Flow-3D会自动设置网格块连接的边界条件为Mesh Block，也可以自己设置。（哪里设置？在General Tab中设置结束时间为0.5s，Flow-3

16、D始终在结束时间到达后停止计算，也可以在Finish condition中指定填充比例作为结束条件，如果不设置Fill fraction而在入口处使用流速边界条件时在注满溶液后求解器可能因为要继续注入溶液而出错退出，这里设置为1.0即全部填充并把Units改为CGS单位。因为空间内开始时没有液体，因此这个例子中也就不存在初始条件。Physics Tab中需要使用的模块包括：模块名称 理由Viscosity 湍流引起的墙面剪应力Wall shear force and turbulent mixing，选择牛顿粘性的Renormalized Group （RNG）模型，此模型健壮性、易用性和精度

17、均较好。（紊流模型都推荐这个吗？）其它取缺省值Energy equation（流体中热传导） 计算金属中的温度 To compute the temperature in the metalHeat Transfer and Conduction in obstacles 砂模和其中流体的热交换Heat exchange between the sand mold and the fluid and within mold 能量对流有2个选项，二阶计算更精确但是更费时，通常在浇铸中使用一阶即可。流体和固体热传导中的无热传导选项用于隔热模拟，后面的2个选项表示温度从流体传递给温度恒定的固体，因为

18、浇铸时间很长，砂模温度变化很大，因此这里需要使用Full Energy Equation选项求解完全的能量方程。Gravity 重力加速度 Acceleration due to gravity CGS单位中Z向加速度是-980。Solidification 模拟注入过程中的金属凝固 Simulate metal freezing during the filling 可能没有充满就凝固了，这里激活此部分并取默认值（模拟没有shrinkage的情况）。Defect Tracking 记录氧化物相对浓度和缺陷的可能位置 Records the relative concentration oxi

19、des and thus the probable location of oxide related defects 此模块计算液态金属暴露在空气中的时间，能表示溶液前锋氧化物含量水平，在氧化物缺陷位置预测时有用。由于氧化物生成速度没有试验确定，因此这里用1表示氧化物相对浓度。在Fluids tab中载入alloy AL 356 in SI属性并在对话框中选择转化为CGS单位。推荐在设置模块后立刻设置材料属性，因为有些模块如Solidification和特定的流体特性有关，这样可以保证需要的特性从数据库中载入了。（材料特性是否必须立刻载入，能否设置完所有模块后再载入？下面在Meshing &

20、 Geometry Tab中的Tools菜单中选择Solids Database，从列表中选择Generic Silica Sand at 1000K，载入，此时thermal conductivity和density*specific Heat等都应该有数值了，然后在左侧Component 1的Initial Conditions中修改温度为293。运行，然后载入flsgrf.casting看结果。先在Analyze Tab右下角的Mesh Block中选择所有网格块，否则只会显示一块。Contour变量用Temperature，Render后可以看到结果。由于速度很快，温度高于熔点，因此溶

21、液没有凝固现象。如果结果中有空洞可能由两个原因，一是Fill Fraction Finish Time条件只检查空网格，而不是部分空网格引起，因此可能在没有空网格但是有部分空网格时计算结束；二是温度图把流体体积比小于0.5（液体比例小于50%）的网格屏蔽了，这个比例可以在Analyze Tab左下角的Advanced对话框中设置。把等值线变量修改为Surface Defect Concentration，Render后可以看出氧化物在金属表面累积情况，可以看出顶部较高，可以设置一个冒口把这部分溶液排出。把计算网格Y方向尺度变为-2.02.0，网格大小为0.5，把模型变为三维情况。计算时间修改为

22、1.0s，边界条件变为随时间变化的情况，线性变化到0.7s时2450，此后压力值保持不变。运行后观看，Iso-Surface变量设置为fraction of fluid，表示溶液表面，Color Variable用选择的变量覆盖流体的iso-surface，选择Temperature，Render。（对这些设置糊涂，哪里有说明？创建热启动模拟凝固，在General tab中选择从0.81s开始，并选中restart time to zero保证初始的非惯性力从0开始（上个例子也是这样，这是什么意思？），其它用缺省选项，结束时间为120s。凝固模拟时假设上部没有流体注入，在Boundaries

23、Tab中Block 1的Z-max边界设置Fluid Fraction为0（应该是F Fraction）。在General tab中修改Finish Condition为Solidified Fluid Fraction允许FLOW-3D在模拟结束时间前停止。模拟120s要很长时间，凝固与热传导有很大关系，在Numerics Tab中缺省是Explicit方式求解，需要很小的时间步，用Implicit方式可以避免这个问题（Explicit/implicit solver options中的Heat transfer选项）。假设在凝固过程中z向流速很小，在Fluid flow solver op

24、tions中选择Use Zero Velocity Field可以大大减少计算时间。并在输出选项中选择Solid fraction。开始模拟。计算完成后在Analyze Tab中载入flsgrf.dat，选择2-D Tab，选择3个网格块和XZ平面，在Data Source中选择Selected，并用Solid Fraction做等值线变量，可以看出凝固过程是从右侧开始的，然后向上凝固。还有其它变量可以检查浇铸过程，在Data Source中选择Solidification，可用的等值线变量包括Solidification Time 网格中温度降低到凝固温度下Solidus Velocity（

25、VEL） 表示solidus front运动的速度，-1表示没有数据，如果两个solidus front碰撞时FLOW-3D无法确定速度或者开始凝固时发生。Cooling Rate（SCR） 在凝固温度时单位时间内能量的变化Local Solidification Time（LST） 网格在液态和固态之间变化使用的时间还有一些其它变量，因不会用到，略过。Flow3D学习5 算例3 Hydraulics TutorialHydraulics Tutorial这个问题中流体从18cm高的坝上流下，根据自由落体运动，水流流速为Velocity = sqrt（2*980*18） = 187.8 cm/

26、sReynold数为 Re = 30cm x 187.8cm/s2 / 10-2cm2/s = 5.6 x 105因此粘性应力可以忽略，不需要用很密的网格描述墙面粘性剪切应力，但是由于存在紊流，流体内会存在粘性剪切应力，因此需要在模拟中使用粘性模型。Bond数和Weber数为Bo = 980cm/s2 * 1 gm/cc * （30cm）2/（73gm/s2） = 1.2 x 104We = 30cm * （187.8 cm/s）2 * 1gm/cc / （73gm/s2） = 1.45 x 104即不需要考虑表面张力。根据问题的性质，采用对称方式模拟。创建workspace和simulati

27、on，在General中设置结束时间为1.0s，单位CGS，在下面Notes中的第一行是标题，将出现在所有的输出文件和图形中，这里写为“Flow over a Weir”。在Meshing&Geometry tab中导入c:flow3dguistl_libweir1.stl，按缺省设置。下面增加上游水库的底面，单击工具栏Box图标，为定义不同的属性（如糙率参数），选择New Component 2，X -101，Y -2020，Z 0.08。选取计算域时上游不能太小，否则会因为突然出现的加速度引起计算不稳定，下游如果计算范围小会使边界条件影响流体流态，但是计算范围太大会使计算时间增加。这个例子

28、先用一块网格计算，然后嵌套加密的网格，网格定义可以用手工或图形方式。设置计算网格范围X -1020，Y 010，Z 018，各方向网格分别为30、10、18，为在对称处加密，在Y方向的Pt（1）选项中设置网格大小为0.5。Boundaries Tab：这个问题中X两端边界被设置为hydrostatic pressure边界条件，如果边界的“fluid height”比计算域内高则会向计算域内流入，否则会流出。在X Min边界上选择Specified Pressure边界条件同时选择Stagnation Pressure，并设置Fluid Height为15.5；在X Max中同样使用Speci

29、fied Pressure和择Stagnation Pressure，并设置Fluid Height为1.7，设置F fraction为0.0，这样可以防止流体从这个边界进入计算域，但是允许流体自由流出。（开边界？由于在X两端使用了静水压力边界，必须激活静水压力选项才能保证边界条件正确设置，在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.