通过平滑颗粒流体动力学 (SPH) 模型可以访问物理模型和设置,以定义 SPH 多相模拟的属性和离散化。
表 1. 平滑颗粒流体动力学 (SPH) 模型参考
| 理论 |
多相流的平滑颗粒流体动力学 (SPH) 法。 |
| 提供方式 |
|
| 节点路径示例 |
|
| 要求 |
前提条件选择:
- 时间:隐式非定常
- 空间:三维
- 材料:多相
- 多相流模型:平滑颗粒流体动力学 (SPH)
|
| 属性 |
关键属性包括:边界公式、内核、高阶边界积分和内核相对尺寸比 请参见平滑颗粒流体动力学 (SPH) 模型属性。 |
| 激活 |
物理模型 |
自动选择的模型:
多相流相互作用(通过选择多相流模型激活)
不可压缩流体 (ISPH)
求解插值
自适应时间步
|
| 区域设置 |
物理条件: 物理值:
- 颗粒基础尺寸
- 直接旋转选项
- 动量源
- 动量源速度导数
- 颗粒生成
- 部件方向
请参见区域设置。 |
| 场函数 |
无。 |
平滑颗粒流体动力学 (SPH) 模型属性
- 边界公式
-
可用于指定用于计算颗粒-壁面相互作用的方法。
| 方法 |
对应属性 |
- 幽灵粒子法
-
选择幽灵粒子法以计算流体-壁面相互作用。
此方法最快且最稳健,尤其适用于在 GPGPU 上运行的模拟。这是默认选项。
|
无。 |
- 边界积分法
选择边界积分法以计算流体-壁面相互作用。 此方法通常更精确,特别是对于间隙小于颗粒尺寸的复杂边界。对于涉及无流体相互作用的区域的应用,这种方法的效率也更高,因为它可以减少初始化期间的 CPU 时间。但是,此方法的稳健性较差,计算成本更高。 使用边界积分法时,颗粒有时可能会与边界过于接近。启用边界保护颗粒修复选项可以减轻这种影响。请参见颗粒修复模型。
|
- 高阶边界积分
-
激活用于计算壁面边界条件的高阶边界积分方法。
请参见高阶边界积分。
|
- 内核
-
指定在 SPH 插值中使用的内核。
| 方法 |
对应属性 |
- 三次
-
选择用于 SPH 插值的三次内核方法。请参见 Eqn. (1072)。
|
无。 |
- Wendland
-
选择 SPH 插值所使用的 Wendland 内核方法。请参见 Eqn. (1073)。
建议尽可能使用此方法,因为它在精度和稳健性之间实现了很好的平衡,尤其是在结合边界积分法对壁面-流体相互作用进行建模时。
|
无。 |
- 边界积分法细分细化度
- 仅当边界公式设为边界积分法时才可用。
- 指定单个表面网格单元中的壁表面单元(Eqn. (1099) 和 Eqn. (1100) 中的
)的细分数。值越大,每个表面单元内的细分数就越多,从而提高模拟中边界条件强制执行的精度。但是,这也会导致计算成本增加。
- 高阶边界积分
- 仅当边界公式设为边界积分法时才可用。
- 激活后,Simcenter STAR-CCM+ 将应用高阶边界积分法,从而提高壁面边界建模的精度。请参见 Eqn. (1098)。
- 内核相对尺寸比
-
指定内核半径与颗粒尺寸之比。如果使用幽灵粒子法对流体-壁面相互作用建模,则此比率固定为 2.0,不可更改。
如果使用边界积分法,则可以调整内核尺寸与颗粒尺寸之比。在计算中加入更多粒子可以增大此比率,从而提高插值精度,但同时也会增加计算成本。
区域设置
以下设置适用于每个表面区域。
- 初始条件选项
- 用于自定义单个区域的初始条件。
请参见设置特定区域的初始条件。
- 动量源选项
- 用于将额外动量添加到流体域。
-
| 动量源选项 |
对应的物理值节点 |
| 无 |
无。 |
| 指定 将指定的动量源添加到动量方程。Eqn. (1060) 中的
。 用于外部加速度(如重力)随时间或空间变化的模拟。
|
- 动量源
- 用于指定动量源的矢量分布。动量源项的量纲为每单位体积的力。Eqn. (1061) 中的
。
- 动量源速度导数
- 对于作为速度函数的动量源,将动量源相对于 x-, y-, [z-] 速度分量的导数设为矢量分布值。Eqn. (1061) 中的
。
- 添加导数的结果取决于是使用显式还是隐式速度求解器时间积分公式,因为它依赖于当前速度或前一个速度。
|
- 运动指定选项
-
| 方法 |
对应的物理值节点 |
- 直接旋转运动值
-
用于直接在区域内定义旋转运动。它可用于将旋转单独添加到区域的不同部分。
|
- 直接旋转运动
- 此选项支持部件子分组,也就是说,可单独为同一区域的每个几何部件定义旋转值。
- 直接旋转运动值
- 指定旋转的属性。
请参见定义直接旋转运动。
|
- 颗粒基础尺寸
- 指定用于在域中间隔流体颗粒的颗粒尺寸(离散化长度)。此值表示颗粒生成时使用的笛卡尔栅格的目标尺寸。
- 颗粒生成
-
在表示液体初始形状和位置的几何部件内生成颗粒。求解器根据初始液体体积和颗粒基础尺寸计算颗粒总数。颗粒均匀分布并提供液体离散化。
-
部件、边界距离比和网格偏移比属性控制流体颗粒的生成。表面颗粒采样因子和表面颗粒采样优化配置幽灵颗粒或边界表面颗粒,以使边界(包括流入边界)与流体颗粒相互作用。
-
- 部件
- 选择几何部件以在其中生成定义初始流体域的粒子。几何部件的表面必须没有错误且完全闭合。
- 边界公式
-
只读属性,用于指示在节点下选择的对颗粒-壁面相互作用建模的方法。
- 边界距离比
- 指定粒子中心与初始流体体积边界表面之间的初始最小距离比。
- 网格偏移比
- 指定初始化时粒子中心与流体原点之间的偏移比。网格偏移比与笛卡尔模式一致,以便适合块状初始流体体积。
-
下图说明了栅格偏移比和边界距离比的基本概念。这些值与初始流体体积的表面(红色)相关,偏移距离确定了实际颗粒体积(绿色)。通过将网格偏移和边界距离除以颗粒尺寸可获得比率。
要增强模拟初始阶段的稳定性,调整这些属性以在流体域内稍微向内移动颗粒。
- 表面颗粒采样因子
-
指定颗粒采样因子,用于控制何时停止生成其他幽灵颗粒和边界表面颗粒。通过 SPH 模型可调整表面颗粒采样,以创建沿边界(包括流入边界)分布的其他幽灵颗粒,以支持与流体颗粒的相互作用并防止出现孔。
采样过程分批(随机点集)继续进行,直到生成的样本总数不再增加到超过指定的改进因子为止。默认值 0.001 对于大多数模拟已足够;增大该值会生成更多颗粒,但会导致计算成本提高。
- 表面颗粒采样优化
- 激活后,Simcenter STAR-CCM+ 支持使用分边细化方法进行附加细化优化,以实现最大表面颗粒采样。建议为长宽比值大的网格激活此选项,以提高计算精度。激活此选项也会增加计算成本。
- 部件方向
-
提供用于反转表面法向矢量相对于流体位置的方向的选项。默认情况下,闭合几何部件的表面法向朝内,即边界的正面朝外。SHP 颗粒仅在正面与边界相互作用。因此,对于闭合几何部件周围的外部流体,默认方向是正确的。但是,对于几何部件内的内部流体,默认方向不正确,必须反转。
-
- 根据部件子组指定
- 为激活允许每个部件值时可用。
- 激活时,可以对分配给同一表面区域的各个几何部件进行子分组。根据部件子组节点显示在模拟树中,同时具有预定义的子组节点。
- 反转方向
- 激活时,它会反转子组部件的表面法线的方向。