多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究.pdf

应用数学和力学编委会，：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究何 树，娄 钦（上海理工大学 能源与动力工程学院，上海；上海理工大学 上海市动力工程多相流动与传热重点实验室，上海）摘要：采用介观相变格子（，）方法，在孔隙尺度下研究了多孔介质的孔隙率对池沸腾换热过程的影响，重点分析了不同孔隙率时气泡的运动过程，并对气泡在多孔介质中的典型状态进行了力平衡分析，进而探究了多孔介质孔隙率影响沸腾传热的机理 结果表明，与无多孔介质的平板表面相比，多孔材料能够有效地降低初始成核的壁面过热度，增强流体的扰动，并且能够显著提升临界热流密度（，）值 在所研究的工况中，孔隙率 时，值提升最大，约为平板的 倍，其余孔隙率的多孔介质最小也可将其 值提升至平板的 倍 研究发现，当孔隙率从 开始逐渐减小时，值逐渐增大，同时沸腾换热曲线向左上方移动，这是因为减小孔隙率能够增大有效换热面积，减小气泡成核的壁面过热度，从而强化沸腾换热 当孔隙率减小到 时，若继续减小孔隙率，热流密度将突然下降，沸腾传热性能显著降低 通过对沸腾过程中气泡的受力进行分析后发现，当孔隙率较小时，过小的孔隙直径显著增大了气泡的逸出阻力，降低了气泡的上升速度，延长了气泡脱离多孔介质的时间，且此时气泡会在蒸发动量力、接触压力以及摩擦力等的共同作用下聚集在加热器上表面，形成气膜，从而恶化沸腾传热 关 键 词：多孔介质；池沸腾换热；气泡动力学；孔隙率中图分类号：文献标志码：，（，；，）：，应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金（）；上海市浦江计划（）作者简介：何树（），男，硕士生（：）；娄钦（），女，教授，博士，博士生导师（通讯作者：）引用格式：何树，娄钦 多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究 应用数学和力学，（）：，（），；，：；引 言池沸腾作为一种高效的换热方式，已经被广泛应用于航空航天器、电子元器件以及电厂等各个领域 对于沸腾换热过程来说，、换热系数以及初始壁面过热度等指标是衡量沸腾换热性能好坏的关键，因此如何提高、有效降低初始壁面过热度，增大换热系数，对优化沸腾换热性能至关重要 近几十年来，为了提高沸腾换热性能，学者们采取了各种措施，包括对加热器表面结构改造，在工作流体中添加活性剂以及在加热器表面施加磁场或者电场进行强化等 多孔介质因具有较大的换热表面积、较高的有效导热系数等优势，成为相变传热领域最具发展前景的冷却技术之一，一些学者采用多种方法研究了多孔材料在不同条件下改善沸腾换热的潜力 作为人们认识复杂物理现象的重要方法，实验方法被用于含有多孔介质的池沸腾换热问题的研究 和 采用实验方法研究了 介质液体在数十至数百微米孔径的多孔石墨中的沸腾过程，并与相同尺寸的光滑铜进行了对比 他们发现在沸腾开始之前，由多孔石墨释放的微小气泡增强了自然对流换热，同时，多孔石墨的存在增加了成核点位，强化了初始核态沸腾 同样采用实验方法，等研究了去离子水在更小直径（孔径大小为 ）多孔介质中的池沸腾换热 他们发现，纳米多孔表面的初始壁面过热度比非涂层的表面要小，而且纳米多孔涂层表面的核态沸腾换热速率比未涂层的表面大 相似地，还有一些学者也发现了不论是微米级孔径还是纳米级孔径，多孔介质的存在都能够有效地降低初始壁面过热度，从而强化沸腾换热 此外，也有一些学者发现，一些特殊的孔隙结构也可以提高沸腾换热性能 和采用实验方法研究了饱和水在多孔蜂窝结构板上的池沸腾过程，发现多孔蜂窝板的毛细作用实现了自动供液，蜂窝状结构使液汽分离，降低了蒸汽逸出的流动阻力，最终使 增大至 ，约为普通表面的 倍 另一方面，等用实验研究了莲花型多孔铜板中饱和水的沸腾传热性能，他们发现采用莲花铜板加热表面，可以增加到 ，几乎是光滑表面情况下的两倍 尽管 等和 等采用的孔隙结构不同，但这些孔隙结构提高 的原理都在于增强毛细作用下的自动供液能力，增加有效成核点位以及增大传热面积 此外，还有学者探究了多孔介质的结构参数对池沸腾换热的影响，等研究了非均匀和均匀多孔涂层表面的池沸腾传热，他们的研究结果表明，周期性地改变颗粒基底厚度的多孔涂层表面可以提供更高的 值（最大可达到普通表面的 倍）等研究了多孔材料的导热率和表面润湿性对沸腾换热的影响，发现多孔材料导热率对沸腾传热的影响比润湿性的影响更大 除了实验方面，一些学者对多孔介质在不同工作流体中的沸腾过程进行了数值模拟研究 和 第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究基于两相混合模型对泡沫铝在水中的流动沸腾过程进行了数值模拟，发现在沸腾过程中，气泡通常成核于金属泡沫韧带和加热器表面接触的角点处，同时也得到了该过程的沸腾曲线，其数值结果与他们的实验结果吻合较好 等对分离相模型进行了改进，研究了多孔介质中冷却剂的相变，同时考虑了液相区和气相区的对流效应 前期的大多数模拟工作采用的是体积平均法 而近年来，在众多数值方法中，格子（）方法这种介观方法在两相流建模方面表现出了巨大的潜力，特别是对于界面复杂的拓扑变化（如气泡变形、聚结、破裂等）或复杂几何形状的实现等，可以认为是孔隙尺度建模的良好方法 而在多孔介质沸腾领域，等采用 方法研究了方形截面实心柱阵列内池沸换热的增强作用 其结果表明，多孔阵列的存在使得受热面在较早的时间步长上有更多的成核位点，从而增加了热流密度（值约为平板表面的 倍），除此之外，他们对比了两种高孔隙率（，）的多孔阵列后发现，孔隙率对传热速率的敏感性显著，在表面过热度较高时尤为明显，并预测了存在最佳孔隙率，但并没有对最佳孔隙率进行量化研究 等采用多松弛 方法模拟了亲水金属泡沫表面的池沸腾传热过程，研究了金属泡沫的孔隙率对沸腾传热的影响，但其研究的孔隙率仅限在 范围内，并没有考虑更小的孔隙率，且没有得到最优的孔隙率结果 除了对孔隙率进行研究之外，等也采用 方法模拟了一种二维多孔结构在水中的沸腾传热过程，在孔隙率一定的情况下得到了泡沫厚度对沸腾传热的影响 以上研究表明，多孔介质对于池沸腾换热的强化作用显著 但由于其内部存在特殊的拓扑结构，因此，在沸腾过程中多孔介质孔隙内的气泡动态行为十分复杂，而目前对于该过程的数值研究仍较少，且对池沸腾换热的影响机理研究尚待完备，对于多孔介质孔隙率的研究还缺乏一些具体的影响机理以及定量的分析 因此，有必要对多孔介质的池沸腾过程及其孔隙率的影响机理进行更详细的孔隙尺度研究 基于此，本文构造了一种二维均匀排列的多孔介质，采用伪势 模型在孔隙尺度下模拟了多孔介质的池沸腾过程，详细分析了不同过热度下，高、中、低孔隙率的多孔介质中气泡的动力学行为特性，并对沸腾过程中典型的气泡状态进行了力平衡分析，最后深入探究了多孔介质的孔隙率对沸腾换热曲线的影响 相变 模型在 模型中，流体密度的分布函数演化方程为（，）（，）（，）（，），（）其中（，）是 时刻 位置的分布函数，是分布函数的方向，速度在 方向上的离散为，是离散速度的个数，为无量纲松弛时间，它与运动黏度 的关系为 （），是一个常数，其中 ，为格子空间步长，为时间步长 在式（）中，平衡态分布函数（，）可由速度 和密度 计算得到：（），（）其中 是权重系数，在 模型中速度的 个方向由下式给出：（，），（），（），（），（），（）权重系数分别为 ；，；，在方程（）中，力项 由下式给出：（，），）（，），），（）其中 为由合力 计算出的速度，合力 定义为 ，（）其中 为粒子间的相互作用力，为合外力，一般为重力或者浮力，为表征流体与固体壁面润湿性的应 用 数 学 和 力 学 年 第 卷力 定义为（）（）（），（）定义为（）（）（），（）其中 为调节接触角的流固相互作用强度；（）为指标函数，当为固体格点时，（）为，当为流体格点时，（）为；（）为有效质量，它与局部密度和局部压力有关，（）（），其中为相互作用强度，局部压力 可以由状态方程得到 在本研究中采用（）状态方程，它具有相当高的稳定性和准确性，其由下式给出：（），（）其中 为气体常数，为流体温度，和分别为修正压力和修正温度 据研究表明，能够以可接受的精度再现饱和水和蒸汽，外力项 由下式计算得到：（）（），（）为重力加速度，为每个计算步下整个计算域的平均密度 温度的演化方程为（，）（，）（，）（，），（）其中（，）为温度分布函数，为温度的无量纲松弛时间，它由热扩散系数 计算得到，（）温度的平衡态分布函数定义为 （），（）其中 为宏观速度 式（）中的源项 由下式给出：，（）为比热容 相关宏观量的计算如下：，（），（），（）（）由于固体加热器以及多孔介质与流体的热物性不同，导致流固交界面温度变化不连续 为了避免这种不连续性，采用 等提出的二阶精确共轭传热边界条件 在该方法中，沿格子速度方向指向流体的加热器或者多孔介质，以及沿格子速度方向指向固体的加热器或者多孔介质周围的流体节点 的温度分布函数分别设为（，）（，）（，），（）（，）（，）（，），（）其中 为碰撞后的分布函数，为指向流体的方向，为指向固体的方向，为固体与流体的热质量比，（）（）注意，在液气界面上，流体性质（如热扩散系数和黏度）由下式给出：第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究 ，（）其中 和 分别为饱和液体和饱和蒸汽的密度 物理模型和计算参数由于均匀多孔介质结构内部流体输运机理更容易刻画，同时也因为随着 打印技术的发展，越来越多的学者将规则多孔结构应用在沸腾过程中以达到强化换热的目的 因此，本文设计了一种简化的均匀二维排列的多孔介质进行研究，不仅可以便于分析多孔介质内的沸腾换热机理，也可以更方便地推广到实际应用中 多孔介质池沸腾的模拟实验如图（）所示，整个计算区域设置为 （格子单位），其中红色部分为加热器，高 ，长 ，加热器上方放置均匀多孔介质 初始时，计算区域内充满温度为，密度为 的静止饱和液体 下边界（）的流体部分采用绝热无滑移边界条件，计算域的两侧为周期边界条件，计算域顶部边界为对流边界条件 加热器底部加热方式采用等温加热，其加热温度为，加热器的左右边界采用绝热无滑移边界条件，加热器上方速度无滑移，温度分布采用流固耦合边界条件 多孔介质内 方向相邻两层多孔介质块之间的传热方式为导热，其处理方式如图（）所示 如图（）所示的任意两个在 方向相邻的多孔介质块，其中、示意部分分别是下方多孔介质临近流体侧的固体部分、上方多孔介质临近流体侧的固体部分、最靠近下方多孔介质的流体层，以及最靠近上方多孔介质的流体层 为了实现多孔介质之间的导热，将下方多孔介质 层固体格点的温度分布赋值给上方多孔介质下最临近的流体格点 层；同时将上方多孔介质 层固体格点的温度分布赋值给 层流体格点 除此之外，层和 层流体格点将与其周围的流体格点保持绝热（即流体中画叉的部分与周围画圆的部分保持绝热），而多孔介质则分别与 层和 层流体进行耦合传热，最终实现多孔介质间的连续热传导（）物理问题示意图（）多孔介质间的导热示意图（）（）图 物理问题以及多孔介质间的导热示意图 注 为了解释图中的颜色，读者可以参考本文的电子网页版本，后同 本文中所采用的模拟参数如下：，加热器与多孔介质内部初始温度设置为，重力加速度 取（，），比热容，动力黏度 ，热扩散系数 ，加热台的热物性参数取为（），而多孔介质的热物性参数取为（），其中，根据其有效导热率的变化而变化 应 用 数 学 和 力 学 年 第 卷下文中所采用的量均为无量纲量，其对应的特征长度、特征速度以及特征时间如下：（），（）则对应的无量纲长度以及无量纲时间为 ，本文所用的固体材料为铜，流体为，模拟工况设计和参数的选取是依据实际固体和流体的性质，相关参数格子单位与物理单位的转换关系详见表 无量纲壁面过热度 以及加热器底部（）的局部温度梯度 ，计算如下：，（），（）其中 为流体的潜热，本模拟中取 ，为格子空间步长 无量纲空间平均热流密度 计算如下：，（）无量纲时间空间平均热流定义为 （）（）表 格子单位与物理单位转换 ，（）（）（）（）模 拟 结 果 不同过热度下多孔介质中气泡的动力学行为本小节研究了不同过热度下气泡在多孔介质中的成核、生长、滑动、上升、到达顶部以及逸出多孔介质等动态特性 为了对比气泡行为在不同沸腾模态下的差异，我们模拟了加热温度分别为，以及 的沸腾过程，研究对象为 个不同孔隙率的多孔介质，它们的孔隙率分别为，以及，除孔隙率外，它们的其余参数设置均相同，具体为孔隙密度 ，高度 ，有效导热率，表面接触角设置为 图 为 的多孔介质在 种加热温度下，池沸腾过程中气泡的运动状态图，其中图（）、（）和（）的加热温度分别为，以及 从图（）中可以发现：在 时刻，加热面中心有两个气泡成核 随着沸腾的持续进行，这两个气泡逐渐生长，在 时刻，由于多孔介质的阻碍发生形变（如图中的红色虚线框所示），随后两个气泡在浮力的作用下脱离加热器表面，此时加热器表面残余了两个气泡 初始的两个气泡在 时刻穿过多孔介质顶部并逸出，上一时刻残余在加热器表面的气泡则开始生长，随后会参与下一阶段的气泡逸出循环过程，显然此时为典型的核态沸腾状态 如图（）所示，随着加热温度提高到，可以发现：多孔介质中气泡的成核位置明显增加，成核时间也随之提前，在 时刻，加热器上表面的 个间隔内均有气泡生成 随着气泡的生长，它们遇到了多孔介质的阻碍，这与低过热度的情况一致，但不同的是，由于过热度的增加，相变更加剧烈，生成的气泡质量明显增加（）随后，这些气泡在浮力的作用下向上升，并到达多孔介质顶部，两侧的气泡率先逸出（），第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究中心的气泡则在 时刻才与下方的气泡分离，并且在分离的过程中由于接触线的拉长，会产生一些小气泡随之分离 继续增大加热温度至，从图（）中可以发现：由于底部加热温度的增加，成核时间进一步提前（），而在第一轮相变产生的初始气泡逸出后（），新一轮的气泡聚集在加热器上表面，形成了一层气膜（），显然此时已到达稳定膜态沸腾状态（）（）（）图 孔隙率为 的多孔介质在加热温度 ，下的各个时刻气泡状态图 ，如图 所示，当多孔介质的孔隙率减小至 时，不同过热度下的气泡状态也发生了变化 图（）的加热温度为，可以发现：在此加热温度下，气泡的成核和脱离状态与 的多孔介质大致相同，均有两个气泡在中心成核，随后在生长过程中受到多孔介质的阻碍，最后在浮力的作用下逸出多孔介质 但有所不同的是，的多孔介质其气泡成核时间略晚（），分别统计了 和 多孔介质内的气泡从成核到与加热器表面脱离的时间段内生成的气泡质量（即图（）中的 至 时刻和图（）中的 至 时刻），发现在相同时间内（），的气泡增加了（格子单位），而 的气泡只增加了（格子单位），这说明孔隙率的减小虽然会使气泡成核时间相对延后，但沸腾的相变速率却会增大 这与 等的研究结果一致，他们发现孔隙率越小，气泡越容易合并，骨架对气泡离开的阻碍作用越强，且在孔隙率较小的金属泡沫中，蒸汽成核体积较大，相变速率增加，这是传热面积增大和气泡释放阻力增大的相互作用结果 如图（）和（）所示，继续增加加热温度至 和 由图（）可以发现：壁面过热度的增加使气泡成核点明显增加，气泡质量也随之增大，但成核与脱离的时间与 的多孔介质相比略晚（和）而当加热温度提升至 后，如图（）的后 幅图所示：在图中的红色虚线框内，气泡脱离之后加热器表面仍然会残余一层气膜，这说明在此过热度下沸腾同样达到了膜态 应 用 数 学 和 力 学 年 第 卷如图 所示，进一步减小孔隙率至 由图（）可以发现，在加热温度为 时，成核时间也相对延后（），同时可以发现生成的气泡在生长过程中更容易聚集在一起（），且气泡在 脱离时尾迹会拖得更长 我们同样统计了气泡从成核到与加热器表面脱离的时间段内生成的气泡质量（即图（）中的 至 时刻，），发现气泡质量增加了 （格子单位），与 的多孔介质相比又有所增加 由图（）可以发现：提高加热温度至 后，气泡在生长过程中很容易在加热台表面聚集形成一层气膜（），且在气泡脱离后残余的气泡仍然会聚集在一起，这应该是由于孔隙直径过小导致的，而这种气泡状态也必然会恶化沸腾传热 继续提高加热温度至，如图（）所示：随着壁面过热度的增加，气泡的聚集行为更加明显（），且此时的气泡阻力过大，不易脱离（）（）（）（）图 孔隙率为 的多孔介质在加热温度 ，下的各个时刻气泡状态图 ，以上研究完整地展现了含多孔介质的加热器的池沸腾从初始核态到完全模态的过程 除此之外，对比 种孔隙率多孔介质的沸腾过程后发现：减小多孔介质的孔隙率能够有效地增加沸腾的相变速率，强化沸腾换热；然而，气泡的成核时间会相对延后，并且减小孔隙率还会使气泡更易聚集，形成气膜，气泡上升阻力会增大，其逸出多孔介质的速率也会随之降低 需要注意的是，从 个沸腾状态的图中均会发现气泡成核于多孔介质间隔内 而这与 等所得到的，气泡通常成核于金属泡沫韧带紧密附着在加热器表面的节点处的结论略有不同，这是因为气泡成核的位置受多孔介质的结构参数（如孔隙密度、有效导热率以及表面润湿性第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究等）影响，我们发现：在一些参数下，气泡成核发生在方肋格子角点处；而在另一些情况下，气泡会成核在方肋格子之间 如对于本文所研究的工况，当接触角为，孔隙率为，加热温度为 时，气泡首先在方肋格子角点处成核（如图 所示）；而本文中图 对应的参数，此时是在方肋格子之间成核，这一结论与 等的模拟结果一致（）（）（）图 孔隙率为 的多孔介质在加热温度 ，下的各个时刻气泡状态图 ，图 孔隙率为、接触角为 的多孔介质在加热温度 下的气泡状态图 应 用 数 学 和 力 学 年 第 卷 沸腾过程中单气泡在多孔介质中的力平衡本小节对沸腾过程中单气泡在多孔介质中的运动行为进行了力平衡分析 根据上一小节的分析，由于多孔介质的存在，气泡的运动行为主要包含成核与穿过多孔介质这两种情况，这两种典型的气泡状态显示在图 中（其中浅蓝色形状代表气泡，红色方格代表二维多孔介质骨架，剩余的箭头代表沸腾过程中作用于气泡上的力）总体来说，如图（）所示，沸腾过程中作用于气泡成核生长的力一般包括重力、浮力、蒸发动量力、接触压力 以及表面张力 除此之外，如图（）所示，当气泡通过多孔介质时，由于气泡在多孔介质表面滑动，故它们之间会产生摩擦阻力 因此，池沸腾过程中多孔介质内的气泡在 方向和 方向的力平衡方程可由如下两式表示：，（）（）（）单气泡成核生长（）单气泡穿过多孔介质（）（）图 沸腾过程中单气泡在多孔介质中的力平衡 表面张力 的定义为 ，（）其中 为气泡的接触直径，为表面张力系数，为气泡与固体的接触角 接触压力 定义为，（）其中 为气泡的等效半径，接触压力的方向垂直于接触壁面，接触直径 可由下式计算：，（）其中 为气泡的基线周长 蒸发动量力 定义为，（）其中 为蒸发热流密度，为界面处的蒸发潜热，为界面面积 当气泡的受力情况满足式（）和（）时，气泡在当前位置成核生长，直至气泡体积达到一定大小时，浮力增大，方向上的平衡被打破，气泡脱离加热面 气泡继续上升后会遇到多孔介质的阻碍，如图（）所示，此时影响气泡上升的阻力主要为气泡与多孔介质间的接触压力 与摩擦力 摩擦力定义为 ，其中 为摩擦因子，摩擦力的定义参考前人的研究，此外，根据 等的研究，摩擦因数 的大小和气泡与固体壁面间微小流体层的厚度相关，而该厚度则与工作流体的黏度、表面张力以及气泡的运动速度相关 因此，接触压力的大小在很大程度上第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究影响着气泡上升的速度，这与 等的实验研究中对接触压力的观点一致 孔隙率对多孔介质池沸腾曲线的影响本小节研究了孔隙率对多孔介质池沸腾换热曲线的影响 根据 小节的研究，继续设置另外 个不同孔隙率的样本进行对比研究，总共 个样本，它们的孔隙率由高到低分别设置为 ，其余参数与 小节一致 图 为不同孔隙率下的多孔介质样本以及不含多孔介质的光滑平板的沸腾曲线 从图中可以观察到众多孔介质样本的沸腾曲线均在光滑平板曲线的左上方，这说明加上多孔介质之后，既相当程度地提高了 值（平均提高了约 倍），也降低了成核的初始壁面过热度（降低了约），显著增强了沸腾换热 而这一结果也与 和 的实验结果一致，他们对厚度为 的氧化铝纳米多孔表面在饱和 介质液中的池沸腾换热进行了实验研究，发现纳米多孔表面比普通表面的初始壁面过热度降低了 他们将引起这个现象的原因归结为纳米多孔结构的存在增加了蒸汽体积和活性成核点密度，从而增强了沸腾换热性能 图 不同孔隙率的多孔介质以及光滑平板的沸腾曲线 从图 中还可以发现：当孔隙率从 开始减小直至 时，沸腾换热曲线随之向左上方移动 我们随机选取了 的 个样本在核态沸腾区域的气泡脱离状态进行分析 图 显示了不同孔隙率的多孔介质样本在加热温度为，时刻的气泡状态以及其所对应的温度场（该图所显示的气泡均为第一次相变成核的气泡脱离情况）由图 可以发现：和 的两个样本其气泡脱离较快，而随着孔隙率降低到 和，气泡脱离速度逐渐变慢，这与 小节所观察到的现象一致，这是由于孔隙率的减小增大了气泡的脱离阻力 与此同时，观察同一时刻下对应的加热器内部的温度分布却发现：大孔隙率的多孔介质内气泡虽然脱离较快，但是其加热器内部的温度场分布却不如小孔隙率的样本分布均匀，且其上表面（与多孔介质相接触的那一层表面）温度较高；而随着孔隙率的降低，加热器内部的温度场分布愈加趋于均匀，其上表面的温度亦随之降低 这是因为随着孔隙率的降低，换热面积随之增大，热量能够更快地传导 由此可见，在本研究中，当孔隙率 时，随着孔隙率减小，换热面积增大，与此同时，气泡自由上升阻力也随之增大，但此时增大的阻力并不足以让气泡聚集在加热面上形成气膜，对沸腾传热的恶化作用有限 因此当孔隙率 时，由于减小孔隙率而增大换热表面积的强化传热作用大于增加气泡上升阻力的恶化传热作用 在 等的金属泡沫表面池沸腾换热与槽阵泡沫表面气泡团簇生长行为的实验中同样发现，在孔密度一定的情况下，孔隙率的大小影响着气泡与金属泡沫间的传热面积以及气泡的上升阻力，从而影响沸腾传热性能，这三者之间存在相互耦合的关系，这与本文的结论一致 应 用 数 学 和 力 学 年 第 卷（）（）（）（）图 不同孔隙率的多孔介质在加热温度为 下的气泡状态及加热器温度场分布（）（）如图 所示，当 时，继续减小孔隙率却发现沸腾传热发生了恶化 为了进一步分析其原因，随机选取了两个加热温度下的某两个时刻，对比分析了 和 两个样本的气泡状态 图 显示了这两个样本在同一加热温度下的气泡行为 横向对比图（）和（）可以发现：当加热温度 时，分别在 和 这两个时刻下，样本沸腾过程中的气泡均聚集在一起，在加热器表面形成了一层气膜（如图中的红色虚线框所示）；而当加热温度增加至 时，在 以及 这两个时刻，样本的气泡也聚集在了一起（），（），（），（），图 两种加热温度下 和 的样本其沸腾过程中气泡的聚集状态 第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究反观图（）和（），的样本在同一时刻下则并没有聚集形成气膜 通过随机选择的两个加热温度分析发现：的样本与 的样本相比，在沸腾过程中其气泡都更容易聚合在一起，除此之外，我们在 和 这两个样本中也观察到同样的现象，小节已经对 的样本进行了研究，而对于 的样本这里也不再赘述 以上研究表明：多孔介质的孔隙率存在一个临界最优值，即 ，使得沸腾传热性能达到最优，一些学者也通过实验研究得到了多孔介质的孔隙率过高或者过低都不利于沸腾换热性能提高的结论 如 等用实验研究了水平导电微孔包覆表面池沸腾过程，发现多孔涂层的孔隙率过高或者过低都不利于提高 值；等用实验对烧结铝粉芯多孔介质内的池沸腾换热性能进行了研究，同样发现烧结多孔芯的孔隙率过大或者过小都不利于传热 此外，如图 所示，我们将模拟所得到的不同孔隙率的 结果与 等的理论预测结果进行了对比 可以明显看出，等的理论预测结果与本文的模拟结果趋势一致，即存在一个最佳孔隙率使得 值达到最优 当然，从图中还可发现，本文的计算结果与 等的预测结果存在一定的差异，这是因为多孔介质内部的结构以及流体或固体的物理性质（如黏度、表面润湿性、导热系数等）所造成的 图 本研究得到的不同孔隙率中的 值与 等的理论预测结果对比 为了探究孔隙率小于的样本内的气泡更容易聚集的原因，以 的样本为例对其气泡上升阻力进行了分析 如图（）所示：加热温度 ，在 时，正中的两个气泡其自由生长遇到了多孔介质的阻碍，未能及时地脱离，而随着体积的不断增大，蒸发动量力 的横向 分量使得气泡横向扩张 与此同时，由于中心气泡的生长引起了流场的扰动，两侧气泡在流场的扰动作用以及 的横向 分量的共同作用下也向中心聚合，总体的运动结果如图（）所示，在时间到达 时，这 个气泡聚集在加热面上形成一层气膜（）气泡聚集（）形成气膜（）（）图 气泡在 样本中的聚集过程的受力分析 应 用 数 学 和 力 学 年 第 卷表 统计了 个孔隙率样本的最大接触压力大小 以及气泡脱离多孔介质前的平均上升速度 可以发现：由于孔隙率的降低，在固定高度以及孔隙密度的情况下，小孔隙率的多孔介质其孔隙直径会随之减小，气泡的接触压力项 中的当量半径变化不大，而气泡的接触直径却随之增大 根据式（）可知，最终导致气泡的接触压力 增加，气泡上升阻力增大，上升速度随之降低 在 等所做的，关于水和庚醇溶液在池沸条件下的梯度金属泡沫池沸传热性能的实验研究中也发现，金属泡沫的孔隙减小，有限的空间会阻碍气泡的膨胀与上升 表 不同孔隙率的多孔介质样本其气泡最大接触压力 以及平均上升速度（格子单位）（）根据本节的研究可以发现：在其余参数设置均相同的情况下，随着孔隙率的减小，换热面积也会增大，但孔径会减小，气泡脱离阻力会逐渐增大，且当多孔介质的孔隙率减小到一定程度后，阻力过大，气泡难以逸出多孔介质，而更多的气泡在受到阻碍后会聚集在一起，最终形成一层气膜覆盖在加热器表面，导致周围的液体不能及时的补充，使沸腾换热恶化 结 论本文采用 方法研究了均匀多孔介质的孔隙率对池沸腾传热性能的影响 首先，以孔隙率 ，的均匀多孔介质为例分析了其在不同过热度下的气泡运动行为，并对气泡成核和受到多孔介质阻碍这两种典型状态进行了力平衡分析 研究发现：减小孔隙率能够增大相变速率，但同时也会增大气泡脱离阻力，使气泡易聚集成气膜，并延后成核与脱离的时间 随后，详细研究了不同孔隙率对池沸腾换热性能的影响机理，重点分析了沸腾过程中的气泡行为，并给出了不同孔隙率下池沸腾从核态沸腾到模态沸腾的热流密度曲线，并且将模拟结果与相关沸腾实验进行了印证 研究发现：虽然多孔介质能够增强流体的扰动，增加换热面积与成核点位，强化导热，且能够降低壁面的初始过热度（平均约）并显著提高 值（平均 倍），但是其提升的性能大小与多孔材料的孔隙率有关，具体研究结论如下：）有多孔介质存在可以降低气泡成核的初始壁面过热度，并提高 与平板相比，有多孔介质存在时气泡成核的初始壁面过热度降低了，同时 最大可提高 倍）当多孔介质的孔隙率 时，减小孔隙率能够增大换热面积，增加相变速率，降低初始壁面过热度，从而强化沸腾换热）当 时，由于过小的孔隙增加了沸腾过程的气泡逸出多孔介质的阻力，延长了气泡在多孔介质内的停留时间，最终在流场与蒸发动量力 的共同作用下，多个气泡更易聚集在加热器表面形成气膜，这不利于沸腾换热的持续进行，最终导致沸腾传热恶化）通过与多孔介质内的气泡受力分析发现当孔隙率过小时，气泡的接触压力 逐渐增大，导致气泡上升阻力增大，最终使得气泡上升速度降低 本文的研究结果表明：多孔介质对于池沸腾换热的强化作用显著，但其孔隙率存在最优值，使得沸腾换热性能提升最大，在实际实验或者应用过程中，针对固定的结构参数，仅改变多孔介质的孔隙率的情况下，可以优先在高过热度下对孔隙率由低到高的多孔介质进行沸腾传热实验测试，若观察到沸腾过程中的气泡从第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究易聚集状态变为自由脱离状态时，即可确定该结构的最佳孔隙率范围 参考文献（）：，（）：李迎雪，王浩原，娄钦 含多个矩形加热器通道内流动沸腾传热性能的介观数值方法研究 应用数学和力学，（）：（，（）：（），：，：，：，：，：，：，（）：，（）：，（）：，（）：，（）：，（）：，：，（）：，：，：，：应 用 数 学 和 力 学 年 第 卷 ，（）：，：，（）：，（）：，：，：，：，（）：，（）：陆威，王婷婷，徐洪涛，等 多孔介质复合方腔双扩散混合对流 模拟 应用数学和力学，（）：（，（）：（），：，（）：，：，（）：，（）：，（）：，（）：，：，：胡卓焕，罗婷，许佳寅，等 毛细芯蒸汽槽道孔径对环路热管（）传热性能影响研究 热能动力工程，（）：（，（）：（），（）：第 期 何树，等：多孔介质孔隙率对池沸腾传热性能影响机理的模拟研究 ，（）：，（）：，：，：，（）：，（）：，：，：，：，（）：，：，：，（）：应 用 数 学 和 力 学 年 第 卷