MATLAB拟合多孔介质参数-fluent计算惯性阻力系数和粘性阻力系数.docxVIP

MATLAB拟合多孔介质参数，fluent计算惯性阻力系数和粘性阻力系数

包涵：1，MATLAB调取NISTrefprop物性参数；

2，MATLAB拟合多孔介质实验参数（二次函数，参数项目为0）

3，按fluent计算公式计算惯性阻力系数和粘性阻力系数

步骤：

一、根据需要安装refprop物性参数和MATLAB调用连接（物性参数也可手动输入不用安装refprop）

1)下载refprop并安装，

2)下载如下4个文件到新建的MATLAB程序保持目录下:

二、建立实验参数数据表excl文件。

实验名称与matlab调用名称一致，实验数据自己根据实验数据填写

三、编写MATLAB程序

没有安装refprop

没有安装refprop物性参数可直接激活输入值

结果显示：

调取的空气物性

拟合的二次函数值：

Fluent软件多孔介质中惯性阻力系数C2和粘性阻力系数D

Fluent软件中输入

程序：

filename1=D:\porous\?à?×?é?êêμ?é2?êy.xlsx

[F1]=importdata(filename1)%êy?Yμ?è????ó?D%

V=F1.data(:,1);%?áè??ù?èm/s%

P=F1.data(:,2);%?áè??1?μPa%

L=F1.data(1,3);%?áè??à?×?é?ê3¤?è%

symsx

f=fittype(a*x*x+b*x,independent,x,coefficients,{a,b});

cfun=fit(V,P,f)%??ê??ao?oˉêy￡?êy?Y±?D??aáD?òá?D?ê?

xi=0:1:40;

yi=cfun(xi);

N=polyfit(xi,yi,2)

figure

plot(V,P,r*,xi,yi,b-);%?ao??ú??oíêμ?é2?êyμ?

title(?ao?oˉêyí?D?);

%figure

%plot(x,sqrt(y-yi(1:1:18,:).^2/18),r*);

%title(±ê×??ó2?í?);

[rho,u,Cp,lamda]=refpropm(DVCL,T,273.15+20,P,101.15,air.mix)%μ÷è?NISTèí?trefprop2?á???D?

%20??￡?±ê????2?á???D?

%rho=input(ê?è??ü?èkg/m3:)

%N=[3.307,70.08]

%L=input(ê?è??à?×??óò3¤?èm:)

%u=input(ê?è??ˉá|?3?è?μêyPa*s:)

C2=N(1)/(rho*L/2)%1?D?×èá|?μêyC2%

%%C2=solve(C2*rho*L/2-N(1)=0)%1?D?×èá|?μêyC2%

D=N(2)/u/L%?3D?×èá|?μêyViscousresistancecoefficient

%%D=solve(D*u*L-N(2)=0)%?3D?×èá|?μêy