Flow Field Effect on the Performance of Direct Formic Acid Membraneless Fuel Cells: A Numerical Study

: The performance of both air-breathing and air-feeding direct formic acid membraneless fuel cells (DFAMFCs) possessing different ﬂow ﬁelds were numerically investigated in this study at given concentration and ﬂow rate for both fuel and electrolyte. Single serpentine, stepwise broadening serpentine, multi-serpentine and parallel channel were tested as liquid ﬂow ﬁeld, while single serpentine, stepwise broadening serpentine, multi-serpentine and pin channel were tested as air ﬂow ﬁeld. The channel width was either 0.8 mm or 1.3 mm. The simulation results showed that the air-breathing DFAMFC having identical ﬂow ﬁeld for both fuel and electrolyte yielded highest cell output. The air-breathing DFAMFC having SBS liquid ﬂow ﬁeld yielded a maximum power density of 10.5 mW/cm 2 , while the air-breathing DFAMFC having S(1.3) liquid ﬂow ﬁeld produced an open circuit voltage of 1.0 V owing to few formic acid penetration into the cathode. Concerning the air-feeding DFAMFCs, the DFAMFC having SBS liquid ﬂow ﬁeld and MS(0.8) air ﬂow ﬁeld yielded highest peak power density, 12 mW/cm 2 , at an airﬂow rate of 500 sccm. Considering the power generated by the DFAMFCs together with the power consumed by the air pump, DFAMFC having SBS liquid ﬂow ﬁeld and Pin(0.8) air ﬂow ﬁeld could be the preferred design. Author Contributions: Software, S.-H.H.; validation, S.-H.H.; formal analysis, J.-C.S. and S.-H.H.; investigation, J.-C.S.; writing—original draft preparation, J.-C.S. and S.-H.H.; writing—review and editing, J.-C.S.; supervision, J.-C.S.; project administration, J.-C.S.


Introduction
Direct liquid fuel cells with polymer electrolyte membrane are becoming highly regarded worldwide because they are favorable to be power sources for portable electronics at room temperature. However, owing to the significant fuel crossover of the polymer electrolyte membrane while feeding liquid fuel [1,2], one of the particular types of fuel cells that generates electricity at room temperature by feeding liquid electrolyte between both electrodes instead of installing a polymer electrolyte membrane have been developed for about a decade. Numerous studies on exploring the characteristics of membraneless fuel cells have been conducted based on both experimental test and numerical simulation [3][4][5][6]. Those membraneless fuel cell usually featured a single, short and straight microchannel, and is also termed as microfluidic fuel cells.
Unlike the traditional DMFCs that possess a solid polymer membrane between two electrodes, a flowing electrolyte-direct methanol fuel cell (FE-DMFC) having a single serpentine, parallel serpentine, triple serpentine and a grid type cathode flow field were simulated by Ouellette et al. [7] to examine the cell performance. The three-dimensional fuel cell model featured an electrolyte stream flowing between two polymer electrolyte membranes, which individually assembled an electrode on one side, for eliminating the fuel crossover. Results showed that the single serpentine flow field provided the best performance due to the better reactant distribution over the cathode catalyst layer. Similar fuel cell configuration using a circulating liquid electrolyte between two MEAs of the PEMFCs to tackle both water and thermal problems was also proposed and tested with dry reactant gases [8].
Besides the fuel crossover issue of the direct liquid fuel cells, flow field design for both electrodes of direct liquid fuel cell has also been considered as another crucial factor that affects the fuel cell performance. An optimal flow field design for fuel cells causes a uniform fuel and oxidant distribution on the electrode surface and ensures the products, H 2 O and CO 2 for example, to be easily removed from the active surface area.
Numerous studies have been conducted to investigate the effect of flow field designs on the DMFCs performance. Among those studies, single serpentine, parallel, and multiple serpentine flow fields are commonly employed as the anodic flow field [9][10][11]. In addition, a few particular anodic flow fields were also proposed. Wang et al. [12] numerically investigated the anodic flow velocity and temperature distributions of the DMFC with four different anodic flow fields, including double serpentine, parallel, helix and single serpentine, using 3-D simulation model. However, the helix anodic flow field yielded a non-uniform inner temperature distribution. Yuan et al. [13] tested an active liquid-feed DMFC with three different types of anodic flow fields, including traditional right-angle serpentine flow field, rounded-corner serpentine flow field, and stepwise broadening serpentine flow field, and investigated the effects of those flow fields on the cell performance, gas bubble behavior and pressure drop characteristics of the DMFC. The results show that the broadening serpentine flow field contributes to the CO 2 emission and uniform distribution of reactants, respectively, thus yielding the highest power output. El-Zoheiry et al. [14] proposed new multi-path spiral flow field designs to improve the under-rib convection mass transport, and consequently the performance of direct methanol fuel cells. Their simulation results indicate a significant increase in fuel cell performance with the enhancement of convection mass transport.
Besides the anodic flow field design, the serpentine, multiple serpentine, and parallel, as well as grid cathode flow fields were also studied for DMFCs in order to examine the cell performance over the past few years [15][16][17][18][19]. Jung et al. [15] found that the gas is well distributed within serpentine flow field because of the favorable mass transport and pressure drop while extremely non-uniform oxygen concentration distribution is observed within parallel flow field. The results obtained by Oliveira et al. [16] show that for the lower value of fuel cell temperature and lower value of methanol concentration, the use of multiple serpentine flow field for cathode leads to a better performance, because the pressure-driven mass flow in the channels ensures a good ability for water removal. Jung et al. [17] found that serpentine flow fields exhibited the highest cell output for low stoichiometric factor operation, while serpentine flow field and grid flow field revealed comparable cell performance for medium stoichiometric factor operation.
Moreover, the performance of DMFCs having bio-inspired flow fields [20,21] was also numerically investigated, in order to achieve both a uniform reactant distribution on the electrode surface and a favorable pressure drop between the inlet and outlet of the reactant.
The flow field effect on the direct liquid fuel cells performance that explored in the aforementioned literature are the fuel cells possessing a polymer electrolyte membrane or multiple polymer electrolyte membranes [8]. In order to eliminate the fuel crossover of direct liquid fuel cells as much as possible, designing direct liquid fuel cells with an electrolyte stream flowing between two gas diffusion electrodes is considered as an effective way. However, such design causes the fuel cells to possess an additional flow field for liquid electrolyte, resulting in more complicated fuel cell design to analyze. According to the literature survey in the previous paragraph, it has been confirmed that the flow fields for both anode and cathode of the direct liquid fuel cells play important roles in the fuel cell performance. Besides, an additional flow field for liquid electrolyte imposes an additional factor on the performance of the direct liquid membraneless fuel cells that has not yet been investigated using three-dimensional numerical simulation, this study aims to numerically investigate the performance of both air-breathing and air-feeding direct formic acid membraneless fuel cells (DFAMFCs) with different flow fields.

Computational Domain
The numerical model of the two different types of DFAMFCs in Figure 1 accounts for the coupled transport of fluid, species and charge in conjunction with the electrochemical reactions. Similar model was also employed in the studies [22,23] using commercially available CFD software. The detailed dimensions of the fuel cell geometry are summarized in Table 1. The following assumptions were made while performing the simulation.
(1) The three-dimensional system is at steady state and isothermal.
(2) Laminar, incompressible fluid flow and body force is negligible.
(3) The physical properties of the electrodes are isotropic and homogeneous. (4) The solutions are dilute and uniformly mixed. (5) Proton transport from anode to cathode is by electromigration only. (6) Electromigration of formate ions is negligible due to the high concentration of supporting electrolyte. (7) The product of carbon dioxide is fully dissolved in the solution. (8) Oxygen transport in the porous gas diffusion electrode is by diffusion only. numerically investigate the performance of both air-breathing and air-feeding direct formic acid membraneless fuel cells (DFAMFCs) with different flow fields.

Computational Domain
The numerical model of the two different types of DFAMFCs in Figure 1 accounts for the coupled transport of fluid, species and charge in conjunction with the electrochemical reactions. Similar model was also employed in the studies [22,23] using commercially available CFD software. The detailed dimensions of the fuel cell geometry are summarized in Table 1. The following assumptions were made while performing the simulation.
(1) The three-dimensional system is at steady state and isothermal.
(2) Laminar, incompressible fluid flow and body force is negligible.
(3) The physical properties of the electrodes are isotropic and homogeneous.
(4) The solutions are dilute and uniformly mixed.
(5) Proton transport from anode to cathode is by electromigration only.

Fluid Flow
The laminar, incompressible fluid flow in the fuel, liquid electrolyte and oxidant flow channels is governed by the continuity and Navier-Stokes equation as follows, where ρ, → u , µ, and p are the fluid density, the velocity vector, the dynamic viscosity of fluid, and the static pressure, respectively. Within the porous media including the catalyst layer and gas diffusion layer, namely carbon paper, the fluid flow is governed by the continuity equation as Equation (1) and the Brinkman equation as follows, where ε p and κ denote the porosity and the permeability of the porous medium, respectively. An identical flow rate of 2.0 mL/min is prescribed at both the fuel and electrolyte inlet boundaries, while the outlet of both the fuel and electrolyte is assumed to be at atmospheric pressure. No-slip and impermeable boundary conditions are applied to all channel walls and continuity is maintained at the interfaces between the flow channel and the porous media.

Species Transport
The species transport in both the flow channel and porous gas diffusion layer by both diffusion and advection is governed as follows, However, since both the formic acid and oxygen are consumed within the catalyst layer, the species transport in the porous catalyst layer should be: where D i,eff and c i denote the effective diffusion coefficient of each species within the porous medium and the local species concentration, and source term S i represents the consumption of reactant within the catalyst layer. The effective diffusion coefficient of the species within the gas diffusion layer and the catalyst layer is, respectively expressed as where τ is the tortuosity of the carbon paper and D i is the species diffusion coefficient in a flow channel. A constant value of the fuel concentration is set at the inlet, and a convective boundary condition is imposed at the outlet. The continuity is maintained at the interfaces between the flow channels and the porous media. All other walls are set to be zero-flux boundaries. In addition, the oxygen concentration at the cathode/air interface is prescribed as a constant value and all other walls are set as zero-flux boundary condition. Table 2 shows the velocity and concentration boundary conditions at all inlets of both types of fuel cells in the study. Besides, the physical and transport properties used in the simulation are listed in Table 3.

Electrochemical Reaction Kinetics
The following two electrochemical reactions indicating the formic acid oxidation at the anode and the oxygen reduction at the cathode are considered here.
In brief, the liquid fuel, HCOOH, is oxidized to produce electrons and protons on the anode following Equation (5a). The former travel along the external circuit of the system through the load and generate the electricity, while the latter move to the cathode through the liquid electrolyte, H 2 SO 4 in the study. Subsequently, the protons, electrons and the oxygen molecules in the cathode react to achieve the reduction reaction indicate as Equation (5b).
The source term in Equation (3b) for anode, S f , and cathode, S o , can be expressed in terms of both the active specific area of the catalyst layer on both electrodes (a v,a and a v,c ), and the local current density of the catalyst layer in both electrodes (i loc,a and i loc,c ), respectively as follows, where i loc is estimated with both the limiting current density, i lim , depending on the maximum reactant transport rate, and the current density generated on the electrode surface, i exp r , as follows, where d is the diffusion distance, n is the number of electrons transported on the electrode, and D i,eff is the effective diffusion coefficient of the reactant within the porous medium. The i expr in Equation (7) is determined by the concentration-dependent Butler-Volmer eqauation for both anode and cathode, respectivley, as follows, where i o,a and i o,c are the exchange current density on the anode and cathode, respectively, c f and c o are the local fuel and oxygen concentrations, respectively, c f ,re f and c o,re f are the reference fuel and oxygen concentrations, respectively, β is the reaction order (β = 1), α a and α c are the charge transfer coefficients, and η is the activation overpotential, which can be estimated by where E eq , φ l and φ s are the equilibrium potential, the local electrolyte and electric potential, respectively. Note that M cross in Equation (9b) denotes the effective fuel crossover flux at the electrolyte/cathode catalyst layer interface. The reaction kinetics used in the abovementioned equations are listed in Table 4.

Electric Current Flow
Both the local electrolyte potential and local electric potential in Equation (10) can be determined by: where σ l and σ s are the conductivity of the liquid electrolyte and effective conductivity of the electrode that is determined by the following equations where ε denotes the porosity of either the catalyst layer or the gas diffusion layer. Once the local electric potential is determined, the electric current flows through the electrolyte, → i l , and the electrode, → i s , is expressed as follows, The anode side wall is set to zero potential and the cathode end face is set to the cell voltage. All other boundaries are zero-flux boundaries.

Solution Procedure
All unknowns including → u , p, c f , c o , φ l and φ s are solved by means of coupling the continuity equation, momentum equation, species transport equation and the electrochemical reactions, as well as charge transfer mentioned above using the commercial software, COMSOL Multiphysics 5.2. Structured grids were established to discretize the flow field, while unstructured grids were established to discretize the rest of the computational domain. Some of the structured grids were converted in order to join the structured grids and unstructured grids at the interface. Grid independence was ensured by testing numerous sets of meshes for the computational domain. Once the differences in both the current density and the average formic acid concentration within the catalyst layer at 0.4 V between two successive sets of meshes was less than 1%, the grid independence was considered to be achieved. According to the test result revealed in Figure 2, the final solutions were obtained with mesh elements of about 100,000 and 240,000 for air-breathing and air-feeding DFAMFCs in this study. Parametric sweep and a direct solution procedure with PARDISO and MUMPS solvers are employed to calculate the current density at different voltages with the convergence criteria for the residual of 10 −4 . Finally, the current density, j, of the fuel cells is calculated using the following equation, where A R is the entire electrode surface, 20 mm × 21.8 mm, and i s,n is the component of the electric current normal to the surface of the catalyst layer. In order to validate the simulation results, the polarization curves of the air-breathing DFAMFCs having S(0.8) as both fuel and electrolyte flow fields obtained by simulation were compared with the measured values in the experiment [26]. Note that the configuration and dimensions of the air-breathing DFAMFCs in that experiment [26] were identical to the model depicted in Figure 1a and the Nafion content in the anode catalyst layer of the DFAMFC was 3.73 mg/cm 2 with catalyst loading of 2 mg-Pd/cm 2 and 2 mg-Pt/cm 2 on the anode and cathode, respectively. Figure 2b shows that both data are in fairly good agreement especially at both high fluid flow rate and concentration, for example, the results presented by purple dashed curve and purple hollow circles in Figure 2b with fluid flow rate of 2.0 mL/min and reactant concentration of 3.0 M, because of the minor effect of gas bubble produced on the catalyst layer on the cell performance that was not considered in the simulation model.

Results and Discussion
Various configurations of the liquid flow fields for transporting the fuel and electrolyte and the gas flow fields for air distribution shown in Table 5 were numerically investigated here. Table 6 shows the channel width, channel depth and rib width of each flow field design. Note that the rib width between two adjacent channels, and the depth of both channel and rib are 0.8 mm and 1.0 mm, respectively, for each flow field. Table 5. The flow fields of the DFAMFCs simulated in this study with the inlet indicated with a pink circle.

Liquid Flow Field Gas Flow Field
Single Serpentine (S) Single Serpentine (S) Figure 2. Results of (a) grid independence test, and (b) the comparison of V-I curves between measured data and the numerical Scheme 0. As both fuel and electrolyte flow fields.
In order to validate the simulation results, the polarization curves of the air-breathing DFAMFCs having S(0.8) as both fuel and electrolyte flow fields obtained by simulation were compared with the measured values in the experiment [26]. Note that the configuration and dimensions of the air-breathing DFAMFCs in that experiment [26] were identical to the model depicted in Figure 1a and the Nafion content in the anode catalyst layer of the DFAMFC was 3.73 mg/cm 2 with catalyst loading of 2 mg-Pd/cm 2 and 2 mg-Pt/cm 2 on the anode and cathode, respectively. Figure 2b shows that both data are in fairly good agreement especially at both high fluid flow rate and concentration, for example, the results presented by purple dashed curve and purple hollow circles in Figure 2b with fluid flow rate of 2.0 mL/min and reactant concentration of 3.0 M, because of the minor effect of gas bubble produced on the catalyst layer on the cell performance that was not considered in the simulation model.

Results and Discussion
Various configurations of the liquid flow fields for transporting the fuel and electrolyte and the gas flow fields for air distribution shown in Table 5 were numerically investigated here. Table 6 shows the channel width, channel depth and rib width of each flow field design. Note that the rib width between two adjacent channels, and the depth of both channel and rib are 0.8 mm and 1.0 mm, respectively, for each flow field. Table 5. The flow fields of the DFAMFCs simulated in this study with the inlet indicated with a pink circle.

Liquid Flow Field Gas Flow Field
Single Serpentine (S)

Results and Discussion
Various configurations of the liquid flow fields for transporting the fuel and electrolyte and the gas flow fields for air distribution shown in Table 5 were numerically investigated here. Table 6 shows the channel width, channel depth and rib width of each flow field design. Note that the rib width between two adjacent channels, and the depth of both channel and rib are 0.8 mm and 1.0 mm, respectively, for each flow field. Table 5. The flow fields of the DFAMFCs simulated in this study with the inlet indicated with a pink circle.

Liquid Flow Field Gas Flow Field
Single Serpentine (S) Single Serpentine (S) Single Serpentine (S)

Results and Discussion
Various configurations of the liquid flow fields for transporting the fuel and electrolyte and the gas flow fields for air distribution shown in Table 5 were numerically investigated here. Table 6 shows the channel width, channel depth and rib width of each flow field design. Note that the rib width between two adjacent channels, and the depth of both channel and rib are 0.8 mm and 1.0 mm, respectively, for each flow field. Table 5. The flow fields of the DFAMFCs simulated in this study with the inlet indicated with a pink circle.

Liquid Flow Field Gas Flow Field
Single Serpentine (S) Single Serpentine (S) Stepwise Broadening Serpentine (SBS)         Figure 3 shows the V-I curves (solid ones) and P-I curves (dashed ones) of the airbreathing membraneless fuel cells having three different flow field configurations, including the 0.8-mm-wide serpentine flow field for both fuel and electrolyte, denoted as S(0.8) × S(0.8), the 0.8-mm-wide multiple serpentine fuel flow field and 0.8-mm-wide serpentine   Figure 3 shows the V-I curves (solid ones) and P-I curves (dashed ones) of the airbreathing membraneless fuel cells having three different flow field configurations, including the 0.8-mm-wide serpentine flow field for both fuel and electrolyte, denoted as S(0.8) × S(0.8), the 0.8-mm-wide multiple serpentine fuel flow field and 0.8-mm-wide serpentine   Figure 3 shows the V-I curves (solid ones) and P-I curves (dashed ones) of the airbreathing membraneless fuel cells having three different flow field configurations, including the 0.8-mm-wide serpentine flow field for both fuel and electrolyte, denoted as S(0.8) × S(0.8), the 0.8-mm-wide multiple serpentine fuel flow field and 0.8-mm-wide serpentine electrolyte flow field, denoted as MS(0.8) × S(0.8), and the 0.8-mm-wide parallel fuel flow  Figure 3 shows the V-I curves (solid ones) and P-I curves (dashed ones) of the airbreathing membraneless fuel cells having three different flow field configurations, including the 0.8-mm-wide serpentine flow field for both fuel and electrolyte, denoted as S(0.8) × S(0.8), the 0.8-mm-wide multiple serpentine fuel flow field and 0.8-mm-wide serpentine electrolyte flow field, denoted as MS(0.8) × S(0.8), and the 0.8-mm-wide parallel fuel flow field and 0.8-mm-wide serpentine electrolyte flow field, denoted as P(0.8) × S(0.8). It can be found in Figure 3 Figure 3. In addition, it is observed in Figure 3 that the DFAMFC having identical flow field for both fuel and electrolyte has highest maximum power density.
Processes 2021, 9, x FOR PEER REVIEW 11 of all because of the significant ohmic loss as shown in Figure 3. In addition, it is observe in Figure 3 that the DFAMFC having identical flow field for both fuel and electrolyte ha highest maximum power density.   Figure 4a is not only significant higher than that in both MS(0.8) and P(0.8) flow fields, but also more uniform along th entire channel, resulting in a highest pressure drop between the inlet and the outlet of fu as shown in Figure 4d. Besides, the fuel stream of extremely low velocity is observed the parallel channel of P(0.8) flow field in Figure 4c.     Figure 5 is the top surface of the anode catalyst layer which contacts the fuel flow field, while the backmost section is the bottom surface of the anode catalyst layer which attaches to the gas diffusion layer. It can be seen in Figure 5a that the uniformity of the formic acid concentration distribution on the top surface of the anode catalyst layer is the best in S(0.8) fuel flow field in Figure 5a, followed by MS(0.8) fuel flow field in Figure  5b, and is the worst in P(0.8) fuel flow field in Figure 5c.   Figure 5 is the top surface of the anode catalyst layer which contacts the fuel flow field, while the backmost section is the bottom surface of the anode catalyst layer which attaches to the gas diffusion layer. It can be seen in Figure 5a that the uniformity of the formic acid concentration distribution on the top surface of the anode catalyst layer is the best in S(0.8) fuel flow field in Figure 5a Figure 7a between the cathode catalyst layer and the electrolyte flow field is plotted in Figure 7b. Due to the lack of the polymer membrane between two GDEs, the formic acid indeed invades the cathode catalyst layer, and the formic acid concentration becomes increased from the first channel to the lower reach for each fluid flow field combination as revealed in Figure 7b. Each peak of those curves in Figure 7b corresponds to the boundary of the fluid channel. Note that the illustrations above Figure  7b denote the position of the fluid channel depicted as blue region and the rib marked with slashes for the three fuel and electrolyte flow field combinations along the centerline in Figure 7a. As observed in Figure 5, as the channel width of the flow field is increased,   Figure 7a between the cathode catalyst layer and the electrolyte flow field is plotted in Figure 7b. Due to the lack of the polymer membrane between two GDEs, the formic acid indeed invades the cathode catalyst layer, and the formic acid concentration becomes increased from the first channel to the lower reach for each fluid flow field combination as revealed in Figure 7b. Each peak of those curves in Figure 7b corresponds to the boundary of the fluid channel. Note that the illustrations above Figure 7b denote the position of the fluid channel depicted as blue region and the rib marked with slashes for the three fuel and electrolyte flow field combinations along the centerline in Figure 7a. As observed in Figure 5, as the channel width of the flow field is increased, the pressure established in the flow field reduces at a given flow rate, resulting in less driving force of the fuel crossover to the cathode catalyst layer and higher open circuit voltage. Therefore, as SBS flow field was employed, whose channel width is stepwise broadened from 0.8 mm to 1.3 mm, the formic acid concentration in Figure 7b is second highest, corresponding to second highest open circuit voltage in Figure 6.
Processes 2021, 9, x FOR PEER REVIEW 14 of 2 the pressure established in the flow field reduces at a given flow rate, resulting in les driving force of the fuel crossover to the cathode catalyst layer and higher open circu voltage. Therefore, as SBS flow field was employed, whose channel width is stepwis broadened from 0.8 mm to 1.3 mm, the formic acid concentration in Figure 7b is secon highest, corresponding to second highest open circuit voltage in Figure 6. As the upstream channel width of SBS fluid flow field is narrower than that of S(1.3 flow field, the pressure of the fuel flow in SBS is higher than that in S(1.3), causing a fa vorable condition for formic acid to penetrate into the anode catalyst layer. In addition t the slight fuel crossover of the DFAMFC possessing SBS fluid flow field, the air-breathin the pressure established in the flow field reduces at a given flow rate, resulting in less driving force of the fuel crossover to the cathode catalyst layer and higher open circuit voltage. Therefore, as SBS flow field was employed, whose channel width is stepwise broadened from 0.8 mm to 1.3 mm, the formic acid concentration in Figure 7b is second highest, corresponding to second highest open circuit voltage in Figure 6. As the upstream channel width of SBS fluid flow field is narrower than that of S(1.3) flow field, the pressure of the fuel flow in SBS is higher than that in S(1.3), causing a favorable condition for formic acid to penetrate into the anode catalyst layer. In addition to the slight fuel crossover of the DFAMFC possessing SBS fluid flow field, the air-breathing DFAMFC possessing SBS fluid flow field generates highest peak power density in Figure  6. As the upstream channel width of SBS fluid flow field is narrower than that of S(1.3) flow field, the pressure of the fuel flow in SBS is higher than that in S(1.3), causing a favorable condition for formic acid to penetrate into the anode catalyst layer. In addition to the slight fuel crossover of the DFAMFC possessing SBS fluid flow field, the air-breathing DFAMFC possessing SBS fluid flow field generates highest peak power density in Figure 6.
Following the optimal fluid flow field configuration presented in Figure 6, the numerical simulation was further performed for air-feeding DFAMFCs with SBS flow field for both fuel and electrolyte incorporating various air flow fields at airflow rates of 200 sccm and 500 sccm to investigate the fuel cell performance. Firstly, both V-I and P-I curves of the air-feeding DFAMFCs with various air flow fields were plotted in Figure 8 along with the V-I and P-I curves of the air-breathing DFAMFC for comparison. It can be observed in Figure 8 that, compared with the performance of the air-breathing DFAMFC, the maximum power density was enhanced whatever air flow field was employed. The fairly high pressure established in the air flow field as shown in Figure 9d-f might possibly not only resist the fuel invasion to the cathode catalyst layer, but also cause the liquid electrolyte to diffuse in a more uniform way within the GDEs because of the counterpart in pressure distribution on both sides of the electrolyte flow field. Therefore, the air-feeding DFAMFCs yielded higher open circuit voltage and lower ohmic overpotential than the air-breathing one. In addition, the effective supply of oxygen to the cathode of the DFAMFCs is responsible for the higher maximum current density. Besides, the slight effect of the airflow rate on the performance of the DFAMFCs having pin type air flow field can be observed in Figure 8. The DFAMFC having MS(0.8) air flow field fed at an airflow rate of 500 sccm yielded the highest peak power density than other cells in Figure 8, 12 mW/cm 2 , while the air-breathing one produced the peak power density of 10.5 mW/cm 2 .
Processes 2021, 9, x FOR PEER REVIEW 15 of 20 Following the optimal fluid flow field configuration presented in Figure 6, the numerical simulation was further performed for air-feeding DFAMFCs with SBS flow field for both fuel and electrolyte incorporating various air flow fields at airflow rates of 200 sccm and 500 sccm to investigate the fuel cell performance. Firstly, both V-I and P-I curves of the air-feeding DFAMFCs with various air flow fields were plotted in Figure 8 along with the V-I and P-I curves of the air-breathing DFAMFC for comparison. It can be observed in Figure 8 that, compared with the performance of the air-breathing DFAMFC, the maximum power density was enhanced whatever air flow field was employed. The fairly high pressure established in the air flow field as shown in Figure 9d-f might possibly not only resist the fuel invasion to the cathode catalyst layer, but also cause the liquid electrolyte to diffuse in a more uniform way within the GDEs because of the counterpart in pressure distribution on both sides of the electrolyte flow field. Therefore, the air-feeding DFAMFCs yielded higher open circuit voltage and lower ohmic overpotential than the air-breathing one. In addition, the effective supply of oxygen to the cathode of the DFAMFCs is responsible for the higher maximum current density. Besides, the slight effect of the airflow rate on the performance of the DFAMFCs having pin type air flow field can be observed in Figure 8. The DFAMFC having MS(0.8) air flow field fed at an airflow rate of 500 sccm yielded the highest peak power density than other cells in Figure 8, 12 mW/cm 2 , while the air-breathing one produced the peak power density of 10.5 mW/cm 2 .   In order to exhibit the air flow field effect on the performance of the air-feeding DFAMFCs, the oxygen concentration distribution within the cathode catalyst layer with air various flow fields, including SBS, Pin(0.8), and MS(0.8), at 0.4 V with airflow rates of 200 sccm and 500 sccm is plotted in Figure 10. It can be seen in Figure 10 that as the airflow rate was increased, the oxygen concentration within the catalyst layer of each DFAMFC become high. This phenomenon is especially obvious on the plane of the catalyst layer in contact with the electrolyte flow field of the DFAMFCs. However, even though the airflow rate is increased, both oxygen and electrolyte seems to maintain uniform distribution within the entire cathode catalyst layer of the DFAMFC having Pin(0.8) air flow field, which could be one of the possible reasons causing the slight increase in the peak power density of the DFAMFC having Pin(0.8) air flow field as observed in Figure 8. In order to exhibit the air flow field effect on the performance of the air-feeding DFAMFCs, the oxygen concentration distribution within the cathode catalyst layer with air various flow fields, including SBS, Pin(0.8), and MS(0.8), at 0.4 V with airflow rates of 200 sccm and 500 sccm is plotted in Figure 10. It can be seen in Figure 10 that as the airflow rate was increased, the oxygen concentration within the catalyst layer of each DFAMFC become high. This phenomenon is especially obvious on the plane of the catalyst layer in contact with the electrolyte flow field of the DFAMFCs. However, even though the airflow rate is increased, both oxygen and electrolyte seems to maintain uniform distribution within the entire cathode catalyst layer of the DFAMFC having Pin(0.8) air flow field, which could be one of the possible reasons causing the slight increase in the peak power density of the DFAMFC having Pin(0.8) air flow field as observed in Figure 8. Since an air pump has to be used to deliver the specific airflow rate in the air flow field for an air-feeding DFAMFC, the air pump consumes additional power. The power consumption of the air pump, Wp,air, can be estimated using the following equation, where p Δ (Pa) is the air pressure drop between the inlet and outlet of the air flow field, Qa (m 3 /s) is the volumetric flow rate of the air, and AR (cm 2 ) is the electrode surface area. In order to evaluate the performance of the DFAMFCs having different air flow field designs in terms of the net power output, the net power density output of the air-feeding DFAMFCs is shown in Figure 11. The net power density output is defined as the difference between the power density produced by the air-feeding DFAMFCs and the power consumption of the air pump as follows,  Since an air pump has to be used to deliver the specific airflow rate in the air flow field for an air-feeding DFAMFC, the air pump consumes additional power. The power consumption of the air pump, W p,air , can be estimated using the following equation, where ∆p (Pa) is the air pressure drop between the inlet and outlet of the air flow field, Q a (m 3 /s) is the volumetric flow rate of the air, and A R (cm 2 ) is the electrode surface area. In order to evaluate the performance of the DFAMFCs having different air flow field designs in terms of the net power output, the net power density output of the air-feeding DFAMFCs is shown in Figure 11. The net power density output is defined as the difference between the power density produced by the air-feeding DFAMFCs and the power consumption of the air pump as follows, W net = W cell − W p,air