FLUENT-第六节传热模型

1、第六节：传热模型概要n能量方程n壁面边界条件n共轭传热n薄壁和双面壁n自然对流n辐射模型n报告-输出能量方程n能量输运方程: 单位质量的能量 E : 对可压缩性流体，或者密度基求解器，总是考虑压力做功和动能。对压力基求解器计算不可压流体，这些项被忽略，可以用下面的命令加入: define/models/energy? ConductionSpeciesDiffusionViscousDissipationConductionUnsteadyEnthalpy Source/Sink固体域的能量方程n能计算固体域的导热n能量方程: h 显焓:n固体域的各向异性导热系数（压力基求解器）壁面边界条件n

2、五类热边界条件 热流量 温度 对流 模拟外部环境的对流（用户定义换热系数） 辐射 模拟外部环境的辐射（用户定义外部发射率和辐射温度） 混合 对流和辐射边界的 结合.n壁面材料和厚度可以定义n为一维或壳导热计算共轭传热nCHT固体域的导热和流体域的对流换热耦合n在流体/固体交界面使用耦合边界条件Coolant Flow Past Heated RodsGridVelocity VectorsTemperature Contours共轭传热例子Circuit board (externally cooled)k = 0.1 W/mKh = 1.5 W/m2KT = 298 KAir inletV

3、= 0.5 m/sT = 298 KElectronic Component(one half is modeled)k = 1.0 W/mKHeat generation rate of 2 watts (each component)Top wall(externally cooled)h = 1.5 W/m2KT = 298 KSymmetryPlanesAir outlet问题设置-热源n在固体域加入热源模拟电子部件的生成热温度分布FlowdirectionConvection Boundary1.5 W/m2 K298 K free stream temp.Convection bo

4、undary1.5 W/m2 K298 K free stream tempFront ViewTop View(image mirrored about symmetry plane)Elect. Component(solid zone)2 Watts sourceBoard(solid zone)Air (fluid zone)298426410394378362346330314Temp.(F)Flowdirection替代的模拟策略n可替代的策略为模拟壁面为一有厚度面 (Thin Wall model).n这时，不需对固体域划分网格对固体板划分网格 vs. 薄壁方法n对固体板划分网格

5、 在固体域求解能量方程l. 板厚度需用网格离散 最精确的方法，但需要多计算网格 由于壁面两侧都有网格，总是应用耦合热边界条件Fluid zoneSolid zoneWall zone(with shadow)Wall thermal resistance directly accounted for in the Energy equation; Through-thickness temperature distribution is calculated.Bidirectional heat conduction is calculated.对固体板划分网格 vs. 薄壁方法n薄壁方法 人

6、工模型模拟壁面热阻 壁面需要必要的数据输入（材料导热系数，厚度） 只有对内部边界用耦合边界条件Fluid zoneWall zone(no shadow)Wall thermal resistance is calculated using artificial wall thickness and material type. Through-thickness temperature distribution is assumed to be linear.Conduction is only calculated in the wall-normal direction unless S

7、hell Conduction is enabled.壳导热模型n壳导热模型处理板内部的导热n求解器创建额外的导热单元，但不能显示，也不能通过UDF获得n固体属性必须是常数，不能和温度相关Static Temperature(cell value)Virtual conduction cells自然对流n当流体加热后密度变化时，发生自然对流n流动是由密度差引起的重力驱动的n有重力存在时，动量方程的压力梯度和体积力项重写为：:其中自然对流 Boussinesq 模型nBoussinesq 模型假设流体密度是不变的，只是改变动量方程沿着重力方向的体积力 适用于密度变化小的情况 (例如，温度在小范围

8、内变化).n对许多自然对流问题，Boussinesq 假设有更好的收敛性 常密度假设减少了非线性. 密度变化较小时适合. 不能和有化学反应的组分输运方程同时使用.n封闭空间的自然对流问题 对稳态问题，必须使用 Boussinesq 模型. 非稳态问题，可以使用 Boussinesq 模型或者理想气体模型自然对流的用户输入n在操作条件面板中定义重力加速度n定义密度模型 Boussinesq 模型 激活重力项. 设置操作温度 T0. 选择 Boussinesq 模型，输入密度值 0. 设置热膨胀系数 . 使用温度变化模型 (ideal gas, Aungier-Redlich-Kwong, pol

9、ynomial): 设置操作密度或 让 FLUENT 从单元平均中计算 0 辐射n当和对流及导热换热相比， 量级相当时，应该考虑辐射效应 ， Stefan-Boltzmann常数, 5.6710-8 W/(m2K4)n要考虑辐射，需求解辐射强度输运方程RTEs 当地流体对辐射能的吸收，以及边界对辐射的吸收，把RTEs 和能量方程耦合起来 这些方程常常和流动方程分离求解，然而，他们也可以和流动耦合n辐射强度, I(r,s),和方向及空间是相关的nFLUENT中有五个辐射模型 离散坐标模型 (DOM) 离散传输辐射模型 (DTRM) P1 模型 Rosseland 模型l Surface-to-S

10、urface (S2S)选择辐射模型n指南： 计算代价 P1 计算代价小，有合理的精度 精度 DTRM 和 DOM 最精确. 光学厚度 DTRM/DOM 适合光学厚度小的模型 (L 1) P1 适合光学厚度大的模型. S2S 适合零厚度模型 散射 只有 P1 和 DO能考虑散射 颗粒辐射 P1 和 DOM 能考虑气体和颗粒间的辐射换热 局部热源 适合用DTRM/DOM 带足够数量的射线/坐标计算附录太阳辐射模型n太阳辐射模型 太阳辐射能量的射线追踪算法，和其他辐射模型兼容 允许并行计算（但射线追踪算法不能并行） 仅适用3Dn特点 太阳方向向量 太阳强度（方向，散射） 使用理论最大或气象条件计算

11、方向和方向强度 瞬态情况 当方向向量是用太阳计算器算出的化，瞬态计算中太阳方向矢量会随时间改变 设置 “time steps per solar load update”能量方程源项 粘性耗散n粘性耗散引起的能量源项: 也称为粘性加热 对粘性剪切力大的流体（如润滑油）和高速可压缩流动比较重要 常常忽略 缺省的压力基求解器不包括. 密度基求解器一般包括. 当 Brinkman 数接近或超过1时重要能量方程源项 组分扩散n多组分流中因为组分扩散引起的能量源项： 包括了由于组分扩散引起的焓输运效应 密度基求解器总包含 在压力基求解器中可以不显示此项能量方程 源项n化学反应流中由于化学反应引起的能量源

12、项 所有组分的生成焓 所有组分的体积生成率n由于辐射引起的能量源项n相间能量源项: 包括连续相和离散相间的传热 DPM, 喷雾, 颗粒薄壁中的温度分布n薄壁模型应用于法向导热，不生成实际的单元n壁面热边界条件应用于外层Thermal boundary condition on wallStatic temperature(cell value)Thin wall(no mesh)Wall temperature(outer surface)Wall temperature(inner surface)薄壁和两侧壁面n薄壁方法中，壁面厚度不需划分网格n在两个区域之间模拟薄层的材料n求解器施加热阻

13、 x/kn边界条件施加在外层面上Thermal boundary conditions are supplied on the inner surface of a thin wallExterior wall(user-specified thickness)Fluid or solid cellsOuter surface(calculated)Inner surface(thermal boundary condition specified here)xInterior wall(user-specified thickness)1k2kInterior wall shadow(use

14、r-specified thickness)22 or Tq11 or TqThermal boundary conditions are supplied on the inner surfaces of uncoupled wall/shadow pairsFluid or solid cellsFluid or solid cells离散坐标模型AbsorptionEmissionScatteringn在有限的离散立体角度s上求解辐射输运方程 :n优势: 守恒方法能保证粗的离散方式上实现热平衡 通过更密的离散方式能提高精度 最综合性的模型: 考虑了散射、半透明介质、镜面以及波长相关的灰体

15、模型n局限性: 求解大数量坐标耗费CPU过多离散传输辐射模型 (DTRM)n主要的假设 特定范围角度的离开表面的辐射能用一束射线近似n使用射线跟踪技术，沿着每条射线积分辐射强度n优势: 相对简单的模型 增加射线数量能提高精度 适用大范围的光学厚度n局限性: 假设所有表面是漫射的. 不包括散射. 求解大数量的射线耗费CPU过多.P-1 模型n主要假设 对RTE积分后，和方向不再相关，导出入射辐射的扩散方程n优势: 辐射传热方程更易求解，耗费资源少 包括散射效应 颗粒、液滴和烟灰的影响 对光学厚度大的应用（如燃烧）较合理n局限性: 假设所有面都是漫射的 如果光学厚度小的话，可能导致精度损失（取决于

16、几何的复杂性） 对局部热源或汇，预测的辐射热过高Surface-to-Surface (S2S) 辐射模型nS2S辐射模型用于模拟介质不参与的辐射 例如，太空飞船的排热系统、太阳能搜集系统、辐射加热器、汽车发动机舱散热等 S2S 是基于角系数的模型 假设没有介质参与n局限性: S2S 模型假设所有面是散射的 假设是灰体辐射 随着表面数量的增加，存储和内存增加很快 可以使用面族来减少内存使用 面族不能和滑移网格及悬节点同时使用 不能使用于周期性或对称边界条件输出 ANSYSn 输出扩展名为 .rfl的 ANSYS 结果文件，读入到ANSYS的顺序为： 1. 在 ANSYS中，到“ General

17、 Postproc Data”及“ File Options”，读入FLUENT生成的文件2. 到“ Results Summary ”，点击第一行，能看到ANSYS_56_OUTPUT窗口显示的几何信息3. 在 ANSYS 输入窗口，键入下面的命令: SET,FIRST/PREP7ET,1,142最后一个命令对应 FLOTRAN 3D 单元，如果你使用二维计算，应改为：ET,1,141. 4. 在 ANSYS MULTIPHYSICS UTITLITY 菜单，选择 Plot 及 Nodes 或Elements, 在 下拉窗口的Results中，选择包括节点输出 ANSYS n通过 GUI 或

18、 TUI输出ANSYS 文件 /file/export/ansys file-namen文件包括坐标、连接关系及下面 的标量： Density, viscosity X, Y, Z velocity, pressure, temperature Turbulence kinetic energy, turbulence dissipation rate, turbulent viscosity, effective viscosity Thermal conductivity (laminar, turbulent, effective) Total pressure and temperature, pressure coefficient, Mach number, stream function, heat flux, heat transfer coefficient, wall shear stress, specific heat输出 ABAQUSn输出