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.     

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.     

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.     

Numerical study on laminar burning velocity of ammonia flame with methanol addition

The combustion characteristics of ammonia/methanol mixtures were investigated numerically in this study. Methanol has a dramatic promotive effect on the laminar burning velocity (LBV) of ammonia. Three mechanisms from literature and another four self-developed mechanisms constructed in this study were evaluated using the measured laminar burning velocities of ammonia/methanol mixtures from Wang et al. (Combust.Flame. 2021). Generally, none of the selected mechanisms can precisely predict the measured laminar burning velocities at all conditions. Aiming to develop a simplified and reliable mechanism for ammonia/methanol mixtures, the constructed mechanism utilized NUI Galway mechanism (Combust.Flame. 2016) as methanol sub-mechanism and the Otomo mechanism (Int. J. Hydrogen. Energy. 2018) as ammonia sub-mechanism was optimized and reduced. The reduced mechanism entitled ‘DNO-NH3’, can accurately reproduce the measured laminar burning velocities of ammonia/methanol mixtures under all conditions. A reaction path analysis of the ammonia/methanol mixtures based on the DNO-NH3 mechanism shows that methanol is not directly involved in ammonia oxidation, instead, the produced methyl radicals from methanol oxidization contribute to the dehydrogenation of ammonia. Besides, NO_x emission analysis demonstrates that 60% methanol addition results in the highest NO_x emissions. The most important reactions dominating the NO_x consumption and production are identified in this study.     

Mechanical Properties of Single-Layer Diamond Reinforced Poly(vinyl alcohol) Nanocomposites through Atomistic Simulation

Low-dimensional carbon nanostructures are ideal nanofillers to reinforce the mechanical performance of polymer nanocomposites due to their excellent mechanical properties. Through molecular dynamics simulations, the mechanical performance of poly(vinyl alchohol) (PVA) nanocomposites reinforced with a single-layer diamond – diamane is investigated. It is found the PVA/diamane exhibits similar interfacial strengths and pull-out characteristics with the PVA/bilayer-graphene counterpart. Specifically, when the nanofiller is fully embedded in the nanocomposite, it is unable to deform simultaneously with the PVA matrix due to the weak interfacial load transfer efficiency, thus the enhancement effect is not significant. In comparison, diamane can effectively promote the tensile properties of the nanocomposite when it has a laminated structure as it deforms simultaneously with the matrix. With this configuration, the interlayer sp³ bonds endows diamane with a much higher resistance under compression and shear tests, thus the nanocomposite can reach very high compressive and shear stress. Overall, enhancement on the mechanical interlocking at the interface as triggered by surface functionalization is only effective for the fully embedded nanofiller. This work provides a fundamental understanding of the mechanical properties of PVA nanocomposites reinforced by diamane, which can shed lights on the design and preparation of next generation high-performance nanocomposites.     

In Silico Investigation of Potential Applications of Gamma Carbonic Anhydrases as Catalysts of CO2 Biomineralization Processes: A Visit to the Thermophilic Bacteria Persephonella hydrogeniphila,Persephonella marina,Thermosulfidibacter takaii,and Thermus thermophilus

Carbonic anhydrases (CAs) have been identified as ideal catalysts for CO₂ sequestration. Here, we report the sequence and structural analyses as well as the molecular dynamics (MD) simulations of four γ-CAs from thermophilic bacteria. Three of these, Persephonella marina, Persephonella hydrogeniphila, and Thermosulfidibacter takaii originate from hydrothermal vents and one, Thermus thermophilus HB8, from hot springs. Protein sequences were retrieved and aligned with previously characterized γ-CAs, revealing differences in the catalytic pocket residues. Further analysis of the structures following homology modeling revealed a hydrophobic patch in the catalytic pocket, presumed important for CO₂ binding. Monitoring of proton shuttling residue His69 (P. marina γ-CA numbering) during MD simulations of P. hydrogeniphila and P. marina’s γ-CAs (γ-PhCA and γ-PmCA), showed a different behavior to that observed in the γ-CA of Escherichia coli, which periodically coordinates Zn²⁺. This work also involved the search for hotspot residues that contribute to interface stability. Some of these residues were further identified as key in protein communication via betweenness centrality metric of dynamic residue network analysis. T. takaii’s γ-CA showed marginally lower thermostability compared to the other three γ-CA proteins with an increase in conformations visited at high temperatures being observed. Hydrogen bond analysis revealed important interactions, some unique and others common in all γ-CAs, which contribute to interface formation and thermostability. The seemingly thermostable γ-CA from T. thermophilus strangely showed increased unsynchronized residue motions at 423 K. γ-PhCA and γ-PmCA were, however, preliminarily considered suitable as prospective thermostable CO₂ sequestration agents.     

Nonlinear dynamic behaviors and control based on simulation of high-purity heat integrated air separation column

In this paper, the dynamic behaviors on the basis of simulation for high-purity heat integrated air separation column (HIASC) are studied. A nonlinear generic model control (GMC) scheme is proposed based on the nonlinear behavior analyses of a HIASC process, and an adaptive generic model control (AGMC) scheme is further presented to correct the model parameters online. Related internal model control (IMC) scheme and multi-loop PID (M-PID) scheme are also developed as the comparative base. The comparative researches are carried out among these linear and nonlinear control schemes in detail. The simulation research results show that the proposed AGMC schemes present advantages in both servo control and regulatory control for the high-purity HIASC.     

The effect of charge on the dihydrogen storage capacity of Sc2–C6H6

The effect of charge on the dihydrogen storage capacity of Sc₂–C₆H₆ has been investigated at B3LYP-D3/6-311G(d,p) level. The neutral system Sc₂–C₆H₆ can store 8H₂ with gravimetric density of 8.76 wt %, and one H₂ dissociates and bonds atomically on the scandium atom. The adsorption of 8H₂ on Sc₂–C₆H₆ is energetically favorable below 155 K. The atom-centered density matrix propagation (ADMP) molecular dynamics simulations show that Sc₂–C₆H₆ can adsorb 3H₂ within 1000 fs at 300K. Compared with Sc₂–C₆H₆, the charged systems can adsorb more hydrogen molecules with higher gravimetric density, and all the H₂ are adsorbed in the molecular form. The gravimetric densities of Sc₂–C₆H₆⁺ and Sc₂–C₆H₆²⁺ are 9.75 and 10.71 wt%. Moreover, the maximum adsorption of charged systems are favorable in wider temperature range. Most importantly, the ADMP-MD simulations indicate that Sc₂–C₆H₆²⁺ can adsorb 6 hydrogen molecules within 1000 fs at 300K. It can be found that the gravimetric density (6.72 wt%) of Sc₂–C₆H₆²⁺ still exceeds the target of US Department of Energy (DOE) under ambient conditions.     

Revisiting the Parallel Strategy for DOACROSS Loops

Journal of Computer Science and Technology - DOACROSS loops are significant parts in many important scientific and engineering applications, which are generally exploited pipeline/wave-front...     

Mixture ratio modeling of dynamic systems

Any knowledge extraction relies (possibly implicitly) on a hypothesis about the modelled-data dependence. The extracted knowledge ultimately serves to a decision-making (DM). DM always faces uncertainty and this makes probabilistic modelling adequate. The inspected black-box modeling deals with “universal” approximators of the relevant probabilistic model. Finite mixtures with components in the exponential family are often exploited. Their attractiveness stems from their flexibility, the cluster interpretability of components and the existence of algorithms for processing high-dimensional data streams. They are even used in dynamic cases with mutually dependent data records while regression and auto-regression mixture components serve to the dependence modeling. These dynamic models, however, mostly assume data-independent component weights, that is, memoryless transitions between dynamic mixture components. Such mixtures are not universal approximators of dynamic probabilistic models. Formally, this follows from the fact that the set of finite probabilistic mixtures is not closed with respect to the conditioning, which is the key estimation and predictive operation. The paper overcomes this drawback by using ratios of finite mixtures as universally approximating dynamic parametric models. The paper motivates them, elaborates their approximate Bayesian recursive estimation and reveals their application potential.