Fluent VOF罐体晃动（二）_Fluent_UDF_Origin-仿真秀干货文章

Fluent VOF罐体晃动（二）

7小时前浏览10

本案例利用Fluent中的VOF模型，对罐体晃动问题进行了仿真计算。

具体的操作与上次推文的罐体晃动（一）一致，只是不再采用命名表达式的方式进行罐体晃动仿真，而是通过UDF编译，本推文主要对UDF的编译和加载进行了介绍。该方法的缺点就是无法进行多工况的快速计算，优点则是在开启能量方程等模型时，能够通过UDF统一编译进行处理，提高计算效率。

UDF

1.1 UDF编译

部分编译如下，

注意：使用时需检查符号，因在录入代码时采用手敲，不能保证输入法正确，因此大家录入时所有代码要确保全是英文字符。

DEFINE_ZONE_MOTION(vel，omega，axis，origin，velocity, time，dtime )

{

real Am=0.2；

real f=1；

vel[0]=Am*f*2*PI*cos(2*PI*time)；

}

1.2 UDF加载

加载操作如下，首先选择Build，待编译成功后，选择加载。加载成功后，在运动区域的设置中将命名表达式去掉，采用UDF定义运动。

其他设置与上一篇文章完全一致，因此不再阐述。

来源：CFD仿真库

Fluent Eulerian 风沙两相流

本案例利用Fluent中的Eulerian(欧拉)模型，对防风挡沙墙问题进行了仿真计算。在防风挡沙墙问题中，需要利用入口风速廓线对入口速度进行定义。本文主要对Flunet命名表达式的方法进行仿真计算。1FLUENT设置1.1导入网格计算域为150m×20m，防风挡沙墙高H为2m，厚度为0.3m，放置于距入口处30m处，如图网格由Icem软件生成的结构化网格。1.2General设置由于是风沙问题，此处设置为瞬态计算，并开启重力。1.3材料定义此处添加材料为自定义材料，具体的材料参数如下图所示。1.4模型设置采用k-wSST湍流模型（实际仿真需要使用合理的湍流模型），并开启欧拉多相流模型。欧拉模型设置如下。2定义命名表达式通过入口风廓线（根据实际需求设置，此处参照文献进行设置）进行入口速度设置，命名表达式的使用可以参照前文FluentVOF罐体晃动（一）中的使用。设置结束后的风廓线云图示意图，可以看到从下到上速度成对数形式逐渐增大。3边界设置上壁面设置为对称面，下壁面设置为WALL，砂石部分设置为滑移壁面。左侧设置为速度进口，右侧设置为压力出口。其中速度进口处分为风沙混合进口与净风进口，以2m的高度为界。所有速度入口选择表达式构建的速度信息。风沙入口处设置砂石体积分数为0.02。4初始化进行初始化。5后处理云图结果对风沙运动情况进行监测，速度分布情况如下。来源：CFD仿真库