Mechanical properties,corrosion resistance,and rubber pad forming of cold differential speed-rolled pure titanium for bipolar plates of proton-exchange membrane fuel cells

Due to problems such as pores on surface-treated coatings, the corrosion resistance of pure titanium bipolar plates for proton-exchange membrane fuel cells can be further improved by increasing the corrosion resistance of pure titanium by using differential speed-rolling (DSR); however, these materials have not yet reached the standard requirements of bipolar plates (corrosion current density i_corr＜10³ nA·cm^?2). In this work, the corrosion resistance of pure titanium was improved by optimizing the DSR process while the strength was maintained. The best corrosion resistance of the DSR pure titanium was achieved when the roller speed ratio was 2, while i_corr was 429 nA·cm^?2 in a solution of 0.5 M H₂SO₄ and 2 mg/L HF at room temperature. The formability of the DSR pure titanium for bipolar plates was verified. The optimal holding pressure range was 6.8–7.0 kN.     

Investigation of La0.6Sr0.4Co1-xNixO3-δ (x=0, 0.2, 0.4, 0.6, 0.8) catalysts on solid oxide fuel cells anode for biogas dry reforming

The introduction of catalyst on anode of solid oxide fuel cell (SOFC) has been an effective way to alleviate the carbon deposition when utilizing biogas as the fuel. A series of La_0.6Sr_0.4Co_1-xNi_xO_3-δ (x = 0, 0.2, 0.4, 0.6, 0.8) oxides are synthesized by sol-gel method and used as catalysts precursors for biogas dry reforming. The phase structure of La_0.6Sr_0.4Co_1-xNi_xO_3-δ oxides before and after reduction are characterized by X-ray diffraction (XRD). The texture properties, carbon deposition, CH₄ and CO₂ conversion rate of La_0.6Sr_0.4Co_1-xNi_xO_3-δ catalysts are evaluated and compared. The peak power density of 739 mW cm^?2 is obtained by a commercial SOFC with La_0.6Sr_0.4Co_0.4Ni_0.6O_3-δ catalyst at 850 °C when using a mixture of CH₄: CO₂ = 2:1 as fuel. This shows a great improvement from the cell without catalyst for internal dry reforming, which is attributed to the formation of NiCo alloy active species after reduction in H₂ atmosphere. The results indicate the benefits of inhibiting the carbon deposition on Ni-based anode through introducing the La_0.6Sr_0.4Co_0.4Ni_0.6O_3-δ catalyst precursor. Additionally, the dry reforming technology will also help to convert part of the exhaust heat into chemical energy and improve the efficiency of SOFC system with biogas fuel.     

Effect of micro-cracks on the in-plane electronic conductivity of proton exchange membrane fuel cell catalyst layers based on lattice Boltzmann method

Micro-cracks commonly occur on the catalyst layers (CLs) during the manufacturing of catalyst coated membranes (CCMs). However, the crack shape parameters effect on CLs in-plane (IP) electronic conductivity λ_s is not clear. In this work, the relationship between crack parameters and the λ_s is obtained based on the two-dimensional (2D) multiple-relaxation time (MRT) lattice Boltzmann method (LBM). The LBM numerical model is validated by the normalized λ_s experiment applied on three different home-made cracked CLs, and the parameter study focus on crack width, length, quantity and phase angle are carried out. The results show that the decrease of λ_s has different sensitivity |k| to the parameters above. The crack width has little effect on λ_s decrease, and the |k_w| is 0.038. However, crack arm length and quantity show more significant impact, which |k_l| and |k_N| are 0.753 and 0.725, respectively. The CLs with different crack propagation directions show significant anisotropy on λ_s, and a 53.53% decrease in λ_s is observed between 0° and 90° crack phase angle change. To manufacture a high electronic conductivity CL, crack initiation and migration mitigation are highly encouraged.     

A novel pre-dilution,swirling jet diffuser to enhance effluent mixing: Hydrodynamics and dilution performance

Diffusers are widely-used to quickly dilute effluents in receiving water bodies. This study proposed a novel diffuser that pre-mixes effluent with ambient water before discharging and that uses the swirling jet to further enhance near-field dilution. The nozzle of the diffuser was examined in two ambient flow conditions: co-flow and counter-flow that are commonly-met in the environment such as oceans due to tidal effect. Physical experiments were first conducted in co-flow on its dilution performance and hydrodynamics, using heated water as the effluent. A 3-D CFD model was developed and calibrated the co-flow scenarios, and then used to investigate the diffuser in counter-flow. The results showed that the nozzle can effectively reduce the maximum temperature rise of the effluent by about 50 % before discharging. The swirling jet from the outlet has a larger shear area, half-width and entrainment rate, enabling the effluent to be rapidly diluted to a minimum of around 10 times at x/D = 6 in co-flow, whereas the dilution for conventional nozzles is about 1 because of the potential core. The flow amplification ratio (α) decreases gradually with increasing velocity ratio in co-flow but increases with increasing velocity ratio in counter-flow. The counter-flow reduces the water drawn into the device; however, the pre-dilution effect at the outlet remains stable. The near-field dilution in counter-flow was significantly enhanced than that in co-flow. Environmental regulations at outfalls and mixing zones can be more easily met using this novel diffuser.     

Three-dimensional and anisotropic numerical analysis of a PEM fuel cell

In this study, the effects of cell temperature and relative humidity on charge transport parameters are numerically analyzed. In order to perform this analysis, three-dimensional and anisotropic numerical models are developed. The numerical models are integrated into the experimental values for anisotropic electrical conductivities, as depending on cell temperature and relative humidity, that were obtained from our previous study. The achieved results indicate that the values of current densities in the in-plane direction increase with increasing cell temperature and relative humidity, while the current densities reach a maximum in the rib regions for both the numerical model at the through-plane direction. The behaviors of electrolyte potentials are similar with changes in the cell temperature and relative humidity. In addition, the cathode electrical potentials in both the in-plane direction and through-plane direction do not change to a considerable amount with increasing cell temperature and relative humidity.     

Simulation of hydrogen and coal gas fueled flat-tubular solid oxide fuel cell (FT-SOFC)

Solid oxide fuel cells (SOFCs) are considered an important technology in terms of high efficiency and clean energy generation. Flat-tubular solid oxide fuel cell (FT-SOFC) which is a combination of tubular and planar cell geometries stands out with its performance values and low costs. In this study, the performance of an FT-SOFC is analyzed numerically by using finite element method-based design as a result of changing parameters by using different fuels which are pure hydrogen and coal gas with various proportions of CO. In addition, cell performance values for different temperatures were analyzed and interpreted. Analyzes have been performed by using COMSOL Multiphysics software. The rates of CO composition used are 10%, 20%, and 40%, respectively. In addition, the air was used as the oxidizer in all cases. The cell voltage and average cell power of the FT-SOFC were examined under the 800 °C operating condition. The maximum power value and current density value were obtained as 710 W/m² and 1420 A/m² for the flat-tubular cell, respectively. As a result of the study, it was observed that the maximum cell power densities increased with increasing temperature. Analysis results showed that FT-SOFCs have suitable properties for different fuel usage and different operating temperatures. High-performance values and design features in different operating conditions are expected to make FT-SOFC the focus of many studies in the future.     

Motion trajectory prediction model of hydrogen leak and diffusion in a stable thermally-stratified environment

The motion trajectory of hydrogen leakage is an essential safe issue for the application of hydrogen energy. A dimensionless fast-running motion trajectory prediction model is proposed to predict the dispersion characteristics of the buoyant jet of hydrogen leakage for the accident. The impact of different leakage angles, leakage velocity and thermal stratification of ambient air on hydrogen leakage behavior was analyzed. The new developed model was verified by experimental results in literatures. Leakage hydrogen can flow upwards freely in a uniform environment. However, it shows an oscillating trajectory at a certain height in a thermally stratified environment, which is so called “locking phenomenon”. The trajectory of hydrogen leakage is upward and hydrogen gathers at the top of the space to form stratification in a uniform environment, while the hydrogen leakage shows an oscillating trajectory at a certain height in a thermal stratification environment. With the increase of Froude number Fr, it shows that the stable height and maximum height of the leakage airflow have a trend of rising first and then falling in a thermally stratified environment. The findings are expected to give guidance in real-world situations, for example, a larger Fr value and a larger temperature gradient can lead to a decrease in the stable height in the thermally stratified environment. It is found that the fitting of the stable height with different temperature gradients satisfies the power function relationship. This work is expected to be helpful for reducing hydrogen leakage accumulation and explosion risk.     

Numerical simulation of water and heat transport in the cathode channel of a PEM fuel cell

Cathode channel of a PEM fuel cell is the critical domain for the transport of water and heat. In this study, a mathematical model of water and heat transport in the cathode channel is established by considering two-phase flow of water and air as well as the phase change between water and vapor. The transport process of the species of air is governed by the convection-diffusion equation. The VOSET (coupled volume-of-fluid and level set method) method is used to track the interface between air and water, and the phase equilibrium method of water and vapor is employed to calculate the mass transfer rate on the two-phase interface. The present model is validated against the results in the literature, then applied to investigate the characteristics of two-phase flow and heat transfer in the cathode channel. The results indicate that in the inlet section, water droplets experience three evolution stages: the growing stage, the coalescence stage and the generation stage of dispersed water drops. However, in the middle and outlet sections of the channel, there are only two stages: the growth of water droplets, and the formation of a water film. The mass transfer rate of phase change in the inlet section of the channel varies over time, exhibiting an initial increase, a decrease followed, and a stabilization finally, with the maximum and stable values of 1.78 × 10^?4 kg/s and 1.52 × 10^?4 kg/s for Part 1, respectively. In the middle and outlet sections, the mass transfer rate increase firstly and then keeps stable gradually. Furthermore, regarding the distribution of the temperature and vapor mass fraction in the channel, near the upper surface of the channel, the temperature and vapor mass fraction first change slightly (x < 0.03 m) and then rapidly decrease with fluctuations (x > 0.03 m). In the middle of the channel, the temperature and vapor mass fraction slowly decrease with fluctuation.     

Analysis and proposition of limiting current density measurement protocol

Limiting current density at different temperatures, backpressures, and balance gases can be used to separate molecular diffusion resistance, Knudsen diffusion resistance and local transport resistance of membrane electrode assembly (MEA). However, the measurement of limiting current density has no unified protocol. The diverse choices in the literature, either in the control of current or voltage or in the atmosphere like relative humidity and O₂ concentrations, make it difficult to compare the results and identify the true bottleneck hindering the mass transport. In this work, the current-voltage curves obtained by current scanning/stepping and voltage scanning/stepping methods under dilute O₂ of different concentrations and a wide range of relative humidity were measured and analyzed systematically. It is found that the voltage stepping method is superior to the other three ways of control for the reliable determination of the limiting current density. Aided with simultaneous electrochemical impedance spectroscopy measurement, the limiting current density can be determined with pinpoint accuracy. When the limiting current density is just used to qualitatively evaluate different MEA, the voltage scanning method can be used instead for its high time efficiency. The selection of the atmosphere also plays an important role in suppressing the distortion from excessive water and reducing the spurious contribution from proton conduction resistance. It is found that O₂ concentrations at 0.5 vol% and relative humidity at 90% can give the best estimation of O₂ transport resistance in membrane electrode assembly.