文档介绍：相变传热与流体流动数值分析作业2
相变传热与流体流动数值分析
作业2
学院(系):能源与动力学院专业:能源与环境工程
学生姓名:
学号:
指导教师:
完成日期:
大连理工大学
Dalian University of Technology
Finite volume method for two-dimensional diffusion problem (Soluted with ADI and TDMA)
Subject:
The plate size are L=1m,H=1m,and the heat transfer coefficient k equals 1000w/m/k.
Solution:
// 王佳琪- : 定义控制台应用程序的入口点。
//
#include<>
#include<iostream>
#include<>
#include<>
#define N 5
#define M 5
#define pi
#define A1
#define
A2
usingnamespace std;
double T[N][M],an[N][M],as[N][M],aw[N][M],ae[N][M],b[N][M],f1[N][M];
/*------定义追赶法数组,from south to north------*/
void tdma1(double f2[N][M],int p)
{int i;
double y[N][M],l[N][M],u[N][M];
for(i=N-2;i>=0;i--)
{
u[N-1][p]=b[N-1][p];
l[N-1][p]=0;
l[i][p]=as[i][p]/u[i+1][p];
u[i][p]=b[i][p]-l[i][p]*an[i+1][p];
}
y[N-1][p]=f2[N-1][p];
for(i=N-2;i>=0;i--)
y[i][p]=f2[i][p]-l[i][p]*y[i+1][p];
T[0][p]=y[0][p]/u[0][p];
for(i=1;i<N;i++)
T[i][p]=(y[i][p]-an[i][p]*T[i-1][p])/u[i][p];
}
/*-------定义追赶法数组,from west to east-------*/
void tdma2( double f2[N][M],int p)
{int j;
double y[N][M],l[N][M],u[N][M];
for(j=1;j<M;j++)
{
u[p][0]=b[p][0];
l[p][0]=0;
l[p][j]=aw[p][j]/u[p][j-1];
u[p][j]=b[p][j]-l[p][j]*ae[p][j-1];
}
y[p][0]=f2[p][0];
for(j=1;j