Investigation of liquid water heterogeneities in large area proton exchange membrane fuel cells using a Darcy two-phase flow model in a multiphysics code

Proper management of the liquid water and heat produced in proton exchange membrane (PEM) fuel cells remains crucial to increase both its performance and durability. In this study, a two-phase flow and multicomponent model, called two-fluid model, is developed in the commercial COMSOL Multiphysics® software to investigate the liquid water heterogeneities in large area PEM fuel cells, considering the real flow fields in the bipolar plate. A macroscopic pseudo-3D multi-layers approach has been chosen and generalized Darcy's relation is used both in the membrane-electrode assembly (MEA) and in the channel. The model considers two-phase flow and gas convection and diffusion coupled with electrochemistry and water transport through the membrane. The numerical results are compared to one-fluid model results and liquid water measurements obtained by neutron imaging for several operating conditions. Finally, according to the good agreement between the two-fluid and experimentation results, the numerical water distribution is examined in each component of the cell, exhibiting very heterogeneous water thickness over the cell surface.     

Fabrication of crescent-shaped ceramic microparticles based on single emulsion microfluidics

Ceramic microparticles have great potentials in various fields such as materials engineering, biotechnology, microelectromechanical systems, etc. Morphology of the microparticle performs an important role on their application. To date, it remains difficult to find an effective and controllable way for fabricating nonspherical ceramic microparticles with 3D features. This work demonstrates a method that combines UV light lithography and single emulsion opaque-droplet-templated microfluidic molding to prepare the crescent-shaped ceramic microparticles. By tailoring the intensity of UV light and flow rate of fluid, the shapes of microparticles are accordingly tuned. Therefore, varieties of crescent-shaped microparticles and their variations have been fabricated. After sintering, the crescent-shaped alumina ceramic microparticles were obtained. Benefitting from the light absorption and scattering behavior of most ceramic nanoparticles, this system can serve as a general platform to produce crescent-shaped microparticles made from different materials, and hold great potentials for applications in microrobotics, structural materials in MEMS, and biotechnology.     

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.     

Numerical investigations on flashback dynamics of premixed methane-hydrogen-air laminar flames

Injecting hydrogen into the natural gas network to reduce CO₂ emissions in the EU residential sector is considered a critical element of the zero CO₂ emissions target for 2050. Burning natural gas and hydrogen mixtures has potential risks, the main one being the flame flashback phenomenon that could occur in home appliances using premixed laminar burners. In the present study, two-dimensional transient computations of laminar CH₄ + air and CH₄ + H₂ + air flames are performed with the open-source CFD code OpenFOAM. A finite rate chemistry based solver is used to compute reaction rates and the laminar reacting flow. Starting from a flame stabilized at the rim of a cylindrical tube burner, the inlet bulk velocity of the premixture is gradually reduced to observe flashback. The results of the present work concern the effects of wall temperature and hydrogen addition on the flashback propensity of laminar premixed methane-hydrogen-air flames. Complete sequences of flame dynamics with gradual increases of premixture velocity are investigated. At the flame flashback velocities, strong oscillations at the flame leading edge emerge, causing broken flame symmetry and finally flame flashback. The numerical results reveal that flashback tendency increase with increasing wall temperature and hydrogen addition rate.     

Dense and core-rim structured B4C-TiB2 ceramics with Mo-Co-WC additive

B₄C-TiB₂ ceramics (TiB₂ ranging 5~70 vol%) with Mo-Co-WC as the sintering additive were prepared by spark plasma sintering. In comparison with B₄C-TiB₂ without additive, the enhanced densification was evident in the sintered specimen with Mo-Co-WC additive. Core-rim structured grain was observed around TiB₂ grains. The interface of the rim between TiB₂ and B₄C phases demonstrated different feature: the inner borderline of the rim exhibited a smooth feature, whereas a sharp curved grain boundary was observed between the rim and the B₄C grain. The formation mechanism is discussed: the epitaxial growth of (Ti,Mo,W)B₂ rim around the TiB₂ core may occur as a result of the solid solution and dissolution-precipitation between TiB₂ phase and the sintering additive. It was revealed that the fracture toughness increased as the content of TiB₂ content increased, alongside the decreased hardness. B₄C-30 vol% TiB₂ specimen demonstrated the optimal combination of mechanical properties, reaching Vickers hardness of 24.3 GPa and fracture toughness of 3.33 MPa·m^1/2.     

基于Abaqus的双裂纹不同间距对裂纹扩展速率的研究

扩展有限元法是近年经过大量运用的,在传统有限元的范围中求解不连续问题一种有效计算方法,它是基于单位分解的思想,在计算不连续问题时加入跳跃函数。以ABAQUS为平台,基于扩展有限元方法 (XFEM),以含双穿透型裂纹的有限宽板受横向拉伸载荷为力学模型,建立相应的裂纹尖端应力的有限元模型,研究焊接接头区域不同间距双裂纹相互作用对裂纹扩展速率的影响。结果表明:双裂纹间距的的大小并没有对裂纹的扩展速率产生影响。     

Study on ice-melting performance of gradient gas diffusion layer in proton exchange membrane fuel cell

In this study, a three-dimensional model was established using the lattice Boltzmann method (LBM) to study the internal ice melting process of the gas diffusion layer (GDL) of the proton exchange membrane fuel cell (PEMFC). The single-point second-order curved boundary condition was adopted. The effects of GDL carbon fiber number, growth slope of the number of carbon fibers and carbon fiber diameter on ice melting were studied. The results were revealed that the temperature in the middle and lower part of the gradient distribution GDL is significantly higher than that of the no-gradient GDL. With the increase of the growth slope of the number of carbon fiber, the temperature and melting rate gradually increase, and the position of the solid-liquid interface gradually decreases. The decrease in the number of carbon fibers has a similar effect as the increase in the growth slope of the number of carbon fibers. In addition, as the diameter of the carbon fiber increases, the position of the solid-liquid interface gradually decreases first and then increases.     

Effects of upholstery materials on the burning behavior of real‐scale residential upholstered furniture mock‐ups

Fire spread and growth on real‐scale four cushion mock‐ups of residential upholstered furniture (RUF) were investigated with the goal of identifying whether changes in five classes of materials (barrier, flexible polyurethane foam, polyester fiber wrap, upholstery fabric, and sewing thread), referred to as factors, resulted in statistically significant changes in burning behavior. A fractional factorial experimental design plus practical considerations yielded a test matrix with 20 material combinations. Experiments were repeated a minimum of two times. Measurements included fire spread rates derived from video recordings and heat release rates (HRRs). A total of 13 experimental parameters (3 based on the videos and 10 on the HRR results), referred to as responses, characterized the measurements. Statistical analyses based on Main Effects Plots (main effects) and Block Plots (main effects and factor interactions) were used. The results showed that three of the factors resulted in statistically significant effects on varying numbers of the 13 responses. The Barrier and Fabric factors had the strongest main effects with roughly comparable magnitudes. Foam was statistically significant for fewer of the responses and its overall strength was weaker than for Barrier and Fabric. No statistically significant main effects were identified for Wrap or Thread. Multiple two‐term interactions between factors were identified as being statistically significant. The Barrier*Fabric interaction resulted in the highest number of and strongest statistically significant effects. The existence of two‐term interactions means that it will be necessary to consider their effects in approaches designed to predict the burning behavior of RUF.     

A practical development of engineering simulation-assisted educational AR environments

Engineering simulations have opened several gates for today’s chemical engineers. They are powerful tools to provide technical content as physics-based numerical solvers. Augmented reality (AR) and virtual reality (VR), on the other hand, are already underway to digitize environments in many fields. The combination of AR/VR environments and simulations in engineering education has been attracting widespread interest. Literature has demonstrated a massive amount of digital educational environments in several contexts as being complementary to conventional educational methods. Nevertheless, hosting technical content produced by engineering simulations with educational AR/VR is still challenging and requires expertise from multiple disciplines throughout the technical development. Present work provides a facile and agile methodology for low-cost hardware but content-wise rich AR software development. Inspired by the Covid-19 pandemic, a case study is developed to teach chemical-engineering concepts using a liquid-soap synthesis process. Accordingly, we assess and conclude the digital development process to guide inexperienced developers for the digitalization of teaching content. The present contribution serves as an example of the power of integrating AR/VR with traditional engineering simulations for educational purposes. The digital tool developed in this work is shared in the online version.     

Organic mesh template-based laminated object manufacturing to fabricate ceramics with regular micron scaled pore structures

A new aqueous slurry-based laminated object manufacturing process for porous ceramics is proposed: firstly, an organic mesh sheet is pre-paved as a pore-forming template before slurry layer scraping; secondly, the 2D pattern is built with laser outline cutting of the dried mesh–ceramic composite layer; finally, the pore structure is formed after degreasing and sintering. Alumina parts with porosities of 51.5 %, round hole diameters of 80 ± 5 μm were fabricated using 70 wt. % solid content slurry and 100 mesh nylon net. Using an organic mesh as the framework and template not only reduces the risk of damage of the green body but also ensures the regularity, uniformity and connectivity of the micron scaled pore network. The layer-by-layer drying method avoids the delamination phenomenon and improves the paving density. The new method can realize the flexible design of the pore structure by using various organic mesh templates.