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.     

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.     

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein''s folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues'' side chains form a tight interaction network within the protein''s core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.     

Does decision quality (always) increase with the size of information samples? Some vicissitudes in applying the law of large numbers.

Adaptive decision making requires that contingencies between decision options and their relative assets be assessed accurately and quickly. The present research addresses the challenging notion that contingencies may be more visible from small than from large samples of observations. An algorithmic account for such a seemingly paradoxical effect is offered within a satisficing-choice framework. Accordingly, a choice is only made when the sample contingency describing the relative evaluation of the 2 options exceeds a critical threshold. Small samples, because of the high dispersion of their sampling distribution, facilitate above-threshold contingencies. Across a broad range of parameters, the resulting small-sample advantage in terms of hits is stronger than their disadvantage in false alarms. Computer simulations and experiments support the model predictions. The relative advantage of small samples is most apparent when information loss is low, when the threshold is high relative to the ecological contingency, and when the sampling process is self-truncated. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.