合作共享

clc

clear all

close all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%气泡实验物性条件,输入参数

stf=0.072;%表面张力

ug=1;%曝气速度

d0=0.001;%孔径

rhog=1;%气相密度

rhol=998;%液相密度

miul=0.001;%液相粘度

g0=9.81;%重力加速度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Q0=pi.*(d0.^2).*ug./4.;%孔口气体流量

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%气泡直径范围与步长

dp=[1.E-5:1.E-5:5.E-2]';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

V=pi./6.*dp.^3;

t=V./Q0;%生成气泡所需时间

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 颈长

lcc=d0;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

z=dp.*0.5+lcc;%气泡中心高度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 虚拟质量力系数cv

cv=11./16;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%气泡上升速度

ub=z./t;%平均速度假设

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%无量纲数

wei=rhol.*d0.*ug.^2./stf;

mo=(rhol-rhog).*g0.*miul.^4./rhol.^2./stf.^3;

%液速

ueq=Q0./(pi.*dp.^2./4);%r为气泡半径

K=(1.5897+0.0971.*wei+0.00000282.*wei.^2+27590.*mo)./(1+0.301.*wei+155403.*mo-6105.*mo.^2);%%%拟合液速模型

ul=K.*ueq;%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%相对速度

us=ub-ul;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%浮力

Fb=(rhol-rhog).*g0.*(pi./6.*dp.^3);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%气体动量力

Fm=rhog.*Q0.*ug;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%表面张力

Fs=pi.*d0.*stf;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%惯性力

Fi=(rhog+rhol.*cv).*pi.*dp.^3./6.*us./t;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%曳力

eo=g0.*dp.^2.*(rhol-rhog)./stf;

redp=dp.*us.*rhol./miul;%雷诺数

cd=max(min(24./redp.*(1+0.15.*redp.^0.687),72./redp),8./3.*eo./(eo+4));%曳力系数

Fd=pi.*dp.^2./4.*cd.*rhol.*us.^2./2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%迭代筛选

po=(Fb+Fm-Fs-Fi-Fd).^2.;%%%%%%Fb=Fs

pre=find(po==min(po));

dpred=dp(pre)%预测气泡直径

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%