
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38538652
57846
10.1038/s41598-024-57846-w
Article
Theoretical analysis of the effect of isotropy on the effective diffusion coefficient in the porous and agglomerated phase of the electrodes of a PEMFC
Pacheco C.
Barbosa Romeli romeli.barbosa@cicy.mx
romelix1@gmail.com

Navarro-Montejo A.
Ordoñez L. C.
grid.418270.8 0000 0004 0428 7635 Centro de Investigación Científica de Yucatán, Parque Científico y Tecnológico de Yucatán, Carretera Sierra Papacal–Chuburna Puerto, km 5, Sierra Papacal, 97302 Mérida, Yucatán Mexico
27 3 2024
27 3 2024
2024
14 72996 2 2024
22 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
In polymer membrane fuel cells (PEMFC), the pore microstructure and the effective diffusion coefficient (Deff) of the catalytic layer have a significant impact on the overall performance of the fuel cell. In this work, numerical methods to simulate PEMFC catalytic layers were used to study the effect of isotropy (Ixy) on the Deff. The proposed methodology studies reconstructed systems by Simulated Annealing imaging with different surface fractions of microstructures composed by two diffusive phases: agglomerates and pores. The Deff is determined numerically by the Finite Volume Method solved for Fick's First Law of Diffusion. The results show that the proposed methodology can effectively quantify the effect of isotropy on the Deff for both diffusion phases. Two trends were obtained in the magnitude of the Deff concerning the change in isotropy: (1) an analytical equation is proposed in this article for Deff≥5%D0 and (2) numerical solutions are determined for Deff<5%D0. In our analytical equation are both a lineal and a logarithmic sweep. When the surface fraction is ∅= 50%, the Deff decreases more linearly than ∅=10% at the beginning of the isotropy change, which indicates that small changes in isotropy in the particulate material modify it drastically; under these conditions the diffusion coefficient in the pore is predominant. (3) When the surface fraction is less than 50%, the Deff decreases more exponentially at the beginning and more linearly at the end of the isotropy change, which shows that small isotropy changes in the bar-aligned material drastically alter it. In this trend, diffusion in the agglomerate is less affected by isotropy. The proposed methodology can be used as a design tool to improve the mass transport in porous PEMFC electrodes.

Keywords

Effective transport coefficient
PEMFC
Fick's diffusion
Catalyst layer
Agglomerates
Subject terms

Fuel cells
Fuel cells
Fluidics
Two-dimensional materials
CONAHCYT4047288 04992 Pacheco C. Navarro-Montejo A. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The urgent need to reduce pollution and carbon dioxide in the atmosphere has increased the interest in cleaner and more efficient systems1–3. Proton Exchange Membrane Fuel Cells (PEMFC) are expected to play a critical role in sustainable energy production in the near future as they promise to be energy conversion devices for applications such as automobiles, trains, portable electronic devices, and others due to their high efficiency, fast start-up, and zero emissions when fueled with hydrogen4–7.

However, this technology still faces challenges; one of them is related to mass transport in the agglomerates and porous media, as the performance of a PEMFC is mainly constrained by mass transport limitations. In addition to the effective transport of electrons, ions, and reactant gases are prerequisites for electrochemical reactions to occur8–12. To improve this mechanism, a detailed understanding of the relationship between the properties of the agglomerates and the porous media is required.

To comprehend mass transport in agglomerates and the porous media, one must be aware of three crucial phenomena: molecular diffusion, convective diffusion, and Knudsen diffusion13,14. Molecular diffusion refers to the transport of chemical species under concentration gradients (Fick Diffusion)15. Convection occurs under pressure gradients16. Knudsen diffusion applies to the displacement of flow or gases on small length scales17,18. These three phenomena can be examined at the agglomerates, pore level, or macro scales. However, pore-scale models must be used to account for the structural details of porous media. Furthermore, transport phenomena are linked to the heterogeneous morphology of the components that comprise the catalytic layer (CL), ionomer, and pores19,20.

The catalytic layer is a crucial component of the PEMFC, where the electrochemical reaction rates are controlled. La CL consists of a void region, ionomer, and platinum carbon-supported catalyst (Pt/C)21. The transport phenomena and characteristics of the CL depend on the composite morphological structures, compositions and pore structure (closed pores and dead-end pores)22,23.

Significant progress has been made in recent years in the CL simulation of PEMFC. However, it remains a challenge to evaluate transport in porous media. Nonetheless, various studies have suggested evaluations of pore diffusion behavior. In a recent study, Shin et al.24 proposed an effective utilization factor of the CL related to the tortuosity of the pathways for reagent mass transport within the CL. To simulate the mass transport of reactants within the CLs, this study utilized the Boltzmann method. Numerous studies have reported that Boltzmann method is useful for simulating mass transport phenomena inside heterogeneous porous media25–27. Lange et al.28 showed that the effect of Fick diffusion on CL proved to be significant in numerical estimation. They used stochastic reconstruction to obtain the CL pore distribution and evaluated the effective Diffusion-coefficient for a range of porosity (∅) from 25 to 50%; they found that the oxygen diffusion coefficient is in the range of 0.0025–0.02 cm2 s−1. They concluded that the diffusion transport limitations are underestimated if the Fick diffusion is not considered. The same effect is reported by Yu & Carter29 their calculations of the diffusion coefficient of the whole electrode (gas diffusing layer and CL) showed that ignoring the Fick diffusion underestimates the diffusion limitations. J.O. Ceballos et al.30 present several numerical simulations of the PEM fuel cell and analyze the effect of tortuosity on its overall performance. They present polarization curves evaluated with different effective diffusion models in the anode and cathode. They reported different current density values when employing the various diffusion models. Congfan Zhao et al.31 studied the effect of different catalyst layer designs used in PEMFC using stochastic reconstruction methods. The generation of the catalytic layer consists of 4 phases, firstly, a skeleton of carbon particles, then the generation of Pt particles on the carbon surface, followed by the addition of ionomer. The last phase consists of growing the microstructure with the complete elements. The results presented demonstrate that the fabrication of the MEA using the conventional process with the consideration of carbon selection, porosity, dispersion and ionomer content is of great benefit for the efficient use of Pt in the CL.

On the other hand, several studies have found that the CL can be analyzed as a random heterogeneous material (RHM)32–35; RHMs can be characterized by an Effective Transport Coefficient (ETC), which measures the material’s ability to transport physical properties, and which is affected by the properties of the phases composing the CL layer. Stochastic reconstruction allows defining the microstructure in a mesh of nodes to study the ETC in more detail34,36.

Strategies for the study of RHM have been previously published by our research group. Rodriguez et al.35 analyzed the effect of size reduction on microstructure ETC´s and implemented reconstruction using the SA method. Barbosa et al.36 developed a scaling method to determine the ETC´s in RHM using stochastic reconstruction techniques and analytical techniques. Pacheco et al.37 reported on the effect of topological entropy on the ETC´s in composites aligned and reconstructed by SA. In this research work, the use of stochastic reconstruction using the SA method is employed to analyze the behavior of oxygen diffusion through the ionomer inside the agglomerate and the void space that forms in the CL microstructures that start with an aligned bar architecture and are reconstructed to obtain a dispersed random material. The novelty of our proposal consists of understanding the influence of isotropy on the oxygen diffusion coefficient in both porous and the agglomerate medium to design or analyze catalyst layers of PEMFC. The use of the numerical o theoretical results allows making decisions about the modification of the geometric structure of the materials. Considering in decision making a change of the isotropy or a change of surface fraction. This also helps to improve the accuracy of models and numerical simulations, allowing the optimization of the overall processes. The limitations of simulated annealing found in this study are the boundary conditions, specifically the edge effects since the pixels at the edges result in unwanted effects due to the lack of neighboring pixels in certain directions. This limitation was solved by implementing strategies in the code that allow preserving the continuity of the microstructures without unwanted border effects. The samples studied comprise domains of 350 × 350 pixels where each pixel is equivalent to 0.1 µm, so the samples studies in this work are 35 µm × 35 µm. The evolution of the isotropy and the Deff in the SA agitation steps are systematically evaluated to allow the design of better CLs. The results are obtained using the working group's computational algorithms, which are programmed in C language and developed using the Dev-C +  + compiler (Company Free Software Foundation, Inc. Version 5.11).

Methods

The scope of this work is constrained to isotropic biphasic agglomerates, specifically those consisting of a randomly distributed phase (Pt/C catalyst) within a matrix phase (ionomer). The digital system is spatially defined by a mesh of control volumes (CV). These CV have recognizable phases and diffusivity coefficient values defined at each node. The SA method35,37,38 is used to stir the initial system into a reference system sequentially. The materials are statistically characterized by intrinsic correlation functions in the SA reconstruction methodology. Because of this, it is possible to determine the isotropy (Ixy) at a specific SA stirring interval, defined by the iteration number of the SA method (ISA). The Deff is established by the FVM and is relates to Ixy of the current system. A block diagram of the methodology proposed in this research is presented in Fig. 1.Figure 1 Block diagram of the implemented methodology.

The source system (initial microstructure) is conditioned to random bars ordered from south to north in the sense of the dominant direction of the mass diffusion. The target system (reference or final microstructure) is conditioned to a dispersed system randomly distributed that is a homogeneous, periodic, and isotropic CL. Our initial and reference microstructures consider only one filter: the one-point correlation function (pore fraction). The numerical solution considers the black color as the pore phase and the white color as the agglomerate phase. Agglomerates are made up of Pt/C + Nafion.

It must be pointed out that the unknown dimension of the three-dimensional catalyst layer is considered the same as the aligned dimension in our bidimensional simulated microstructure. In this work the aligned dimension is the vertical south-north direction. On this geometrical point of view, the 3D image is extruded along the unknown direction of our 2D images. On this geometrical assumption, and with the knowledge that the oxygen diffusion coefficient is a directional and intensive property, we can define that our two-dimensional values of the effective diffusion coefficient are the same than the extruded three-dimensional structure. There are three sections in Fig. 1: (1) Top section is for user-defined input that includes the initial and final microstructures of the SA method, CL microstructures are fixed to 10%<∅<50% (with steps of 10%). Surface fractions higher than 50% were not analyzed because experimentally a higher porosity implies a weak microstructure. (2) Left section describes the SA reconstruction strategy for generating the proposed universe of microstructure meshes. (3) Right section shows the numerical outputs that are reported below.

Simulated annealing

The SA methodology is an optimization method used to determine local or global optima in complex problems39,40. In the iteration of the SA method, a pixel of a different phase is exchanged, and the exchange is evaluated and accepted based on a cooling probability based on a fictitious temperature. This approach transforms the original system into a target system. This technique has been extensively employed in reconstructing stochastic heterogeneous materials (SHM)32–34,41–43, demonstrating its versatility and effectiveness. However, it is essential to define some fundamental bases.

The statistical analysis considers the generation of an ensemble (Ω) composed of W realizations ꞷ; in this work, W = 15. Statistically, the realization ꞷ represents the domain V of the studied system (CL) within the measurable space ℜd, partitioned by a j phase. In this context, the region in the space Vjω, the surface fraction ϕj(ω), and the index function Ijω can be defined.1 Ijω=1,ifthepositionx,y∈Vjω.0,otherwise,

The realization ω can be interpreted as a CL of the set Ω of CL´s. All these are microstructures with a finite number of nodes on the x-axis (Nx) and the y-axis (Ny); the specific phase of each node can be identified. This work studies two-phase isometric materials in a square domain of Nx=Ny pixels. In this way, it is possible to represent the agglomerate matrix phase (white pixels) and the pore dispersion phase (black pixels). When examining the white phase, the index function 1 can be rewritten as:2 Ijω=1,ifthepositionpixelx,yiswhiteagglomeratematrix.0,otherwiseitisblackporedispertion.

Once the domain has been defined, applying different correlation functions to obtain statistical information to characterize the material is possible. In this work, the two-point correlation functions (Sj2r) and the linear path correlation function (Pj) were employed in both phases of SA reconstruction.

Considering a material with domain V⊆ℜd, volume V, and composed of j = 2 disjoint random phases. The Sj2r function is the probability that the initial point and the endpoint of a line of length r fall in the same phase j. Considering an isotropic and homogeneous medium. The Sj2r can be defined as a function of distance r28:3 Sj2r=IjxIjx+r,

where 〈〉 refers to the statistical average from evaluating the whole domain ꞷ. IJ is the index function of the computational domain (Eqs. 1 or 2). The function Sj2r provides a measure of how the endpoints of a vector r in the studied phase are correlated. On the other hand, the Pj function is defined by:4 Pj=∫Ijx+arda

Both correlation functions are used for the SA reconstruction; they are considered for both phases (j0 and j1). Equation (5) defines the error between the source system and the current system (ESA):5 ESA=∑rF′r-Fr2

where,6 Fr=14(∑Sj1r+∑Sj2r+∑Pj1r+∑Pj2r)

Fr characterizes the source system, and F′r characterizes the SA current system. The pixel phase shift is accepted if the probability PΔESA is satisfied,7 PΔESA=1,ΔESA≤0,exp-ΔESAT,ΔESA>0

In this work, the SA methodology is applied pixel by pixel, simplifying the theoretical study of the evolution of an initially bar-ordered system as its transitions to a random dispersion system. It should be noted that the user defines both the initial system (ordered bars) and the final system (random dispersion).

Isotropy

The isotropy of a material can be expressed as a function of the orientation of the tensor describing the transport phenomenon44 or as a geometric constraint influencing the covariance matrix in a random field45. The relationship between the sum of covariance models in factor spaces and zonal isotropy and the relationship between geometric isotropy and an isotropic covariance function are relevant aspects in this context46. It is important to note that geometric and zonal isotropy are interrelated.

The concept of geometric isotropy is a valuable tool for describing the behavior of microstructures based on both accurate models, such as images acquired through scanning electron microscopy and synthetic CLs. This approach for describing microstructures is applicable to CL structures of varying dimensions. However, one of the fundamental challenges is the selection of appropriate scaling factors, which is a significant hurdle to overcome in the field of research. In this study, we will use the Sj2r to generate a vertical axis representation Sj,y2r versus the horizontal deviation Sj,x2r,8 Exy=1n∑Sj,yr-Sj,xr2

where, Exy is the mean square error of the correlation functions. Since the user defines the initial and final microstructure of the study, it is feasible to normalize the difference for a generalized comparison. The proposed isotropy (Ixy) is shown in Eq. (9).9 Ixy=1-Exy,actualExy,initial

Exy,initial corresponds to the source system of ordered bars (higher anisotropy) and Exy,actual is the mean of the difference between squares Exy of the actual SA iteration system. The bar microstructure is represented by Ixy=0, as SA shakes the system Ixy tends to 138.

Effective diffusion coefficient

The solution of the FVM is iterative and requires the implementation of the physical model in the nodal point46,47. The molar flux due to diffusion can be calculated from the concentration gradient and a diffusivity coefficient, as are indicated by the first Fick's law described below48:10 JF=-DΔCΔx

where JF is the flux density of the species diffusing through a unit area, D is the diffusion coefficient, ΔC is the differential of the species concentration, and Δx is the diffusion length. Considering the general equation of transport phenomena described in classical fluid mechanics and that the study is limited to a stationary diffusion state, without a convective and transient phenomenon, as well as without generation or consumption of species, the diffusion in a volume can be expressed as in Eq. (11):11 ∫SCn·D∇CdA

Using the Gaussian divergence theorem, we have that the integration over the CV is equal to the integration over the surface bounding that CV. For two-dimension mesh, the diffusion can be expressed as follows:12 De∂xPEAe+Dw∂xWPAw+Ds∂yPSAs+Dn∂xPNAnCP-De∂xPEAeCE-Dw∂xWPAwCW-Ds∂xPSAsCS-Dn∂xPNAnCN=0

where A is the CV cross-sectional area, ∂x is the distance of the volume or control element between boundaries, C is the species concentration, D the phase diffusivity coefficient as a function of the flow direction, the subindices e, w, s, and n are the east, west, south, and north direction, respectively. As well as P is the VC central position. The tri-diagonal matrix algorithm (TDMA) is used to determine the pattern of the complete spatial concentration distribution (Cx,y).

In this study, we define the effective flux (Jeff) as the sum of the local fluxes across horizontal mesh layers. We then utilize the numerical Jeff value and modify the proportionality coefficient in the fundamental Eq. (10), but renaming it to Deff, as it is shown in Eq. (13),13 Deff=JeffNyΔCBoundarycondition

The material properties of the CL mesh define, in a deterministic way, the ETC (Deff). Figure 2 shows FVM boundary conditions for Jeff numerical determination.Figure 2 FVM boundary conditions and graphical index of the local fluxes used for the global JEff.

In this work, we simulate both pore (empty medium) and agglomerate phases (Pt/C + Nafion). The diffusion coefficient of oxygen in an empty medium is DO2,Pore=3.23×10-5m2s-1 and the diffusion coefficient of oxygen in the agglomerate is the same that in the ionomer DO2,Agg=8.45×10-9m2s-149. The exact analytical solution for a parallel diffusion resistance is given by the Eq. (14).14 D0=∅DO2,Pore+1-∅DO2,Agg

where D0 is the exact value of the diffusion coefficient when the both phases are vertically aligned (Ixy=0.0).

Results

The process of reconstruction of the catalytic layers starts with the generation of the synthetic microstructures comprised by the pore phase (black color phase) and the agglomerate (white color phase), the microstructures are bars aligned vertically from north to south and have the following microstructural characteristics: surface fraction (ϕ) of 50% pore—50% agglomerate, 40% pore—60% agglomerate, 30% pore—70% agglomerate, 20% pore—80% agglomerate and 10% pore—90% agglomerate. The domain of all samples is 350 × 350 pixels where each pixel is equivalent to 0.1 µm, so the samples are 35 µm × 35 µm.

The SA method is employed to obtain an ensemble Ω composed of W = 15 random repetitions for each studied CL. The results present the CL isotropy and the Deff determined for both pore and agglomerate phases. Figure 3 shows representative images of the studied CL microstructures. The first column shows porosity, and the C1 column shows the aligned bars when ISA=0 iteration, C2 column is for ISA=5000 iteration, C3 column is for ISA=20000 iteration, and the C4 column shows the dispersed random material (reconstructed image).Figure 3 Representative images of the SA process of the studied CL microstructures.

It is important to highlight that all samples' SA reconstruction temperature was the same ( T=1x10-6) and the reconstructions were stirring until the convergence of the annealing error ESA=1x10-6. The C1 column indicates the bar-ordered CL system, Ixy=0.00±0.00, while the C4 column represents the homogeneously agitated CL system Ixy=1.00±0.00. C2 and C3 columns were selected to analyze their behavior in the results shown in Fig. 6. Each image in Fig. 3 represents a Ω ensemble, but the values are the averages and standard deviation obtained for W = 15.

Figure 4 presents plots based on ISA iterations: (A) for the error variation of the reconstruction process ESA (Eq. 5), (B) for the isotropy variation Ixy (Eq. 9) and (C) for Deff (Eq. 13), all graphs based on the ISA statistical iteration. The values presented in this figure include all the universe of ω realizations of the total samples studied.Figure 4 ESA error response (A), Ixy isotropy response (B), and Deff coefficient (C), the three graphs as a function of ISA iteration.

Figure 4A shows that the ESA value tends to zero, which indicates that the actual microstructure of the CL is converging to the values of the objective microstructure of the CL, as observed in the representative images of Fig. 3. As the ISA iterations increase, the overall error between the reconstructed microstructure and the target microstructure decreases, the reconstruction can be seen in Fig. 3, where in column 1 we have the initial microstructure arranged in aligned bars and in the following columns 2 to 4 the dispersion of the pore phase and the agglomerated phase is observed.

The Ixy values in Fig. 4B tend to unity, indicating that the error between the correlation functions Sj,xr and Sj,yr are decreasing, as indicated in Eq. (9). It is observed that Ixy is more dependent of Φ than ESA, this trend is similar to that reported by38. Figure 4C shows the response of Deff coefficient depending on the reconstruction process iteration ISA. It is observed that the Deff variation is higher for the first ISA iterations.

Figure 5 shows the response of Deff compared to the evolution of isotropy Ixy. It must be pointed out that the numerical solution considers the black color as the pore phase and the white color as the agglomerate phase, both for a synthesized microstructure of a catalytic layer. The values presented in this figure includes all possibilities (points), average values (solid lines) and the curves of Deff after solving the new proposed Eq. (15) (black dotted lines).Figure 5 Deff coefficient as a function of the isotropy evolution Ixy: (A) for all possibilities and analytical solutions magnitudes @ Ixy=0 and (B) for a zoom section and numerical results @ Ixy=1.

Figure 5A shows the evolution of Deff as a function of isotropy Ixy. The results from all examined samples (Ω∈W=15ω) are presented with markers and solid lines denoting their averages. It can be observed that when the Ixy=0 the microstructure corresponds to the pore phase and the agglomerate are vertically aligned. In this condition, it is observed that Deff values are equal to the analytical solution of a parallel diffusion resistances; these values their equations are highlighted in the gray shading of Fig. 5A. The values of D0 obtained by exact analytical solution of Eq. (14) are: D0∅=50=1.6×10-5m2s-1@ (50%DO2,Pore—50%DO2,Agg) and D0∅=10=3.2×10-6m2s-1@ (10%DO2,Pore—90%DO2,Agg). The theoretically based Eq. (15) is proposed for the range Deff≥5%D0. This range is selected by the authors to reduce the variance between the numerical average values and the analytical solution values.15 DeffD0=DeffD0=1-m1Ixy-m2lnIxyIxy0,ifDeff≥5%D0Numericalsolution,ifDeff<5%D0

where D0 is the exact solution of Eq. (14), Ixy0 is the lowest numerical Ixy value different to cero (in this work Ixy0=5.40×10-3) and m1 y m2 are analytical slopes obtained from numerical data, related by the porosity ∅ in Eqs. (16) and (17).16 m1=1.9∅-0.19

17 m2=-0.48∅+0.29

Table 1 shows the solution of Eqs. (14, 16 and 17). These values are applied in Eq. (15) to obtain the black dotted lines in Fig. 5A.Table 1 Numerical values used in equation 15.

Porosity ∅ (%)	D0(m2 s−1)	m1	m2	
10	3.24E−06	0.00	0.24	
20	6.47E−06	0.19	0.19	
30	9.70E−06	0.38	0.15	
40	1.29E−05	0.57	0.10	
50	1.62E−05	0.76	0.05	

Ceballos et al.30 reports effective oxygen diffusion values for a catalytic layer of a PEMFC with 30% pore phase surface fraction; the results reported in their work are obtained by solving different effective diffusion models, they report oxygen Deff in a range of 1.87×10-6m2s-1 for Nam and Kabiany model and 8.91×10-6m2s-1 for Tomadakis and Storirchos model; In this work, for the same pore phase surface fraction condition an analytical result of 9.7×10-6m2s-1 is obtained. The variations in the reported values are attributed to the different microstructures of the evaluated electrodes.

As isotropy increases, Deff values exhibit an exponential decrease, the exponential decrease being particularly pronounced in the sample ϕ=10% (green color) and less pronounced in the sample ϕ=50% (red color). Figure 5B shows a zoom section for details of the range Deff≤1.5×10-7m2s-1 and Ixy≥0.3. In this figure, the diffusion coefficient of oxygen in the agglomerate (data input: DO2,Agg=8.45×10-9m2s-1) is noted by the dashed blue line.

Zhao et al.50 report experimentally obtained magnitudes of the effective diffusivity coefficient for catalytic layers. The reported averages of Deff for catalytic layers range from 4.2±0.9×10-7m2s-1 @ 25 °C to 4.6±0.5×10-7m2s-1 @ 75 °C. Experimental samples exhibited a surface fraction of the pore phase within a range of ϕ = 30% to ϕ = 40%. In this study, the reported values for different surface fractions of the pore phase and agglomerates are in the range of magnitudes of the order of 10-6 when Ixy is equal to zero and 10-8 when Ixy is equal to one. Our numerical model takes into account the diffusive coefficient of oxygen within the agglomerate, attributing that as the pore phase decreases the Deff values approach the reported value of DO2,Agg=8.45×10-9m2s-1. Similarly, Shen et al.51 report experimentally obtained magnitudes of Deff for a catalytic layer with different thicknesses; the average value obtained is 1.47±0.05×10-7m2s-1, a figure that also falls within the range of results obtained in this study. Chen et al.52 presents results of the effective diffusivity coefficient of oxygen for a sample with a surface fraction of 17% for the pore phase using the lattice Boltzmann method, reporting a Deff of 7.71×10-8m2s-1. In our work, the result obtained for a similar surface fraction of the pore phase is 1.4×10-8m2s-1 when it is 20% DO2,Pore- 80% DO2,Agg. The variation in results is attributed to the difference in the isotropy value between synthetic and experimental materials. Microstructural isotropy directly influences the behavior of the effective diffusivity coefficient. The results show that the exponential tendency is preserved at higher isotropy ranges. The box highlighted in the gray shading of Fig. 5B presents Deff@Ixy=1.

As evidenced by the results, the impact of oxygen diffusion within the agglomerate exhibits insignificance concerning the effective diffusion coefficient, Deff, particularly when the surface fraction of the porous phase is predominant. Conversely, in scenarios where the surface fraction of the agglomerate phase predominates, it becomes apparent that the diffusion within the agglomerate significantly influences the ultimate determination of the diffusion coefficient. This observation is elucidated in the zoom section in Fig. 5.

Figure 6 shows Deff as a function of pore surface fraction ϕ, for three selected isotropies referred to in columns of Figs. 3 and 5: the solid black line for Ixy=0.0, the gay line for Ixy=1.0 and dashed black line for Ixy=0.5.Figure 6 Deff coefficient as a function of agglomerate surface fraction ϕ, for different values of isotropy Ixy.

Figure 6 shows the results obtained for Deff as a function of the superficial fraction of the agglomerate phase. The results from all examined samples (Ω∈W=15ω) are presented with markers and solid lines denoting their averages. For Ixy = 0.0 and Ixy = 0.5 the Deff ranges are on the left side of the graph; for Ixy = 1.0 the Deff ranges are on the right side of the graph.

The results are presented within ϕ = 10% to ϕ = 50%, with 10% intervals. The solid black line, labeled as Ixy = 0.0, corresponds to the average of all samples when the microstructure is vertically aligned. Numerical results indicate that the effective diffusion coefficient is directly related to the porosity of the catalyst layer microstructure. Higher porosity has a positive impact on Deff enhancement, this trend agrees with empirical correlation models of the effective diffusion coefficient in porous materials (Bruggeman53, Neale and Nader54, Tomadakis and Sotirchos55, Mezedur et al.56, Zamel et al.57, Das et al.58). The results suggest that the isotropy of the microstructure also have a significant effect on the Deff as shown in Fig. 6. An increasing trend is observed when the porosity of the microstructure is more significant. When the porosity is minor, the diffusion of the agglomerate becomes more relevant in the effective diffusion coefficient, because the higher the surface fraction of the agglomerate, the higher the value of the predominant oxygen diffusion is that of the diffusion in the ionomer, which can be seen in Fig. 650,59 Deff values range from 3.26 × 10–6 m2 s−1@ ϕ = 10% to 1.60 × 10–5 m2 s−1@ ϕ = 50%. This behavior is attributed to the reduced diffusive resistance when the phases are vertically aligned. As isotropy in the microstructure increases, Deff diminishes due to the agitation within the microstructure. When Ixy = 0.5, Deff values are within the range of 3.74 × 10–8 m2 s−1@ ϕ = 10% and 6.53 × 10–6 m2 s−1@ ϕ = 50%. In the case of a randomly monodisperse microstructure (Ixy = 1.0), the value of Deff varies between 1.08 × 10–8 m2 s−1@ ϕ = 10% and 6.22 × 10–8 m2 s−1@ ϕ = 50%. Accordingly, microstructures with lower isotropy (aligned bars) and higher porosity have lower diffusion resistance, which translates as higher Deff.

Conclusions

In this study, a theorical investigation has been carried out to determine the effect of the isotropy on the effective diffusion coefficient of PEMFC catalyst layers that change from aligned to dispersed pattern microstructure. The CL´s studied in this work are composed of two diffusive phases, the porosity, and agglomerated phases. We apply two advanced numerical methodologies: (1) Finite volume method to define a specific oxygen diffusion coefficient into each mesh node: DO2,Pore for pores and DO2,Agg for agglomerates. (2) Simulated Annealing to simulate the evolution of unidirectional aligned bars to random mono-dispersed CL´s.

The analysis confirms that the aligned bar is the better structure for improving the CL mass transport phenomena when the diffusion is made in the parallel direction to the bars because there is less resistance to the transport of the oxygen molecule. Further, this statistical strategy allows us to determine the effect of geometric isotropy (Ixy) on the effective diffusion coefficient (Deff). When the superficial fraction of the pore and the agglomerate are both equal to 50% in a monodisperse random configuration (Ixy = 1.0 + 0.0), Deff exhibits values within a range with a maximum of 6.20 × 10–8 m2 s−1 and a minimum close to the oxygen diffusion coefficient in the ionomer. As the superficial fraction of the pore decreases, the value of Deff decreases significantly. For instance, when the superficial fraction of the pore phase is 10% values close to the diffusion coefficient of oxygen in the ionomer are obtained. As the surface fraction of the agglomerate is predominant, the Deff is mostly influenced by the DO2,Agg, which impacts the overall performance of the fuel cell.

In conclusion, it has been observed that in samples with a surface fraction of the porous phase equal or higher than 50% at isotropy = 0, anisotropic material, the effective diffusion coefficient presents a higher value than when the isotropy is equal to 1. This decrease is attributed to the dispersion of the phases by the SA agitation process, the diffusion coefficient decreases exponentially until reaching a maximum of isotropy. At surface fractions greater than or equal to 50% of the pore phase the diffusion coefficient is influenced by the diffusion of oxygen in the pore, since it is the predominant phase.

On the other hand, as isotropic values approach unity, a change in the behavior of the effective diffusion is observed, with oxygen diffusion in the agglomerate being the predominant factor. This change becomes more evident in samples with porous phase surface fractions below 50%, where the effective diffusion coefficient shows an exponential decrease as the isotropy process progresses.

Ultimately, as the isotropic values approach unity, the effective diffusion coefficient falls into ranges close to the diffusion of oxygen in the ionomer, indicating that this phase becomes predominant in the microstructure. These findings are relevant to understanding and controlling diffusion processes in anisotropic materials, providing valuable information for the design and optimization of devices and systems that rely on the diffusion species.

The results obtained in this study provide information on the behavior of the Diffusion coefficients in different porosity conditions and how this directly affects the effective transport coefficient. In turn, it provides a methodology that is aimed at a better design and manufacture of catalytic layers of PEM fuel cells, leading to a better overall performance of the cell.

Acknowledgements

C. Pacheco acknowledges CONAHCYT for the postdoctoral research fellowship 2022(3), application number 4047288. A. Navarro-Montejo acknowledges CONAHCYT for the doctoral research fellowship, 2022-000018-02NACF-04992.

Author contributions

C.P. conceived and designed the experiments, performed the experiments, data analysis, interpreted, and wrote. R.B. conceptualized the idea of this research project, responsible for research, analysis, interpretation of data, and writing. A.N. and L.C.O. performed the experiments, analysis tools, and writing.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Stern AG A new sustainable hydrogen clean energy paradigm Int. J. Hydrogen Energy 2018 43 4244 4255 10.1016/j.ijhydene.2017.12.180
Stern, A. G. A new sustainable hydrogen clean energy paradigm. Int. J. Hydrogen Energy 43, 4244–4255 (2018).10.1016/j.ijhydene.2017.12.180
2. Nemmour, A., Inayat, A., Janajreh, I. & Ghenai, C. Green hydrogen-based E-fuels (E-methane, E-methanol, E-ammonia) to support clean energy transition: A literature review. Int. J. Hydrogen Energy (2023).
3. Aykut, Y. & Yurtcan, A. B. The role of the EHC system in the transition to a sustainable energy future: A review. Int. J. Hydrogen Energy (2023).
4. Sergi F PEM fuel cells analysis for grid connected applications Int. J. Hydrogen Energy 2011 36 10908 10916 10.1016/j.ijhydene.2011.05.161
Sergi, F. et al. PEM fuel cells analysis for grid connected applications. Int. J. Hydrogen Energy 36, 10908–10916 (2011).10.1016/j.ijhydene.2011.05.161
5. Torreglosa JP Jurado F García P Fernández LM PEM fuel cell modeling using system identification methods for urban transportation applications Int. J. Hydrogen Energy 2011 36 7628 7640 10.1016/j.ijhydene.2011.03.133
Torreglosa, J. P., Jurado, F., García, P. & Fernández, L. M. PEM fuel cell modeling using system identification methods for urban transportation applications. Int. J. Hydrogen Energy 36, 7628–7640 (2011).10.1016/j.ijhydene.2011.03.133
6. Gencoglu MT Ural Z Design of a PEM fuel cell system for residential application Int. J. Hydrogen Energy 2009 34 5242 5248 10.1016/j.ijhydene.2008.09.038
Gencoglu, M. T. & Ural, Z. Design of a PEM fuel cell system for residential application. Int. J. Hydrogen Energy 34, 5242–5248 (2009).10.1016/j.ijhydene.2008.09.038
7. Yuan X A novel feature susceptibility approach for a PEMFC control system based on an improved XGBoost-Boruta algorithm Energy AI 2023 12 100229 10.1016/j.egyai.2023.100229
Yuan, X. et al. A novel feature susceptibility approach for a PEMFC control system based on an improved XGBoost-Boruta algorithm. Energy AI 12, 100229 (2023).10.1016/j.egyai.2023.100229
8. Higareda A Hernández-Arellano DL Ordoñez LC Barbosa R Alonso-Vante N Advanced electrocatalysts for the oxygen evolution reaction: From single-to multielement materials Catalysts 2023 13 1346 10.3390/catal13101346
Higareda, A., Hernández-Arellano, D. L., Ordoñez, L. C., Barbosa, R. & Alonso-Vante, N. Advanced electrocatalysts for the oxygen evolution reaction: From single-to multielement materials. Catalysts 13, 1346 (2023).10.3390/catal13101346
9. Pan W Chen X Wang F Dai G Mass transfer enhancement of PEM fuel cells with optimized flow channel dimensions Int. J. Hydrogen Energy 2021 46 29541 29555 10.1016/j.ijhydene.2020.09.105
Pan, W., Chen, X., Wang, F. & Dai, G. Mass transfer enhancement of PEM fuel cells with optimized flow channel dimensions. Int. J. Hydrogen Energy 46, 29541–29555 (2021).10.1016/j.ijhydene.2020.09.105
10. Carcadea E Effects of geometrical dimensions of flow channels of a large-active-area PEM fuel cell: A CFD study Int. J. Hydrogen Energy 2021 46 13572 13582 10.1016/j.ijhydene.2020.08.150
Carcadea, E. et al. Effects of geometrical dimensions of flow channels of a large-active-area PEM fuel cell: A CFD study. Int. J. Hydrogen Energy 46, 13572–13582 (2021).10.1016/j.ijhydene.2020.08.150
11. Li S Sundén B Effects of gas diffusion layer deformation on the transport phenomena and performance of PEM fuel cells with interdigitated flow fields Int. J. Hydrogen Energy 2018 43 16279 16292 10.1016/j.ijhydene.2018.07.064
Li, S. & Sundén, B. Effects of gas diffusion layer deformation on the transport phenomena and performance of PEM fuel cells with interdigitated flow fields. Int. J. Hydrogen Energy 43, 16279–16292 (2018).10.1016/j.ijhydene.2018.07.064
12. Athanasaki G Jayakumar A Kannan AM Gas diffusion layers for PEM fuel cells: Materials, properties and manufacturing–A review Int. J. Hydrogen Energy 2023 48 2294 2313 10.1016/j.ijhydene.2022.10.058
Athanasaki, G., Jayakumar, A. & Kannan, A. M. Gas diffusion layers for PEM fuel cells: Materials, properties and manufacturing–A review. Int. J. Hydrogen Energy 48, 2294–2313 (2023).10.1016/j.ijhydene.2022.10.058
13. Obayopo SO Bello-Ochende T Meyer JP Modelling and optimization of reactant gas transport in a PEM fuel cell with a transverse pin fin insert in channel flow Int. J. Hydrogen Energy 2012 37 10286 10298 10.1016/j.ijhydene.2012.03.150
Obayopo, S. O., Bello-Ochende, T. & Meyer, J. P. Modelling and optimization of reactant gas transport in a PEM fuel cell with a transverse pin fin insert in channel flow. Int. J. Hydrogen Energy 37, 10286–10298 (2012).10.1016/j.ijhydene.2012.03.150
14. Nguyen VT Pham NH Papavassiliou DV Prediction of the aggregation rate of nanoparticles in porous media in the diffusion-controlled regime Sci. Rep. 2024 14 1916 10.1038/s41598-023-50643-x 38253573
Nguyen, V. T., Pham, N. H. & Papavassiliou, D. V. Prediction of the aggregation rate of nanoparticles in porous media in the diffusion-controlled regime. Sci. Rep. 14, 1916 (2024).38253573 10.1038/s41598-023-50643-x
15. Krishna R Wesselingh JA Taylor R The Maxwell-Stefan approach to mass transfer Chem. Eng. J. 1995 57 861
Krishna, R., Wesselingh, J. A. & Taylor, R. The Maxwell-Stefan approach to mass transfer. Chem. Eng. J. 57, 861 (1995).
16. Chen F Chang M-H Hsieh P-T Two-phase transport in the cathode gas diffusion layer of PEM fuel cell with a gradient in porosity Int. J. Hydrogen Energy 2008 33 2525 2529 10.1016/j.ijhydene.2008.02.077
Chen, F., Chang, M.-H. & Hsieh, P.-T. Two-phase transport in the cathode gas diffusion layer of PEM fuel cell with a gradient in porosity. Int. J. Hydrogen Energy 33, 2525–2529 (2008).10.1016/j.ijhydene.2008.02.077
17. Agarwal RK Yun K-Y Balakrishnan R Erratum:“Beyond Navier–Stokes: Burnett equations for flows in the continuum-transition regime” [Phys. Fluids 13, 3061 (2001)] Phys. Fluids 2002 14 1818 10.1063/1.1466827
Agarwal, R. K., Yun, K.-Y. & Balakrishnan, R. Erratum:“Beyond Navier–Stokes: Burnett equations for flows in the continuum-transition regime” [Phys. Fluids 13, 3061 (2001)]. Phys. Fluids 14, 1818 (2002).10.1063/1.1466827
18. Inoue G Park K So M Kimura N Tsuge Y Microscale simulations of reaction and mass transport in cathode catalyst layer of polymer electrolyte fuel cell Int. J. Hydrogen Energy 2022 47 12665 12683 10.1016/j.ijhydene.2022.02.021
Inoue, G., Park, K., So, M., Kimura, N. & Tsuge, Y. Microscale simulations of reaction and mass transport in cathode catalyst layer of polymer electrolyte fuel cell. Int. J. Hydrogen Energy 47, 12665–12683 (2022).10.1016/j.ijhydene.2022.02.021
19. Carcadea E Influence of catalyst structure on PEM fuel cell performance–A numerical investigation Int. J. Hydrogen Energy 2019 44 12829 12841 10.1016/j.ijhydene.2018.12.155
Carcadea, E. et al. Influence of catalyst structure on PEM fuel cell performance–A numerical investigation. Int. J. Hydrogen Energy 44, 12829–12841 (2019).10.1016/j.ijhydene.2018.12.155
20. Carcadea E PEM fuel cell performance improvement through numerical optimization of the parameters of the porous layers Int. J. Hydrogen Energy 2020 45 7968 7980 10.1016/j.ijhydene.2019.08.219
Carcadea, E. et al. PEM fuel cell performance improvement through numerical optimization of the parameters of the porous layers. Int. J. Hydrogen Energy 45, 7968–7980 (2020).10.1016/j.ijhydene.2019.08.219
21. Molaeimanesh GR Akbari MH Agglomerate modeling of cathode catalyst layer of a PEM fuel cell by the lattice Boltzmann method Int. J. Hydrogen Energy 2015 40 5169 5185 10.1016/j.ijhydene.2015.02.097
Molaeimanesh, G. R. & Akbari, M. H. Agglomerate modeling of cathode catalyst layer of a PEM fuel cell by the lattice Boltzmann method. Int. J. Hydrogen Energy 40, 5169–5185 (2015).10.1016/j.ijhydene.2015.02.097
22. Shojaeefard MH Molaeimanesh GR Nazemian M Moqaddari MR A review on microstructure reconstruction of PEM fuel cells porous electrodes for pore scale simulation Int. J. Hydrogen Energy 2016 41 20276 20293 10.1016/j.ijhydene.2016.08.179
Shojaeefard, M. H., Molaeimanesh, G. R., Nazemian, M. & Moqaddari, M. R. A review on microstructure reconstruction of PEM fuel cells porous electrodes for pore scale simulation. Int. J. Hydrogen Energy 41, 20276–20293 (2016).10.1016/j.ijhydene.2016.08.179
23. Dou S Hao L Liu H Effects of liquid water on the pore structure and transport coefficients in the cathode catalyst layer of PEM fuel cells Int. J. Hydrogen Energy 2022 47 41138 41153 10.1016/j.ijhydene.2022.09.187
Dou, S., Hao, L. & Liu, H. Effects of liquid water on the pore structure and transport coefficients in the cathode catalyst layer of PEM fuel cells. Int. J. Hydrogen Energy 47, 41138–41153 (2022).10.1016/j.ijhydene.2022.09.187
24. Shin S Kim A-R Um S Computational prediction of nanoscale transport characteristics and catalyst utilization in fuel cell catalyst layers by the lattice Boltzmann method Electrochim. Acta 2018 275 87 99 10.1016/j.electacta.2018.04.138
Shin, S., Kim, A.-R. & Um, S. Computational prediction of nanoscale transport characteristics and catalyst utilization in fuel cell catalyst layers by the lattice Boltzmann method. Electrochim. Acta 275, 87–99 (2018).10.1016/j.electacta.2018.04.138
25. Succi S Foti E Higuera F Three-dimensional flows in complex geometries with the lattice Boltzmann method EPL Europhys. Lett. 1989 10 433 10.1209/0295-5075/10/5/008
Succi, S., Foti, E. & Higuera, F. Three-dimensional flows in complex geometries with the lattice Boltzmann method. EPL Europhys. Lett. 10, 433 (1989).10.1209/0295-5075/10/5/008
26. Kim SH Pitsch H Reconstruction and effective transport properties of the catalyst layer in PEM fuel cells J. Electrochem. Soc. 2009 156 B673 B681 10.1149/1.3106136
Kim, S. H. & Pitsch, H. Reconstruction and effective transport properties of the catalyst layer in PEM fuel cells. J. Electrochem. Soc. 156, B673–B681 (2009).10.1149/1.3106136
27. Gao Y Using MRT lattice Boltzmann method to simulate gas flow in simplified catalyst layer for different inlet–outlet pressure ratio Int. J. Heat Mass Transf. 2015 88 122 132 10.1016/j.ijheatmasstransfer.2015.04.031
Gao, Y. Using MRT lattice Boltzmann method to simulate gas flow in simplified catalyst layer for different inlet–outlet pressure ratio. Int. J. Heat Mass Transf. 88, 122–132 (2015).10.1016/j.ijheatmasstransfer.2015.04.031
28. Lange KJ Sui P-C Djilali N Pore scale modeling of a proton exchange membrane fuel cell catalyst layer: Effects of water vapor and temperature J. Power Sources 2011 196 3195 3203 10.1016/j.jpowsour.2010.11.118
Lange, K. J., Sui, P.-C. & Djilali, N. Pore scale modeling of a proton exchange membrane fuel cell catalyst layer: Effects of water vapor and temperature. J. Power Sources 196, 3195–3203 (2011).10.1016/j.jpowsour.2010.11.118
29. Yu Z Carter R Measurements of effective oxygen diffusivity, pore size distribution, and porosity in PEM fuel cell electrodes ECS Trans. 2009 19 1 10.1149/1.3242364
Yu, Z. & Carter, R. Measurements of effective oxygen diffusivity, pore size distribution, and porosity in PEM fuel cell electrodes. ECS Trans. 19, 1 (2009).10.1149/1.3242364
30. Ceballos JO Ordoñez LC Sierra JM Numerical simulation of a PEM fuel cell: Effect of tortuosity parameters on the construction of polarization curves Int. J. Hydrogen Energy 2022 47 30291 30302 10.1016/j.ijhydene.2022.03.112
Ceballos, J. O., Ordoñez, L. C. & Sierra, J. M. Numerical simulation of a PEM fuel cell: Effect of tortuosity parameters on the construction of polarization curves. Int. J. Hydrogen Energy 47, 30291–30302 (2022).10.1016/j.ijhydene.2022.03.112
31. Zhao C The effect of catalyst layer design on catalyst utilization in PEMFC studied via stochastic reconstruction method Energy AI 2023 13 100245 10.1016/j.egyai.2023.100245
Zhao, C. et al. The effect of catalyst layer design on catalyst utilization in PEMFC studied via stochastic reconstruction method. Energy AI 13, 100245 (2023).10.1016/j.egyai.2023.100245
32. Torquato S Random Heterogeneous Materials: Microstructure and Macroscopic Properties 2013 Springer
Torquato, S. Random Heterogeneous Materials: Microstructure and Macroscopic Properties Vol. 16 (Springer, 2013).
33. Torquato, S. Theory of random heterogeneous materials. in Handbook of Materials Modeling 1333–1357 (Springer, 2005).
34. Torquato S Haslach HW Jr Random heterogeneous materials: Microstructure and macroscopic properties Appl. Mech. Rev. 2002 55 B62 B63 10.1115/1.1483342
Torquato, S. & Haslach, H. W. Jr. Random heterogeneous materials: Microstructure and macroscopic properties. Appl. Mech. Rev. 55, B62–B63 (2002).10.1115/1.1483342
35. Rodriguez A Effect of an image resolution change on the effective transport coefficient of heterogeneous materials Materials 2019 12 3757 10.3390/ma12223757 31731587
Rodriguez, A. et al. Effect of an image resolution change on the effective transport coefficient of heterogeneous materials. Materials 12, 3757 (2019).31731587 10.3390/ma12223757
36. Barbosa R Andaverde J Escobar B Cano U Stochastic reconstruction and a scaling method to determine effective transport coefficients of a proton exchange membrane fuel cell catalyst layer J. Power Sources 2011 196 1248 1257 10.1016/j.jpowsour.2010.08.033
Barbosa, R., Andaverde, J., Escobar, B. & Cano, U. Stochastic reconstruction and a scaling method to determine effective transport coefficients of a proton exchange membrane fuel cell catalyst layer. J. Power Sources 196, 1248–1257 (2011).10.1016/j.jpowsour.2010.08.033
37. Pacheco C Numerical simulation to determine the effect of topological entropy on the effective transport coefficient of unidirectional composites Crystals Basel 2020 10 423 10.3390/cryst10060423
Pacheco, C. et al. Numerical simulation to determine the effect of topological entropy on the effective transport coefficient of unidirectional composites. Crystals Basel 10, 423 (2020).10.3390/cryst10060423
38. Escobar B Simulated annealing and finite volume method to study the microstructure isotropy effect on the effective transport coefficient of a 2D unidirectional composite Mater. Today Commun. 2020 24 101343 10.1016/j.mtcomm.2020.101343
Escobar, B. et al. Simulated annealing and finite volume method to study the microstructure isotropy effect on the effective transport coefficient of a 2D unidirectional composite. Mater. Today Commun. 24, 101343 (2020).10.1016/j.mtcomm.2020.101343
39. Kirkpatrick S Gelatt CD Vecchi MP Optimization by simulated annealing Science 1983 1979 220 671 680 10.1126/science.220.4598.671
Kirkpatrick, S., Gelatt, C. D. & Vecchi, M. P. Optimization by simulated annealing. Science 1979(220), 671–680 (1983).10.1126/science.220.4598.671
40. Adachi D Tsujimoto N Akashi R Todo S Tsuneyuki S Search for common minima in joint optimization of multiple cost functions Comput. Phys. Commun. 2019 241 92 97 10.1016/j.cpc.2019.02.004
Adachi, D., Tsujimoto, N., Akashi, R., Todo, S. & Tsuneyuki, S. Search for common minima in joint optimization of multiple cost functions. Comput. Phys. Commun. 241, 92–97 (2019).10.1016/j.cpc.2019.02.004
41. Martins JS Moura CS Vargas RMF Image reconstruction using simulated annealing in electrical impedance tomography: A new approach Inverse Probl. Sci. Eng. 2018 26 834 854 10.1080/17415977.2017.1361420
Martins, J. S., Moura, C. S. & Vargas, R. M. F. Image reconstruction using simulated annealing in electrical impedance tomography: A new approach. Inverse Probl. Sci. Eng. 26, 834–854 (2018).10.1080/17415977.2017.1361420
42. de Castro Martins T Interval simulated annealing applied to electrical impedance tomography image reconstruction with fast objective function evaluation Comput. Math. Appl. 2016 72 1230 1243 10.1016/j.camwa.2016.06.021
de Castro Martins, T. et al. Interval simulated annealing applied to electrical impedance tomography image reconstruction with fast objective function evaluation. Comput. Math. Appl. 72, 1230–1243 (2016).10.1016/j.camwa.2016.06.021
43. Wu W Jiang F Simulated annealing reconstruction and characterization of the three-dimensional microstructure of a LiCoO2 lithium-ion battery cathode Mater. Charact. 2013 80 62 68 10.1016/j.matchar.2013.03.011
Wu, W. & Jiang, F. Simulated annealing reconstruction and characterization of the three-dimensional microstructure of a LiCoO2 lithium-ion battery cathode. Mater. Charact. 80, 62–68 (2013).10.1016/j.matchar.2013.03.011
44. García-Pelayo R Stadler PF Correlation length, isotropy and meta-stable states Physica D 1997 107 240 254 10.1016/S0167-2789(97)00091-2
García-Pelayo, R. & Stadler, P. F. Correlation length, isotropy and meta-stable states. Physica D 107, 240–254 (1997).10.1016/S0167-2789(97)00091-2
45. De Iaco S Posa D Cappello C Maggio S Isotropy, symmetry, separability and strict positive definiteness for covariance functions: A critical review Spat. Stat. 2019 29 89 108 10.1016/j.spasta.2018.09.003
De Iaco, S., Posa, D., Cappello, C. & Maggio, S. Isotropy, symmetry, separability and strict positive definiteness for covariance functions: A critical review. Spat. Stat. 29, 89–108 (2019).10.1016/j.spasta.2018.09.003
46. Patankar S Numerical Heat Transfer and Fluid Flow 2018 Taylor & Francis
Patankar, S. Numerical Heat Transfer and Fluid Flow (Taylor & Francis, 2018).
47. Versteeg HK Malalasekera W An Introduction to Computational Fluid Dynamics: The Finite Volume Method 2007 Pearson education
Versteeg, H. K. & Malalasekera, W. An Introduction to Computational Fluid Dynamics: The Finite Volume Method (Pearson education, 2007).
48. Shin S Liu J Akbar A Um S Nanoscale transport characteristics and catalyst utilization of vertically aligned carbon nanotube catalyst layers for fuel cell applications: Comprehensive stochastic modeling of composite morphological structures J. Catal. 2019 377 465 479 10.1016/j.jcat.2019.07.053
Shin, S., Liu, J., Akbar, A. & Um, S. Nanoscale transport characteristics and catalyst utilization of vertically aligned carbon nanotube catalyst layers for fuel cell applications: Comprehensive stochastic modeling of composite morphological structures. J. Catal. 377, 465–479 (2019).10.1016/j.jcat.2019.07.053
49. Jung C-Y Park C-H Lee Y-M Kim W-J Yi S-C Numerical analysis of catalyst agglomerates and liquid water transport in proton exchange membrane fuel cells Int. J. Hydrogen Energy 2010 35 8433 8445 10.1016/j.ijhydene.2010.05.035
Jung, C.-Y., Park, C.-H., Lee, Y.-M., Kim, W.-J. & Yi, S.-C. Numerical analysis of catalyst agglomerates and liquid water transport in proton exchange membrane fuel cells. Int. J. Hydrogen Energy 35, 8433–8445 (2010).10.1016/j.ijhydene.2010.05.035
50. Zhao J Shahgaldi S Alaefour I Yang S Li X Pore structure and effective diffusion coefficient of catalyzed electrodes in polymer electrolyte membrane fuel cells Int. J. Hydrogen Energy 2018 43 3776 3785 10.1016/j.ijhydene.2018.01.019
Zhao, J., Shahgaldi, S., Alaefour, I., Yang, S. & Li, X. Pore structure and effective diffusion coefficient of catalyzed electrodes in polymer electrolyte membrane fuel cells. Int. J. Hydrogen Energy 43, 3776–3785 (2018).10.1016/j.ijhydene.2018.01.019
51. Shen J Measurement of effective gas diffusion coefficients of catalyst layers of PEM fuel cells with a Loschmidt diffusion cell J. Power Sources 2011 196 674 678 10.1016/j.jpowsour.2010.07.086
Shen, J. et al. Measurement of effective gas diffusion coefficients of catalyst layers of PEM fuel cells with a Loschmidt diffusion cell. J. Power Sources 196, 674–678 (2011).10.1016/j.jpowsour.2010.07.086
52. Chen L Nanoscale simulation of shale transport properties using the lattice Boltzmann method: Permeability and diffusivity Sci. Rep. 2015 5 8089 10.1038/srep08089 25627247
Chen, L. et al. Nanoscale simulation of shale transport properties using the lattice Boltzmann method: Permeability and diffusivity. Sci. Rep. 5, 8089 (2015).25627247 10.1038/srep08089
53. von Bruggeman DAG Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen: III: Die elastischen Konstanten der quasiisotropen Mischkörper aus isotropen Substanzen Ann. Phys. 1937 421 160 178 10.1002/andp.19374210205
von Bruggeman, D. A. G. Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen: III: Die elastischen Konstanten der quasiisotropen Mischkörper aus isotropen Substanzen. Ann. Phys. 421, 160–178 (1937).10.1002/andp.19374210205
54. Neale GH Nader WK Prediction of transport processes within porous media: Diffusive flow processes within an homogeneous swarm of spherical particles AIChE J. 1973 19 112 119 10.1002/aic.690190116
Neale, G. H. & Nader, W. K. Prediction of transport processes within porous media: Diffusive flow processes within an homogeneous swarm of spherical particles. AIChE J. 19, 112–119 (1973).10.1002/aic.690190116
55. Tomadakis MM Sotirchos SV Ordinary and transition regime diffusion in random fiber structures AIChE J. 1993 39 397 412 10.1002/aic.690390304
Tomadakis, M. M. & Sotirchos, S. V. Ordinary and transition regime diffusion in random fiber structures. AIChE J. 39, 397–412 (1993).10.1002/aic.690390304
56. Mezedur MM Kaviany M Moore W Effect of pore structure, randomness and size on effective mass diffusivity AIChE J. 2002 48 15 24 10.1002/aic.690480104
Mezedur, M. M., Kaviany, M. & Moore, W. Effect of pore structure, randomness and size on effective mass diffusivity. AIChE J. 48, 15–24 (2002).10.1002/aic.690480104
57. Zamel N Li X Shen J Correlation for the effective gas diffusion coefficient in carbon paper diffusion media Energy Fuels 2009 23 6070 6078 10.1021/ef900653x
Zamel, N., Li, X. & Shen, J. Correlation for the effective gas diffusion coefficient in carbon paper diffusion media. Energy Fuels 23, 6070–6078 (2009).10.1021/ef900653x
58. Das PK Li X Liu Z-S Effective transport coefficients in PEM fuel cell catalyst and gas diffusion layers: Beyond Bruggeman approximation Appl. Energy 2010 87 2785 2796 10.1016/j.apenergy.2009.05.006
Das, P. K., Li, X. & Liu, Z.-S. Effective transport coefficients in PEM fuel cell catalyst and gas diffusion layers: Beyond Bruggeman approximation. Appl. Energy 87, 2785–2796 (2010).10.1016/j.apenergy.2009.05.006
59. Cunningham RE Williams RJJ Diffusion in Gases and Porous Media 1980 Springer
Cunningham, R. E. & Williams, R. J. J. Diffusion in Gases and Porous Media Vol. 1 (Springer, 1980).
