CN106159296B - A kind of calculation method of solid oxide fuel cell combination electrode property - Google Patents
A kind of calculation method of solid oxide fuel cell combination electrode property Download PDFInfo
- Publication number
- CN106159296B CN106159296B CN201610270337.3A CN201610270337A CN106159296B CN 106159296 B CN106159296 B CN 106159296B CN 201610270337 A CN201610270337 A CN 201610270337A CN 106159296 B CN106159296 B CN 106159296B
- Authority
- CN
- China
- Prior art keywords
- combination electrode
- particle
- fuel cell
- solid oxide
- oxide fuel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000446 fuel Substances 0.000 title claims abstract description 15
- 238000004364 calculation method Methods 0.000 title claims abstract description 12
- 239000007787 solid Substances 0.000 title claims abstract description 12
- 239000000463 material Substances 0.000 claims abstract description 33
- 150000002500 ions Chemical class 0.000 claims abstract description 12
- 238000009826 distribution Methods 0.000 claims abstract description 9
- 239000008187 granular material Substances 0.000 claims abstract description 7
- 239000011148 porous material Substances 0.000 claims abstract description 3
- 239000000203 mixture Substances 0.000 claims description 44
- 239000002245 particle Substances 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 2
- RVPVRDXYQKGNMQ-UHFFFAOYSA-N lead(2+) Chemical compound [Pb+2] RVPVRDXYQKGNMQ-UHFFFAOYSA-N 0.000 claims description 2
- 230000035699 permeability Effects 0.000 abstract description 3
- 230000000875 corresponding Effects 0.000 abstract description 2
- 229910001233 yttria-stabilized zirconia Inorganic materials 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 4
- 238000010884 ion-beam technique Methods 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000001218 confocal laser scanning microscopy Methods 0.000 description 1
- 238000010192 crystallographic characterization Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006011 modification reaction Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Abstract
The invention discloses a kind of calculation method of solid oxide fuel cell (SOFC) combination electrode property, this method object uses the obtained various random morphologies of SOFC combination electrode of spheric granules random distribution method or electron-microscope scanning method.Regular grid is constructed directly against the combination electrode random structure entirely reconstructed, and assigns corresponding material properties to each grid lattice point by comparison mesh point.Finally by related electrodes properties such as the effective electron conductivity of these lattice point information importing finite element analysis software numerical value calculating combination electrode, effective ion conductivity, permeability, effective three-phase line length, pore radius distributions.To achieve the purpose that being studied different microstructure appearances using the Regular grid with random structure characteristic information acts on combination electrode property.It effectively avoids carrying out subnetting and numerical solution bring technical problem directly against random structure.
Description
Technical field
The present invention relates to the property analysis calculating field of a variety of combination electrodes of fuel cell, specifically a kind of soild oxides
The calculation method of combined electrode of fuel cell property.
Background technique
Solid oxide fuel cell (Solid Oxide Fuel Cell:SOFC) is used as a kind of energy conversion devices, with
Outstanding advantages of it is efficiently, cleaning and fuel forms are various is paid close attention to by people.One typical SOFC is mainly by cathode, electricity
Xie Zhi, anode and bipolar plates composition.With the continuous development of combination electrode function, ingredient and manufacture craft, simple homogeneous hypothesis
Numerical simulation cannot fully meet the needs of SOFC Composites Analysis and structure designs.Because for different
And thus for method for preparing composite electrode and preparation condition, direct result is to result in different electrode microstructure patterns,
Lead to different electrode performances.Therefore, deeply understand that specific random structure pattern is to obtain to the mechanism of action of electrode performance
Optimize the key factor of electrode design.In recent years, focused ion beam flying-spot microscope (Focused Ion-Beam Scanning
Electron Microscopy) and x-ray tomography (X-Ray Tomography) be used to obtain high-resolution
The combination electrode sectional view of rate.By the judgement of different materials phase gray scale in pair cross-section figure, to realize to microscopic electrode structure
Three-dimensionalreconstruction.In this way, we more careful can be accurately obtained the shape characteristic inside combination electrode.But due to high power
The reason of amplification, this method are only used for reconstructing the three-dimensional structure in combination electrode very small part volume.Kee is small within 2010
Group constructs LSM/YSZ combination electrode morphology based on the method for spheric granules random distribution, and attempts using 4 face bodies
Grid directly carries out subnetting to the random structure pattern that YSZ particle is formed.And then in the region charging electric charge conservation equation meter
Calculate the effective ion conductivity of combination electrode.Similar, Mitra group by focused ion beam flying-spot microscope also for being obtained
Combination electrode microstructure solves the effective ion conductivity of combination electrode by way of direct subnetting.This method is people
Good trial has been done in the random structure analysis for further carrying out SOFC electrode.However have two o'clock it is noted that (a) by
In the irregular structure the characteristics of, complicated network results in big calculation amount and numerical instability, therefore this method is only
Can be calculated for one region of very little, and be difficult further to be coupled more equations come it is anti-to the microscopic electrochemical of SOFC
Mechanism is answered to do careful analysis；(b) this method is only used for calculating limited structure and property, is not suitable for by mixing conductance
Material (LSCM of such as anode or the LSCF of cathode) and the pure performance for leading the electrode structure that ionic material (such as YSZ) mixes
Prediction.And the property of the combined electrode structure mixed by three kinds or more materials is predicted.Because while LSCF and YSZ have
There is ionic conductivity, but the two has different intrinsic conductivities.
Summary of the invention
The present invention be exactly in view of a) it is careful consider random morphology for develop the importance of combination electrode with
B) two factors of the defect of existing SOFC combination electrode emulational computation propose a kind of solid oxide fuel cell compound electric
The calculation method of polarity matter.
In order to solve the above technical problems, solution technical solution of the invention is: a kind of solid oxide fuel cell
The calculation method of combination electrode property, includes the following steps:
(1) the porous of LSM/YSZ combination electrode is reconstructed using the method for spheric granules random distribution or electron-microscope scanning figure
Random microstructure, and record the coordinate, radius and attribute of each particle.
(2) according to the length ruler of the porous random microstructure of the LSM/YSZ combination electrode of step (1) reconstruct
It is very little, a regular matrix grid is constructed, and close by compareing each mesh point and reconstructing the mapping of obtained electrode structure
System determines that the material phase attribute of the mesh point is indicated that the x, y, z respectively represent electronic material particle, ion material by x, y, z
Expect particle, gas phase.
(3) by material phase the attribute x, y, z of mesh point in step (2), electronic material particle, ionic material particle, gas phase
Information imports finite element software COMSOL Multiphysics, loaded ion charge conservation equation as input file, and is tying
Potential boundary V1, V2 is arranged in structure both ends, and the electric current I by cross section is obtained by way of numerical solution.
(4) the effective ion conductivity for calculating the material isWherein, A is the combination electrode cross section
Product, L are combination electrode geometric thickness.
The beneficial effects of the invention are that: the method that the present invention uses spheric granules random distribution first reconstructs combination electrode
Random microscopic appearance, and explore the Regular grid characterization combination electrode random structure with different materials phase character information
Pattern, and the method based on Regular grid prediction combination electrode relevant nature, for further research random structure pattern
SOFC combination electrode is laid a good foundation.
Detailed description of the invention
Fig. 1 is a kind of LSM/YSZ combination electrode micro-structure diagram of embodiment reconstruct of the method for the present invention.
Fig. 2 is a kind of LSM/YSZ combination electrode microstructure matrix grid figure of embodiment reconstruct of the method for the present invention.
Fig. 3 is that a kind of LSM/YSZ combination electrode microstructure matrix grid of embodiment reconstruct of the method for the present invention passes through it
The solution schematic diagram of the electric current of cross section.
Specific embodiment
The present invention is done further by taking typical LSM/YSZ binary SOFC combination electrode as an example below with reference to attached drawing 1-3
Explanation, but should not limited the scope of the invention with this.
For convenience of explanation and understand technical solution of the present invention, illustrates the used noun of locality with attached drawing institute below
Subject to the orientation of displaying.
A kind of calculation method of solid oxide fuel cell combination electrode property (abbreviation method) of the present invention, feature exists
In including the following steps:
(1) the porous random microcosmic knot of LSM/YSZ combination electrode is reconstructed using the method for spheric granules random distribution
Structure (referring to Fig. 1), and record the coordinate, radius and attribute of each particle；The attribute of the particle is to lead electric particles or lead
Ion particles.
(2) according to the length ruler of the porous random microstructure of the LSM/YSZ combination electrode of step (1) reconstruct
It is very little, a regular matrix grid (referring to fig. 2) is constructed, and by compareing each mesh point and reconstructing obtained electrode structure
Enantiomorphic relationship, determine that the material phase attribute of the mesh point is indicated by x, y, z, wherein the x, y, z respectively indicate electronics material
Expect particle, ionic material particle, stomata.For example, can table if mesh point coordinate is in the radius for leading ion particles
Be shown as (x, y, z, 0,0).WhereinFor the intrinsic conductivity of YSZ material.In this way, we can be by one
Irregular structure is reduced to the Regular grid input file with different materials attribute information.
(3) by material phase the attribute x, y, z of mesh point in step (2), electronic material particle, ionic material particle, stomata)
Information imports finite element software Comsol Multiphysics, loaded ion charge conservation equation as input file, and is tying
Potential boundary V1, V2 is arranged in structure both ends, and the electric current I by cross section is obtained by way of numerical solution (referring to Fig. 3).
(4) the effective ion conductivity for calculating the material isWherein, A is that the combination electrode geometry is horizontal
Sectional area, L are combination electrode geometric thickness.
May further calculate LSM/YSZ combination electrode some other such as electronic conductivity, porous media permeability,
The nature parameters such as Diffusion of gas stream coefficient, thermal conductivity coefficient and pore-size distribution.
The method of the present invention, first using the random of the method reconstruct multicomponent SOFC combination electrode of spheric granules random distribution
Then morphology.Then Regular grid is constructed directly against the combined electrode structure entirely reconstructed, and passes through comparison grid
Point assigns corresponding material properties to each grid lattice point.Both not using the Regular grid research with random structure characteristic information
Influence with microstructure appearance to combination electrode property.To avoid carrying out subnetting and numerical solution directly against random structure
Bring a series of problems.These lattice point information are finally imported into effective electricity that finite element analysis software numerical value calculates combination electrode
The related electrodes properties such as electron conductivity, effective ion conductivity, permeability, effective three-phase line length, pore radius distribution, thus
Achieve the purpose that carry out exact value solution to random structure using regular grid.The method of the present invention is that further research is random
Then structure and morphology SOFC combination electrode is laid a good foundation.According to the disclosure and teachings of the above specification, of the art
Technical staff can also change and modify the above embodiment.Therefore, the invention is not limited to be disclosed above and retouch
The specific embodiment stated should also be as falling into scope of protection of the claims of the invention to some modifications and changes of the invention
It is interior.In addition, these terms are merely for convenience of description although using some specific terms in this specification, it is not right
The present invention constitutes any restrictions.
The present invention does not address place and is suitable for the prior art.
Claims (4)
1. a kind of calculation method of solid oxide fuel cell combination electrode property, which is characterized in that include the following steps:
(1) 3 dimensions for reconstructing multicomponent SOFC combination electrode using spheric granules random distribution method or electron-microscope scanning figure are porous
Random morphology, and record the coordinate, radius and material properties of each particle；
(2) according to the length size of the porous random microstructure of step (1) reconstruct, a regular matrix net is constructed
Lattice, and by compareing each mesh point and reconstructing the enantiomorphic relationship of obtained electrode structure, determine the material phase of the mesh point
Attribute, the material phase attribute are indicated that x is electronic material particle by x, y, z, and y is ionic material particle, and z is stomata；
(3) finite element software Comsol is imported using the material phase attribute information of mesh point in step (2) as input file
Multiphysics, loaded ion charge conservation equation, and potential boundary V1, V2 are set at structure both ends, pass through numerical solution
Mode obtain the electric current I by cross section；
(4) the effective ion conductivity of calculating multicomponent SOFC combination electrode isWherein, A is the compound electric
Pole geometrical surface, L are combination electrode geometric thickness.
2. a kind of calculation method of solid oxide fuel cell combination electrode property according to claim 1, feature
It is, the material properties of particle are to lead electric particles or lead ion particles in the step (1).
3. a kind of calculation method of solid oxide fuel cell combination electrode property according to claim 1, feature
Be, in the step (1) material properties of particle be mix electric particle i.e. can conductive son also can diversion.
4. a kind of calculation method of solid oxide fuel cell combination electrode property according to claim 1, feature
It is, the material properties of particle are that pore creating material particle can conduction gas in the step (1).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610270337.3A CN106159296B (en) | 2016-04-28 | A kind of calculation method of solid oxide fuel cell combination electrode property |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610270337.3A CN106159296B (en) | 2016-04-28 | A kind of calculation method of solid oxide fuel cell combination electrode property |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106159296A CN106159296A (en) | 2016-11-23 |
| CN106159296B true CN106159296B (en) | 2019-07-16 |
Family
ID=
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013080608A (en) * | 2011-10-04 | 2013-05-02 | Toyota Motor Corp | Electrode conductivity calculation system and electrode conductivity calculation method |
| CN104091960A (en) * | 2014-07-25 | 2014-10-08 | 哈尔滨工业大学 | Method for regulating and controlling microstructure of oxide electrode in working period of fuel cell |
| CN104156539A (en) * | 2014-08-22 | 2014-11-19 | 江苏科技大学 | Method for predicting gas transmission curvature factors of solid oxide fuel cell electrode |
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013080608A (en) * | 2011-10-04 | 2013-05-02 | Toyota Motor Corp | Electrode conductivity calculation system and electrode conductivity calculation method |
| CN104091960A (en) * | 2014-07-25 | 2014-10-08 | 哈尔滨工业大学 | Method for regulating and controlling microstructure of oxide electrode in working period of fuel cell |
| CN104156539A (en) * | 2014-08-22 | 2014-11-19 | 江苏科技大学 | Method for predicting gas transmission curvature factors of solid oxide fuel cell electrode |
Non-Patent Citations (1)
| Title |
|---|
| "Multiscale Model for Solid Oxide Fuel cell with electrode containing mixed conducting material";daifen chen et al.;《AICHE JOURNAL》;20151130;第61卷(第11期);全文 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wilson et al. | Effect of composition of (La0. 8Sr0. 2MnO3–Y2O3-stabilized ZrO2) cathodes: Correlating three-dimensional microstructure and polarization resistance | |
| Iwai et al. | Quantification of SOFC anode microstructure based on dual beam FIB-SEM technique | |
| Cronin et al. | Three-dimensional reconstruction and analysis of an entire solid oxide fuel cell by full-field transmission X-ray microscopy | |
| Lu et al. | Improving intermediate temperature performance of Ni-YSZ cermet anodes for solid oxide fuel cells by liquid infiltration of nickel nanoparticles | |
| Zhu et al. | Enhancement in three-phase boundary of SOFC electrodes by an ion impregnation method: a modeling comparison | |
| Bertei et al. | The fractal nature of the three-phase boundary: a heuristic approach to the degradation of nanostructured solid oxide fuel cell anodes | |
| Wilson et al. | Three-dimensional analysis of solid oxide fuel cell Ni-YSZ anode interconnectivity | |
| Zhang et al. | High-throughput, super-resolution 3D reconstruction of nano-structured solid oxide fuel cell electrodes and quantification of microstructure-property relationships | |
| Zhang et al. | Parametric study for electrode microstructure influence on SOFC performance | |
| Zhao et al. | The effect of catalyst layer design on catalyst utilization in PEMFC studied via stochastic reconstruction method | |
| Zheng et al. | On the tortuosity factor of solid phase in solid oxide fuel cell electrodes | |
| Jankovic et al. | Multiscale imaging and transport modeling for fuel cell electrodes | |
| Zhang et al. | Thermal stress analysis at the interface of cathode and electrolyte in solid oxide fuel cells | |
| Budiman et al. | Evaluation of electrochemical properties of LaNi0. 6Fe0. 4O3− δ-Ce0. 9Gd0. 1O1. 95 composite as air electrode for SOFC | |
| Sar et al. | Three dimensional analysis of Ce0. 9Gd0. 1O1. 95–La0. 6Sr0. 4Co0. 2Fe0. 8O3− δ oxygen electrode for solid oxide cells | |
| Timurkutluk et al. | Synthetical designing of solid oxide fuel cell electrodes: Effect of particle size and volume fraction | |
| Chen et al. | Controllable synthesis of cobalt oxide nanoflakes on three-dimensional porous cobalt networks as high-performance cathode for alkaline hybrid batteries | |
| Riazat et al. | The effect of cathode porosity on solid oxide fuel cell performance | |
| Budáč et al. | Prediction of electrical conductivity of porous composites using a simplified Monte Carlo 3D equivalent electronic circuit network model: LSM-YSZ case study | |
| CN106159296B (en) | A kind of calculation method of solid oxide fuel cell combination electrode property | |
| Bertei et al. | Microstructural modeling and effective properties of infiltrated SOFC electrodes | |
| Setevich et al. | 3D resistor-network modeling of infiltrated SOFC electrodes | |
| Wan et al. | A numerical study of infiltrated solid oxide fuel cell electrode with dual‐phase backbone | |
| Hu et al. | Ultrathin nickel cobalt nitride nanoflowers embedded in Cu3N porous nanorod arrays for ultrahigh capacitance energy storage | |
| Zhang et al. | Tortuosity factor of three-dimensional infiltrate network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201021 Address after: 215500 No.13, Caotang Road, Changshu, Suzhou, Jiangsu Province Patentee after: Changshu intellectual property operation center Co.,Ltd. Address before: 212003 Jiangsu University of Science and Technology, Zhenjiang, Jiangsu province 2 Patentee before: JIANGSU University OF SCIENCE AND TECHNOLOGY |
|
| CP02 | Change in the address of a patent holder |
Address after: 215500 5th floor, building 4, 68 Lianfeng Road, Changfu street, Changshu City, Suzhou City, Jiangsu Province Patentee after: Changshu intellectual property operation center Co.,Ltd. Address before: No.13 caodang Road, Changshu City, Suzhou City, Jiangsu Province Patentee before: Changshu intellectual property operation center Co.,Ltd. |