几个月前,在上微反应器课的时候,看到老师讲微通道反应器的二氧化碳吸收
突然灵光一闪,感觉这个过程可以用MATLAB复现啊
于是就写了一个代码,大家感兴趣可以看看
主要是最近好几个项目写不下去了闲的无聊
具体模拟的结果和动画效果如下
这不是我的主要方向,所以感觉有很多不足的地方,如果代码有什么不对,请多多帮忙
%这是一个微反应器模拟的主函数
%气体系统为氮气和二氧化碳
%溶液为乙二醇胺溶液
clc
clear
%尺寸输入这里采用了m为单位
L1=4*10^-4;
L2=L1;
L0=0.03;
%为了直观方便打印的时候放大了10^6倍
x_equipmennt=10^6*[-8*L1,-8*L1,2*L1,2*L1,3*L1,3*L1,3*L1+L0,3*L1+L0,-8*L1];
y_equipmennt=10^6*[0,L1,L1,3*L1,3*L1,L1,L1,0,0];
%这是二氧化碳体积分数的定义
fain2=0.80;
faico2=1-fain2;
%这是气液相流量单位ml/h
QL=10;
QG=50;
%这是黏度 单位mP·s
nd=1.633*10^-3;
%这是表面张力 单位mN/m
zl=(25.28*10^-3)*0.75;
uc=QG*10^-6/3600/L1^2;
Ca=nd*uc