板翅式换热器数值模拟

第一章 绪论1.1 板翅式换热器的发展1.1.1 板翅式换热器总体发展概述 板翅式换热器首先应用于汽车和航空工业中，20世纪30年代英国马尔斯顿艾克歇息=尔瑟公司首先生产了铜质钎焊的板翅式换热器，20世纪40年代中期出现了更轻巧的铝质钎焊板翅式换热器，随后研究、生产了更多结构形式的翅片，使其趋于更加紧凑、轻巧，20世纪50年代开始在空气分离设备中应用板翅式换热器，随着冶金、石化工业对空分设备的大量需求以及空分设备大型化的发展趋势，板翅式换热器的研究、实验、设计与制造也得到有力的推进。目前它已广泛应用于航空、汽车、内燃机车、工程机械、空分、石化、制冷、空调、深低温等领域。 板翅式换热器在其初期发展阶段，由于人们对于这种传热表面的传热机理及设计数据缺乏认识，加上结构与制造工艺方面存在许多问题，因此在相当长的一段时间内，仅处于小规模的试验阶段。1942年，美国的诺利斯首先进行了平直形翅片、锯齿翅片、波纹翅片、钉状翅片的传热机理研究，找出几种主要翅片的摩擦因子、传热因子与雷诺数的关系，为后来的研究、设计与应用奠定了基础。1947年美国海军研究署、船舶局、航空局合作在斯坦福大学拟定了系统的研究计划并扩大了研究范围，在19481954年间，美国海军研究署公布了22篇关于紧凑式换热表面的实验研究报告。后来凯斯和伦敦两人编著了紧凑式换热器，较系统的总结了研究成果，在目前这已成为研究、设计与应用板翅式换热器的基本参考文献。-引自换热器设计手册1.1.2 板翅式换热器在我国的发展 我国板翅式换热器的研究开始于20世纪60年代中期，由杭州制氧机厂、营口通用机械厂、开封空分设备厂、上海第一五金厂等单位研制并相互协作，先后研制了6000、3200、10000等空分设备上用的板翅式换热器，在制造工艺上虽然取得一些经验，但也存在一些问题，后来机械工业不组织了攻关小组，重点研究制造工艺中的某些关键性问题。1972年，机械工业部系统的总结了研制与攻关的成果，为我国的板翅式换热器的发展打下一个良好的基础。 现在我国已经掌握了成熟的板翅式换热器的设计、制造技术，我们已经拥有自己的设计软件，根据这些软件设计计算的结果，经过实际考核都达到了预期的结果。我们不仅掌握了盐浴炉的钎焊技术，而且自行设计、制造了大型真空钎焊炉，成熟而批量的生产力大、中、小型各种真空钎接的板翅式换热器，我国生产的板翅式换热器不仅完全满足国内各种型号空分设备的需要，而且还有一定数量的出口或满足世界其他厂商的订货。对生产单元尺寸在1000mmX1000mmX6000mm左右，承压能力在23Mpa以下的板翅式换热器已不存在问题。但对于石化工业中的某些大型、高压、有相变的多股流板翅式换热器的设计、制造，由于过去这些石化设备基本上都是进口，国内行业没有机遇进行研制，所以这方面有待于开发。-引自换热器设计手册1.2 板翅式换热器数值模拟的发展状况随着计算流体动力学(Computationall Fluid Dynamics简称CFD)和计算传热学的蓬勃发展，使得采用数值模拟的方法进行传热、传质、动量传递及燃烧、多相流和化学反应等的研究成为可能。CFD在最近20年中的飞速发展，主要因为板翅式换热器无论分析的方法或实验的方法都有较大的限制，CFD通过引入计算流体力学技术，能使设计思想迅速可视化并以较低成本进行测试，从而大大缩短其从概念到最终实现的时间。 因此，国内的部分学者通过利用CFD对板翅式换热器进行了大量的研究。祝银海针对平直形和锯齿形两种不同翅片类型进行了数值模拟研究。王武林应用SMPLE算法对错列翅片板翅式换热器换热表面的流动及传热进行数值模拟。漆波对低雷诺数下白叶窗翅片内的传热和流动特性进行了数值模拟。 板翅式换热器的传热与流动阻力性能主要决定于翅片的表面特性。 李媛通过稳态试验法，对锯齿形、波纹形、百叶窗形铝翅片表面性能进行了研究，取得了相应的Re-j和Re-f曲线。指出随着雷诺数的增大，翅片表面的传热因了下降，摩擦因了也下降。指出平直翅片的高度对其传热性能有直接影响，翅片间距越大，传热效果越好;锯齿翅片的切开长度越短，传热性能越好;波纹翅片的波幅越大，翅片间距越大其传热效果越好。庞铭采用计算机编程对错位翅片型板翅换热器的传热因子、摩擦因子及翅片结构对传热性能影响进行了数值模拟，并指出在需要较大传热系数的场合，选用高度小、翅片厚度大、有效长度短、间距小的翅片;相反，在需要传热系数较小场合以选用高度大、翅片厚度小、有效长度长、间距大的翅片为宜。李建军通过对油水板翅式换热器进行的性能试验，得到了低雷诺数流动下板翅式换热器翅片侧传热与阻力特性的数据，在此基础上获得了错位翅片传热因子与摩擦系数的准则关系式，传热因子和摩擦系数的最大计算误差分别为62%和1.44%。根据这些准则关系式提出了一个衡量翅片质量的经济系数。-板翅式换热器的研究进展1.3 板翅式换热器的结构与传热机理简述 板翅式换热器由两块板片和翅片钎焊在一起，两端以封条紧固形成的流道单元组成。冷热流体在流道中流动，流道间隔布置达到换热的目的。翅片是板翅式换热器的主要结构，也是换热的主要结构。翅片作为换热器的二次表面扩充了换热表面积，同时扰动流体流动，使边界层不断破裂和再生从而达到增强换热的目的。另外，翅片在板片间的支撑和加固作用，提高了板翅式换热器的强度和承压能力。如图1-1（a）图1-1 板翅式换热器结构及翅片型式板翅式换热器的翅片形式经过多年的发展，主要的形式有平直型翅片、锯齿型翅片、打孔型翅片和波纹型翅片。如图1-1（b）（c）（d）（e） 平直型翅片由薄金属片冲压或滚轧而成，其换热和流动特性与管内流动相似。相对其它结构形式的翅片，其特点是传热系数和流动阻力系数都比较小。这种翅片一般用于流动阻力要求较小而其自身的传热系数又比较大的场合。平直翅片具有较高的承压强度。 锯齿型翅片可看作是由平直翅片切成许多短小的片段，并且互相错开一定间隔而形成的间断式翅片。这种翅片对促进流体的湍动、破坏热阻边界层十分有效，属于高效能翅片，但流动阻力也相应增大。锯齿型翅片多用于需要强化换热的场合。 打孔型翅片是先在金属片上冲孔，然后再冲压或滚轧成型。翅片上密布的小孔使热阻边界层不断破裂，从而提高了传热性能。打孔有利于流体分布，但同时也使翅片的传热面积减小，翅片强度降低。打孔型翅片多用于导流片及流体中夹杂着颗粒或相变换热的场合。 波纹型翅片是将金属片冲压或滚轧成一定的波形，形成弯曲流道，通过不断改变流体的流动方向，促进流体的湍动、分离和破坏热阻边界层，其效果相当于翅片的折断。波纹愈密、波幅愈大，越能强化传热。 板翅式换热的流动形式又大致可以分为逆流、叉流和叉逆流。如图1-2图1-2 a. 逆流式 b.叉流式 c. 叉逆流式此外，板翅式换热器还有一个重要的结构是导流片。导流片位于流道的入口处和出口处。其作用就是为了引导由进口管经封头流入流道的流体，使之均匀的分布于流道之中，或是汇集从流道流出的流体使之经过封头由出口管排出。导流片尚有保护翅片，壁面通道堵塞的作用。常见的导流片形式如图1-3 。图1-3 导流片的形式第二章 板翅式换热器的设计理论jc posted 2010年5月25日 06:39 in 毕业论文（不再更新） , 520 阅读 2.1 板翅式换热器结构参数定义2.1.1 水力半径和水力直径水力半径被定义为最小自由流通面积和湿周之比 (2-1)其中： -总传热面积 -流体流动长度水力直径等于4倍的水力半径： (2-2)值得注意的是，水力直径和当量直径不是相同的概念。水力直径是非圆形截面等效为圆形截面管道的一个几何尺寸，用于计算雷诺数和压力损失，判断管道内流体是层流还是湍流状态。当量直径是指，在压力损失相等的前提下，非圆形截面管道与圆形截面管道的一个等效直径。通常是靠实验总结出的经验公式获得。在本文中，由于主要计算为准则数，所以统一使用水力直径。2.1.2 传热面积密度、翅片面积比和孔度表面密度被定义为换热器一侧（即热流体或冷流体侧）的总传热面积与该侧板间体积之比： (2-3)式中的为孔度，被定义为最小自由流通面积与迎风面积 的比： (2-4)翅片面积比被定义为换热器一侧的翅片表面积与总传热表面积之比： (2-5)需要指出的是，对于换热器一侧的总换热表面积A，由一次传热面积（山和二次传热面积（翅片面积）两部分组成： (2-6)2.1.3 翅片效率和表面效率等截面直翅的翅片效率： (2-7)其中： (2-8)-翅片导热系数- 表面换热系数（对流换热系数）对于一般的平直矩形，三角形和锯齿形翅片均可近似为等截面直翅。 具体的结构参数如下图2-1、图2-2所例示：图 2-1图 2-2-翅片间距-翅片厚度-隔板厚度-板间距-热流体流动长度-冷流体流动长度-换热器高度（包括侧板） -侧板厚度-封条宽度-翅片高度-翅片宽度-热流体流道数-冷流体流道数2.2 板翅式换热器的热流量方程和效率在板翅式换热器中，定义热流体进出口温度分别为：，定义冷流体的进出口温度分别为：那么热流体和冷流体的换热量分别为：其中，为流体的热容量，单位。即是对应单位温度变流动流体的能量储存速率。那么换热器的换热量可以定义为：(2-9)在换热器的不同部位，壁面两侧的流体温度是不一样的。同样，在换热器的不同部位换热系数也是不一样的。因此，此公式中的和是微元面积处的温差和换热系数。取平均温度来简化换热热流量方程 (2-10)此时的乃是板翅式换热器的换热系数的一个总体平均值，与式(2-8)中的意义不同。是翅片的换热系数，其表征翅片表面的对流和流体本身的传热效果。是以换热器及出口温度为基准算得的流体在换热器流道中的温度分布平均值，其计算方法有对数和算术两种方法。对数平均温差因流动方式而异，分别为：顺流：(2-11)逆流：(2-12)对于其它的流动方式，可采用基准于逆流的修正因子来近似表示。的意义为接近逆流方式平均温度的程度，一般根据无因次量来查表。(2-13)(2-14)在工程计算中，在热容量较小时可以以算术平均值简化平均温度的求解。当时，可取算术平均温度作为其平均温度：当时，侧取算术平均温度：侧的平均温度为：若为热侧时则取加号。令式(2-10)：其中表示和的较大者。采用的原因是热容量小的流体温差较大，其误差相对较小。由(1-2)式定义传热单元数NTU为：(2-15)换热器的效率定义为换热器的实际换热量和理论上最大换热量的比值，即(2-16)当时，有当时，有换热器中理论最大的换热量只有在换热表面积无限大的逆流中出现。在理想化的情况下，热流体可以被冷却到,或者冷流体可以被加热到,那么冷热流体的最大温差为对于,则对于,则在各种布置方式的换热器中，其效率和单元数之间存在着一定的函数关系。逆流单流程流动：(2-17)(2-18)对于两种流体均不混合的单流程叉流流动，Mason应用拉普拉斯变换法，得出效率、单元数和热容比的关系式是一个无穷级数解。工程计算中目前采用Drake提出的一个近似关系：(2-19)2.3 板翅式换热器的物性参数 板翅式换热器内部流体的温度在其空间位置中是变化的，因此其物性参数都随温度或压力等基本参数而发生变化。在工程计算中，一般取入口处的温度、压力值来作为求解物性参数的参考值。 在气体温度为常温，压力不高的情况下（温度低于200K，压力高于100atm），作为理想气体处理，其误差不大。理想气体的气体状态方程： (2-20)其中为比体积，为摩尔气体常数，空气为为摩尔质量，空气为 气体的粘度主要由分子动量交换强度决定，当温度升高时，分子运动加剧，动量交换剧烈，表现切应力增大，使粘度也相应增大。在一般情况下，压强变化对粘度几乎没有什么影响，只有发生几百个大气压变化时，粘度才有明显变化。因此可认为粘度至于温度相关，用Sutherland关系式计算： (2-21)其中为气体在时的动力粘度，为Sutherland常数，为参考温度。空气： , ,。 物体的导热系数的定义式由傅里叶定律的数学表达式给出，数值上它等于在单位温度梯度作用下物体内热流密度矢量的模：第三章 物理模型的建立jc posted 2010年5月26日 21:13 in 毕业论文（不再更新） , 238 阅读 为了模拟板翅式换热器的换热和流动特性，首先设计一台空气空气板翅式换热器，以作为数值模拟的模型基础和参考。3.1 原始数据表3-1 原始数据 热空气 冷空气入口温度 180 40入口压力 0.2648 0.1013入口质量流量 3200 5600换热器效率 压降 流动长度 480 2703.2 平均温度及物性参数计算首先假设换热器的效率为75%，则：假设，则：由于取算术平均温度作为其平均温度：根据平均温度，求冷热空气的物性参数。同理同理同理同理表3-2 冷热空气的物性参数参数 热空气 冷空气 3.3 结构参数的计算表3-3 产品结构尺寸表 热侧 冷侧型面 平直型 平直型翅片节距 2.0 2.1板间距 4.7 6.5翅片层数 7 8翅片厚度 0.3 0.3隔板厚度 0.5侧板厚度 2封条厚度 6换热器高度：翅片高：翅片宽：水力直径：翅片面积比：传热面积密度：空气流通面积：空气迎风面积：孔度：一次传热面积：二次传热面积：3.4 传热计算质量流量：雷诺数：由雷诺数可知，流体处于充分发展的湍流状态，采用Gnielinsk公式：翅片参数：有效传热面积：两侧总有效传热面积：两侧表面效率：壁面热阻：忽略污垢热阻，则：传热单元数：德思克近似公式：换热器传热热流量：流体出口温度：与假设效率误差：3.5 阻力计算热侧阻力：设：由，得，查图2-4得设局部阻力损失系数为，则局部损失为：压降符合设计要求。冷侧预设，采用与热侧相同的计算方式可得，也符合设计要求。第四章 数值模型的建立jc posted 2010年5月29日 18:18 in 毕业论文（不再更新） , 421 阅读 4.1 计算流体动力学简介 计算流体动力学（Computatonal Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所作的分析。CFD的基本思想可以归结为：把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。（周雪漪 计算水力学 清华大学出版社 1995 陶文铨 数值传热学（第二版）西安交通大学出版社 2001 郭鸿志 传输过程数值模拟 冶金工业出版社 1998） 采用CFD的方法对流体流动进行数值模拟，通常包括以下步骤： （1） 建立反映工程问题或物理问题本质的数学模型。具体说来就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 （2）寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是CFD的核心。 （3）编制程序和进行计算。这部分工作包括计算网格划分、初始化和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Navier-Stokes方程就是一个十分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正式从这个意义上讲，数值模拟又叫数值试验。 （4）显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。 以下将以以上4点为纲领，详细叙述模型的建立和计算过程。 本文使用商业CFD软件Fluent及其前置软件Gambit来进行模拟，使用压力基解算器求解。压力基求法原称为分离求解方法，即分别求解各个控制方程的方法。由于控制方程是非线性的，因此求解必须经过多次迭代才能获得收敛解。图4-1为求解流程图，其过程可概述如下： （1）流场变量更新。在第一次计算时，变量由初始化过程更新。在随后的计算中，没迭代一部既得到一个更新的解。 （2）用当前压强和质量通量的值求解动量方程，以得到新的速度场。 （3）因为（2）中得到的速度场的数值解无法完全满足连续方程，于是再求解压强修正方程。压强修正方程是由连续方程导出的泊松型方程，求解这个方程可以得到对压强场、速度场和质量通量的修正，进而使连续方程得到满足。 （4）利用前面求出的解，求解湍流方程、能量方程和组元方程。 （5）在多项流计算中那个如果考虑相间干扰，则需要通过求解弥散相轨迹计算得到连续相方程中的源项解。 （6）检验收敛条件是否满足。如果收敛条件满足，则停止计算。如果计算没有收敛，则继续迭代过程。图 4-14.2 流体的控制方程及其离散化设为流体系统在时刻所具有的某种物理量，表示单位质量流体所具有的这种物理量。由此推导出输运方程为： (4-1)其中为速度在微元面法线方向的投影。该式表明，流体某一物理量时间的全变化率等于控制体内这种物理量的时间变化率与经过控制面的这种物理量的净通量之和。单位质量的流体的质量，质量为因为质量守恒，所以有：由(4-1)式可得： (4-2)此即是连续方程的积分形式，其微分形式为： (4-3)在直角坐标系中表示为： (4-4)单位质量的流体动量为，流体系统动量为，由(4-1)式得： (4-5)根据质点系动量定理，流体系统动量的时间变化率等于作用在系统上的外力之和。作用在流体流体微元上的外力有质量力和表面力之分。定义为作用在单位质量流体上的质量力分布函数，定义为作用在单位质量流体上的表面力分布函数。那么运动方程的完整积分表述方式为：(4-6)又根据奥高定理：于是式(4-6)变为： (4-7)此即是运动方程的微分形式在直角坐标系中，变为： (4-8a) (4-8b) (4-8c)由本构方程：代入式（4-8）得到N-S方程： (4-9a) (4-9b) (4-9c)当流体为均质不可压缩常粘度时，N-S方程简化为： (4-10a) (4-10b) (4-10c)单位质量的流体的能量为，其中为单位质量流体的热力学能。则，由（4-1）式得： (4-11)根据能量守恒和转换定律，流体系统中能量的时间全变化率等于作用在系统上的质量力和表面力所作的功率和与外界换热率之和。则得： (4-12)则能量守恒方程积分形式为： (4-13) 其微分形式为： (4-14)对于不可压缩，常物性无内热源可简化为： (4-15) 比较基本的控制方程，可以看出，尽管这些方程中因变量各不相同，但它们均反映了单位时间单位体积内物理量的守恒性质。用 表示通用变量，则上述控制方程可以表示成以下通用形式： (4-16)其中是变量的扩散系数 是的单位体积源项转化为微分形式为： (4-17)式中从左到右依次代表瞬态项、对流项、扩散项和源项。对反映物理问题的控制方程离散化是数值模拟里核心的 步骤，离散化方式的选择决定着计算的稳定性和精度。目前成熟和应用较多的离散化方式有：有限差分法（Finite Difference Method）、有限元法（Finite Element Method）和有限体积法（Finite Volume Method）。 有限体积法的基本思路是：将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；将待解的微分方程对每一个控制体积积分，从而得到一组离散方程。其中的未知数是网格点上的因变量。为了求出控制体积的积分，必须假定值在网格点之间的变化规律。在有限体积法的网格中，计算节点（node）被看成是控制体积（cell）的代表。在离散化过程中，将一个控制体积上的物理量定义并储存在该节点处。在有限体积法中，计算节点的值即因变量处于控制体积的中心，有异于网格交点。其数值的确定方式和在节点间的数值传递（变化）方式即是离散化方式。 Fluent对于时间上的偏导数，有一阶和二阶的显式或隐式格式可选。由于本文讨论的问题为定常流动，因此控制方程中没有瞬态项。对于空间上的离散方式即对对流项和扩散项的值的离散，通常有一阶迎风、二阶迎风、中心差分和QUICK等格式。一阶迎风格式将控制体积界面上的未知量恒取上游的计算节点的值，即：而二阶迎风格式在一阶迎风的基础上考虑了物理量在节点间分布曲线的曲率影响，其表示为：其中分别为计算节点的值和其在上游节点的梯度 为上游节点质心道界面质心的方向向量中心差分格式计算控制体积界面上的值为：其中指共面于的控制体的计算节点值 本文所建立的模型流动方向与网格方向基本一致，且为结构化网格，因此选用一阶迎风格式。一阶迎风格式具有稳定性高，计算速度快的优点。对于的离散化，在Fluent中使用Green-Gauss理论，分别有Cell-Based和Node-Based两种方法。对于控制体中心的值，其梯度的离散化形式为： (4-18)Cell-Based方法： (4-19)Node-Based方法： (4-20)其中界面上的节点数 包围节点的单元的值的加权平均Node-Based方法比默认的Cell-Based方法精度更高，因此本模型选用Node-Based方法计算。4.3 几何模型的建立及网格化 本文用Gambit建立两个平直型翅片的模型，其尺寸参数分别为： 模型一 模型二水力直径mm 2.45 2.79翅片间距mm 2.0 2.1翅片厚度mm 0.3 0.3翅片高度mm 4.7 6.5隔板厚度mm 0.5 0.5流动长度mm 50 50 由于翅片很薄，其与板片接触部分以钎焊连接，换热量少。为了简化模型，减少网格数量节约计算资源将此部分去除。此外，参考1947年美国海军研究署的实验装置，取板翅式换热器一个完整流道和其相邻的两个流道的一半，建立模型结构如图4-2。图 4-2 图中1-2-3部分为流体区域，上下板片及翅片被分割为12部分。注意在Gambit中对体的切割有连接（connected）和非连接（unconnected）之分，非连接的体在其交接面出产生连个重叠的面，在载入Fluent时默认为两个壁面（wall）。如果两个壁面间有参数的传递，比如流体/固体内部或流固耦合处。则这种设置将阻断参数的传递，使其彼此成为孤立的体。若要使其之间有参数的传递，即使其网格结构不同或属于流固两种区域，则需要使其连接或分别定义为interface边界条件。interface边界条件为不同网格结构的区域提供了一种插值的参数传递方式，然而其设置繁杂，会增加计算量且存在一定误差。此处体的分割均为连接。即两个分割的体公用一面，在载入Fluent时会自动生成wall和wall.shadow两个壁面。其中一个为虚拟壁面，Fluent自动耦合两壁面间参数。 由于平直型翅片几何结构较为规整，因此使用六边形的网格单元以map的方式分别网格化图 4-2中的各区域。考虑到网格密度对计算精度的影响，先行对各待网格化的区域在流动方向上的面和流体入口面的边网格化，设置合适的网格数量。得到整体网格如下图 4-3 图 4-3 为了减少入口效应和出口效应，在流体的进口处和出口处分别建立延长段。进口延长段长度为进口水力直径的1.5倍，出口延长段长度为出口水力直径的5倍，如上图。4.4 解的设置4.4.1 边界条件 此模型考虑到计算资源有限，采用压力基隐式求解方式。因此将流体近似认为不可压缩流体，密度取为常数。因为速度入口条件应用于不可压缩流体，且比质量流量入口条件稳定，因此流体入口定义为速度入口边界条件。入口流体的温度保持为300K，以考察各模型对流体的加热情况。速度方向设置为垂直于边界（Normal to Boundary），所以忽略了导流片对流动分配的影响，近似认为入口流动是均匀的。模型一预期取得十个实验数据点，雷诺数从5009500以1000递增。模型二取雷诺数500和9500两个数据点与模型一对比。模型一速度入口的速度值与雷诺数对应为：速度m/s 雷诺数1.45004.2 15006.95 25009.7350012.5 450015.3 550018 650020.87 750023.65 850026.4 9500模型二速度与雷诺数对应为：速度m/s 雷诺数2.2 50039.36 9500 对于雷诺数由25009500的情形，虽然实际上流体处于过度流动阶段。然而在模拟中使用湍流模型，调整湍流强度来模拟其流动状况。界面上的湍流强度定义为： (4-21)其中是速度脉动量，是平均速度。 湍流强度小于1%时，可以认为湍流强度是比较低的，而在湍流强度大于10%时，则可以认为湍流强度是比较高的。内流问题的湍流强度取决于上游的流动状态。如果上游是没有充分发展的未受扰动的流动，则进口处可以使用低湍流强度。如果上流是充分发展的湍流，则进口处湍流强度可以达到几个百分点。鉴于此，考虑到本模型的实际情况，选用较小的湍流强度。其数值在1%到10%浮动，以水力直径为特征长度。 压力出口边界条件指定出口边界处的静压，而静压值仅在流场为亚音速时使用。如果在出口边界上流场达到超音速，则边界上的压强将从流场内部通过插值得到。其他流场变量均从流场内部插值获得。定义压力出口的静压值均为1000Pa（表压），无回流，回流温度均为零。 定义板片的上下两个面为定热流密度（热通量）的壁面，热流密度均设置为。定义模型的左右各个面为对称边界条件。其余的面均定义为绝热壁面。4.4.2 物理材质的物性参数定义 对于Fluent中的物理材质，还需根据本模型的具体情况进行调整。 空气的密度设置为常数，其数值为。定压比热容以温度的多项式计算，设置为式（2-24）。导热系数以运动理论计算，设置为式（2-22）。运动粘度系数以Sutherland三系数公式计算，系数取默认值。 固体材质为铝，其密度、比热容和导热系数均认为常数。其值分别设置为。4.4.3 湍流模型及壁面函数定义 本文模型选用Standard湍流模型，标准的模型通过解两个独立的输运方程来求取湍动粘度第五章 计算及结果分析jc posted 2010年6月01日 19:09 in 毕业论文（不再更新） , 249 阅读 5.1 模型计算调整模型在上一章的设置后，基本可以进行迭代求解了。只是在迭代之前，还需对影响收敛的一些参数进行调整，以加速收敛的速度。 在控制方程以给定的方法离散化后，还必须对离散方程进行某种调整，并且对各未知量（速度、压力、温度）的求解顺序以及方式进行特殊处理。这些需要依赖流场计算方法（SIMPLE和SIMPLEC）来进行调整。 SIMPLE是英文Semi-Implicit Method for Pressure-Linked Equation的缩写，意为”求解压力耦合方程组的半隐式方法“。SIMPLE算法的基本思想可描述如下：对于给定的压力场（它可以是假定的，或是上次迭代计算所得的结果），求解离散形式的动量方程，得出速度场。因为压力场是假定的或不精确的，这样，由此得到的速度场一般不满足连续方程。因此，必须对给定的压力场加以修正。修正的原则是：与修正后的压力场相对应的速度场能满足这一迭代层次上的连续方程。据此原则，我们把由动量方程的离散形式所规定的压力与速度的关系代入连续方程的离散形式，从而得到压力修正方程，由压力修正方程得出压力修正值。接着，根据修正后的压力场，求得新的速度场。然后检查速度场是否收敛。若不收敛，用修正后的压力值作为给定的压力场，开始下一层的计算，如此反复，直到获得收敛的解。 SIMPLEC是SIMPLE的改进算法之一，由Van Doormal和Raithby提出。其稳定性较好，在计算中可以将亚松驰因子适当放大，特别是在层流计算时，如果没有在计算中使用辐射模型等辅助方程，用SIMPLEC可以大大加速计算速度。因此，本模型在层流计算时选用SIMPLEC算法，在湍流时选用SIMPLE算法。 因为Fluent中解的方程为非线性，所以有必要对变量的变化进行控制。通常的方法便是调整亚松驰因子，以减小在迭代过程中的变化。一般情况下，一个单元（cell）中的新值依赖于上一个旧值。亚松驰因子因此定义为： (5-1)由此公式可见，调节亚松驰因子直接影响着计算收敛的速度。在计算过程中，先以默认亚松驰因子迭代，当残差曲线下降过于缓慢时可适当减小亚松驰因子以加速收敛速度。 初始化时迭代计算前最重要的一步，初始化数值是否恰当是影响计算收敛的最重要因素之一。本模型经过不断的摸索尝试，总结出恰当的初始化方式和数值。全局初始化以压力出口为参考，压力值设为1000Pa。待全局初始化完毕，再对其进行适当的修改。在Path面板中，设置Fluid zone的压力为1000Pa，Z方向的速度值为速度入口的值，温度为300K。 收敛条件以默认值为准，即残差值分别为：Continuity 1e-3；X-Velocity 1e-3；Y-Velocity 1e-3；Z-Velocity 1e-3；k 1e-3；epsilon 1e-3；Energy 1e-6。绘制残差曲线以监视计算收敛情况。 另外，在进行迭代计算前还需建立等值面（ISO-Surface）用以生成等值线图和矢量图。以建立Surface Zone，即板翅式换热器模型在其高度一半处沿流动方向的一个分割面。 计算残差图如下：图 5-1 模型一在Re=500时的残差图图 5-2 模型一在Re=9500时的残差图5.2 计算结果分析5.2.1 换热分析当流体从空间进入管道时，流动的边界层会有一个从零到汇合于流道中心线的过程。相同的，热边界层也如此。当热边界层汇合于中心线后，换热已经充分发展，此后的换热强度将不再保持不变。如图5-3所示。图5-3（a）Re=500时温度分布图图5-3（b）Re=9500时温度分布入口段的热边界层较薄，局部换热系数比充分发展段的高，且沿流动方向逐渐降低