求解流体力学中有限体积（FVM）法的matlab代码

求解流体力学中有限体积（FVM）法的matlab代码

function datagen
clc;
clear all;
format long;

h=0.01;    %%每隔h秒传感器测一次温度，滤波器工作一次
nc=1000;    %%温度传感器测量次数
time=linspace（0.0110.0nc）;
%%%%%%%%%%%%%%%%%初边界条件设置%%%%%%%%%%%%%%%%%%%%%
cq=1.0;%%热流
rcp=1.0;     %%密度和比热的乘积
SL=1.0;    %%传热长度
%%%%%%%%%%%%%%%%%初边界条件设置完毕%%%%%%%%%%%%%%%%%

nx=50;   %%控制体积的个数--nx
nt=10;   %%每个滤波步长内积分步数--nt

T0=0.0*ones（1nx）;   %%每个控制体积的初始温度

%%%%%%%%%%%%%%%%%%热传导系数确定%%%%%%%%%%%%%%%%%%%
theta=zeros（13）;
theta（1）=1.0;
theta（2）=0.0e-003;
theta（3）=0.0e-005;
%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%%%

Tout=zeros（1nx）;  %%t+dt时刻的各控制体积的温度分布
TC=zeros（1nc）;   %%各时刻的温度测量值
CK=zeros（1nc）;   %%各刻传系数的记录值
realTC=zeros（1nc）;

M_ran=random（‘norm‘0.01.01nc）;
%%%%%%%%%%%%%%%%%%%%开始温度测量%%%%%%%%%%%%%%%%%%%%
Tin=T0;

for i=1:nc
    [ToutTCK] = fvm（TinthetanxSLnthrcpcq）;       %%[ToutTCK] = feval（‘fvm‘TinthetanxSLnthrcpcq）;    
    display（i）;
    Tin=Tout;
    TC（i）=TCK+0.001*M_ran（i）;
    realTC（i）=TCK;
    CK（i）=theta（1）+theta（2）*TCK+theta（3）*TCK*TCK;
end
%%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%解析解%%%%%%%%%%%%%%%%%%%%%%%
Tjx=zeros（1nc）;
for i=1:nc
    vn=0.0;
    for j=1:1000
        vn=vn+2.0*exp（-j*j*pi*pi*time（i））*cos（j*pi）/（j*j*pi*pi）;
    end
    Tjx（i）=T0（nx）+time（i）-1.0/6.0-vn;
end
%%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%%

figure;
plot（timeTC‘b‘timeTjx‘r‘）;
figure;
plot（timeTC-Tjx‘r‘）;

%%%%%%%%%%%%%%%%%%%%测量值存储%%%%%%%%%%%%%%%%%%%%%%

save TC.dat TC -ascii;
save CK.dat CK -ascii;
save realTC.dat realTC -ascii;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      16002  2014-04-03 09:03  有限体积法data_gen\data_gen\CK.dat

     文件       1418  2009-09-07 09:51  有限体积法data_gen\data_gen\datagen.asv

     文件       1745  2009-09-21 22:13  有限体积法data_gen\data_gen\datagen.m

     文件       3336  2009-09-08 11:29  有限体积法data_gen\data_gen\fvm.m

     文件      16002  2014-04-03 09:03  有限体积法data_gen\data_gen\realTC.dat

     文件      16002  2014-04-03 09:03  有限体积法data_gen\data_gen\TC.dat

     目录          0  2014-04-04 10:50  有限体积法data_gen\data_gen

     目录          0  2014-04-04 10:50  有限体积法data_gen

----------- ---------  ---------- -----  ----

                54505                    8