基于Fluent石蜡相变材料模拟.docx

基于Fluent石蜡相变材料模拟.docx

《基于Fluent石蜡相变材料模拟.docx》由会员分享，可在线阅读，更多相关《基于Fluent石蜡相变材料模拟.docx（16页珍藏版）》请在冰豆网上搜索。

基于Fluent石蜡相变材料模拟

基于Fluent石蜡相变材料模拟

一、问题背景：

为了解决日益严重的能源短缺问题，如何更加充分地利用现有的化石能源，开发利用绿色能源成为世界各国关注的重要议题。

随着现有化石能源的逐步开采利用，世界各国已经普遍认识到降低对传统能源（如煤炭、石油、天然气能源等）的依赖性，以及对绿色能源（如太阳能、风能、地热能等）实现充分开发利用的重要性，使用再生类能源并通过提高能源利用效率的方式成为应对能源枯竭现状的重要手段。

蓄热技术就是这类能够提高能源利用效率的典型技术手段，蓄热技术通过将间歇性或者不稳定的热量通过蓄能介质暂时储存，在有使用要求时释放能量，解决能源利用高峰阶段造成的能源匹配不足的问题。

经过多年应用发展，蓄热技术已经在太阳能、地热能、风能、工业废热、电网系统的“移峰填谷”等领域有了一定程度的应用，并表现出强劲的发展势头。

二、蓄热技术

蓄热技术一般通过利用蓄热介质的比热容、潜热等物理特性实现对采集能源多余热量的暂时储存，主要分为显热蓄热、潜热蓄热和化学热反应蓄热，其中前两种技术属于物理蓄热范畴。

●显热蓄热

显热蓄热通过提升蓄热介质材料的温度进而提高物质内能的方式实现储热，储热能力取决于材料的比热容（提升物质单位温升所需要的外部能量）等物理参数。

显热蓄热的突出弱点在于由于蓄能介质（例如水）在多余能量的储存过程中伴随着物质温度的上升，不能满足部分设备对于恒温放热的要求。

同时，显热蓄热材料蓄热能力有限，储能密度较低，往往需要较大容积的容器提供储能保证，限制了显热蓄热技术的大规模推广应用。

●潜热蓄热

潜热蓄热技术充分利用了相变介质在相态改变时会吸收或释放巨大的能量，并藉此实现对多余能量的储存和释放。

这类利用相变过程实现能量吸收释放的材料被称为相变材料。

相比显热蓄热材料，相变材料的相变潜热与蓄热材料的比热特性相比在储热能力方面有了极大的提升，同等质量的储热介质能够实现对更多能量的储存，降低了对储热设备容积、质量等方面的要求，降低了整体设备成本。

同时，潜热蓄热材料在相变吸热放热过程中近似等温过程，方便了实际工程控制。

尽管如此，由于相变蓄热材质在导热、传热特性方面的不足，以及某些无机盐相变材料存在的相分离和过冷现象以及某些固有的化学反应，会严重影响储热设备的使用寿命。

●化学热反应蓄热

化学热反应在反应过程中会释放大量热量，与外部环境进行热交换，化学热反应蓄热即是利用这一热交互过程实现能量和吸收和释放。

基本原理如下：

利用正化学反应吸热，将能量转化为化学能储存，之后利用负反应放热，将通过正反应储存的能量通过热量的方式释放出来。

化学热反应蓄热与相变潜热蓄热技术有相似的优势特点，但是由于化学热反应存在反应过程复杂、反应速度过快、反应过程复杂难于控制等问题，对于热反应发生装置要求较高，目前仅仅在一些特殊要求应用场合得到了一定程度的应用。

综上所述，在热能综合利用和蓄能技术对比中，采用相变材质的潜热蓄热技术具有更为突出的应用前景，如图1所示为一种典型的通过相变材料实现蓄热的蓄热系统构成（相变蓄能热水器系统）。

该系统首先通过吸收外部太阳辐射热量对水进行加热，并导入热水箱结构进行储存，对于多余的能量，系统通过热水将热能再次转化为相变材料的相变能量，实现对多余能量的储存。

在夜间缺少太阳能的时候，可以通过相变蓄热模块对冷水进行初步加热，并将初步加热的温水导入热水箱通过传统的辅助加热设备加热到固定的使用温度。

上述系统能够在连续多天没有阳光或者日照不足的情况下，实现对所储存能量的更高效的利用，降低对传统能源的依赖程度。

图1相变材料实现蓄热的蓄热系统构成

三、相变材料介绍

根据蓄热材料的组成分为无机类、有机类和复合相变材料三个大的类别。

通过相变材料的相变过程温度变化范围分为高温相变材料（相变温度≥250℃），中低温相变材料（相变温度：

100℃～250℃），低温相变材料（相变温度≤100℃）。

如图2所示，对相关相变材料分类及类别内的物质种类进行了简要划分。

图2相变材料划分及典型代表

四、基于Fluent的石蜡相变材料蓄热模拟

本次毕业设计采用ANSYSFluent计算流体力学仿真软件对考虑对流状态下的石蜡相变材料融化和凝固过程进行了数值模拟，研究石蜡材料在热交互过程中的储热机理和传热效果，为基于潜热蓄热激励的储能设备的设计提供参考数据。

4.1固-液相变传热相关机理

相变材料的凝固和融化过程实际上完全是一个物理变化过程，在热力学中往往利用潜热的概念（KJ/kg）对单位质量的相变物质相态转化所需要的能量进行定义，并称相变材料在相态转化过程中的潜热为溶解热（或者凝固热）。

相变材料相态转变过程主要分为三个阶段：

1.相变材料的显热存储阶段：

相变材料的温度低于外部热源温度，材料与热源发生热传导实现材料之间的能量交互，相变材料温度逐步上升并达到物质熔点温度。

2.相变材料的潜热存储阶段：

相变材料达到材料熔点后，依旧持续从外部热源吸收热量并转化为自身能量，但在存储阶段相变材料的温度不变，吸收的热量转化为融化热能量，并促使相变材料由固态转变为液态。

3.相变材料的第二显热存储阶段：

当全部相变材料转换为液态形式，持续从外部热源吸收的热量将促使液态相变材料温度进一步升高。

相变传热过程是一个复杂的固-液两相转化过程，相变过程中固-液交接面的移动依赖于外部边界条件，界面移动的速度取决于相变材料结构、属性以及外部热力学边界条件，涉及了固-液两相界面的移动，在数学上属于高度非线性的复杂问题，无法单独对系统方程进行求解，需要采用单独的求解器对系统进行分别求解。

4.2相变传热的数学模型和求解机理

由于相变材料的相变过程中，对于单一的相变材料，由于有确定的相变温度，存在不同状态条件下的相变材料共存的现象。

不同状态的相变材料由于具有不同的物理属性（导热性质和比热容），并且在固态相变材料中仅仅具有热传导方式的热量转化而在流体状态的相变材料中还包括了热对流的传热方式。

目前，对于固-液相变材料传热模拟的数值求解思路主要包括以下两种：

1.温度法模型

温度法模型适合于单一成分有固定相变温度的相变材料传热过程模拟，以温度为系统因变量，分别对不同状态下的相变材料进行控制能量方程的建立，并进一步进行温度求解。

2.焓法模型

相比前述温度法模型，焓法模型不仅适用于单一组分的相变传热模拟过程，对于模糊的边界移动等传热问题同样适用，该方法将焓作为系统因变量，在固、液以及混合模糊区域建立统一的能量守恒控制方程。

图3固-液两相系统示意图

如图3所示，为典型的相变材料固-液两相系统示意图，选择对应控制体V作为研究对象，经过相界面的移动，在时刻t，系统分为固体相

和液体相

，将相间由于温度变化产生的循环流动速度定义为v,控制体内材料的变化以及相界面移动取决于交界面S处的传热边界条件，初始材料温度以及材料的热力学物理属性。

对于，焓法模型，应用积分形式（如公式1）：

式中，以温度和焓参数作为系统因变量，相互关系可以通过公式2实现转换：

公式2中，参数

和

分别表示相变材料固相比热和液相比热，参数

和

分别表示相变材料的液体相与固体相焓值，单位为

。

将上述方程带入分别表示固体相和液体相的能量控制方程（公式3和公式4）：

式中，

和

分别表示固体相和液体相的密度参数，

和

表示对应相态的导热系数，

和

分别表示对应相态的温度，

表示梯度算子。

想要求解温度控制方程，需要首先对流体场的速度方程及相关系数进行求解。

将公式1,2带入公式3,4后，可以得到如下方程：

通过上述公式可以将活动区域以及相态界面区域的能量控制方程统一为一个共同的控制方程组，实现了相变传热模拟的数值计算。

对于上述问题的求解方法，对于一维问题求解主要包括：

Neumann法，Lightfoot积分法，Paterson法，对于多维问题求解主要包括有限差分法，有限元法和有限体积法。

4.3基于Fluent的石蜡材料相变传热模拟

4.3.1Fluent计算流体力学数值模拟软件介绍

Fluent是一款较为常用的集成流体、热力学领域数值模拟技术的数值模拟软件，采用了多重网格技术，具有求解收敛稳定、收敛速度快的特点，软件提供了丰富的物理模型。

依托ANSYS有限元仿真软件多场耦合平台，ANSYSMesh模块为Fluent提供了较好的前处理功能，可以高效集成多种CAD辅助建模软件，能够生成二维数值模拟所需要的三角形、四边形网格，生成三维数值模拟所需要的四面体、六面体网格，并能够方便实现对模型细节网格的处理，包括网格细化、节点耦合等处理方式。

同时，ANSYS平台为Fluent软件提供了更人性化的结果后处理窗口CFD-Post。

ANSYS平台下基于Fluent的流体数值模拟一般流程如图4所示。

图4Fluent数值模拟基本流程

4.3.2FluentSolidification/Melting模型介绍

为了实现对相变材料传热条件下的相变传热过程进行数值模拟，Fluent提供了专业的融化模型（Solidification/Melting模型），该模型采用前述焓法模型作为系统控制方程内核，适用于带有模糊区域的相变传热过程模拟。

在数值计算模拟过程中，由于相变材料由固态吸热变为液态，材料孔隙率由0逐渐增加为1，相反的，当材料由液相逐渐转变为固相状态时，材料孔隙率相应的由1逐渐转变为0,。

同时，Solidification/Melting模型引入了液相率参数

，用温度表示为如下形式：

式中，对于单一成分的相变材料模拟有

，对于混合成分的相变材料模型有

切相变率参数满足

。

同时在焓变计算过程中，Fluent对能量方程和动量方程中的源项进行了修正，表示为如下形式：

能量方程源项修正形式：

动量方程源项修正形式：

其中：

表示模糊区域流体设计参数，一般取值为

，

为小于0.0001的固定参数。

4.3.3基于Fluent的石蜡凝固和融化设置模拟研究

本章将基于Fluent流体仿真软件的凝固-融化模型对圆管外侧石蜡蓄能材料的融化和凝固过程进行数值模拟，在此基础上对不同结构形式，不同热源（恒定进水温度以及太阳能辐射热源）条件下的液相分数以及对应相变材料的温度场数据进行了提取，以此为进一步研究石蜡材料的融化过程提供初始理论数据。

如图5为本次设计所使用的石蜡相变材料融化模拟模型结构示意图（本次设计相变蓄热箱体结构总体容积为50L），本次模拟所采用的石蜡材料的相关物理特性如表1所示，同时为简化模拟，做出如下合理假设：

1.石蜡相变材料分布均匀，切各物理特性属各项同性；

2.石蜡液相区域属于牛顿不可压缩流体区域，不考虑可压缩性；

3.石蜡液相区域热对流属于二维层流流动模型；

4．模拟过程中仅仅考虑密度随温度的变化；

5.管道厚度与热阻系数在模拟过程中忽略不计。

图5圆管外石蜡相变材料相变过程模拟模型

表1石蜡材料物理特性参数表

名称

密度

（kg/m3）

比热

相变潜热

（kJ/kg）

固体相温度点

（K）

液体相温度点

（K）

动力粘度

（kg/（m.s））

导热系数（W/（m.K））

石蜡

760

2100

170

324.15

330.15

0.00324

0.25

4.3.3.1Fluent石蜡凝固和融化参数设置

1）进入ANSYS协同仿真平台，如图6所示，在工作窗口内建立Fluent流体仿真模块。

图6基于ANSYS的Fluent流体仿真模块

2）进入ANSYS仿真平台的建模模块（Geometry），如图7所示，在建模窗口内建立石蜡模型（容积50L，内部管道直径

，内部管道直径

，蓄热箱体高度

）。

图7石蜡箱体模拟模型

3）进入ANSYS仿真平台的网格划分模块，对石蜡模型进行网格划分，网格划分效果如图8所示（包括了58680个单元，63240个单元节点，网格平均质量0.93）。

图8石蜡模型网格划分效果

4）进入ANSYS仿真平台的Fluent仿真模块，设置求解器类型等相关设置，如图9，并依次对分析类型、模型材料属性、边界条件等相关参数进行设置。

图9Fluent基本设置窗口

●选择求解类型为瞬态求解，并设置求解器为压力基求解器，设置重力方向为+Z方向（大小为

），如图10所示;

图10求解设置

●激活Fluent融化-凝固模型，如图11所示；

图11Fluent凝固-融化模型激活

●在材料定义窗口对石蜡材料的相关物理属性进行设置（如表1，包括石蜡材料的密度、比热容参数、导入系数、潜热参数、动力粘度以及热膨胀参数等物理参数），如图12所示。

图12石蜡物理材料属性设置

●设置模型边界条件（热源温度、热流密度、热对流等相关参数），如图13所示（设置内壁温度为350K）；

图13模型边界条件设置

●选择求解器，并对求解存储子步等相关参数进行设置，最终求解计算。

不考虑石蜡箱体结构与外部环境的散热情况，石蜡相变材料随时间变化温度变化及液相体积分数变化整理如下：

1.如图14为入口水温80℃条件下，石蜡结构不同时刻（50min～500min）相变温度分布云图，对应液相体积分数分布云图如图15所示。

从图中可以看出，石蜡材料温度从靠近热水管管壁位置开始得到提升，当加热时间为500min时，整体最大相变材料温度基本达到50℃左右，同时液相体积分数与温度变化基本保持一致。

选取特点结果提取坐标点，提取温度变化曲线如图16所示。

图14不同时刻相变材料温度分布云图

图15不同时刻相变材料液相体积分数分布云图

图16特定位置点温度变化曲线（点坐标：

（60,0,0））

图17特定位置点温度变化曲线（点坐标：

（70,0,0））

图18特定位置点温度变化曲线（点坐标：

（80,0,0））