基于Fluent石蜡相变材料模拟_第1页
基于Fluent石蜡相变材料模拟_第2页
基于Fluent石蜡相变材料模拟_第3页
基于Fluent石蜡相变材料模拟_第4页
基于Fluent石蜡相变材料模拟_第5页
已阅读5页,还剩14页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

一、问题背景:为了解决日益严重的能源短缺问题,如何更加充分地利用现有的化石能源,开发利用绿色能源成为世界各国关注的重要议题。随着现有化石能源的逐步开采利用,世界各国已经普遍认识到降低对传统能源(如煤炭、石油、天然气能源等)的依赖性,以及对绿色能源(如太阳能、风能、地热能等)实现充分开发利用的重要性,使用再生类能源并通过提高能源利用效率的方式成为应对能源枯竭现状的重要手段。蓄热技术就是这类能够提高能源利用效率的典型技术手段,蓄热技术通过将间歇性或者不稳定的热量通过蓄能介质暂时储存,在有使用要求时释放能量,解决能源利用高峰阶段造成的能源匹配不足的问题。经过多年应用发展,蓄热技术已经在太阳能、地热能、风能、工业废热、电网系统的“移峰填谷”等领域有了一定程度的应用,并表现出强劲的发展势头。二、蓄热技术蓄热技术一般通过利用蓄热介质的比热容、潜热等物理特性实现对采集能源多余热量的暂时储存,主要分为显热蓄热、潜热蓄热和化学热反应蓄热,其中前两种技术属于物理蓄热范畴。显热蓄热显热蓄热通过提升蓄热介质材料的温度进而提高物质内能的方式实现储热,储热能力取决于材料的比热容(提升物质单位温升所需要的外部能量)等物理参数。显热蓄热的突出弱点在于由于蓄能介质(例如水)在多余能量的储存过程中伴随着物质温度的上升,不能满足部分设备对于恒温放热的要求。同时,显热蓄热材料蓄热能力有限,储能密度较低,往往需要较大容积的容器提供储能保证,限制了显热蓄热技术的大规模推广应用。潜热蓄热潜热蓄热技术充分利用了相变介质在相态改变时会吸收或释放巨大的能量,并藉此实现对多余能量的储存和释放。这类利用相变过程实现能量吸收释放的材料被称为相变材料。相比显热蓄热材料,相变材料的相变潜热与蓄热材料的比热特性相比在储热能力方面有了极大的提升,同等质量的储热介质能够实现对更多能量的储存,降低了对储热设备容积、质量等方面的要求,降低了整体设备成本。同时,潜热蓄热材料在相变吸热放热过程中近似等温过程,方便了实际工程控制。尽管如此,由于相变蓄热材质在导热、传热特性方面的不足,以及某些无机盐相变材料存在的相分离和过冷现象以及某些固有的化学反应,会严重影响储热设备的使用寿命。化学热反应蓄热化学热反应在反应过程中会释放大量热量,与外部环境进行热交换,化学热反应蓄热即是利用这一热交互过程实现能量和吸收和释放。基本原理如下:利用正化学反应吸热,将能量转化为化学能储存,之后利用负反应放热,将通过正反应储存的能量通过热量的方式释放出来。化学热反应蓄热与相变潜热蓄热技术有相似的优势特点,但是由于化学热反应存在反应过程复杂、反应速度过快、反应过程复杂难于温度法模型适合于单一成分有固定相变温度的相变材料传热过程模拟,以温度为系统因变量,分别对不同状态下的相变材料进行控制能量方程的建立,并进一步进行温度求解。焓法模型相比前述温度法模型,焓法模型不仅适用于单一组分的相变传热模拟过程,对于模糊的边界移动等传热问题同样适用,该方法将焓作为系统因变量,在固、液以及混合模糊区域建立统一的能量守恒控制方程。图SEQ图\*ARABIC3固-液两相系统示意图如图3所示,为典型的相变材料固-液两相系统示意图,选择对应控制体V作为研究对象,经过相界面的移动,在时刻t,系统分为固体相VS和液体相Vl,将相间由于温度变化产生的循环流动速度定义为v,d式中,以温度和焓参数作为系统因变量,相互关系可以通过公式2实现转换:T公式2中,参数cs和cl分别表示相变材料固相比热和液相比热,参数hs和hl分别表示相变材料的液体相与固体相焓值,单位为J/ρρ式中,ρs和ρl分别表示固体相和液体相的密度参数,ks和kl表示对应相态的导热系数,Ts和Tl分别表示对应相态的温度,ρρ=ρs通过上述公式可以将活动区域以及相态界面区域的能量控制方程统一为一个共同的控制方程组,实现了相变传热模拟的数值计算。对于上述问题的求解方法,对于一维问题求解主要包括:Neumann法,Lightfoot积分法,Paterson法,对于多维问题求解主要包括有限差分法,有限元法和有限体积法。4.3基于Fluent的石蜡材料相变传热模拟4.3.1Fluent计算流体力学数值模拟软件介绍Fluent是一款较为常用的集成流体、热力学领域数值模拟技术的数值模拟软件,采用了多重网格技术,具有求解收敛稳定、收敛速度快的特点,软件提供了丰富的物理模型。依托ANSYS有限元仿真软件多场耦合平台,ANSYSMesh模块为Fluent提供了较好的前处理功能,可以高效集成多种CAD辅助建模软件,能够生成二维数值模拟所需要的三角形、四边形网格,生成三维数值模拟所需要的四面体、六面体网格,并能够方便实现对模型细节网格的处理,包括网格细化、节点耦合等处理方式。同时,ANSYS平台为Fluent软件提供了更人性化的结果后处理窗口CFD-Post。ANSYS平台下基于Fluent的流体数值模拟一般流程如图4所示。图SEQ图\*ARABIC4Fluent数值模拟基本流程4.3.2FluentSolidification/Melting模型介绍为了实现对相变材料传热条件下的相变传热过程进行数值模拟,Fluent提供了专业的融化模型(Solidification/Melting模型),该模型采用前述焓法模型作为系统控制方程内核,适用于带有模糊区域的相变传热过程模拟。在数值计算模拟过程中,由于相变材料由固态吸热变为液态,材料孔隙率由0逐渐增加为1,相反的,当材料由液相逐渐转变为固相状态时,材料孔隙率相应的由1逐渐转变为0,。同时,Solidification/Melting模型引入了液相率参数β,用温度表示为如下形式:β=式中,对于单一成分的相变材料模拟有Ts=Tl,对于混合成分的相变材料模型有Ts能量方程源项修正形式:S动量方程源项修正形式:SS其中:Amush表示模糊区域流体设计参数,一般取值为104~4.3.3基于Fluent的石蜡凝固和融化设置模拟研究本章将基于Fluent流体仿真软件的凝固-融化模型对圆管外侧石蜡蓄能材料的融化和凝固过程进行数值模拟,在此基础上对不同结构形式,不同热源(恒定进水温度以及太阳能辐射热源)条件下的液相分数以及对应相变材料的温度场数据进行了提取,以此为进一步研究石蜡材料的融化过程提供初始理论数据。如图5为本次设计所使用的石蜡相变材料融化模拟模型结构示意图(本次设计相变蓄热箱体结构总体容积为50L),本次模拟所采用的石蜡材料的相关物理特性如表1所示,同时为简化模拟,做出如下合理假设:1.石蜡相变材料分布均匀,切各物理特性属各项同性;2.石蜡液相区域属于牛顿不可压缩流体区域,不考虑可压缩性;3.石蜡液相区域热对流属于二维层流流动模型;4.模拟过程中仅仅考虑密度随温度的变化;5.管道厚度与热阻系数在模拟过程中忽略不计。图SEQ图\*ARABIC5圆管外石蜡相变材料相变过程模拟模型表SEQ表\*ARABIC1石蜡材料物理特性参数表名称密度(kg/m3)比热(相变潜热(kJ/kg)固体相温度点(K)液体相温度点(K)动力粘度(kg/(m.s))导热系数(W/(m.K))石蜡7602100170324.15330.150.003240.254.3.3.1Fluent石蜡凝固和融化参数设置1)进入ANSYS协同仿真平台,如图6所示,在工作窗口内建立Fluent流体仿真模块。图SEQ图\*ARABIC6基于ANSYS的Fluent流体仿真模块2)进入ANSYS仿真平台的建模模块(Geometry),如图7所示,在建模窗口内建立石蜡模型(容积50L,内部管道直径D1=100mm,内部管道直径D图SEQ图\*ARABIC7石蜡箱体模拟模型3)进入ANSYS仿真平台的网格划分模块,对石蜡模型进行网格划分,网格划分效果如图8所示(包括了58680个单元,63240个单元节点,网格平均质量0.93)。图SEQ图\*ARABIC8石蜡模型网格划分效果4)进入ANSYS仿真平台的Fluent仿真模块,设置求解器类型等相关设置,如图9,并依次对分析类型、模型材料属性、边界条件等相关参数进行设置。图SEQ图\*ARABIC9Fluent基本设置窗口选择求解类型为瞬态求解,并设置求解器为压力基求解器,设置重力方向为+Z方向(大小为9.8m/s图SEQ图\*ARABIC10求解设置激活Fluent融化-凝固模型,如图11所示;图SEQ图\*ARABIC11Fluent凝固-融化模型激活在材料定义窗口对石蜡材料的相关物理属性进行设置(如表1,包括石蜡材料的密度、比热容参数、导入系数、潜热参数、动力粘度以及热膨胀参数等物理参数),如图12所示。图SEQ图\*ARABIC12石蜡物理材料属性设置设置模型边界条件(热源温度、热流密度、热对流等相关参数),如图13所示(设置内壁温度为350K);图SEQ图\*ARABIC13模型边界条件设置选择求解器,并对求解存储子步等相关参数进行设置,最终求解计算。不考虑石蜡箱体结构与外部环境的散热情况,石蜡相变材料随时间变化温度变化及液相体积分数变化整理如下:1.如图14为入口水温80℃条件下,石蜡结构不同时刻(50min~500min)相变温度分布云图,对应液相体积分数分布云图如图15所示。从图中可以看出,石蜡材料温度从靠近热水管管壁位置开始得到提升,当加热时间为500min时,整体最大相变材料温度基本达到50℃左右,同时液相体积分数与温度变化基本保持一致。选取特点结果提取坐标点,提取温度变化曲线如图16所示。图SE

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论