Numerical simulation of two-phase flow in a multi-gas channel of a proton exchange membrane fuel cell

Understanding the two-phase distribution characteristics within the multi-gas channel of a fuel cell is important for improving fuel cell performance. In the paper, the volume of fluid model is used to predict the dynamic behaviour of water in the multi-gas channel, analyze the pressure drop, velocity distribution, and flow resistance coefficient between different channels, and investigate the influence of operating conditions, surface wettability and channel structure on the two-phase distribution characteristics in the channel. The results show that water undergoes the processes of growth, separation, single droplet transport, wall impact, droplet collision, liquid film formation, and liquid film transport in the multi-gas channel. Inlet velocity and surface wettability significantly affect the pressure drop, water saturation, and surface water coverage. As the inlet velocity and gas diffusion layer surface wettability increase, the flow resistance coefficient and unevenness of the distribution decrease, indicating that the in-channel flow distribution homogeneity is enhanced. The rectangular channel has better water removal and flow distribution uniformity than the tapered channel, and the unevenness of distribution decreases significantly with decreasing rectangular width, from 0.15715 to 0.00315. The research work is a guide to understanding water transport in multi-gas channels, accelerating water removal, and improving inter-channel flow distribution uniformity.     

An investigation of two-phase flow pressure drop in helical rectangular channel

Numerical simulations have been carried out to evaluate the two-phase frictional pressure drop for air-water two-phase flow in horizontal helical rectangular channels by varying configurations, inlet velocity and inlet sectional liquid holdup. The investigations performed using eight coils, five different inlet velocity and four different inlet sectional liquid holdups. The effects of curvature, torsion, fluid velocity and inlet sectional liquid holdup on two-phase frictional pressure drop have been illustrated. It is found that the two-phase frictional pressure drop relates strongly to the superficial velocities of air or water, and that the curvature and torsion have some effect on the pressure drop for higher Reynolds number flows in large-scale helical rectangular channel; the inlet sectional liquid holdup only increases the magnitude of pressure drop in helical channel and has no effect on the development of pressure drop. The correlation developed predicts the two-phase frictional pressure drop in helical rectangular channel with acceptable statistical accuracy.     

Water management studies in PEM fuel cells,Part II: Ex situ investigation of flow maldistribution,pressure drop and two-phase flow pattern in gas channels

Two-phase flow of water and reactant gases in the gas distribution channels of proton exchange membrane fuel cells (PEMFCs) plays a critical role in proper water management. In this work, the two-phase flow in PEMFC cathode parallel channels is studied over a wide range of superficial air velocity (air stoichiometry) and superficial water velocity in a specially designed ex situ experimental setup, which enables the measurement of instantaneous flow rates in individual gas channels and simultaneous visualization of the water flow structure. It is found that the two-phase flow at low superficial air velocities (air stoichiometry below 5) is dominated by slugs or semi-slugs, leading to severe flow maldistribution and large fluctuations in the pressure drop. Slug residence time, measured from the video observation and the instantaneous flow rate data, is found to be a new parameter to describe the slug flow. At higher air velocities, a water film is formed on the channel walls if they are hydrophilic. The pressure drop for the film flow is characterized by smaller but frequent fluctuations, which are found to result from the water buildup at the channel-exit manifold interface. As the superficial air velocity increases further, mist flow is obtained where little water buildup is observed. The water buildup in the gas channels at the two-phase flow is well described by the two-phase friction multiplier, defined as the ratio of the two-phase pressure drop to the single gas phase pressure drop. It is found that the two-phase friction multiplier increases with increasing water flow rate. A flow pattern map is developed using superficial water and air velocities with clearly defined transition regions.     

Numerical investigation of liquid water distribution in the cathode side of proton exchange membrane fuel cell and its effects on cell performance

A three-dimensional unsteady two-phase model for the cathode side of proton exchange membrane fuel cell (PEMFC) consisting of gas diffusion layer (GDL) with hybrid structural model is developed to investigate liquid water behaviors under different operating and geometrical conditions and to quantitatively evaluate effects of liquid water distribution on reactant transport and current density distribution. Simulation results reveal that liquid water transport processes and distributions are significantly affected by inlet air velocity, wall wettability and water inlet position, which in turn play a prominent role on local reactant transport and cause considerable disturbances of the current density. Liquid water film spreading on the gas channel (GC) top wall is identified as the most desirable flow pattern in the GC based on overall evaluations of current density magnitude, uniformity of current density distribution and pressure drop in the GC. Modification to GDL structure is proposed to promote the formation of the desirable flow pattern.     

Characteristics of Two-Phase Flows in a Rectangular Microchannel with a T-Junction Type Gas-Liquid Mixer

In this study, gas–liquid two-phase flows in a horizontal rectangular microchannel have been investigated. The rectangular microchannel has a hydraulic diameter of 0.235 mm, and a width and depth of 0.24 mm and 0.23 mm, respectively. A T-junction-type gas–liquid mixer was used to introduce gas and liquid in the channel. In order to know the effects of liquid properties, distilled water, ethanol, and HFE7200 were used as the test liquids, with nitrogen gas was used as the test gas. The flow pattern, the bubble length, the liquid slug length, and the bubble velocity in two-phase flow were measured with a high-speed video camera, and the void fraction was determined from the bubble velocity data and the superficial gas velocity data. In addition, the pressure drop was also measured with a calibrated differential pressure transducer. The bubble length data were compared with the calculation by the scaling law proposed by Garstecki et al. [7]. The bubble velocity data and/or the void fraction data were well correlated with the well-known drift flux model [12] with a new distribution parameter correlation developed in this study. The frictional pressure drop data were also well correlated with the Lockhart-Martinelli method with a correlation of the two-phase friction multiplier.     

Three-dimension simulation of two-phase flows in a thin gas flow channel of PEM fuel cell using a volume of fluid method

This study investigates the two-phase flow in a thin gas flow channel of PEM fuel cells and wall contact angle's impact using the volume of fluid (VOF) method with tracked two-phase interface. The VOF results are compared with experimental data, theoretical solution and analytical data in terms of flow pattern, pressure drop and water fraction. Stable film flow is predicted, as observed experimentally, for the contact angle ranging from 5° to 40° including varying contact angles at different walls of a channel. The contact angle is found to have small impact on the gas pressure drop for the stratified flow regime, but it determines the meniscus of the two-phase interface, which affects the optical detection of the liquid thickness in experiment. The work is important to study of two-phase flow dynamics, multichannel design, experimental design and control of two-phase flows in thin gas flow channels for PEM fuel cells.     

Fundamental data on the gas–liquid two-phase flow in minichannels

We report on the results of investigations into the characteristics of an air–water isothermal two-phase flow in minichannels, that is, in capillary tubes with inner diameters of 1 mm, 2.4 mm, and 4.9 mm, also in capillary rectangular channels with an aspect ratio of 1 to 9. The directions of flow were vertical upward, horizontal and vertical downward. Based on the authors 15 years of fundamental research into the gas–liquid two-phase flows in circular tubes and rectangular channels, we summarized the characteristics of the flow phenomena in a minichannel with special attention on the flow patterns, the time varying holdup and the pressure loss. The effects of the tube diameters and aspect ratios of the channels on these flow parameters and the flow patterns were investigated. Also the correlations of the holdup and the frictional pressure drop were proposed.     

Modeling gas flow in PEMFC channels: Part I – Flow pattern transitions and pressure drop in a simulated ex situ channel with uniform water injection through the GDL

The two-phase flow in the gas channels of a proton exchange membrane fuel cell (PEMFC) is studied with an ex situ setup using a gas diffusion layer (GDL) as the sidewall of the channels. Air is supplied at the channel inlet manifold and water is supplied continuously and uniformly through the GDL along the length of the channel. This is different from the simultaneous air and water introduction at the inlet of the channel as studied by previous two-phase flow researchers. The GDL is compressed between the gas channels and the water chambers to simulate PEMFC conditions. The superficial velocity for air and water ranged from 0.25 to 34.5 m/s and 1.54 × 10⁻⁵ to 1.54 × 10⁻⁴ m/s, respectively. The ex situ setup was run in both vertical and horizontal orientations with two GDLs, – Baseline (Mitsubishi Rayon Co. MRC 105 with 5 wt.% PTFE and coated with an in-house MPL by General Motors) and SGL 25 BC – and three channel treatments – hydrophobic, hydrophilic, and untreated Lexan, with contact angles of 116°, 11° and 86°, respectively. No appreciable effect was noted because of the orientation, GDL type or channel coatings. The flow regime is observed at different locations along the channel and is expressed as a function of the superficial air and water velocities. Flow regime criteria are developed and validated against the range of ex situ data observations. A new variable water flow rate pressure drop model is developed in order to account for the variation of water entering the channel at multiple locations along the flow length. Pressure drop models are developed for specific flow regimes and validated against experimental data. The models are able to predict the experimental pressure drop data with a mean error of less than 14%.     

Flow Boiling Heat Transfer Characteristics in Minitubes With and Without Hydrophobicity Coating

Abstract

Wettability plays an important role during flow boiling inside micro and mini channels. The present work focuses on the flow boiling heat transfer characteristics inside copper minitube (inner diameter of 3?mm) coated internally to render the inside surface nearly hydrophobic. Electroless Galvanic Deposition technique is employed for hydrophobic coating inside the copper tube. Both single phase heat transfer and two-phase flow boiling heat transfer and pressure drop characteristics were investigated in regular and internally coated hydrophobic copper minitubes. The experiments were performed with deionized water as a working fluid and the mass flux was varied from 100 to 650?kg/m²s. The two-phase heat transfer characteristics was observed to be both functions of mass flux as well as heat flux. The two phase heat transfer has been observed to be augmented due to the wettability within the tubes. The two-phase pressure drop has also been observed to increase when compared to the regular, uncoated tube; however, the proportional increment is lower than the augmentation achieved in two-phase heat transfer. The enhanced heat transfer effects observed have been explained on the basis of wetting physics.     

Three dimensional numerical simulation of gas-liquid two-phase flow patterns in a polymer-electrolyte membrane fuel cells gas flow channel

Water management in polymer-electrolyte membrane fuel cells (PEMFCs) has a major impact on fuel cell performance and durability. To investigate the two-phase flow patterns in PEMFC gas flow channels, the volume of fluid (VOF) method was employed to simulate the air-water flow in a 3D cuboid channel with a 1.0 mm × 1.0 mm square cross section and a 100 mm in length. The microstructure of gas diffusion layers (GDLs) was simplified by a number of representative opening pores on the 2D GDL surface. Water was injected from those pores to simulate water generation by the electrochemical reaction at the cathode side. Operating conditions and material properties were selected according to realistic fuel cell operating conditions. The water injection rate was also amplified 10 times, 100 times and 1000 times to study the flow pattern formation and transition in the channel. Simulation results show that, as the flow develops, the flow pattern evolves from corner droplet flow to top wall film flow, then annular flow, and finally slug flow. The total pressure drop increases exponentially with the increase in water volume fraction, which suggests that water accumulation should be avoided to reduce parasitic energy loss. The effect of material wettability was also studied by changing the contact angle of the GDL surface and channel walls, separately. It is shown that using a more hydrophobic GDL surface is helpful to expel water from the GDL surface, but increases the pressure drop. Using a more hydrophilic channel wall reduces the pressure drop, but increases the water residence time and water coverage of the GDL surface.     

Forchheimer's inertial effect on liquid water removal in proton exchange membrane fuel cells with baffled flow channels

In this study, a two-dimensional, two-phase, non-isothermal and steady-state modified model of proton exchange membrane fuel cells is developed. The Forchheimer's effect (Non-Darcy effect) is coupled in the model, and its impact on liquid water removing process in flow channels with baffles having different shapes is discussed. Simulation results show that the liquid water is able to be removed more at the regions around baffles. At the same time, the baffle shapes reform the liquid water distribution. When using the baffles having larger dimensions (e.g. using rectangular baffles or trapezoidal baffles), the flow spaces around baffles decrease more and the liquid water is removed more because of the increase in local flow velocity. As a result, the concentration polarization is weakened and cell performance is improved more. Moreover, a streamline baffled flow channel that is designed for the purpose of both increasing the cell performance and avoiding excessive increase in pressure drops is discussed. Simulation results show that this flow channel design can both avoid too much increase in pressure drop and facilitate the liquid water removing out from the fuel cell.     

Gas–liquid two-phase flow patterns in parallel channels for fuel cells

Two-phase flow in horizontal parallel channels has been experimentally investigated under fuel cell related operating conditions. Pronounced hysteresis is observed in the pressure drop versus flow characteristic curve when starting from either flooded or dry conditions. When gas is introduced into channels initially filled with water (flooded initial condition), both gas and liquid tend to flow predominantly in one channel at low gas or liquid flow velocities. As the gas flow velocity increases, even distribution of gas and liquid flow in both channels is observed, accompanied with a sudden decrease in the pressure drop. On the other hand, even gas and liquid flow distribution between both channels is found at comparatively lower gas flow velocities when starting with dry-gas flow conditions with liquid introduced into channels filled with gas (stratified flow regime). The flow regimes of this system are visualized in plots of the pressure drop against gas and liquid flow velocities. However, this phenomenon tends to vanish at high gas and liquid flow velocities, suggesting that high gas and liquid flow velocities are required to ensure even flow distribution in parallel channels. The hysteresis points appear at the same level of the pressure drop, reflecting intrinsic characteristics of the parallel channels used in this study. These results have important implications for PEM fuel cell operational strategies. In order to avoid reactant mal-distribution in parallel flow channels in the flow field in the two-phase flow regime, fuel cells should be operated at sufficiently high gas flow velocities.     

Experimental Investigation and RELAP5 Modeling of Two-Phase Flow in Horizontal Rectangular Channel

The investigation of steam, water, and air flow characteristics in horizontal channel is a part of major investigations program at the Lithuanian Energy Institute. The objective of this program is to identify condensation effects on two-phase flow stability and to predict conditions when rapid condensation could be induced in two-phase condensable flow. This article presents investigation of steam–water and air multiphase flow in nearly horizontal rectangular channel. The experimental data for pressure drop and interfacial and wall shear stresses in the channel with uniform distribution of void fraction are presented in this paper. Overall channel dimensions are length = 1.2 m, width = 0.02 m, height = 0.1 m; however, the test section was about 0.84 m in length. Three different flow types were analyzed at atmospheric pressure: (1) single-phase air flow (height of the channel was reduced to 0.075 m); (2) non-condensable air–water two-phase flow at void fraction of 0.75; (3) two-phase steam–water flow at almost saturation conditions, and void fraction of 0.75. RELAP5 Mod3.3 code was selected to model test cases. Modeling results and experimental measurements show good agreement with each other. The developed model will be used for calculating different cases of the process.     

Gas flow rate distributions in parallel minichannels for polymer electrolyte membrane fuel cells: Experiments and theoretical analysis

In the present work, instantaneous gas flow rates in each of two parallel channels of gas-liquid two-phase flow systems were investigated through measurements of the pressure drop across the entrance region. Liquid flow rates in two branches were pre-determined through liquid injection independently into each channel. Experiments were conducted in two different manners, i.e., the gas flow rate was varied in both ascending and descending paths. Flow hysteresis was observed in both gas flow rate distributions and the overall pressure drop of two-phase flow systems. Effects of liquid flow rates on gas flow distributions were examined experimentally. The presence of flow hysteresis was found to be associated with different flow patterns at different combinations of gas and liquid flow rates and flow instability conditions. A new and simple method was developed to predict gas flow distributions based on flow regime-specific pressure drop models for different experimental approaches and flow patterns. In particular, two different two-phase pressure drop models were used for slug flow and annular flow, separately. Good agreement was achieved between theoretical predictions and our experimental data. The developed new method can be potentially applied to predict gas flow distributions in parallel channels for fuel cells.     

Two-phase flow and maldistribution in gas channels of a polymer electrolyte fuel cell

Liquid water transport in a polymer electrolyte fuel cell (PEFC) is a major issue for automotive applications. Mist flow with tiny droplets suspended in gas has been commonly assumed for channel flow while two-phase flow has been modeled in other cell components. However, experimental studies have found that two-phase flow in the channels has a profound effect on PEFC performance, stability and durability. Therefore, a complete two-phase flow model is developed in this work for PEFC including two-phase flow in both anode and cathode channels. The model is validated against experimental data of the wetted area ratio and pressure drop in the cathode side. Due to the intrusion of soft gas diffusion layer (GDL) material in the channels, flow resistance is higher in some channels than in others. The resulting flow maldistribution among PEFC channels is of great concern because non-uniform distributions of fuel and oxidizer result in non-uniform reaction rates and thus adversely affect PEFC performance and durability. The two-phase flow maldistribution among the parallel channels in an operating PEFC is explored in detail.     

Numerical study of water management in the air flow channel of a PEM fuel cell cathode

The water management in the air flow channel of a proton exchange membrane (PEM) fuel cell cathode is numerically investigated using the FLUENT software package. By enabling the volume of fraction (VOF) model, the air–water two-phase flow can be simulated under different operating conditions. The effects of channel surface hydrophilicity, channel geometry, and air inlet velocity on water behavior, water content inside the channel, and two-phase pressure drop are discussed in detail. The results of the quasi-steady-state simulations show that: (1) the hydrophilicity of reactant flow channel surface is critical for water management in order to facilitate water transport along channel surfaces or edges; (2) hydrophilic surfaces also increase pressure drop due to liquid water spreading; (3) a sharp corner channel design could benefit water management because it facilitates water accumulation and provides paths for water transport along channel surface opposite to gas diffusion layer; (4) the two-phase pressure drop inside the air flow channel increases almost linearly with increasing air inlet velocity.     

General characteristics of two-phase flow distribution in a compact heat exchanger

An experimental loop representing a compact plate heat exchanger was built up to study the two-phase distribution in the different header channels. The test section consists of a cylindrical horizontal header and eight rectangular channels in which the liquid and vapour flow rates are evaluated and the flow inside the header can be visualized. Several geometrical and functional parameters to study the two-phase distribution were tested using “HFE 7100” at a temperature close to 57 °C and a pressure close to 100 kPa. A flow pattern map in the header was built up using the different entry parameters on which a quantitative understanding of the two-phase distribution could be deduced.     

A critical review of two-phase flow in gas flow channels of proton exchange membrane fuel cells

Water management in PEM fuel cells has received extensive attention due to its key role in fuel cell performance. The unavoidable water, from humidified gas streams and electrochemical reaction, leads to gas-liquid two-phase flow in the flow channels of fuel cells. The presence of two-phase flow increases the complexity in water management in PEM fuel cells, which remains a challenging hurdle in the commercialization of this technology. Unique water emergence from the gas diffusion layer, which is different from conventional gas-liquid two-phase flow where water is introduced from the inlet together with the gas, leads to different gas-liquid flow behaviors, including pressure drop, flow pattern, and liquid holdup along flow field channels. These parameters are critical in flow field design and fuel cell operation and therefore two-phase flow has received increasing attention in recent years. This review emphasizes gas-liquid two-phase flow in minichannels or microchannels related to PEM fuel cell applications. In situ and ex situ experimental setups have been utilized to visualize and quantify two-phase flow phenomena in terms of flow regime maps, flow maldistribution, and pressure drop measurements. Work should continue to make the results more relevant for operating PEM fuel cells. Numerical simulations have progressed greatly, but conditions relevant to the length scales and time scales experienced by an operating fuel cell have not been realized. Several mitigation strategies exist to deal with two-phase flow, but often at the expense of overall cell performance due to parasitic power losses. Thus, experimentation and simulation must continue to progress in order to develop a full understanding of two-phase flow phenomena so that meaningful mitigation strategies can be implemented.     

Modeling two-phase flow in PEM fuel cell channels

This paper is concerned with the simultaneous flow of liquid water and gaseous reactants in mini-channels of a proton exchange membrane (PEM) fuel cell. Envisaging the mini-channels as structured and ordered porous media, we develop a continuum model of two-phase channel flow based on two-phase Darcy's law and the M² formalism, which allow estimate of the parameters key to fuel cell operation such as overall pressure drop and liquid saturation profiles along the axial flow direction. Analytical solutions of liquid water saturation and species concentrations along the channel are derived to explore the dependences of these physical variables vital to cell performance on operating parameters such as flow stoichiometric ratio and relative humility. The two-phase channel model is further implemented for three-dimensional numerical simulations of two-phase, multi-component transport in a single fuel-cell channel. Three issues critical to optimizing channel design and mitigating channel flooding in PEM fuel cells are fully discussed: liquid water buildup towards the fuel cell outlet, saturation spike in the vicinity of flow cross-sectional heterogeneity, and two-phase pressure drop. Both the two-phase model and analytical solutions presented in this paper may be applicable to more general two-phase flow phenomena through mini- and micro-channels.     

A novel solution for pressure drop in singly connected microchannels of arbitrary cross-section

This paper outlines a novel approximate solution for determining the pressure drop of fully developed, laminar, single-phase flow in singly connected microchannels of arbitrary cross-section. Using a “bottom-up” approach, it is shown that for constant fluid properties and flow rate in fixed cross-section channels, the Poiseuille number is only a function of geometrical parameters of the cross-section, i.e., perimeter, area, and polar moment of inertia. The proposed model is validated with experimental data for rectangular, trapezoidal, and triangular microchannels. The model is also compared against numerical results for a wide variety of channel cross-sections including: hyperellipse, trapezoid, sine, square duct with two adjacent round corners, rhombic, circular sector, circular segment, annular sector, rectangular with semi-circular ends, and moon-shaped channels. The model predicts the pressure drop for the cross-sections listed within 8% of the values published.