计算热物理精要.docx

计算热物理大作业  PAGE \* MERGEFORMAT 29  目录  HYPERLINK \l a1 题目 2  HYPERLINK \l a2 1.方程和定解条件的无量纲化 4  HYPERLINK \l a3 2.控制方程的离散 4  HYPERLINK \l a4 3.稳定性分析 5  HYPERLINK \l a5 4.温度场的计算求解 6  HYPERLINK \l a6 5.温度场的等值线图 9  HYPERLINK \l a7 6.分析比较 12  HYPERLINK \l a8 7.个人心得体会 14  HYPERLINK \l a9 附录 16 题目：定义在正方形区域（如图所示）H*H=3m*3m 的二维非稳态导热问题的控制方程及其定解条件为 其中ρ，c，λ分别为导热体的密度、比热和导热系数，且ρ=7820kg/m3，c=460J/(kg?K)，λ=15W/(m?K)。将定义域各方向均匀分成10个区块，取点中心网格分割方式，各边界共 11个节点; 或者取块中心网格分割，各边界共 10个节点，（两种网格，至少做一种），按以下要求，数值求解该问题。 取（X，Y）=（x，y）H，θ=T-T0T1-T0，τ=tρcH2/λ，将方程和定解条件无量纲化，并在无量纲化条件下求解。 采用以下方法对控制方程