相变传热与流体流动数值分析作业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 王佳琪 作业 cpp 定义控制台应用程序的入口点 include include include include define N 5 define M 5 define pi 3 1415926 define A1 0 1 define A2 0 1 using namespace 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 0 j T p j y p j ae p j T p j 1 u p j void main double Sp N M f2 N M k x y H L Tn M Ts Tw Te int i j h 边界条件 k 1000 H 1 L 1 y H N x L M Ts 0 Tw 0 Te 0 for j 0 j M j Tn j 100 sin 2 j 1 pi 2 M cout Tn j Tn j 离散追赶法函数系数 for i 0 i N i for j 0 j M j T i j 0 if i 0 an i j 0 else an i j k A1 y if i N 1 as i j 0 else as i j k A1 y if j 0 aw i j 0 else aw i j k A2 x if j M 1 ae i j 0 else ae i j k A2 x if i 0 else if i N 1 else if i 0 else if i N 1 else if i 0 else if i N 1 else if j 0 else if j M 1 else Sp i j 0 b i j Sp i j aw i j ae i j as i j an i j 初始化边界热量传递 for i 0 i N i for j 0 j M j if i 0 else if i 0 else if i N 1 else if i N 1 else if i 0 else if i N 1 else if i 0 else if i 0 else f1 i j 0 for h 1 h 100 h 不考虑东西方向 只南北循环计算 for j 0 j M j for i 0 i N i if j 0 f2 i j f1 i j ae i j T i j 1 else if j M 1 f2 i j f1 i j aw i j T i j 1 else f2 i j f1 i j aw i j T i j 1 ae i j T i j 1 tdma1 f2 j 不考虑南北方向 只东西循环计算 for i 0 i N i for j 0 j M j if i 0 f2 i j f1 i j as i j T i 1 j else if i N 1 f2 i j f1 i j an i j T i 1 j else f2 i j f1 i j an i j T i 1 j as i j T i 1 j tdma2 f2 i cout h endl 输出结果 for i 0 i N i for j 0 j M j cout T i j T i j cout endl The result Discussion From the figure we can see the temperature distribution is the same as the physical analysis Because of the grid number limitation we can impro