东南大学传热学二维稳态差分接点计算论文

二维稳态计算实验报告

王## 能源与环境学院03012###一、题目要求

二、各节点离散化代数方程

取研究节点为t(i,j)，其上部节点温度用T表示，下部节点温度用B表示，左侧节点温度用L表示，右侧节点温度用R表示。该题目各节点离散化代数方程可分为四个区域：

区域一：

t(i,j)=

区域二：

t(i,j)=t i-1,j+t i,j-1+t i+1,j+t i,j+1

4

区域三：

t(i,j)=

区域四：

t(i,j)=

三、源程序

Jacobi迭代计算代码（C#）：

for (int n = 0; n < 6; n++)

{

mg0[m, n] = new Button();

mg0[m, n].Text = mg[m, n].Text.ToString ();

}

for (int m = 1; m < 4; m++)

{

int n = 1;

double t = Convert.ToDouble(mg0[m - 1, n].Text.ToString()); double b = Convert.ToDouble(mg0[m + 1, n].Text.ToString()); double l = Convert.ToDouble(mg0[m, n - 1].Text.ToString()); double r = Convert.ToDouble(mg0[m, n + 1].Text.ToString()); double center = (200 + t + 2*r + b) / 24;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

for (int m = 1; m < 4; m++)

{

for (int n = 2; n < 5; n++)

{

double t = Convert.ToDouble(mg0[m - 1,n].Text.ToString()); double b = Convert.ToDouble(mg0[m + 1, n].Text.ToString()); double l = Convert.ToDouble(mg0[m, n - 1].Text.ToString()); double r = Convert.ToDouble(mg0[m, n + 1].Text.ToString()); double center = (l + t + r + b) / 4;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

}

for (int m = 4; m < 5; m++)

{

int n = 1;

double t = Convert.ToDouble(mg0[m - 1, n].Text.ToString());

double l = Convert.ToDouble(mg0[m, n - 1].Text.ToString()); double r = Convert.ToDouble(mg0[m, n + 1].Text.ToString()); double center = (100 + t + r) / 12;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

for (int n = 2; n < 5; n++)

{

int m = 4;

double t = Convert.ToDouble(mg0[m - 1, n].Text.ToString());

double r = Convert.ToDouble(mg0[m, n + 1].Text.ToString());

double center = (2*t + r + l) / 4;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

j = Convert.ToDouble(mg[1, 1].Text.ToString()); for (int m = 0; m < 6; m++)

{

for (int n = 0; n < 6; n++)

{

mg0[m, n] = new Button();

mg0[m, n].Text = mg[m, n].Text.ToString();

}

}

高斯（Gauss）-赛德尔（Seidel）迭代计算代码（C#）：

for (int m = 1; m < 4; m++)

{

int n = 1;

double t = Convert.ToDouble(mg[m - 1, n].Text.ToString());

double b = Convert.ToDouble(mg[m + 1, n].Text.ToString());

double l = Convert.ToDouble(mg[m, n - 1].Text.ToString());

double r = Convert.ToDouble(mg[m, n + 1].Text.ToString());

double center = (200 + t + 2*r + b) / 24;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

for (int m = 1; m < 4; m++)

{

for (int n = 2; n < 5; n++)

{

double t = Convert.ToDouble(mg[m - 1, n].Text.ToString());

double b = Convert.ToDouble(mg[m + 1, n].Text.ToString());

double l = Convert.ToDouble(mg[m, n - 1].Text.ToString());

double r = Convert.ToDouble(mg[m, n + 1].Text.ToString());

double center = (l + t + r + b) / 4;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

}

for (int m = 4; m < 5; m++)

{

int n = 1;

double l = Convert.ToDouble(mg[m, n - 1].Text.ToString());

double r = Convert.ToDouble(mg[m, n + 1].Text.ToString());

double center = (100 + t + r) / 12;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

for (int n = 2; n < 5; n++)

{

int m = 4;

double t = Convert.ToDouble(mg[m - 1, n].Text.ToString());

double l = Convert.ToDouble(mg[m, n - 1].Text.ToString());

double r = Convert.ToDouble(mg[m, n + 1].Text.ToString());

double center = (2 * t + r + l) / 4;

mg[m, n].Text =center.ToString();

mg[m, n].BackColor = Color.LightBlue;

}

四、不同初值时的收敛快慢

将从0~100的不同初值代入，均已绝对误差小于等于0.001为收敛条件。以迭代次数为纵坐标，初值为横坐标作图。

1.采用Jacobi迭代计算时的最佳初值

可看到当初值=77℃时，最小迭代次数为7次。

2.高斯（Gauss）-赛德尔（Seidel）迭代计算时的最佳初值

可看到当初值=88℃时，最小迭代次数为7次。

五、上下边界的热流量

上边界热流量：

=

解得，=358.99W。

下边界热流量：

∵下边界绝热∴=0。

六、计算结果的等温图

七、计算小结

本次实验的目的是研究二维稳态导热问题计算效率的影响因素，运用MATLAB，C#等计算工具，从迭代方法，初值设置等方面研究计算效率。得出以下结论：

1.迭代计算方法

G-S迭代计算方法优于Jacobi迭代计算方法。因为G-S迭代方法在每次计算中都用到周围节点的最新计算值，该值比初值更接近于稳态值，所以一定程度上减少了计算步骤。

2.初值设置

通过将不同初值代入并比较运算次数，结合图线我们可以找到迭代计算的最佳初值，在该初值下迭代计算最少的次数便可达到收敛条件。G-S迭代法与Jacobi的最佳初值不同。研究最佳初值的设置可以帮助我们在工程中提高实际传热问题的计算效率。

3.稳态温度等值线

由稳态温度等值线我们可以发现二维稳态导热稳态温度分布的特点：

（1）在外界温度均匀稳定非绝热面，温度等值线平行于该平面。在绝热面，温度等值线垂直于该平面。在各个边界面交界处，温度等值线倾斜分布。

（2）当相邻两边界外界温差大时，等温线在两边界处分布稠密。

通过本次实验，我加深了对运用有限差分方程求解二维稳态导热问题的理解，并掌握了使用计算机软件辅助解题，将课堂知识与实践相结合，对今后更好的学习传热学起到很大的帮助。