Two-phase flow dynamics in the gas diffusion layer of proton exchange membrane fuel cells: Volume of fluid modeling and comparison with experiment

This paper proposes a three-dimensional (3D) volume of fluid (VOF) study to investigate two-phase flow in the gas diffusion layer (GDL) of proton exchange membrane (PEM) fuel cells and liquid water distribution. A stochastic model was adopted to reconstruct the 3D microstructures of Toray carbon papers and incorporate the experimentally-determined varying porosity. The VOF predictions were compared with the water profiles obtained by the X-ray tomographic microscopy (XTM) and the Leverett correlation. It was found local water profiles are similar in the sample’s sub-regions under the pressure difference p = 1000 Pa between the two GDL surfaces, but may vary significantly under p = 6000 Pa. The water-air interfaces inside the GDL structure were presented to show water distribution and breakthrough

Investigating the in-/through-plane effective diffusivities of dry and partially-saturated gas diffusion layers

In this study, the effective oxygen diffusivity in the dry or partially-saturated gas diffusion layer (GDL) is numerically investigated by an oxygen diffusion model in GDLs reconstructed by a stochastic method. The predicted effective diffusivity in dry GDLs is compared with various diffusivity models from literatures. Reasonable agreements with other models were obtained. The effect of the PTFE loading in the dry Toray carbon paper is also investigated and compared with recent experimental data. It is found that the effective diffusivity becomes lower under higher PTFE loading due to the decreased pore volume, as expected. The relative effective oxygen diffusivity in partially-saturated GDLs is calculated using the two-phase volume of fluid (VOF) model and an oxygen diffusion model. The effects of different local water profiles and porosity distribution on the effective oxygen diffusivity in both the through-plane (TP) and in-plane (IP) directions are investigated and compared with a lattice Boltzmann model and experimental data. The present results are in good agreement with other studies. It is found that local water profile has significant impacts on the effective diffusivity in partially-saturated GDLs and the diffusivity in the TP direction is more sensitive to the water distribution than the IP direction

Topological Classification of Crystalline Insulators with Point Group Symmetry

We show that in crystalline insulators point group symmetry alone gives rise to a topological classification based on the quantization of electric polarization. Using C3 rotational symmetry as an example, we first prove that the polarization is quantized and can only take three inequivalent values. Therefore, a Z3 topological classification exists. A concrete tight-binding model is derived to demonstrate the Z3 topological phase transition. Using first-principles calculations, we identify graphene on BN substrate as a possible candidate to realize the Z3 topological states. To complete our analysis we extend the classification of band structures to all 17 two-dimensional space groups. This work will contribute to a complete theory of symmetry conserved topological phases and also elucidate topological properties of graphene like systems

Evolution of Threats in the Global Risk Network

With a steadily growing population and rapid advancements in technology, the global economy is increasing in size and complexity. This growth exacerbates global vulnerabilities and may lead to unforeseen consequences such as global pandemics fueled by air travel, cyberspace attacks, and cascading failures caused by the weakest link in a supply chain. Hence, a quantitative understanding of the mechanisms driving global network vulnerabilities is urgently needed. Developing methods for efficiently monitoring evolution of the global economy is essential to such understanding. Each year the World Economic Forum publishes an authoritative report on the state of the global economy and identifies risks that are likely to be active, impactful or contagious. Using a Cascading Alternating Renewal Process approach to model the dynamics of the global risk network, we are able to answer critical questions regarding the evolution of this network. To fully trace the evolution of the network we analyze the asymptotic state of risks (risk levels which would be reached in the long term if the risks were left unabated) given a snapshot in time, this elucidates the various challenges faced by the world community at each point in time. We also investigate the influence exerted by each risk on others. Results presented here are obtained through either quantitative analysis or computational simulations.Comment: 27 pages, 15 figure

Numerical simulation of two-phase cross flow in the gas diffusion layer microstructure of proton exchange membrane fuel cells

The cross flow in the under-land gas diffusion layer (GDL) between 2 adjacent channels plays an important role on water transport in proton exchange membrane fuel cell. A 3-dimensional (3D) two-phase model that is based on volume of fluid is developed to study the liquid water-air cross flow within the GDL between 2 adjacent channels. By considering the detailed GDL microstructures, various types of air-water cross flows are investigated by 3D numerical simulation. Liquid water at 4 locations is studied, including droplets at the GDL surface and liquid at the GDL-catalyst layer interface. It is found that the water droplet at the higher-pressure channel corner is easier to be removed by cross flow compared with droplets at other locations. Large pressure difference Δp facilitates the faster water removal from the higher-pressure channel. The contact angle of the GDL fiber is the key parameter that determines the cross flow of the droplet in the higher-pressure channel. It is observed that the droplet in the higher-pressure channel is difficult to flow through the hydrophobic GDL. Numerical simulations are also performed to investigate the water emerging process from different pores of the GDL bottom. It is found that the amount of liquid water removed by cross flow mainly depends on the pore's location, and the water under the land is removed entirely into the lower-pressure channel by cross flow