纳米流体沸腾传热研究进展姚寿广;董招生【摘要】综述近年来纳米流体在池沸腾和流动沸腾传热领域的实验研究和数值模拟研究。

实验包括纳米流体沸腾传热的临界热流密度、沸腾换热系数以及换热机理方面的研究，并简要分析纳米流体强化或弱化沸腾传热的主要原因。

数值模拟主要介绍格子玻尔兹曼方法（Lattice Boltzmann Method，LBM）在纳米流体沸腾传热领域的最新研究进展，包括基于不同的LBM模型模拟气泡产生、成长到脱离壁面的过程，对气泡的脱离直径和频率进行分析。

最后展望纳米流体沸腾传热的发展方向。

%This paper presents the latest researches on the heat transfer of pool boiling and flow boiling in the re-cent years.The studies on the critical heat flux,heat transfer coefficient and mechanism are included in the ex-perimental researches on the heat transfer of boiling.Moreover,the main reasons are briefly analyzed why the nanofluids enhance or deteriorate the heat transfer in comparison to base fluid.Then,this paper presents a re-view of recent numerical simulation researches on the boiling heat transfer behavior of nanofluids based on LBM. With different LBMmodels,the nucleation and growth of bubbles and departure from the wall are simulated and investigated as well as the departure diameter and release period.Finally,several topics worthy of attention for future investigations are identified.【期刊名称】《江苏科技大学学报（自然科学版）》【年(卷),期】2017(031)001【总页数】6页(P49-54)【关键词】纳米流体;强化传热;沸腾换热;数值模拟【作者】姚寿广;董招生【作者单位】江苏科技大学能源与动力工程学院，镇江 212003;江苏科技大学能源与动力工程学院，镇江 212003【正文语种】中文【中图分类】TK172.4随着科学技术的飞速发展,工业界对研制出紧凑、轻量和高效的热交换设备提出了迫切的要求,发展高热导率和传热性能良好的换热工质已经成为换热设备中强化传热技术的研究热点之一[1]．纳米材料因具有小尺寸、大比表面积以及独特的光、磁、电、化学和机械特性,为强化传热提供了新的有效途径．自1995年美国Argonne国家实验室提出纳米流体作为强化传热介质以来,将纳米流体作为流动工质应用于各类紧凑型换热冷板或换热器，从而提高其传热性能,在电子设备散热及紧凑式换热器上已展示出良好的应用效果[2]．文中主要研究近期纳米流体在强化沸腾换热领域实验研究和数值模拟的进展．纳米流体的导热系数高于基液并能提高单相流体换热能力,而沸腾作为有效的强化换热手段,纳米流体沸腾传热的机理十分复杂,纳米流体沸腾是否强化换热,不同研究者得出了不同甚至完全相反的结果[3]．1.1 沸腾换热的实验研究文献[4]研究了Al2O3/H2O纳米流体的池沸腾传热特性,采用NiCr丝作为加热面,纳米颗粒的体积分数为0.01 %～0.1 %．结果发现，在每一个体积分数下,纳米流体的临界热流密度(critical heat flux,CHF)都比纯水有所提高,而且随着体积分数的增加而增加,最高可提高48%．文中研究得出纳米流体CHF的提高主要是由于池沸腾时纳米颗粒沉积在加热表面形成一层多孔的涂层,多孔涂层增加汽化核心所需要的微穴，并且纳米颗粒的沉积形成多孔介质,从而进一步强化换热,文中通过测量加热表面的粗糙度和扫描电子显微镜(scanning electron microscope,SEM)图像证实此观点．文献[5]对γ-Al2O3/CMC非牛顿纳米流体的池沸腾传热特性进行实验研究．结果表明：CMC溶液的核态沸腾的传热系数随着CMC浓度的升高而降低,但是在CMC较低浓度情况下,由于溶液的流变特性的变化,沸腾换热系数会有小幅提升．在CMC溶液中添加纳米粒子后,沸腾表面的温度降低,沸腾传热系数升高,而且随着纳米粒子质量浓度的升高,强化传热效果更加明显．文中认为影响纳米流体沸腾传热效果的因素中存在不同和对立的因素,如溶液的粘度,纳米颗粒与加热表面和气泡的冲突．对比与非牛顿流体基液,这些因素的不同组合作用导致非牛顿纳米流体有更好的强化传热表现．文献[6]研究低浓度(≤1 g/L)Al2O3/H2O纳米流体的池沸腾特性,结果发现：当纳米粒子的浓度低于0.025 g/L时,纳米流体的CHF较纯水都有提高,而且CHF随着纳米粒子浓度的升高而升高；当浓度为0.025 g/L时,纳米流体的CHF较纯水有80%的提高；当浓度大于0.025 g/L时，纳米流体的CHF大小不变,但是壁面的过热温度增大,说明池沸腾传热系数降低．文献[7]研究TiO2纳米流体的浓度、加热表面的材料和粗糙度对沸腾的影响．Cu、Al制成水平的圆盘作为加热面,表面粗糙度为0.2 μm和4 μm,纳米颗粒的体积浓度分别为0.000 05、0.000 1、0.000 5、0.005、0.01．结果发现,铜圆盘作为加热表面，浓度为0.000 1的纳米流体的传热系数较基液有所提高,且粗糙度为0.2μm时提高约为15%,粗糙度为4 μm时提高约为4%;当浓度超过0.000 1时,传热系数在两种粗糙度下均小于基液的传热系数,文中得出浓度小于0.000 1时,沸腾产生的向上微对流导致纳米粒子无法沉积进而强化传热,当浓度大于0.000 1时,纳米离子沉积在表面,增加了热阻,所以恶化传热过程．另外，同样的粗糙度,Al比Cu有更高的沸腾传热系数,而且无论加热表面哪种材质，粗糙度为4 μm时的传热系数要比0.2 μm时的大,这是由于粗糙的表面增加了沸腾的汽化核心,强化传热．文献[8]研究δ-Al2O3/R141b纳米流体的池内核态沸特性,将平均粒径为20 nm 的δ-Al2O3纳米粉添加到制冷工质R141b中制备体积分数为0.001%、0.01%、0.1 %的纳米流体,沸腾表面为2000#砂纸打磨光滑的铜表面．实验表明：纳米流体强化了沸腾传热特性,且强化倍数随着纳米流体浓度的增加而增大．体积分数0.1 %纳米流体的沸腾传热系数比基液增大了50.2 %,分析认为表面颗粒沉积是强化换热的主要因素．沉积颗粒的粒径可能比原沸腾表面的粗糙度大,覆盖在表面上增大了表面局部粗糙度,从而强化了沸腾传热系数．文献[9]研究Al2O3/H2O纳米流体的池沸腾传热性能和现象．纳米颗粒的平均直径为47 nm,沸腾表面为用2000#砂纸打磨过的水平矩形光滑的铜板,纳米流体的体积浓度分别为0.5%、1%、2%、4%,并用超声波振荡以确保稳定性,实验结果发现,纳米流体在自然对流和核态沸腾阶段的换热系数比纯水有所降低,而且随着纳米粒子浓度的升高,降幅增大．同时还对纳米流体的CHF进行测量,结果发现在水平加热面上的CHF提升32 %．分析认为纳米颗粒的沉积改变了表面的粗糙度,在沸腾表面形成一层纳米膜,导致沸腾表面的汽化核心减少,所以传热效果恶化．文献[10]对Al2O3/乙二醇纳米流体(质量分数为0.1%～0.3%)进行池沸腾换热特性实验研究,分析热流密度、沸腾表面的粗糙度、纳米流体的浓度以及污垢的热阻分别对纳米流体池沸腾换热系数的影响．实验结果表明：① 纳米流体在自然对流和核态沸腾阶段,沸腾传热系数随着热流密度的升高而增大,但是核态沸腾阶段传热系数增加更明显,文中认为主要的原因是核态沸腾产生的气泡将大部分热量从加热壁面传递给液体,而液体的湍动和局部扰动则是次要因素;② 纳米流体的传热系数随着纳米颗粒的浓度增加而减小,文中认为,随着浓度的增加,纳米流体颗粒在沸腾表面沉积,导致了沸腾表面的粗糙度降低,汽化核心的数量减少,还有气泡的导热系数低,导致热量传递恶化,所以沸腾传热系数降低．文献[11]实验研究了R-134a和R-134a/POE混合液分别与SiO2、CuO制成的纳米流体的流动沸腾的换热情况．结果发现，SiO2/R-134a纳米流体流动沸腾换热系数相比于基液有所减小,由于沸腾表面形成的纳米颗粒膜和纳米颗粒不稳定性导致纳米颗粒团聚是阻碍热交换的原因．但是在CuO/R-134a/POE纳米流体中,换热系数随着纳米颗粒体积分数的增加而增大,而且纳米颗粒对纳米流体流动沸腾压降的影响很小,几乎可以忽略,因此，纳米颗粒的浓度太低,不至于影响流体的动力粘度．文献[12-13]研究以制冷剂R113为基液的纳米流体水平管内流动沸腾的传热特性,纳米粒子CuO平均粒径为40 nm,质量分数为0～0.5 %．实验结果表明,纳米流体在光滑管内的流动沸腾传热系数比基液大,最大可以提高29.7 %,分析认为纳米粒子的扰动减薄了边界层的厚度和纳米粒子表面生成分子吸附膜是强化传热的主要原因．但是与此同时管内的压降比纯制冷剂增大,而且随着粒子浓度的增加而增加．文献[14]研究不同纳米粒子(Cu、Al、Al2O3、CuO)与制冷剂R141b制备的纳米流体(质量分数为0.1%、0.2%、0.3%)在水平热管内的流动沸腾换热特性,以及质量流量对沸腾的影响．结果表明：随着纳米颗粒浓度的增加,沸腾换热系数增大,而且其中0.3%的Cu/R141b纳米流体的强化换热效果较好,达到49%,这是由于Cu 的导热系数最大造成的．上述实验结果表明：纳米流体是否强化沸腾传热与纳米流体的性质(纳米颗粒的种类、大小、浓度以及稳定性等),基液的性质(粘度、表面张力等)和沸腾表面特性(粗糙度、润湿性能)等因素有关;这些因素的相互作用,不同配合也决定着纳米流体是否强化传热以及强化强度的大小．但是,目前的实验研究并未考虑到纳米颗粒不同的制备方法对其传热性能的影响,所以有必要结合纳米颗粒的制备、纳米颗粒的稳定性进行实验研究．1.2 纳米流体沸腾换热机理的研究由于影响因素复杂性,加上实验研究自身的局限,目前还未能清楚阐述纳米流体强化沸腾换热的机理,所以需要进一步展开大量的实验研究工作．文献[15]对Al2O3/H2O纳米流体在加热壁面的强制对流沸腾换热的CHF进行实验研究,结果表明，纳米流体的CHF相比于纯水提高,并且随着液流速度的提高而提高．研究认为纳米颗粒沉积导致壁面的湿润性发生变化,并通过SEM图像得以证明．文献[16]通过对不同类型纳米流体的过冷流动沸腾传热研究也得出相同的结论．文献[17]对纳米流体的核态沸腾中纳米颗粒的尺寸相对于沸腾表面的粗糙度进行综述并研究,重点探讨了纳米粒子的直径相对表面粗糙度对沸腾传热的影响．当加热壁面的粗糙度和纳米颗粒的尺寸差不多时,核态沸腾汽化核心减少,传热恶化;当壁面粗糙度和纳米颗粒的尺寸相差甚远的时候,由于汽化核心增加或者减小数量非常少,强化传热．文献[18]研究纳米颗粒悬浮液池内泡状沸腾机理,研究表明,纳米颗粒的加入对基液池内泡状沸腾换热过程有正反两方面的影响:一方面,纳米颗粒悬浮液有效导热系数较基液有了很大的提高,颗粒的加入减小了液体的表面张力,减少了沸腾过热度,从而强化液体内部换热过程;另一方面,颗粒的团聚以及液体粘度的变大易造成颗粒的沉积,颗粒的沉积和对凹坑的充填减小原加热表面的凹坑口径和数量,沸腾过热度变大,削弱沸腾换热．对于不同性质颗粒、不同浓度、不同热流密度的情况下,纳米流体沸腾换热的强化或恶化,是这些因素综合作用的结果．文献[19]从纳米流体沸腾核化的热力学角度探讨纳米流体更容易核化的主要原因是颗粒在界面处的聚集,有效降低了气液界面自由能．文献[20]从两个方面阐述纳米流体导热系数的强化机理:一是纳米颗粒改变液体的结构,使液体变成液-固两相悬浮液,影响了纳米流体内部的能量传递过程;二是由悬浮于基液中的纳米粒子随机布朗运动引起并促进了流体内部的微扰动,从而增强了纳米颗粒与基液之间的能量传递速率．综上实验研究,目前对于纳米流体沸腾换热机理的研究主要在以下几个方面:① 纳米颗粒沉积表面改变润湿性能;② 纳米颗粒的粒径与表面粗糙度的匹配;③ 纳米颗粒导热特性和纳米流体的粘度以及表面张力;④ 两相界面的自由能;⑤ 纳米粒子的布朗运动．但是目前还不清楚传热的机制,对于探究纳米流体强化传热的机理,还需要从纳米颗粒的稳定性、沉积表面的润湿性能和毛细结构、纳米流体的导热系数、表面张力和粘度、纳米流体的运动以及气泡动力学这几个方面进行研究．数值模拟作为研究手段之一,在深入研究多相流机理方面扮演着重要的角色．目前多相流换热的数值模拟方法主要分为界面追踪和界面捕捉两类,界面追踪是将每相分开考虑,该方法精度高,稳定性好,但是该方法难以处理激烈的、不规则变形，并且用该方法模拟3D情况比较复杂;界面捕捉的方法有VOF、level-set和phase-field,其中phase-field是基于扩散界面模型而提出,主要用于晶体生长和凝固．在众多的方法中,LBM由于具有处理、分析复杂流体的潜力而广受关注,一些学者已经采用这种方法进行了多相流的相关研究．文献[21]通过LBM的数值模拟研究微尺度池沸腾中热流密度、湿润性和重力对气泡脱离3种不同结构壁面的直径和释放频率的影响．研究发现气泡脱离的直径随着接触角的增加而线性增加,而且壁面的粗糙度会减缓增加的幅度．气泡直径可以大致看成是热流密度的线性函数,加热壁面的洞穴影响这个函数的系数大小．如果考虑壁面的导热和洞穴,发现气泡脱离频率是接触角的单调升函数,可以用抛物线近似表达．文献[22]研究了气液间的相变换热,通过在伪势模型中添加新的能源项来提高数值稳定性,并采用该模型模拟水平壁面上气泡的生成和分离过程．文献[23]基于相变传热模型,用LBM模拟光滑亲水表面添加疏水点后的整个沸腾过程,研究发现疏水点的添加增加了核化气泡,减少了成核时间,强化了沸腾传热,同时提高了沸腾的CHF,其原因是疏水点调节气泡在壁面的伸展．至于热流密度和温度的变化,发现生成气泡下面的疏水区域的热流密度很高,而且在三相交接处的热流密度最大但是温度最低．文献[24]采用多相流稳定模型模拟过热壁面的气泡生成和脱离过程,讨论了密度、重力、表面张力、接触角和壁面过热度对气泡脱离直径和频率的影响．模拟结果和实验关联式很好的吻合,结果表明接触角对气泡脱离直径几乎没影响,但是随着密度比的增加,气泡直径将减小,气泡的释放周期随着密度比的增加而延长．气泡生成的过程中,气泡的过热度分布不均,但气泡脱离之后,气泡温度下降到略低于饱和温度,气泡直径随着气泡上升逐渐增大．文献[25]对微通道流动沸腾换热的气泡进行研究,将多相模型与热格子LBM结合，提出一个混合模型来自动捕捉界面变化,模拟结果表明气泡的成长阻碍了流动的发展,导致流动阻力增加,直至气泡脱离壁面;气泡周围的相变潜热的消耗强化了传热过程,但是正在成核的气泡与脱离壁面的气泡结合会导致气泡过早脱离壁面,在某种程度上限制了传热效率．文献[26]提出一种不可融合、大密度比的两相流LBM模型,基于这种方法模拟了气泡流的复杂结构和速度场,模拟结果与理论解符合较好．文献[27]耦合多相模型和Inamuro 的LBM模型模拟了重力加速度对气泡的成长和脱离过程的影响．文献[28]用多相模型结合相变模型模拟了重力加速度、接触角等对气泡生长过程的影响．文献[29]提出将Inamuro模型和文献[30]所提模型耦合来模拟气泡的形成、上升和脱离过程,研究了接触角、热流密度和加热壁面的长度对气泡离开壁面时间和直径的影响,发现接触角的增大、重力的减小以及热流密度的增加可以增加气泡脱离壁面的直径,伸展壁面的长度也能够明显增加气泡脱离时的直径,并指出当小气泡上升时受到很大的重力,则气泡将消失．文献[31]将LBM和有限体积法(finite volume method,FVM)结合,认为该模型存在三大优势:① 从宏观参数到LBM分布函数,形成了统一的重构算子;② 该模型可以解决LBM中密度变化和FVM中不连续的问题;③ LBM和FVM之间时间的耦合可以解决模拟不稳定性的问题．文中用该方法做了3个案例,分别是高斯脉冲的对流扩散、流体在圆柱体表面的流动和方腔内不稳定的自然对流,证明了该方法的灵活性与可行性．目前研究纳米流体强化换热的LBM模拟多集中于单相流体换热,如文献[32]采用LBM中双分布函数模型,模拟矩形腔内Al2O3/H2O纳米流体的Rayleigh-Benard 对流,但是进一步考虑纳米流体沸腾相变换热的LBM模拟研究还很少,仅有文献[33-35]基于伪势模型,通过加入纳米颗粒的布朗力和范德瓦尔力,构建适合模拟纳米流体等温相分离的LBM模型，并以此纳米流体LBM模型为基础耦合相变方程,建立一种能够描述纳米流体气液相变的多组分多相LBM模型;在大密度比模型基础上,耦合Inamuro的传热模型,根据纳米流体物性参数修改平衡态方程,建立描述纳米流体沸腾换热的单组分多相格子Boltzmann模型，模拟纳米流体沸腾过程中气泡成长及脱离过程．从已有研究工作看,基于LBM模拟纳米流体沸腾相变换热的基础研究较少,当前的研究大多把纳米流体看成是一种单组分单相的均质流体,这与实际的纳米流体有所出入,所以还应从以下几个方面进行大量的研究工作:① 将纳米流体做多组分处理,把纳米颗粒看作一个组分,基液看作另一种组分,考虑纳米颗粒所受到的布朗力、范德瓦尔力等;② 由于沸腾换热的机理复杂性,还应考虑相变过程,所以纳米流体的沸腾换热是多组分多相的,建立多组分多相LBM模型对研究其沸腾传热机理必不可少;③ 结合气泡动力学、沸腾表面的润湿性能以及纳米流体的粘度和表面张力等因素,模拟出纳米流体沸腾换热过程,并探讨各因素对换热的影响．从目前的实验研究和数值模拟来看,纳米流体能否强化沸腾传热还不能一概而论,但是基于单相纳米流体强化传热的共识,纳米流体强化沸腾传热的潜力得到大多学者的认同．虽然对纳米流体强化沸腾传热机理的理解是不充分的,但对于纳米流体沸腾CHF的研究结果却一致是增大的,目前认为主要原因是纳米颗粒沉积改变沸腾表面的微观形态结构,进而增强壁面的湿润性、粗糙度和虹吸作用．而对于提高纳米流体沸腾传热系数的研究结果却是对立的,由于纳米颗粒的种类、尺寸、浓度、沸腾表面的特性、热物性的改变等多种因素共同影响纳米流体的换热性能．强化沸腾传热的主要原因有:添加纳米颗粒增加了汽化核心;纳米颗粒的扰动降低了边界层厚度,减小了热阻;弱化传热的主要集中在纳米颗粒沉积表面形成纳米膜,增大热阻,减小沸腾的汽化核心等．由于纳米流体的复杂性,至今没有彻底清楚纳米流体强化传热的机理．为了理解纳米流体沸腾传热的机理,还可从以下几个方面进行实验或数值模拟的研究:① 制备能够长时稳定的纳米流体;② 结合沸腾表面特性(润湿度、粗糙度)和纳米颗粒(尺寸、浓度)相互作用的研究;③ 结合气泡动力学和流体的粘度、表面张力研究对沸腾的影响;④ 建立起可以预测纳米流体沸腾换热系数、CHF、压降、热阻的实验关联式,为工程应用提供参考依据．【相关文献】[1] 杨波,王姣,刘军. 碳纳米流体强化传热研究[J]. 强激光与粒子束,2014,26(5):1-3. YANGBo,WANG Jiao,LIU Jun. Heat transfer enhancement of carbon nanofluids[J]. High Power Laser and Particle Beams,2014,26(5):1-3.(in Chinese)[2] 张正国,燕志鹏,方晓明,等. 纳米技术在强化传热中应用的研究进展[J]. 化工进展,2011,309(1):34-39. ZHANG Zhengguo,YAN Zhipeng,FANG Xiaoming,et al. Research development of applications of nanotechnology in heat transfer enhancement[J]. Chemical Industry and Engineering Progress,2011,309(1):34-39. (in Chinese)[3] TAYLOR R A,PHELAN P E. Pool boiling of nanofluids:comprehensive review of existing data and limited new data [J]. International Journal of Heat and Mass Transfer,2009,52(23-24):5339-5347.[4] HEGDE R N, RAOB S S,REDDY R P. Flow visualization and study of critical heat flux enhancement in pool boiling with Al2O3-water nanofluids[J]. ThermalScience,2012,16(2):445-453.[5] SOLTANI S, ETEMAD S G,THIBAULT J. Pool boiling heat transfer of non-Newtonian Nanofluids [J]. International Communications in Heat and Mass Transfer,2010,37(1):29-33.[6] KWARK S M, KUMAR R, MORENO G,et al. Pool boiling characteristics of low concentration nanofluids[J]. International Journal of Heat and Mass Transfer,2010,53(5-6):972-981.[7] SURIYAWONG A, WONGWISES S. Nucleate pool boiling heat transfer characteristics of TiO2-water nanofluids at very low concentrations [J]. Experimental Thermal and Fluid Science,2010,34(8):992-999.[8] 唐潇,刁彦华,赵耀华,等. δ-Al2O3-R141b纳米流体的池内核态沸腾传热特性[J]. 化工学报,2012,63(1):64-70. TANG Xiao,DIAO Yanhua,ZHAO Yaohua,et al. Nucleate pool boiling heat transfer of δ-Al2O3-R141b nanofluid on horizontal plate[J]. CIESCJournal,2012,63(1):64-70. (in Chinese)[9] BANG I C,CHANG S H. Boiling heat transfer performance and phenomena of Al2O3-water nano-fluids from a plain surface in a pool[J]. International Journal of Heat and Mass Transfer,2005,48:2407-2419.[10] SARAFRAZ M M, HORMOZI F. Nucleate pool boiling heat transfer characteristics of dilute Al2O3-ethyleneglycol nanofluids[J]. International Communications in Heat and Mass Transfer，2014,58:96-104.[11] HENDERSON K,PARK Y G,LIU L P,et al. Flow-boiling heat transfer of R-134a-based nanofluids in a horizontal tube [J]. International Journal of Heat and MassTransfer,2010,53(5-6):944-951.[12] PENG Hao,DING Guoliang,JIANG Weiting,et al. Heat transfer characteristics of refrigerant-based nanofluid flow boiling inside a horizontal smooth tube [J]. International Journal of Refrigeration,2009,32(6):1259-1270.[13] PENG Hao,DING Guoliang,JIANG Weiting,et al. Measurement and correlation of frictional pressure drop of refrigerant-based nanofluid flow boiling inside a horizontal smooth tube [J]. International Journal of Refrigeration,2009,32(7):1756-1764.[14] SUN Bin,YANG Di. Flow boiling heat transfer characteristics of nano-refrigerants in a horizontal tube[J]. International Journal of Refrigeration,2014,38:206-214.[15] AHN H S,KIM H,JO H J,et al. Experimental study of critical heat flux enhancement during forced convective flow boiling of nanofluid on a short heated surface [J]. International Journal of Multiphase Flow,2010,36(5):375-384.[16] KIM S J,MCKRELL T,BUONGIORNO J,et al. Experimental study of flow critical heat flux in alumina-water,zinc-oxide-water,and diamond-water nanofluids [J]. Journal of Heat Transfer,2009,131(4):1-7.[17] DAS S K,PRAKASH N G,BABY A K. Survey on nucleate pool boiling of nanofluids:the effect of particle size relative to roughness [J].Journal of NanoparticleResearch,2008,10(7):1099-1108.[18] 帅美琴, 施明恒, 李强. 纳米颗粒悬浮液池内泡状沸腾机理[J]. 东南大学学报(自然科学版), 2006, 36(5):785-789. SHUAI Meiqin,SHI Mingheng,LI Qiang. Pool boiling heat transferm echanism of nano-particle suspensions[J]. Journal of Southeast University(Natural Science Edition),2006,36(5):785-789. (in Chinese)[19] 李启明, 王补宣, 段远源. 纳米流体沸腾核化的热力学探讨[J]. 化工学报,2011,62(11): 2987-2992. LI Qiming,WANG Buxuan,DUAN Yuanyuan. Thermodynamics of nanofluid boiling nucleation[J]. CIESC Journal,2011,62(11):2987-2992. (in Chinese)[20] 宣益民. 纳米流体能量传递理论与应用[J]. 中国科学: 技术科学,2014,44(3):269-278. XUAN Yimin. An overview on nanofluids and applications[J]. Scientia SinicaTechnologica,2014,44(3):269-278. (in Chinese)[21] MRKUS A,HZI G. On pool boiling at microscale level: the effect of a cavity and heat conduction in the heated wall[J]. Nuclear Engineering and Design,2012,248:238-247. [22] GONG Shuai,CHENG Ping. A lattice Boltzmann method for simulation of liquid-vapor phase-change heat transfer[J]. International Journal of Heat and MassTransfer,2012,55:4923-4927.[23] GONG Shuai,CHENG Ping. Numerical simulation of pool boiling heat transfer on smooth surfaces with mixed wettability by lattice Boltzmann method[J]. International Journal of Heat and Mass Transfer,2015,80:206-216.[24] BEGMOHAMMADI A, FARHADZADCH M, RAHIMIAN M H. Simulation of pool boiling and periodic bubble release at high density ratio using lattice Boltzmann method[J]. International Communications in Heat and Mass Transfer,2015,61:78-87.[25] DONG Zhiqiang,XU Jinliang,JIANG Fangming,et al. Numerical study of vapor bubble effect onflow and heat transfer in microchannel[J]. International Journal of Thermal Sciences,2012,54:22-32.[26] INAMURO T,OGATA T,TAJIMA S,et al. A lattice Boltzmann method for incompressible two-phase flows with large density differences[J]. Journal of ComputationalPhysics,2004,198:628-644.[27] SUN Tao,LI Weizhong,YANG Shuai. Numerical simulation of bubble growth and departure during flow boiling period by lattice Boltzmann method[J]. International Journal of Heat and Fluid Flow,2013,44:120-129.[28] RYU Seungyeob,KO Sungho. Direct numerical simulation of nucleate pool boiling using a two-dimensional lattice Boltzmann method[J]. Nuclear Engineering and Design,2012,248: 248-262.[29] SATTARI E,DELAVAR M A,FATTAHI E, et al. Numerical investigation the effects of working parameters on nucleate pool boiling[J]. International Communications in Heat and Mass Transfer,2014,59:106-113.[30] TANAKA Y,WASHIO Y,YOSHINO M,et al. Numerical simulation of dynamic behavior of droplet on solid surface by the two-phase Lattice Boltzmann method[J]. Computers & Fluids，2011,40 (1):68-78.[31] TONG Zixiang,HE Yaling. A unified coupling scheme between lattice Boltzmann method and finite volume method for unsteady fluid flow and heat transfer[J]. International Journal of Heat and Mass Transfer,2015,80: 812-824.[32] 郭亚丽,徐鹤函,沈胜强,等. 利用格子Boltzmann方法模拟矩形腔内纳米流体Raleigh-Benard 对流[J].物理学报,2014,62(14):1-5. GUO Yali,XU Hehan,SHEN Shengqiang,et al. Nanofluid Raleigh-Benard convection in rectangular cavity: simulation with lattice Boltzmann method[J]. Acta Physica Sinica,2014,62(14):1-5. (in Chinese)[33] YAO Shouguang,ZHOU Changjiang,NIE Yuhong. Numerical simulation of the nanofluid phase separation by lattice Boltzmann method[J]. Advanced Materials Research,2014:619-622.[34] 姚寿广,贾新旺,胡安杰,等. 基于LBM的管内纳米流体相变过程的分析[C]∥中国工程热物理学会学术会议论文集,西安,2014.[35] YAO Shouguang,JIA Xinwang,HUANG Tao,et al. Numerical simulation of bubble motion in boiling nanofluids based on Lattice Boltzmann method[J].International Journal of Heat and Technology,2015,33(1):71-76.。