1. 三种空化模型介绍

恒温下的液体当内部压力较小时低于饱和蒸汽压,液体可能会气化形成空腔,这一过程称为空化cavitation。在这种过程中,空化区会发生非常陡峭的密度变化。

Fluent中提供了三种空化模型:

Singhal et al. model:适用多相流模型为mixture。这种空化模型需要使用TUI命令打开,TUI命令:solve/set/expert

输入上面的TUI命令,在控制台弹出下面选项,其中Singhal et al. model选择yes

在Phase Interaction界面会出现Singhal-et-al Cavitation Model

Zwart-Gerber-Belamri model:多相流模型可选择 mixture 和Eulerian

Schnerr and Sauer model:默认的空化模型,多相流模型可选择 mixture 和Eulerian

2. 空化模型的特点

特点1:系统必须包含液相和气相,液相和气相之间发生传质。在空化模型中,同时考虑了气泡的形成(蒸发)和坍塌(凝结)。

特点2:对于空化模型,Fluent定义从液体到蒸汽的传质为正值。由于fluent提供的空化公式按照这种方式定义传质的正负,因此空化模型的传质必须从液相到气相。

特点3:空化模型是基于Rayleigh-Plesset 方程,描述了液体中单个蒸汽气泡的生长。

特点4:Singhal et al.引入了不可凝气体,假设不可凝气体的质量分数是已知常数。Singhal et al模型要求主相为液体,次相为蒸汽,只适用于mixture模型;Singhal et al.模型只能用于单一的空化过程。

也就是说,当定义传质时,只有一个液体可以发生空化。这种液体可以是一种相,也可以是混合相中的一种物质。

特点5:Zwart-Gerber-Belamri模型和Schnerr and Sauer模型不考虑不可凝气体的影响,能够与Fluent中所有可用的湍流模型兼容。Singhal et al.模型与LES湍流模型不兼容。

特点6:基于密度求解器density-based不能使用空化模型

特点7:多相流模型使用Explicit VOF模型时,不推荐使用空化模型

特点8:当使用隐式VOF公式时,默认添加湍流效应引起的数值扩散。这种扩散会增强计算的稳定性,但会造成界面精度下降。

3. 三种空化模型理论

两相流连续方程

其中,源项R表示蒸气蒸发或冷凝的速率。空化理论的关键点就是给出相变率的表达式,不同的理论给出的R表达式不同。

3.1 Singhal et al. Model理论

Singhal et al. 空化模型基于”full cavitation model”,提出了一个以蒸汽质量分数为输运方程中因变量的模型

其中,fv表示蒸气质量分数;fg表示不凝性气体的质量分数;Г为扩散系数;Re为蒸发速率,Rc为冷凝速率

如果压力小于饱和蒸汽压,则蒸发:

如果压力大于饱和蒸汽压,则冷凝:

其中,Fvap为蒸发系数,Fcond为冷凝系数,Fvap=0.02,Fcond=0.01

此模型需要设置不凝性气体的质量分数,即上面的参数fg

3.2 Zwart-Gerber-Belamri Model模型

假设系统中所有的气泡具有相同的大小,空化率R可用气泡数密度和单个气泡的质量变化率相乘得到

经过代入推导可得到最终表达式:

如果压力小于饱和蒸汽压,则蒸发:

如果压力大于饱和蒸汽压,则冷凝:

其中,

ℜB=1e-6,为气泡半径bubble radius

αnuc=5e-4,为成核位点体积分数nucleation site volume fraction

Fvap=50,为蒸发系数evaporation coefficient

Fcond=0.01,为冷凝系数 condensation coefficient

此模型设置如下:

Bubble Diameter:气泡直径,即ℜB的2倍

nucleation site volume fraction:成核位点体积分数,即αnuc=5e-4

evaporation coefficient:蒸发系数Fvap=50

condensation coefficient:冷凝系数Fcond=0.01

3.3 Schnerr and Sauer Model模型

Schnerr and Sauer 提出的模型只需要确定气泡的数量密度,气体的参数如气泡的直径、成核位点体积分数可通过此模型自动推导出,不必设置。

空化率R表达式如下:

如果压力小于饱和蒸汽压,则蒸发:

如果压力大于饱和蒸汽压,则冷凝:

其中,

Fvap=1,为蒸发系数evaporation coefficient

Fcond=0.2,为冷凝系数 condensation coefficient

这两个值为模型的默认值,在模型设置里无法更改,想要更改需要使用TUI命令

修改冷凝系数TUI命令:

solve/set/multiphase-numerics/heat-mass-transfer/cavitation/schnerr-cond-coeff

修改蒸发系数TUI命令:

solve/set/multiphase-numerics/heat-mass-transfer/cavitation/schnerr-evap-coeff

使用以上命令,控制台会出现如下文本:

4. 湍流因子Turbulence Factor

对于Schnerr-Sauer和Zwart-Gerber-Belamri模型,可以考虑湍流对于饱和压力的影响

也就是虽然我们设置了空化饱和压力,但是由于湍流的影响,这个压力值不一定准确,需要进行一定的修正。修正公式如下:

其中,

kl为湍动能;

coff为湍流系数,即上图设置中的Turbulent Coefficient;推荐值为0.39

5. 空化模型的使用依据

a. 模型的选择

Zwart-Gerber-Belamri和Schnerr and Sauer 模型的稳定性更强,收敛速度更快,Fluent强烈推荐使用这两个模型。Singhal et al模型收敛性较差,不推荐使用。

b. 求解器的选择

Fluent中无论分离式求解器(SIMPLE, SIMPLEC, and PISO)还是耦合式求解器coupled 都适用空化模型。

耦合求解器更稳定,收敛速度更快,特别是对于旋转机械中的空化流动,如液体泵、诱导器、叶轮等,因此推荐使用耦合式求解器Coupled

对于Singhal模型,耦合求解器没有任何显著的优势,建议使用分离求解器。

c. 初始条件

对于Zwart-Gerber-Belamri和Schnerr and Sauer 模型,初始条件对收敛性影响不大,因此任何初始条件都可以

对于Singhal模型,由于收敛性不好,因此对初始条件比较敏感,Fluent建议初始化时,蒸汽分数设置为入口值。初始化压力设置为入口和出口之间的最高压力

d. 压力离散格式

推荐使用如下格式

PRESTO!

body force weighted

second order

而standard和linear尽量避免使用

五十六、Fluent空化模型理论恒温下的液体当内部压力较小时低于饱和蒸汽压,液体可能会气化形成空腔,这一过程称为空化cavitation。在这种过程中,空化区会发生非常陡峭的密度变化。https://mp.weixin.qq.com/s/9LVw5rbobHBU-whxaK32MQ

以上案例的cas和dat文件均可以免费获取,需要的朋友,只需要点赞关注收藏一键三连后私信我即可哦

五十六、Fluent空化模型理论相关推荐

  1. Python编程基础:第五十六节 filter函数Filter

    第五十六节 filter函数Filter 前言 实践 前言 filter函数的作用是对可迭代对象内部的元素按照特定条件进行过滤,其书写方式为filter(function, iterable),第一个 ...

  2. OpenCV学习笔记(五十六)——InputArray和OutputArray的那些事core OpenCV学习笔记(五十七)——在同一窗口显示多幅图片 OpenCV学习笔记(五十八)——读《Mast

    OpenCV学习笔记(五十六)--InputArray和OutputArray的那些事core 看过OpenCV源代码的朋友,肯定都知道很多函数的接口都是InputArray或者OutputArray ...

  3. JavaScript学习(五十六)—寄生式继承(临时构造器的使用)

    JavaScript学习(五十六)-寄生式继承(临时构造器的使用) 寄生式继承(临时构造器的使用): 之所以要使用临时构造器就是要在实现原型链继承时获取上级构造方法原型中的内容,而不要它的私有属性 默 ...

  4. 合肥工业大学机器人技术五十六题

    合肥工业大学机器人技术五十六题 题目要求 //拿球后行为,利用已有 Worldmodel(21) (1)在 playOn 模式下,拿到球以后朝前方快速带球. (2)在 PlayOn 模式下,拿到球以后 ...

  5. MFC Windows 程序设计[五十六]之开启子对话框(附源码)

    MFC Windows 程序设计[五十六]之开启子对话框 程序之美 前言 主体 运行效果 核心代码 逻辑分析 结束语 程序之美 前言 MFC是微软公司提供的一个类库(class libraries), ...

  6. 达芬奇密码 第五十五章 第五十六章

    达芬奇密码 第五十五章 第五十六章[@more@] 第五十五章 索菲靠着兰登坐在长沙发上,喝着茶吃着烤饼,享受着食物的美味.雷·提彬爵士微笑着,在炉火前面笨拙地踱来踱去.假肢敲在地面上,发出" ...

  7. 【Visual C++】游戏开发五十六 浅墨DirectX教程二十三 打造游戏GUI界面(一)

    本系列文章由zhmxy555(毛星云)编写,转载请注明出处. 文章链接: http://blog.csdn.net/poem_qianmo/article/details/16384009 作者:毛星 ...

  8. 中国五十六个民族列表

    以后要养成一个习惯,把常用的代码保存到自己的 "代码库",虽说资料网上都有,但是有些经过自己验证的代码,如果能自己建立自己的库的话,将来也是一笔财富.所以,今天给自己的博客添加了 ...

  9. [Python从零到壹] 五十六.图像增强及运算篇之图像平滑(中值滤波、双边滤波)

    欢迎大家来到"Python从零到壹",在这里我将分享约200篇Python系列文章,带大家一起去学习和玩耍,看看Python这个有趣的世界.所有文章都将结合案例.代码和作者的经验讲 ...

最新文章

  1. Thymeleaf中each标签遍历list如何获取index
  2. SpringBoot应用和PostgreSQL数据库部署到Kubernetes上的一个例子
  3. mysql group原理_MySQL Group By 实现原理分析
  4. sql server 里面怎么支持数字使用双引号_国查:用中文编写SQL
  5. Linux查看文件的首个字母 文件属性字段
  6. 求1+2!+3!+......+20!的和
  7. C/C++基础面试-Const的全面理解(C部分)
  8. python怎么创建函数_Python创建与调用函数
  9. 反卷积,上采样,上池化的理解
  10. 疯狂讲义java_《疯狂Java讲义》 1-概述
  11. matlab 加权回归估计_Matlab:地理加权回归基本操作
  12. Ubuntu 13.10 安装后配置
  13. css找某个元素的下个子元素,使用CSS获取特定位置的子元素
  14. 把计算机知识列表合为一列,怎么把相同表格的数据合并
  15. android 强制退出当前应用程序,android 强制退出应用
  16. 理解COM的线程套件(转)
  17. Python---Excel文件xls格式转为xlsx格式
  18. xp系统无打印机服务器,WinXP局域网共享及网络打印机无任何网络提供程序怎么办?...
  19. VBoxGuestAdditions下载地址
  20. 计算机counta函数怎么用,counta函数的使用方法及实例 看完你就知道了

热门文章

  1. 拿好!52个实用的数据可视化工具!
  2. 含金量较高的数据分析证书有哪些
  3. 北京、上海、福州新基建三年规划比较: 区块链地位、作用大不同
  4. python b站 排行_用python爬虫追踪知乎/B站大V排行
  5. 凝汽式汽轮机调节系统
  6. 浅谈PE(私募股权)业务及系统建设:02. 系统建设
  7. 【狮子数学】chapter3-07-积分中值定理(第66讲)
  8. 计算机军中绿花谱子,军中绿花的曲谱_军中绿花简谱
  9. 大数据开发平台:数仓组件Apache Kylin详细安装暨使用教程
  10. C语言实现四舍五入?