MATLAB擬合多孔介質(zhì)參數(shù)-fluent計(jì)算慣性阻力系數(shù)和粘性阻力系數(shù)

MATLAB擬合多孔介質(zhì)參數(shù)，fluent計(jì)算慣性阻力系數(shù)和粘性阻力系數(shù)包涵：1，MATLAB調(diào)取NISTrefprop物性參數(shù)；2，MATLAB擬合多孔介質(zhì)實(shí)驗(yàn)參數(shù)（二次函數(shù)，參數(shù)項(xiàng)目為0）3，按fluent計(jì)算公式計(jì)算慣性阻力系數(shù)和粘性阻力系數(shù)步驟：一、根據(jù)需要安裝refprop物性參數(shù)和MATLAB調(diào)用連接（物性參數(shù)也可手動(dòng)輸入不用安裝refprop）1)下載refprop并安裝，2)下載如下4個(gè)文件到新建的MATLAB程序保持目錄下:二、建立實(shí)驗(yàn)參數(shù)數(shù)據(jù)表excl文件。實(shí)驗(yàn)名稱與matlab調(diào)用名稱一致，實(shí)驗(yàn)數(shù)據(jù)自己根據(jù)實(shí)驗(yàn)數(shù)據(jù)填寫(xiě)三、編寫(xiě)MATLAB程序沒(méi)有安裝refprop沒(méi)有安裝refprop物性參數(shù)可直接激活輸入值結(jié)果顯示：調(diào)取的空氣物性擬合的二次函數(shù)值：Fluent軟件多孔介質(zhì)中慣性阻力系數(shù)C2和粘性阻力系數(shù)DFluent軟件中輸入程序：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¤?è%symsxf=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)figureplot(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?×èá|?μêy