周期性扩缩微通道：单相与两相流动传热特性的深入剖析

周期性扩缩微通道：单相与两相流动传热特性的深入剖析一、绪论1.1研究背景与意义随着现代科技的迅猛发展，电子设备、能源系统、生物医学等众多领域都在朝着微型化、集成化和高性能化的方向迈进，这使得设备单位体积内的功率密度急剧增加，产生的热量也大幅增多。高效散热技术已然成为制约这些领域进一步发展的关键因素。微通道技术作为一种极具潜力的散热解决方案，因其具有比表面积大、传热效率高、结构紧凑等诸多优点，在过去几十年间受到了学术界和工业界的广泛关注，并在多个领域得到了成功应用。在电子领域，微通道热沉被广泛应用于芯片散热。随着芯片集成度的不断提高，如英特尔酷睿i7处理器中晶体管数量已达到7×10⁸个，散热功率高达106W/m²，传统的散热方式已难以满足需求。微通道热沉能够有效降低芯片温度，提高电子设备的性能和可靠性。在能源领域，微通道换热器可应用于太阳能集热器、燃料电池等设备，提高能源转换效率，促进能源的高效利用。在生物医学领域，微通道可用于生物芯片的构建，实现对生物分子的快速分离和检测，其精确的温度控制能力有助于生物反应的进行。然而，常规的等截面微通道在传热性能和流动阻力方面存在一定的局限性。周期性扩缩微通道作为一种新型的微通道结构，通过改变通道截面的大小，使流体在通道内产生复杂的流动现象，如漩涡、二次流等，从而增强流体与壁面之间的传热传质过程，有望在降低流动阻力的同时提高传热效率。对周期性扩缩微通道单相及两相流动传热特性的深入研究，不仅能够丰富微尺度下流体力学和传热学的理论体系，还可为新型散热设备的设计和优化提供坚实的理论依据，对推动散热技术的发展具有重要的现实意义。1.2国内外研究现状1.2.1微通道内单相液体流动与传热研究微通道内单相液体的流动与传热特性一直是研究的热点。早期，Tuckerman和Pease在1981年率先提出了微通道热沉的概念，并通过实验与理论分析，对矩形微通道内的强制对流换热进行了研究，为后续的研究奠定了基础。此后，众多学者围绕微通道内单相液体的流动与传热开展了大量的研究工作。在流动特性方面，研究发现微通道内的流体流动存在明显的尺寸效应。由于微通道的特征尺度通常在微米至毫米量级，表面积与体积比大幅增加，表面效应变得尤为重要。壁面粗糙度、表面极性等因素会对流体的流动产生显著影响，使得微通道内的流动现象比常规尺度下更为复杂。一些实验研究表明，微通道内的流动可能会出现非牛顿流体行为，即使是牛顿流体，其流动特性也可能偏离传统的Navier-Stokes方程预测结果。在传热特性方面，学者们通过实验和数值模拟等方法，研究了不同因素对微通道内传热的影响。Kandlikar和Steinke运用CFD软件对微通道内的单相流传热进行了数值模拟，详细分析了不同工况下的传热特性，发现传热系数会随着雷诺数、普朗特数等参数的变化而发生改变。国内学者李震等利用数值模拟方法研究了不同截面形状微通道内的流动与传热特性，结果表明矩形微通道在一定条件下具有较好的综合性能。此外，微通道内的传热还受到流体物性、通道表面性质等因素的影响。然而，目前对于微通道内单相液体流动与传热的研究仍存在一些不足之处。部分研究仅考虑了单一因素对流动和传热的影响，而实际应用中，多种因素往往相互作用，使得流动和传热现象更加复杂。对于微通道内复杂流动现象的机理，如漩涡、二次流等的形成和发展机制，尚未完全明确，需要进一步深入研究。1.2.2微通道内两相流型与压降特性研究微通道内的气液两相流在化学、化工、能源等领域具有广泛的应用背景，因此受到了众多学者的关注。微通道内的气液两相流动特性与常规尺度下有很大区别，主要是因为在微通道内流体间的表面张力、润湿性等占主导地位，使得气泡的生成、运动、聚并和破裂等过程变得更加复杂。在流型研究方面，国内外学者通过实验观测和数值模拟等手段，对微通道内的气液两相流型进行了大量研究。Thome等对微通道内的气液两相流流动与传热进行了深入研究，通过高速摄像等可视化技术，观察到微通道内存在泡状流、弹状流、环状流等多种流型，并揭示了流型转变的规律。国内学者柴磊和夏国栋采用高速摄像仪和宽视场体式显微镜相结合的可视化观测技术，实验研究了三角凹穴型微通道热沉内氮气/水气液两相流型，发现微通道热沉内为气液两相和单相液体周期性地交替流过，且流型完全不同于等截面微通道。在压降特性研究方面，压降是衡量微通道内气液两相流动性能的重要参数之一。曹传超、冯龙龙和贾洪伟通过实验研究了制冷剂R134a在两种微通道尺寸下（0.5mm和1mm）的两相压降特性，发现摩擦压力梯度随干度的升高呈现出先增大后减小的趋势，且不同尺寸微通道的压降特性存在差异。然而，由于微通道内气液两相流的复杂性，目前已有的压降模型在预测微通道内的压降时，往往存在一定的误差，尤其是对于复杂流型和特殊工况下的压降预测，准确性有待提高。当前研究虽然取得了一定的成果，但仍存在一些问题。微通道内气液两相流的流型转变机制尚未完全清晰，不同流型之间的过渡条件和影响因素还需要进一步深入研究。对于微通道内两相压降特性的研究，缺乏统一的理论模型，已有的模型大多是基于特定的实验条件和流型建立的，适用范围有限。1.3研究目标与内容本研究旨在深入了解周期性扩缩微通道的单相及两相流动传热特性，揭示其内在的物理机制，为新型散热设备的设计和优化提供理论依据。具体研究内容如下：周期性扩缩微通道单相流动传热特性研究：通过数值模拟和实验研究相结合的方法，研究不同雷诺数、普朗特数以及通道几何参数（如扩缩比、周期长度等）对周期性扩缩微通道内单相液体流动和传热特性的影响。分析流体在通道内的速度分布、温度分布以及传热系数的变化规律，探讨漩涡、二次流等复杂流动现象对传热的强化机制。周期性扩缩微通道两相流型与压降特性研究：采用高速摄像等可视化技术，实验观测周期性扩缩微通道内气液两相流的流型，并与数值模拟结果进行对比分析。研究不同工况（如质量通量、热通量、饱和温度等）下的流型转变规律，确定流型转变的临界条件。同时，测量周期性扩缩微通道内气液两相流的压降，分析压降与流型、工况参数之间的关系，建立适用于周期性扩缩微通道的两相压降模型。周期性扩缩微通道传热性能优化研究：基于上述研究结果，对周期性扩缩微通道的结构进行优化设计，以提高其传热性能。通过改变通道的几何形状、尺寸以及表面性质等参数，寻找最优的设计方案。采用多目标优化算法，综合考虑传热效率和流动阻力等因素，实现周期性扩缩微通道性能的最大化。二、周期性扩缩微通道流动传热实验装置及方法2.1实验装置2.1.1单相液体流动与传热特性测试实验平台为了深入研究周期性扩缩微通道内单相液体的流动与传热特性，搭建了一套高精度的实验平台，该平台主要由液体循环系统、实验段、数据采集系统和温度控制系统等部分组成。液体循环系统采用高精度的齿轮泵，能够精确控制液体的流量，流量调节范围为0.1-10L/min，精度可达±0.01L/min，以满足不同实验工况下的流量需求。实验段由周期性扩缩微通道芯片和加热模块组成，微通道芯片采用MEMS加工工艺制作，确保通道尺寸的精确性和表面的光滑度。加热模块采用薄膜电阻加热器，能够实现对微通道壁面的均匀加热，加热功率可在0-100W范围内调节，精度为±0.1W，通过控制加热功率来模拟不同的热流密度工况。数据采集系统配备了高精度的压力传感器和温度传感器。压力传感器用于测量微通道进出口的压力差，精度为±0.01kPa，能够准确获取流体在通道内的流动阻力。温度传感器采用T型热电偶，测量精度为±0.1℃，分别布置在微通道的进出口以及壁面上，用于实时监测流体和壁面的温度变化。所有传感器采集的数据通过数据采集卡传输至计算机，利用专业的数据采集软件进行实时记录和分析。温度控制系统采用恒温槽，能够将液体的入口温度精确控制在设定值，温度波动范围不超过±0.1℃，为实验提供稳定的入口温度条件。在实验过程中，通过调节恒温槽的温度和液体流量，改变雷诺数和热流密度等实验参数，研究不同工况下周期性扩缩微通道内单相液体的流动与传热特性。2.1.2并行硅基微通道内气液两相流型与压降测试实验平台该实验平台主要用于研究并行硅基微通道内气液两相流的流型和压降特性，由气液供应系统、实验段、可视化观测系统和数据采集系统等组成。气液供应系统分别采用气体质量流量控制器和液体注射泵来精确控制气体和液体的流量。气体质量流量控制器的量程为0-1000sccm，精度为±1%FS，可提供稳定的气体流量。液体注射泵的流量范围为0-10mL/min，精度可达±0.001mL/min，确保液体流量的精确控制。实验段由并行硅基微通道芯片组成，微通道采用MEMS加工工艺制作，通道尺寸精确，表面质量高。芯片上下表面分别设置了进液口和进气口，使气液两相能够在微通道内充分混合流动。可视化观测系统采用高速摄像仪，帧率可达10000fps，配合高分辨率显微镜镜头，能够清晰捕捉微通道内气液两相流的流型变化。通过将高速摄像仪与计算机相连，实时记录流型图像，并利用图像处理软件对图像进行分析，确定不同流型的特征参数和转变条件。数据采集系统与单相液体实验平台类似，配备了高精度的压力传感器，用于测量微通道进出口以及不同位置处的压力，精度为±0.01kPa，以获取气液两相流的压降数据。通过改变气液流量比、质量通量等实验参数，结合可视化观测和压降测量，深入研究并行硅基微通道内气液两相流型与压降的特性。2.1.3PDMS单通道内气液两相流型及其演变特性测试实验平台此实验平台专门用于研究PDMS单通道内气液两相流型及其演变特性，主要由气液供应装置、PDMS微通道实验件、高速摄像系统和数据处理单元等构成。气液供应装置同样采用高精度的气体质量流量控制器和液体注射泵，能够精确调节气体和液体的流量，为实验提供稳定的气液两相流条件。PDMS微通道实验件通过PDMS加工工艺制作而成，具有良好的柔韧性和透光性，便于可视化观测。微通道的尺寸根据实验需求进行设计和制作，通道表面经过特殊处理，以控制表面的润湿性，研究润湿性对气液两相流型的影响。高速摄像系统由高速摄像机和照明光源组成，高速摄像机的帧率最高可达20000fps，分辨率为1280×1024像素，能够清晰拍摄到微通道内气液两相流的细微变化。照明光源采用高亮度的LED冷光源，通过光纤传导至实验件下方，为高速摄像提供充足的照明，确保拍摄的图像质量清晰、对比度高。数据处理单元主要负责对高速摄像获取的图像数据进行处理和分析。利用专业的图像处理软件，对图像进行降噪、增强、分割等操作，提取气液两相流的流型信息，如气泡尺寸、形状、分布等参数。通过对不同时刻的图像进行对比分析，研究气液两相流型的演变规律，以及流型演变与气液流量比、表面润湿性等因素之间的关系。2.2实验件加工工艺2.2.1MEMS加工工艺MEMS加工工艺是制作硅基实验件的关键技术，其具有高精度、可重复性好等优点，能够精确控制微通道的尺寸和形状。本研究中，利用MEMS加工工艺制作周期性扩缩微通道的主要流程如下：光刻：首先在硅片表面涂覆一层光刻胶，通过光刻掩模版将设计好的微通道图案曝光在光刻胶上。光刻掩模版采用高分辨率的石英玻璃制作，图案精度可达亚微米级。曝光过程中，利用紫外线照射光刻胶，使曝光区域的光刻胶发生化学反应，溶解性发生改变。然后通过显影工艺，去除未曝光区域的光刻胶，从而在硅片表面形成与掩模版图案一致的光刻胶图案。刻蚀：采用反应离子刻蚀（RIE）技术对硅片进行刻蚀，去除未被光刻胶保护的硅材料，形成微通道结构。RIE刻蚀过程中，将硅片放置在真空腔室内，通入反应气体（如CF4、O2等），在射频电源的作用下，反应气体被电离形成等离子体。等离子体中的离子在电场作用下加速撞击硅片表面，与硅材料发生化学反应，从而实现对硅的刻蚀。通过精确控制刻蚀时间、刻蚀功率和反应气体流量等参数，能够精确控制微通道的深度和侧壁垂直度。掺杂：为了改变硅材料的电学性能，提高微通道的传热性能，对微通道壁面进行掺杂处理。采用离子注入或扩散工艺，将杂质原子（如硼、磷等）引入硅材料中。离子注入过程中，利用高能离子束将杂质原子注入硅片表面，通过控制离子能量和注入剂量，精确控制杂质原子的分布深度和浓度。扩散工艺则是将硅片放置在高温炉中，在含有杂质原子的气氛中进行扩散，使杂质原子向硅片内部扩散，形成一定的浓度分布。清洗与封装：刻蚀和掺杂完成后，对硅片进行清洗，去除表面的光刻胶、杂质和反应物等。清洗过程采用去离子水、丙酮、乙醇等有机溶剂进行多次清洗，确保硅片表面干净无污染。最后，对硅基实验件进行封装，将硅片与玻璃盖板通过阳极键合或热压键合等方式进行密封，形成完整的微通道器件。封装过程中，要确保键合质量良好，避免出现漏液、漏气等问题。2.2.2PDMS加工工艺PDMS加工工艺具有制作成本低、工艺简单、灵活性高等优点，适合制作用于气液两相流实验的微通道实验件。PDMS加工工艺的具体操作过程如下：模具制作：首先根据微通道的设计图案，利用光刻、刻蚀等MEMS加工工艺在硅片上制作出微通道的模具。模具的表面粗糙度和尺寸精度对PDMS微通道的质量有重要影响，因此要确保模具制作的精度和质量。PDMS混合与浇铸：将PDMS预聚体和固化剂按照一定比例（通常为10:1）在烧杯中混合均匀，然后将混合液倒入真空干燥箱中进行脱气处理，去除混合液中的气泡。脱气完成后，将PDMS混合液缓慢浇铸在硅片模具上，使混合液充分填充模具的微通道图案。为了确保浇铸质量，可在浇铸过程中轻轻晃动硅片，使混合液均匀分布。固化与脱模：将浇铸好PDMS混合液的硅片放入烘箱中进行固化，固化温度一般为60-80℃，固化时间为2-4小时。固化完成后，将硅片从烘箱中取出，冷却至室温，然后小心地将PDMS微通道从模具上脱模。脱模过程中要注意避免对微通道造成损伤，可使用镊子等工具辅助脱模。打孔与清洗：根据实验需求，在PDMS微通道的两端打孔，作为气液进出口。打孔过程可使用打孔器或激光打孔设备，确保孔的位置和尺寸准确。打孔完成后，用去离子水和乙醇对PDMS微通道进行清洗，去除表面的杂质和残留的固化剂。清洗完成后，将PDMS微通道晾干或放入烘箱中烘干。键合与老化：将清洗后的PDMS微通道与玻璃基底或其他PDMS片材进行键合，形成完整的实验件。键合过程可采用等离子体处理键合或化学粘接剂键合等方法。等离子体处理键合是将PDMS微通道和基底分别放入等离子体清洗机中进行处理，使表面产生活性基团，然后将两者紧密贴合，在室温下放置一段时间即可实现键合。化学粘接剂键合则是在PDMS微通道和基底表面涂抹适量的粘接剂，然后将两者贴合，在一定温度和压力下固化，实现键合。键合完成后，将实验件放入烘箱中进行老化处理，进一步提高键合强度和实验件的稳定性。老化温度一般为60-80℃，老化时间为1-2小时。2.3实验件结构及尺寸2.3.1MEMS加工的硅基实验件结构MEMS加工的硅基实验件结构设计是实现周期性扩缩微通道流动传热特性研究的关键。硅基实验件主要由硅片主体、周期性扩缩微通道、进出液口和加热电阻等部分组成，如图1所示。图1硅基实验件结构示意图周期性扩缩微通道采用折线形设计，这种设计能够使流体在通道内产生更复杂的流动形态，增强传热传质效果。微通道的扩大部分和缩大部分交替排列，形成周期性变化的截面。扩大部分的宽度为W1，长度为L1；缩大部分的宽度为W2，长度为L2。通道的深度为H，整个微通道的周期长度为T=L1+L2。通过精确控制这些尺寸参数，可以研究不同几何结构对流动传热特性的影响。进出液口位于硅片的两端，采用圆形设计，直径为D，便于与外部管路连接。加热电阻采用薄膜电阻制作，均匀分布在微通道的底部，通过施加电流来实现对微通道壁面的加热，模拟不同的热流密度工况。加热电阻的材料为铂，具有良好的电阻温度系数和稳定性，能够精确控制加热功率。2.3.2MEMS加工的硅基实验件的封装硅基实验件的封装对于保证实验件在实验中的性能稳定至关重要。本研究采用玻璃盖板与硅片通过阳极键合的方式进行封装。在封装前，首先对硅片和玻璃盖板进行清洗和预处理，去除表面的杂质和氧化物，提高键合质量。清洗过程采用去离子水、丙酮、乙醇等有机溶剂依次进行超声清洗，然后用氮气吹干。预处理过程中，对硅片表面进行亲水处理，使其表面带有羟基，便于与玻璃盖板进行键合。阳极键合过程中，将硅片和玻璃盖板紧密贴合，放入键合设备中。在高温（约300-400℃）和高电压（约500-1000V）的作用下，玻璃中的钠离子向阴极（硅片）迁移，在界面处形成一层很薄的氧化硅层，实现硅片与玻璃盖板的紧密键合。键合完成后，对封装好的实验件进行气密性检测，确保无漏液、漏气现象。通过封装，不仅能够保护硅基实验件的微通道结构不受外界环境的影响，还能够为实验提供稳定的边界条件，保证实验数据的准确性和可靠性。2.3.3PDMS实验件结构PDMS实验件主要用于气液两相流型及其演变特性的研究，其结构设计具有独特的优势。PDMS实验件由PDMS微通道、进出液口和观察窗口等部分组成，如图2所示。图2PDMS实验件结构示意图PDMS微通道采用直通道设计，通道的宽度为W，深度为H。直通道设计便于观察气液两相在通道内的流动形态和流型演变过程。进出液口位于PDMS实验件的两端，采用圆形通孔设计，直径为D，与外部气液供应系统连接。观察窗口位于PDMS微通道的上方，采用透明的PDMS材料制作，厚度为t。观察窗口的尺寸较大，能够提供足够的观察视野，便于利用高速摄像系统对微通道内的气液两相流进行可视化观测。为了提高观察效果，观察窗口的表面经过抛光处理，减少光线的散射和折射。PDMS实验件具有良好的柔韧性和透光性，能够实时观察气液两相在微通道内的流动情况，并且对气液两相的流动干扰较小，有利于准确研究气液两相流型及其演变特性。2.4实验步骤2.4.1单相流体流动传热实验步骤实验准备：检查实验装置各部分连接是否正确、牢固，确保液体循环系统、数据采集系统和温度控制系统等设备正常运行。将实验所需的去离子水或其他实验液体加入到液体循环系统的储液罐中，并开启恒温槽，将液体的入口温度调节至设定值。实验段安装：将制作好的周期性扩缩微通道硅基实验件安装到实验装置的实验段中，确保实验件与加热模块、压力传感器和温度传感器等设备连接紧密，无漏液现象。初始数据测量：在不加热的情况下，启动齿轮泵，调节液体流量至设定值，稳定运行一段时间后，记录微通道进出口的压力、温度以及液体的流量等初始数据。加热实验：根据实验需求，调节加热模块的功率，使微通道壁面达到设定的热流密度。待系统稳定后，每隔一定时间记录一次微通道进出口的压力、温度、液体流量以及壁面温度等数据，确保数据的准确性和可靠性。改变工况：依次改变液体流量和热流密度等实验工况，重复步骤4，测量不同工况下的实验数据。每个工况下至少测量3组数据，取平均值作为该工况下的实验结果，以减小实验误差。实验结束：完成所有实验工况的测量后，关闭加热模块和齿轮泵，停止实验。将实验液体从储液罐中排出，清洗实验装置和实验件，整理实验设备。2.4.2气液两相流型及压降测试实验步骤实验准备：检查气液供应系统、实验段、可视化观测系统和数据采集系统等设备是否正常工作。确保气体质量流量控制器和液体注射泵的精度满足实验要求，高速摄像仪和压力传感器等设备校准准确。准备好实验所需的气体（如氮气）和液体（如水或其他实验液体），并将其分别接入气液供应系统。实验段安装：将并行硅基微通道芯片或PDMS微通道实验件安装到实验段中，确保气液进出口与气液供应系统连接紧密，无漏气、漏液现象。调整高速摄像仪的位置和焦距，使其能够清晰拍摄到微通道内气液两相流的流型。初始条件设定：设定气体和液体的初始流量，一般先将气体流量设置为较小值，液体流量设置为较大值。开启气体质量流量控制器和液体注射泵，使气液两相开始流入微通道。流型观测与数据采集：利用高速摄像仪实时观察微通道内气液两相流的流型，并记录流型图像。同时，数据采集系统实时采集微通道进出口以及不同位置处的压力数据。每隔一定时间，调整气体和液体的流量比，观察流型的变化，并记录相应的流型图像和压力数据。流型转变研究：逐渐改变气液流量比，从低气液流量比逐渐增加到高气液流量比，观察流型的转变过程。记录不同流型转变时的气液流量、质量通量等参数，确定流型转变的临界条件。实验结束：完成所有实验工况的测量后，先关闭气体质量流量控制器和液体注射泵，停止气液供应。然后关闭高速摄像仪和数据采集系统，拆除实验件，清洗实验装置和设备，整理实验现场。2.5实验数据的分析与处理2.5.1单相液体流动传热实验数据分析与处理流量计算：根据齿轮泵的流量调节参数和校准曲线，计算实验过程中液体的实际流量Q（单位：m³/s）。流量计算公式为：Q=k\times\DeltaP其中，k为齿轮泵的流量系数，通过校准实验确定；\DeltaP为齿轮泵进出口的压力差。雷诺数计算：三、周期性扩缩微通道内单相液体流动传热特性3.1流动与传热实验数据分析3.1.1流动阻力实验数据分析在周期性扩缩微通道内，流动阻力的特性十分关键，它与流量、通道结构等因素紧密相关。对实验数据进行深入分析后，能够揭示这些因素之间的内在联系。流量的变化会直接影响流体在通道内的流速和流态，进而改变流动阻力。当流量较小时，流体处于层流状态，此时粘性力起主导作用，流动阻力与流量近似呈线性关系。随着流量的增加，流速增大，流体逐渐过渡到湍流状态，惯性力的影响逐渐增强，流动阻力与流量的关系变得更加复杂，不再是简单的线性关系，而是呈现出非线性的增长趋势。通道结构是影响流动阻力的重要因素之一。扩缩比，即通道扩大部分与缩大部分的尺寸比例，对流动阻力有着显著影响。较大的扩缩比会使流体在通道内的流速变化更加剧烈，在扩大部分，流体流速降低，压力升高；在缩大部分，流体流速增大，压力降低。这种流速和压力的大幅变化会导致流动阻力增加。周期长度，也就是扩缩结构重复出现的间隔距离，也会影响流动阻力。较短的周期长度会使流体频繁地经历扩缩过程，加剧流体的扰动，从而增加流动阻力；而较长的周期长度则相对减少了流体的扰动，流动阻力相对较小。通过对不同流量和通道结构下的实验数据进行回归分析，得到了流动阻力与各因素之间的经验关联式。该关联式能够较为准确地描述周期性扩缩微通道内的流动阻力特性，为工程应用提供了重要的参考依据。例如，在设计微通道热沉时，可以根据实际需求，利用该关联式优化通道结构参数，以降低流动阻力，提高散热效率。3.1.2传热特性实验数据分析传热特性是周期性扩缩微通道的另一个重要研究内容，它与温度、流速等参数密切相关。传热系数是衡量传热性能的关键指标，它反映了单位时间内单位面积上传递的热量。通过对实验数据的详细分析，研究传热系数与各参数之间的关系，对于理解微通道内的传热机制具有重要意义。温度对传热系数有着显著影响。随着流体温度的升高，流体的物性参数，如导热系数、粘度等会发生变化。导热系数的增大有利于热量的传递，从而提高传热系数；而粘度的变化则会影响流体的流动状态，进而间接影响传热系数。在一定温度范围内，传热系数可能会随着温度的升高而增大，但当温度过高时，可能会出现其他因素的影响，导致传热系数的变化趋势发生改变。流速是影响传热系数的另一个重要因素。当流速较低时，流体的对流换热较弱，传热主要依靠导热进行，此时传热系数较小。随着流速的增加，对流换热逐渐增强，流体与壁面之间的热量传递加快，传热系数显著增大。当流速进一步增加，达到一定程度后，传热系数的增长趋势会逐渐变缓，这是因为此时流体的湍流程度已经较高，再增加流速对传热的强化效果有限。通过对实验数据的拟合，得到了传热系数与温度、流速等参数之间的经验公式。该公式能够定量地描述传热系数与各参数之间的关系，为工程设计和分析提供了有力的工具。在实际应用中，可以根据该公式预测不同工况下的传热性能，从而优化微通道的设计，提高传热效率。例如，在设计微通道换热器时，可以根据实际的工作温度和流速要求，利用该公式选择合适的通道尺寸和结构，以实现高效的热量传递。3.2微通道内单相流动与传热的数值模型3.2.1基本物理模型为了深入研究周期性扩缩微通道内单相液体的流动与传热特性，建立了合理的物理模型。该模型基于以下假设：流体为不可压缩的牛顿流体，这意味着流体的密度不随压力和温度的变化而改变，且其粘性符合牛顿粘性定律，即切应力与速度梯度成正比；流动为稳态层流，在稳态条件下，流体的速度、压力等物理量不随时间变化，层流状态下流体的流动较为规则，没有明显的湍流脉动；忽略重力等体积力的作用，在微通道内，流体的流动主要受粘性力和惯性力的影响，重力等体积力相对较小，可以忽略不计；流体的物性参数，如密度、粘度、导热系数等，均为常数，不随空间位置和温度的变化而改变。在建立物理模型时，充分考虑了微通道的周期性扩缩结构。微通道的扩大部分和缩大部分交替排列，形成周期性变化的截面。这种结构会使流体在通道内产生复杂的流动现象，如漩涡、二次流等，从而增强传热传质效果。通过对微通道结构的精确描述，能够更准确地模拟流体在通道内的流动与传热过程。该物理模型适用于研究微通道内单相液体在一定工况下的流动与传热特性，为后续的数值模拟提供了基础。在实际应用中，需要根据具体的研究对象和工况条件，对模型进行适当的调整和验证，以确保模型的准确性和可靠性。3.2.2边界条件设置在数值模拟中，边界条件的设置对计算结果的准确性至关重要。对于周期性扩缩微通道内单相液体的流动与传热问题，主要设置以下边界条件：入口边界条件采用速度入口，给定流体的入口速度和温度。入口速度根据实验条件或实际工况确定，温度则根据实验测量值或热管理需求设定。速度入口条件能够准确地模拟流体进入微通道时的初始状态，为后续的流动与传热计算提供准确的起始条件。出口边界条件采用压力出口，设定出口压力为环境压力。在实际工程中，微通道的出口通常与大气相通，压力接近环境压力。采用压力出口边界条件可以合理地模拟流体从微通道流出时的状态，确保计算结果的真实性。壁面边界条件采用无滑移边界条件和恒定热流边界条件。无滑移边界条件意味着流体在壁面上的速度为零，这符合实际情况，因为流体与壁面之间存在粘性力，使得流体在壁面处无法滑动。恒定热流边界条件则给定壁面的热流密度，模拟微通道壁面的加热或冷却过程。通过设置合适的热流密度，可以研究不同热负荷下微通道内的传热特性。这些边界条件的设置是基于对实际物理过程的理解和简化，能够准确地反映微通道内单相液体流动与传热的实际情况。在数值模拟过程中，需要严格按照这些边界条件进行设置，以确保计算结果的准确性和可靠性。同时，对于一些复杂的实际工况，可能需要根据具体情况对边界条件进行适当的调整和改进。3.2.3控制方程离散格式及迭代精度的控制在数值模拟中，选择合适的控制方程离散格式对于保证计算结果的准确性和稳定性至关重要。对于周期性扩缩微通道内单相液体的流动与传热问题，采用有限体积法对控制方程进行离散。有限体积法的基本思想是将计算区域划分为一系列的控制体积，然后在每个控制体积上对控制方程进行积分，将偏微分方程转化为代数方程。这种方法具有守恒性好、计算精度高、适应性强等优点，能够有效地处理复杂的几何形状和边界条件。在有限体积法中，对流项采用二阶迎风格式进行离散。二阶迎风格式考虑了流速方向上的梯度变化，能够更准确地模拟对流过程，减少数值扩散误差。扩散项采用中心差分格式进行离散，中心差分格式具有较高的精度，能够准确地计算扩散通量。通过合理选择对流项和扩散项的离散格式，能够提高计算结果的精度和稳定性。为了保证计算结果的收敛性，需要对迭代精度进行严格控制。在迭代过程中，设置了残差收敛标准，当各物理量的残差小于设定的收敛阈值时，认为计算结果已经收敛。残差是指迭代过程中当前迭代步与上一迭代步之间物理量的差值，通过监测残差的变化可以判断计算是否收敛。通常，收敛阈值设置为10⁻⁶~10⁻⁸之间，具体数值根据计算问题的复杂程度和精度要求进行调整。在实际计算中，还需要密切关注迭代过程中各物理量的变化情况，如速度、压力、温度等，确保计算结果的合理性。如果发现计算结果出现异常波动或不收敛的情况，需要检查离散格式、边界条件等设置是否正确，并进行相应的调整。3.2.4网格划分对微通道模型进行网格划分是数值模拟的重要步骤，它直接影响到计算效率和精度。在网格划分过程中，采用结构化网格对微通道模型进行离散。结构化网格具有网格质量高、计算效率高、数据存储量小等优点，能够有效地提高数值模拟的精度和速度。对于周期性扩缩微通道这种复杂的几何结构，在网格划分时，对扩缩部分和拐角处进行了局部加密。扩缩部分和拐角处是流体流动和传热变化较为剧烈的区域，局部加密网格可以更准确地捕捉这些区域的流动和传热细节，提高计算精度。在加密过程中，根据几何形状和物理量的变化梯度，合理确定网格的疏密程度，避免出现网格质量过低或计算量过大的问题。为了验证网格的独立性，进行了网格无关性测试。通过逐渐加密网格，比较不同网格数量下的计算结果。当网格数量增加到一定程度后，计算结果的变化小于设定的误差范围，此时认为网格已经足够细密，计算结果与网格数量无关。经过网格无关性测试，确定了合适的网格数量，既保证了计算精度，又不会过多增加计算时间和内存消耗。在实际应用中，网格划分还需要考虑计算资源的限制和计算效率的要求。对于大规模的数值模拟问题，需要在保证计算精度的前提下，合理优化网格划分，以提高计算效率，降低计算成本。3.2.5沿程计算公式在数值模拟中，需要给出计算沿程阻力和传热的公式，为模拟提供理论依据。沿程阻力是指流体在微通道内流动时，由于粘性力的作用而产生的阻力。根据流体力学理论，沿程阻力可以通过达西-韦斯巴赫公式进行计算：\DeltaP_f=f\frac{L}{D_h}\frac{\rhov^2}{2}其中，\DeltaP_f为沿程阻力，单位为Pa；f为摩擦因数，它与雷诺数和通道壁面粗糙度等因素有关；L为微通道的长度，单位为m；D_h为水力直径，单位为m，对于非圆形截面的微通道，水力直径的计算公式为D_h=\frac{4A}{P}，其中A为通道截面积，P为湿周；\rho为流体密度，单位为kg/m^3；v为流体平均流速，单位为m/s。传热的计算公式基于牛顿冷却定律：q=h\DeltaT其中，q为热流密度，单位为W/m^2；h为传热系数，单位为W/(m^2\cdotK)，它与流体的物性、流速、通道壁面条件等因素有关；\DeltaT为流体与壁面之间的温差，单位为K。传热系数可以通过经验公式或数值模拟方法得到，在周期性扩缩微通道中，传热系数的计算较为复杂，需要考虑通道结构、流动状态等多种因素的影响。这些沿程计算公式是数值模拟的基础，通过准确计算沿程阻力和传热，可以深入研究周期性扩缩微通道内单相液体的流动与传热特性。在实际应用中，需要根据具体的工况条件和实验数据，合理选择和修正这些公式，以提高计算结果的准确性。3.3压降特性分析3.3.1总体压降特性周期性扩缩微通道内的总体压降特性是研究其流动性能的重要指标，它与流量、通道几何形状密切相关。随着流量的增加，流体在微通道内的流速增大，惯性力增强，导致总体压降呈现出上升的趋势。这是因为流速的增加使得流体与壁面之间的摩擦加剧，同时在扩缩部分和拐角处，流体的流动方向和速度发生剧烈变化，产生了额外的压力损失。通道几何形状对总体压降也有着显著影响。扩缩比的增大使得流体在通道内的流速变化更加剧烈，在扩大部分，流体速度降低，压力升高；在缩大部分，流体速度增大，压力降低。这种大幅度的速度和压力变化会导致压力损失增加，从而使总体压降增大。周期长度的变化也会影响总体压降，较短的周期长度意味着流体更频繁地经历扩缩过程，加剧了流体的扰动，增加了压力损失，使得总体压降升高；而较长的周期长度相对减少了流体的扰动，总体压降相对较小。通过数值模拟和实验研究，深入分析了不同流量和通道几何形状下的总体压降特性。结果表明，总体压降与流量之间存在非线性关系，且随着扩缩比和周期长度的变化而变化。这些研究结果为微通道的设计和优化提供了重要的参考依据，在实际应用中，可以根据具体的流量需求和压降限制，合理选择通道的几何形状，以降低总体压降，提高微通道的流动性能。3.3.2局部压降特性除了总体压降特性，微通道局部区域的压降变化也不容忽视，它对理解微通道内的流动机制具有重要意义。在周期性扩缩微通道中，局部区域的压降变化主要集中在扩缩部分和拐角处。在扩大部分，流体的流速降低，根据伯努利方程，流速的降低会导致压力升高，形成局部高压区。这是因为流体的动能部分转化为压力能，使得压力升高。在缩大部分，流体的流速增大，压力降低，形成局部低压区，此时压力能转化为动能。拐角处也是局部压降变化较大的区域。当流体流经拐角时，由于流动方向的突然改变，流体受到强烈的扰动，形成漩涡和二次流。这些复杂的流动现象导致流体的能量损失增加，从而使拐角处的压降显著增大。漩涡和二次流的形成会消耗流体的动能，使得压力降低，同时它们还会加剧流体与壁面之间的摩擦，进一步增加压力损失。通过数值模拟和实验测量，详细分析了局部区域的压降变化规律。结果表明，局部压降与局部结构密切相关，不同的扩缩比和拐角形状会导致局部压降的差异。这些研究结果对于优化微通道的局部结构具有重要指导意义，在设计微通道时，可以通过合理设计扩缩部分和拐角的形状和尺寸，减少局部压降，提高微通道的整体性能。3.4传热特性分析3.4.1总体传热特性周期性扩缩微通道的总体传热性能是评估其热管理能力的关键指标，它与雷诺数、普朗特数等参数密切相关。雷诺数是衡量流体惯性力与粘性力相对大小的无量纲数，它反映了流体的流动状态。当雷诺数增大时，流体的流速增加，对流换热增强，使得总体传热系数增大。这是因为流速的增加促进了流体与壁面之间的热量传递，增强了对流换热效果。普朗特数是流体动量扩散系数与热量扩散系数的比值，它反映了流体的物性对传热的影响。对于不同的流体，普朗特数不同，其传热性能也会有所差异。在周期性扩缩微通道中，普朗特数的变化会影响传热系数的大小。当普朗特数较大时，流体的粘性较大，热量扩散相对较慢，传热系数相对较小；而当普朗特数较小时，流体的粘性较小，热量扩散较快，传热系数相对较大。通过数值模拟和实验研究，深入分析了总体传热系数与雷诺数、普朗特数之间的关系。结果表明，总体传热系数随着雷诺数的增大而增大，且在一定范围内，与雷诺数的对数呈线性关系。同时，总体传热系数还与普朗特数有关，不同的普朗特数会导致传热系数的变化。这些研究结果为微通道的热设计提供了重要的理论依据，在实际应用中，可以根据具体的工况条件，选择合适的流体和通道结构，以提高微通道的总体传热性能。3.4.2局部传热特性研究微通道局部区域的传热特性对于深入理解其传热机制和优化传热性能具有重要意义。在周期性扩缩微通道中，局部区域的传热特性受到局部结构的显著影响。在扩缩部分，由于流体的流速和温度分布发生变化，传热特性也会相应改变。在扩大部分，流体流速降低，停留时间增加，使得流体与壁面之间的热量传递更加充分，传热系数增大；而在缩大部分，流体流速增大，传热系数相对较小。这是因为流速的变化会影响流体与壁面之间的对流换热强度，进而影响传热系数。拐角处也是局部传热特性变化较大的区域。由于拐角处存在漩涡和二次流，这些复杂的流动现象增强了流体的混合和扰动，使得热量传递更加均匀，传热系数增大。漩涡和二次流的存在使得流体在拐角处的流动更加紊乱，增加了流体与壁面之间的接触面积和换热机会，从而提高了传热系数。通过数值模拟和实验测量，详细分析了局部区域的传热特性。结果表明，局部结构对传热具有强化或削弱作用，合理设计局部结构可以有效提高微通道的传热性能。在实际应用中，可以通过优化扩缩部分和拐角的形状和尺寸，利用局部结构对传热的强化作用，提高微通道的整体传热效率。3.5热阻特性分析热阻是衡量微通道传热性能的重要指标，它反映了热量传递过程中的阻力大小。热阻与流量、传热系数等因素密切相关，通过分析热阻特性，可以评估微通道的传热性能四、周期型扩缩微通道结构参数优化4.1扇形凹穴型微通道结构尺寸优化4.1.1基本计算单元及优化的尺寸扇形凹穴型微通道的基本计算单元如图3所示，主要由主流通道和周期性布置的扇形凹穴组成。主流通道宽度为W_1，长度为L_1，扇形凹穴的宽度为W_2，长度为L_2，凹穴深度为H。图3扇形凹穴型微通道基本计算单元示意图通过数值模拟和实验研究相结合的方法对其尺寸进行优化。数值模拟采用有限体积法对控制方程进行离散，选用合适的湍流模型（如k-\epsilon模型）来模拟流体的湍流流动，设置入口速度边界条件和出口压力边界条件，壁面采用无滑移边界条件。实验中，利用微加工工艺制作不同尺寸的扇形凹穴型微通道实验件，搭建实验平台测量不同工况下的流动阻力和传热性能。研究发现，当W_1=0.5mm，W_2=0.2mm，L_1=2mm，L_2=1mm，H=0.3mm时，在保证一定传热性能的前提下，微通道的流动阻力相对较小，综合性能较优。此时，流体在凹穴内形成的漩涡能够有效地增强流体与壁面之间的传热传质，同时又不会导致过大的压力损失。4.1.2L_1/L_2参数优化L_1/L_2参数对微通道内的流动与传热性能有着显著影响。当L_1/L_2较小时，扇形凹穴的数量相对较多，凹穴之间的间隔较小，流体在凹穴内的停留时间较短，形成的漩涡强度较弱，对边界层的破坏作用有限，传热性能提升不明显。随着L_1/L_2的增大，扇形凹穴的数量减少，间隔增大，流体在凹穴内能够充分发展，形成较强的漩涡，边界层被有效破坏，传热性能得到显著提升。然而，当L_1/L_2过大时，凹穴之间的主流通道长度过长，流体在主流通道内的边界层会逐渐恢复，导致传热性能下降，同时流动阻力也会因通道长度的增加而增大。通过数值模拟和实验测量，分析不同L_1/L_2值下微通道的努塞尔数（Nu）和摩擦系数（f）。以Nu与f的综合性能评价指标PEC=(Nu/Nu_0)/(f/f_0)^{1/3}（其中Nu_0和f_0为基准通道的努塞尔数和摩擦系数）为优化目标，确定最优的L_1/L_2参数范围。结果表明，当L_1/L_2在2-3之间时，微通道具有较好的综合性能，此时传热性能得到明显增强，而流动阻力的增加在可接受范围内。4.1.3W_2/W_1参数优化W_2/W_1参数反映了扇形凹穴与主流通道的宽度比例关系，对微通道的性能有着重要影响。当W_2/W_1较小时，扇形凹穴相对较窄，流体在凹穴内的流速较高，形成的漩涡较为集中，但漩涡的作用范围有限，对传热的强化效果不显著。随着W_2/W_1的增大，扇形凹穴变宽，流体在凹穴内的流速降低，漩涡的作用范围扩大，能够更有效地增强流体的混合和传热传质，传热性能得到明显提升。然而，当W_2/W_1过大时，凹穴宽度过大可能会导致流体在凹穴内形成较大的回流区域，增加流动阻力，同时也可能会影响主流通道内的流体流动稳定性，导致传热性能下降。通过数值模拟和实验研究，分析不同W_2/W_1值下微通道的传热系数和压降。以传热系数与压降的比值作为优化目标，寻找使传热性能最佳的W_2/W_1参数值。研究结果表明，当W_2/W_1在0.3-0.4之间时，微通道的传热性能最佳，此时在一定的压降条件下能够实现较高的传热系数，有效地提高了微通道的换热效率。4.2三角凹穴型微通道结构尺寸优化4.2.1基本计算单元及优化的尺寸三角凹穴型微通道的基本计算单元主要由直通道和周期性布置的三角凹穴组成，如图4所示。直通道宽度为W，长度为L，三角凹穴的底边长度为a，高度为h，相邻凹穴之间的间距为s。图4三角凹穴型微通道基本计算单元示意图在确定优化尺寸时，采用数值模拟与实验相结合的方式。数值模拟基于CFD软件，采用有限体积法离散控制方程，选用合适的湍流模型（如k-\omegaSST模型）以准确模拟复杂的流动情况。实验中，利用光刻、刻蚀等微加工技术制作不同尺寸的三角凹穴型微通道芯片，搭建实验平台测量不同工况下的流动与传热参数。通过大量的模拟和实验数据分析，当W=0.6mm，L=3mm，a=0.3mm，h=0.15mm，s=0.5mm时，微通道在传热性能和流动阻力之间取得较好的平衡，综合性能较为理想。此时，三角凹穴能够有效地破坏边界层，增强流体的扰动，提高传热效率，同时流动阻力也在合理范围内，有利于实际应用。4.2.2\alpha参数优化\alpha参数定义为三角凹穴的顶角角度，它对微通道内的流动和传热有着重要影响。当\alpha较小时，三角凹穴相对较尖锐，流体流经凹穴时，在凹穴附近形成的漩涡较小且强度较弱，对边界层的破坏作用有限，传热性能提升不明显。随着\alpha的增大，凹穴逐渐变宽，流体在凹穴内的流动更加复杂，形成的漩涡强度和范围增大，能够更有效地破坏边界层，增强流体与壁面之间的传热传质，传热性能得到显著提升。然而，当\alpha过大时，凹穴过于宽阔，流体在凹穴内的流动阻力会显著增加，同时可能会导致主流通道内的流体流动出现较大的波动，影响传热性能的稳定性。通过数值模拟和实验测量，分析不同\alpha值下微通道的努塞尔数和摩擦系数。以努塞尔数与摩擦系数的综合性能评价指标PEC=(Nu/Nu_0)/(f/f_0)^{1/3}为优化目标，确定最优的\alpha参数值。结果表明，当\alpha=60^{\circ}时，微通道具有较好的综合性能，此时传热性能得到明显改善，而流动阻力的增加在可接受范围内，能够实现高效的传热和稳定的流动。4.2.3\beta参数优化\beta参数为三角凹穴底边与直通道宽度的比值，即\beta=a/W，它对微通道的性能有着重要影响。当\beta较小时，三角凹穴的底边相对较窄，凹穴对主流通道内流体的扰动较小，边界层的破坏程度有限，传热性能提升有限。随着\beta的增大，凹穴底边变宽，对主流通道内流体的扰动增强，边界层被更有效地破坏，流体的混合和传热传质得到促进，传热性能得到明显提升。然而，当\beta过大时，凹穴底边过宽可能会导致流体在凹穴内形成较大的回流区域，增加流动阻力，同时也可能会影响主流通道内的流体流动稳定性，导致传热性能下降。通过数值模拟和实验研究，分析不同\beta值下微通道的传热系数和压降。以传热系数与压降的比值作为优化目标，寻找使传热性能最佳的\beta参数值。研究结果表明，当\beta=0.5时，微通道的传热性能最佳，此时在一定的压降条件下能够实现较高的传热系数，有效地提高了微通道的换热效率。在该参数下，凹穴对流体的扰动恰到好处，既能充分增强传热，又能保持较低的流动阻力。4.2.4\gamma参数优化\gamma参数定义为三角凹穴高度与直通道宽度的比值，即\gamma=h/W，它对微通道性能有着显著影响。当\gamma较小时，三角凹穴相对较浅，流体在凹穴内的流动较为平缓，形成的漩涡较弱，对边界层的破坏作用不明显，传热性能提升有限。随着\gamma的增大，凹穴深度增加，流体在凹穴内的流速和流动方向变化更加剧烈，形成的漩涡强度和范围增大，能够更有效地破坏边界层，增强传热传质，传热性能得到显著提升。然而，当\gamma过大时，凹穴过深可能会导致流体在凹穴内的流动阻力急剧增加，同时可能会影响主流通道内的流体流动，使流动变得不稳定，从而降低传热性能。通过数值模拟和实验测量，分析不同\gamma值下微通道的努塞尔数和摩擦系数。以努塞尔数与摩擦系数的综合性能评价指标PEC=(Nu/Nu_0)/(f/f_0)^{1/3}为优化目标，确定最优的\gamma参数值。结果表明，当\gamma=0.25时，微通道具有较好的综合性能，此时传热性能得到明显改善，而流动阻力的增加在可接受范围内，能够在保证一定流动稳定性的前提下实现高效传热。4.2.5\delta参数优化\delta参数表示相邻三角凹穴之间的间距与直通道宽度的比值，即\delta=s/W，它对微通道性能有着重要影响。当\delta较小时，相邻凹穴之间的间距较小，凹穴对流体的扰动相互叠加，可能会导致流体流动过于复杂，流动阻力增大，同时传热性能的提升也可能会受到限制，因为凹穴之间的流体来不及充分发展。随着\delta的增大，凹穴之间的间距增大，流体在相邻凹穴之间有足够的空间进行发展，能够形成稳定的流动和传热过程，传热性能得到改善。然而，当\delta过大时，凹穴对流体的扰动作用减弱，边界层的破坏程度减小，传热性能会下降。通过数值模拟和实验研究，分析不同\delta值下微通道的传热系数和压降。以传热系数与压降的比值作为优化目标，寻找使传热性能最佳的\delta参数值。研究结果表明，当\delta=0.8时，微通道的传热性能最佳，此时在一定的压降条件下能够实现较高的传热系数。在该参数下，凹穴之间的间距既能保证凹穴对流体的有效扰动，又能使流体在相邻凹穴之间充分发展，从而实现高效的传热和稳定的流动。4.3横断扰流微通道热沉结构尺寸优化4.3.1基本计算单元及优化的尺寸横断扰流微通道热沉的基本计算单元由直通道和横断区内的扰流元组成，如图5所示。直通道宽度为W，长度为L，横断区宽度为W_1，扰流元的形状为矩形，长度为l，宽度为w，扰流元与直通道壁面的距离为d，相邻扰流元之间的间距为s。图5横断扰流微通道热沉基本计算单元示意图为确定优化尺寸，采用数值模拟与实验相结合的方法。数值模拟采用有限体积法离散控制方程，选用适当的湍流模型（如雷诺应力模型）来准确模拟复杂的湍流流动。实验中，通过微加工工艺制作不同尺寸的横断扰流微通道热沉实验件，搭建实验平台测量不同工况下的流动与传热性能。经过大量的模拟和实验数据分析，当W=0.8mm，L=5mm，W_1=0.3mm，l=0.2mm，w=0.1mm，d=0.05mm，s=0.15mm时，微通道热沉在传热性能和流动阻力之间达到较好的平衡，综合性能较为理想。此时，扰流元能够有效地扰动横断区的流体，增强流体的混合和传热传质，同时流动阻力也在可接受范围内，有利于实际应用。4.3.2横断扰流微通道内速度场、压力场及温度场分布通过数值模拟和实验测量，对横断扰流微通道内的速度场、压力场及温度场分布进行深入分析。在数值模拟中，利用CFD软件对微通道内的流体流动和传热过程进行模拟，设置合适的边界条件，如入口速度边界条件、出口压力边界条件和壁面无滑移边界条件等。在实验测量中，采用粒子图像测速（PIV）技术测量速度场，利用压力传感器测量压力场，通过红外热成像仪测量温度场。速度场分布显示，在直通道内，流体速度分布较为均匀，而在横断区，由于扰流元的存在，流体速度发生明显变化。扰流元前方的流体速度降低，形成滞止区；扰流元后方的流体速度增大，形成高速射流区，同时在扰流元两侧形成漩涡，漩涡的存在增强了流体的混合和传热传质。压力场分布表明，在扰流元前方，压力升高，形成高压区；在扰流元后方，压力降低，形成低压区，这种压力分布有利于流体的流动和传热。温度场分布显示，在直通道内，流体温度沿流动方向逐渐升高，而在横断区，由于扰流元对流体的扰动，温度分布变得更加均匀，有效提高了传热效率。通过对速度场、压力场和温度场分布的分析，深入了解了横断扰流微通道内的流动和传热机理，为进一步优化微通道结构提供了理论依据。4.3.3矩形扰流元尺寸对流动与换热的影响矩形扰流元的尺寸（长度l和宽度w）对微通道内的流动与换热有着显著影响。当扰流元长度l较小时，扰流元对横断区流体的扰动范围较小，形成的漩涡强度较弱，对传热的强化作用有限。随着l的增大，扰流元对流体的扰动范围扩大，形成的漩涡强度和范围增大，能够更有效地破坏边界层，增强流体与壁面之间的传热传质，传热性能得到显著提升。然而，当l过大时，扰流元对流体的阻力也会增大，导致流动阻力急剧增加，同时可能会影响主流通道内的流体流动稳定性，导致传热性能下降。扰流元宽度w的变化也会对流动与换热产生影响。当w较小时，扰流元对流体的扰动程度较小，传热性能提升不明显。随着w的增大，扰流元对流体的扰动增强，传热性能得到提高。但当w过大时，流动阻力会显著增加，同时可能会导致横断区的流体流动出现堵塞现象，降低传热性能。通过数值模拟和实验研究，分析不同扰流元尺寸下微通道的努塞尔数和摩擦系数。以努塞尔数与摩擦系数的综合性能评价指标PEC=(Nu/Nu_0)/(f/f_0)^{1/3}为优化目标，确定最优的扰流元尺寸。结果表明，当扰流元长度l=0.2mm，宽度w=0.1mm时，微通道具有较好的综合性能，此时传热性能得到明显改善，而流动阻力的增加在可接受范围内，能够实现高效的传热和稳定的流动。4.3.4矩形扰流元位置对流动与换热的影响矩形扰流元在横断区的位置对微通道的流动与换热性能有着重要作用。当扰流元靠近直通道入口时，流体在进入横断区后立即受到扰流元的扰动，能够迅速增强流体的混合和传热传质，但可能会导致入口处的流动阻力较大。随着扰流元向出口方向移动，流体在直通道内有一定的发展距离，进入横断区五、并联硅基周期性扩缩微通道内气液两相流型与压降特性5.1实验参数测量在实验中，对气液两相流参数的精确测量至关重要，它直接关系到对气液两相流型和压降特性的准确理解。对于气体和液体的流量测量，分别采用了高精度的气体质量流量控制器和液体注射泵。气体质量流量控制器利用热式质量流量测量原理，通过测量气体对加热元件的冷却效应来确定气体的质量流量。其量程范围为0-1000sccm，精度高达±1%FS，能够提供稳定且精确的气体流量控制。液体注射泵则采用容积式计量原理，通过高精度的步进电机驱动活塞，精确控制液体的注射量，从而实现对液体流量的精确调节。其流量范围为0-10mL/min，精度可达±0.001mL/min，确保了液体流量的高精度测量。压力测量选用了精度为±0.01kPa的高精度压力传感器。该传感器基于压阻效应，通过测量压力作用下敏感元件的电阻变化来确定压力值。在微通道的进出口以及不同位置处布置压力传感器，能够实时监测气液两相流在微通道内的压力变化，获取准确的压降数据。为了保证压力测量的准确性，在实验前对压力传感器进行了严格的校准，采用高精度的压力校准仪对传感器进行标定，确保其测量误差在允许范围内。为了直观地观察气液两相流的流型，采用了高速摄像仪。该高速摄像仪帧率可达10000fps，配合高分辨率显微镜镜头，能够清晰捕捉微通道内气液两相流的瞬间流型变化。在实验过程中，通过合理调整高速摄像仪的拍摄角度、焦距和曝光时间等参数，确保拍摄到的流型图像清晰、准确。同时，利用专业的图像处理软件对拍摄的图像进行分析，能够准确识别不同的流型，并测量流型的相关参数，如气泡尺寸、形状、分布等。通过这些精确的实验参数测量方法和仪器，为深入研究并联硅基周期性扩缩微通道内气液两相流型与压降特性提供了可靠的数据支持。5.2并联硅基周期性扩缩微通道热沉内氮气/水气液两相流型5.2.1气液两相流型描述在并联硅基周期性扩缩微通道热沉内，通过高速摄像仪观察到了丰富多样的氮气/水气液两相流型，主要包括泡状流、弹状流、块状流和环状流等，每种流型都具有独特的特征。泡状流是在气相体积分数较低时出现的一种流型。此时，气相以离散的小气泡形式均匀分布在连续的液相中，气泡的尺寸相对较小且分布较为均匀。由于气泡的存在，液相的连续性受到一定程度的干扰，但整体流动相对较为稳定。气泡在液相中受到浮力、粘性力和表面张力等多种力的作用，其运动轨迹呈现出一定的随机性。在微通道的不同位置，气泡的分布和运动状态可能会有所不同，靠近壁面处的气泡由于受到壁面的影响，运动速度相对较慢，而在通道中心处的气泡则运动速度较快。随着气相体积分数的增加，流型逐渐转变为弹状流。弹状流的特征是较大的气泡（弹状气泡）在液相中依次排列，气泡之间夹杂着少量的液相。弹状气泡的长度通常与微通道的水力直径相当，其头部呈球形，尾部扁平，形如炮弹。在弹状流中，弹状气泡的运动对液相的流动产生较大的影响，导致液相的流速分布不均匀。弹状气泡在上升过程中，会带动周围的液相一起运动，形成局部的流速增大区域。同时，弹状气泡与壁面之间的液层也会受到气泡运动的影响，产生一定的剪切应力。当气相体积分数进一步增加时，块状流出现。在块状流中，气相形成较大的块状结构，液相则分布在块状气相的周围。块状流的流动稳定性较差，相界面复杂，容易出现局部的漩涡和湍流现象。由于相界面的复杂性和流动的不稳定性，块状流中的气液混合更加剧烈，传热传质效率相对较高。然而，这种复杂的流动状态也给流动阻力的计算和压降特性的研究带来了较大的困难。在高气速和较高气相体积分数下，环状流形成。环状流的特点是液相在微通道壁面上形成一层连续的液膜，气相则在通道中心形成核心流。在环状流中，液膜的厚度和稳定性对流动特性有着重要影响。液膜的厚度会随着气相流速和液相流量的变化而变化，当气相流速较大时，液膜会被吹薄，甚至出现破裂的情况。液膜与气相之间的相互作用也会影响流动的稳定性，液膜的波动和破裂可能会导致局部的压力波动和流动阻力的增加。5.2.2气液两相流型分布图为了更直观地展示不同工况下流型的变化规律，绘制了气液两相流型分布图，如图6所示。该图以气相表观速度和液相表观速度为坐标轴，将不同流型的出现区域清晰地划分出来。图6气液两相流型分布图在低气相表观速度和低液相表观速度区域，主要出现泡状流。随着气相表观速度的增加，当达到一定值时，流型开始从泡状流转变为弹状流，转变的临界条件与微通道的几何尺寸、流体物性等因素有关。继续增加气相表观速度，弹状流逐渐过渡为块状流，块状流区域位于弹状流和环状流之间。当气相表观速度和液相表观速度都较高时，环状流出现。从流型分布图中可以看出，流型的转变是一个连续的过程，不同流型之间的过渡区域并非截然分明，而是存在一定的重叠。这是因为流型的转变受到多种因素的综合影响，在过渡区域内，不同流型的特征可能会同时存在。通过流型分布图，能够直观地了解不同工况下流型的变化趋势，为进一步研究气液两相流的流动特性和压降特性提供了重要的参考依据。在实际应用中，可以根据流型分布图，合理选择气液流量，以获得所需的流型，从而优化微通道热沉的性能。5.3并联硅基周期性扩缩微通道热沉内氮气/水气液压降特性分析5.3.1微通道内气液两相压降特性微通道内气液两相的压降特性是研究气液两相流的重要内容，它与流量、流型等因素密切相关。随着气液流量的增加，微通道内的流速增大，惯性力增强，导致压降呈现上升的趋势。在不同的流型下，压降特性也存在明显差异。在泡状流中，由于气泡尺寸较小且分布均匀，对流体的扰动相对较小，压降主要由液相的粘性阻力和气泡与液相之间的摩擦阻力引起。随着气相体积分数的增加，气泡数量增多，气泡与液相之间的相互作用增强，导致压降逐渐增大。弹状流中，弹状气泡的存在使得流体的流动变得更加复杂，弹状气泡的运动不仅会带动周围液相的流动，还会在气泡与壁面之间形成局部的流速梯度，从而增加了流动阻力，使得压降明显增大。弹状气泡的长度、速度以及它们之间的间距等因素都会影响压降的大小。块状流中，由于气相形成较大的块状结构，相界面复杂，气液混合剧烈，流动阻力进一步增大，压降也随之增大。在块状流中，局部的漩涡和湍流现象会消耗更多的能量，导致压降急剧增加。环状流中，液相在壁面上形成液膜，气相在中心形成核心流。液膜的厚度和稳定性对压降有重要影响。当液膜较薄且稳定时，压降相对较小；但当液膜出现波动或破裂时，会增加气液之间的摩擦和能量损失，导致压降增大。5.3.2实验结果同均相模型计算值的比较为了评估均相模型在预测并联硅基周期性扩缩微通道内气液两相压降方面的适用性，将实验结果与均相模型计算值进行了对比。均相模型假设气液两相均匀混合，将气液混合物视为一种具有平均物性的单相流体来处理。通过计算发现，均相模型在某些工况下能够较好地预测压降趋势，但在具体数值上与实验结果存在一定的偏差。在低气相体积分数和低流速的工况下，均相模型的计算值与实验结果较为接近，这是因为此时气液两相的混合较为均匀，均相模型的假设基本成立。然而，随着气相体积分数的增加和流速的增大，气液两相的相互作用变得更加复杂，均相模型无法准确描述这种复杂的流动现象，导致计算值与实验结果的偏差逐渐增大。进一步分析偏差产生的原因，主要是均相模型忽略了气液两相间的滑移现象和相界面的影响。在实际的气液两相流中，气相和液相的速度并不相同，存在滑移现象，这会导致动量传递和能量损失的变化，而均相模型未能考虑这一因素。此外，相界面的存在使得气液之间的相互作用更加复杂，均相模型无法准确描述相界面处的物理过程，从而影响了压降预测的准确性。5.3.3实验结果同分相模型计算值的比较分相模型将气液两相分别进行考虑，分别计算气相和液相的压降，然后将两者相加得到总的压降。为了评估分相模型在本研究中的准确性，将实验结果与分相模型计算值进行了比较。实验结果表明，分相模型在一定程度上能够更准确地预测气液两相的压降。在不同的流型下，分相模型的计算值与实验结果的吻合度相对较高，尤其是在高气相体积分数和高流速的工况下，分相模型的优势更加明显。这是因为分相模型考虑了气液两相间的滑移现象和相界面的影响，能够更真实地描述气液两相的流动特性。然而，分相模型也存在一定的局限性。在某些复杂的流型下，如块状流和环状流，由于相界面的复杂性和流动的不稳定性，分相模型仍然难以准确预测压降。此外，分相模型中一些参数的确定较为困难，如气液相间的摩擦系数等，这些参数的不确定性也会影响模型的预测精度。5.4小结本部分通过实验研究了并联硅基周期性扩缩微通道热沉内氮气/水气液两相流型与压降特性。采用高精度的测量仪器，准确测量了气液两相流的参数。观察到了泡状流、弹状流、块状流和环状流等多种流型，并对其特征进行了详细描述，绘制了气液两相流型分布图，清晰展示了不同工况下流型的变化规律。分析了微通道内气液两相的压降特性，发现压降与流量、流型密切相关。将实验结果与均相模型和分相模型计算值进行比较，结果表明均相模型在某些工况下存在较大偏差，分相模型虽然在一定程度上能够更准确地预测压降，但在复杂流型下仍存在局限性。这些研究结果为深入理解并联硅基周期性扩缩微通道内气液两相流的特性提供了重要的实验依据，也为相关理论模型的改进和完善提供了参考。六、PDMS单微通道内气液两相流型及其演变特性6.1PDMS单微通道内气液两相流型及其演变特性的实验研究6.1.1PDMS单通道内气液两相流型描述在PDMS单通道内，通过高速摄像仪对气液两相流型进行了细致的观察和记录，识别出了多种典型的流型，每种流型都具有独特的形态和流动特征。泡状流是在低气相流速和高液相流速条件下常见的流型。此时，气相以微小气泡的形式均匀分散在连续的液相中，气泡的尺寸相对较小，一般在几十微米到几百微米之间。这些小气泡在液相的带动下，沿着通道轴线方向流动，其运动轨迹相对较为规则，主要受到液相的粘性力和浮力的作用。由于气泡的存在，液相的连续性受到一定程度的干扰，但整体流动仍然较为稳定，流场中的速度分布相对均匀。随着气相流速的逐渐增加，流型逐渐转变为弹状流。弹状流的显著特征是气泡尺寸增大，形成较大的弹状气泡，其长度通常与微通道的水力直径相当，甚至更长。这些弹状气泡在液相中依次排列，气泡之间夹杂着一定量的液相。弹状气泡的头部呈球形，尾部扁平，在流动过程中，弹状气泡会带动周围的液相一起运动，形成局部的流速增大区域，导致液相的流速分布变得不均匀。同时，弹状气泡与壁面之间的液层也会受到气泡运动的影响，产生一定的剪切应力，使得液层的流动状态变得复杂。当气相流速进一步增大时，块状流出现。在块状流中，气相聚集成较大的块状结构，液相则分布在块状气相的周围。这种流型的流动稳定性较差，相界面复杂且不规则，容易出现局部的漩涡和湍流现象。由于相界面的复杂性和流动的不稳定性，块状流中的气液混合更加剧烈，热量和质量传递效率相对较高。然而，这种复杂的流动状态也给流动阻力的计算和压降特性的研究带来了较大的困难，因为在块状流中，气液之间的相互作用更加复杂，难以用简单的理论模型进行描述。在高气速和较高气相体积分数的条件下，环状流形成。环状流的特点是液相在微通道壁面上形成一层连续的液膜，气相则在通道中心形成核心流。液膜的厚度和稳定性对环状流的流动特性有着重要影响。液膜的厚度会随着气相流速和液相流量的变化而变化，当气相流速较大时，液膜会被吹薄，甚至出现破裂的情况。液膜与气相之间的相互作用也会影响流动的稳定性，液膜的波动和破裂可能会导致局部的压力波动和流动阻力的增加。此外，在环状流中，气相核心流中的速度分布也不均匀，中心区域的速度较高，靠近液膜的区域速度较低。6.1.2PDMS单通道内气液两相流型分布图为了直观地展示不同工况下气液两相流型的分布和演变规律，绘制了PDMS单通道内气液两相流型分布图，如图7所示。该图以气相表观速度（j_g）和液相表观速度（j_l）为坐标轴，将不同流型的出现区域清晰地划分出来。图7PDMS单通道内气液两相流型分布图在低气相表观速度和高液相表观速度的区域，主要出现泡状流。随着气相表观速度的逐渐增加，当达到一定阈值时，流型开始从泡状流向弹状流转变，这个转变的临界条件与微通道的几何尺寸、流体物性以及表面润湿性等因素密切相关。继续增加气相表观速度，弹状流逐渐过渡为块状流，块状流区域位于弹状流和环状流之间，其范围受到多种因素的影响，如气相和液相的流量比、微通道的粗糙度等。当气相表观速度和液相表观速度都较高时，环状流出现，环状流的形成需要一定的气相流速和气相体积分数条件，以确保气相能够在通道中心形成稳定的核心流，同时液相能够在壁面上形成连续的液膜。从流型分布图中可以看出，流型的转变是一个连续的过程，不同流型之间的过渡区域并非截然分明，而是存在一定的重叠。这是因为流型的转变受到多种因素的综合影响，在过渡区域内，不同流型的特征可能会同时存在。例如，在泡状流向弹状流转变的过渡区域，可能会同时观察到小气泡和较大的弹状气泡；在弹状流向块状流转变的过渡区域，弹状气泡可能会逐渐聚集成块状结构。通过流型分布图，能够直观地了解不同工况下流型的变化趋势，为进一步研究气液两相流的流动特性和传热特性提供了重要的参考依据。在实际应用中，可以根据流型分布图，合理选择气液流量，以获得所需的流型，从而优化微通道的性能。6.2PDMS单微通道内气液两相流型及其演变特性的数值模拟6.2.1多相流模型选择在对PDMS单微通道内气液两相流型及其演变特性进行数值模拟时，选择合适的多相流模型至关重要。经过综合考虑，选用了VOF（VolumeofFluid）模型。VOF模型是一种处理多相流问题的有效方法，特别适用于模拟没有相互穿插的多相流情况，在PDMS单微通道内气液两相流中，气液之间存在明显的分界面，且不会相互穿透，VOF模型的这一特点与实际情况相符。该模型通过求解各相的体积分数方程来确定相界面的位置和形状。在每个计算单元中，定义气相和液相的体积分数\alpha_g和\alpha_l，满足\alpha_g+\alpha_l=1。当\alpha_g=0时，表示该单元中只有液相；当\alpha_g=1时，表示该单元中只有气相；当0<\alpha_g<1时，表示该单元中存在气液交界面。通过追踪体积分数的变化，可以准确地捕捉气液两相的界面位置和形状，从而模拟不同流型下的气液两相流动特性。与其他多相流模型相比，VOF模型具有计算效率高、内存需求小等优点。在处理复杂的几何形状和动态界面问题时，VOF模型能够有效地减少计算量，提高模拟的速度和精度。同时，VOF模型在处理分层流动或自由表面流方面具有独特的优势，这与PDMS单微通道内气液两相流的特点相契合，能够准确地模拟气液两相在微通道内的流动行为和流型演变过程。6.2.2流体动力学建模为了准确描述PDMS单微通道内气液两相流的流体动力学行为，建立了相应的数学模型。该模型基于质量守恒、动量守恒和能量守恒方程，考虑了气液两相间的相互作用以及壁面的影响。质量守恒方程描述了气液两相在微通道内的质量变化情况，对于气相和液相分别有：\frac{\partial(\alpha_g\rho_g)}{\partialt}+\nabla\cdot(\alpha_g\rho_g\vec{v}_g)=0\frac{\partial(\alpha_l\rho_l)}{\partialt}