传热学第讲非稳态导热.pptVIP

* * 习题解： 有一房间的砖墙厚0.3m，导热系数0.85W/(m·K)，热容量 1.05×106 J/(m3·K)，室内温度20℃保持不变，对流换热系数6 W/(m2·K)。起初该墙的温度处于稳态，内墙表面温度为15℃。后寒潮入侵，室外温度下降到零下10℃，外墙表面换热系数为35 W/(m2·K)。问寒潮入侵后多少时间内墙表面方可感到外界气温的变化？ 模型：常物性无内热源一维非稳态导热 解： （外墙表面） （内墙表面） 2010年10月24日 墙内温度分布示意 假定内外墙表面对流换热系数不变 初始条件的确定 热流密度： 初始温度分布： 外墙表面温度： 数值解（自编程序）： 下面分两种解，自编程序求解和利用FLUET软件求解： 将墙体划分8等分，区域离散后的标记见图 共9个节点 一维非稳态导热网格划分及记号 计算结果 右侧温度开始变化点在时间3492秒 当时间达到2900分时达到新的平衡分布（稳态导热状态） 此时墙体散热热流密度为： 时间为60s 时间为600s 时间为3000s 时间为6000s 时间为3492s 时间为18000s 时间为12000s 时间为2900min 利用FLUENT解： GAMBIT建模： 点击Update钮确定产生工作区图形 物理问题的几何形状 物理问题的几何面 工作区图形 几何形状 几何面 建模坐标系统 x最小值0,最大值30。 y最小值0,最大值10。 incrememt 各取10 划分计算所需要的网格 x向分15等分 y向分4等分 网格点 生成计算域内的网格 网格 设置边界名称和类型 为方便先将网格关闭，点击 在mesh项点Off，再点Apply 在Specify Boundary Types中 对4个边界取名并各设置为Wall类型 设置好的边界面列表如下 关闭对话框 保存数据 点选择菜单中的File，选Export，选Mesh…， 在对话框图中，选定Export2-D(X-Y) Mesh项，File Name框中填入文件名unsteady后，点击Accept，点击Close保存网格文件。 FLUENT计算 启动FLUENT程序 读入网格文件：File→Read→Case 找到文件D:\heatconduction\onedim\unsteady.msh，读入。 检查网格以确保网格质量 确定几何区域的长度单位 Grid Scale Grid Was Created In的下拉菜单中选择cm 注意，此处可选的单位是米，厘米和毫米三种，故在GAMBIT建模时要预先考虑到。 结果： 确定求解方程 对于导热问题，只有一个能量方程（导热微分方程） Define→Models→Energy… 选中 （先计算初始时稳态的温度分布值） 设置导热体的物性参数： Define→Materias…，在材料项选择固体（Solid） 取名brick，导热系数0.85 W/m·K，热容量 1.05×106 J/(m3·K)，（密度1900kg/m3, 比热容552.6J/kg·K） 边界条件设置： Define→Boundary Conditions…， 在边界条件设置对话框内点击Zone区域中的fluid，再点击Type区域中的solid后按Yes按钮、点击OK按钮，计算域为固体。 接下来依次对bottom、left、right和top四个边界设定边界条件 bottom和top分别是绝热条件，设置热流密度为0。 left、right分别是对流换热条件，按题数据设置。 初时核算出外界温度为3.55℃。 解的设置 Solve→Control→Solution…，在Equation框内选择Energy项后点击OK。 初始化： Solve→Initialize→Initialize…，无须设置点击Init。 设置迭代计算前后两次的误差监视窗： Solve→Monitors→Residential…，在对话框中选Plot，点击OK。 迭代计算： Solve→Iterate…，在对话框中填入迭代次数，点击Iterate进行计算。 右图是迭代情况 先按稳态导热进行计算 Define→Models→Solver →选Unsteady，选implicity（隐格式） 按初时的稳态条件计算得到下列温度分布 这就是下面计算所需要的初始温度分布条件 下面开始非稳态问题的计算 更改新的边界条件： 将左侧边界条件中的流体温度更改为263K Define→Boundary Conditions…，left:中设置温度为263K Solve→Iterate…，在对话框中time step size(s)中填入时间步长1s，在Number of time steps中填入6