fdm之一维动态热传导两层不同导热系数K
来源:互联网 发布:js控制标签隐藏 编辑:程序博客网 时间:2024/04/28 21:05
function heat_transit_1d_varK
q=0.065;
Z=20;K=20; % conductivity
TT=4;
times=100;
dt=TT/times;
n=200;
delz=Z/(n-1);
z=0:delz:Z;
A=zeros(n,n);
K2=0.3*K;
for ii=2:n-1
if ii<n/4
alpha=K*dt/(delz*delz);
A(ii,ii)=1+2*alpha;
A(ii,ii-1)=-alpha;
A(ii,ii+1)=-alpha;
elseif ii>n/4
alpha3=K2*dt/((delz*delz));
A(ii,ii)=1+2*alpha3;
A(ii,ii-1)=-alpha3;
A(ii,ii+1)=-alpha3;
elseif ii==n/4
alpha2=(K+K2)*dt/((delz*delz)*2);
% alpha2=2*K*K2*dt/((delz*delz)*(K+K2));
A(ii,ii)=1+2*alpha2;
A(ii,ii-1)=-K*dt/(delz*delz);
A(ii,ii+1)=-K2*dt/(delz*delz);
end
end
% a=1+ones(1,n)*2*alpha;
% b=-ones(1,n-1)*alpha;
% A=diag(a,0)+diag(b,-1)+diag(b,1);
A(1,1)=1;
A(1,2)=0;
A(n,n-1)=0;
A(n,n)=1;
upperBC = 20; % set temp BCs, in degrees
lowerBC = 20;
B=zeros(n,1);
B(1:n) = 10;
%line(T,z);
axis ij;
xlabel('temperature');
ylabel('Depth');
title('Geothermal Temps');
for ii=1:dt:TT
B(1) = 10*(1+sin(2*pi*ii));
%BC
T=inv(A)*B;
h=plot(B,z);
axis ij;
hold on;
pause(0.1);
% delete(h);
B=T;
end
0 0
- fdm之一维动态热传导两层不同导热系数K
- fdm之一维静态热传导不同导热系数K
- fdm之一维动态热传导
- FDM之二维静态热传导--含有不同传导系数K
- fdm之一维静态热传导
- fdm之二维静态热传导
- FEM之一维静态热传导
- 泡沫金属导热系数测试
- FDM
- 导热系数/热阻测试--美信检测
- fdm之一维弹性波的右边无反射边界处理
- 两层交换机不同vlan通讯
- 一维热传导问题与C++描述
- 皮尔逊系数皮尔逊系数---K近算法之皮尔逊系数
- 对齐系数 #pragma pack(k)
- 动态规划 计算二项式系数
- 动态规划之二项式系数
- 动态规划 计算二项式系数
- iOS #include与#import的区别
- UML基本架构建模--通用机制的通用建模技术(三)
- 《STL源码解析》读书笔记之序列式容器(2)
- [小技巧] Excel 里搜索所有的 Sheet
- UML基本架构建模--给通用机制建模时的注意事项
- fdm之一维动态热传导两层不同导热系数K
- 关于内存对齐的那些事
- ios Crash Log 分析汇总
- jquery拉拽div
- CSS 元素选择器的基本概念(仅理论知识总结)
- android声音播放
- 新浪微博客户端开发
- 触发事件传递
- Eclipse的快捷键