CN114564866A - 一种热仿真网格划分方法 - Google Patents
一种热仿真网格划分方法 Download PDFInfo
- Publication number
- CN114564866A CN114564866A CN202210205222.1A CN202210205222A CN114564866A CN 114564866 A CN114564866 A CN 114564866A CN 202210205222 A CN202210205222 A CN 202210205222A CN 114564866 A CN114564866 A CN 114564866A
- Authority
- CN
- China
- Prior art keywords
- grid
- grids
- size
- simulation
- mesh
- 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.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 35
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 claims abstract description 17
- 229910001416 lithium ion Inorganic materials 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 description 10
- 238000004146 energy storage Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000000638 solvent extraction Methods 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
-
- 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/10—Energy storage using batteries
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Computer Graphics (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种热仿真网格划分方法,将锂离子电池模组仿真模型划分为电池、连接片、流场以及构件;确定电池网格、连接片网格、连接片网格以及构件尺寸;根据确定的尺寸对锂离子电池模组通过软件进行网格划分,可以有效降低网格数量,提升仿真效率,为后期仿真结果的应用提供便利。
Description
技术领域
本发明属于电池热管理技术领域,具体涉及一种热仿真网格划分方法。
背景技术
在锂离子电池模组应用过程中,电池热管理是储能系统稳定运行的重要保障。
目前电池热管理策略主要是通过热仿真结果制定,因此仿真结果的准确性关系到热管理策略的有效性及系统安全性。若要获得储能模组准确的热仿真结果,需要如下条件:1、准确的3D物理模型;2、准确的物理参数;3、正确的网格划分;4、适合的方程组。网格划分是其中的难点,网格划分是否适当,直接影响到仿真结果。网格划分与网格划分参数密切关系。不同的划分方法,可以获得不同的网格数量,网格数量直接影响到仿真结果的准确性,理论上网格数量越大,仿真结构越准确,但网格数量越大,所需要的计算资源越大,计算时间越长,因此如何在网格数量和计算资源之间寻找平衡点,一直是仿真工作的难点。
但是现有方法(如fluent mesh),在网格划分时,给出的最小网格尺寸偏小,并且在网格生成时,经常只是用一套网格标准,这套标准用于整个模型时,模型整体网数量偏大,计算量巨大,计算时间变长,影响仿真效率。
因此,适宜的网格划分方法将有效提高仿真精度,提升仿真精度,为后期仿真结果的应用提供便利。
发明内容
为克服现有技术中的问题,本发明的目的是提供一种热仿真网格划分方法。
为实现上述目的,本发明的采用的技术方案为:
一种热仿真网格划分方法,包括以下步骤:
将锂离子电池模组仿真模型划分为电池、连接片、流场以及构件;
确定电池网格、连接片网格、连接片网格以及构件网格尺寸;
根据所述电池网格、连接片网格、连接片网格以及构件网格尺寸对锂离子电池模组进行网格划分;
根据网格划分结果,实现对锂离子电池模组的热仿真。
进一步的,电池厚度为αmm时,电池网格尺寸最小值为α/N mm,最大值为αmm,N为正整数。
进一步的,连接片厚度为σmm时,连接片网格尺寸最小值为σ/M mm,最大值为σmm,M为正整数。
进一步的,电池间缝隙为βmm时,流场网格尺寸最小值为β/L mm,最大值为βmm,L为正整数。
进一步的,构件最小尺寸为θ时,构件网格尺寸最小值为θ/L mm,最大值为θmm,L为正整数。
进一步的,N+M+L=6。
进一步的,N=1,M=2,L=3,或N=1,M=3,L=2。
进一步的,N=2,M=1,L=3,或N=2,M=3,L=1。
进一步的,N=3,M=2,L=1,或N=3,M=1,L=2。
进一步的,根据所述电池网格、连接片网格、连接片网格以及构件网格尺寸对锂离子电池模组进行网格划分是通过fluent mesh软件实现的。
与现有技术相比,本发明具有的有益效果:
本发明通过对锂离子电池模组仿真模型划分,确定电池网格、连接片网格、连接片网格以及构件尺寸;根据确定的尺寸特征参量,通过软件实现模型的网格划分,实现储能模组网格数量的优化。采用本发明中的方法进行网格划分,可减少热仿真过程中的计算量,提升仿真效率。本发明的方法能够大幅降低网格数量,同时本发明的划分方法,能够将网格划分的时间,从119min,缩短到18min,节省了时间。
附图说明
图1为锂离子电池模组划分示意图。
图2为采用现有技术划分的网格示意图。
图3为采用本发明方法划分的网格示意图。
图4为不同网格划分方法仿真温度曲线。
图中,1为电池,2为连接片。
具体实施方式
下面结合附图对本发明进行详细说明。
本发明中涉及的术语含义,解释如下:
Fluent mesh:一款网格划分软件。
网格划分:通过有限元方法将物理模型划分为若干区域。
通过分析锂离子电池模组仿真模型,提取模型中特征参量,通过计算确定网格的尺寸,并通过软件实现模型的网格划分。
具体的,将锂离子电池模组仿真模型划分为电池1、连接片2、流场与其他部件(即构件,除电池、连接片与流场之外的部件的总称)4部分,其中,空间网格特征值主要与电池间隙密切相关,当电池1厚度αmm,连接片2厚度σmm,电池间缝隙βmm为已知量。
由此可以确定这三部分的网格尺寸最小值与最大值:
电池网格尺寸最小值为α/N mm,最大值为αmm;
连接片网格尺寸最小值为σ/M mm,最大值为σmm;
流场网格尺寸最小值为β/L mm,最大值为βmm;
其中,N,M,L为正整数。
在确定网格尺寸的过程中,本发明提供了一组计算公式:
N+M+L=6 (1)
f(N,M,L)=1 (2)
g(N,M,L)=0 (3)
通过上述公式可以计算出N,M和L数值,具体的,N=1,M=2,L=3,或N=1,M=3,L=2,或N=2,M=1,L=3,或N=2,M=3,L=1,或N=3,M=2,L=1,或N=3,M=1,L=2。
进而得到电池1,连接片2,流场的网格尺寸范围。
其他部件网格与流场尺寸关系较大,因此其他部件网格范围可与流场网格使用相同的设置。即其他部件网格尺寸最小值为θ/L mm,最大值为θmm;θ为其他部件(即构件)最小尺寸,可以根据实际情况确定。
由此确定了锂离子电池模组中4部分网格尺寸范围,以此设置通过软件fluentmesh进行网格划分,可以有效降低网格数量,提升仿真效率。
本发明通过确定储能模组中各部分的网格划分尺寸,由于各分部的网格尺寸相互关联,因此通过计算能够得到一套相对合理的网格尺寸数据,根据网格划分结果,实现对锂离子电池模组的热仿真。
参见图2,现有技术用软件默认参数划分网格,网格数量约2000万左右。
参见图3,在下面的算例里N=4,M=1,L=1,α、σ和β以及θ具体值分别为25mm,3mm,4mm,3mm,采用软件Fluent mesh划分网格,网格数量约为200万左右。通过图2与图3对比,可以看出,本发明的方法能够大幅降低网格数量,同时本发明的划分方法,将网格划分的时间,从119min,缩短到18min,节省了时间。
本发明对不同网格生成方法进行了对比,对比结果参见表1。
表1不同方法的网格数量和网格生成时间
仿真流程是建模、网格划分、边界条件设置、最后计算得到结果。对同一算例,模型、边界条件、计算都是相同的,只有网格可以修改,不同的网格划分会有不同的仿真结果。
参见图4,从曲线数据可以看出该网格划分方法仿真结果,与其他网格方划分方法仿真结果最大误差约0.6%。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种热仿真网格划分方法,其特征在于,包括以下步骤:
将锂离子电池模组仿真模型划分为电池、连接片、流场以及构件;
确定电池网格、连接片网格、连接片网格以及构件网格尺寸;
根据所述电池网格、连接片网格、连接片网格以及构件网格尺寸对锂离子电池模组进行网格划分;
根据网格划分结果,实现对锂离子电池模组的热仿真。
2.根据权利要求1所述的一种热仿真网格划分方法,其特征在于,电池厚度为αmm时,电池网格尺寸最小值为α/N mm,最大值为αmm,N为正整数。
3.根据权利要求2所述的一种热仿真网格划分方法,其特征在于,连接片厚度为σmm时,连接片网格尺寸最小值为σ/M mm,最大值为σmm,M为正整数。
4.根据权利要求3所述的一种热仿真网格划分方法,其特征在于,电池间缝隙为βmm时,流场网格尺寸最小值为β/L mm,最大值为βmm,L为正整数。
5.根据权利要求3所述的一种热仿真网格划分方法,其特征在于,构件最小尺寸为θ时,构件网格尺寸最小值为θ/L mm,最大值为θmm,L为正整数。
6.根据权利要求5所述的一种热仿真网格划分方法,其特征在于,N+M+L=6。
7.根据权利要求5或6所述的一种热仿真网格划分方法,其特征在于,N=1,M=2,L=3,或N=1,M=3,L=2。
8.根据权利要求5或6所述的一种热仿真网格划分方法,其特征在于,N=2,M=1,L=3,或N=2,M=3,L=1。
9.根据权利要求5或6所述的一种热仿真网格划分方法,其特征在于,N=3,M=2,L=1,或N=3,M=1,L=2。
10.根据权利要求1所述的一种热仿真网格划分方法,其特征在于,根据所述电池网格、连接片网格、连接片网格以及构件网格尺寸对锂离子电池模组进行网格划分是通过fluentmesh软件实现的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210205222.1A CN114564866A (zh) | 2022-03-02 | 2022-03-02 | 一种热仿真网格划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210205222.1A CN114564866A (zh) | 2022-03-02 | 2022-03-02 | 一种热仿真网格划分方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114564866A true CN114564866A (zh) | 2022-05-31 |
Family
ID=81717605
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210205222.1A Pending CN114564866A (zh) | 2022-03-02 | 2022-03-02 | 一种热仿真网格划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114564866A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117272763A (zh) * | 2023-11-22 | 2023-12-22 | 苏州易来科得科技有限公司 | 电池热力学仿真网格的划分方法,装置,存储介质和设备 |
CN117350088A (zh) * | 2023-12-06 | 2024-01-05 | 苏州易来科得科技有限公司 | 生成电池极片仿真网格的方法,装置,存储介质和设备 |
-
2022
- 2022-03-02 CN CN202210205222.1A patent/CN114564866A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117272763A (zh) * | 2023-11-22 | 2023-12-22 | 苏州易来科得科技有限公司 | 电池热力学仿真网格的划分方法,装置,存储介质和设备 |
CN117272763B (zh) * | 2023-11-22 | 2024-02-13 | 苏州易来科得科技有限公司 | 电池热力学仿真网格的划分方法,装置,存储介质和设备 |
CN117350088A (zh) * | 2023-12-06 | 2024-01-05 | 苏州易来科得科技有限公司 | 生成电池极片仿真网格的方法,装置,存储介质和设备 |
CN117350088B (zh) * | 2023-12-06 | 2024-02-23 | 苏州易来科得科技有限公司 | 生成电池极片仿真网格的方法,装置,存储介质和设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114564866A (zh) | 一种热仿真网格划分方法 | |
CN111398833B (zh) | 一种电池健康状态评估方法 | |
WO2017161786A1 (zh) | 基于光储系统运行优化的混合储能配比计算方法 | |
CN111276973A (zh) | 一种计及风功率波动的电力系统惯量需求评估的方法 | |
CN106128518A (zh) | 一种获取高精度的快中子反应堆组件少群截面的方法 | |
CN103106344A (zh) | 一种建立电力系统聚类负荷模型的方法 | |
CN114240019A (zh) | 适用新能源电力系统的灵活资源价值评估方法及装置 | |
CN113285472A (zh) | 一种储能系统的优化配置方法和装置 | |
WO2020010830A1 (zh) | 动力电池梯次利用的全寿命周期成本的估算方法 | |
CN113033136B (zh) | 简化的光伏电池物理参数提取优化方法及系统 | |
CN105262146A (zh) | 含风电的电力系统备用容量计算方法和系统 | |
CN106056305A (zh) | 一种基于状态聚类的发电系统可靠性快速评估方法 | |
CN113094878A (zh) | 一种用于铸造流动场数值模拟的并行计算方法及系统 | |
CN112736944A (zh) | 一种电化学储能电站有功功率调度方法和系统 | |
CN112564160A (zh) | 一种计及风电不确定性的储能系统随机配置方法、终端及存储介质 | |
CN113609686B (zh) | 一种新能源置信容量分析方法及系统 | |
CN114142470A (zh) | 电力系统新能源电站配套储能合理规模测算方法及系统 | |
CN114400660A (zh) | 一种风电光伏制氢系统配比确定方法及系统 | |
CN114188986A (zh) | 一种区域光伏最大消纳电量计算方法 | |
CN114243740B (zh) | 储能电站替代输电通道容量的分析方法及系统 | |
CN110031769A (zh) | 锂电池的电池包容量测算方法 | |
Xi et al. | A pseudo-analytical mix sampling strategy for reliability assessment of power girds | |
CN112100790B (zh) | 一种基于近似函数的高斯混合模型参数获取方法 | |
CN113805084B (zh) | 电池容量衰减的计算方法、装置、计算机设备及存储介质 | |
CN110912117B (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 |