FLUENT操作过程及全参数选择

振动流化床模拟工艺及参数选择1生成流化床模型。根据陈海波论文中提供的试验机参数构建流化床模型。流化床直148毫米，高1米，开口率9%，光圈2毫米。在屏幕上铺两层帆布，使气流均匀。实验机是圆形的流化床，所以只能用二维模型简化。在实际实验中，中间化高度明显小于1m或500mm，因此，为了提高计算时间，可以将模型高度减少到500mm。为了使气流均匀分布，屏幕上铺了两层帆布，所以我认为沿整个屏幕进口风速均匀。最终简化模型如下图所示。上图显示流动流化床模型表明，流化床底部末端的网格比较贴合在顶部。流态化行为主要发生在流化床底部，由于加快计算时间，因此采用了这种紧密划分上梁的方式。这里的收入是velocity inlet设置为。出口流出；设定为；左右分隔设置为wall。在GAMBIT中进行设置，然后导出2d模型vfb.msh。Outflow边界条件不必指定入口的物理条件，但其应用限制在以下四个方面：1.仅适用于不可压缩流2.出口的流向完全发展了3.不能与压力边界条件(压力入口、压力出口)一起使用4.不能用于计算流量分配问题(例如，具有多个出口的问题)2打开FLUENT 6.3.26，然后导入模型vfb.mshGRID单击CHECK以查看网格信息和模型中设置的信息，并验证其是否正确。特别是检查是否有负区域和负区域(例如，立即更正)。完成后，单击GRID-SCALE显示SCALE GRID窗口，将单位设置为mm，然后单击change length unit按钮。具体设置包括：3设定解析器实际液化过程会随时间变化，因此请将其他设置保留为默认值，并将TIME更改为unsteady。(1)在求解不可压缩流体时，在动量方程和连续性方程中不连续求解代数方程，可以直接获得角速度分量及其压力值，但由于使用了大量计算内存，因此可用在Fluent6.3中，所需内存是分离算法的1.5-2倍。Density based解决方案是针对压力流体而设计的，因此更适合于将速度分量、密度(基于密度)作为主变量的压力场计算。压力由状态方程求解。基于预应力的求解器使用压力补偿算法，解决的控制方程是标量形式，可以很好地解决不可压缩流，可以解决可压缩流。Fluent 6.3之前的Solver版本只有Segregated Solver和Coupled Solver，这是基于Pressure的Solver的两种处理方法。基于Density的解算器是在Fluent 6.3中新开发的，所解释的控制方程是矢量形式，主要离散格式Roe、AUSM基本上设计为Fluent可以更好地解决可压缩流动性，但目前格式中没有添加限制器，因此还不完善。只有Coupled的算法。对于低速问题，使用Preconditioning方法进行处理，以便计算低速问题。Density-Based Solver下显然没有SIMPLEC，PISO这些选项。因为它们是这种类型求解器中未出现的压力修正算法；通常，使用基于Pressure的解决方案进行故障诊断。(2)在下面的“渐变选项”(GRADIENT OPTION)选项组中，指定用于计算控制方程导数的压力拔模。CELL-BASED(按单元格的压力梯度计算)和NODE-BASED(按节点的案例渐变计算)。“Porous formulation”选项组用于开发多孔介质的速度。(3)选择UNSTEADY后，将显示UNSTEASDY FORMULATION选项组，从而使用户可以基于顶级时间相关主题的计算公式和方法。对大部分巨大的问题做一次隐式选择就足够了。只有在对精度有特殊要求的情况下，才选择二次隐式。设置4多相流模型。如果设置为Euler模型，将相位数设置为2，则为两相流，如下所示：Fluent有三个Euler-Euler多相流模型：Volume Of Fluid(VOF)模型、Mixture模型和Euler(Euler Ian)模型。(1) VOF模型。VOF模型是固定Euler栅格下的曲面跟踪方法。如果需要一个或多个不相互融合的流体之间的连接面，则可以使用此模型。在VOF模型中，其他流体元件共用一组动量方程式，并记录计算时每个流体元件在整个流场的每个计算单位内所占的体积百分比。VOF模型的应用示例包括层流、自由表面流、灌注、晃动、液体中大气泡的流、坝破裂时的水流、求任意液-气界面的稳态或瞬时分割界面。(2)混合物模型。混合物模型可用于两相流或多相流(流体或粒子)。在欧拉模型中，每个相位都被视为互相穿透的连续体，因此混合物模型求解混合物的动量方程，并以相对速度描述不连续的相位。混合物模型的应用包括低负荷颗粒载荷流、气泡流、沉降和旋流器。混合物模型也可以用于没有离散相相对速度的均匀多相流动。(3)欧拉伊恩模型。Fluent中最复杂的多相流模型。建立n个动量方程式和一组连续方程式来解决每个拓朴。压力和界面交换系数相互结合。合并方法依赖于粒子流(流-固体)和非粒子流(流-流)不同的包含阶段。欧拉模型应用气泡柱、浮动、粒子悬浮和流化床。根据振动流化床的实际情况，本论文采用欧拉模型进行仿真。5设置具有粘性的模型。第一步，依次单击DEFINE-MODELS-VISCOUS、VISCOUS MODEL对话框弹出菜单、选择K-EPSILO模型和确定。在第二步中，在任务窗格中输入以下命令：define/models/viscous/turbulence-expert/low-re-k屏幕显示：/define/models/viscous/turbulence-expert low-re-kenable the low-re k-Isilon tor bulence model？否输入y以在“选择模型”面板上查看低射线模型low-re-ke模型。默认情况下，使用第0低雷诺数模型。第三步是在Fluent中提供六个低雷诺数模型，然后使用low-re-ke-index命令设置一个。Low-re-ke-index索引模特0Abid1Lam-Bremhorst2装载机-沙尔马3yang市4Abe-Kondoh-Nagano5Chang-Hsieh-Chen默认情况下，此仿真使用第0个低雷诺数模型。标准k-epsilo模型是使用非常充分的湍流流开发湍流的模型创建的，它是高雷诺数的湍流计算模型，与0方程模型和1方程模型相比，有很大的改进，但是在强涡流、弯曲壁面流或弯曲流流中会发生扭曲。与标准模型相比，RNG k-模型修改了湍流粘度，考虑了平均流的旋转和旋涡流，可以更好地处理高应变速率和流线型弯曲成都的大流。此外，它还针对完全发展的湍流，即高雷诺数模型。Realizable k-模型通常应用于具有射流和混合流的自由流、管道中的流、边界层流等。实际计算的雷诺数较少，因此与上述三个湍流模型不匹配。FLUENT提供了标准k-的一些专家模型，对其进行修改以适合低雷诺数(低雷诺数k-epsilo模型)。6定义材料特性。DEFINE-MATERIALS、“材料弹出对话框”、“创建”按钮，首先选择“空气”作为“气象”。然后单击FLUENT DATABASE MATERIALS按钮，从材质库中选择流体，然后单击COPY按钮。然后将材质的密度和名称更改为所需材质的材质特性，如下所示，然后单击“更改”。7正义奖。DEFINE-PHASE。首先，定义空气的主要阶段，如下所示：然后将二次相设置为固相米勒。订购Millet后，单击SET按钮弹出secondary phase对话框，继续进行下一设置。首先，将材质定义为颗粒，即粒子，以定义粒子大小。包装的bed是填充床，与实际不符，因此未选择。选择粒子温度模型PHASE PROPERTY相属性。Partical differential equation是偏微分方程。固体剪切粘度包括碰撞和动力部分、摩擦部分。其中，“增强”部分提供了两种表达。默认值为sya mlal表示和GIDASPOW ET AL表示，通过实验选择sya mlal表示。固体体积粘度被解释为对粒子压缩和膨胀的抵抗力，对于这个项目一般没有争议，目前学术界广泛采用Lun et al这个说法。这篇论文的模拟忽略了摩擦粘度。填充限制设置为0.6。也就是说，初始固体上的体积分数最高为0.6。设置气固封闭关系：在“裁切”对话框中单击INTERACTION，然后设置气固相互作用的抗力函数一般为Wen-Yu、GIDs pow、SYAMLAL-OBRIEN三种，实验结果表明sya mlal-obrien更真实选择Syamlal-Obrien恒力函数模型。(1)Syamlal-OBrien型号234(20.4.31)以下是Dalla Valle47提供的阻力函数：(20.4.32)此模型是根据流化床或沉降床粒子的结束速度测量的，体积分数与相对雷诺数的函数关系193:(20.4.33)这里下标是第一液体相，第一固体相，第一固体相粒子的直径。液固交换系数为(20.4.34)这里是与固体相关的端速度73:(20.4.35)其中(20.4.36)是，(20.4.37)是，(20.4.38)此模型适用于根据Syamlal et al定义固体阶段的剪应力时的235(方程式20.4.52)。(2)对于温家宝和俞模型262，液固交换系数的形式如下：(20.4.39)在这里，(20.4.40)数量由方程式20.4.33定义。该模型适用于稀释系统。(3) GIDs pow模型76是Wen and Yu模型262和Ergun方程62的组合。当时的液固交换系数为：(20.4.41)这里(20.4.42)当时，(20.4.43)对于高密度流化床，建议使用此模型。由于此流化床中粒子的直径大于粒子之间的距离，因此不适合将正向力包含在几乎充满的粒子中，请忽略正向力的影响，然后从“左”选项中选择“无”。在“恢复系数”选项下，保留默认设置值0.9。虚拟质量力可以忽略，因为第二相密度大于第一相。具体设置包括：8编译UDF程序。define-user-defined-function-compiled，导入程序。(1) voiddefine _ CG _ motion (UDF name，dynamic _ thread * dt，realvel ，real omega ，real omega此函数界面用于控制刚体的移动，允许动态边界控制，方法是将刚体中心运动速度和每个速度分别指定给vel和omega，fluent根据相应的值自动计算边界下一步的位置。刚体质心的位置可以在“函数界面”对话框中定义。Dynamic Zones中的dwall是要控制的动态边界，Motion UDF/Profile中的stc1sta010a0ph0是UDFname。在此开发的目的是控制dwall，理论上，FLUEN T可以控制无限数量的动态边界。C.G.Location用于设置初始位置的质心，C.G.Orientation用于设置刚体的初始角度。通常适用于刚体本身不变形的运动。2) void define _ geom (char name，domain * d，dynamic _ thread * dt，real * position)。此函数界面用于控制变形实体的边界行为，变形实体是运动边界中网格节点的位置值，可以指定此值以实现控制效果。position 0对应于边界节点的x坐标，position 1对应于y坐标，position 2对应于z坐标。FLUENT适用于规则运动边界，可由函数描述，因为FLUENT会自动通过所有边界节点。3) void define _ grid _ motion (name、d、dt、time、dtime)。此函数界面还用于控制形状的边界行为。UDF编程比前面两种情况