CN112487735A - 一种流体锂氧电池功率平衡模型算法 - Google Patents
一种流体锂氧电池功率平衡模型算法 Download PDFInfo
- Publication number
- CN112487735A CN112487735A CN202011390131.7A CN202011390131A CN112487735A CN 112487735 A CN112487735 A CN 112487735A CN 202011390131 A CN202011390131 A CN 202011390131A CN 112487735 A CN112487735 A CN 112487735A
- Authority
- CN
- China
- Prior art keywords
- lithium
- oxygen
- battery
- glassy carbon
- carbon electrode
- 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
Links
- QTJOIXXDCCFVFV-UHFFFAOYSA-N [Li].[O] Chemical compound [Li].[O] QTJOIXXDCCFVFV-UHFFFAOYSA-N 0.000 title claims abstract description 137
- 239000012530 fluid Substances 0.000 title claims abstract description 85
- 229910021397 glassy carbon Inorganic materials 0.000 claims abstract description 93
- 238000012360 testing method Methods 0.000 claims abstract description 27
- 150000002500 ions Chemical class 0.000 claims description 39
- 230000002572 peristaltic effect Effects 0.000 claims description 17
- 239000011244 liquid electrolyte Substances 0.000 claims description 13
- 239000007789 gas Substances 0.000 claims description 12
- 230000035699 permeability Effects 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 8
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 4
- 229910001416 lithium ion Inorganic materials 0.000 claims description 4
- 229910052760 oxygen Inorganic materials 0.000 claims description 4
- 239000001301 oxygen Substances 0.000 claims description 4
- 125000004432 carbon atom Chemical group C* 0.000 claims description 3
- 239000003365 glass fiber Substances 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 1
- 230000002708 enhancing effect Effects 0.000 claims 1
- 239000000126 substance Substances 0.000 abstract description 6
- 239000000203 mixture Substances 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000001311 chemical methods and process Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 7
- 229910052744 lithium Inorganic materials 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005086 pumping Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910006715 Li—O Inorganic materials 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000003792 electrolyte Substances 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- -1 redox mediators Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 229910001323 Li2O2 Inorganic materials 0.000 description 1
- OUUQCZGPVNCOIJ-UHFFFAOYSA-M Superoxide Chemical class [O-][O] OUUQCZGPVNCOIJ-UHFFFAOYSA-M 0.000 description 1
- JDZCKJOXGCMJGS-UHFFFAOYSA-N [Li].[S] Chemical compound [Li].[S] JDZCKJOXGCMJGS-UHFFFAOYSA-N 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009969 flowable effect Effects 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 230000002687 intercalation Effects 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 239000002086 nanomaterial Substances 0.000 description 1
- 239000005486 organic electrolyte Substances 0.000 description 1
- 150000002978 peroxides Chemical class 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Abstract
本发明涉及一种流体锂氧电池功率平衡模型算法，包括：步骤1、设计测试单体电池的测试装置；步骤2、建立锂氧流单体电池的输出功率的半经验模型。本发明的有益效果是：本发明提出的流体锂氧电池功率平衡模型针对复杂电池系统，在不需要获得动态数据的前提下，描述流体锂氧电池的特殊化学组成和反应过程；本发明设计一种合理的多孔玻璃态碳电极及测试该单体电池的测试装置，最大化输出功率且最小化功率损耗；半经验模型将拦截阴极流体的三维多孔阴极的几何结构与净功率输出进行函数关联，给出了净功率输出与锂氧流流速的关系方程，同时考虑单体电池的特殊化学过程和结构组成；充分开发了流体锂氧电池的化学特性。
Description
技术领域
本发明属于流体锂氧电池功率平衡领域，具体涉及一种流体锂氧电池功率平衡模型算法。
背景技术
随着可再生能源技术的发展，氧化还原流体锂氧电池(RFBs)作为电能存储系统(EES)越来越受到关注，其主要优点为设计灵活性，能量/功率解耦可行性，同时可以减少不常用电池组件的数量。
尽管目前人们使用不同的化学物质(催化剂、氧化还原介质、有机电解液、锂插层流动电极、微/纳米结构材料)来提高动力电池的性能，但多数研究成果主要集中在材料科学上。Arenas等人提出对流体锂氧电池的结构和实际性能进行工程研究，评价电池的整体功率平衡。材料科学应用中，电化学模型方法将推动锂离子电池、锂硫电池和流体锂氧电池等新型储能技术的发展。锂氧离子(Li2O2)是一种绝缘材料，在放电过程中，随着Li2O2在阴极上的形成，电池的过电位增大。Li2O2限制了电池的放电容量，导致充电电压上限高于3.7V。同时，在高化学反应性化合物(如过氧化物和超氧物)存在的情况下，高压充电会降低电解液和碳电极的化学稳定性，减少电池的循环寿命。Laoire和peng深入研究电解液组成对流体锂氧电池性能的影响程度。Monaco等人利用电池流体力学设计原理，显著提高电池的效率响应。
实际应用中，流动循环的能量消耗降低了锂流电池的整体效率。相应的功率损耗(Ploss)为锂电池的体积流量速率(V)和压降(Δp)的乘积，即：
上式中，η为流速效率。在设计良好的流动性RFB中，Ploss可以降到2％～4％。但是，使用粘性和半固态流体会增加功率损耗，其值不可忽略。Zeng等人证明了平均速度为v的流体通过长度为L的玻璃态碳(RVC)电极，Δp可以由Forchheimer方程模拟：
Δp/L＝(μ/k)·v+(ρβ/2)·v2
式中，μ和ρ分别是液体电解质的动态粘度和密度。k表示RVC磁导率，即对穿过多孔介质的流动粘性阻力。β为流体的惯性因子，表示路径弯曲引起的流体加速/减速的微观惯性效应。
发明内容
本发明的目的是克服现有技术中的不足，提供一种流体锂氧电池功率平衡模型算法。
这种流体锂氧电池功率平衡模型算法，包括以下步骤：
步骤1、设计测试单体电池的测试装置；设置测试压降条件，并测量不同阴极离子流速VRVC下，流体锂氧电池的实际输出功率；
步骤2、建立锂氧流单体电池的输出功率Pout的半经验模型，利用测试数据计算锂氧流单体电池的输出功率Pout和功率损耗Ploss；
步骤3、设计长度L、宽度W和厚度σ的棱柱状堆栈流体锂氧电池结构模型；
步骤4、计算棱柱状堆栈流体锂氧电池模型的输出功率Pout，建立棱柱状堆栈流体锂氧电池的输出功率Pout与体积的半经验拟合函数，电流大小、锂氧流单体电池的输出功率Pout和功率损耗Ploss与玻璃态碳电极(RVC)面积相关；计算功率损失Ploss和净功率PN；
步骤5、设计最大限度地提高锂氧流单体电池的输出功率Pout且减少锂氧流单体电池的功率损失Ploss的多孔玻璃态碳电极(RVC)：利用两片GF/F-Whatman玻璃纤维过滤器，分离多孔玻璃态碳电极(RVC)和锂离子电极，为了避免锂表面被溶解在阴极中的氧气污染，允许锂/多孔玻璃态碳电极(RVC)堆积层存在不对称位移。
作为优选，步骤1具体包括如下步骤：
步骤1.1、设置测试压降条件：增大测试装置中锂氧流单体电池(FLOB电池)内多孔玻璃态碳电极(RVC)的入口面积，减小多孔玻璃态碳电极(RVC)的出口面积，使流体锂氧电池的流量通道压降分布均匀；
步骤1.2、在多孔玻璃态碳电极(RVC)中不同阴极离子流速VRVC下，测量流体锂氧电池的实际输出功率和与流体锂氧电池的阴极离子流速VRVC相关的功率损失；对玻璃态碳电极(RVC)的磁导率k和惯性因子β进行实验评估。
作为优选，步骤2具体包括如下步骤：
步骤2.1、计算锂氧流单体电池的输出功率Pout，建立锂氧流单体电池的输出功率Pout的半经验模型：
将流体锂氧电池的实际输出功率与流体锂氧电池的阴极离子流速VRVC关联，设锂氧流单体电池的输出功率Pout只取决于阴极面积，阴极面积与玻璃态碳电极(RVC)总孔表面和体积相关；记录不同充放电电压条件下锂氧流单体电池的输出功率Pout与阴极离子流速VRVC，并将记录的锂氧流单体电池的输出功率Pout与阴极离子流速VRVC拟合得到关系曲线；利用插值算法建立锂氧流单体电池的输出功率Pout的半经验模型：
上式中，
为锂氧流单体电池的输出功率密度，A和P0为拟合常数，vτ为与时间常数相关的流速，vRVC为玻璃态碳电极(RVC)中的阴极离子流速；
步骤2.2、测得不同阴极体积流量速率V下通过流体锂氧电池的压降Δp，并记录该值；计算锂氧流单体电池的与阴极离子相关的功率损耗Ploss：
上式中，η为流速效率，Δp为压降，V为体积流量速率；则功率损失Ploss取决于沿流动方向穿过多孔玻璃态碳电极(RVC)的压降Δp和体积流量速率V；由Forchheimer方程得压降Δp取决于玻璃态碳电极(RVC)的磁导率k、惯性因子β和多孔玻璃态碳电极(RVC)的长度L；Forchheimer方程为：
Δp/L＝(μ/k)·v+(ρβ/2)·v2 (3)
上式中，μ为液体电解质的动态粘度；ρ是液体电解质的密度；k表示玻璃态碳电极(RVC)的磁导率，磁导率k为对穿过多孔介质的流动粘性阻力；β为流体的惯性因子，表示路径弯曲引起的流体加速/减速的微观惯性效应。
作为优选，步骤4具体包括如下步骤：
步骤4.1、建立棱柱状堆栈流体锂氧电池的输出功率Pout与体积的半经验拟合函数，在给定阴极离子流速VRVC条件下，将棱柱状堆栈流体锂氧电池的输出功率Pout与玻璃态碳电极(RVC)的体积关联：
上式中，Pout为棱柱状堆栈流体锂氧电池的输出功率，W为玻璃态碳电极(RVC)的宽度，L为玻璃态碳电极(RVC)的长度，σ为玻璃态碳电极(RVC)的厚度，A和P0为拟合常数，vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速，vτ为与时间常数相关的流速；
为锂氧流单体电池的输出功率密度，计算方式为：
步骤4.2、建立体积流量速率
与多孔玻璃态碳电极(RVC)尺寸的关系方程：
上式中，vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速，ε为多孔玻璃态碳电极(RVC)的空隙率；σ为玻璃态碳电极(RVC)的厚度，W为玻璃态碳电极(RVC)的宽度；
步骤4.3、结合
Forchheimer方程和式(5)计算功率损失Ploss和净功率PN：
上式中，A和P0为拟合常数；
为拟合常数P0的平均值；Δp为压降；
为多孔玻璃态碳电极的体积流量速率；μ为液体电解质的动态粘度；ρ是液体电解质的密度；β为流体的惯性因子；vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速；W为玻璃态碳电极(RVC)的宽度，L为玻璃态碳电极(RVC)的长度，σ为玻璃态碳电极(RVC)的厚度，ε为多孔玻璃态碳电极(RVC)的空隙率；Pout为锂氧流单体电池的输出功率；
上式(7)表明净功率PN随着多孔玻璃态碳电极(RVC)的体积增大而增加；在恒定多孔玻璃态碳电极(RVC)体积条件下改变锂氧流单体电池(FLOB电池)的几何结构，对增强锂氧流单体电池的输出功率Pout与功率损失Ploss没有影响，因为两个参数具有相同的线性依赖性(W·L·σ)；上式(4)至上式(6)表明锂氧流单体电池的功率损失Ploss和阴极离子流速vRVC在指数和多项式方面存在不同的增长速率。
作为优选，步骤1中测试装置包括储气罐、流体锂氧电池、压力计和蠕动泵；储气罐一端连接蠕动泵一端，蠕动泵另一端连接流体锂氧电池一端，流体锂氧电池另一端连接储气罐另一端，形成回路；压力计的两端并联在流体锂氧电池两端上。
作为优选，步骤1中单体电池的充放电电压范围为2.3V～2.7V。
作为优选，步骤2.1中计算锂氧流单体电池的输出功率Pout的方式为：Pout＝电流*电池的充放电电压。
本发明的有益效果是：本发明提出的流体锂氧电池功率平衡模型针对复杂电池系统，在不需要获得动态数据的前提下，描述流体锂氧电池的特殊化学组成和反应过程；本发明设计一种合理的多孔玻璃态碳电极(RVC)及测试该单体电池的测试装置，最大化输出功率且最小化功率损耗；半经验模型将拦截阴极流体的三维多孔阴极的几何结构与净功率输出进行函数关联，给出了净功率输出与锂氧流流速的关系方程，同时考虑单体电池的特殊化学过程和结构组成；充分开发了流体锂氧电池的化学特性。
附图说明
图1为测试单体电池的实验装置图；
图2为输出功率密度与阴极离子流速关系曲线图；
图3为长度L、宽度W和厚度σ的棱柱状堆栈电池结构模型图；
图4为在蠕动泵泵速为200rpm时(阴极离子流速不同)，不同电池电压(2.6V、2.5V、2.4V和2.3V)FLOB放电电流分布曲线图；
图5为输出功率Pout与单体电池电压的关系曲线图；
图6为不同泵速条件下输出功率Pout和损耗功率Ploss的比较曲线图；
图7中(a)为流体速度vRVC与蠕动泵泵速关系图；(b)为不同几何形状(g1,g2和g3)条件下，Pout与蠕动泵泵速的关系曲线图；(c)为不同几何形状(g1,g2和g3)条件下，Ploss与蠕动泵泵速的关系曲线图；(d)为净功率PN随蠕动泵泵速变化曲线图。
附图标记说明：储气罐1、流体锂氧电池2、压力计3、蠕动泵4、外部气体5、内部气体6。
具体实施方式
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
本发明提出流体锂氧电池实际功率性能的实验评估方法和建模方法。设计一种新型半经验模型，建立电池净功率输出与锂氧流流速的关系方程。
作为一种实施例，一种流体锂氧电池功率平衡模型算法，具体为：
步骤1、如图1所示，设计测试单体电池的测试装置，该测试装置包括储气罐1、锂氧流单体电池(FLOB电池)2、压力计3和蠕动泵4；储气罐1一端连接蠕动泵4一端，蠕动泵4另一端连接锂氧流单体电池2一端，锂氧流单体电池2另一端连接储气罐1另一端，形成回路；压力计3的两端并联在锂氧流单体电池2两端上；
步骤1.1、设置测试压降条件：增大测试装置中锂氧流单体电池(FLOB电池)2内多孔玻璃态碳电极(RVC)的入口面积，减小多孔玻璃态碳电极(RVC)的出口面积，使锂氧流单体电池(FLOB电池)2的流量通道压降分布均匀；
步骤1.2、在多孔玻璃态碳电极(RVC)中不同阴极离子流速VRVC下，测量流体锂氧电池的实际输出功率和与流体锂氧电池的阴极离子流速VRVC相关的功率损失；对玻璃态碳电极(RVC)的磁导率k和惯性因子β进行实验评估；
步骤2、利用测试数据估计锂氧流单体电池的输出功率Pout和功率损耗Ploss；
步骤2.1、计算锂氧流单体电池的输出功率Pout，建立锂氧流单体电池的输出功率Pout的半经验模型：
将流体锂氧电池的实际输出功率与流体锂氧电池的阴极离子流速VRVC关联，设锂氧流单体电池的输出功率Pout只取决于阴极面积，阴极面积与玻璃态碳电极(RVC)总孔表面和体积相关；记录不同充放电电压条件下锂氧流单体电池的输出功率Pout与阴极离子流速VRVC，并将记录的锂氧流单体电池的输出功率Pout与阴极离子流速VRVC拟合得到关系曲线，图2为不同充放电电压条件下的实验数据与拟合数据(实点为实验数据，实线为拟合数据的拟合曲线)；利用插值算法建立锂氧流单体电池的输出功率Pout的半经验模型：
上式中，
为锂氧流单体电池的输出功率密度，A和P0为拟合常数，vτ为与时间常数相关的流速，vRVC为玻璃态碳电极(RVC)中的阴极离子流速；
步骤2.2、由压力计3测得不同阴极体积流量速率V下通过流体锂氧电池的压降Δp，并记录该值；计算锂氧流单体电池的与阴极离子相关的功率损耗Ploss：
上式中，η为流速效率，Δp为压降，V为体积流量速率；则功率损失Ploss取决于沿流动方向穿过多孔玻璃态碳电极(RVC)的压降Δp和体积流量速率V；由Forchheimer方程得压降Δp取决于玻璃态碳电极(RVC)的孔隙度磁导率k、惯性因子β和多孔玻璃态碳电极(RVC)的长度L；Forchheimer方程为：
Δp/L＝(μ/k)·v+(ρβ/2)·v2 (3)
上式中，μ为液体电解质的动态粘度；ρ是液体电解质的密度；k表示玻璃态碳电极(RVC)的磁导率，磁导率k为对穿过多孔介质的流动粘性阻力；β为流体的惯性因子，表示路径弯曲引起的流体加速/减速的微观惯性效应；
步骤3、如图3所示，设计长度L、宽度W和厚度σ的棱柱状堆栈流体锂氧电池结构模型；设定三种RVC几何形状记为g1、g2和g3，三种几何形状具有相同的厚度、体积和横截面积，根据不同长宽比定义g1为细长几何体，g2为正方形几何体，g3为大几何体。三种几何形状的长宽比分别为，g1：L＝10cm和W＝1cm(细长几何体)，g2：L＝W＝3.16cm(正方形几何体)，g3：L＝1cm和W＝10cm(大几何体)；
步骤4、计算棱柱状堆栈流体锂氧电池模型的功率值，电流大小、锂氧流单体电池的输出功率Pout和功率损耗Ploss与玻璃态碳电极(RVC)面积相关；
步骤4.1、建立棱柱状堆栈流体锂氧电池的输出功率Pout与体积的半经验拟合函数，在给定阴极离子流速VRVC条件下，将棱柱状堆栈流体锂氧电池的输出功率Pout与玻璃态碳电极(RVC)的体积关联：
上式中，Pout为棱柱状堆栈流体锂氧电池的输出功率，W为玻璃态碳电极(RVC)的宽度，L为玻璃态碳电极(RVC)的长度，σ为玻璃态碳电极(RVC)的厚度，A和P0为拟合常数，vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速，vτ为与时间常数相关的流速；
为锂氧流单体电池的输出功率密度，计算方式见式(1)；
步骤4.2、建立体积流量加速度
与多孔玻璃态碳电极(RVC)尺寸的关系方程：
上式中，vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速，ε为多孔玻璃态碳电极(RVC)的空隙率；σ为玻璃态碳电极(RVC)的厚度，W为玻璃态碳电极(RVC)的宽度；
步骤4.3、结合式(2)、式(3)和式(5)计算功率损失Ploss和净功率PN：
上式中，A和P0为拟合常数；
为拟合常数P0的平均值；Δp为压降；μ为液体电解质的动态粘度；ρ是液体电解质的密度；β为流体的惯性因子；vRVC为多孔玻璃态碳电极(RVC)中的阴极离子流速；W为玻璃态碳电极(RVC)的宽度，L为玻璃态碳电极(RVC)的长度，σ为玻璃态碳电极(RVC)的厚度，ε为多孔玻璃态碳电极(RVC)的空隙率；Pout为锂氧流单体电池的输出功率；
上式(7)表明净功率PN随着多孔玻璃态碳电极(RVC)的体积增大而增加；在恒定多孔玻璃态碳电极(RVC)体积条件下改变锂氧流单体电池(FLOB电池)的几何结构，对增强锂氧流单体电池的输出功率Pout与功率损失Ploss没有影响，因为两个参数具有相同的线性依赖性(W·L·σ)；上式(4)至上式(6)表明锂氧流单体电池的功率损失Ploss和阴极离子流速vRVC在指数和多项式方面存在不同的增长速率；
步骤5、设计一种最大限度地提高锂氧流单体电池的输出功率Pout且减少锂氧流单体电池的功率损失Ploss的多孔玻璃态碳电极(RVC)：利用两片GF/F-Whatman玻璃纤维过滤器，分离多孔玻璃态碳电极(RVC)和锂离子电极，为了避免锂表面被溶解在阴极中的氧气污染，允许锂/多孔玻璃态碳电极(RVC)堆积层存在不对称位移。
实验结果1：
通过蠕动泵4控制泵速为200rpm时，电池电压分别为2.6V、2.5V、2.4V和2.3V的锂氧流单体电池(FLOB电池)的电流分布曲线如图4所示；由于氧在电极上的传输质量得到改善，流量的增加能够改善不同电压数值下的电流值；在给定的流速下，电流随电池极化而增加。
实验结果2：
控制泵速分别为0rmp、50rmp、100rmp和150rmp，通过将电流乘以电池电压来计算锂氧流单体电池的输出功率Pout，得到如图5所示锂氧流单体电池的输出功率Pout与电池电压的关系图；由图可知：
1)泵速越大，锂氧流单体电池的输出功率Pout越大；
2)在相同泵速条件下，输出功率Pout越大，电池电压越小；
3)若泵速增大，则输出功率随着电压增大减小更快。
实验结果3：
得到不同泵速条件下输出功率Pout和损耗功率Ploss的比较曲线如图6所示，图6中实折线为功率损耗Ploss，4条虚曲线为不同电池电压的输出功率Pout。结果表明，本实施例的单体电池具有正能量平衡(Pout>Ploss)的可行性区域，该可行性区域在0到30rpm之间。
实验结果4：
考虑流体速度与蠕动泵速率关系，图7(a)表示对应的流体速度vRVC，图7(b)和图6(c)表示不同的几何形状(g1,g2和g3)条件下，Pout和Ploss与泵速的关系曲线。由图可知：
1)流速增加(体积从g3变到g1)，Pout和Ploss同时增加；
2)几何形状g1的输出功率Pout最高，可以达到20MW，几何形状g3的功率损耗Ploss最小；
3)Pout和Ploss的变化趋势随着锂氧流体的流速变化不同而不同；(4)图7(d)为净功率PN随流速变化曲线。由图可知，PN在vRVC＝0m/s时达到最大值，并随RVC的体积大小不同而不同。
4)实验证明，在蠕动泵速率范围内，最大化净功率增益的最佳几何结构是具有最大流量横截面积的几何结构g3。
Claims (7)
1.一种流体锂氧电池功率平衡模型算法，其特征在于，包括以下步骤：
步骤1、设计测试单体电池的测试装置；设置测试压降条件，并测量不同阴极离子流速VRVC下，流体锂氧电池的实际输出功率；
步骤2、建立锂氧流单体电池的输出功率Pout的半经验模型，利用测试数据计算锂氧流单体电池的输出功率Pout和功率损耗Ploss；
步骤3、设计长度L、宽度W和厚度σ的棱柱状堆栈流体锂氧电池结构模型；
步骤4、计算棱柱状堆栈流体锂氧电池模型的输出功率Pout，建立棱柱状堆栈流体锂氧电池的输出功率Pout与体积的半经验拟合函数；计算功率损失Ploss和净功率PN；
步骤5、设计最大限度地提高锂氧流单体电池的输出功率Pout且减少锂氧流单体电池的功率损失Ploss的多孔玻璃态碳电极：利用两片GF/F-Whatman玻璃纤维过滤器，分离多孔玻璃态碳电极和锂离子电极。
2.根据权利要求1所述流体锂氧电池功率平衡模型算法，其特征在于，步骤1具体包括如下步骤：
步骤1.1、设置测试压降条件：增大测试装置中锂氧流单体电池内多孔玻璃态碳电极的入口面积，减小多孔玻璃态碳电极的出口面积；
步骤1.2、在多孔玻璃态碳电极中不同阴极离子流速VRVC下，测量流体锂氧电池的实际输出功率。
3.根据权利要求1所述流体锂氧电池功率平衡模型算法，其特征在于：步骤2具体包括如下步骤：
步骤2.1、计算锂氧流单体电池的输出功率Pout，建立锂氧流单体电池的输出功率Pout的半经验模型：
将流体锂氧电池的实际输出功率与流体锂氧电池的阴极离子流速VRVC关联，设锂氧流单体电池的输出功率Pout只取决于阴极面积，阴极面积与玻璃态碳电极总孔表面和体积相关；记录不同充放电电压条件下锂氧流单体电池的输出功率Pout与阴极离子流速VRVC，并将记录的锂氧流单体电池的输出功率Pout与阴极离子流速VRVC拟合得到关系曲线；利用插值算法建立锂氧流单体电池的输出功率Pout的半经验模型：
上式中，
为锂氧流单体电池的输出功率密度，A和P0为拟合常数，vτ为与时间常数相关的流速，vRVC为玻璃态碳电极中的阴极离子流速；
步骤2.2、测得不同阴极体积流量速率V下通过流体锂氧电池的压降Δp，并记录该值；计算锂氧流单体电池的与阴极离子相关的功率损耗Ploss：
上式中，η为流速效率，Δp为压降，V为体积流量速率；则功率损失Ploss取决于沿流动方向穿过多孔玻璃态碳电极的压降Δp和体积流量速率V；由Forchheimer方程得压降Δp取决于玻璃态碳电极的磁导率k、惯性因子β和多孔玻璃态碳电极的长度L；Forchheimer方程为：
Δp/L＝(μ/k)·v+(ρβ/2)·v2 (3)
上式中，μ为液体电解质的动态粘度；ρ是液体电解质的密度；k表示玻璃态碳电极的磁导率，磁导率k为对穿过多孔介质的流动粘性阻力；β为流体的惯性因子。
4.根据权利要求1所述流体锂氧电池功率平衡模型算法，其特征在于：步骤4具体包括如下步骤：
步骤4.1、建立棱柱状堆栈流体锂氧电池的输出功率Pout与体积的半经验拟合函数，在给定阴极离子流速VRVC条件下，将棱柱状堆栈流体锂氧电池的输出功率Pout与玻璃态碳电极的体积关联：
上式中，Pout为棱柱状堆栈流体锂氧电池的输出功率，W为玻璃态碳电极的宽度，L为玻璃态碳电极的长度，σ为玻璃态碳电极的厚度，A和P0为拟合常数，vRVC为多孔玻璃态碳电极中的阴极离子流速，vτ为与时间常数相关的流速；
为锂氧流单体电池的输出功率密度，计算方式为：
步骤4.2、建立体积流量速率
与多孔玻璃态碳电极尺寸的关系方程：
上式中，vRVC为多孔玻璃态碳电极中的阴极离子流速，ε为多孔玻璃态碳电极的空隙率；σ为玻璃态碳电极的厚度，W为玻璃态碳电极的宽度；
步骤4.3、结合
Forchheimer方程和式(5)计算功率损失Ploss和净功率PN：
上式中，A和P0为拟合常数；
为拟合常数P0的平均值；Δp为压降；μ为液体电解质的动态粘度；ρ是液体电解质的密度；β为流体的惯性因子；vRVC为多孔玻璃态碳电极中的阴极离子流速；W为玻璃态碳电极的宽度，L为玻璃态碳电极的长度，σ为玻璃态碳电极的厚度，ε为多孔玻璃态碳电极的空隙率；Pout为锂氧流单体电池的输出功率；
上式(7)表明净功率PN随着多孔玻璃态碳电极的体积增大而增加；在恒定多孔玻璃态碳电极体积条件下改变锂氧流单体电池的几何结构，对增强锂氧流单体电池的输出功率Pout与功率损失Ploss没有影响上式(4)至上式(6)表明锂氧流单体电池的功率损失Ploss和阴极离子流速vRVC在指数和多项式方面存在不同的增长速率。
5.根据权利要求1所述流体锂氧电池功率平衡模型算法，其特征在于：步骤1中测试装置包括储气罐(1)、流体锂氧电池(2)、压力计(3)和蠕动泵(4)；储气罐(1)一端连接蠕动泵(4)一端，蠕动泵(4)另一端连接流体锂氧电池(2)一端，流体锂氧电池(2)另一端连接储气罐(1)另一端，形成回路；压力计(3)的两端并联在流体锂氧电池(2)两端上。
6.根据权利要求1所述流体锂氧电池功率平衡模型算法，其特征在于：步骤1中单体电池的充放电电压范围为2.3V～2.7V。
7.根据权利要求3所述流体锂氧电池功率平衡模型算法，其特征在于：步骤2.1中计算锂氧流单体电池的输出功率Pout的方式为：Pout＝电流*电池的充放电电压。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011390131.7A CN112487735B (zh) | 2020-12-02 | 2020-12-02 | 一种流体锂氧电池功率平衡模型方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011390131.7A CN112487735B (zh) | 2020-12-02 | 2020-12-02 | 一种流体锂氧电池功率平衡模型方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112487735A true CN112487735A (zh) | 2021-03-12 |
| CN112487735B CN112487735B (zh) | 2024-03-15 |
Family
ID=74938751
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011390131.7A Active CN112487735B (zh) | 2020-12-02 | 2020-12-02 | 一种流体锂氧电池功率平衡模型方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112487735B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20010038943A1 (en) * | 1997-04-04 | 2001-11-08 | Spillman David M. | Electrochemical cell having multiplate electrodes with differing discharge rate regions |
| CN103345561A (zh) * | 2013-07-16 | 2013-10-09 | 上海航天电源技术有限责任公司 | 一种具有多层片状结构的锂离子电池热仿真分析方法 |
| CN103460469A (zh) * | 2011-04-05 | 2013-12-18 | 布莱克光电有限公司 | 以水为主的电化学氢催化剂动力系统 |
| US20150099199A1 (en) * | 2013-10-03 | 2015-04-09 | Massachusetts Institute Of Technology | Battery with Heterogeneous Flow-Through Porous Electrodes |
| CN109033619A (zh) * | 2018-07-24 | 2018-12-18 | 浙江大学城市学院 | 一种18650型锂电池放电循环的瞬态温度模型建模方法 |
| CN111276960A (zh) * | 2019-05-13 | 2020-06-12 | 中国矿业大学 | 一种光-储直流微电网系统中储能模块预测控制方法 |
-
2020
- 2020-12-02 CN CN202011390131.7A patent/CN112487735B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20010038943A1 (en) * | 1997-04-04 | 2001-11-08 | Spillman David M. | Electrochemical cell having multiplate electrodes with differing discharge rate regions |
| CN103460469A (zh) * | 2011-04-05 | 2013-12-18 | 布莱克光电有限公司 | 以水为主的电化学氢催化剂动力系统 |
| CN103345561A (zh) * | 2013-07-16 | 2013-10-09 | 上海航天电源技术有限责任公司 | 一种具有多层片状结构的锂离子电池热仿真分析方法 |
| US20150099199A1 (en) * | 2013-10-03 | 2015-04-09 | Massachusetts Institute Of Technology | Battery with Heterogeneous Flow-Through Porous Electrodes |
| CN109033619A (zh) * | 2018-07-24 | 2018-12-18 | 浙江大学城市学院 | 一种18650型锂电池放电循环的瞬态温度模型建模方法 |
| CN111276960A (zh) * | 2019-05-13 | 2020-06-12 | 中国矿业大学 | 一种光-储直流微电网系统中储能模块预测控制方法 |
Non-Patent Citations (4)
| Title |
|---|
| 张彦辉;: "锂氟化碳-二氧化锰电池组热分析", 船电技术, no. 04, 15 April 2018 (2018-04-15) * |
| 汪秋婷;姜银珠;陆豪;: "基于UKF的18650锂离子电池健康状况估计", 电源技术, no. 03, 20 March 2016 (2016-03-20) * |
| 汪秋婷;戚伟;肖铎;: "18650型锂电池放电过程的瞬态热模型与电压模拟", 科技通报, no. 02, 29 February 2020 (2020-02-29) * |
| 王求生;王乘;任波;李利军;: "固体氧化物燃料电池三维模拟与性能分析", 能源工程, no. 01, 25 February 2007 (2007-02-25) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112487735B (zh) | 2024-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zeng et al. | A hierarchical interdigitated flow field design for scale-up of high-performance redox flow batteries | |
| Reed et al. | Performance of Nafion® N115, Nafion® NR-212, and Nafion® NR-211 in a 1 kW class all vanadium mixed acid redox flow battery | |
| Huang et al. | Vanadium redox flow batteries: Flow field design and flow rate optimization | |
| CN102324570B (zh) | 一种锂离子电池及其化成方法、制备方法 | |
| CN109446619B (zh) | 一种锂离子电池电极设计参数的优化方法 | |
| US9991545B2 (en) | Metal-air flow batteries using oxygen enriched electrolyte | |
| CN106505200A (zh) | 碳纳米管/石墨烯/硅复合锂电池负极材料及其制备方法 | |
| CN112909232A (zh) | 一种氟化钠浸渍包覆钒掺杂多孔结构焦磷酸铁钠正极材料及其制备方法 | |
| CN109659640A (zh) | 一种锂离子电池快速化成工艺 | |
| CN108155347A (zh) | 提升锂离子电池含镍正极材料首次库仑效率方法及其应用 | |
| CN112086975B (zh) | 一种协调多储能单元参与二次调频的优化调度方法 | |
| CN108110295B (zh) | 液流电池装置的能效优化控制方法及系统 | |
| CN107732200A (zh) | 一种利用光伏产业废料制备锂离子电池负极材料的方法 | |
| CN108493427A (zh) | 用于锂离子电池电极材料的微纳米Nb2O5粉体的制备方法 | |
| Yan et al. | CFD numerical study of a new crossed inverse Z flow field for PEMFC | |
| CN107204480A (zh) | 液流电池电解液参数确定方法及其系统、液流电池 | |
| CN112924867A (zh) | 多场耦合下的锂电池容量衰减计算方法 | |
| CN109713339B (zh) | 一种基于电流优化策略的液流电池系统控制方法 | |
| CN206301881U (zh) | 静态钒电池系统 | |
| CN112487735B (zh) | 一种流体锂氧电池功率平衡模型方法 | |
| Ji et al. | Electrochemomechanical performance of porous electrode incorporating binder network | |
| Xiao et al. | Experimental study on charge/discharge characteristics of zinc-nickel single-flow battery | |
| CN107887566B (zh) | 一种锂离子液流电池系统的驱动控制方法 | |
| CN116125285A (zh) | 基于单体间流量不均匀性的液流电池性能预测方法 | |
| CN106876813A (zh) | 一种锂离子电池预充方法 |
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 |