CN110096774A - Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect - Google Patents

Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect Download PDF

Info

Publication number
CN110096774A
CN110096774A CN201910319320.6A CN201910319320A CN110096774A CN 110096774 A CN110096774 A CN 110096774A CN 201910319320 A CN201910319320 A CN 201910319320A CN 110096774 A CN110096774 A CN 110096774A
Authority
CN
China
Prior art keywords
supercapacitor
model
electrolyte
energy storage
ion
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
Application number
CN201910319320.6A
Other languages
Chinese (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201910319320.6A priority Critical patent/CN110096774A/en
Publication of CN110096774A publication Critical patent/CN110096774A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

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)
  • Electric Double-Layer Capacitors Or The Like (AREA)

Abstract

本发明公开了一种力学冲击下超级电容器电解液离子重分布效应的仿真方法,该方法包括以下步骤:建立超级电容器的储能动力学模型;在储能动力学模型中,嵌入力学冲击下电解液离子的迁移重分布模型,以构成基于超级电容器的储能‑冲击敏感耦合模型体系;利用有限元软件对储能‑冲击敏感耦合模型体系进行仿真计算,以实现力学冲击下超级电容器电解液离子重分布效应的有效仿真。该方法可以对冲击过程中的电解液离子浓度变化、超级电容器输出电压变化进行数值计算,还可以电容器的重要工作参数进行扫描,分析其对于离子重分布效应显著性的影响规律,可有效促进对超级电容器力学敏感效应的抑制或增强,满足其在不同工作环境的应用需求。

The invention discloses a method for simulating the ion redistribution effect of a supercapacitor electrolyte under mechanical impact. The method comprises the following steps: establishing an energy storage kinetic model of a supercapacitor; The migration and redistribution model of liquid ions is used to form a supercapacitor-based energy storage-shock sensitive coupling model system; the finite element software is used to simulate and calculate the energy storage-shock sensitive coupling model system to realize the supercapacitor electrolyte ion under mechanical impact. Efficient simulation of redistribution effects. This method can numerically calculate the change of the ion concentration of the electrolyte and the output voltage of the supercapacitor during the impact process, and can also scan the important working parameters of the capacitor to analyze its influence on the significance of the ion redistribution effect, which can effectively promote the study of The suppression or enhancement of the mechanical sensitivity effect of the supercapacitor meets its application requirements in different working environments.

Description

力学冲击下超级电容器电解液离子重分布效应的仿真方法Simulation Method of Electrolyte Ion Redistribution Effect in Supercapacitor Under Mechanical Shock

技术领域technical field

本发明涉及超级电容器技术领域,特别涉及一种力学冲击下超级电容器电解液离子重分布效应的仿真方法。The invention relates to the technical field of supercapacitors, in particular to a method for simulating the redistribution effect of electrolyte ions in a supercapacitor under mechanical impact.

背景技术Background technique

超级电容器是一种应用广泛的电化学储能器件,它通过双电层效应和法拉第反应两种机理实现能量的存储和释放,具有容量密度大,高低温性能好、循环寿命长等优势,已经广泛应用于电动汽车、工业检测和可穿戴电子器件等领域。Supercapacitor is a widely used electrochemical energy storage device. It realizes energy storage and release through two mechanisms of electric double layer effect and Faraday reaction. It has the advantages of large capacity density, good high and low temperature performance, and long cycle life. It is widely used in fields such as electric vehicles, industrial detection and wearable electronic devices.

建模仿真方法对于超级电容器的设计和开发具有重要作用,可以通过参数扫描进行数值计算,实现对各主要设计参数的优化设计,部分替代传统设计方法中所需的繁琐实验测试。但是,相关技术中的超级电容器建模仿真方法,均仅考虑了其在常温常压状态下电化学特性,而忽略了温度、压力等环境因素的干扰效应,限制了这类仿真方法的实际应用。Modeling and simulation methods play an important role in the design and development of supercapacitors. Numerical calculations can be performed through parameter sweeps to optimize the design of each main design parameter, and partially replace the tedious experimental tests required in traditional design methods. However, the supercapacitor modeling and simulation methods in the related art only consider their electrochemical characteristics at normal temperature and pressure, and ignore the interference effects of environmental factors such as temperature and pressure, which limits the practical application of such simulation methods .

事实上,随着超级电容器在一些工业检测领域的逐步应用,其在工作状态下不可避免会受到外部力学冲击作用。而超级电容器内部储存有液态电解质,在冲击作用下具有离子重分布效应。因此,迫切的需要一种新型仿真方法,实现对力学冲击下超级电容器电解液离子重分布效应的有效仿真。In fact, with the gradual application of supercapacitors in some industrial detection fields, it is inevitable that they will be subjected to external mechanical impacts under working conditions. The liquid electrolyte is stored inside the supercapacitor, which has an ion redistribution effect under the impact. Therefore, there is an urgent need for a new simulation method to realize the effective simulation of the ion redistribution effect of the supercapacitor electrolyte under mechanical impact.

发明内容Contents of the invention

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。The present invention aims to solve one of the technical problems in the related art at least to a certain extent.

为此,本发明的目的在于提出一种力学冲击下超级电容器电解液离子重分布效应的仿真方法,该方法通过comsol仿真软件平台具体实现,并与matlab仿真平台兼容,具有便于移植、拓展的可行性。For this reason, the purpose of the present invention is to propose a simulation method of the ion redistribution effect of the supercapacitor electrolyte under mechanical impact, the method is specifically realized by the comsol simulation software platform, and is compatible with the matlab simulation platform, and has the feasibility of being easy to transplant and expand sex.

为达到上述目的,本发明提出了超级电容器电解液离子重分布效应的仿真方法,包括以下步骤:建立超级电容器的储能动力学模型;在所述储能动力学模型中,嵌入力学冲击下电解液离子的迁移重分布模型,以构成基于超级电容器的储能-冲击敏感耦合模型体系;利用有限元软件对所述储能-冲击敏感耦合模型体系进行仿真计算,以实现力学冲击下超级电容器电解液离子重分布效应的有效仿真。In order to achieve the above object, the present invention proposes a simulation method for the ion redistribution effect of the supercapacitor electrolyte, comprising the following steps: establishing an energy storage kinetic model of the supercapacitor; The migration and redistribution model of liquid ions to form a supercapacitor-based energy storage-shock sensitive coupling model system; use finite element software to perform simulation calculations on the energy storage-shock sensitive coupling model system to realize supercapacitor electrolysis under mechanical impact Efficient simulation of liquid ion redistribution effects.

本发明实施例的力学冲击下超级电容器电解液离子重分布效应的仿真方法,描述了超级电容器在充放电过程中电解液离子浓度不均衡现象的产生机制以及力学冲击作用产生的离子浓度重分布机制,并考虑了离子浓度重分布效应对超级电容器电化学反应的影响作用,通过COMSOL仿真软件平台具体实现,并与matlab仿真平台兼容,具有便于移植、拓展的可行性;进而可以实现对力学冲击过程中超级电容器电解液离子重分布效应的动力学仿真,包括对电解液离子浓度、器件输出电压的精确仿真计算;可以对放电电流等器件工作参数进行扫描,分析其对于超级电容器力学冲击敏感特性的影响规律;可以对冲击时刻的放电程度参数进行扫描,分析其对于超级电容器力学冲击敏感特性的影响规律。The simulation method of the ion redistribution effect of the supercapacitor electrolyte under the mechanical impact of the embodiment of the present invention describes the generation mechanism of the unbalanced ion concentration of the electrolyte in the charging and discharging process of the supercapacitor and the ion concentration redistribution mechanism caused by the mechanical impact , and considering the influence of the ion concentration redistribution effect on the electrochemical reaction of the supercapacitor, it is realized through the COMSOL simulation software platform, and is compatible with the matlab simulation platform, which is easy to transplant and expand; and then can realize the mechanical impact process Dynamic simulation of electrolyte ion redistribution effect in supercapacitors, including accurate simulation calculation of electrolyte ion concentration and device output voltage; can scan device operating parameters such as discharge current, and analyze its influence on supercapacitor mechanical shock sensitive characteristics The law of influence; the parameters of the discharge degree at the time of impact can be scanned to analyze its influence law on the mechanical shock sensitive characteristics of the supercapacitor.

另外,根据本发明上述实施例的力学冲击下超级电容器电解液离子重分布效应的仿真方法还可以具有以下附加的技术特征:In addition, the simulation method for the ion redistribution effect of the supercapacitor electrolyte under mechanical impact according to the above-mentioned embodiments of the present invention may also have the following additional technical features:

进一步地,在本发明的一个实施例中,所述储能动力学模型包括双电层效应模型和法拉第反应模型。Further, in an embodiment of the present invention, the energy storage kinetic model includes an electric double layer effect model and a Faradaic reaction model.

进一步地,在本发明的一个实施例中,所述储能-冲击敏感耦合模型体系包括双电层效应的电极动力学模型,法拉第反应的电极动力学模型,电解液中的离子扩散模型和力学冲击下的电解液离子重分布模型。Further, in one embodiment of the present invention, the energy storage-shock sensitive coupling model system includes the electrode kinetic model of the electric double layer effect, the electrode kinetic model of the Faraday reaction, the ion diffusion model in the electrolyte, and the mechanical Electrolyte ion redistribution model under shock.

进一步地,在本发明的一个实施例中,所述储能动力学模型通过电极活性材料和电解液材料之间发生的高度可逆的氧化还原反应,使得电荷发生转移产生并存储能量,其中,所述法拉第反应模型储能占主导地位。Further, in one embodiment of the present invention, the energy storage kinetic model generates and stores energy through the highly reversible redox reaction between the electrode active material and the electrolyte material, so that the charge transfer occurs, wherein, The Faradaic reaction model described above dominates energy storage.

进一步地,在本发明的一个实施例中,通过所述储能动力学模型对双电层效应模型进行法拉第电流修正,使得所述超级电容器在充电和放电过程中的电势场、离子浓度场在法拉第反应模型和双电层效应模型两种机制共同驱动下的变化趋势得以确定。Further, in one embodiment of the present invention, the electric double layer effect model is corrected by Faraday current through the energy storage dynamics model, so that the electric potential field and ion concentration field of the supercapacitor during charging and discharging are at The change trend driven by the Faraday reaction model and the electric double layer effect model was determined.

进一步地,在本发明的一个实施例中,在嵌入力学冲击瞬间,电解液将受惯性力驱动产生定向流动,其流场变化的动力学过程由Navier-Stokes方程进行描述:Further, in one embodiment of the present invention, at the moment of embedded mechanical impact, the electrolyte will be driven by inertial force to generate directional flow, and the dynamic process of its flow field change is described by the Navier-Stokes equation:

其中,u1,ρ1,μ分别表示电解液流体的流场速度,密度,动力学粘度。Among them, u 1 , ρ 1 , and μ represent the flow field velocity, density, and dynamic viscosity of the electrolyte fluid, respectively.

进一步地,在本发明的一个实施例中,通过COMSOL弱形式方程和通用形式偏微分方程的复杂模型构建方法,对锂离子电池模块内嵌模型进行修正,使得所述储能动力学模型符合构建的所述储能-冲击敏感耦合模型体系。Further, in one embodiment of the present invention, the built-in model of the lithium-ion battery module is corrected through the complex model construction method of the COMSOL weak form equation and the general form partial differential equation, so that the energy storage kinetic model conforms to the construction The energy storage-shock sensitive coupling model system of .

进一步地,在本发明的一个实施例中,所述通用形式偏微分方程架构如下:Further, in one embodiment of the present invention, the general-form partial differential equation framework is as follows:

其中,u是偏微分方程因变量,t是自变量,其他参数为通用系数。Among them, u is the dependent variable of the partial differential equation, t is the independent variable, and other parameters are general coefficients.

进一步地,在本发明的一个实施例中,利用COMSOL软件对超级电容器冲击敏感特性进行仿真,所述超级电容器在不同的工作模式下出现两种的冲击敏感现象,其中,第一种情况为当所述超级电容器进行大电流充电时,反应电流得失离子的速率远远大于溶液中离子扩散的速率,第二种情况与所述第一种情况相反,由于所述超级电容器的自放电效应,离子浓度分布已变均衡。Further, in one embodiment of the present invention, COMSOL software is used to simulate the shock sensitivity characteristics of the supercapacitor, and the supercapacitor has two kinds of shock sensitivity phenomena in different operating modes, wherein the first case is when When the supercapacitor is charged with a large current, the rate at which the reaction current gains and loses ions is far greater than the rate at which the ions diffuse in the solution. The second case is contrary to the first case. Due to the self-discharge effect of the supercapacitor, ions The concentration distribution has become balanced.

本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.

附图说明Description of drawings

本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present invention will become apparent and easy to understand from the following description of the embodiments in conjunction with the accompanying drawings, wherein:

图1为根据本发明实施例的力学冲击下超级电容器电解液离子重分布效应的仿真方法流程图;Fig. 1 is the simulation method flow chart of supercapacitor electrolyte ion redistribution effect under the mechanical impact according to the embodiment of the present invention;

图2为根据本发明实施例中超级电容器结构示意图;Fig. 2 is a schematic structural diagram of a supercapacitor according to an embodiment of the present invention;

图3为根据本发明实施例中超级电容器双电层效应原理图;Fig. 3 is a schematic diagram of the electric double layer effect of a supercapacitor according to an embodiment of the present invention;

图4为根据本发明实施例中超级电容器的法拉第反应原理图;Fig. 4 is according to the Faraday reaction schematic diagram of supercapacitor in the embodiment of the present invention;

图5为根据本发明实施例中超级电容器在冲击作用下的离子重分布原理图;5 is a schematic diagram of ion redistribution of a supercapacitor under impact according to an embodiment of the present invention;

图6为根据本发明实施例中COMSOL软件的弱形式方程入口示意图;Fig. 6 is according to the weak form equation entry schematic diagram of COMSOL software in the embodiment of the present invention;

图7为根据本发明实施例中超级电容器的两种冲击敏感现象仿真图,其中,(a)为第一种敏感现象:电压向下波动,(b)为第二种敏感现象:电压向上波动;Figure 7 is a simulation diagram of two impact sensitive phenomena of a supercapacitor according to an embodiment of the present invention, wherein (a) is the first sensitive phenomenon: voltage fluctuates downward, and (b) is the second sensitive phenomenon: voltage fluctuates upward ;

图8为根据本发明实施例中超级电容器冲击过程中的离子分布仿真图,其中,(a)为第一种敏感现象,(b)为第二种敏感现象,图中ts代表冲击的持续时间;Figure 8 is a simulation diagram of ion distribution during the shock process of a supercapacitor according to an embodiment of the present invention, wherein (a) is the first sensitive phenomenon, (b) is the second sensitive phenomenon, and ts represents the duration of the shock in the figure ;

图9为根据本发明实施例中不同放电电流和放电程度下的超级电容器电压峰幅值仿真图;FIG. 9 is a simulation diagram of a supercapacitor voltage peak amplitude under different discharge currents and discharge degrees according to an embodiment of the present invention;

图10为根据本发明实施例中不同放电电流下超级电容器的离子浓度场仿真图,其中, (a)为放电电流5mA,(b)为放电电流20mA;Fig. 10 is the ion concentration field simulation diagram of the supercapacitor under different discharge currents according to the embodiment of the present invention, wherein, (a) is discharge current 5mA, (b) is discharge current 20mA;

图11为根据本发明实施例中超级电容器的实测冲击敏感现象变化图,其中,(a)为电压向下波动的冲击敏感现象,(b)为电压向上波动的冲击敏感现象。Fig. 11 is a graph showing changes in the measured impact sensitivity phenomenon of a supercapacitor according to an embodiment of the present invention, wherein (a) is the impact sensitivity phenomenon of downward voltage fluctuation, and (b) is the impact sensitivity phenomenon of upward voltage fluctuation.

具体实施方式Detailed ways

下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。Embodiments of the present invention are described in detail below, examples of which are shown in the drawings, wherein the same or similar reference numerals designate the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the figures are exemplary and are intended to explain the present invention and should not be construed as limiting the present invention.

下面参照附图描述根据本发明实施例提出的超级电容器电解液离子重分布效应的仿真方法及系统,首先将参照附图描述根据本发明实施例提出的超级电容器电解液离子重分布效应的仿真方法。The simulation method and system for the ion redistribution effect of the supercapacitor electrolyte proposed according to the embodiments of the present invention will be described below with reference to the accompanying drawings. First, the simulation method for the ion redistribution effect of the supercapacitor electrolyte proposed according to the embodiments of the present invention will be described with reference to the accompanying drawings .

图1是本发明一个实施例的力学冲击下超级电容器电解液离子重分布效应的仿真方法流程图。Fig. 1 is a flow chart of a simulation method for the ion redistribution effect of the supercapacitor electrolyte under mechanical impact according to an embodiment of the present invention.

如图1所示,该力学冲击下超级电容器电解液离子重分布效应的仿真方法包括以下步骤:As shown in Figure 1, the simulation method for the ion redistribution effect of the supercapacitor electrolyte under mechanical impact includes the following steps:

在步骤S101中,建立超级电容器的储能动力学模型。In step S101, an energy storage dynamics model of the supercapacitor is established.

其中,储能动力学模型包括双电层效应模型和法拉第反应模型。Among them, the energy storage kinetic model includes the electric double layer effect model and the Faraday reaction model.

具体而言,储能动力学模型通过电极活性材料和电解液材料之间发生的高度可逆的氧化还原反应,使得电荷发生转移产生并存储能量,其中,法拉第反应模型储能占主导地位。Specifically, the energy storage kinetic model generates and stores energy through the highly reversible redox reaction between the electrode active material and the electrolyte material, in which the Faraday reaction model is dominant.

进一步地,在本发明的一个实施例中,通过储能动力学模型对双电层效应模型进行法拉第电流修正,使得超级电容器在充电和放电过程中的电势场、离子浓度场在法拉第反应模型和双电层效应模型两种机制共同驱动下的变化趋势得以确定。Further, in one embodiment of the present invention, the electric double layer effect model is corrected by faraday current through the energy storage dynamics model, so that the electric potential field and ion concentration field of the supercapacitor in the charging and discharging process are in the Faraday reaction model and The change trend driven by the two mechanisms of the electric double layer effect model was determined.

具体地,如图2所示,超级电容器的正负极均为具有高比表面积的多孔薄膜材料,中间由隔膜隔开,电极和隔膜的孔隙中充满液态电解液。超级电容器的储能机理包括双电层效应和法拉第赝电容效应两种机制,本发明实施例将根据两种机制分别进行机理建模研究。Specifically, as shown in Figure 2, the positive and negative electrodes of the supercapacitor are both porous film materials with high specific surface area, separated by a separator in the middle, and the pores of the electrodes and the separator are filled with liquid electrolyte. The energy storage mechanism of a supercapacitor includes two mechanisms, the electric double layer effect and the Faraday pseudocapacitance effect, and the embodiments of the present invention will carry out mechanism modeling studies based on the two mechanisms.

如图3所示,多孔电极具有微观孔隙结构,与电解液之间形成了丰富的固-液界面。超级电容器的双电层效应储能机制正是基于固-液界面处正负离子吸附所形成的双电层结构。As shown in Figure 3, the porous electrode has a microscopic pore structure, forming a rich solid-liquid interface with the electrolyte. The energy storage mechanism of the electric double layer effect of supercapacitors is based on the electric double layer structure formed by the adsorption of positive and negative ions at the solid-liquid interface.

超级电容器的充放电过程的本质是微观双电层结构的形成与解体,以及由此产生的电极和电解液中电势、离子浓度的耦合变化。针对超级电容器的电极和隔膜两个部分,本发明实施例分别构建电场和电解液浓度场模型,进行动力学建模研究。The essence of the charging and discharging process of a supercapacitor is the formation and disintegration of the microscopic electric double layer structure, and the resulting coupled changes in the potential and ion concentration of the electrode and electrolyte. For the two parts of the electrode and the diaphragm of the supercapacitor, the embodiment of the present invention respectively constructs electric field and electrolyte concentration field models, and conducts dynamic modeling research.

首先,在多孔电极中,从电极基体相流向电极电解液相的电流可以由双电层电流描述,First, in a porous electrode, the current flowing from the electrode matrix phase to the electrode electrolyte phase can be described by the electric double layer current,

其中,i1表示电极电解液相的表面电流密度,iDL表示双电层效应所产生的电流密度。Among them, i 1 represents the surface current density of the electrolyte phase of the electrode, and i DL represents the current density generated by the electric double layer effect.

对于电极基体相,电流和电势之间的关系可以由欧姆定律给出,For the electrode matrix phase, the relationship between current and potential can be given by Ohm's law,

其中,is为电极基体相的表面电流密度,σs,eff电极基体相的等效电导率和Ψs为电极基体相的电势。where is is the surface current density of the electrode matrix phase, σs ,eff is the equivalent conductivity of the electrode matrix phase and Ψs is the potential of the electrode matrix phase.

但是,对于电极电解液相,需要对欧姆定律进行浓度场修正,However, for the electrode electrolyte phase, a concentration field correction to Ohm's law is required,

其中,i1为电极电解液相的表面电流密度,σ1,eff为电极电解液相的等效电导率,Ψ1为电极电解液相的电势,t+为电解液溶液的转移系数,c1为电解液溶液的离子浓度,F,R和 T分别表示法拉第常数,通用气体常数和开尔文温度。Among them, i 1 is the surface current density of the electrolyte phase of the electrode, σ 1,eff is the equivalent conductivity of the electrolyte phase of the electrode, Ψ 1 is the potential of the electrolyte phase of the electrode, t + is the transfer coefficient of the electrolyte solution, c 1 is the ion concentration of the electrolyte solution, F, R and T represent Faraday's constant, universal gas constant and Kelvin temperature, respectively.

在孔隙介质中,电解液相的等效电导率σl,eff服从布拉格曼关系,也即σ1,eff=σ1ε1 1.5,其中,σ1为电解液离子导电率,ε1表示电极的孔隙率。In porous media, the equivalent conductivity σ l,eff of the electrolyte phase obeys the Bragman relation, that is, σ 1,eff1 ε 1 1.5 , where σ 1 is the ionic conductivity of the electrolyte, and ε 1 represents the electrode porosity.

电解液相的电势方程可以由方程(1)和(3)联立获得The potential equation of the electrolyte phase can be obtained by combining equations (1) and (3)

由于电荷守恒,始终成立。将此电荷守恒方程与方程(2)联立:Due to the conservation of charge, always holds. Combine this charge conservation equation with equation (2):

离子浓度场同时受到反映电流和离子迁移的影响,如方程(6)所示The ion concentration field is affected by both reflected current and ion migration, as shown in equation (6)

其中,ε1和D1,eff分别表示电极的孔隙率和电解液溶液的离子扩散系数,t表示时间。Among them, ε 1 and D 1,eff represent the porosity of the electrode and the ion diffusion coefficient of the electrolyte solution, respectively, and t represents the time.

在隔膜部分,由于隔膜基体不导电,is=0。所以可以得到进而电解液相的电势方程可以表示为:In the diaphragm part, since the diaphragm matrix is not conductive, i s =0. so you can get Then the potential equation of the electrolyte phase can be expressed as:

与方程(6)类似,由于隔膜部分没有双电层反应,电解液相的离子浓度场方程可以由下式给出,Similar to equation (6), since there is no electric double layer reaction in the diaphragm part, the ion concentration field equation of the electrolyte phase can be given by the following formula,

根据双电层效应的机理,双电层电流密度iDL可以由下式计算[85–87] According to the mechanism of the electric double layer effect, the electric double layer current density i DL can be calculated by the following formula [85–87]

其中,av表示多孔电极的比表面积,Cdl表示与双电层效应有关的电容参数。Among them, a v represents the specific surface area of the porous electrode, and C dl represents the capacitance parameter related to the electric double layer effect.

除此此外,部分超级电容器还通过法拉第反应进行能量的储存和释放。法拉第反应是电极活性材料和电解液材料之间发生的高度可逆的氧化还原反应(例如氧化钌电极与硫酸电解液之间的氧化还原反应)。如图4所示,氧化还原反应的得失电子效应使电荷发生转移,从而实现能量的储存。In addition, some supercapacitors also store and release energy through Faraday reactions. Faradaic reactions are highly reversible redox reactions that occur between electrode active materials and electrolyte materials (such as the redox reaction between a ruthenium oxide electrode and a sulfuric acid electrolyte). As shown in Figure 4, the electron gain and loss effect of the redox reaction makes the charge transfer, thereby realizing the energy storage.

法拉第反应储能占主导作用的超级电容器,其电化学体系中的法拉第电流要远大于双电层电流。因此,此类超级电容器储能动力学模型需要在双电层储能模型方程(1)~(9) 的基础上进行一系列修正。The faraday current in the electrochemical system of the supercapacitor in which the faraday reaction plays a leading role in energy storage is much larger than the electric double layer current. Therefore, this type of supercapacitor energy storage kinetic model needs a series of corrections based on the electric double layer energy storage model equations (1)-(9).

首先,在多孔电极中,从电极基体相流向电极电解液相的电流由双电层电流和法拉第电流共同构成,原方程(1)需要加入法拉第电流项:First, in a porous electrode, the current flowing from the electrode matrix phase to the electrode electrolyte phase is composed of the electric double layer current and the Faraday current. The original equation (1) needs to add the Faraday current term:

其中,iF表示法拉第反应所产生的电流密度。Among them, i F represents the current density generated by the Faradaic reaction.

相应的,方程(4)的电解液相电势方程和方程(5)的电极基体相电势方程也需要加入法拉第电流修正项:Correspondingly, the electrolyte phase potential equation of Equation (4) and the electrode matrix phase potential equation of Equation (5) also need to add the Faraday current correction term:

对于方程(6)所描述的离子浓度场,反应电流项应由双电层电流和法拉第电流共同组成,For the ion concentration field described in Equation (6), the reaction current term should be composed of the electric double layer current and the Faraday current,

而根据电化学理论,法拉第电流密度可以由下式计算:According to the electrochemical theory, the Faraday current density can be calculated by the following formula:

iF=avj1oc (14)i F =a v j 1oc (14)

其中,j1oc表示法拉第转移电流密度,其取值可以通过Butler-Volmer方程加以计算:Among them, j 1oc represents the Faraday transfer current density, and its value can be calculated by the Butler-Volmer equation:

其中,i0表示法拉第交换电流密度,αa和αc分别表示阳极和阴极的转移系数,Uoc表示法拉第过程的开路电势。Among them, i 0 represents the Faradaic exchange current density, α a and α c represent the transfer coefficients of the anode and cathode, respectively, and U oc represents the open circuit potential of the Faradaic process.

通过方程(10)~(15)对双电层储能模型进行修正,本发明实施例得到了普适的超级电容器的储能动力学模型。通过这个模型,超级电容器在充电和放电过程中的电势场、离子浓度场在法拉第反应和双电层效应两种机制共同驱动下的变化过程可以被完全确定。By modifying the electric double layer energy storage model through equations (10) to (15), the embodiment of the present invention obtains a universal energy storage dynamics model of supercapacitors. Through this model, the change process of the electric potential field and ion concentration field of the supercapacitor during the charging and discharging process driven by the two mechanisms of Faraday reaction and electric double layer effect can be completely determined.

在步骤S102中,在储能动力学模型中,嵌入力学冲击下电解液离子的迁移重分布模型,以构成基于超级电容器的储能-冲击敏感耦合模型体系。In step S102 , in the kinetic model of energy storage, a migration and redistribution model of electrolyte ions under mechanical impact is embedded to form a supercapacitor-based energy storage-shock sensitive coupling model system.

进一步地,在本发明的一个实施例中,储能-冲击敏感耦合模型体系包括双电层效应的电极动力学模型,法拉第反应的电极动力学模型,电解液中的离子扩散模型和力学冲击下的电解液离子重分布模型。Further, in one embodiment of the present invention, the energy storage-shock sensitive coupling model system includes the electrode dynamics model of the electric double layer effect, the electrode kinetics model of the Faraday reaction, the ion diffusion model in the electrolyte, and the Electrolyte ion redistribution model.

需要说明的是,储能动力学建模是研究微纳能源器件的有效方法,在超级电容器研究领域已有诸多应用。本发明实施例综合考虑超级电容器的储能效应和冲击敏感效应,进行一体化耦合建模。It should be noted that energy storage dynamics modeling is an effective method for studying micro-nano energy devices, and has been applied in many fields in the field of supercapacitor research. The embodiment of the present invention comprehensively considers the energy storage effect and impact sensitivity effect of the supercapacitor, and performs integrated coupling modeling.

具体而言,在超级电容器的电化学体系中,电化学反应的进行将不断地在一个电极产生离子,同时在另一个电极消耗离子,总的趋势是使得电解液中离子的正负电极的浓度分布产生不均衡,并不断积累。而另一方面,离子在电解液中具有扩散效应,总的趋势是使得电解液中的离子浓度均衡化。当超级电容器的充放电电流较大、电化学反应较强烈时,电解液的扩散效应将远不足以抵消反应过程导致的离子快速产生和消耗。因此,电解液中的离子浓度不均衡将可能较为显著,如图5所示。Specifically, in the electrochemical system of a supercapacitor, the electrochemical reaction will continuously generate ions at one electrode and consume ions at the other electrode at the same time. The general trend is to make the concentration of the positive and negative electrodes of ions in the electrolyte Distribution creates imbalances that accumulate. On the other hand, ions have a diffusion effect in the electrolyte, and the general trend is to equalize the ion concentration in the electrolyte. When the charge and discharge current of the supercapacitor is large and the electrochemical reaction is strong, the diffusion effect of the electrolyte will be far from enough to counteract the rapid production and consumption of ions caused by the reaction process. Therefore, the ion concentration imbalance in the electrolyte may be more significant, as shown in FIG. 5 .

下面具体说明充放电过程中的离子浓度差的产生过程及其导致的冲击敏感效应。在充放电开始之前,反应电流为0,电解液中的离子经过长期的扩散效应,处于均衡状态,The generation process of the ion concentration difference in the charging and discharging process and the impact sensitivity effect caused by it will be described in detail below. Before charging and discharging starts, the reaction current is 0, and the ions in the electrolyte are in a balanced state after a long-term diffusion effect.

当充放电过程开始后,超级电容器正负电极的反应电流符号相反,两电极分别得到和失去离子。When the charging and discharging process starts, the signs of the reaction currents of the positive and negative electrodes of the supercapacitor are opposite, and the two electrodes gain and lose ions respectively.

对于充电过程:For the charging process:

对于放电过程:For the discharge process:

因此,当充电或放电过程开始后,正负电极的离子浓度产生相反的变化,离子浓度的均衡状态将被打破,正负极之间将会产生离子浓度差。这种离子浓度差效应体现为离子浓度场具有非零的梯度值:Therefore, when the charge or discharge process starts, the ion concentration of the positive and negative electrodes will change oppositely, the equilibrium state of the ion concentration will be broken, and a difference in ion concentration will be generated between the positive and negative electrodes. This ion concentration difference effect is manifested by the ion concentration field having a non-zero gradient value:

在加速度冲击瞬间,电解液将受惯性力驱动产生定向流动,其流场变化的动力学过程可以由Navier-Stokes方程进行描述:At the moment of acceleration shock, the electrolyte will be driven by inertial force to generate directional flow, and the dynamic process of its flow field change can be described by the Navier-Stokes equation:

其中,u1,ρ1,μ分别表示电解液流体的流场速度,密度,动力学粘度。Among them, u 1 , ρ 1 , and μ represent the flow field velocity, density, and dynamic viscosity of the electrolyte fluid, respectively.

根据流体力学理论可知,在存在离子浓度梯度的情况下,电解液流动将在短时间内显著改变离子浓度场分布,因此需要对方程(13)的电极电解液相离子浓度场动力学模型添加流场修正项,According to the theory of fluid mechanics, in the presence of an ion concentration gradient, the electrolyte flow will significantly change the ion concentration field distribution in a short time, so it is necessary to add flow field correction term,

类似地,对于方程(8)所描述的隔膜部分的离子浓度场方程,也需要进行相应的流场修正,Similarly, for the ion concentration field equation of the diaphragm part described by equation (8), the corresponding flow field correction is also required,

方程(20)~(23)构成了超级电容器电解液流场和离子浓度场加速度冲击响应模型。将其与储能动力学模型相结合,超级电容器在冲击瞬间由电解液离子冲击敏感效应对电化学体系产生的复杂耦合变化过程可以被定量描述。从而对是本发明实施例后续对电解液离子冲击敏感效应进行仿真研究奠定了基础。Equations (20)-(23) constitute the supercapacitor electrolyte flow field and ion concentration field acceleration shock response model. Combining it with the energy storage dynamics model, the complex coupling change process of the supercapacitor at the moment of impact caused by the impact sensitive effect of electrolyte ions on the electrochemical system can be quantitatively described. Thus, a foundation is laid for the subsequent simulation research on the ion impact sensitivity effect of the electrolyte according to the embodiment of the present invention.

在步骤S103中,利用有限元软件对储能-冲击敏感耦合模型体系进行仿真计算,以实现力学冲击下超级电容器电解液离子重分布效应的有效仿真。In step S103, the finite element software is used to simulate the energy storage-shock sensitive coupling model system, so as to realize the effective simulation of the ion redistribution effect of the supercapacitor electrolyte under mechanical shock.

其中,本发明实施例利用的有限元软件为COMSOL软件,COMSOL软件是一款功能强大的多物理场仿真软件,其中的锂离子电池模块的主体模型方程与本发明实施例所提出的超级电容器冲击敏感动力学模型较为接近,因此在此模块的基础上进行修正,实现对超级电容器冲击敏感特性的仿真计算。Wherein, the finite element software that the embodiment of the present invention utilizes is COMSOL software, and COMSOL software is a powerful multi-physics field simulation software, and the main body model equation of the lithium-ion battery module wherein and the supercapacitor impact that the embodiment of the present invention proposes The sensitive dynamic model is relatively close, so it is modified on the basis of this module to realize the simulation calculation of the impact sensitive characteristics of supercapacitors.

锂离子模块的内嵌模型如下:The embedded model of the Li-ion module is as follows:

但是,该内嵌模型仅包含了电极储能动力学方程,没有考虑本发明实施例所提出的超级电容器的离子冲击敏感效应。However, the embedded model only includes the electrode energy storage kinetic equation, and does not consider the ion impact sensitivity effect of the supercapacitor proposed by the embodiment of the present invention.

因此,本发明实施例通过COMSOL弱形式方程和通用形式偏微分方程的复杂模型构建方法,对锂离子电池模块内嵌模型进行修正,使得储能动力学模型符合构建建的电化学储能-冲击敏感耦合模型体系。Therefore, the embodiment of the present invention corrects the embedded model of the lithium-ion battery module through the complex model construction method of COMSOL weak form equation and general form partial differential equation, so that the energy storage kinetic model conforms to the established electrochemical energy storage-shock Sensitive coupling model system.

一方面,内嵌模型忽略了电解液的流动及流场变化对离子浓度场、电势场的影响,可以通过对内嵌方程的弱形式代码进行编辑,实现对这一缺陷的修正。所谓弱形式,是有限元计算中将原始偏微分方程组进行变分处理后的一种形式。如图6所示,在COMSOL软件系统设定中选择“方程视图”选项,即可进入弱形式编辑模式。将方程(22)和方程(23) 所描述的离子浓度场修正模型进行弱形式转换后,替代模块中的内嵌的原始代码,就实现了对模块浓度场方程的流场修正。On the one hand, the embedded model ignores the influence of electrolyte flow and flow field changes on the ion concentration field and electric potential field. This defect can be corrected by editing the weak form code of the embedded equation. The so-called weak form is a form after the original partial differential equations are subjected to variational processing in the finite element calculation. As shown in Figure 6, select the "equation view" option in the COMSOL software system settings to enter the weak form editing mode. After transforming the ion concentration field correction model described by Equation (22) and Equation (23) into a weak form, the flow field correction of the concentration field equation of the module is realized by replacing the embedded original code in the module.

另一方面,本发明实施例的模型体系中冲击加速度与电解液流场的关系由方程(20) 和方程(21)描述,可以通过在模块中添加新方程的方式实现。COMSOL提供的通用形式偏微分方程架构如下:On the other hand, the relationship between shock acceleration and electrolyte flow field in the model system of the embodiment of the present invention is described by Equation (20) and Equation (21), which can be realized by adding new equations to the module. The general form partial differential equation framework provided by COMSOL is as follows:

对照如方程(20)和方程(21)所示的模型方程,将各项系数对应输入通用形式方程(25)中,即可完成仿真模块的补充修正。Comparing with the model equations shown in Equation (20) and Equation (21), the corresponding coefficients are input into the general form Equation (25), and the supplementary correction of the simulation module can be completed.

故可以利用COMSOL软件对超级电容器冲击敏感特性进行仿真。仿真结果表明,超级电容器在不同的工作模式下会出现两种典型的冲击敏感现象,分别如图7所示。Therefore, COMSOL software can be used to simulate the impact sensitive characteristics of supercapacitors. The simulation results show that two typical shock-sensitive phenomena will appear in different working modes of the supercapacitor, as shown in Fig. 7 respectively.

如果超级电容器在大电流充电完成后立即开始放电,并在开始放电后不久后受到加速度冲击,将发生第一种敏感现象,如图7(a)所示。此时输出电压在冲击瞬间向下波动,记为ΔV<0。ΔV表示冲击过程中输出电压的上升值。If the supercapacitor starts to discharge immediately after the high-current charging is completed, and receives an acceleration shock shortly after starting to discharge, the first sensitive phenomenon will occur, as shown in Fig. 7(a). At this time, the output voltage fluctuates downward at the moment of impact, which is recorded as ΔV<0. ΔV represents the rising value of the output voltage during the impact.

如果超级电容器在充满电后进行了长时间的自放电过程,在自放电完成以后开始进行放电,放电过程中发生冲击,将发生第二种敏感现象,如图7(b)所示。与图7(a)相反,此时输出电压在冲击瞬间向上波动,记为ΔV>0。If the supercapacitor undergoes a long self-discharge process after being fully charged, and starts to discharge after the self-discharge is completed, and an impact occurs during the discharge process, the second sensitive phenomenon will occur, as shown in Figure 7(b). Contrary to Figure 7(a), the output voltage fluctuates upward at the moment of impact, which is recorded as ΔV>0.

对冲击过程中离子浓度场的仿真可以更好地解释这两种敏感现象的本质。The simulation of the ion concentration field during the impact can better explain the nature of these two sensitive phenomena.

对于第一种情况,当超级电容器进行大电流充电时,反应电流得失离子的速率远远大于溶液中离子扩散的速率。因此,在正极,离子大量堆积,在负极,离子被大量消耗。在冲击开始前,正极的离子浓度远大于负极的离子浓度,其离子浓度场分布如图8(a)中的方形标志线所示。For the first case, when the supercapacitor is charged with a large current, the rate of ion gain and loss by the reaction current is much greater than the rate of ion diffusion in the solution. Therefore, at the positive electrode, ions are accumulated in large quantities, and at the negative electrode, ions are consumed in large quantities. Before the impact begins, the ion concentration of the positive electrode is much greater than that of the negative electrode, and its ion concentration field distribution is shown in the square mark line in Fig. 8(a).

当冲击发生后,电解液溶液的流动迅速的改变了离子浓度场分布,在极短时间内溶液中的离子浓度分布恢复到近似均衡状态。这个离子再分布过程如图8(a)所示,ts代表冲击的持续时间。根据超级电容器的冲击敏感模型,离子浓度场c1的变化将直接影响电解液电势场Φ1,进而影响其输出电压。When the impact occurs, the flow of the electrolyte solution rapidly changes the ion concentration field distribution, and the ion concentration distribution in the solution returns to an approximately balanced state in a very short time. This ion redistribution process is shown in Fig. 8(a), where t s represents the duration of the shock. According to the impact sensitive model of the supercapacitor, the change of the ion concentration field c 1 will directly affect the electrolyte potential field Φ 1 , and then affect its output voltage.

进一步地,离子浓度场变化产生的电压波动现象可以由如下的电化学原理加以解释。在放电过程中,离子不断与正极发生结合,冲击使得离子浓度分布趋向均衡,也即减少了正极的离子。正极离子的减少使得超级电容器正极放电容量降低,从而使其输出电压下降。Furthermore, the voltage fluctuation phenomenon caused by the change of ion concentration field can be explained by the following electrochemical principle. During the discharge process, the ions are continuously combined with the positive electrode, and the impact makes the ion concentration distribution tend to be balanced, that is, the ion of the positive electrode is reduced. The reduction of the positive electrode ions reduces the discharge capacity of the positive electrode of the supercapacitor, thereby reducing its output voltage.

第二种敏感现象与之相反,由于超级电容器自放电效应的作用,在经历了长时间的自放电过程以后,其离子浓度分布已经变得均衡。所以,当电容器开始大电流放电后,负极的离子浓度开始显著的超过正极,如图8(b)中的方形标志线所示。在冲击过程中,离子浓度场同样在极短时间内恢复近似均衡。而根据方程(11)对电解液相电势场的描述,离子浓度场c1对电势场的影响取决于由于图8(a)和图8(b)中的取值的符号相反,所以第二种现象对应的电压波动也与第一种现象相反,表现为冲击过程中电压向上波动。基于电解液离子再分布效应,冲击使得电解液离子浓度场分布更加均衡,也即增加了正极的离子。因此,电解液可以为正极的法拉第反应提供更多的离子,增大了超级电容器正极的放电容量,使输出电压上升。The second sensitive phenomenon is the opposite. Due to the self-discharge effect of the supercapacitor, after a long period of self-discharge, the ion concentration distribution has become balanced. Therefore, when the capacitor starts to discharge with a large current, the ion concentration of the negative electrode begins to significantly exceed that of the positive electrode, as shown by the square mark line in Figure 8(b). During the impact process, the ion concentration field also restores approximately equilibrium in a very short time. According to the description of the potential field of the electrolyte phase by equation (11), the influence of the ion concentration field c 1 on the potential field depends on Since in Figure 8(a) and Figure 8(b) The signs of the values are opposite, so the voltage fluctuation corresponding to the second phenomenon is also opposite to the first phenomenon, which is manifested as an upward fluctuation of the voltage during the impact process. Based on the redistribution effect of electrolyte ions, the impact makes the concentration field distribution of electrolyte ions more balanced, that is, increases the positive ions. Therefore, the electrolyte can provide more ions for the Faraday reaction of the positive electrode, which increases the discharge capacity of the positive electrode of the supercapacitor and increases the output voltage.

总体而言,超级电容器电压波动的情况取决于冲击发生时的离子浓度分布。如图7(a) 和图8(a)所示,如果冲击初始时刻离子集中在正极(c1|positive>c1|negative),则冲击过程中电压将向下波动。如图7(b)和图8(b)所示,相反地,如果冲击初始时刻离子集中在负极(c1|positive<c1|negative),则过程中电压将向上波动。Overall, the voltage fluctuation of the supercapacitor depends on the ion concentration distribution when the shock occurs. As shown in Fig. 7(a) and Fig. 8(a), if the ions are concentrated at the positive electrode (c 1 |positive>c 1 |negative) at the beginning of the shock, the voltage will fluctuate downward during the shock. As shown in Fig. 7(b) and Fig. 8(b), on the contrary, if the ions are concentrated on the negative electrode (c 1 |positive<c 1 |negative) at the initial moment of impact, the voltage will fluctuate upward during the process.

放电电流和放电程度是影响超级电容器离子浓度场分布的两个重要因素(定义放电程度为1减去输出电压与最大充电电压的比值),本发明实施例提出的仿真方法可以分析参数对于力学冲击下离子重分布效应的影响。如图9所示,其展示了不同放电电流和放电程度下,冲击产生的电压向上波动的幅值,由此可以分析超级电容器的冲击响应特性。Discharge current and discharge degree are two important factors affecting the distribution of supercapacitor ion concentration field (the definition of discharge degree is 1 minus the ratio of output voltage and maximum charging voltage), and the simulation method proposed in the embodiment of the present invention can analyze parameters for mechanical impact Under the influence of ion redistribution effect. As shown in Figure 9, it shows the amplitude of the upward fluctuation of the voltage generated by the impact under different discharge currents and discharge degrees, so that the impact response characteristics of the supercapacitor can be analyzed.

首先,当放电电流相同时,超级电容器的放电程度越大,冲击产生的电压向上波动的幅值就越大;其次,当超级电容器的放电程度相同时,放电电流越大,冲击产生的电压向上波动的幅值就越大;第三,当放电电流较小时,电压波动幅值与放电程度的关系如图9 中的方形标志线和圆形标志线所示。此时,当放电程度较小时,电压波动幅值与放电程度为近似线性关系,而当放电程度较大时,随着放电程度继续增大,电压波动幅值的增加越来越缓慢,最终趋于收敛。第四,当放电电流较大时,电压波动幅值与放电程度的关系能够始终保持近似线性,如图9中的三角形标志线所示。First, when the discharge current is the same, the greater the discharge degree of the supercapacitor, the greater the amplitude of the upward fluctuation of the voltage generated by the shock; secondly, when the discharge degree of the supercapacitor is the same, the greater the discharge current, the greater the upward fluctuation of the voltage generated by the shock. The larger the fluctuation amplitude is; third, when the discharge current is small, the relationship between the voltage fluctuation amplitude and the discharge degree is shown as the square mark line and the circle mark line in Figure 9. At this time, when the discharge degree is small, the voltage fluctuation amplitude has an approximately linear relationship with the discharge degree, and when the discharge degree is large, as the discharge degree continues to increase, the voltage fluctuation amplitude increases more and more slowly, and eventually tends to in convergence. Fourth, when the discharge current is large, the relationship between the voltage fluctuation amplitude and the discharge degree can always be approximately linear, as shown by the triangle mark line in FIG. 9 .

图9中超级电容器冲击电压波动幅值与放电电流和放电程度的关系可以由图10的离子浓度场仿真结果进一步加以解释。The relationship between the shock voltage fluctuation amplitude of the supercapacitor in Fig. 9 and the discharge current and discharge degree can be further explained by the ion concentration field simulation results in Fig. 10 .

首先,离子浓度梯度的产生来源于法拉第反应电流得失离子效应的积累。所以,放电程度越大,负极积累的离子就越多,离子浓度梯度就越大。因此,冲击产生的离子浓度均衡化再分布效应对离子浓度场的改变就越大,进而也将对输出电压产生更加显著的影响。First, the generation of the ion concentration gradient comes from the accumulation of the Faraday reaction current gain and loss ion effect. Therefore, the greater the degree of discharge, the more ions accumulated in the negative electrode, and the greater the ion concentration gradient. Therefore, the ion concentration equalization redistribution effect generated by the impact will change the ion concentration field more, and then will have a more significant impact on the output voltage.

第二,离子浓度梯度产生的本质是法拉第反应得失离子的速率超过了离子扩散的速率,因此放电电流越大,法拉第反应越剧烈,离子浓度梯度也将越大。Second, the essence of the ion concentration gradient is that the rate of Faraday reaction gain and loss of ions exceeds the rate of ion diffusion, so the greater the discharge current, the more intense the Faraday reaction and the greater the ion concentration gradient.

第三,根据方程(13)所描述的离子浓度场动力学机制,离子扩散的速率取决于离子浓度梯度的大小,因此随着放电程度增大,离子浓度梯度不断积累。当离子浓度梯度积累到足够大时,离子扩散效应将足以抵消法拉第反应对离子浓度场的影响,进而离子浓度场保持在稳定的状态,如图10所示。所以当放电程度足够大以后,尽管放电程度继续增大,但冲击产生的电压波动幅值基本趋于稳定。Third, according to the ion concentration field dynamics mechanism described by Equation (13), the rate of ion diffusion depends on the size of the ion concentration gradient, so as the degree of discharge increases, the ion concentration gradient continues to accumulate. When the ion concentration gradient accumulates to a large enough size, the ion diffusion effect will be sufficient to offset the influence of the Faradaic reaction on the ion concentration field, and then the ion concentration field will remain in a stable state, as shown in Figure 10. Therefore, when the discharge degree is large enough, although the discharge degree continues to increase, the amplitude of the voltage fluctuation caused by the impact basically tends to be stable.

因此,本发明实施例的仿真方法已经可以较为全面的分析超级电容器在力学冲击下的离子浓度重分布效应。为了进一步验证仿真分析结论的正确性,本发明实施例利用马歇特锤击实验系统对超级电容器进行冲击响应特性测试,并利用CHI660电化学工作站记录器件在冲击瞬间的电压响应信号。Therefore, the simulation method of the embodiment of the present invention has been able to comprehensively analyze the ion concentration redistribution effect of the supercapacitor under mechanical impact. In order to further verify the correctness of the conclusion of the simulation analysis, the embodiment of the present invention uses the Marshall hammer test system to test the impact response characteristics of the supercapacitor, and uses the CHI660 electrochemical workstation to record the voltage response signal of the device at the moment of impact.

下面的实验测试证实了仿真分析中论述的两种冲击敏感现象。The following experimental tests confirm the two shock-sensitivity phenomena discussed in the simulation analysis.

对应于第一种冲击敏感现象,实验测试以20mA电流对超级电容器进行充电,在充电完成后立即开始放电并进行马歇特锤击测试。实验记录的器件输出电压如图11(a)所示,在冲击瞬间出现了电压向下波动(ΔV<0)的现象。Corresponding to the first impact sensitive phenomenon, the experimental test charges the supercapacitor with a current of 20mA, and starts discharging immediately after the charging is completed, and performs a Marshall impact test. The output voltage of the device recorded in the experiment is shown in Figure 11(a), and the voltage fluctuates downward (ΔV<0) at the moment of impact.

对应于第二种冲击敏感现象,超级电容器充满电压后对其进行长时间恒压保持以使得其离子浓度场恢复均衡,随后开始以20mA电流放电,在放电过程中进行马歇特锤击测试。实验记录的器件输出电压如图11(b)所示,在冲击瞬间出现了电压向上波动(ΔV>0)的现象。Corresponding to the second impact sensitive phenomenon, after the supercapacitor is fully charged, it is kept at a constant voltage for a long time to restore its ion concentration field to balance, and then it starts to discharge at a current of 20mA, and a Marshall impact test is performed during the discharge process. The output voltage of the device recorded in the experiment is shown in Figure 11(b), and the voltage fluctuates upward (ΔV>0) at the moment of impact.

因此,根据本发明实施例提出的力学冲击下超级电容器电解液离子重分布效应的仿真方法,可实现对超级电容器在力学冲击下的离子重分布现象的有效仿真,仿真结果在定性上与实验结果一致,可有效促进对超级电容器力学冲击特性的研究,为应用于特殊力学环境的超级电容器研发提供技术支持。Therefore, according to the simulation method of the supercapacitor electrolyte ion redistribution effect under the mechanical impact proposed by the embodiment of the present invention, the effective simulation of the ion redistribution phenomenon of the supercapacitor under the mechanical impact can be realized, and the simulation results are qualitatively consistent with the experimental results. Consistent, can effectively promote the research on the mechanical impact characteristics of supercapacitors, and provide technical support for the research and development of supercapacitors used in special mechanical environments.

此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。In addition, the terms "first" and "second" are used for descriptive purposes only, and cannot be interpreted as indicating or implying relative importance or implicitly specifying the quantity of indicated technical features. Thus, the features defined as "first" and "second" may explicitly or implicitly include at least one of these features. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.

在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。In the present invention, unless otherwise clearly specified and limited, terms such as "installation", "connection", "connection" and "fixation" should be understood in a broad sense, for example, it can be a fixed connection or a detachable connection , or integrated; it may be mechanically connected or electrically connected; it may be directly connected or indirectly connected through an intermediary, and it may be the internal communication of two components or the interaction relationship between two components, unless otherwise specified limit. Those of ordinary skill in the art can understand the specific meanings of the above terms in the present invention according to specific situations.

在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。In the present invention, unless otherwise clearly specified and limited, the first feature may be in direct contact with the first feature or the first and second feature may be in direct contact with the second feature through an intermediary. touch. Moreover, "above", "above" and "above" the first feature on the second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature. "Below", "beneath" and "beneath" the first feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is less horizontally than the second feature.

在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。In the description of this specification, descriptions referring to the terms "one embodiment", "some embodiments", "example", "specific examples", or "some examples" mean that specific features described in connection with the embodiment or example , structure, material or characteristic is included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the described specific features, structures, materials or characteristics may be combined in any suitable manner in any one or more embodiments or examples. In addition, those skilled in the art can combine and combine different embodiments or examples and features of different embodiments or examples described in this specification without conflicting with each other.

尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。Although the embodiments of the present invention have been shown and described above, it can be understood that the above embodiments are exemplary and should not be construed as limiting the present invention, those skilled in the art can make the above-mentioned The embodiments are subject to changes, modifications, substitutions and variations.

Claims (9)

1.一种力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,包括以下步骤:1. A simulation method of supercapacitor electrolyte ion redistribution effect under mechanical impact, it is characterized in that, comprises the following steps: 建立超级电容器的储能动力学模型;Establish the energy storage dynamics model of supercapacitor; 在所述储能动力学模型中,嵌入力学冲击下电解液离子的迁移重分布模型,以构成基于超级电容器的储能-冲击敏感耦合模型体系;以及In the energy storage dynamics model, a migration and redistribution model of electrolyte ions under mechanical impact is embedded to form a supercapacitor-based energy storage-shock sensitive coupling model system; and 利用有限元软件对所述储能-冲击敏感耦合模型体系进行仿真计算,以实现力学冲击下超级电容器电解液离子重分布效应的有效仿真。The finite element software is used to simulate and calculate the energy storage-shock sensitive coupling model system, so as to realize the effective simulation of ion redistribution effect of supercapacitor electrolyte under mechanical impact. 2.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,所述储能动力学模型包括双电层效应模型和法拉第反应模型。2. the simulation method of supercapacitor electrolyte ion redistribution effect under mechanical impact according to claim 1, is characterized in that, described energy storage dynamics model comprises electric double layer effect model and Faraday reaction model. 3.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,所述储能动力学模型是通过电极活性材料和电解液材料之间发生的高度可逆的氧化还原反应,使得电荷发生转移产生并存储能量,其中,所述法拉第反应模型储能占主导地位。3. the simulation method of supercapacitor electrolyte ion redistribution effect under mechanical impact according to claim 1, is characterized in that, described energy storage dynamics model is by the highly reversible that takes place between electrode active material and electrolyte material The oxidation-reduction reaction allows charge transfer to generate and store energy, where the Faradaic reaction model dominates energy storage. 4.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,所述储能-冲击敏感耦合模型体系包括双电层效应的电极动力学模型,法拉第反应的电极动力学模型,电解液中的离子扩散模型和力学冲击下的电解液离子重分布模型。4. the simulation method of supercapacitor electrolyte ion redistribution effect under mechanical shock according to claim 1, is characterized in that, described energy storage-shock sensitive coupling model system comprises the electrode dynamics model of electric double layer effect, Faraday The electrode kinetic model of the reaction, the ion diffusion model in the electrolyte and the electrolyte ion redistribution model under mechanical impact. 5.根据权利要求1或4所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,通过所述储能动力学模型对双电层效应模型进行法拉第电流修正,使得所述超级电容器在充电和放电过程中的电势场、离子浓度场在法拉第反应模型和双电层效应模型两种机制共同驱动下的变化趋势得以确定。5. according to the simulation method of supercapacitor electrolyte ion redistribution effect under the described mechanical shock of claim 1 or 4, it is characterized in that, carry out Faraday current correction to electric double layer effect model by described energy storage dynamics model, make The variation trends of the electric potential field and the ion concentration field of the supercapacitor in the process of charging and discharging under the joint drive of the Faraday reaction model and the electric double layer effect model are determined. 6.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,在嵌入力学冲击时,电解液将受惯性力驱动产生定向流动,其流场变化的动力学过程由Navier-Stokes方程进行描述:6. the simulation method of ion redistribution effect of supercapacitor electrolyte under mechanical impact according to claim 1, it is characterized in that, when embedding mechanical impact, electrolyte will be driven by inertial force to produce directional flow, and the change of its flow field The kinetic process is described by the Navier-Stokes equation: 其中,u1,ρ1,μ分别表示电解液流体的流场速度,密度,动力学粘度。Among them, u 1 , ρ 1 , and μ represent the flow field velocity, density, and dynamic viscosity of the electrolyte fluid, respectively. 7.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,通过COMSOL弱形式方程和通用形式偏微分方程的复杂模型构建方法,对锂离子电池模块内嵌模型进行修正,使得所述储能动力学模型符合构建的所述储能-冲击敏感耦合模型体系。7. the simulation method of supercapacitor electrolyte ion redistribution effect under mechanical impact according to claim 1, it is characterized in that, through the complex model building method of COMSOL weak form equation and general form partial differential equation, to lithium ion battery module The embedded model is corrected so that the energy storage dynamic model conforms to the constructed energy storage-shock sensitive coupling model system. 8.根据权利要求7所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,所述通用形式偏微分方程架构如下:8. the simulation method of supercapacitor electrolyte ion redistribution effect under the mechanical shock according to claim 7, is characterized in that, described general form partial differential equation framework is as follows: 其中,u是偏微分方程因变量,t是自变量,其他参数为通用系数。Among them, u is the dependent variable of the partial differential equation, t is the independent variable, and other parameters are general coefficients. 9.根据权利要求1所述的力学冲击下超级电容器电解液离子重分布效应的仿真方法,其特征在于,利用COMSOL软件对所述超级电容器冲击敏感特性进行仿真,所述超级电容器在不同的工作模式下出现两种的冲击敏感现象,其中,第一种情况为当所述超级电容器进行大电流充电时,反应电流得失离子的速率远远大于溶液中离子扩散的速率,第二种情况与所述第一种情况相反,由于所述超级电容器的自放电效应,离子浓度分布已变均衡。9. the simulation method of supercapacitor electrolyte ion redistribution effect under the mechanical impact according to claim 1, it is characterized in that, utilize COMSOL software to simulate the impact sensitive characteristic of described supercapacitor, described supercapacitor is in different work There are two kinds of impact sensitive phenomena in the mode, wherein, the first case is that when the supercapacitor is charged with a large current, the rate of gaining and losing ions by the reaction current is far greater than the rate of ion diffusion in the solution, and the second case is the same as that of the supercapacitor Contrary to the first case described above, the ion concentration distribution has become equalized due to the self-discharge effect of the supercapacitor.
CN201910319320.6A 2019-04-19 2019-04-19 Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect Pending CN110096774A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910319320.6A CN110096774A (en) 2019-04-19 2019-04-19 Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910319320.6A CN110096774A (en) 2019-04-19 2019-04-19 Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect

Publications (1)

Publication Number Publication Date
CN110096774A true CN110096774A (en) 2019-08-06

Family

ID=67445287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910319320.6A Pending CN110096774A (en) 2019-04-19 2019-04-19 Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect

Country Status (1)

Country Link
CN (1) CN110096774A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110264381A1 (en) * 2010-04-22 2011-10-27 Battelle Energy Alliance, Llc Systems, methods and computer readable media to model kinetic performance of rechargeable electrochemical devices
CN105807092A (en) * 2016-03-21 2016-07-27 吉林大学 Method for designing molecule-electron induction type accelerometer elastic film element
CN106940385A (en) * 2016-11-21 2017-07-11 哈尔滨盈江科技有限公司 A kind of temperature compensation for molecular electronic induction type accelerometer

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110264381A1 (en) * 2010-04-22 2011-10-27 Battelle Energy Alliance, Llc Systems, methods and computer readable media to model kinetic performance of rechargeable electrochemical devices
CN105807092A (en) * 2016-03-21 2016-07-27 吉林大学 Method for designing molecule-electron induction type accelerometer elastic film element
CN106940385A (en) * 2016-11-21 2017-07-11 哈尔滨盈江科技有限公司 A kind of temperature compensation for molecular electronic induction type accelerometer

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王晓峰 等: "混合式MEMS超级电容器的动力学建模与仿真方法", 《天津大学学报(自然科学与工程技术版)》 *
郑美娜 等: "超级电容器的热电化学耦合研究", 《电源技术》 *

Similar Documents

Publication Publication Date Title
CN104991980B (en) The electrochemical mechanism modeling method of lithium ion battery
CN105723231B (en) Devices and methods therefor for the state for estimating mixing secondary cell
Gu et al. numerical modeling of coupled electrochemical and transport processes in lead‐acid batteries
CN101629992B (en) Method for estimating residual capacity of iron-lithium phosphate power cell
CN113779794B (en) Lithium-ion battery SOP estimation method and system taking into account microscopic constraints
CN110457742A (en) A lithium-ion battery electrochemical-thermo-mechanical coupling model and modeling method
CN111177924B (en) Establishment method of solid-state battery performance model and prediction method of solid-state battery performance
CN108761341A (en) A kind of lithium ion battery battery chemical modeling parameter acquisition methods
CN111665451B (en) A Lithium-ion battery aging test method under time-varying cycle conditions
US6656234B2 (en) Tuning battery electrode porosity technical field
CN115017781B (en) Lithium ion battery electrochemistry-heat-force-short circuit-thermal runaway coupling model
CN113868934A (en) Parallel lithium ion battery electrochemical parameter identification method
CN110442901A (en) A kind of acquisition methods of lithium ion battery electrochemistry simplified model and its parameter
CN109033619A (en) A kind of transient temperature model modelling approach of 18650 type lithium battery discharge cycles
Fărcaş et al. Modeling and simulation of supercapacitors
CN115189047A (en) Lithium ion battery lithium analysis-free control method and system
CN112924867B (en) Lithium battery capacity attenuation calculation method under multi-field coupling
CN114280480A (en) A method for decomposing the DC internal resistance of lithium-ion batteries based on a numerical model
CN110580389A (en) An optimization method for battery tabs
Liu et al. Design and application: Simplified electrochemical modeling for Lithium-ion batteries in whole lifespan
CN110096774A (en) Mechanics impacts the emulation mode of lower electrolytic solution for super capacitor ion redistribution effect
Pavlov et al. Influence of H2SO4 concentration on lead-acid battery performance: H-type and P-type batteries
CN109164396B (en) Real-time lithium iron phosphate battery life loss assessment method
Ma et al. Multifunctional NMC-Si batteries with self-actuation and self-sensing
CN116663371A (en) Lithium ion battery aging simulation method, system and storage medium

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190806