CN110210089A - 考虑空气溶解度的液体活塞等效库容模型及其仿真算法 - Google Patents

考虑空气溶解度的液体活塞等效库容模型及其仿真算法 Download PDF

Info

Publication number
CN110210089A
CN110210089A CN201910422161.2A CN201910422161A CN110210089A CN 110210089 A CN110210089 A CN 110210089A CN 201910422161 A CN201910422161 A CN 201910422161A CN 110210089 A CN110210089 A CN 110210089A
Authority
CN
China
Prior art keywords
phase
liquid
storage capacity
gas
air
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.)
Granted
Application number
CN201910422161.2A
Other languages
English (en)
Other versions
CN110210089B (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.)
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
North China Electric Power University
Original Assignee
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
North China Electric Power University
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 State Grid Corp of China SGCC, State Grid Hebei Electric Power Co Ltd, North China Electric Power University filed Critical State Grid Corp of China SGCC
Priority to CN201910422161.2A priority Critical patent/CN110210089B/zh
Publication of CN110210089A publication Critical patent/CN110210089A/zh
Application granted granted Critical
Publication of CN110210089B publication Critical patent/CN110210089B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种考虑空气溶解度的液体活塞等效库容模型及其仿真算法。考虑空气溶解度的液体活塞等效库容模型是由气相区、液相的不溶气液体区与液相的库容区组成的热力学等值模型，并且用传质泵表示溶解速率。基于该模型的仿真算法包括以下步骤：A、确定运行工况，输入两相的初始热力学参数，设定仿真时间步长；B、计算液体活塞等效库容模型参数；C、利用等值模型和当前时刻的两相热力学参数进行热力学计算，得到下一时刻的两相热力学参数；D、修正两相热力学参数；E、仿真运算结束并输出结果。

Description

考虑空气溶解度的液体活塞等效库容模型及其仿真算法
技术领域
本发明属于压缩空气储能领域，涉及一种考虑空气溶解度的液体活塞等效库容模型及其仿真算法。
背景技术
液体活塞技术是一种新兴的等温压缩空气储能技术，它利用比热容大的液体提供近似恒温环境，并且通过气液两相直接接触以提高传热效率，实现空气的储能与发电过程接近于等温过程，从而降低热损失，提高***的储能发电效率。常用的气液两相直接接触的方法是液体喷雾、活塞底部注气和填料换热。直接接触换热的方式不仅强化相间传热速率，而且强化相间传质速率。由于液体活塞内压缩空气的压强变化范围大，因而空气溶解度对装置运行的影响更加显著。但是，现有的文献一般将液相假设成不溶气的液体，其推导的热力学模型没有考虑空气溶解度的影响，导致高压下的仿真结果与实际数据存在一定偏差。
发明内容
本发明的目的是解决背景技术中所提到的液体活塞热力学模型的仿真精度不足的问题，提出一种考虑空气溶解度的液体活塞等效库容模型及其仿真算法。本发明将可以溶解压缩空气的液体空隙等效为库容，考虑了影响库容的各种因素，并且利用库容对等值模型的气液两相热力学参数进行修正，提高了模型的仿真精度。
本发明提出的考虑空气溶解度的液体活塞等效库容模型是由气相区、液相的不溶气液体区与液相的库容区组成的热力学等值模型，并且用传质泵表示溶解速率。
本发明提出的基于考虑空气溶解度的液体活塞等效库容模型的仿真算法包括以下步骤：
A、确定运行工况，输入两相的初始热力学参数，设定仿真时间步长；
B、计算液体活塞等效库容模型参数；
C、利用等值模型和当前时刻的两相热力学参数进行热力学计算，得到下一时刻的两相热力学参数；
D、修正两相热力学参数；
E、仿真运算结束并输出结果。
所述液体活塞是一种等温压缩空气储能装置，其通过比热容大的液体提供近似恒定的温度环境，并且增大气液接触面积和接触时间，使空气在压缩和膨胀过程中无限接近于等温过程。
所述气相区是没有溶解于液相的空气体积总和，所述不溶气液体区是空气无法进入的液体体积总和，所述库容区是相同压强下的空气可以进入的液体空隙体积总和，液相的体积是不溶气液体区与库容区之和。
所述步骤A中的运行工况分为储能压缩过程和发电膨胀过程。
所述步骤B中的液体活塞等效库容模型参数是溶解速率、库容区体积；空气的溶解速率正比相间传热速率；库容区体积的计算方法为：
其中，Veq为库容区体积，rs为空气在该压强下的溶解度，气体i为氮气或者氧气，Vw为装置内液相体积，Kpi为气体i的溶解平衡常数，为气体i的有效间隙度，bi为气体i的VanderWaals体积，R为理想气体常数，T为热力学温度，p为压强。
所述步骤C中的热力学计算包括以下步骤：
(1)计算装置对空气做功量
Wa(k)＝p(k)dV
其中，p(k)为k时刻压强，dV为气相体积在单个时间步的变化量，储能压缩过程为正，发电膨胀过程为负。
(2)计算空气的内能变化量
其中，Ca,V为气相定容比热，ma为装置内气相总质量。
(3)计算相间传热量
Q(k)＝h(k)·s(k)·(Ta(k)+dta(k)-Tw(k))·dt
其中，h(k)为相间换热系数，s(k)为换热面积，Ta(k)、Tw(k)分别为k时刻的气相、液相修正温度，dt为时间步长。
(4)计算两相温度变化量
其中，dTa(k)、dTw(k)分别为k时刻的气相、液相温度变化量，maq(k)为k时刻气相区剩余的气体质量，Cw为液相比热，mw(k)为k时刻装置内液相总质量，Δmw为单个时间步内液体质量的变化量。
(5)计算两相温度
储能压缩过程：
发电膨胀过程：
其中，ta(k+1)、tw(k+1)分别为k+1时刻的气相与液相未修正温度，Tw0为注入液体活塞的液体温度。
(6)更新热力学参数
其中，Va(k+1)、Vw(k+1)为k+1时刻气相区与液相区体积。
所述步骤D中的修正两相热力学参数包括以下步骤：
(1)更新气相区剩余气体质量
其中，Veq(k+1)为k+1时刻液相的库容区体积。
(2)两相温度修正
当ma(k)＞ma(k+1)时：
当ma(k)＜ma(k+1)时：
附图说明
图1是液体活塞等效库容模型的结构图；
图2是仿真算法流程图；
图3是储能压缩过程结果图；
图4是发电膨胀过程结果图。
具体实施方式
下面结合具体实施例对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
如图2所示，采用本发明的考虑空气溶解度的液体活塞等效库容模型仿真算法，对表1所示的运行过程进行仿真模拟，时间步为9000，将得到的仿真结果与实际***运行数据进行对比，验证所述方法的适用性。
本发明的具体实施步骤如下：
A、确定运行工况，输入两相的初始热力学参数，设定仿真时间步长
B、计算液体活塞等效库容模型参数
C、利用等值模型和当前时刻的两相热力学参数进行热力学计算，得到下一时刻的两相热力学参数
1.计算装置对空气做功量；
2.计算空气的内能变化量；
3.计算相间传热量；
4.计算两相温度变化量；
5.计算两相温度；
6.更新热力学参数；
D、修正两相热力学参数
1.更新溶气空穴的库容体积与气相区剩余气体质量；
2.两相温度修正；
3.若迭代完毕，仿真运算结束，否则返回步骤三进行下一时刻的热力学计算；
E、仿真运算结束并输出结果
本发明实施例的计算结果如图3、4所示。

Claims (7)

1.考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述考虑空气溶解度的液体活塞等效库容模型是由气相区、液相的不溶气液体区与液相的库容区组成的热力学等值模型，并且用传质泵表示溶解速率；所述仿真算法包括以下步骤：
A、确定运行工况，输入两相的初始热力学参数，设定仿真时间步长；
B、计算液体活塞等效库容模型参数；
C、利用等值模型和当前时刻的两相热力学参数进行热力学计算，得到下一时刻的两相热力学参数；
D、修正两相热力学参数；
E、仿真运算结束并输出结果。
2.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述液体活塞是一种等温压缩空气储能装置，其通过比热容大的液体提供近似恒定的温度环境，并且增大气液接触面积和接触时间，使空气在压缩和膨胀过程中无限接近于等温过程。
3.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述气相区是没有溶解于液相的空气体积总和，所述不溶气液体区是空气无法进入的液体体积总和，所述库容区是相同压强下的空气可以进入的液体空隙体积总和，液相的体积是不溶气液体区与库容区之和。
4.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述步骤A中的运行工况分为储能压缩过程和发电膨胀过程。
5.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述步骤B中的液体活塞等效库容模型参数是溶解速率、库容区体积；空气的溶解速率正比相间传热速率；库容区体积的计算方法为：
其中：Veq为库容区体积，rs为空气在该压强下的溶解度，气体i为氮气或者氧气，Vw为装置内液相体积，Kpi为气体i的溶解平衡常数，为气体i的有效间隙度，bi为气体i的Van derWaals体积，R为理想气体常数，T为热力学温度，p为压强。
6.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述步骤C中的热力学计算包括以下步骤：
(1)计算装置对空气做功量
Wa(k)＝p(k)dV
其中：p(k)为k时刻压强，dV为气相体积在单个时间步的变化量，储能压缩过程为正，发电膨胀过程为负；
(2)计算空气的内能变化量
其中：Ca,V为气相定容比热，ma为装置内气相总质量；
(3)计算相间传热量
Q(k)＝h(k)·s(k)·(Ta(k)+dta(k)-Tw(k))·dt
其中：h(k)为相间换热系数，s(k)为换热面积，Ta(k)、Tw(k)分别为k时刻的气相、液相修正温度，dt为时间步长；
(4)计算两相温度变化量
其中：dTa(k)、dTw(k)分别为k时刻的气相、液相温度变化量，maq(k)为k时刻气相区剩余的气体质量，Cw为液相比热，mw(k)为k时刻装置内液相总质量，Δmw为单个时间步内液体质量的变化量；
(5)计算两相温度
储能压缩过程：
发电膨胀过程：
其中：ta(k+1)、tw(k+1)分别为k+1时刻的气相与液相未修正温度，Tw0为注入液体活塞的液体温度；
(6)更新热力学参数
其中：Va(k+1)、Vw(k+1)为k+1时刻气相区与液相区体积。
7.根据权利要求1所述的考虑空气溶解度的液体活塞等效库容模型及其仿真算法，其特征在于：所述步骤D中的修正两相热力学参数包括以下步骤：
(1)更新气相区剩余气体质量
其中：Veq(k+1)为k+1时刻液相的库容区体积；
(2)两相温度修正
当ma(k)＞ma(k+1)时：
当ma(k)＜ma(k+1)时：
。
CN201910422161.2A 2019-05-21 2019-05-21 考虑空气溶解度的液体活塞等效库容模型及其仿真算法 Expired - Fee Related CN110210089B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910422161.2A CN110210089B (zh) 2019-05-21 2019-05-21 考虑空气溶解度的液体活塞等效库容模型及其仿真算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910422161.2A CN110210089B (zh) 2019-05-21 2019-05-21 考虑空气溶解度的液体活塞等效库容模型及其仿真算法

Publications (2)

Publication Number Publication Date
CN110210089A true CN110210089A (zh) 2019-09-06
CN110210089B CN110210089B (zh) 2021-02-09

Family

ID=67787909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910422161.2A Expired - Fee Related CN110210089B (zh) 2019-05-21 2019-05-21 考虑空气溶解度的液体活塞等效库容模型及其仿真算法

Country Status (1)

Country Link
CN (1) CN110210089B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112364576A (zh) * 2020-11-09 2021-02-12 中冶赛迪技术研究中心有限公司 一种喷淋塔气液接触效果评价及优化设计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000041645A1 (en) * 1999-01-13 2000-07-20 Koninklijke Philips Electronics N.V. Fluid-dispensing and refilling system for a power toothbrush
CN1945275A (zh) * 2006-10-12 2007-04-11 浙江大学 一种液体中的气体溶解度测试仪
CN104156541A (zh) * 2014-08-26 2014-11-19 武汉大学 一种特高压断路器用液压操动机构的仿真平台***及方法
EP3374003A1 (en) * 2015-11-13 2018-09-19 Fresenius Vial SAS Method for detecting an occlusion in an infusion line
CN108644095A (zh) * 2018-04-18 2018-10-12 华北电力大学 基于分级压缩空气储能***的功率倍增运行策略
CN109340079A (zh) * 2018-09-17 2019-02-15 华北电力大学 一种等温压缩空气储能***及高效发电方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000041645A1 (en) * 1999-01-13 2000-07-20 Koninklijke Philips Electronics N.V. Fluid-dispensing and refilling system for a power toothbrush
CN1945275A (zh) * 2006-10-12 2007-04-11 浙江大学 一种液体中的气体溶解度测试仪
CN104156541A (zh) * 2014-08-26 2014-11-19 武汉大学 一种特高压断路器用液压操动机构的仿真平台***及方法
EP3374003A1 (en) * 2015-11-13 2018-09-19 Fresenius Vial SAS Method for detecting an occlusion in an infusion line
CN108644095A (zh) * 2018-04-18 2018-10-12 华北电力大学 基于分级压缩空气储能***的功率倍增运行策略
CN109340079A (zh) * 2018-09-17 2019-02-15 华北电力大学 一种等温压缩空气储能***及高效发电方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RAMAKRISHNAN: "Temperature Abatement Using Hollow Spheres in Liquid Piston Compressor for Ocean Compressed Air Energy Storage System", 《A THESIS SUBMITTED TO THE GRADUATE FACULTY OF NORTH CAROLINA STATE UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE》 *
吴云翔: "气液相变热声发动机声电类比分析与实验研究", 《浙江大学硕士学位论文》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112364576A (zh) * 2020-11-09 2021-02-12 中冶赛迪技术研究中心有限公司 一种喷淋塔气液接触效果评价及优化设计方法
CN112364576B (zh) * 2020-11-09 2024-02-27 中冶赛迪技术研究中心有限公司 一种喷淋塔气液接触效果评价及优化设计方法

Also Published As

Publication number Publication date
CN110210089B (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
Zhang et al. A near-isothermal expander for isothermal compressed air energy storage system
Houssainy et al. Thermodynamic performance and cost optimization of a novel hybrid thermal-compressed air energy storage system design
CN104279012B (zh) 一种基于深冷储能的核电调峰***
Mozayeni et al. Performance study of an advanced adiabatic compressed air energy storage system
Zhang et al. Combined cooling, heating, and power generation performance of pumped thermal electricity storage system based on Brayton cycle
CN104989473A (zh) 一种发电***以及基于该***的发电方法
CN108930627A (zh) 一种定压抽水压缩气体储能***及储能方法
Zhang et al. Comparative analysis of compressed carbon dioxide energy storage system and compressed air energy storage system under low-temperature conditions based on conventional and advanced exergy methods
CN110210089A (zh) 考虑空气溶解度的液体活塞等效库容模型及其仿真算法
Li et al. Experimental and thermodynamic investigation on isothermal performance of large-scaled liquid piston
CN108953099B (zh) 一种闭式等压压缩空气储能***及方法
Odukomaiya et al. Transient thermofluids analysis of a ground-Level integrated diverse energy storage (GLIDES) system
Li et al. Comprehensive thermo-exploration of a near-isothermal compressed air energy storage system with a pre-compressing process and heat pump discharging
Tong et al. Accumulated and transient exergy analyses of pneumatic systems with isochoric and isobaric compressed air storage tanks
US11391415B1 (en) Method for minimizing power demand for hydrogen refueling station
CN110032826A (zh) 一种联合Matlab与AMESim的间冷循环燃气轮机建模仿真方法
CN205618209U (zh) 一种用于燃气轮机试验台的综合气源***
He et al. Thermodynamic investigation of variable‐speed compression unit in near‐isothermal compressed air energy storage
CN216345433U (zh) 一种站内制氢用撬装增压、储存及加氢装置
CN105756727B (zh) 一种用于燃气轮机试验台的综合气源***
CN115263467A (zh) 单抽热电联产抽凝机组运行功率上下限确定方法及***
CN105019960B (zh) 电力***低谷电存储以及调峰控制***
Zhang et al. Performance analysis of advanced adiabatic compressed air energy storage system with constant-volume and adiabatic gas storage model under multi-cycle operation
CN113757072A (zh) 一种离子液体氢气压缩机及其使用方法
Kiliç et al. Mathematical modeling of a small scale compressed air energy storage system

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210209

Termination date: 20210521