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

Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media

This paper introduces a new discrete fracture model accounting for non-isothermal compositional multiphase Darcy flows and complex networks of fractures with intersecting, immersed and non immersed fractures. The so called hybrid-dimensional model using a 2D model in the fractures coupled with a 3D model in the matrix is first derived rigorously starting from the equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully implicit time integration combined with the Vertex Approximate Gradient (VAG) finite volume scheme which is adapted to polyhedral meshes and anisotropic heterogeneous media. The fully coupled systems are assembled and solved in parallel using the Single Program Multiple Data (SPMD) paradigm with one layer of ghost cells. This strategy allows for a local assembly of the discrete systems. An efficient preconditioner is implemented to solve the linear systems at each time step and each Newton type iteration of the simulation. The numerical efficiency of our approach is assessed on different meshes, fracture networks, and physical settings in terms of parallel scalability, nonlinear convergence and linear convergence

3D particle tracking velocimetry using dynamic discrete tomography

Particle tracking velocimetry in 3D is becoming an increasingly important imaging tool in the study of fluid dynamics, combustion as well as plasmas. We introduce a dynamic discrete tomography algorithm for reconstructing particle trajectories from projections. The algorithm is efficient for data from two projection directions and exact in the sense that it finds a solution consistent with the experimental data. Non-uniqueness of solutions can be detected and solutions can be tracked individually

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified