Formulation, Implementation and Validation of a Two-Fluid model in a Fuel Cell CFD Code

Conference · Mon Dec 01 04:00:00 EST 2008

OSTI ID:1050035

Jain, Kunal ^[1]; Cole, J. Vernon ^[2]; Kumar, Sanjiv ^[3]; Gidwani, Ashok ^[2]; Vaidya, N. ^[4]

---

1. ESI US R&D, San Diego, CA (United States); CFD Research Corporation
2. CFD Research Corporation, Huntsville, AL (United States)
3. Ballard Power Systems, Burnaby, BC (Canada)
4. ESI US R&D, San Diego, CA (United States)

Water management is one of the main challenges in PEM Fuel Cells. While water is essential for membrane electrical conductivity, excess liquid water leads to flooding of catalyst layers. Despite the fact that accurate prediction of two-phase transport is key for optimal water management, understanding of the two-phase transport in fuel cells is relatively poor. Wang et. al. have studied the two-phase transport in the channel and diffusion layer separately using a multiphase mixture model. The model fails to accurately predict saturation values for high humidity inlet streams. Nguyen et. al. developed a two-dimensional, two-phase, isothermal, isobaric, steady state model of the catalyst and gas diffusion layers. The model neglects any liquid in the channel. Djilali et. al. developed a three-dimensional two-phase multicomponent model. The model is an improvement over previous models, but neglects drag between the liquid and the gas phases in the channel. In this work, we present a comprehensive two-fluid model relevant to fuel cells. Models for two-phase transport through Channel, Gas Diffusion Layer (GDL) and Channel-GDL interface, are discussed. In the channel, the gas and liquid pressures are assumed to be same. The surface tension effects in the channel are incorporated using the continuum surface force (CSF) model. The force at the surface is expressed as a volumetric body force and added as a source to the momentum equation. In the GDL, the gas and liquid are assumed to be at different pressures. The difference in the pressures (capillary pressure) is calculated using an empirical correlations. At the Channel-GDL interface, the wall adhesion affects need to be taken into account. SIMPLE-type methods recast the continuity equation into a pressure-correction equation, the solution of which then provides corrections for velocities and pressures. However, in the two-fluid model, the presence of two phasic continuity equations gives more freedom and more complications. A general approach would be to form a mixture continuity equation by linearly combining the phasic continuity equations using appropriate weighting factors. Analogous to mixture equation for pressure correction, a difference equation is used for the volume/phase fraction by taking the difference between the phasic continuity equations. The relative advantages of the above mentioned algorithmic variants for computing pressure correction and volume fractions are discussed and quantitatively assessed. Preliminary model validation is done for each component of the fuel cell. The two-phase transport in the channel is validated using empirical correlations. Transport in the GDL is validated against results obtained from LBM and VOF simulation techniques. The Channel-GDL interface transport will be validated against experiment and empirical correlation of droplet detachment at the interface.

• Research Organization: CFD Research Corporation, Huntsville, AL (United States)
• Sponsoring Organization: USDOE Office of Energy Efficiency and Renewable Energy (EERE)
• Contributing Organization: ESI US R&D, San Diego, CA (United States); Ballard Power Systems, Burnaby, BC (Canada)
• DOE Contract Number: FG36-07GO17010
• OSTI ID: 1050035
• Report Number(s): Doe/GO--17010-4
• Country of Publication: United States
• Language: English

Related Subjects

30 DIRECT ENERGY CONVERSION
fuel cells
gas diffusion layer
model
multicomponent
water management
water transport