转自:微信公众号 cfd之道

用来自记!!!!!
转载链接:https://mp.weixin.qq.com/s/nhDoHNsTwenXZNBYkDVoVg

在液体喷雾过程中,常利用Rosin-Rammler方法来描述液滴粒径分布。

本文内容来自Fluent UserGuide 25.3.14。

Rosin-Rammler方法利用以下方式描述粒径与质量分数之间的函数关系:

齐总\bar{d}为平均粒径(Mean Diameter),n为尺寸分布指数(Spread Parameter)。

通过将粒径分布数据拟合到Rosin-Rammler方程中,可以很容易地定义粒度分布。在这种方法中,完整的粒径范围被划分为一组离散的粒度范围。例如,假设粒径数据服从以下分布:

粒径范围(微米) 质量分数
0~70 0.05
70~100 0.1
100~120 0.35
120~150 0.3
150~180 0.15
180~200 0.05

首先需要处理表中的数据,以累积质量分数的形式显示:

粒径(微米) 超过粒径的质量分数
70 0.95
100 0.85
120 0.50
150 0.20
180 0.05
200 0

将表显示成散点图,如下图所示。

利用Rosin-rammler函数拟合上面的数据。这种非线性估计极为麻烦,曲线拟合常常失败。可以采用python对上表中的数据进行拟合。

import matplotlib.pyplot as plt
from scipy.optimize import curve_fit
import numpy as np
x = np.array([70,100,120,150,180,200])
y = np.array([0.95,0.85,0.5,0.2,0.05,0])# 防止对0取对数,故去掉最后一个元素
xx = np.array([70,100,120,150,180])
yp = np.array([0.95,0.85,0.5,0.2,0.05])yy= np.log(yp)'''
指定的公式,两边取对数然后拟合
'''
def func(x,a,b):return -np.power(x/a,b)
popt, pcov = curve_fit(func, xx, yy)#函数拟合# 绘制图形
xx = np.linspace(70,200,num=50)
yvals=np.exp(func(xx,popt[0],popt[1]))
plot1=plt.plot(x, y, 'o',c='r',label='original values')
plot2=plt.plot(xx, yvals, 'b',linewidth=2,label='$fit: Y_d = e^{-(\\frac{d}{%5.5f})^{%5.5f}}$' % tuple(popt))
plt.xlabel('diameter(μm)')
plt.ylabel('mass fraction')
plt.legend(loc=1)
plt.show()

拟合结果如图所示:

可看到拟合的系数a=134.247,b=3.7794。这里为了防止对0取对数而丢失了一个点的信息。

Fluent中采用另外一种估算方案。

先估算平均粒径。当粒径为平均粒径时,此时质量分数:

此时线性插值得到平均粒径:

可得到平均粒径d=133.2μm。

有了平均粒径,即可代入公式:

将表中的粒径d代入公式中,求解得到多个n,再计算其平均值即可得到分布指数。

粒径 质量分数 n
70 0.95 4.682585
100 0.85 6.5445
120 0.5 3.845475
150 0.2 3.722698
180 0.05 3.53755
200 0
平均值 4.466562

可得到平均分布指数n=4.466562.

    |

| | 平均值 | 4.466562 |

可得到平均分布指数n=4.466562.

[外链图片转存中…(img-SGpXEVvs-1601943704493)]

看图形拟合得还不错。当然如果想要硬生生的拟合,也并不是不可以,不过搞起来麻烦一点罢了,不管从哪个角度来讲,非线性拟合的复杂程度总是要大于代数计算。

Rosin-Rammler液滴粒径分布相关推荐

  1. 实验型喷雾干燥机粒径分布测量方法

    实验型喷雾干燥机粒径分布测量方法 有许多方法可以测量小型实验型喷雾干燥机器的液滴尺寸和分布. 根据其测量原理,可以分为三类. 实验喷雾干燥机的粒度数据对于控制干粉产品的质量是必不可少的. 但是,粒度并 ...

  2. 绿色荧光ps微球/聚苯乙烯彩色荧光微球具有荧光强度高,性能稳定,粒径分布窄

    绿色荧光ps微球/聚苯乙烯彩色荧光微球具有荧光强度高,性能稳定,粒径分布窄 绿色荧光微球 聚苯乙烯荧光微球具有荧光强度高,性能稳定,粒径分布窄等特性, 可广泛用于医疗诊断.血流测定.示踪.体内成像,以 ...

  3. 【Android探索】基于OpenCV的微液滴粒径分析APP

    前言:这个App是之前<数字图像处理>课程的一次课程设计中的产物,现在整理一下记录下来,里面涉及到了比较多的控件以及拓展包,功能不是很丰富但是也比较算齐全,其中使用的技术原理包括在安卓上使 ...

  4. cv曲线面积的意义_关于粒径及粒度分布的定义及其在粉末涂料生产中的指导意义...

    摘要:粒径和粒径分布式粉末涂料的重要参数,文章介绍了粉体材料的粒径和粒径分布定义,表示方法,以及在粉末涂料生产中的指导意义. 1.粒径的定义 当被测颗粒的某种物理特性或物理行为与某一直径的同质球体(或 ...

  5. 含稀土配合物粒径小/发色性能好的聚苯乙烯荧光微球/交联羧基聚苯乙烯微球/单分散荧光微球/红色荧光微球

    含稀土配合物粒径小/发色性能好的聚苯乙烯荧光微球/交联羧基聚苯乙烯微球/单分散荧光微球/红色荧光微球 含稀土配合物粒径小/发色性能好的聚苯乙烯荧光微球/交联羧基聚苯乙烯微球/单分散荧光微球/红色荧光微 ...

  6. 磺酸基修饰的聚苯乙烯核壳结构微球,磁核粒径50-80nm,聚苯乙烯粒径30-50nm

    磺酸基修饰的聚苯乙烯核壳结构微球,磁核粒径50-80nm,聚苯乙烯粒径30-50nm 磺酸基修饰磁性PS微球(20μm) 产品名称:磺酸基修饰的聚苯乙烯微球 别名:磺酸基标识的聚苯乙烯微球:磺酸基标记 ...

  7. 粒径测量中的D10 D50 D90的含义

    D10.D50.D90的的粒径大小的参数,代表的含义是10%.50%.90%的颗粒尺寸在所测得的尺寸值. D10:颗粒累积分布为10%的粒径,即小于此粒径的颗粒体积含量占全部颗粒的10%. D50:颗 ...

  8. 不同粒径大小的金纳米粒子|球形金纳米粒子|单分散金纳米颗粒Glucose modified Gold Nanoparticles(粒径5nm)

    不同粒径大小的金纳米粒子|球形金纳米粒子|单分散金纳米颗粒Glucose modified Gold Nanoparticles(粒径5nm) 金属纳米颗粒是尺寸在1-100纳米的金属原子聚集体,比光 ...

  9. 粒径2-5nm黄色金纳米颗粒|粒径10-20nm酒红色金纳米颗粒|粒径30-80nm紫红色金纳米颗粒|粒径520-530nm黑色金纳米颗粒

    粒径2-5nm黄色金纳米颗粒|粒径10-20nm酒红色金纳米颗粒|粒径30-80nm紫红色金纳米颗粒|粒径520-530nm黑色金纳米颗粒 金属纳米颗粒是尺寸在1-100纳米的金属原子聚集体,比光的波 ...

最新文章

  1. mysql insert 数据_MySQL-插入数据(INSERT)
  2. linux回到桌面的命令符_Linux命令行环境与桌面环境护切换
  3. CodeCombat编程游戏
  4. java swing鼠标事件监听_java swing鼠标监听问题
  5. Python中最重要的知识点:切片的操作,学习python一定要会的操作!
  6. Vue实现仿音乐播放器14-实现搜索页面以及功能
  7. Golang之envconfig的使用
  8. “No operations defined in spec!”一文教你swagger如何扫描多个controller
  9. android实现日志列表,【Android开发】listview+popupwindow实践:日志列表
  10. linux动态可执行文件,Linux中ELF格式 可执行文件+动态链接器 的加载
  11. python 创建二维list,Python创建二维数组(关于list的一个小坑)
  12. 桌面计算机图标变黑块,电脑中的文件夹图标变成黑色的方块的四种解决方法
  13. C语言怎么播放背景音乐
  14. -- 99美金和299美金的开发者证书的区别 --
  15. js-入门(字符串-运算符) html常用命令代码行
  16. Cao Le的 Emacs org mode 学习笔记
  17. PBRT学习笔记:在单位圆内部均匀采样
  18. Push to origin/master was rejected异常解决办法
  19. SQlite跨库查询
  20. 怎么在不重启服务器情况下查看管理口-IPMI的ip是多少?

热门文章

  1. ncode2019安装教程_nCode DesignLife|ANSYS 17.0 nCode DesignLife下载 v12.0免费版 附使用教程 - 121下载站...
  2. 逻辑思维-03取水问题
  3. 超级推荐,Mac 端ssh连接工具termius
  4. 我的PhoneGap安装配置经历
  5. java描绘条形码_java生成条形码工具类
  6. 计算机网络学习笔记之数据链路层的介质访问控制
  7. STC15单片机-通过PWM调整灯亮度
  8. 修改Docker镜像仓库为国内163镜像源
  9. 我是一个假的acmer
  10. 优秀的内存规划方法——环形缓冲区(ring buffer)