【物理应用】Matlab 二维对流扩散温度场

function A=secdim2% 本程序可实现两种可选边界条件的对流扩散实时动态模拟
% 采用乘方格式
% 初始温度场为10，速度可选
% 时间差分采用全隐式格式
clc
dt=0.001;
bc=input('南北边界条件：1-第一类，2-第二类:  ');
q=0;
ddt=input('迭代次数： ');
L1=0.02;r=1000;k=0.017;sp=0;L2=0.01;
u=input('横向速度： ');if bc==2v=0;
elsev=input('纵向速度 ： ');
end
n=100;m=100;
dx=L1/n;dy=L2/m;F1=r*u*dy;F2=r*v*dx;
ap0=r*dx*dy/dt;
D1=k*dy/dx;pe1=F1/D1;D2=k*dx/dy;pe2=F2/D2;tp0=10*ones(n*m,1);%初场温度tw=10;te=10; %边界条件A=sparse(n*m);
b=ap0*tp0;
%二维影响系数采用乘方格式
ae=D1*max(0,(1-0.1*abs(pe1^5)))+max(-F1,0);
aw=D1*max(0,(1-0.1*abs(pe1^5)))+max(0,F1);
as=D2*max(0,(1-0.1*abs(pe2^5)))+max(0,F2);
an=D2*max(0,(1-0.1*abs(pe2^5)))+max(0,-F2);ap=aw+ae+as+an+ap0-sp;spw=-1*2*k*dy/dx;
spe=-1*2*k*dy/dx;
%边界条件源项处理
if bc==2q=input('热流为: ');spn=0;sps=0;tn=0;ts=0;te=0;spe=0;
elsetn=input('北侧温度： ');ts=input('南侧温度： ');
spn=-1*2*k*dx/dy;
sps=-1*2*k*dx/dy;
end
sun=zeros(n,1);
sus=zeros(n,1);
sun(5:64)=2*k*tn*dx/dy+q*dx;
sus(5:64)=2*k*ts*dx/dy+q*dx;
suw=2*k*tw*dy/dx;
sue=2*k*te*dy/dx;%系数矩阵，常数向量赋值
for i=1:n*m
%南s边界点赋值  if(i<n&&i>1)A(i,i)=ae+an+aw+ap0-sps;A(i,i+1)=-ae; A(i,i-1)=-aw;A(i,i+n)=-an;b(i)=b(i)+sus(i);
%内节点inner赋值else if(i>n&&rem(i,n)~=1&&rem(i,n)~=0&&i<n*(m-1)+1)   A(i,i)=ap;A(i,i+n)=-an; A(i,i-n)=-as; A(i,i+1)=-ae; A(i,i-1)=-aw;
%北n边界点赋值else if(i>n*m-n+1&&i<n*m)          A(i,i)=ae+aw+as+ap0-spn;A(i,i+1)=-ae; A(i,i-1)=-aw; A(i,i-n)=-as;b(i)=b(i)+sun(i-n*m+n);%西w边界点赋值    else if(i>n&&i<n*m-n+1&&rem(i,n)==1)A(i,i)=ae+an+as+ap0-spw;A(i,i+n)=-an; A(i,i-n)=-as;  A(i,i+1)=-ae;b(i)=b(i)+suw;%东e边界点赋值else if(i>n&&i<n*m&&rem(i,n)==0)A(i,i)=aw+an+as+ap0-spe;A(i,i+n)=-an; A(i,i-n)=-as;  A(i,i-1)=-aw;b(i)=b(i)+sue;elseendendendendend
end
%四个角点赋值%西南角点A(1,1)=1; %东南角点  A(n,n)=1;  %西北角点    A(n*m-n+1,n*m-n+1)=1;%东北角点A(n*m,n*m)=1;%
end

【物理应用】Matlab 二维对流扩散温度场相关推荐

二维对流方程matlab求解,二维对流扩散方程的有限元计算方法
冯立伟+张成+屈福志 " " " 摘要:针对二维对流扩散方程边值问题,采用三角形剖分,使用二维线性有限元进行计算分析.采用matlab编写了计算程序,使用算例进行了数值实 ...
[计算流体力学][Matlab] 使用 A,B,C 格式与蛙跳格式求解二维对流问题
1. 题目 2. 转述原题目要求可以简化为: 对于二维对流方程: ∂u/∂t+∂u/∂x+∂u/∂y=0 u(x,y,0)={█(1,when-4≤x≤4,-4≤y≤4@0,other)┤ 取计算范 ...
MATLAB二维绘图(二)向图中添加标题，坐标轴，图标和文字信息
MATLAB二维绘图(二)向图中添加标题,坐标轴,图标和文字信息 1.添加标题.图例.x轴信息和y轴信息,示例: %% 添加标题 clear; clc; close all; x = 0:0.1:2* ...
Matlab 二维网格图pcolor和imagesc区别
Matlab 二维网格图pcolor和imagesc区别在速度或者衰减层析成像反演之后会得到地下介质每个网格点的速度/衰减值,在进行呈现的时候可以使用maltab的pcolor和imagesc进行二 ...
matlab 二维矩阵变成一维矩阵
matlab 二维矩阵变成一维矩阵 1.一维变二维: https://blog.csdn.net/qq_40584593/article/details/90691276 reshape 2.a(:) ...
Matlab二维图形绘制与图形处理
Matlab二维图形绘制与图形处理一.二维图形绘制 1.极坐标图 2.散点图 3. 平面等值线图二.图形处理 1.添加格栅,图例和标注 2.定制坐标 3.在之前基础上继续作图 4.新建图形置于当 ...
matlab二维绘图部分
matlab二维绘图部分 X,Y是向量,分别表示点集的横坐标和纵坐标 PLOT(X,Y,S) 符号函数(显函数.隐函数和参数方程)画图 (1) ezplot ezplot('f(x)',[a,b]) ...
matlab二维三维图形绘制和坐标轴范围设置
matlab二维绘图一．二维图形(Two dimensional plotting) 1. 基本绘图函数(Basic plotting function):Plot, semilogx, se ...
MATLAB二维图形绘制
MATLAB二维图形绘制数据点标记数据点颜色线型曲线图一.plot函数 (一)最简单的plot函数调用格式: plot(x) (1)plot函数的参数 X 为普通向量, (2)plot函数 ...

【物理应用】Matlab 二维对流扩散温度场

【物理应用】Matlab 二维对流扩散温度场相关推荐

最新文章

热门文章