凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究.pdf

第 卷第 期 年 月能 源 环 境 保 护 ，卢浩，新疆大学西北能源碳中和教育部工程研究中心副主任，教授，博士生导师。入选中组部国家海外高层次人才计划，新疆杰出青年科学基金获得者，新疆大学双一流建设“能源碳中和”创新团队负责人，新疆大学“新能源”研究中心首席专家。卢浩教授在国内外学术期刊及会议上发表论文 余篇，其中本领域国际知名期刊 、等发表 论文 余篇（第一 通讯作者 余篇，国际 一区论文 篇），累计 影响因子总和达到，论文引用次数 余次，个人 因子为，申请 授权国家发明专利 余项，参与撰写学术专著 部。卢浩教授作为项目负责人或主研人员先后主持或参与了国家重点研发计划项目、国家自然科学基金项目、国家社科基金重大项目、新疆杰出青年科学基金项目、新疆重大科技专项等 余项课题。移动扫码阅读韩尊师，卢浩，赵文君，等 凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究 能源环境保护，（）：，（）：收稿日期：基金项目：国家海外高层次人才计划、国家自然科学基金（）；国家社会科学基金重大项目（）作者简介：韩尊师（），男，山东济南人，硕士在读，主要研究方向为能源多相流及传热传质。：通讯作者：卢 浩（），男，湖北荆州人，教授，主要研究方向为能源多相流与传热传质。：凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究韩尊师，卢 浩，赵文君，王泽宇，肖志博（新疆大学 电气工程学院 能源碳中和实验室，新疆 乌鲁木齐；新疆大学 西北能源碳中和教育部工程研究中心，新疆 乌鲁木齐；新疆大学 未来技术学院 新能源研究中心，新疆 乌鲁木齐）摘要：颗粒物在换热通道内的沉积会显著降低换热效率。为了提高换热器的性能，更好地理解颗粒物在换热通道中的沉积特性对于提高换热器性能至关重要。利用雷诺应力模型（）和离散颗粒模型（）对三维凹凸酒窝通道中颗粒的沉积特性进行了模拟。在网格无关性测试和数值验证之后，对凹凸酒窝通道中的流动结构、二次流、温度分布和颗粒沉积分布进行了分析和研究。研究结果表明，影响颗粒在凹凸酒窝通道中沉积特性的主要机制是重力、热泳力、湍流结构、二次流和颗粒惯性，这些机制的影响随颗粒大小的变化而显著改变。对于小颗粒（），凹凸酒窝结构的存在会促进沉积；但对于大颗粒（），凹凸酒窝结构的存在只会对颗粒沉积位置的分布产生很小的影响。关键词：颗粒沉积；凹凸酒窝结构；离散颗粒模型；数值模拟中图分类号：文献标识码：文章编号：（），（，；，；，）能 源 环 境 保 护第 卷第 期：，（）（），（），；，（），：；引 言热交换器广泛应用于能源利用领域，包括加热和冷却系统在内的各种设备都需要在气流和热交换器表面之间进行热传递。然而，在长期使用过程中，湍流携带的颗粒会沉积在热交换器表面。这将导致传热效率的大大降低，同时还会增加热交换器的压降。因此，了解颗粒如何附着在热交换器表面非常重要。研究方法在数值模拟研究中非常重要，研究方法的不同往往会对研究结果产生一定的影响。等对比了拉格朗日技术和欧拉技术在相同运行环境下预报湍流中颗粒沉积的效果。结果表明，欧拉技术的计算时间低于拉格朗日方法。等通过数值和实验研究，帮助理解了紧凑型热交换器的热力学特征，从而进行了有效的设计，并对热交换器管束进行热力、水力和空气动力测试，提出了具有不同翅片比的各种全尺寸管试样。热交换器的结构会影响其传热性能。研究发现，可以通过优化和改进热交换器结构来实现热交换器的优化。等通过数值模拟的方法分别对凹酒窝换热管道和凸酒窝换热管道进行了研究。研究者采用了 方法对两种换热管到进行了仿真，模型被用于流场的计算，结果表明与光滑管道相比凹酒窝和凸酒窝换热管道均可强化管道的换热且可以促进大颗粒的沉积。等对 形弯曲带肋通道中的颗粒沉积进行了离散模拟。气相计算采用雷诺平均纳维斯托克斯技术，结果表明后排肋条受到前排肋条的保护，较高的雷诺数和较大的颗粒直径促进了颗粒在壁上的沉积。本研究的重点是分析热交换通道内结构的改变对通道中的流场和颗粒沉积情况的影响。模拟中研究了一种规律性凹凸酒窝热交换管道的换热和颗粒沉积情况。通过比较管道中有无凹凸酒窝存在以及改变入射颗粒大小，研究不同条件对管道中热交换特性和颗粒沉积特性的影响。数值方法 气相湍流模型雷诺应力模型（）考虑了湍流的各向异性，可以计算出管道中由于应力引起的二次流，是 中提供的最精细的 湍流模型。采用求解雷诺应力输运方程以及耗散率方程而不是涡黏性假设来封闭方程。因此 模型在三维求解中需要额外求解七个方程，气流的质量、动量和能量方程可以描述为：（）（）()()（）（）()（）式中，为流体密度，；为压力，；为焓值，；，分别表示笛卡尔坐标分量、速度分量以及质量力的分量；为分子粘性系数；为韩尊师等 凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究湍流粘性系数，由湍流模型确定。由于粗糙结构的存在，在壁面附近会出现一些复杂的湍流结构。本文为了准确地模拟粗糙流道的湍流流动，采用了 模型来对气相部分进行模拟，模型相较于，这两种模型来说，模型更多地考虑了湍流涡粘度的各向异性，而且它大大提高了计算的精度。湍流雷诺应力模型（）的方程表达式为：()时间变化率()对流项扩散项压力产生项应变项耗散项（）式中，是经验常数，分别为，和。湍流耗散率 的计算公式为：（）其中，。离散相模型颗粒在粗糙流道壁面沉积主要包括流场中的运动、与粗糙壁面的接触粘附或反弹、在流场的出口逃逸这三个过程，颗粒在流场中流动时主要受接触力（碰撞力和接触粘附力）、非接触力（非接触粘附力）、流体对颗粒的作用力（涡团区域的浮力、布朗力、升力和热泳力）以及颗粒自身重力的作用。因为空气密度与粒子密度之比非常小（本研究中为 ），所以与粒子上的外力相比，压力梯度力和虚假质量力都很小，因此忽略这些力对颗粒的作用。对于 模拟，研究采用了单向耦合。模型中考虑了颗粒相的传热过程，其中使用的颗粒物类型为惰性，其传热过程为惰性加热 冷却，离散相中还考虑了热泳力对颗粒物运动的影响。根据拉格朗日轨道追踪法求解了颗粒运动方程得到其运动轨迹，颗粒的运动满足牛顿第二定律，平衡方程为：（）其中 为颗粒质量，；为运动速度，；等式右边依次为颗粒所受曳力、重力、布朗力、升力和热泳力。在 模型中，离散随机漫步模型（）考虑了颗粒的离散旋涡与流体之间的相互作用，从而能够更好地预测颗粒的湍流扩散。本研究采用雷诺应力模型（），其脉动速度考虑了雷诺应力的各向异性：（）（）（）在湍流涡团中，流体脉动速度、和 满足高斯概率密度分布函数，。是正态分布随机数。研究发现，如果直接应用 模型和 模型，会导致对颗粒沉积速率的过度预测，。因为垂直于壁面的速度波动对于计算颗粒沉积速度是很重要的，等的 数据校正近壁区域的壁面法向湍流波动可以准确预测颗粒沉积速度。计算公式为：()，（）其中 为摩擦速度，；表示脉动速度，；常数，它可以通过下式计算得出：（）颗粒物的无量纲弛豫时间 为：（）其中 为运动粘度，；为坎宁安滑移修正系数（）；为粒子与空气的密度比；为颗粒物粒径，。与墙壁的无量纲距离 的定义如下：（）颗粒的沉积速度 由下式得出：（）式中 是沉积的颗粒数量；是平均颗粒浓度，（体积分数）；颗粒的无量纲沉积速度为：（）能 源 环 境 保 护第 卷第 期颗粒沉积率的计算方法为：（）其中 和 分别为沉积颗粒的数目和管道入口释放的颗粒数目。工况描述及求解方法图（）是凹凸酒窝管道 截面结构化示意图，图（，）为 截面结构化示意图。通道宽度 为 ，长度 为 。通道的下表面排列着规律性的凹凸酒窝。通道下表面的凹凸酒窝为，共 个。入口区位于通道前部，以确保凹凸酒窝结构附近的气流均匀，而出口位于通道后部，并在通道后部延伸，以避免空气倒流。根据表，研究共模拟了 个实例。入口风速设为 ，采用出口受压条件，通道壁面和凹凸面均设为无滑动条件。采用了第一类边界条件，规定通道壁面和凹凸酒窝面温度为 ，入口热空气温度为 。图 凹凸酒窝通道结构化网格示意图 研究中所采用的颗粒物为无烟煤颗粒，即“”，密度为 。个颗粒在入口处以圆锥面（）方式入射，通道内地所有粒子都被追踪，直到它们沉积在管道中或离开管道。为了关注凹凸酒窝的存在对通道中颗粒沉积的影响，将凹凸酒窝和下表面的壁面边界条件设置为“”，即粒子一旦接触到上述两个壁面，即视 为 沉 积。通 道 的 出 口 是“”边 界条件。采用有限体积法（）求解湍流风场的控制方程。对流项和扩散项分别采用二阶迎风格式和二阶中心差分格式离散。采用 算法对压力场和速度场进行解耦。采用龙格库塔法求解尘埃粒子的运动方程。表 仿真工况 工况空气速度（）颗粒粒径 光滑管道，粗糙管道，模型验证 网格无关性验证利用网格划分软件 划分了凹凸酒窝通道的六面体结构化网格，并对边界层区域进行了加密处理，网格总数为 ，壁面最近网格中心距离壁面，增长因子为，如图 所示。为探究网格数量对仿真模拟的影响，对模型进行了网格无关性验证。分别挑选了网格数量为：、和 地 种网格数量的模型进行了验证，计算了颗粒粒径分别在、时通道内的颗粒沉积率，如图 所示。结果表明当通道网格数量大于 时，网格数量对颗粒沉积的结果影响 较 小，因 此 本 研 究 选 用 了 网 格 数 量 为 的网格。流场验证选择光滑管道中 处的湍流流速剖韩尊师等 凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究面，是因为此时流场已经完全演化。如图 所示，本文 模型的计算结果与 等的直接数值模拟（）结果非常吻合。颗粒沉积验证图 描述了颗粒在光滑通道中的无量纲沉积速度分布与无量纲弛豫时间之间的联系，并与过去的实验数据和数值模拟进行了比较，。随着无量纲弛豫时间的增加，无量纲沉积速率先增加，然后趋于稳定。图 表明，颗粒的无量纲临界松弛周期约为。当 时，沉积特征主要由流体的湍流涡动决定，但当 时，颗粒的惯性占主导地位。无量纲沉积速率与之前的实验结果一致。因此，这种数值技术可用于预测颗粒沉积的特征。图 网格无关性验证折线图 图 光滑通道速度数值验证图 图 光滑通道颗粒沉积数值验证图 结果分析 通道流场分析热交换通道中，凹凸酒窝的存在会极大地影响通道中的速度场和温度场，以及通道壁面处气流的流动状态。如图（）所示，从 处开始，可以明显观察到沿着凹凸酒窝的周期性速度边界层。在凹酒窝空腔内部气流速度较小，凸酒窝后方存在低流速区。在凸酒窝后面的凹酒窝空腔内的速度要大于左侧第一个凹酒窝空腔内速度，这是因为凸酒窝的存在影响了通道内的流场结构，为其后面的凹酒窝起到了引流的作用。凹凸酒窝的存在还影响了通道内的温度场。如图（）所示，凹酒窝空腔内部是低温区，通道中部位置是高温区。温度由通道中部靠近壁面递减，此时则会存在由高温区域指向低温区域的热泳力，从而促进颗粒趋壁沉积。另外通过观察通道 截面的速度矢量图，在通道壁面处和凹酒窝空腔内部存在二次流动，可以促进颗粒物在此处的沉积。图 展示了通道中 截面的 值（）和流场流线图，高 值在凹酒窝的右上位置，凸酒窝的迎风面上方和背风面的后方。凹酒窝空腔内和凸酒窝后方存在明显的大涡团，高 值和大涡团的存在会极大地影响通道内的颗粒沉积特性。颗粒沉积分析通道中凹凸酒窝的存在会影响通道内的流场进而影响颗粒物在通道内的运动轨迹，图 为流场均匀时颗粒轨迹分布图，在凹凸酒窝上方颗粒物较密集而在表面光滑区域则颗粒物较分散。能 源 环 境 保 护第 卷第 期图 （）凹凸酒窝速度云图，（，）截面温度场及二次流动示意图 （），（，）图 截面 值及流场流线图 图 凹凸酒窝通道颗粒轨迹图 在光滑通道内进行了颗粒物沉积实验，选用的是美国 灰尘，沉积情况如图 所示，同时跟光滑通道仿真对比，沉积形貌几乎一致，证明了本仿真模拟颗粒沉积模型的正确性。图 展示了不同粒径的颗粒物在凹凸酒窝通道中的三维沉积示意图。图（）展示了粒径为 的颗粒部分沉积在凹酒窝空腔内部和凸酒窝迎风面上，但是光滑表面处也有部分颗粒沉积。此时颗粒物质量较小，受流场影响较大，流场均匀后颗粒物沉积分布也较为均匀。图（）中 颗粒集中沉积在凹酒窝空腔内、凸酒窝迎风面以及凹凸酒窝之间的光滑壁面。这是因为颗粒粒径增大，重力和颗粒自身的惯性对颗粒物的沉积影响变大，因此凸酒窝的迎风面沉积颗粒数目增多。图（）的颗粒物主要沉积在凸酒窝迎风面以及第一个和最后一个凸酒窝两侧光滑壁面处。此时颗粒粒径进一步增大，受流场、重力以及颗粒物自己惯性的影响，颗粒物集中沉积于凸酒窝迎风面，凹酒窝右上部以及光滑壁面。当颗粒粒径较小时，凹凸酒窝通道内颗粒物沉积率大于表面光滑通道内。随颗粒粒径增大，总体颗粒物在通道中的沉积数目越多，然而两种通道中的沉积率趋于一致，当粒径增大到 时凹凸酒窝通道中沉积率略低于表面光滑通道，如图 所示。韩尊师等 凹凸酒窝表面换热管道内颗粒沉积与传热特性的数值模拟研究图 光滑通道颗粒沉积实验仿真对比图 图 （，）、及 颗粒沉积示意图 （，），图 通道颗粒物沉积率对比 结 论研究了颗粒物在具有凹凸酒窝结构的三维传热通道中的沉积情况。采用了结构化网格进行了网格划分，流体相采用雷诺应力模型（），颗粒相采用离散颗粒模型。研究并讨论了气流场、粒子沉积特性和沉积效率。研究了壁面结构和颗粒物对流场和沉积特性的影响。结论如下：（）凹凸酒窝的存在会影响酒窝通道中的流场，并强化通道中的传热。它还会影响颗粒在通道中的运动轨迹。更多的颗粒聚集在凹凸酒窝上方，酒窝的存在促进了小颗粒（）的沉积。（）颗粒粒径大小的差异会明显影响颗粒在通道中的沉积。颗粒沉积速率随颗粒直径的增大而增大，凹凸酒窝结构的存在对小颗粒（）的沉积有很好的促进作用。参考文献（）：，（）：，：，：，：能 源 环 境 保 护第 卷第 期 ，：，：，（）：，：，（）：，（）：，：，（）：，（）：，：（），（）：，：，：，（）：，（）：（上接 页），：，：，：，：，：，：，（）：，：，（）：，：，：，：，：，（）：，：，：，：，：，：