Polymer electrolyte fuel cell (PEFC) is promising as a next generation power source for vehicles due to the high efficiency. To make the PEFC commercially competitive, the reduction of the cost is necessary. For the reduction of the cost of PEFC, it is required to increase the maximum current density. Under the high current density operation, the water produced in a catalyst layer (CL) accumulates in a gas diffusion layer (GDL) and the supply of oxygen is suppressed. Therefore, the water transport in the GDL is important to increase the maximum current density of PEFC. To improve the water transport in the GDL, the wettability distribution is important factor. Utaka el al. showed that better cell performance is obtained by forming the hydrophilic region in the hydrophobic GDL because the water transport in the GDL is controlled by the hydrophilic region (1). Similarly, Forner-Cuenca et al. showed that the cell performance is improved by the GDL with the hydrophilic region (2). There are some reports which showed the water transport is improved by the control of the wettability of the GDL. However, these reports are focused on the planer distribution in the GDL and the effect of the wettability distribution in the thickness direction is not clear. For the optimization of GDL structure, it is necessary to clarify the effect of the wettability distribution in the thickness direction. Therefore, in this study, the effect of the wettability distribution in the thickness direction was investigated. To evaluate the water transport in the GDL with the wettability distribution in the thickness direction, a lattice Boltzmann method (LBM) was used (3). Figure 1 shows the simulation domain and boundaries. The size of GDL is 80μm × 80μm × 80μm and the porosity of GDL is 70%. The bottom of the simulation domain is assumed the micro-porous layer (MPL) and it is treated as a solid wall. MPL crack as inlet boundary is set in the MPL. The top of the simulation domain is outlet boundary. The four sides of the simulation domain are solid walls. The contact angle of GDL, MPL and side walls is set 130º. The density of the two-phase in the simulation is set identical to water, ρ = 978kg/m3. The interfacial tension between the two-phase is σ = 6.40×10-2N/m, and liquid and gas viscosities are μL = 4.04×10-4Pa∙s and μG = 2.30×10-5Pa∙s respectively. Figure 2 shows simulation results of water transport in GDL. Figure 2 (a) shows the result for the GDL without the wettability distribution. The contact angle of the GDL is uniform (130 º). Figure 2 (b) shows the result for the GDL with the wettability distribution in the thickness direction. The contact angle from the GDL surface to the depth of 16 μm is 50º and the other regions of the GDL is 130º. In these figures, the water is shown in red. With the water transport in the GDL without the wettability distribution (Fig. 2 (a)), the water was inflowed to the large pore in the center of the GDL (0.32ms). The water was spread in the large pore (0.96ms). After the much water was accumulated in the GDL, the water was reached at the GDL surface (1.9ms). With the water transport in the GDL with the wettability distribution in the thickness direction (Fig. 2 (b)), the water was inflowed to the large pore in the center of the GDL and the water accumulated in the GDL as well as the GDL without the wettability distribution (0.32ms, 0.96ms). When the water was reached at the hydrophilic region in the top of GDL, the water in the GDL was pulled up to the GDL surface and the amount of the water in the center of the GDL was reduced (1.2ms). Compared with the GDL without wettability distribution, the time when the water reached the GDL surface was earlier in the GDL with wettability distribution. Therefore, the water transport in the thickness direction is improved by the wettability distribution in the thickness direction. In the final full paper, the large-scale simulation including the rib and channel will be carried out and the effect of the wettability in the thickness direction will be discussed. Reference (1) Y. Utaka, R. Koresawa, J. Power Sources, 363, 227-233 (2017). (2) A. Forner-Cuenca, J. Biesdorf, V. Manzi-Orezzoli, L. Gubler, T. J. Schmidt, P. Boillat, J. Electrochem. Soc., 163 (13), F1389 (2016). (3) S. Sakaida, Y. Tabe, and T. Chikahisa, J. Power Sources, 361, 133-143 (2017). Figure 1
Read full abstract