本算例演示利用STAR CCM+中的Resolved VOF-Lagrangian Transition 模型模型模拟喷泉。其中涉及到液体介质从从连续相转化为离散相的过程模拟。

计算结果如下图所示。

注：本算例取自STAR CCM+的随机案例文档。

”

在STAR-CCM+ 模拟中，混合多相流是指使用多个不同多相流模型的组合对多相介质的多个流态进行模拟。这种模拟功能拓宽了现实生活中多相流的模拟范围。

在混合多相流模拟中，每个多相模型都涵盖特定的流态。通过考虑局部流动条件和网格单元尺寸的相相互作用来处理流态之间的相的转变。此类模拟有多种应用，例如车辆水管理（挡风玻璃雨刷、后视镜）、电动机中的喷油冷却或多相泵等。

1 问题描述

在本算例中，将模拟开放环境中简单喷泉的流动。以下屏幕截图显示了模拟中显示的几何形状和不同的流态：

自由表面水射流以 5 m/s 的速度离开圆柱形喷泉几何体进入空气域中，从而在水和空气之间形成大尺度的界面。射流的初级破裂产生大的连续水滴和水带。这两种流态都使用VOF 模型进行建模。

大水滴和液带分解成较小的水滴。当球形水滴的直径低于某个阈值（本算例为 5 mm）时，其会转变为使用拉格朗日多相模型跟踪其路径的拉格朗日水滴。VOF 到 Lagrangian 的转换由 Resolved VOF-Lagrangian Transition 模型处理。当拉格朗日水滴和 VOF 斑点撞击地面时，会在地面会形成薄膜水坑，可以积聚到水池中。薄水膜可以通过液膜模型以非解析的方式进行建模，或者通过 VOF 模型以解析的方式进行建模。自适应网格细化 (AMR) 用于沿水-空气自由表面局部细化网格。一旦 VOF 水滴转变为拉格朗日液滴，AMR 就会粗化网格。

本教程侧重于演示设置混合多相仿真的工作流程。

2 STAR CCM+设置

• 启动STAR CCM+2021.2并新建Simulation
• 利用菜单 File → Load加载文件 fountain_start.sim

2.1 生成计算网格

• 右键选择模型树节点 Operations → Automated Mesh，选择弹出菜单项 Execute 执行网格划分过程

计算网格如下图所示。

2.2 选择物理模型

• 右键选择模型树节点 Physics 1，点击弹出菜单项 Select Models… 打开物理模型选择对话框

选择下表所示的物理模型。

注：最核心的选项在于自适应网格与自适应时间步长。

”

选择完毕后的对话框如下图所示。

2.3 自适应网格与自适应时间步长

• 选择模型树节点 Physics 1 → Adaptive Mesh ，属性窗口中激活选项 Limit Cell Size，设置参数 Min Adaption Cell Size 为 0.003 m

• 鼠标右键选择模型树节点 Models → Adaptive Mesh → Adaptive Mesh Criteria，点击弹出菜单项 New → Free Surface Mesh Refinement 创建新的节点Free Surface Mesh Refinement

• 选中上一步创建的新节点 Free Surface Mesh Refinement ，属性窗口中设置如下参数
• 设置 Max Refinement Level 为 4
• 设置 Resolution Criterion for Interface Detection 为 10
• 其它参数保持默认设置

• 鼠标右键选择模型树节点 Adaptive Time-Step → Time-Step Providers，点击弹出菜单项 New → Free Surface CFL Condition ，其它参数保持默认设置

2.4 定义相

1、创建欧拉相Air

• 鼠标右键选择模型树节点 Multiphase → Eulerian Phases，点击弹出菜单项 New 新建流体相，重新命名为 Air

• 右键选择节点 Air ，点击弹出菜单项 Select Models… 打开模型选择对话框

• 选择如下表所示的选项

设置完毕的对话框如下图所示。

2、创建欧拉相water

• 相同的方法创建第二种相，并命名为 water，设置其Material为 Liquid，指定其Equations of State为 Constant Density

3、创建颗粒相Droplets

• 右键选择模型树节点 Lagrangian Multiphase ，点击弹出菜单项 New 新建颗粒相，修改其名称为 Droplets

• 右键选择节点 Droplets ，点击弹出菜单项 Select Models… 打开模型选择对话框

• 如下表所示选择颗粒相的模型参数

选择完毕后的对话框如下图所示。

4、创建液膜相Water Film

• 右键选择模型树节点 Fluid Film → Fluid Film Phases ，点击弹出菜单项 New 创建液膜相，命名为Water Film

• 右键选择节点 Water Film ，点击弹出菜单项 Select Models… 打开模型选择对话框，选择如下表所示的模型参数

设置完毕后的对话框如下图所示。

注：这里定义了四种相，包括Air、Water、Droplets、Water Film。

”

2.5 设置相间作用及流型转换

相间相互作用描述了空气和水两相之间的相互影响。在混合多相流模型设置中，可以使用相间相互作用来控制水向其他不同流型的转变。本算例需要考虑5种相间相互作用。

1、创建Air-Water相互作用

• 右键选择模型树节点 Physics 1 > Models > Multiphase Interaction > Phase Interaction ，点击弹出菜单项 New → Air → water
• 将新创建的节点命名为 Air-Water

• 右键选择模型树节点 Air-Water ，点击弹出菜单项 Select Models… 打开模型选择对话框

• 按下表所示选择模型参数

选择完毕后的模型对话框如下图所示。

2、创建Water-Droplet相互作用

• 按下图所示创建水和液滴之间的相间作用 Water-Droplet

• 如下图所示打开模型选择对话框

• 选择下表所示模型参数

选择完毕后的模型对话框如下图所示。

这里利用Field Function定义连续相water破碎为液滴的转化判据。大液滴（VOF滴）和液带通过VOF模型计算，液带随后分裂成液滴。只有小于指定转化长度的球形VOF液相才会转化为拉格朗日液滴。

转化判据定义为：

其中 为转化长度尺度， 为水滴 的球形度。球形度是指在相同的圆球体积下，理想球体的表面积与圆球实际表面积的比值。转化判据是指只有球形度大于0.7，直径小于5 mm的VOF液滴才能转化到拉格朗日液滴。利用VOF-Lagrangian模型处理这个转化过程。

• 右键选择模型树节点 Tools → Field Functions ，点击弹出菜单项 New → Scalar 创建函数 Blob Sphericity Criterion

• 如下图所示设置 Function Name 为 blobSphericityCriterion，设置 Definition 为 ${BlobSphericity} > 0.7 ? 0.005 : 0

• 如下图所示，选择模型树节点 VOF-Lagrangian Impingement ，设置属性 Scalar Function 为 Blob Sphericity Criterion

3、创建Droplet-Film相互作用

• 如下图所示创建液滴与液膜相间相互作用并重命名为 Droplet-Film

• 右键选择节点 Droplet-Film → Models ，点击弹出菜单项**Select Models…**打开模型选择对话框

选择以下模型选项。

选择完毕后的对话框如下图所示。

4、创建Water-Water Film相互作用

• 如下图所示创建水相与液膜相之间的相互作用，并命名为 Water-Water Film

• 如下图所示为相间作用添加物理模型

• 添加模型选项 Multiphase Material及 Resolved Fluid Film，添加完毕后模型对话框如下图所示

5、创建Air-Water Film相互作用

• 如下图所示创建空气与液膜之间相互作用，命名为 Air-Water Film

• 如下图所示为相间作用添加物理模型

• 添加模型选项 Multiphase Material及 Surface Tension，添加完毕后模型对话框如下图所示

2.5 设置边界条件

1、为液膜创建壳区域

• 右键选择模型树节点 Regions → Fountain → Boundaries → Floor ，点击弹出菜单项 Create Shell Region 创建壳体区域，软件会自动创建新的节点 Floor Shell

• 选择模型树节点 Floor shell，属性窗口中设置参数 Continuum 为 Water Film ，如下图所示

2、设置Inflow边界

• 如下图所示，设置边界 Inflow的速度为 5 m/s，设置体积分数为**[0,1]**，表示进入区域的全部为Water

3、压力出口边界

• 同时选择 Sides节点与 Sky节点下的 Volume Fraction ，指定其值为**[1.0, 0.0]** ，表示这两个边界回流介质为空气

4、颗粒边界条件

• 展开节点 Continua > Physics 1 > Models > Lagrangian Multiphase，选中节点 Lagrangian Phases > Water Droplets > Boundary Conditions > Fluid-Film Boundary > Physics Conditions > Mode ，属性窗口中指定 Active Mode 为 Fluid Film

• 选中节点 Boundary Conditions > Wall > Physics Conditions > Mode node ，属性窗口中指定 Active Mode 为 VOF Conversion

2.6 创建入射器

• 右键选择模型树节点 Injectors ，点击弹出菜单项 New 新建入射器Injector 1

• 选中节点 Injector 1，在属性窗口中设置如下参数

设置完毕后如下图所示。

2.7 显示多流型多相流流动

1、创建等值面

创建等值面用于显示空气-水分界面。

• 右键选择模型树节点 Derived Parts ，点击弹出菜单项 New Part → Isosurface…打开设置对话框

• 按下图所示顺序及参数进行设置

2、创建求解历程文件

• 右键选择模型树节点 Solution Histories ，点击弹出菜单项 New… 打开文件保存对话框，保存文件名为 fountain

• 选择模型树节点 Solution Histories > fountain，设置如下参数

设置完毕后如下图所示。

• 选择节点 fountain > Update > Time-Step Frequency ，设置参数 Frequency 为 50

3、设置Scene

• 展开模型树节点 Scenes > Multi-Regime Flow
• 选择模型树节点 Multi-Regime Flow > Water > Parts ，设置属性框参数 Parts 为 Derived Parts > Isosurface
• 设置模型树节点 Droplets ，按下表所示参数进行设置

• 设置模型树节点 Water Film ，设置下面的参数

2.8 设置计算参数

• 选择模型树节点 Solvers → Implicit Unsteady ，设置 Time-Step 为 3.0e-4 s

• 选中模型树节点 Solvers > Adaptive Mesh > Trigger > Time-Step Frequency ，设置 Frequency 为 50

• 选中节点 Lagrangian Multiphase > Two-Way Coupling > Volume Source Smoothing Method ，指定 Method 为 Cell Cluster

• 选中节点 Stopping Criteria > Maximum Physical Time ，指定 Maximum Physical Time 为 1.5 s

• 选中节点 Maximum Steps ，取消选中选项 Enable

• 选择模型树节点 Continua → Physics 1 → Initial Conditions → Volume Fraction ，设置 Value 为**[1.0 ,0.0]**，表示初始时刻计算域内全部为空气

• 点击菜单 Solution → Run 开始计算

3 计算结果

• 喷水过程如下图所示

STAR CCM+中的Resolved VOF-Lagrangian Transition 模型等同于Fluent中的VOF-to-DPM模型，是最近几个版本中才提供的新模型。

案例相关文件：

链接：https://pan.baidu.com/s/1_AzAkLgViIlIf3yQ3LCVkg 提取码：821a

文章来源：CFD之道

” 、