微热管及其传热理论分析

1、微热管及其传热理论分析摘要：随着微电子制造技术的快速发展，微热管在航天器热控系统、微电子元器件散热等领域中有着广泛的应用。微热管是利用密封在管内工质相变进行热量传输的器件，具有体积小、重量轻、传热效率高、成本低、易于集成、无需外加动力等显著优点，能有效解决目前微小型器件和芯片的散热问题，具有广泛的应用前景。作者综述了微热管的发展与当前研究现状，详细介绍了微热管的工作原理，并指出微热管与常规微热管的区别，对槽道式平板微热管进行理论分析，最后展望了该领域的未来研究方向。关键词：微热管，工作原理，平板微热管，引言随着电子科技技术的进步，许多电子产品向着高性能化、高功率化和小型化方向发展，同时产品的高

2、集成度使其散热空间更为狭小，导致了电子元器件单位面积的热量急剧上升，如高性能微处理器的热流密度已达到100W/cm21。元器件的温度每升高10C,系统的可靠性降低50%2,所以必须采用高效的传热技术对电子元器件进行散热。微热管是一种利用相变传热的高效传热元件，其导热能力大大超越了铜、铝材料的空气强制对流散热方案3-4,因此，具有高导热率、良好的等温性，以及结构简单等优点1,5的微热管成为微电子散热领域的关键元件，并广泛应用于各种电子产品。其中平板微热管由于其良好的蒸发吸热特性和形状易于与芯片贴合等优点被越来越多地应用于高效散热中。而微热管或热管内微结构具有强化传热传质的作用，引起研究者越来越多

3、的关注。1 .微热管的发展与国内外研究现状微热管是利用密封在管内工质相变进行热量传输的器件，具有体积小、重量轻、传热效率高、成本低、易于集成、无需外加动力等显著优点，能有效解决目前微小型器件和芯片的散热问题，具有广泛的应用前景。1944年Gaugler第一次提出了热管的工作原理；1963年美国应用物理杂志报道了世界上的第一根热管；1984年Cotter等人提出了热管微型化的设想，为微热管的研究开辟了道路；1984年，T.PCotter在第五届国际热管会议上首次提出了微热管的概念，并指出微热管在用于电子芯片冷却散热领域具有广阔的应用前景。关于微热管的研究，最初集中于几个厘米长，工质通道横截面为带

4、有尖角区域的图形，通道的水力半径在10科m100科m的单根微热管。工质回流主要靠的是横截面尖角区域所形成的毛细力。这种单根微热管主要应用在传输热量不是很大，但要求温度分布均匀稳定的领域。随后微热管的研究分别从实验研究和理论研究两方面逐步展开，研究结果均体现出这一传热元件相比其它传热手段具有效率高而无需外加动力的优点。而关于微热管结构的研究也从单根微热管逐步发展到微热管阵列，即在固体基板上开出一簇簇微型槽道，这样的方式大大提高了微热管的传热能力，但这只是单根微热管的一种简单的并列组合。进一步的改进是具有连通蒸汽腔的平板微热管。平板微热管通过连通蒸汽腔降低了气液界面高速对流产生的界面摩擦力，使热管

5、的传热能力进一步提高，从而成为目前微热管领域的研究热点。2 .微热管工作原理图l所示为微热管工作原理示意图。根据微热管内部蒸汽流动情况，沿其轴向可分为蒸发段、绝热段和冷凝段。从结构上分析，微热管包括管壳、毛细吸液芯和工作介质（液流）。为降低热阻和工作介质沸点，提高微热管工作效率，管壳内部需保持一定的真空度。在微热管工作时，工作介质在蒸发段吸收热源热量发生相变，蒸汽流经过绝热段到达冷凝段释放热量并凝结为液体，冷凝液流在毛细吸液芯的毛细作用下回流到蒸发段，如此循环下去，微热管不断将热量从热源带走。tn，汽流jii一蒸发段绝热段冷凝段图1微热管工作原理示意图热量由一端传向另一端，在整个热量传递的过程

6、中主要包含以下六个相互关联的过程:(1)热量从热源通过管壁和充满工作液体的管芯传递到液汽分界面；(2)液体在蒸发段内的液汽分界面上蒸发；(3)蒸汽腔内的蒸汽从蒸发段流到冷凝段；(4)蒸汽在冷凝段内的汽液分界面上凝结；(5)热量从汽液分界面通过管芯、液体和管壁传给冷源；(6)冷凝后的工作介质液体在管芯内毛细力作用下回到蒸发段。3 .微热管与常规热管的区别微热管的结构及工作原理与常规热管类似，最大的区别是常规热管内部通常存在专门提供毛细力以供工质回流的毛细吸液芯；而微热管则主要是通过沟道尖角区完成工质的回流6。微热管稳态工作时，微热管蒸发段受到外界加热，热量通过热管管壁及液态工质传递到汽液分界面，

7、使工作液体在蒸发段内的汽液分界面上蒸发。由于液体蒸发，蒸汽腔内产生压差，而蒸汽正是在这压差的作用下由蒸发段流向冷凝段并在冷凝段内的汽液分界面上凝结，释放出热量。热量通过液态工质和管壁传给冷源，最后由于热管内腔尖角区域的毛细作用使冷凝后的工作液体回流到蒸发段7。图2给出了常见的微热管沟道截面结构。如前所述，微热管与常规热管最大的差别在于微热管依靠沟道尖角区提供毛细力，理论分析发现只要是非圆形的截面都能提供或大或小的毛细力。微热管发展之初，常用的沟道结构主要是简单的三角形和矩形结构，如2图中（a）、（b）、（c）所示。随着机械加工和其他各种加工技术的发展，一些比较不规则的结构也进入研究之列，如图2

8、-2中（d）、（e）、（f）所示。在实际应用中采用哪种沟道结构主要由应用场合、热管材料、加工方法等共同决定。4 .微热管的传热极限微平板热管通过相变换热的工作机理使其具有很强的传热能力，但也不能无限制的增大其热载荷，热管工作过程中也受到很多因素影响存在一定的传热极限。传热极限的大小是由热管的内部流道结构、外部形状特征、管壳材料、加工方法、工质特性、工作温度等共同决定的。而多种影响因素分别产生的传热极限中最小值即决定了微热管的最大传热能力8。通过国内外大量研究发现目前对微热管性能影响最大的几种常见传热极限为：工质回流极限、工质沸腾极限、气体冷凝极限、对流携带极限、蒸汽压力极限等。现对其逐一介绍：

9、（1）工质回流极限随着热管蒸发段热载荷的不断增加，液体由冷凝段流向蒸发段的速度不能无限制增大，会受到毛细结构产生的驱动力的限制，为此称其为工质回流极限。此时当热量超过此极限值，由毛细力作用回流的冷凝工质便满足不了蒸发段的需求，结果会在热管蒸发段发生干涸现象，蒸发段温度急剧升高，严重时就会出现烧损等情况。微热管与常规热管相似，工质回流极限依然是影响其传热能力的主要问题，而且对于微热管而言，靠沟道尖角区域回流工质时更容易产生工质回流不足的现象，因此工质回流极限对微热管的影响更为突出。（2）工质沸腾极限随着热管蒸发段热量的增加，蒸发段径向热流密度或热管内壁液体温度过高，会在与管壁接触液体的核化中心产

10、生沸腾气泡，而热管的吸液芯结构会阻碍气泡脱离吸液芯结构，因此气泡停留在液态工质回流的路径上会阻碍工质的回流。若产生气泡较多时，会在工质和管壁之间形成蒸汽层，使气体和管壁之间的热阻迅速增大，也会导致管壁温度快速增高，造成蒸发段干涸或者热管烧损的现象。通过研究发生工质沸腾极限时气体的临界热流和临界过热度可以发现微热管的工质沸腾极限由热管工作条件和沸腾气泡产生的临界半径共同决定。（3）气体冷凝极限微热管冷凝段的冷却能力可能会受到两种情况的限制：一种是充液量过大，微热管的实际工况与最佳工况有偏差，过量的工质在微热管冷凝段堆积，使换热面积变小，从而使微热管的性能降低或工作失效；另一种是不凝气体的存在，影

11、响了冷凝段的冷凝效率。(4)对流携带极限随着微热管热载荷的不断增大，工质流道中的液态和气态工质的流速都会不断增大，会在气液交界面上产生对热管影响较大对流剪切力。当这个剪切力增大到一定程度时，气流会将部分液态工质卷入气流中，从而减少了液态工质回流到蒸发段的量。这种现象发生时会大大降低热管的传热能力，使其过早达到传热极限。(5)蒸汽压力极限又称之为粘滞极限，当热管的工作温度较低时，内部气态工质较为稀薄，会出现分子不连续现象，而由于管壁与气体分子的粘滞力作用会使气态工质在冷凝段末端出现压降变小的状况，使其内部整蒸汽的运动受到限制，导致热管的传热能力受到一定的影响。5 .梢道式平板微热管理论分析本节将

12、以微槽道平板微热管为例，如图3所示，利用一维稳态模型，对影响其性能的几个结构性问题做定量分析，指出其性能可进一步提高的潜在范围，并提出通过结构改进的方案和意义。图？三曾形截面平板旅缺管咕掏示意图5.1 微槽道平板微热管的一维稳态模型分析其内部工作过程的一维稳态模型是目前最为常用的用于预测平板微热管传热能力、一种理论分析方法，其计算结果能准确反应微热管的工作状态平板微热管的一维稳态模型通常利用有限容积法建立，即将所研究对象。通过某种方法划分为若干个有限容积单元，在每个有限容积单元内选取控制节点；再分别对每个有限容积单元列出守恒方程，并对其进行离散化处理，得出离散化方程组；最后通过数值方法将其求解

13、。5.2模型的假设条件平板微热管工作时内部发生复杂的流动、相变、传热过程，一维稳态模型只能将其物理过程简化处理。因此，参照国内外相关理论研究，在保证模型计算精度的前提下，为使模型简化需做出如下假设：(1)热管工作在稳定的条件下，不考虑微热管的瞬态和不稳定状态，因为只有在热管处于稳态时，后文的质量守恒方程才是准确有意义的；(2)气体和液体工质的密度为恒定值，均为不可压缩流体，因为在对流体建立动量守恒方程时要求流体不可压缩，只有这样动量守恒方程才准确有意义；(3)液态工质达到稳定工作状态时其液体表面弯月面半径仅沿轴向变化，因为是一维模型，无法考虑其在其它方向的变化；(4)同理工质的物理参量也只能沿

14、轴向发生变化，例如压强、流速等均取其在流体横街面上的平均值。5.2.1模型的控制方程在建立控制方程之前对微热管进行有限元划分，如图4所示为微热管内部某槽道的剖面图，其近似给出了液体与气体稳态工作时的分布图，沿热管轴向将划分为若干个长度为dz的有限元体积，图4中也给出了单个有限元体积的示意图，下文中将对单个有限元体积建立相应的控制方程。图4平板放热管有限元划分示意图(1)质量守恒方程对于有限容积单元dz,在其处于稳定状态时满足质量恒定，即单位时间人有限容积单兀dz和流出dz的质量是相等的，如图5所示。>+-1q(iniBBBirniBriBiiKBaBiIt4491!4-V9B11-11I

15、9)¥HVapor$-<:pvj4rur>不uj-t<Pv&n.：.：.1JQ1巧4乜-Lr-i-4)ud-1*dz|图5质量守恒示意图如图5所不，某时刻流入有限容积单兀的泊志工质质量为。.dnA%dz1dzd(A4%)i+dzdi(5-1)某时刻流出该有限容积单元的液态工底质量是；%=/<5-2）单位时间内该有限容积单元增加的质量是：A-=d（4%）次（5.3）az式中制融在某一叶刻流入的液态工质质量一在某一时刻演出的液态工质质量；Am在某一时刻有限容积单元内工质质量的变化量Jr一液态工质的密度?%涌态工后在流体某一横截面位置的平均流速.；4在某一位

16、置液态工质流体的横截面积。因为热管工作在稳定状态，且对于任意一个有限容积单元都应处在稳定状态所以每一个有限容积单元内都应满足质量守恒定律，即单元内工质质量为恒定值，且随着时间的变化该恒定值是不变的。因此对于某一有限容积单元流入质量与流出质量的差值为发生相变的质量。相变的质量可以由m=QXfgh得出，其中Q为吸收或放出的热量，fgh为汽化潜热，m为相变液态工质的质量。当微热管槽道数为N时，则Q=NXmXfgh,那么单根槽道内发生相变的液态工质质量是：C5K)(5-5)=Nxh&八fs因为Am=mc,所以得到平板微热管内部工质的质量守恒方程为："（.4%）_Q工dz-Nxh对于质

17、量守恒方程式（5-5）无论在蒸发段、绝热段、冷凝段其表达式不变，只是在不同位置，输入热量Q的取值不同而已“如果假设无论是蒸发段还是冷海段输入热量都是均匀分布则输入热量的表达式可以写为：£5-6)再对式(5-6)两侧分别对上求导可得:dz04QinTc(5乂)因为Am=mt所以得到平板微热管内部工后的质量守恒方程为：(5-5)对于质量守恒方程式C5-5)无论在蒸发段、绝热段、冷凝段箕表达式不变，只是在不同位置，输入熟量Q的取值不同而已“如果假谩无论是蒸发段还是冷凝股输入热量都是均匀分布则输入热量的表达式可以写为：(2)动量守恒方程本文仍以液态工质流动为主要研究对象。在此前我们已经阐述了

18、微槽道内液相工质流动的工作原理：由于蒸发段与冷凝段毛细半径值的不同所产生的毛细力的作用，使得微槽道内的液体会从冷凝段流回蒸发段，微槽道所提供的毛细力的大小与这两个毛细半径之间的差值是成正比关系。在毛细力的作用下液态工质体由冷凝段流向蒸发段的过程中还受到其它作用力的影响，如重力、压力和摩擦力等。在微槽道内液态工质流动的动量守恒方程式中，关于摩擦力项的部分争议较大，而关于重力和压力的动量项定义都不存在任何分歧。实际上摩擦力项应由槽道壁面产生的摩擦力和反向流动蒸汽与液态工质之间长生的高速对流剪切力两部分组成。而在气态工质的流动过程中，也存在这两种摩擦力影响，其中不可以忽略液态工质流动对蒸汽作用的摩擦

19、力。图6对整个动量关系有详细的描述。Pv4北+吆电空dzdzAM“I空dzdzp同哈粤dz图6有限容积单元内动量关系示意图以z方向为正，对各动量项符号确定之后，列出液态工质流动的动量守恒方程："竺.)+A孚+A°gsin夕=窃A+A.(5-10)dzdz式中:片液态工质在某一位置的压强平为值3方一界面对滂摩擦系数；A某一位置上气液界面的长度3*3一液体与槽道管壁间的摩擦系数；At某一位置液体与槽道管壁间的接触长度。式5-10中各项表达式的含义从左到右分别为：动量项、压力项、重力项、界面剪切力项和管壁摩擦力项。同理，我们也可以列出气态工质流动的动量守恒方程：d(nAu2).d

20、P.a、与-+A会+Apgsin0=-已A-小A“，(5-n)dzdz式中1一气态工质在某一位警的压强平均值j%一方程调节系数,参照以往研究取L3力Q,一界面对流摩擦系数5J气体流动与管壁间的摩擦系数JAv某一位置气体与管壁的接触长度因为管壁摩擦力相对于界面再切力来说是很小的，所以在实际计算中我们将其忽略，那么(5/，(5-11)两式就可以改写成：勺丁,+4塞+3mJ=0A(5-12)dz+A户gsin8=心A(5-H)(3)Laplace-Young方程微热管在工作时，蒸汽在冷凝段冷凝，液体在冷凝段增多，使弯月面半径变大；而液体在蒸发段受热蒸发，液体向尖角区流动，使毛细半径变小。槽道内部的毛

21、细作用就是由这种冷凝段和蒸发段之间的毛细半径之差产生的。弯月面的产生是由于工质表面张力、气提压强和液体压强之间共同作用的结果，Laplace-Young方程给出了弯月面半径和气液压强之间的关系：(5-14)(5-15)a以一必二一rLaphce-Young方程的激分形式如下：肛马(7drdz一r2也式中仃一表面张力心(4)边界条件为求解上述徵分方程，需要确定徜分方程所涉及相关参量的边界条件,微分方程组在蒸发段揣部(的边界条件如下:厂二仁=凡凡=7皿二凸二。血一07%式中rD理月面半径在蒸发段端处的初值;P在工作温度下的饱和蓦汽压口5.2.2模型的求解过程式(5-8),(5-9),(5-12),

22、(5-13),(5-15)和(5-16)给出的边界条件，建立了微热管一维稳态的毛细流动模型。上述模型的方程组建立完毕后可以在MATLAB环境下采用四阶Runga-Kutta值积分算法建立迭代程序，对其进行数值求解。在实际求解过程中，边界条件0r是未知的，因此与之对应的r和lp的初值也是未知的，为能继续求解需要给定0r一个假设的初值，然后通过迭代收敛求解其真实值。在固定某一输入热量且热管正常工作情况下，即没出发生蒸发段干涸时，其蒸发端处最小弯月面半径minr是由微槽道尖角区域的几何结构决定的。而在冷凝端，液态工质的最大弯月面半径maxr,则是由槽道的几何结构参数和工质与管壁间接触角共同决定的，可

23、以由下式给出：因此在蒸发端液态工质的弯月面半径计算方法如下：首先给定蒸发端处的弯月面半径0r一个任意小的初值；然后通过模型计算得出输入热量为inq时的弯月面半径的轴向分布；再不断增加0r并不断的通过模型计算出冷凝端处的弯月面半径，直到其恰好等于maxr,则此时蒸发端处的弯月面半径0r即为真实值。当弯月面初始值0r确定后，各参量沿微热管槽道的轴向分布就可以通过耦合方程组计算得到。需要特别说明的是，不同的inq所对应的0r的值是不相同的。微型多槽道平板热管的最大传热量Q,需要如下步骤计算：首先假设蒸发端的弯月半径0r为minr;通过模型计算得出输入热量为inq时的弯月面半径的轴向分布；然后不断增大

24、inq并通过模型计算出冷凝端处的弯月面半径，直到其恰好等于maxr,则此刻的输入热量inq即为微型多槽道平板热管的最大传热量Q。6.结语本章以常规热管的工作原理为出发点，给出了微热管的定义、工作原理，并对微热管与常规热管的工作原理差异进行了较为深入的说明。此外，本章还分别阐述了单根微热管、微热管阵列以及微平板热管的结构特点以及一些常见的沟道形状，并对微热管工作时可能发生的各种传热极限进行了简要的说明。-尽管在微热管领域已经有了一定的研究成果，但是该领域毕竟是一个全新领域，尚有很多问题有待解决。展望未来，主要的研究领域有：1）发展更精确、更可靠的微热管的机械加工技术以解决细微结构机械加工费用昂贵的问题，进而扩大微热管的使用领域；2）根据微热管的具体运用环境，设计出与热源表面结构相匹配的异截面形状热管；3）展开对热管结构可靠性研究以及热管在变热流密度、变热源位置条件下间歇性