青岛大学传热学上机报告-@孟得明

右边界为第二类边界条件，边界节点 N 的向后差分得： 迭代形式，得： θ N = θ N −1
θ N − θ N −1 = 0 ，将此式整理为 x
1.3.4 最终离散格式
θ1 = θ w = 1 θi =
1 (θi +1 + θi −1 ) 2 + m2 x 2
(i=2,3……,N-1)
θ N = θ N −1
程序 #include"math.h" #include"stdio.h" #define N 5 #define M 5 void main() { printf("\t 二维稳态导热----------孟得明 学号：201040502074\n"); printf("矩形区域，边长 L=W=1，假设区域内无内热源，导热系数为常熟，三个边温度为 T1=0，一个边温度为 T2=1。\n"); double T[N][M],m=0,x,y; int i,j; for(i=0;i<N;i++)
(i=1,2,3,…,N-1)
T1K +1 = T2K +1
TNK +1 =
1 h x 1 (TNK−+ T∞ ) 1 + h x λ +1 λ
其中 K=0，1，2……。 隐式： Ti = T0
0
T1K +1 = T2K +1
Ti K +1 = TNK +1 =
1 K +1 K +1 K ( fTi+ 1 + fTi −1 + Ti ) 1+ 2 f 1 h x 1 (TNK−+ T∞ ) 1 + h x λ +1 λ
K K +1 K K
⎛ ∂ 2T ⎞ = a⎜ 2 ⎟ ⎝ ∂x ⎠ i
K +1
TiK +1 − Ti K ⎛ ∂T ⎞ 将上式的左端温度对时间的偏导数进行差分离散为： ⎜ ⎟ = τ ⎝ ∂τ⎠i ⎛ ∂T ⎞ ⎜ ⎟ ⎝ ∂τ⎠i
K +1
=
TiK +1 − Ti K τ
大学那些事儿：微信账号：college007 欢迎添加！
1.3.5 代数方程组的求解及其程序
假定一个温度场的初始发布，给出各节点的温度初值： θ1 ， θ 2 ，….， θ N 。将这些初 值代入离散格式方程组进行迭代计算，直至收敛。假设第 K 步迭代完成，则 K+1 次迭代计算 式为： θ1
K +1 0 0 0
= θw
θiK +#43;1 −1 2+m x
其中 K=0，1，2……
3.4 程序及截图
程序： #include"math.h" #include"stdio.h" #define N 10 #define M 11 #define f 6.2824114e-2 #define Bi 6.6781312e-2 void main()
大学那些事儿：微信账号：college007 欢迎添加！
∅=
hp (tw − t∞ )th( mH ) m
1.3 数值离散
1.3.1 区域离散
计算区域总节点数取 N。
1.3.2 微分方程的离散
对任一借点 i 有： ⎜
⎛ d 2θ 2 ⎝ dx
⎞ 2 ⎟ − m θi = 0 ⎠
用θ在节点 i 的二阶差分代替θ在节点 i 的二阶导数，得： 整理成迭代形式： θi =
(
)
(i=2,3……,N-1)
K +1 K +1 θN =θN −1
程序: #include"stdio.h" #include"math.h" #define N 11 #define p 8.1818182e-3 void main() { printf("\t 一维稳态导热----------孟得明 学号：201040502074\n"); printf("已知：h=45，t∞=80, tw=300, λ=110, δ=0.01, H=0.1\n"); double t1=300.0,t2=80.0,T[N+1]={1,1,1,1,1,1,1,1,1,1,1,1},k=0,b=0; int i,j; for(i=0;;i++) { k=T[N-1]; for(j=2;j<=N-1;j++) T[j]=(T[j+1]+T[j-1])/(2+p); T[N]=T[N-1]/(1+0.5*p); if(fabs(T[N-1]-k)<=1e-6) break; } printf("坐标\t\t 温度\n"); for(i=1;i<=N;i++) printf("%.2f\t\t%.2f\n",b+=0.01,T[i]*(t1-t2)+t2); }
大学那些事儿：微信账号：college007 欢迎添加！
一维稳态导热的数值计算 1.1 物理问题
一个等截面直肋，处于温度=80 的流体中。肋表面与流体之间的对流换热系数为,肋基处温 度，肋端绝热。肋片由铝合金制成，其导热系数为，肋片厚度为，高度为 H=0.1m。试计算 肋内的温度分布及肋的总换热量。
(K=0,1,…,i=2,3…,N-1)
aτ x2 3.3.3 边界条件的离散
其中 f = 边界节点的差分方程： T1K +1 = T2K +1
TNK +1 =
3.3.4 最终离散格式
显式：Ti=T0 (i=1,2,3,…,N)
1 h x 1 (TNK−+ T∞ ) 1 + h x λ +1 λ
K K K Ti K +1 = ( fTi+ 1 + fTi −1 + (1 − 2 f )Ti )
1.2 数学描述及其解析解
引入无量纲过余温度，则无量纲温度 描述的肋片导热微分方程及其边界条件：
d 2θ − m 2θ = 0 2 dx
x=0,θ=θw=1 x=H,
∂θ =0 ∂x
其中 m =
hp λA ch[ m( x − H )] ch(mH )
上述数学模型的解析解为： t − t∞ = (tw − t∞ ) ⋅
大学那些事儿：微信账号：college007 欢迎添加！
大学那些事儿：微信账号：college007 欢迎添加！
一维非稳态导热的数值计算
非稳态导热问题由于有时间变量，其数值计算出现了一些新的特点。
3.1 问题
一块无限大平板，其一半厚度为 L=0.1m ，初始温度 T0=1000 ，突然将其插入温度
用 I,j 节 点 的 二 阶 中 心 差 分 代 替 上 式 中 的 二 阶 导 数 ， 得 ：
Ti +1, j − 2Ti , j +Ti -1, j x
2
+
Ti , j +1 − 2Ti , j +Ti , j -1 y2
=0
y2 x2 上式整理成迭代形式： Ti , j = (Ti+1, j + Ti-1, j ) + 2( x 2 + y 2 ) (Ti, j +1 + Ti, j -1 ) 2( x 2 + y 2 )
∞
n
2.3 数值离散
2.3.1 区域离散
区域离散 x 方向总节点数为 N，y 方向总节点数为 M，区域内任一节点用 I,j 表示。
2.3.2 方程的离散
对于图中所有的内部节点方程可写为： ⎜
⎛ ∂ 2t ⎞ ⎛ ∂ 2t ⎞ + ⎜ 2 ⎟ =0 2 ⎟ ⎝ ∂x ⎠i , j ⎝ ∂y ⎠i , j
(i=2,3……,N-1)，(j=2,3……,M-1) 补充四个边界上的第一类边界条件得： T1, j = T1 (j=1,2,3……,M) (j=1,2,3……,M) (i=1,2,3……,N) (i=1,2,3……,N)
TN , j = T1 Ti , j = T1 Ti , M = T2
2.4 程序及截图
大学那些事儿：微信账号：college007 欢迎添加！
K K 显式差分格式： Ti K +1 = fTi+1 + (1 − 2 f )Ti K + fTi-1
(K=0,1,…,i=2,3…,N-1)
全隐式差分格式： Ti K +1 =
1 K +1 K +1 ( fTi+1 + fTi-1 + Ti K） 1+ 2 f
T =20 的流体介质中。平板的导热系数 =34.89W/（m ） ，密度 =7800kg/m3,比热 c=712J/ （kg ） ，平板与介质的对流换热系数为 h=233W/（m2 ） ，求平板内各点的温度分布。
3.2 数学描述
由于平板换热关于中心线是对称的， 仅对平板一半区域进行计算即可。 坐标 x 的原点选 在平板中心线上，因而一半区域的非稳态导热的数学描述为：
大学那些事儿：微信账号：college007 欢迎添加！
大学那些事儿：微信账号：college007 欢迎添加！
截图：
二维稳态导热的数值计算 2.1 物理问题
一矩形区域，其边长 L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为 T1=0，一个边温度为 T2=1，求该矩形区域内的温度分布。
2.2 数学描述
对上述问题的微分方程及其边界条件为：
∂2T ∂2T + =0 ∂x 2 ∂y 2