CFD丨多相流案例丨03VOF：空化

Shmily89

2月前浏览2985

本案例为STAR-CCM+帮助文档案例学习，只是笔记记录！

使用的网格同CFD丨多相流案例丨01VOF：重力驱动流体

只是此例中的相关物理维度均缩小 1000 倍。因此，喷嘴宽度约1mm。相关几何体左侧边界的规定压力为 5e+07 Pa，而右侧边界处于大气压力下。底部的边界是一个对称平面，并且所有其他边界都是壁面。

最初，左腔室注满水，而求解域的其余部分则注满空气。使用默认的 K-Epsilon 模型为湍流建模。窄喷嘴入口处的尖角与流体加速结合产生一个可发生空化的低压区域。

1 导入网格文件

① 启动STAR-CCM+。

② 选择文件→导入→导入体网格。

③ 保存文件。

2 转换为2D网格

① 选择网格→转换为2D。

② 确保已激活转换后删除3D区域选项。

③ 删除连续体→物理1。

注意：三维网格转换成二维网格要求：

a. 网格必须在X-Y平面对齐。

b. 网格必须在Z=0位置有一个边界平面。

3 缩放网格

初始网格并未按正确比例构建，因此需要缩小 1000 倍。

① 从菜单中选择网格→比例缩放网格。

②在比例缩放网格对话框中，选择 Default_Fluid 2D 区域，比例因子设置为0.001。

③ 单击应用后，网格区域尺寸会缩小。

④ 单击可视化工具栏中的重置视图。

4 选择物理模型

物理模型定义模拟的主变量，包括压力、温度、速度和用于生成求解的数学公式。

在该模拟中，流体是湍流且问题涉及多相流体和空化。

此次分析需要三种流体（空气、水和水蒸气）。

但是，由于这些流体占据相同的域，所以仅需要一个连续体和一个区域即可设置模拟。

① 将physical 1 2D 连续体节点重命名为喷射器。

② 选择物理模型如下图所示。

5 设置材料特性

在欧拉相定义各混合物成分所对应的材料。

① 在喷射连续体中，右键单击模型→多相→欧拉相，然后创建一个新相。

② 将相 1 节点重命名为 H2O。

③ 对H2O相，选择模型。

④ 对于 H2O 相，请选择下列模型。

⑤ 创建第2相，然后将其重命名为air，air空气模型选择。

⑥ 创建第3相，然后将其重命名为H2O（g），H2O（g）模型选择。

注意：

如果您要设置标准自由表面问题，则仅需两相。

但由于此问题涉及空化，所以必须定义空化产生的第3相，即水蒸气。

因此，在STAR-CCM+ 材料数据库中将第3相中的空气替换为水蒸气。

⑦ 右键单击 H2O (g) → 模型→ 气体→ 空气节点，然后选择替换为。

⑧ 在替换材料对话框中，展开材料数据库 → 标准 → 气体→ H2O(Water)。

6 定义相间相互作用

使用多相交互作用模型可定义空气和液相之间的相互作用。

要建模空化，可使用 VOF 多相流相互作用模型和 Schnerr-Sauer 空化模型。

注意：

当液体内特定位置的静压低于液体的饱和蒸汽压时，液体就会经历一种称为空化的相变。这种相变产生了充满液体蒸汽的空化气泡。

创建相间相互作用时，首先选择主相，然后选择次相。在此模拟中，主相为 H2O，次相为 H2O (G)。

① 在喷射器连续体中，右键单击模型→多相流相互作用 →相间相互作用节点，然后选择新建→H2O→ H2O (g)。

定义适当的 Schnerr-Sauer 空化模型参数。在相关的液相材料特性节点设置饱和压力。

② 选择模型→多相 →欧拉相→ H2O→模型→ 液体→H2O→材料属性→饱和压力→ 常数，设为 2338.0 Pa。

7 设置初始条件

物理连续体中的两个流体空间分布的初始条件是：只在左腔室中注入水，在右腔室和通道内注入空气。指定这种分布的便捷方法是创建并使用场函数。

定义场函数，用于指定连续体中两个流体的空间分布的初始条件：

① 工具→场函数→新建→标量。

② 将User Field Function 1重命名为初始分布（空气）。

③ 定义初始分布，设置属性。

④ 重复以上所列步骤，使用称为初始分布（H2O）的第二个场函数指定 H2O 相的初始分布，定义如下：

注意：

由于这两个函数互相可逆，因此所有区域的水蒸气浓度都为零。

因为静止流体中不存在空化，所以这种情况下此值是合适的。

设置物理连续体的初始条件和基准值。

物理连续体→初始条件→体积分数，进行体积分数设置。

8 设置边界条件

定义边界类型并指定适当的属性值。本教程中所用的几何体有六个边界，其中三个有为其分配的无滑移壁面条件。剩下的三个边界中，两个定义为压力出口边界，而另一个则定义为对称平面。

① 将 Default_Fluid 2D 节点重命名为流体。

② 流体→边界。

③ left→物理值，体积分数为1:0：0，压力为5E+07。

④ right→物理值，体积分数为0:1：0。

9 设置求解器参数和停止条件

如果要解算非稳态问题，必须指定时间步长和消耗的模拟时间。

使用时间步长 10e-7 s 运行此计算 5 x 10e-5 s，所以需要 500 个时间步。

① 求解器→隐式非定常→时间步→1e-07s。

要设置求解器参数和停止条件。

② Mmaximum physical time

10 可视化和初始化求解

创建标量场景显示模拟结果，查看空气和水的分布情况。

① 新建场景→标量。

② 选择体积分数 > H2O。

和预期一样，在运行初期，整个左腔室注满水。将标量函数更改为空气体积分数显示右腔室和连接通道都完全注满空气。在两个流体之间的交界面处可以看到有一块小区域明显存在两种流体，这种效果是由于网格的粗糙度造成的。

11 运行模拟

模拟准备现已结束，可以运行模拟。

单击求解工具栏中的运行。

输出窗口中会显示求解进度。图形窗口中会自动创建残差屏幕，并在其中显示求解器的进度。您也可单击图形窗口顶部的标量场景 1 选项查看模拟进度。

创建场函数绘图。

① 工具右键单击→场函数，选择新建→标量。

② 将用User Field Function 1重命名为体积分数（所有相）。

③ 选择体积分数（所有相）。

a. 将函数名更改为体积分数（所有相）。

b. 要打开自定义编辑器，请单击定义框。

c. 在编辑器窗口中，输入下列语句：

${VolumeFractionH2O}+2*${VolumeFractionH2O(g)}

④ 在标量场景 1 ，选择体积分数（所有相）。

合理的近似分布情况为，蓝色和蓝绿色区域主要是空气，绿色区域主要是液态水，而黄色、橙色和红色区域则主要是水蒸气。

以这种方式处理结果可以定性目测各相的空间范围。

但是，在空气和水蒸气一起出现的区域，所得到的绘图可能不精确。液态水体积分数稍大于水蒸气体积分数的区域也可能不精确。在此类情况下，后者会导致水蒸气高浓度区域的大小预测值过大。

来源：认真的假装VS假装的认真

好书推荐10余本

第1类 WB学习第1本：周炬，苏金英. ANSYS Workbench有限元分析实例详解（静力学）+（动力学）这两本书，已经推荐太多次，抽奖两次奖品都有它，足以说明我对这本书的喜爱。热分析据说也快出了，十分期待！第2类学科基础类第2本：Yunus A Gengel，John M. Cimbala.流体力学基础及其工程应用（上册+下册）第四版如果，没学过流体力学，想快速入门，一定要看这本。写的太具体太形象了！第3本：符松，王亮.湍流模式理论这本，很贵，很薄，很凝练，整体行文脉络极好，如果看那些经典的CFD理论书籍看崩了，翻这本会不会能缓口气？第4本：Yunus A Gengel.传热学（heat transfer：a practical approach）如果，看陶文铨那本，看蒙了，可以借助这本。和流体力学那本一样形象，才发现作者是同一人？第5本：张雄.有限元法基础这本，刚看了个开头，有限元法准备上路的朋友可以翻来瞧一瞧。（其实我是更风买的，一时半会看不上）。第6本：于靖军，毕树生等著.柔性设计柔性机构的分析与综合是关于柔性设计的基础理论书籍，内容非常详细，写作脉络非常的高教社。第3类光机结构类第7本：凯斯.B.道尔.光机集成分析（第2版）如果每个写软件的人，都写出这样一本书，真的是软件学习者的福音！学习sigfit之前，一定要看它！第8本：John E Greivenkamp. Field Guide to Geometrical Optics这本，全是精华，好书中的好书！第9本：Katie Schwerz,James H Burge.Field guide to Optomechanical Design and Analysis推荐理由，同上。好书中的好书！第10本：Alson E.Hatheway.The Optomechanical Constraint Equations Theory and Applications.这本只有137页，绝对是今年遇到最爱的专业书籍，精品！搞光机结构分析一定要看。（你会来感谢我的）！过段时间会认认真真，详详细细，一字一字扣着写一下这本读书笔记，因为真的太爱了。爱死国防工业出版社，引进了，翻译质量还不错！今天是世界读书日，愿，我们，都可以在书中找到想要的答案，或者可能，被点亮一下！来源：认真的假装VS假装的认真