计算热物理(共29页).docx

目录244569121416题目：定义在正方形区域（如图所示）H*H=3m*3m 的二维非稳态导热问题的控制方程及其定解条件为其中，c，分别为导热体的密度、比热和导热系数，且=7820kg/m3，c=460J/(kgK)，=15W/(mK)。将定义域各方向均匀分成10个区块，取点中心网格分割方式，各边界共 11个节点; 或者取块中心网格分割，各边界共 10个节点，（两种网格，至少做一种），按以下要求，数值求解该问题。 （1） 取（X，Y）=（x，y）H，=T-T0T1-T0，=tcH2/，将方程和定解条件无量纲化，并在无量纲化条件下求解。（2） 采用以下方法对控制方程离散：（a）用控制容积法离散成显式、全隐式格式；（b）用 ADI格式离散。（3）用 Fourier分析方法对以上3种格式的稳定性进行分析。（4）嵌入初、边条件，分别对以上3种格式构成的定解问题编制计算程序，作数值求解，直到达到稳态温度分布能直接求解 显式格式和 AD