CN108846205B - 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法 - Google Patents

一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法 Download PDF

Info

Publication number
CN108846205B
CN108846205B CN201810617324.8A CN201810617324A CN108846205B CN 108846205 B CN108846205 B CN 108846205B CN 201810617324 A CN201810617324 A CN 201810617324A CN 108846205 B CN108846205 B CN 108846205B
Authority
CN
China
Prior art keywords
heat
phase change
total
reservoir
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810617324.8A
Other languages
English (en)
Other versions
CN108846205A (zh
Inventor
徐阳
郑章靖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201810617324.8A priority Critical patent/CN108846205B/zh
Publication of CN108846205A publication Critical patent/CN108846205A/zh
Application granted granted Critical
Publication of CN108846205B publication Critical patent/CN108846205B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Abstract

本发明公开了一种固‑液相变储热器的快速性能计算、预测方法及储热器设计方法，针对固‑液相变储热器的动态、非线性难点，提出无量纲充/放热时间，即实际充/放热时间与基准充/放热时间的比值。基准充/放热时间通过解析方式建立显示计算公式，可表征相变储热器的非线性特性。提出了时均载热流体侧出口温度和时均壁面温度，表征固‑液相变储热器的动态特性。在此基础上，通过实验和数值模拟方法建立无量纲充/放热时间与各影响参数的修正公式。进而可以开展固‑液相储热器的两类快速设计问题(已知总换热面积设计充/放热时间和已知充/放热时间设计总换热面积)。该方法具有较高的准确性，且具有通用性，易于推广到不同类型的结构形式。

Description

一种固-液相变储热器的快速性能计算、预测方法及储热器设 计方法
技术领域
本发明属于固-液相变储热器技术领域，具体涉及一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法。
背景技术
随着能源危机和环境问题的日益加剧，储热得到了越来越多的关注。储热不仅可以节约能源消耗，同时也可以解决能源供需之间的不匹配问题。目前，储热方式主要有3种：显热储热、潜热储热及化学储热。其中，相变储热具有比显热高数倍的储热密度；且相变储热过程处于恒温状态，具有比化学储热更好的稳定性。目前常见的相变储热方式是指固-液相变储热。
固-液相变储热器是固-液相变储热系统的关键设备之一，其性能好坏以及匹配问题是影响整个储热系统性能及成本的关键因素之一。常见的固-液相变储热器是通过载热流体(HTF)将热量带进或带出相变储热器；即在充热的时候，在储热器内实现热量从载热流体侧传递给相变材料侧(PCM)；而在放热的时候，在储热器内实现热量从相变材料侧传递给载热流体侧。因此，相变储热器与常规换热器具有类似的结构。两者最大的不同之处在于：常规换热器中一般有冷热两种载热流体，流体通过壁面进行换热；而相变储热器内的相变材料侧一般是密封在相变储热器内，所以不能流动。正因为如此，相变储热器的相变材料侧的换热效果一般不太理想。截至目前，针对固-液相变储热器内的传热机制和性能提升方法已经开展了大量的研究工作。例如，为了改善固-液相变材料的传热能力弱的问题，研究工作包括制备一些高导材料(如纳米颗粒等)与相变材料的复合相变材料、在相变材料中填充高导翅片、以及设计优化相变储热器的结构等。这些研究工作为促进相变储热技术的发展起到了很大的作用，因此相变储热技术正处于由实验室向实际工程推进的阶段。
在将相变储热技术推向实际工程应用的阶段，如何快速预测相变储热器性能以及快速设计相变储热器成为了制约的关键因素之一。普通换热器性能的快速预测和设计技术已经十分成熟。然而，普通换热器预测与设计方法并能不够直接应用于相变储热器，这是因为相变储热器的固-液相变过程是一个非线性的动态过程。这种非线性的动态过程给相变储热器的快速预测和设计提出了挑战。
发明内容
本发明提供一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法，能够解决上述固-液相变过程的非线性动态特性对固-液相变储热器性能预测与设计造成的问题。
为实现上述目的，本发明采用如下技术方案：
一种固-液相变储热器的快速性能计算方法，建立无量纲充/放热时间计算方法，计算公式如下：
Figure GDA0002432664940000021
(储热器结构、载热流体侧参数、初始温度) (1)
其中：ttotal为实际储热器的充/放热时间；
ttotal,0是实际储热器的简化模型所对应的基准充/放热时间，所述简化模型是具有解析解或近似解析解的模型，其表达式为：
ttotal,0＝f(Ste,δPCM,λPCM,ρPCM,cp,PCM) (2)
其中，Ste为相变材料的Stefan数，ρPCM为相变材料的密度，cp,PCM为相变材料的比热容，λPCM为相变材料的导热系数，δPCM为相变材料的厚度；
f(储热器结构、载热流体侧参数、初始温度)是与储热器结构、载热流体侧参数以及整个储热器的初始温度等因素有关的公式，其与储热器结构、载热流体侧参数以及初始温度等因素之间的关系通过数值模拟或实验方法建立。
优选的，
Figure GDA0002432664940000022
Figure GDA0002432664940000023
其中，Ste0为无量纲初始时间，Tinitial为整个储热器的初始温度，Tmelting是相变材料的熔点，β是相变材料的相变潜热，L为储热器的长度，δHTF为载热流体的厚度，a、b、c、m以及n是待定常数，通过数值模拟或实验方法进行确定。
优选的，所述简化模型为单相Stefan模型，
Figure GDA0002432664940000031
Figure GDA0002432664940000032
Figure GDA0002432664940000033
其中，Th是热源温度或冷源温度，对于充热过程为热源温度，对于放热过程为冷源温度。
优选的，Th＝Tinlet，Tinlet为载热流体入口温度。
优选的，Th＝0.5(Tinlet+Toutlet)，其中Toutlet为载热流体出口温度。
优选的，Th＝Twall,ave，
Figure GDA0002432664940000034
Figure GDA0002432664940000035
其中，Twall,ave为储热器的时均壁温，Toutlet,ave为载热流体时均出口温度，qave为储热器的时均热流密度，qm,HTF为载热流体的质量流量，cp,HTF为载热流体的比热容，hHTF为载热流体与壁面的对流换热系数，A为储热器的换热面积；熔化充热过程时，式(8)右边最后一项取负号，凝固放热过程时，取正号；
Figure GDA0002432664940000036
Q＝Qs+Ql (11)
Ql＝mPCMβ (12)
Figure GDA0002432664940000037
储热量Q为充/放热过程中相变材料吸收或释放出来的总热量，由相变潜热Ql和显热Qs组成，mPCM为相变材料的质量。
一种基于上述计算方法的固-液相变储热器性能预测方法，包括如下步骤：
(a)确定基本参数：确定相变材料及其热物性，包括密度ρPCM、比热容cp,PCM、导热系数λPCM、熔点Tmelting以及相变潜热β等；确定相变材料的质量mPCM，确定载热流体及其物性，包括密度ρHTF、比热容cp,HTF以及导热系数λHTF等；确定载热流体侧的工况参数，包括进口温度Tinlet以及质量流量qm,HTF或流速uHTF等；确定储热器换热面积A和储热器的长度L、载热流体侧通道高度δHTF以及相变材料厚度δPCM等；
(b)计算储热量：根据公式(12)计算潜热量，如果需要考虑显热，则先给定一个初始假定的显热Qs,init；
(c)给定初始假定的充/放热时间ttotal,init；
(d)根据估计的显热Qs,init和充/放热时间ttotal,init，按照公式(9)和(10)计算时均出口温度Toutlet,ave；
(e)根据公式(8)计算时均壁面温度Twall,ave；
(f)核算显热量：根据时均壁面温度Twall,ave，按照例如公式(13)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)；
(g)计算基准充/放热时间ttotal,0：根据例如公式(5)、(6)、(7)计算基准充/放热时间；
(h)计算ttotal：根据公式(3)计算ttotal，将计算得到的ttotal和假定的ttotal,init进行对比，如果满足差异小于规定值则预测完成；如果差异大于规定值，则将计算得到的ttotal作为新的估计值代入步骤(c)，重复步骤(c)到(h)直到两者的差异小于规定值，预测完成。
一种基于上述计算方法的固-液相变储热器设计方法，包括如下步骤：
(a)确定基本参数：选择相变材料，确定其基本热物性，包括密度ρPCM、比热容cp,PCM、导热系数λPCM、熔点Tmelting以及相变潜热β等；确定载热流体及其物性，包括密度ρHTF，比热容cp,HTF、以及导热系数λHTF等；确定载热流体侧的工况参数包括Tinlet以及流量qm,HTF或流速uHTF等，确定充/放热时间ttotal；
(b)计算储热量：根据公式(12)计算潜热量，如果需要考虑显热，则先给定一个初始假定的潜热Qs,init；
(c)给定初始假定的总换热面积Ainit；
(d)根据已知的充/放热时间ttotal及假定的总换热面积Ainit，按照公式(9)、(10)计算时均出口温度Toutlet,ave；
(e)根据给定的Ainit，并按照公式(8)计算时均壁面温度Twall,ave；
(f)核算显热量：根据时均壁面温度Twall,ave，按照例如公式(13)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)；
(g)计算基准充/放热时间ttotal,0：根据例如公式(5)、(6)、(7)计算基准充/放热时间；
(h)根据例如公式(3)计算t’total：将计算得到的t’total和给定的充/放热时间ttotal进行对比，如果满足差异小于规定值则认为设计完成；如果差异大于规定值，则更改Ainit，重复步骤(d)到(h)直到两者的差异小于规定值，设计完成。
有益效果：与现有方法相比，本发明具有如下有益技术效果：1、通过引入无量纲充/放热时间的概念解决了动态、非线性特性对固-液相变储热器性能快速预测的影响；2、提出了与无量纲充/放热时间相匹配的固-液相变储热器快速设计方法；3、所提出的相变储热器快速设计方法具有通用性，易于推广到不同类型的结构形式。
附图说明
图1为本发明的一种根据相变储热器换热面积来校核充/放热时间的流程示意图；
图2为本发明的一种根据相变储热器充/放热时间来设计换热面积的流程示意图；
图3为本发明涉及的一种矩形腔式相变储热器结构示意图；
图4为更改参数对矩形腔式相变储热器快速预测与设计结果的影响；
图5为一种矩形腔式相变储热器快速预测与设计结果的验证；
图中标号：1、相变材料侧；2、载热流体侧；3、壁面。
具体实施方式
下面针对相变储热器的两类设计问题，对本发明做进一步的详细说明，所述是对本发明的解释而不是限定。
对于固-液相变储热器而言，其性能主要涉及3个参数，分别是：储热量Q、充/放热时间ttotal和换热面积A。其中储热量Q一般是已知的，因为相变潜热比显热高几个数量级，所以相变储热量主要取决于相变材料的多少。所以相变储热器一般涉及两类设计问题：第一类是已知换热面积A来校核充/放热时间ttotal，第二类是已知充/放热时间ttotal来设计换热面积A。其中，第二类可以转换为第一类，即第二类可以通过假定A来校核ttotal。
因此，相变储热器快速设计的基础是可以对给定A情况下的相变储热器的性能(充/放热时间)进行快速预测。本专利提出了一种无量纲充/放热时间的快速预测方法，并在该快速预测方法的基础上提出了一种时均快速设计方法。无量纲充/放热时间定义为实际充/放热时间ttotal和基准充/放热时间ttotal,0的比值，即ttotal/ttotal,0。其中：ttotal,0是实际储热器的简化模型所对应的基准充/放热时间。简化模型是一些具有解析解或近似解析解的模型，例如单相Stefan模型。因此，通过解析或近似解析方式可以获得ttotal,0的显式计算公式。但是，ttotal,0对应的是简化模型，其估计的结果和实际值之间还有一定的差异。因此需要通过数值模拟或实验的方法对ttotal,0进行修正，最后得到ttotal/ttotal,0的修正关系式。由于ttotal,0是通过理论方式获得的结果，ttotal,0的计算式已经能够很大程度反映相变储热过程的动态非线性特性，因此建立ttotal/ttotal,0与各影响因素的关系式肯定比直接建立ttotal和不同影响参数之间的非线性关系式要容易很多。
在建立无量纲充/放热时间的快速预测方法时，ttotal,0的建立很重要。一般而言，ttotal,0与相变材料的物性(密度ρPCM、比热容cp,PCM以及导热系数λPCM)、相变材料的厚度δPCM以及Stefan数(Ste)等有关系，其表达式为：
ttotal,0＝f(Ste,δPCM,λPCM,ρPCM,cp,PCM) (1)
对于单相Stefan模型，ttotal,0的计算公式如下：
Figure GDA0002432664940000061
Figure GDA0002432664940000062
ρPCM为相变材料的密度，cp,PCM为相变材料的比热容；
在储热器的设计过程中，相变储热材料一般是事先选定的，所以物性是确定的。而对于厚度δPCM而言，当储热器的换热面积A已知时，即储热器结构已确定时，δPCM也是已知的；即使当换热面积A未知时，也可以先假设一个δPCM后进行校核。所以，Ste数的计算成为了决定ttotal,0的关键。Ste数的定义式如下：
Figure GDA0002432664940000063
其中，Tmelting是相变材料的熔点，β是相变潜热。Tmelting，β，cp,PCM这三个参数对于给定的相变材料是确定值。而Th是热源温度或冷源温度(充热过程为热源温度，放热过程为冷源温度)，其取值关系到ttotal,0的准确度。Th的计算方式可以有三种：(1)以载热流体入口温度Tinlet作为Th，即Th＝Tinlet。这种估算方式最简单，因为载热流体的入口温度可以认为是恒定值，不随时间而变化。但是，这种方式估计的Ste数计算得到的ttotal,0与实际ttotal之间的差异较大。(2)以载热流体入口Tinlet和出口温度Toutlet的平均值作为Th，即Th＝0.5(Tinlet+Toutlet)。这种计算结果比第一种方式准确，但是问题在于载热流体的出口温度Toutlet是时变的，因为相变储热器的换热性能是动态变化的；(3)最精准的Th应该是取为载热流体和相变材料之间的壁面温度Twall，但是该壁温也是随时间而变化的。为了解决载热流体出口温度和壁面温度随时间的变化对计算造成的难度，本专利提出了时均方法。
载热流体时均出口温度Toutlet,ave的计算公式为：
Figure GDA0002432664940000071
其中，qave为时均热流密度，qm,HTF为载热流体的质量流量，cp,HTF为载热流体的比热容。
时均壁温Twall,ave的计算公式为：
Figure GDA0002432664940000072
其中，hHTF为载热流体与壁面的对流换热系数，其值取决于载热流体物性和流速等。当熔化充热过程时，式(6)右边最后一项取负号，当凝固放热过程时，取正号。
在计算时均出口温度Toutlet,ave和时均壁面温度Twall,ave的时候，都涉及时均热流密度qave的概念，其计算公式为：
Figure GDA0002432664940000073
其中，储热量Q为充/放热过程中相变材料吸收或释放出来的总热量。其由相变潜热Ql和显热Qs两部分构成。一般而言，在相变储热器内，相变潜热比显热高1-2个数量级，因此在近似计算中，可以用相变潜热近似替代储热量，即：
Q＝Qs+Ql≈Ql (8)
相变潜热Ql的计算比较简单，取决于相变材料的质量mPCM和相变潜热β。对于储热器设计而言，相变材料的用量和种类是已知的，所以相变潜热是已知的。Ql的计算公式如下：
Ql＝mPCMβ (9)
显热Qs的计算方法则比较复杂，除了与相变材料的物性有关系外，还与相变储热器的结构和其壁面温度有关系，而壁面温度由前面的介绍可知是随时间变化的。如果采用时均估计方法，对于矩形腔相变储热器结构，Qs的计算公式如下：
Figure GDA0002432664940000081
至此，已经建立了基准充/放热时间ttotal,0的计算方法。下面介绍如何构建无量纲充/放热时间ttotal/ttotal,0的修正关系式。该关系式与储热器结构、载热流体侧参数以及初始温度等因素有关。以矩形腔式储热器为例，可以构建如公式(11)所示的计算公式。其中，Ste与载热流体侧的入口温度有关，Ste0与初始温度有关，L/δHTF与储热器的长度有关，而δPCM/δHTF与储热器的厚度有关。
Figure GDA0002432664940000082
其中，a、b、c、m以及n是待定常数，需要通过数值模拟或实验方法进行确定。具体而言，就是保证Ste0、L/δHTF以及δPCM/δHTF三个参数不变，先通过更改Ste的值，来获得对应的ttotal/ttotal,0，进而通过拟合方法确定ttotal/ttotal,0和Ste之间的关系，从而确定常数b的取值。同理可以确定c、m以及n的取值。最后确定a的取值。当然，公式(11)可以涵盖更多的影响参数；考虑的参数越多，公式(11)的准确度越高。另外，针对不同的储热器结构，例如管壳式储热器也可以建立与公式(11)类似的关系式。因此，无量纲充/放热时间的修正关系式(11)具有通用性。
至此，建立了完整的无量纲充/放热时间快速计算方法，根据该方法就可以开展相变储热器的两类设计方法。
结合图1，针对换热面积已知的情况下的相变储热器设计进行说明：
(a)确定基本参数。选择合适的相变材料，确定其基本热物性(包括密度、比热容、导热系数、熔点以及相变潜热等)，确定载热流体及其物性(包括密度，比热容以及导热系数等)，确定载热流体侧的工况参数(包括进口温度、流量等)，确定换热面积和储热器的其他几何参数。
(b)计算储热量。根据公式(9)计算潜热量，如果需要考虑显热，则需要先给定一个初始假定的显热Qs,init。
(c)给定初始假定的充/放热时间ttotal,init。
(d)根据估计的充/放热时间ttotal,init，按照公式(5)计算时均出口温度Toutlet,ave。
(e)根据公式(6)计算时均壁面温度Twall,ave。
(f)核算显热量。根据时均壁面温度Twall,ave，按照例如公式(10)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值(例如5％)，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)。
(g)计算基准充/放热时间ttotal,0。根据例如公式(2)计算基准充/放热时间。
(h)计算ttotal。根据例如公式(11)计算ttotal。
(i)将计算得到的ttotal和假定的ttotal,init进行对比，如果满足差异小于规定值(例如5％)则认为设计完成；如果差异大于规定值，则将计算得到的ttotal作为新的估计值代入步骤(c)，直到两者的差异小于规定值，设计完成。
结合图2，针对充/放热时间已知情况下的相变储热器设计进行说明：
(a)确定基本参数。选择合适的相变材料，确定其基本热物性(包括密度、比热容、导热系数、熔点以及相变潜热等)，确定载热流体及其物性(包括密度，比热容以及导热系数等)，确定载热流体侧的工况参数(包括进口温度、流量等)。
(b)计算储热量。根据公式(9)计算潜热量，如果需要考虑显热，则需要先给定一个初始假定的显热Qs,init。
(c)根据已知的充/放热时间ttotal，按照公式(5)计算时均出口温度Toutlet,ave。
(d)给定初始假定的总换热面积Ainit。
(e)计算时均壁面温度Twall,ave。根据给定的Ainit，按照公式(6)计算时均壁面温度Twall,ave。
(f)核算显热量。根据时均壁面温度Twall,ave，按照例如公式(10)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值(例如5％)，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)。
(g)计算基准充/放热时间ttotal,0。根据例如公式(2)计算基准充/放热时间。
(h)根据例如公式(11)计算t’total。
将计算得到的t’total和给定的充/放热时间ttotal进行对比，如果满足差异小于规定值(例如5％)则认为设计完成；如果差异大于规定值，则更改Ainit，重复步骤(d)到(h)直到两者的差异小于规定值，设计完成。
图3显示了一个给定换热面积的矩形相变储热器，利用本专利提供的方法对不考虑自然对流情况下的熔化过程进行快速预测与设计，并将快速设计结果与数值模拟结果进行对比。
(1)选定相变储热材料为石蜡，其基本物性为：密度为820kg·m-3，比热为2500J·kg-1·K-1，导热系数为0.195W·m-1·K-1，熔点为321.66K；选定载热流体为水，其基本物性为：密度为998.2kg·m-3，比热为4182J·kg-1·K-1，导热系数为0.6W·m-1·K-1；储热器宽度W已知为1m，载热流体侧通道高度δHTF＝0.02m。
图4显示了更改不同影响参数(储热器长度L，相变材料厚度δPCM，载热流体流速uHTF，以及载热流体入口温度Tinlet)时的无量纲时间差(定义为(ttotal,0-ttotal)/ttotal)与载热流体HTF进出口温差(Tinlet-Toutlet,ave)之间的变化关系。可以看出，此时，这些影响参数对无量纲时间的影响可以归纳为Tinlet-Toutlet,ave对无量纲时间的影响，所以，公式(11)可以简化为：
Figure GDA0002432664940000101
对图4中的(ttotal,0-ttotal)/ttotal随Tinlet-Toutlet,ave的变化关系进行线性拟合得到：
Figure GDA0002432664940000102
对公式(13)进行变换，可以得到无量纲充/放热时间的计算公式为：
Figure GDA0002432664940000103
图5显示了不同储热器长度L时，利用公式(14)快速预测的充热(熔化)时间与数值模拟的对比情况，可以看出，在无量纲长度L/δHTF＝10～100的变化范围内，快速预测结果与数值模拟结果的误差远小于5％，说明本专利提出的针对固-液相变储热器的快速性能预测与设计方法是完全可行，且具有较高的准确性。

Claims (7)

1.一种固-液相变储热器的快速性能计算方法，其特征在于，建立无量纲充/放热时间计算方法，计算公式如下：
Figure FDA0002447820740000011
其中：ttotal为实际储热器的充/放热时间；
ttotal,0是实际储热器的简化模型所对应的基准充/放热时间，所述简化模型是具有解析解或近似解析解的模型，其表达式为：
ttotal,0＝f(Ste,δPCM,λPCM,ρPCM,cp,PCM) (2)
其中，Ste为相变材料的Stefan数，ρPCM为相变材料的密度，cp,PCM为相变材料的比热容，λPCM为相变材料的导热系数，δPCM为相变材料的厚度；
f(储热器结构、载热流体侧参数、初始温度)是与储热器结构、载热流体侧参数以及整个储热器的初始温度等因素有关的公式，其与储热器结构、载热流体侧参数以及初始温度等因素之间的关系通过数值模拟或实验方法建立；
Figure FDA0002447820740000012
Figure FDA0002447820740000013
其中，Ste0为无量纲初始时间，Tinitial为整个储热器的初始温度，Tmelting是相变材料的熔点，β是相变材料的相变潜热，L为储热器的长度，δHTF为载热流体的厚度，a、b、c、m以及n是待定常数，通过数值模拟或实验方法进行确定。
2.根据权利要求1所述的一种固-液相变储热器的快速性能计算方法，其特征在于，所述简化模型为单相Stefan模型，
Figure FDA0002447820740000014
Figure FDA0002447820740000015
Figure FDA0002447820740000016
其中，Th是热源温度或冷源温度，对于充热过程为热源温度，对于放热过程为冷源温度。
3.根据权利要求2所述的一种固-液相变储热器的快速性能计算方法，其特征在于，Th＝Tinlet，Tinlet为载热流体入口温度。
4.根据权利要求2所述的一种固-液相变储热器的快速性能计算方法，其特征在于，Th＝0.5(Tinlet+Toutlet)，其中Tinlet为载热流体入口温度，Toutlet为载热流体出口温度。
5.根据权利要求2所述的一种固-液相变储热器的快速性能计算方法，其特征在于，Th＝Twall,ave，
Figure FDA0002447820740000021
Figure FDA0002447820740000022
其中，Tinlet为载热流体入口温度，Twall,ave为储热器的时均壁温，Toutlet,ave为载热流体时均出口温度，qave为储热器的时均热流密度，qm,HTF为载热流体的质量流量，cp,HTF为载热流体的比热容，hHTF为载热流体与壁面的对流换热系数，A为储热器的换热面积；熔化充热过程时，式(8)右边最后一项取负号，凝固放热过程时，取正号；
Figure FDA0002447820740000023
Q＝Qs+Ql (11)
Ql＝mPCMβ (12)
Figure FDA0002447820740000024
储热量Q为充/放热过程中相变材料吸收或释放出来的总热量，由相变潜热Ql和显热Qs组成，mPCM为相变材料的质量。
6.一种基于权利要求5所述计算方法的固-液相变储热器性能预测方法，其特征在于，包括如下步骤：
(a)确定基本参数：确定相变材料及其热物性，确定载热流体及其物性，确定载热流体侧的工况参数，确定储热器换热面积A及储热器的结构；
(b)计算储热量：根据公式(12)计算潜热量，如果需要考虑显热，则先给定一个初始假定的显热Qs,init；
(c)给定初始假定的充/放热时间ttotal,init；
(d)根据估计的显热Qs,init和充/放热时间ttotal,init，按照公式(9)和(10)计算时均出口温度Toutlet,ave；
(e)根据公式(8)计算时均壁面温度Twall,ave；
(f)核算显热量：根据时均壁面温度Twall,ave，按照公式(13)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)；
(g)计算基准充/放热时间ttotal,0：根据公式(5)、(6)、(7)计算基准充/放热时间；
(h)计算ttotal：根据公式(3)计算ttotal，将计算得到的ttotal和假定的ttotal,init进行对比，如果满足差异小于规定值则预测完成；如果差异大于规定值，则将计算得到的ttotal作为新的估计值代入步骤(c)，重复步骤(c)到(h)直到两者的差异小于规定值，预测完成。
7.一种基于权利要求5所述计算方法的固-液相变储热器设计方法，其特征在于，包括如下步骤：
(a)确定基本参数：确定相变材料及其热物性，确定载热流体及其物性，确定载热流体侧的工况参数，确定充/放热时间ttotal；
(b)计算储热量：根据公式(12)计算潜热量，如果需要考虑显热，则先给定一个初始假定的潜热Qs,init；
(c)给定初始假定的总换热面积Ainit；
(d)根据已知的充/放热时间ttotal及假定的总换热面积Ainit，按照公式(9)、(10)计算时均出口温度Toutlet,ave；
(e)根据给定的Ainit，并按照公式(8)计算时均壁面温度Twall,ave；
(f)核算显热量：根据时均壁面温度Twall,ave，按照公式(13)进行显热计算，并将计算值与Qs,init进行对比，如果两者的差异大于规定值，则重新(b)到(f)步骤，直到两者的误差小于规定值；如果不考虑显热则可以跳过该步骤直接进入步骤(g)；
(g)计算基准充/放热时间ttotal,0：根据公式(5)、(6)、(7)计算基准充/放热时间；
(h)根据公式(3)计算t’total：将计算得到的t’total和给定的充/放热时间ttotal进行对比，如果满足差异小于规定值则认为设计完成；如果差异大于规定值，则更改Ainit，重复步骤(d)到(h)直到两者的差异小于规定值，设计完成。
CN201810617324.8A 2018-06-15 2018-06-15 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法 Active CN108846205B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810617324.8A CN108846205B (zh) 2018-06-15 2018-06-15 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810617324.8A CN108846205B (zh) 2018-06-15 2018-06-15 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法

Publications (2)

Publication Number Publication Date
CN108846205A CN108846205A (zh) 2018-11-20
CN108846205B true CN108846205B (zh) 2020-08-11

Family

ID=64203009

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810617324.8A Active CN108846205B (zh) 2018-06-15 2018-06-15 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法

Country Status (1)

Country Link
CN (1) CN108846205B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782602B (zh) * 2019-02-02 2021-03-23 中国科学院过程工程研究所 一种相变储热系统运行的智能动态优化方法
CN111581742B (zh) * 2020-04-30 2022-12-13 西安交通大学 一种宽应用领域固体储热器的综合设计方法
CN111907936B (zh) * 2020-08-03 2022-07-29 松冷（武汉）科技有限公司 一种实现温度时间智能可视化的保温箱及方法
CN113239541B (zh) * 2021-05-12 2023-08-29 中国矿业大学 用于相变储热器性能强化的纵向翅片结构快速优化方法
CN113283199B (zh) * 2021-06-28 2022-04-22 中国人民解放军国防科技大学 含相变的空气预冷器设计方法、装置、计算机系统及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593503A (zh) * 2013-10-16 2014-02-19 中国运载火箭技术研究院 一种肋片式相变蓄热装置优化设计方法
JP2014152947A (ja) * 2013-02-05 2014-08-25 Taisei Corp 蓄熱量算出方法および蓄熱量算出装置
CN105004749A (zh) * 2015-07-03 2015-10-28 浙江大学 固-液相变材料熔化传热性能参数测试系统及其方法
CN107368952A (zh) * 2017-06-28 2017-11-21 西安交通大学 一种相变储热的性能强化及经济性评价方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10890383B2 (en) * 2014-01-21 2021-01-12 Drexel University Systems and methods of using phase change material in power plants

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014152947A (ja) * 2013-02-05 2014-08-25 Taisei Corp 蓄熱量算出方法および蓄熱量算出装置
CN103593503A (zh) * 2013-10-16 2014-02-19 中国运载火箭技术研究院 一种肋片式相变蓄热装置优化设计方法
CN105004749A (zh) * 2015-07-03 2015-10-28 浙江大学 固-液相变材料熔化传热性能参数测试系统及其方法
CN107368952A (zh) * 2017-06-28 2017-11-21 西安交通大学 一种相变储热的性能强化及经济性评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Development of Correlation for Melting Time of Phase Change Material in Latent Heat Storage Unit;Rathod ManishK et al.;《 Energy Procedia 》;20150831;第75卷;第2125-2130页 *
Lattice-Boltzmann方法计算多孔介质内固液相变问题;钱吉裕等;《自然科学进展》;20060430;第16卷(第4期);全文 *
Total solidification time of a liquid phase change material enclosed in cylindrical/spherical containers;Levent Bilir et al.;《Applied Thermal Engineering》;20041231;第27卷(第17-18期);全文 *

Also Published As

Publication number Publication date
CN108846205A (zh) 2018-11-20

Similar Documents

Publication Publication Date Title
CN108846205B (zh) 一种固-液相变储热器的快速性能计算、预测方法及储热器设计方法
Liu et al. Design of sensible and latent heat thermal energy storage systems for concentrated solar power plants: Thermal performance analysis
Daghigh et al. Theoretical and experimental analysis of thermal performance of a solar water heating system with evacuated tube heat pipe collector
Miyara Thermal performance investigation of several types of vertical ground heat exchangers with different operation mode
El Qarnia Numerical analysis of a coupled solar collector latent heat storage unit using various phase change materials for heating the water
Wang et al. Numerical optimization of the cooling effect of the bionic spider-web channel cold plate on a pouch lithium-ion battery
Tang et al. Film condensation heat transfer on a horizontal tube in presence of a noncondensable gas
Chiu et al. Submerged finned heat exchanger latent heat storage design and its experimental verification
Wang et al. Thermal performance of solar air collection-storage system with phase change material based on flat micro-heat pipe arrays
Cui et al. A new model for predicting performance of fin-and-tube heat exchanger under frost condition
Wang et al. A combined experimental and simulation study on charging process of Erythritol–HTO direct-blending based energy storage system
Raccanello et al. Simplified dynamic modeling of single-tank thermal energy storage systems
d'Entremont et al. Modeling of a thermal energy storage system based on coupled metal hydrides (magnesium iron–sodium alanate) for concentrating solar power plants
Yang et al. Criteria for performance improvement of a molten salt thermocline storage system
Shirazi et al. Thermal capacity effects in borehole ground heat exchangers
Wang et al. Influence of inlet structure on thermal stratification in a heat storage tank with PCMs: CFD and experimental study
Yang et al. Numerical analysis on the thermal behavior of high temperature latent heat thermal energy storage system
Zhu et al. Experimental validation of a semi-dynamic simplified model of active pipe-embedded building envelope
Zhang Modified computational methods using effective heat capacity model for the thermal evaluation of PCM outfitted walls
Yu et al. Study on flow and heat transfer characteristics of composite porous material and its performance analysis by FSP and EDEP
Ha et al. A numerical study of thermo-fluid phenomena in metal hydride beds in the hydriding process
Liu et al. Experimental and numerical studies on melting/solidification of PCM in a horizontal tank filled with graded metal foam
Chen et al. Experimental and numerical investigations of a lauric acid‐multichannel flat tube latent thermal storage unit
Zhu et al. Thermal performance revival of composite PCM for hybrid BTMSs by architecture and formula integrated optimization
Bie et al. A novel empirical heat transfer model for a solar thermal storage process using phase change materials

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Xu Yang

Inventor after: Zheng Zhangjing

Inventor before: Zheng Zhangjing

Inventor before: Xu Yang

GR01 Patent grant
GR01 Patent grant