换热器简化模型仿真.docVIP

本次汇报的主要内容为：根据一篇论文中换热器的物理简化模型，和已有的参数数据，对模型进行了重现（应用了欧拉迭代法）。 1.换热器模型简化 （1） 式中—单位长度管壁的金属质量；—金属比热；—管外放热介质对单位管长的管壁金属在单位时间内的放热量；—管壁金属向管内介质在单位管长和单位时间内的放热量。 管内介质在单位管长和单位时间内的放热量可表示： （2） 式中—管子内径；—单位长度管子的内表面积；—分别表示金属和管内介质的温度；—放热系数，可表示为（3），式中是常数，是管内介质的流量。 管外介质在单位管长和单位时间内的放热量可表示： （4） 式中—介质向单位管长管壁的对流放热量；—介质向单位管长管壁的辐射放热量。且有 （5） （6） 式中、—分别为对流和辐射放热系数；—管外介质温度；—管子外径；—单位长度管子外表面积。把方程（2）-（6）代入方程（1）中 （7） 令，，，，（7）式可以： （8） 令（时间常数）,，（8）式 （9） （9）式存在解析解 (10) 解析解中有指数的形式，占用机时较多。仿真计算中大多不用解析算法，而是从（9）式出发，构造差分格式，用一定精度的近似解代替解析解，达到实时仿真的目的。 2. 差分格式算法（前进欧拉法） 2.1基本原理方法 前进欧拉法是基于向前微分近似变形为 其中，以n循环得到 2.2具体构造的欧拉法迭代公式 %显示欧拉法程序代码：（p=0.02;q=10;h=1;T0=20;） %T = dsolve(DT +p*T-q=0,T(0) =T0,t) %T =(q - (q - T0*p)/exp(p*t))/p %微分方程的函数文件： function z=mTf4(x,T) p=0.02; q=10; z=q-p*T; %准确解的函数文件（解析解）： function fz=fz(x) p=0.02; q=10; T0=20; fz=q/p-(1/p)*(q-T0*p)./exp(p*x); %欧拉格式算法： %初始值和步长假定 p=0.02; q=10; h=1; T0=20; x=zeros(1,01); T=zeros(1,301); x(1)=0; T(1)=20; for n=1:300 x(n+1)=x(n)+h; T(n+1)=T(n)+h*feval(@mTf4,x(n),T(n)); end E=[x,T]; T1=fz(x); plot(x,T,:og,x,T1,--b); title(Euler格式与准确解比较图); 从图中可以看出欧拉迭代法的近似解和微分方程的精确解的误差很小，从而证明了欧拉迭代算法满足了题目要求。 3.存在问题 1）迭代的差分格式算法替代微分方程的精确解，仅仅是因为计算时间问题？，还是有其他原因？ 2） 原论文中采取了其他迭代算法（改进的欧拉法、梯形法和四阶龙格库塔法），这里只是选取了一种算法进行重现，也没有进行算法稳定性的对比分析。 汇报人： 杨鹏志