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.     

Process Intensification of Living Cationic Polymerization of Isobutylene by a Flow Reaction System

Increasing the reaction temperature of the living cationic polymerization of isobutylene is crucial for industrial production due to the cost of refrigeration. The reaction temperature increase was achieved with an accelerated reaction rate using a flow reaction system. The polymerization conditions, including the flow reactor design, were based on the results of kinetic studies. Utilizing a milli‐scale flow reactor, polyisobutylene, which has a narrow molecular weight distribution, was obtained within a considerably short residence time at a high temperature. Furthermore, it was confirmed that the value of M_w/M_n correlates with the product of the Reynolds number and the angle of collision.     

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.     

Hybrid hydrate processes for CO2/H2 mixture purification: A techno-economic analysis

In this work, hydrate based separation technique was combined with membrane separation and amine-absorption separation technologies to design hybrid processes for separation of CO₂/H₂ mixture. Hybrid processes are designed in the presence of different types of hydrate promoters. The conceptual processes have been developed using Aspen HYSYS. Proposed processes were simulated at different flow rates for the feed stream. A comprehensive cost model was developed for economic analysis of novel processes proposed in this study. Based on the results from process simulation and equipment sizing, the amount of total energy consumption, fixed cost, variable cost, and total cost were calculated per unit weight of captured CO₂ for various flow rates of feed stream and in the presence of different hydrate promoters. Results showed that combination of hydrate formation separation technique with membrane separation technology results in a CO₂ capture process with lowest energy consumption and total cost per unit weight of captured CO₂. As split fraction and heat of hydrate formation increases, the share of hydrate formation section in total energy consumption increases. When TBAB is applied as hydrate promoter, due to its higher hydrate separation efficiency, more amount of CO₂ is captured in hydrate formation section and consequently the total cost for process decreases considerably. Hybrid hydrate-membrane process in the presence of TBAB as hydrate promoter with 29.47 US$/ton CO₂ total cost is the best scheme for hybrid hydrate CO₂ capture process. Total cost for this process is lower than total cost for single MDEA-based absorption process as the mature technology for CO₂ capture.     

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.     

Theoretical framework for formulating MIP scheduling models with multiple and non-uniform discrete-time grids

We present a framework for the formulation of MIP scheduling models based on multiple and nonuniform discrete time grids. In a previous work we showed that it is possible to use different (possibly non-uniform) time grids for each task, unit, and material. Here, we generalize these ideas to account for general resources, and a range of processing characteristics such as limited intermediate storage and changeovers. Each resource has its own grid based on resource consumption and availability allowing resource constraints to be modeled more accurately without increasing the number of binary variables. We develop algorithms to define the unit-, task-, material-, and resource-specific grids directly from problem data. Importantly, we prove that the multi-grid formulation is able to find a schedule with the same optimal objective as the discrete-time single-grid model with an arbitrarily fine grid. The proposed framework leads to the formulation of models with reduced number of binary variables and constraints, which are able to find good solutions faster than existing models.     

A semicontinuous approach for heterogeneous azeotropic distillation processes

The separation of azeotropes has substantial energy and investment costs, and the available methods require high capital costs for reconstruction of process plants. As an alternative, a semicontinuous configuration that utilizes an existing plant with minor modifications has been explored. In this paper, a semicontinuous, heterogeneous azeotropic distillation process is proposed and acetic acid dehydration process is used as a case study. To carry out the simulation work, Aspen HYSYS^® simulation software is used along with MATLAB^® and an interface program to handle the mode-transition of the semicontinuous process. Sensitivity analyses on operating parameters are performed to identify the process limits. Comparisons are made to conventional heterogeneous azeotropic distillation, and dividing-wall distillation column on the annual cost. The results proved that the semicontinuous system is the best setup in terms of total annual costs and energy requirements.