基于LBM的多孔骨架热物性对固液相变的影响研究.docx

该【基于LBM的多孔骨架热物性对固液相变的影响研究 】是由【科技星球】上传分享，文档一共【27】页，该文档可以免费在线阅读，需要了解更多关于【基于LBM的多孔骨架热物性对固液相变的影响研究 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。基于LBM的多孔骨架热物性对固液相变的影响研究??宋林泉,陈宝明,2,3,*,郜凯凯(,山东济南250101;,山东济南250101;,山东济南250101)基于LBM的多孔骨架热物性对固液相变的影响研究宋林泉1,陈宝明1,2,3,*,郜凯凯1(,山东济南250101;,山东济南250101;,山东济南250101)对不同物性骨架对固液相变过程的影响研究可为中低温相变储能技术的应用和发展奠定理论基础。文章基于格子玻尔兹曼方法(LBM),采用两区域焓—多孔介质模型研究了方腔内无填充多孔介质骨架固液相变过程,从孔隙尺度分析了相变过程的流动和传热机理,探讨了方腔内填充不同导热系数的骨架对于相变过程的影响。结果表明:在无填充多孔介质骨架方腔内固液相变过程中传热方式由热传导逐渐向自然对流换热转变,形成向右倾斜的糊状区;它的存在导致相变材料不能完全融化,且在方腔的左侧壁面处存在上窄下宽的固相相变材料;在填充多孔介质骨架方腔内,融化的初始阶段,高导热系数多孔骨架的相变材料融化速率较大,对相变换热起到了明显的促进作用,而当相变过程发展至准稳态阶段,受到右壁面处的低温影响和糊状区的综合作用,相变过程受到明显的抑制,且骨架的导热系数越大,其融化率越低。固液相变;格子玻尔兹曼方法;糊状区;多孔骨架;两区域焓—多孔介质模型Abstract:Changesofheattransferandflowfilledwithdifferentthermalconductivityskeletononthesolid--zoneenthalpy-porousmodeltostudythesolid-liquidphasetransitionprocessbasedontheLatticeBoltamannMethod(LBM).Onthisbasis,thesolid-,,,whenthephasetransitiondevelopstothequasisteadystate,:solid-liquidphasetransition;LatticeBoltamannMethod;mushyregion;porousskeleton;two-zoneenthalpy-porousmodel0引言多孔介质内固液相变传热流动研究广泛存在于各行各业中,例如建筑外围护结构内添加相变材料[1-4]可以有效地减弱室外不利环境对室内温湿度的影响。在相变过程中,相变材料特性(融化温度、相变半径)影响方腔内的流动、传热以及相变界面的移动。对于单质(如铝、铁等金属物质),其固液相变过程中有显著的固液交界面;而相变物质(如石蜡、羧酸、多元醇等有机物),其融化温度往往发生在一个区间,此时就会存在一个固相和液相共存的区域(糊状区)。由于两相界面随融化过程的进行而移动,在界面处的能量守恒条件是非线性的,同时伴随着热质交换的进行,因此相变过程是一个非稳态、非线性的问题。Kazmierczak首次基于Darcy模型进行了多孔介质中相变材料融化的自然对流现象的研究[5]。Jany等从孔隙尺度分析了填充多孔介质方腔内相变材料的融化过程[6],并总结出了自然对流换热强度对壁面平均Nu数的影响。之后,Bejan采用相同的方法研究了Ste数的影响[7]。近年来,利用金属泡沫和骨架的高导热性对蓄能换热器的相变传热效率强化进行逐渐重视。Mancin对具有不同相变温度的三种石蜡在同孔隙率、孔密度的铜泡沫中的固液相变过程进行了对比实验[8]。何叶从等基于焓法分析了相变墙板传热性能的影响因素[9]。崔娜等采用焓法分析了石膏基—石蜡相变储能构件在被动式建筑中的应用效果[10]。固液相变过程中存在着界面能量守恒的非线性特征、变化过程的非稳定性以及孔介质结构的复杂性,现阶段多孔介质内固液相变过程的研究主要采用的是数值模拟方法。而格子玻尔兹曼方法LBM(LatticeBoltamannMethod)pntationalFluidDynamics)方法相比,具有处理复杂边界、容易编程、后处理简单和并行性高等优点。并且在填充多孔介质方腔内的传热流动以及多相流的数值模拟等方面取得了许多成功的应用。许多学者将LBM应用在固液相变传热过程研究中。Wen等在前人对纯导热LBM模型研究的基础上,建立了固液相变过程中的的纯导热LBM模型[11]。夏莉等采用VOF(VolumeofFluid)模型研究了圆柱形容器中具有自由液面的石蜡相变过程,分析了多孔介质模型中形态系数对石蜡熔化过程的影响[12]。Gao等考虑局部非热平衡效应,使用REV(RepresentativeElementaryVolume)和孔隙尺度上的LBM模型对泡沫金属中的固液相变过程进行了模拟研究[13]。虽然目前大量对无填充多孔介质的固液相变过程研究表明,相变过程中固液共存的区域(糊状区)会显著影响固液相变的过程[14-15],但对于填充多孔骨架方腔的储能材料固液相界面的传递机理缺乏进一步的研究,尤其是相界面糊状区的传热和流动的模型缺少深入的探讨。文章拟以前期的研究为基础,根据非平衡热力学和复杂介质中的相变机理,建立糊状区的两区域模型,揭示无填充多孔介质内固液相变过程的流动和热质交换的规律,以及填充不同导热系数的骨架对固液相变过程的传热和流动的变化规律,为实际工程应用提供必要的理论依据,为中低温相变储能技术的研究和应用奠定理论基础。1多孔骨架的固液相变两区域模型构建多孔骨架模型如图1所示,方腔的边长为H,腔体内填充摆放边长为1/15H的正方形多孔介质骨架共5×5个,。其中黑色部分代表固体骨架,而白色部分为填充均匀的相变材料。方腔内部填充相变材料,,。相变材料在融化过程中当温度达到θms=θm-θR时开始融化;而完成融化时,对应的温度达到θml=θm+θR,方腔上下壁面为绝热壁面,整个方腔内,初始温度θc=0。融化的初始时刻,左侧的高温壁面温度θh为1,高于相变材料的融化温度。基于糊状区的焓—多孔介质模型,文章提出两区域模型,即对糊状区中高液相率区(rfl>rtr)采用液—固两相流,将其看作不同组分均匀混合的单相介质,将原来适用单质流体的宏观输运方程直接用于液—固两相流相流区域,而涉及到的相关物性参数则使用表征参数;而低含液率区域(rfl>rtr)采用适用范围更广的Brinkmann-Forchheimer-Darcy多孔介质渗流模型,其中渗透率、形态系数由液相率计算得到,从而建立更加准确的数学模型。相变材料的融化过程可分为3个区域,分别为液相区(红色)、固相区(蓝色)以及固液共存区域——糊状区(介于红色和蓝色之间的区域),如图2所示。=。为了便于建立数学模型,作出了5项假设:(1)相变过程中液相区流体(含液率rfl=1)视为不可压缩牛顿流体;(2)方腔内液相相变材料的流动为层流;(3)相变材料与填充腔体内的多孔介质骨架物性参数均视为常数;(4)考虑液相区内的自然对流,流体中的粘性耗散忽略不计,除密度外其它物性为常数;(5)相变过程中相变材料的体积变化忽略不计。在假设的基础上,Brinkmann-Forchheimer-Darcy模型可表示为广义N—S方程的形式,在一定条件下广义N—S方程可以简化为标准的N—S方程,而基于孔隙尺度的流体流动仍遵循标准N—S方程,因此REV尺度的糊状区模型与孔隙尺度的固液相变过程可统一为一区域方程。基于焓—多孔介质模型,糊状区内多孔介质区域的孔隙率ε由该区域的液相率rfl获得,因此,孔隙尺度下含糊状区的固液相变的控制方程分别由式(1)~(3)分别表示为式中:fl为相变材料液相;fs为相变材料固相;s为多孔骨架;v为有效运动粘性系数,m2/s;u为渗流速度矢量,m/s;p为表观压力,Pa;ε为多孔介质孔隙率,ε=0为固相区,ε=1为液相区,0<ε<1为糊状区;k为导热系数,W/(m·K);ktotal为有效导热系数,W/(m·K);vfl表示液相流体的动力粘性系数,m2/s;β表示液相流体的热膨胀系数,1/K;F为外力源项由式(4)表示为式中:K为多孔介质的渗透率;Fε为形状因子。在糊状区的低含液率区域(rfl<rtr),对应的渗透率K、有效导热系数ktotal分别由式(5)、(6)表示为在液相区rfl=1(即ε=1),控制方程演化为无多孔介质自然对流条件下固液相变的流动和传热控制方程;在固相区rfl=0(即ε=0),进而速度u=0,控制方程则演化为导热的控制方程;在糊状区的低液相率区域0<rfl<rtr,即为多孔介质区域,该区域内流动、传热满足上述控制方程。在高液相率1>rfl>rtr,液固两相流的表征导热系数、表征运动粘度由式(7)、(8)[16-17]表示为式中:φ为相变材料固相体积分数,即φ=1-rfl。多孔介质骨架传热为纯导热过程能量方程由式(9)表示为式中:ks表示多孔骨架的导热系数;Ts表示多孔骨架的温度。采用焓法[18]求解液相率,相变材料在某温度时的焓值由式(10)表示为相变材料的焓值En和温度Tf之间的关系由式(11)表示为式中:Ens、Enl分别为相变开始时相应温度Tfs的焓值和相变完成时对应温度Tfl相应的焓值;TR为相变半径,TR=(Tfl-Tfs)/2。根据相变材料的液相率和焓值之间的关系,计算得到式(12)为由式(10)、(11),结合式(12)可知,多孔介质骨架传热的纯导热方程(9)中的温度场与含液率是相互耦合的,可以通过数值迭代求解。为减少控制方程中的变量,探讨相变过程中糊状区对相变过程影响的机理,引入了无量纲参数