许可优化
许可优化
产品
产品
解决方案
解决方案
服务支持
服务支持
关于
关于
软件库
当前位置:服务支持 >  软件文章 >  Fluent PBM | 由于成核和颗粒生长引起的相间传质

Fluent PBM | 由于成核和颗粒生长引起的相间传质

阅读数 3
点赞 0
article_banner


1.设置

在涉及颗粒产生、溶解或生长(如结晶)的应用中,由于这些现象,颗粒相的总体积分数方程将具有源项。由于质量的增加,颗粒相的动量方程也会有源项。在Ansys Fluent中,可以使用UDF宏DEFINE_HET_RXN_RATE(参见异相反应部分该宏介绍)或使用相位交互选项卡指定质量源项。


    Fluent相间作用|多相反应
   


     专栏
   

b占余文乐

8755阅读44回复

例如,在结晶过程中,粒子是通过成核(n0)产生的,也可以指定生长速率(G)。所有尺寸颗粒的传质速率(in kg/m³-s)为

%5Cbegin%7Baligned%7D%0A%5Cdot%7Bm%7D%20%26%20%3D3%20%5Crho%20K_v%20%5Cint_0%5E%7B%5Cinfty%7D%20L%5E2%20G%20n(L)%20d%20L%20%5C%5C%0A%26%20%3D%5Cfrac%7B1%7D%7B2%7D%20%5Crho%20K_a%20%5Cint_0%5E%7B%5Cinfty%7D%20L%5E2%20G%20n(L)%20d%20L%0A%5Cend%7Baligned%7D

对于离散方法,由于生长导致的传质速率可以写成:

%5Cbegin%7Baligned%7D%0A%5Cdot%7Bm%7D%20%26%20%3D%5Crho%20%5Cint_0%5E%7B%5Cinfty%7D%20G_v%20n(L)%20d%20L%20%5C%5C%0A%26%20%3D%5Crho%20%5Cint_0%5E%7B%5Cinfty%7D%20G_v%20n(V)%20d%20V%20%5C%5C%0A%26%20%3D%5Crho%20%5Csum_i%20G_%7BV%2C%20i%7D%20N_i%0A%5Cend%7Baligned%7D

如果总传质过程中包含成核速率,则传质过程为:

%5Cdot%7Bm%7D%3D%5Crho%20V_0%20%5Cdot%7Bn%7D_0%2B%5Csum_i%20%5Crho%20G_%7BV%2C%20i%7D%20N_i

重点:对于离散方法,总体平衡方程的源必须和为总传质速率。要访问源代码,可以使用宏C_PB_DISCI_PS (cell, thread, i)

对于SMM,只有一个与尺寸无关的增长率是可用的。因此,传质速率可以写成:

%5Cdot%7Bm%7D%3D%5Cfrac%7B1%7D%7B2%7D%20%5Crho%20K_a%20G%20m_2

对于QMOM,传质速率可以写成

%5Cdot%7Bm%7D%3D%5Cfrac%7B1%7D%7B2%7D%20%5Crho%20K_a%20%5Csum_i%20L_i%5E2%20w_i%20G%5Cleft(L_i%5Cright)

对于SMMQMOM,由于成核的传质可以忽略不计。

重要提示:注意,对于结晶,初级相有多个组分;至少,有溶质和溶剂。要定义多组分多相系统,您需要在激活多相模型后,在物种模型对话框中激活主要阶段的物种传输。建立物种迁移问题的其余步骤与建立单一阶段的物种迁移问题相同。多相反应定义为: liquid(solvent)j→crystal(solute)

当激活种群平衡模型时,可以自动完成非反应物质(如沸腾)和异相反应(如结晶)相之间的质量传递,而不是挂钩 UDF。

由于非反应物质的成核和生长现象,对于主和次相之间的简单单向传质,配置以下设置:

  1. 在Multiphase Model对话框中,打开Phase Interaction > Heat, Mass, Reactions > Mass
  2. 指定案例Number of Mass Transfer Mechanisms
  3. 对于每个机制,在 From Phase 下指定源材料的阶段,在 To Phase 下指定目标材料阶段的阶段。
  4. 从机制下拉列表中,选择一种机制:
  • none 如果您不想在相之间进行任何传质
  • constant-rate 对于固定的、用户指定的速率
  • user-defined 如果你挂载了一个描述传质机制的 UDF
  • population-balance 用于质量转移的自动化方法,不涉及 UDF。由种群平衡核计算的成核和增长率用于传质。
注:对于非均匀离散种群平衡模型,如果存在一个以上的次级相,则可以选择种群平衡作为溶剂相(例如结晶)和非均匀离散群体平衡模型下定义的每个溶质相之间的传质机制。

 5.点击apply

对于异相反应,请配置以下设置:

  1. 对于主相,激活Species Transport模型
  2. 在Multiphase Model对话框中,打开Phase Interaction > Heat, Mass, Reactions > Reactions
  3. 在“Reactions”选项卡下,指定反应物和产物的化学计量。
  4. 选择“population-balance”作为 Reaction Rate Function.
  5. 点击OK保存

这种方法或使用DEFINE_HET_RXN_RATE宏中描述的UDF都会产生相同的结果。

对于涉及成核和生长的非均匀离散平衡模型,您可以选择 population-balance 作为已设置的每个非均相反应的反应速率函数。

2.理论

2.1 种群平衡方程(PBE)

假设 Φ是粒子体积,数密度函数的输运方程为:

%5Cbegin%7Baligned%7D%0A%26%20%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%20t%7D%5Bn(V%2C%20t)%5D%2B%5Cnabla%20%5Ccdot%5B%5Cvec%7Bu%7D%20n(V%2C%20t)%5D%2B%5Cunderbrace%7B%5Cnabla_v%20%5Ccdot%5Cleft%5BG_v%20n(V%2C%20t)%5Cright%5D%7D_%7B%5Ctext%20%7BGrowth%20term%20%7D%7D%3D%5Cunderbrace%7B%5Cfrac%7B1%7D%7B2%7D%20%5Cint_0%5EV%20a%5Cleft(V-V%5E%7B%5Cprime%7D%2C%20V%5E%7B%5Cprime%7D%5Cright)%20n%5Cleft(V-V%5E%7B%5Cprime%7D%2C%20t%5Cright)%20n%5Cleft(V%5E%7B%5Cprime%7D%20t%5Cright)%20d%20V%5E%7B%5Cprime%7D%7D_%7B%5Ctext%20%7BBirth%20due%20to%20Aggregation%20%7D%7D%20%5C%5C%0A%26%20%5Cunderbrace%7B-%5Cint_0%5E%7B%5Cinfty%7D%20a%5Cleft(V%2C%20V%5E%7B%5Cprime%7D%5Cright)%20n(V%2C%20t)%20n%5Cleft(V%5E%7B%5Cprime%7D%2C%20t%5Cright)%20d%20V%5E%7B%5Cprime%7D%7D_%7B%5Ctext%20%7BDeath%20due%20to%20Aggregation%20%7D%7D%20%5C%5C%0A%26%20%5Cunderbrace%7B%2B%5Cint_%7B%5COmega_r%7D%20p%20g%5Cleft(V%5E%7B%5Cprime%7D%5Cright)%20%5Cbeta%5Cleft(V%20%5Cmid%20V%5E%7B%5Cprime%7D%5Cright)%20n%5Cleft(V%5E%7B%5Cprime%7D%20t%5Cright)%20d%20V%5E%7B%5Cprime%7D%7D_%7B%5Ctext%20%7BBirth%20due%20to%20Breakgge%20%7D%7D%20%5C%5C%0A%26%20-%5Cunderbrace%7Bg(V)%20n(V%2C%20t)%7D_%7B%5Ctext%20%7BDeath%20due%20to%20Breakage%20%7D%7D%20%5C%5C%0A%26%0A%5Cend%7Baligned%7D

边界和初始条件由下式给出:

n(V%2C%20t%3D0)%3Dn_v%20%3B%20%5Cquad%20n(V%3D0%2C%20t)%20G_V%3D%5Cdot%7Bn%7D_0

其中 n0是以particles/m³-s 为单位的成核率。

2.2 颗粒生长

基于颗粒体积的生长速率,Gv,(m³/s)定义为:

G_V%3D%5Cfrac%7B%5Cpartial%20V%7D%7B%5Cpartial%20t%7D

基于颗粒直径(或长度)的生长速率定义为:

G%3D%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20t%7D

单个粒子的体积V定义为KvL³,因此GvG之间的关系为:

G_v%3D3%20K_v%20L%5E2%20G

单个粒子的表面积A定义为KaL²。因此,对于立方体或球体Ka=6Kv

重要:

  • 在颗粒溶解的情况下,由于颗粒体积减少,颗粒体积增长率为负。
  • 在密度(压力)变化导致气泡膨胀的情况下,数密度函数的输运方程中的增长率Gv为气泡体积膨胀率。
  • 考虑到体积V和密度ρ的气泡的质量m,

%5Cbegin%7Baligned%7D%0Am%20%26%20%3D%5Crho%20V%20%5C%5C%0A%5Cfrac%7Bd%20m%7D%7Bd%20t%7D%20%26%20%3D%5Crho%20%5Cfrac%7Bd%20V%7D%7Bd%20t%7D%2BV%20%5Cfrac%7Bd%20%5Crho%7D%7Bd%20t%7D%0A%5Cend%7Baligned%7D

对于恒定m(无传质的气泡膨胀),dm/dt=0,因此

G_V%3D-%5Cfrac%7BV%7D%7B%5Crho%7D%20%5Cfrac%7Bd%20%5Crho%7D%7Bd%20t%7D

米尔斯的论文中有更详细的解释[1]

参考文献:

[1]M. Millies, D. Mewes. “Interfacial area density in bubbly flow”. Chemical Engineering and Processing. 38. 4-6. 307–319. 1999.  


免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

相关文章
技术文档
QR Code
微信扫一扫,欢迎咨询~
customer

online

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 board-phone 155-2731-8020
close1
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

姓名不为空

姓名不为空
手机不正确

手机不正确

手机不正确
公司不为空

公司不为空

公司不为空