【旧代码】传热过程数值模拟（《传热学》实验指导书第四部分第一题，第一，第二类边界条件）

最新推荐文章于 2023-12-26 15:53:07 发布

原创

最新推荐文章于 2023-12-26 15:53:07 发布 · 4.5k 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#stream #up #网格 #lambda #bi #delete

这篇博客分享了作者2010年使用C++编写的一段旧代码，用于实现传热过程的数值模拟。代码针对《传热学》实验指导书第四部分第一题，涉及第一类（给定边界温度）和第二类（对流边界）边界条件。文章介绍了模拟中涉及的网格形状，并提到使用C++进行编程的原因是个人兴趣。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

2010年十月写的旧代码。

第一类边界条件是给定边界温度。

第二类是对流边界。

区域都是如下形状的:

--------------------------------

| |

| ----------------------

| |

--------

用C++纯属蛋疼。

第一类边界条件：

/*
 * 等温边界
 */
#include<iostream>
#include<cmath>
using namespace std;
const double out_temp=30.0;//外边界温度
const double in_temp=0.0;//内边界温度
const double accuracy=0.00000001;//精度
const double lambda=0.53;//导热系数

const int width=15;//上部点阵宽度
const int width_bottom=4;//下部点阵宽度
const int height=4;//上部点阵高度
const int height_bottom=7;//下部点阵高度
/*
//另一组参数
const int width=16;
const int width_bottom=6;
const int height=6;
const int height_bottom=6;
*/
//总共点数
const int num_of_points=width*height+width_bottom*height_bottom;
//单个点
class point
{
		public:
				double temp;//温度
				int up;//数组下标
				int down;
				int left;
				int right;

				point(){
						temp=0.0;//初始化成0摄氏度
						up=down=right=left=0;
				}
};

ostream & operator<<(ostream & src_stream,const point & src){
		src_stream<<"temp="<<src.temp;
		src_stream<<" up="<<src.up<<" down="<<src.down;
		src_stream<<" left="<<src.left<<" right="<<src.right;
		return src_stream;
}

void print_grid(point * points){//输出网格各点的温度
		cout<<endl;
		for(int position=0;position>=0;position=points[position].down){
				for(int tmp=position;tmp>=0;tmp=points[tmp].right){
						cout.width(10);
						cout<<points[tmp].temp;
				}
				cout<<endl<<endl<<endl<<endl;
		}
}

double calc_direction(point * points,int direction,double &opposit_weight){
		//计算给定方向邻点的温度和反方向的权重。
		switch (direction) {
				case -1:
						return out_temp;
						break;
				case -2:
						return in_temp;
						break;
				case -3:
						opposit_weight*=2;
						return 0.0;
						break;
				default:
						return points[direction].temp;
		}
}

void compu_point(point * points,int now){
		//根据周围四个点算出指定点的温度
		double left_temp;
		double right_temp;
		double up_temp;
		double down_temp;
		double up_wight=0.25;//上方权重,默认0.25
		double down_wight=0.25;
		double left_wight=0.25;
		double right_wight=0.25;
		left_temp=calc_direction(points,points[now].left,right_w