热工数值设计报告.docVIP

《热工数值计算》程序 设计报告 1 设计题目 利用试射法计算一维等截面圆柱形直肋的温度分布和散热量。并与分析解进行对比。截面面积为5mm2，材料的导热系数为14W/(m•℃)，肋高100mm,肋根温度为320℃，肋端绝热，环境温度为20℃，总表面传热系数为20 W/(m2•℃)。 2．物理与数学模型 2．1物理模型 如图1示温度为320℃的壁面上有一圆柱形直肋，肋片的形状和对象的各参数如图1中所标示。 图1 一维等截面圆柱形直肋 2．2数学模型 为了建立描述该一维等截面圆柱形直肋的温度分布的数学描述，建立如图1示的坐标系，并从肋片中距肋基x处取出长为dx微无体，如图2所示 图2 研究对象的微元体模型 在稳态的条件下，微元体满足如下的能量守恒， Qx= Qt+ Qx+dx 式中，Qx和Qx+dx 分别为在x处导入微元体的热量和在x+dx处导出的微元体的热量。Qt是通过微元体侧面散布失到环境中的热量。 根据Fourier定律，Qx和Qx+dx可表示为， 式中λ为肋片材料的导热系数；A为肋片的截面积。 通过肋片侧表面积的散热量Qt可表示为， 式中，h为肋片表面传热系数，F为微元体的侧表面积，它应等于 式中，U为圆柱形微元体的周长，dx为微元体的长度。 将Qx、Qx+dx和Qt的表达式代入到能量守恒方程式中，有 （1） 相应的边界条件为 （2） （3） 此即描述该一维等截面圆柱形直肋的温度分布的数学模型。 3 数值处理与程序设计 3.1 将数学模型改写为试射法的形式 所建数学模型为边值问题，通过降阶处理，方程（1）可以转换为一阶常微分方程组，然后利用试射法进行数值计算，即， 令 则有 这样有 由，有 由，有 3．2程序的编写及基本算法 根据上述分析编写用户子程序，利用试射法主体程序确定m值，将边值问题转化为它所对应的初值问题，然后再调用Runge-kutta方法计算各点的温度分布及温度梯度，其中试射法的基本算法如下， 1．先任意给，求它们对应的初值问题，得； 2．若时，则按下式插值 转向步骤（4）；否则，转向步骤（3）； 3．若,则 否则， 然后针对新的再解初值问题得，直到满足 4．利用求得的m解初值问题得,计算; 5．若|,求得m，转（7）；否则转（6）； 6.若,则 否则， 再插值求得m,转(4); 7.结束。 而用户子程序需要根据本题分析编写，具体如附录2所示。 4 模型与程序验证 4．1分析解计算 根据方程（1）（2）（3），利用传热学有关知识，计算分析解如下， 引入过余温度， 则原方程（1）（2）（3）可变为， （4） （5） （6） 式中为一常数。 解得该圆柱形直肋的温度分布为， 也即， 也可设计分析解程序，来得到分析解结果。 4．2数值解与分析解结果比较验证 利用试射法程序，得到数值解的结果，提取数据，再利用分析解程序，解得分析解结果，提取数据，作表格如附录3表1数据比较结果，根据表格数据，在同一坐标系中做出数值解与分析解的该肋温度分布曲线如下图3， 图3 温度分布曲线比较 不论是从表中的数据还是温度曲线的分布曲线，都可以看出数值解与分析解之间在误差允许的范围内的吻合，以此证明了数值解计算的正确性。 5 计算结果与分析 5．1温度分布的数值结果 根据试射法程序，可以得到各