Silver_777的博客

一个形状如下图所示的同心圆套管，两圆内外径分别为17.8mm，46.25mm（比为2.6），内部流体为不可压缩理想气体，使用ANSYS Fluent求其对流换热问题。

章熙民的第六版《传热学》里，较为简单的介绍了非稳态导热的数值计算，本文根据此书，以计算一个可视为无限大平壁的复合墙体传热过程为例，讨论一维非稳态导热问题数值求解的问题。

使用有限单元法计算一维常物性无内热源非稳态导热的温度场，以及内外壁面的热流密度，并进行温度场和热流的特点分析。室内温度恒定为20℃，无限大平壁的比热为1.0kJ/kg℃，室外气温周期性变化，周期为1天（即24小时）。二、网格划分与节点方程式由问题所给条件可知，边界条件为对流边界条件，则该一维常物性无内热源非稳态导热问题的有限单元节点方程如下（注这里不是有限差分法，是有限单元法，离散格式为隐式）：.

牛顿迭代法的具体内容不赘述它的核心算法是：k = 1;x = x0;x0 = x + e*2; % 为了让初启动时满足循环条件while (abs(x-x0))>e && (k<=N) % 同时限定误差和最大循环次数 x0 = x; x = x0 - f(x0)/df(x0); % 牛顿迭代法式子 k = k+1;end这是一个非常简单的牛顿法实现，根据上学期学的数值分析和这学期上的Matlab，把它用GUI简单包装了一下。