Search

The series of metal–organic frameworks M-MOF-74 gained popularity in the field of capture and separation of CO2 due to the presence of numerous, highly reactive open-metal sites. The description of effective interactions between guest molecules and open-metal sites without accounting for polarization effects is challenging but it can significantly reduce the computational cost of simulations. In this study, we propose a non-polarizable force field for CO2, and H2 adsorption in M-MOF-74 (M = Ni, Cu, Co, Fe, Mn, Zn) by scaling the Coulombic interactions of M-MOF-74 atoms, and Lennard-Jones interaction potentials between the center of mass of H2 and the open-metal centers. The presented force field is based on UFF and DREIDING parameters, characterized by high transferability and efficiency. The quantum behavior of H2 at cryogenic temperatures is considered by incorporating Feynman–Hibbs quantum corrections. To validate the force field, the experimental isotherms of CO2 at 298 K and 10−1 – 102kPa, the isotherms of H2 at 77 K and 10−5 – 102kPa, the corresponding enthalpy of adsorption, and the binding geometries in the M-MOF-74 series were reproduced using Monte Carlo simulations in the grand-canonical ensemble. The computed loadings, heats of CO2 and H2 adsorption, and binding geometries in M-MOF-74 are in very good agreement with the experimental values. The temperature transferability of the force field from 77 K to 87 K, and 298 K was shown for adsorption of H2. The validated force field was used to study the adsorption and separation of CO2/H2 mixtures at 298 K. The adsorption of H2 practically does not occur when CO2 is present in the mixture. As indicated from simulated breakthrough curves, the breakthrough time of CO2 in M-MOF-74 follows the same order a, Engineering Thermodynamics

The use of reactive working fluids in thermodynamic cycles is currently being considered as an alternative to inert working fluids, because of the preliminarily attested higher energy-efficiency potential. The current needs to simulate their use in thermodynamic cycles, which may operate in liquid, vapour or vapour-liquid state, are an accurate real-fluid equation of state and ideal gas thermochemical properties of each molecule constituting the mixture, to calculate the equilibrium constant. To this end, the appeal to a multi-scale theoretical methodology is paramount and its definition represents the objective of the present work. This methodology is applied and validated on the system N2O4 ⇌ 2NO2. Firstly, the equations solved for simultaneous two-phase and reaction equilibrium are presented. Secondly, ideal gas thermochemical properties of N2O4 and NO2 are computed at atomic scale by quantum mechanics simulations. Then, to apply the selected cubic equation of state, pure-component properties of the species forming the reactive mixture (critical point coordinates and acentric factor) are required as input. However, these properties are not measurable, since NO2 and N2O4 do not exist in nature as pure components. To get around this difficulty, the methodology relies on molecular Monte Carlo simulations of the pure N2O4 and NO2, as well as on the reactive N2O4 ⇌ 2NO2, enabling the determination of those missing pure-component properties and thus the calculation, on a macroscopic scale, of the reactive mixture properties. Finally, the comparison of calculated mixture properties with available experimental data leads to validate the accuracy of the proposed methodology., Engineering Thermodynamics, Process and Energy

Thermodynamic factors for diffusion connect the Fick and Maxwell-Stefan diffusion coefficients used to quantify mass transfer. Activity coefficient models or equations of state can be fitted to experimental or simulation data, from which thermodynamic factors can be obtained by differentiation. The accuracy of thermodynamic factors determined using indirect routes is dictated by the specific choice of an activity coefficient model or an equation of state. The Permuted Widom’s Test Particle Insertion (PWTPI) method developed by Balaji et al. enables direct determination of thermodynamic factors in binary and multicomponent systems. For highly dense systems, for example, typical liquids, it is well known that molecular test insertion methods fail. In this article, we use the Continuous Fractional Component Monte Carlo (CFCMC) method to directly calculate thermodynamic factors by adopting the PWTPI method. The CFCMC method uses fractional molecules whose interactions with their surrounding molecules are modulated by a coupling parameter. Even in highly dense systems, the CFCMC method efficiently handles molecule insertions and removals, overcoming the limitations of the PWTPI method. We show excellent agreement between the results of the PWTPI and CFCMC methods for the calculation of thermodynamic factors in binary systems of Lennard-Jones molecules and ternary systems of Weeks-Chandler-Andersen molecules. The CFCMC method applied to calculate the thermodynamic factors of realistic molecular systems consisting of binary mixtures of carbon dioxide and hydrogen agrees well with the NIST REFPROP database. Our study highlights the effectiveness of the CFCMC method in determining thermodynamic factors for diffusion, even in densely packed systems, using relatively small numbers of molecules., Reservoir Engineering, Engineering Thermodynamics

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties., Engineering Thermodynamics, Team Poulumi Dey

Hydrogen peroxide plays a key role in many environmental and industrial chemical processes. We performed classical Molecular Dynamics and Continuous Fractional Component Monte Carlo simulations to calculate thermodynamic properties of H2O2 in aqueous solutions. The quality of the available force fields for H2O2 developed by Orabi and English (2018) [67], and by Cordeiro (2014) [69] was systematically evaluated. To assess which water force field is suitable for predicting properties of H2O2 in aqueous solutions, four widely used water force fields were used, namely the TIP3P, TIP4P/2005, TIP5P-E, and a modified TIP3P force field. While the computed densities of pure H2O2 in the temperature range of 253 - 353 K using the force field by Orabi & English are in excellent agreement with experimental results, the densities using the force field by Cordeiro are underestimated by 3%. The TIP4P/2005 force field in combination with the H2O2 force field developed by Orabi & English can predict the densities of H2O2 aqueous solution for the whole range of H2O2 mole fractions in very good agreement with experimental results. The TIP4P/2005 force field in combination with either of the H2O2 force fields can predict the viscosities of H2O2 aqueous solutions for the whole range of H2O2 mole fractions in reasonably good agreement with experimental results. The computed diffusion coefficients for H2O2 and water molecules using the TIP4P/2005 force field with either of the H2O2 force fields are almost constant for the whole range of H2O2 mole fractions. Hydrogen bond analysis showed a steady increase in the number of hydrogen bonds with the solute concentrations in H2O2 aq, Engineering Thermodynamics

Water transport is one of the major factors responsible for moisture damage in asphalt pavements. To study the thermodynamics and kinetics of water transport in bitumen and to uncover microscale mechanisms of moisture-induced damage, molecular dynamics simulations were performed for up to 600 ns for water–bitumen systems with realistic water contents that varied from 0 to 1.76 wt%. Hydrogen bonding interactions and clustering of water molecules at various combinations of temperature and water content were investigated, and their effects on the self-diffusion coefficient of water and bitumen properties are computed and discussed. It is shown that the saturated water concentration in bitumen is small, especially at low temperatures, and additional water molecules tend to form large water clusters via hydrogen bonding, indicating micro-phase separation of the water and bitumen phases inside the simulation box. Hydrogen bonding and water clustering play a crucial role on the magnitude of the self-diffusion coefficient of water. Physical properties of bitumen that include viscosity and cohesive energy are affected by water. The presence of large water clusters is indicative of how degradation in cohesion is observed on the microscale., Pavement Engineering, Process and Energy, Engineering Thermodynamics

Due to the intermittency of renewable energy sources, alkaline water electrolyzers are typically operated at partial load compared to the nominal design value. It is well-known that gas crossover is dominant at low current densities leading to higher anodic hydrogen content and higher cathodic oxygen content in the separator tanks. High anodic hydrogen content is tantamount to loss of product hydrogen which results in an explosive atmosphere in the gas phase if the volumetric hydrogen content in oxygen exceeds 4%. We have developed a transient model of a multi-cell stack which can describe the operation of the electrolyzer with mixed electrolyte flows (anolyte and catholyte), separated flows, or a combination thereof (dynamic switching). This is a major extension of the steady-state model developed by Haug et al. (International Journal of Hydrogen Energy, 2017, 42, 15,689–15707). In sharp contrast to the steady-state model by Haug et al., the transient model can calculate the gas crossover as the operating conditions (e.g. electrolyte flow cycles) dynamically change in time. Depending on the size of the stack and the separator tanks, the model estimates different rates for impurities to build up. The transient model is validated using independent experimental results by Haug et al. and Brauns et al. (Electrochimica Acta, 2022, 404, 139,715) The results show that the dynamic model can follow experimental results for fluctuating current densities for a period of several days. We found that the dynamic response and transition time to steady state depend significantly on the geometrical volume of the gas separators with respect to the single-cell stack. For a multi-cell stack, we find that the impurities build-up faster when increasing the number of cells in the stack. This model serves as a tool for sizing and process management of the electrolyzer system and the separator tanks especially with respect to explosion safety., Engineering Thermodynamics, Process and Energy

The removal of acid gas impurities from synthesis gas or natural gas can be achieved using several physical solvents. Examples of solvents applied on a commercial scale include methanol (Rectisol), poly(ethylene glycol) dimethyl ethers (Selexol), n-methyl-2-pyrrolidone (Purisol), and propylene carbonate (Fluor solvent). Continuous Fractional Component Monte Carlo (CFCMC) simulations in the osmotic ensemble were used to compute the Henry coefficients of the pure gases CO (Formula presented.), CH (Formula presented.), C (Formula presented.) H (Formula presented.), CO, H (Formula presented.), N (Formula presented.), N (Formula presented.) O, and H (Formula presented.) S in the aforementioned solvents. The predicted Henry coefficients are in good agreement with the experimental results. The Monte Carlo method correctly predicts the gas solubility trend in these physical solvents, which obeys the following order: H (Formula presented.) S > CO (Formula presented.) > C (Formula presented.) H (Formula presented.) > CH (Formula presented.) > CO > N (Formula presented.) > H (Formula presented.). The gas separation selectivities for the precombustion process and the natural gas sweetening process are calculated from the pure gas Henry coefficients. The CO (Formula presented.) /N (Formula presented.) O analogy is verified for the solubility in these solvents., Engineering Thermodynamics

Electrochemical reduction of carbon dioxide (CO2) to useful products is an emerging power-to-X concept, which aims to produce chemicals and fuels with renewable electricity instead of fossil fuels. Depending on the catalyst, a range of chemicals can be produced from CO2 electrolysis at industrial-scale current densities, high Faraday efficiencies, and relatively low cell voltages. One of the main challenges for up-scaling the process is related to (bi)carbonate formation (carbonation), which is a consequence of performing the reaction in alkaline media to suppress the competing hydrogen evolution reaction. The parasitic reactions of CO2 with the alkaline electrolytes result in (bi)carbonate precipitation and flooding in gas diffusion electrodes, CO2 crossover to the anode, low carbon utilization efficiencies, electrolyte carbonation, pH-drift in time, and additional cost for CO2 and electrolyte recycling. We present a critical review of the causes, consequences, and possible solutions for the carbonation effect in CO2 electrolyzers. The mechanism of (bi)carbonate crossover in different cell configurations, its effect on the overall process design, and the economics of CO2 and electrolyte recovery are presented. The aim is to provide a better understanding of the (bi)carbonate problem and guide research directions to overcome the challenges related to low-temperature CO2 electrolysis in alkaline media., Process and Energy, Team Peyman Taheri, Engineering Thermodynamics

Formic acid production from CO2 allows the reduction of carbon dioxide emissions while synthesizing a product with a wide range of applications. CO2 hydrogenation is challenging due to the cost of transition metal catalysts and the toxicity of the transition elements. In this work, the thermodynamic confinement effects of the metal–organic framework UiO-66 on the CO2 hydrogenation to formic acid were studied by force field-based molecular simulations. The confinement effects of UiO-66 and the metal–organic frameworks Cu-BTC, and IRMOF-1 were compared, to assess the impact of different pore size and metal centers on the production of HCOOH. Monte Carlo simulations in the grand-canonical ensemble were performed in the frameworks, using gas phase mole fractions of CO2, H2, and HCOOH at chemical equilibrium, obtained from Continuous Fractional Component Monte Carlo simulations in the Reaction Ensemble. The adsorption isobars of the components in metal–organic frameworks were computed at 298.15 – 800 K, 1 – 60 bar. The enhancement of HCOOH production due to preferential adsorption of HCOOH in metal–organic frameworks was calculated for all studied conditions. UiO-66, Cu-BTC, and IRMOF-1 affect CO2 hydrogenation reaction, shifting the thermodynamical equilibrium toward HCOOH formation. The prevailing factor is the type of metal center in the metal–organic framework. The confinement effect of Cu-BTC turns out to exceed the enhancement caused by UiO-66, and IRMOF-1. The resulting mole fraction of HCOOH increased by ca. 2000 times compared to the gas phase at 298.15 K, 60 bar. Cu-BTC can be considered as an alternative to improve the production of HCOOH due to elimination of the high-cost temperature elevation, cost reduction of downstream processing methods, and comparable final concentration of HCOOH to the reported concentrations of formate obtained using transition metal catalysts., Process and Energy

We present the RUPTURA code (https://github.com/iraspa/ruptura) as a free and open-source software package (MIT license) for (1) the simulation of gas adsorption breakthrough curves, (2) mixture prediction using methods like the Ideal Adsorption Solution Theory (IAST), segregated-IAST and explicit isotherm models, and (3) fitting of isotherm models on computed or measured adsorption isotherm data. The combination with the RASPA software enables computation of breakthrough curves directly from adsorption simulations in the grand-canonical ensemble. RUPTURA and RASPA have similar input styles. IAST is implemented near machine precision but we also provide several explicit mixture prediction methods that are non-iterative and potentially faster than IAST. The code supports a wide variety of isotherm models like Langmuir, Anti-Langmuir, BET, Henry, Freundlich, Sips, Langmuir-Freundlich, Redlich-Peterson, Toth, Unilan, O'Brian & Myers, Asymptotic Temkin, and Bingel & Walton. The isotherm model parameters can easily be obtained by the fitting module. Breakthrough plots and animations of the column properties are automatically generated. In addition to highlighting the code, we also review all the developed techniques from literature for mixture prediction, breakthrough simulations, and isotherm model fitting, and provide a tutorial discussing the workflows., Engineering Thermodynamics, Process and Energy

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions., Engineering Thermodynamics, Team Poulumi Dey

Experimentally measuring the diffusivities of CO2 and H2S in aqueous alkanolamine solutions presents an extremely challenging task. To overcome this challenge, we performed Molecular Dynamics (MD) simulations to study the effects of temperature and N-methyldiethanolamine (MDEA) concentration on self-diffusivities of CO2 (DCO2) and H2S (DH2S) in aqueous MDEA solutions. We compute the densities and viscosities of aqueous MDEA solutions for an MDEA concentration range of 10–50 wt% and a temperature range of 288–333 K showing an excellent agreement with experimental data from literature. We compute the self-diffusivity of MDEA (DMDEA) in aqueous MDEA solutions and our findings show that the computed values of DMDEA are in excellent agreement with experimental and simulation results from literature. The self-diffusivities DCO2 and DH2S in aqueous MDEA solutions are computed for a wide range of temperatures and MDEA concentrations and our results show that both DCO2 and DH2S depend significantly on temperature and MDEA concentration. We also show that both CO2 and H2S diffuse slower in aqueous MDEA solutions than in aqueous MEA solutions. By comparing the radial distribution functions of CO2, H2S, water, and MDEA, we show that H2S has stronger interactions with the surrounding molecules than CO2, which makes H2S diffuse slower in aqueous MDEA solutions. We also investigate the densities and viscosities of acid gas loaded aqueous MDEA solutions and self-diffusivities of the reaction products of CO2 and H2S with aqueous MDEA solutions. We show that the self-diffusivities of CO2-loaded solutions significantly decrease with increasing CO2 loading while the sel, Engineering Thermodynamics

It has been reported that lipid droplets (LDs), called oleosomes, have an inherent ability to inflate or shrink when absorbing or fueling lipids in the cells, showing that their phospholipid/protein membrane is dilatable. This property is not that common for membranes stabilizing oil droplets and when well understood, it could be exploited for the design of responsive and metastable droplets. To investigate the nature of the dilatable properties of the oleosomes, we extracted them from rapeseeds to obtain an oil-in-water emulsion. Initially, we added an excess of rapeseed oil in the dispersion and applied high-pressure homogenization, resulting in a stable oil-in-water emulsion, showing the ability of the molecules on the oleosome membrane to rearrange and reach a new equilibrium when more surface was available. To confirm the rearrangement of the phospholipids on the droplet surface, we used molecular dynamics simulations and showed that the fatty acids of the phospholipids are solubilized in the oil core and are homogeneously spread on the liquid-like membrane, avoiding clustering with neighbouring phospholipids. The weak lateral interactions on the oleosome membrane were also confirmed experimentally, using interfacial rheology. Finally, to investigate whether the weak lateral interactions on the oleosome membrane can be used to have a triggered change of conformation by an external force, we placed the oleosomes on a solid hydrophobic surface and found that they destabilise, allowing the oil to leak out, probably due to a reorganisation of the membrane phospholipids after their interaction with the hydrophobic surface. The weak lateral interactions on the LD membrane and their triggered destabilisation present a unique property that can be used for a targeted release in foods, pharmaceuticals and cosmetics., Engineering Thermodynamics

Sodium borohydride (NaBH4) has a high hydrogen (H2 ) gravimetric capacity of 10.7 wt %. NaBH4 releases H2 through a hydrolysis reaction in which aqueous NaB(OH)4 is formed as a byproduct. NaB(OH)4 strongly influences the thermophysical properties of aqueous solutions (i.e., densities, viscosities, and electrical conductivities) and the hydrolysis reaction kinetics and conversion of NaBH4. Here, molecular dynamics (MD) simulations are performed to compute viscosities, electrical conductivities, and self-diffusivities of H2 , Na+, and B(OH)4- for a temperature and concentration range of 298-353 K and 0-5 mol NaB(OH)4/kg water, respectively. Continuous fractional component Monte Carlo (CFCMC) simulations are used to compute the solubilities of H2 and activities of water in aqueous NaB(OH)4 solutions for the same temperature and concentration range. A new force field is developed (Delft force field of B(OH)4-: DFF/B(OH)4-) in which B(OH)4- is modeled as a tetrahedral structure with a scaled charge of −0.85. The OH group in B(OH)4- is modeled as a single interaction site. This force field is based on TIP4P/2005 water and the Madrid-2019 Na+ force field. The MD simulations can accurately capture the densities and viscosities within 2.5% deviation from available experimental data at 298 K up to a concentration of 5 mol NaB(OH)4/kg water. The computed electrical conductivities deviate by ca. 10% from experimental data at 298 K for the same concentration range. Based on the molecular simulations results, engineering equations are developed for shear viscosities, self-diffusivities of H2, Na+, and B(OH)4-, and solubilities of H2, which can be used to design an, Engineering Thermodynamics, Complex Fluid Processing, Team Poulumi Dey

Ideal Adsorbed Solution Theory (IAST) is a common method for modelling mixture adsorption isotherms based on pure component isotherms. When the adsorbent has distinct adsorption sites, the segregated version of IAST (SIAST) provides improved adsorbed loadings compared to IAST. We have adopted the concept of SIAST and applied it to an explicit isotherm model which takes into account the different sizes of the adsorbates: the so called Segregated Explicit Isotherm (SEI). The purpose of SEI is to have an explicit adsorption model that can consider both size-effects of the co-adsorbed molecules and surface heterogeneities. In sharp contrast to IAST and SIAST, no iterative scheme is required in case of SEI, which leads to much faster simulations. A comparative study has been performed to analyse the adsorption isotherms calculated using these three methods. The adsorbed loadings predicted by SEI and SIAST are in excellent agreement with the Grand-Canonical Monte Carlo (GCMC) simulation data. The loadings estimated by IAST show considerable deviations from the GCMC data at high pressures. Breakthrough curve modelling is used to compare the effects of these three models at dynamic conditions. The explicit model (SEI) leads to the fastest simulation run time, followed by SIAST., Engineering Thermodynamics

Hydrate dissociation is often accompanied by the formation of nanobubbles. Knowledge of the effects of nanobubbles on hydrate dissociation is essential for understanding the dynamic behavior of the hydrate phase change and improving the gas production efficiency. Here, molecular dynamics simulations were performed to study the methane hydrate dissociation kinetics with and without a pre-existing methane nanobubble. The results show that the hydrate cluster in the liquid phase dissociates layer-by-layer. This process is shown to be independent of the temperature and nanobubble presence at the simulation conditions. Hydrate dissociation does not always lead to nanobubble formation because the supersaturated methane solution can be stable for a long time. A steep methane concentration gradient was observed between the hydrate cluster surface and the methane nanobubble, which can enhance the directional migration of methane and effectively minimize the methane concentration in the liquid phase, thereby increasing the driving force for the hydrate dissociation. Our findings indicate that the presence of a nanobubble near the hydrate surface does not decrease the activation energy of hydrate dissociation, but it can increase the intrinsic decomposition rate. The average hydrate dissociation rate is linearly correlated with the mass flow rate towards the nanobubble. The mass flow rate is determined by the nanobubble size and hydrate-nanobubble distance. Our findings contribute to the fundamental understanding of the dissociation mechanism of gas hydrates in the liquid phase, which is crucial for the design and optimization of efficient gas hydrate production techniques., Engineering Thermodynamics

Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%., Geoscience and Engineering, Engineering Thermodynamics, Team Poulumi Dey

Vapor-liquid equilibrium (VLE) data for the binary systems tetrahydrofuran (THF) + acetic acid (AA) and THF + trichloroethylene (TCE) were measured under isobaric conditions using an ebulliometer. The boiling temperatures for the systems (THF + AA/THF + TCE) are reported for 13/15 compositions and five/six different pressures ranging from 50.2/60.0 to 101.1/101.3 kPa, respectively. The THF + AA system shows simple phase behavior with no azeotrope formation. The THF + TCE system does not exhibit azeotrope formation but seems to have a pinch point close to the pure end of TCE. The nonrandom two-liquid (NRTL) and universal quasichemical (UNIQUAC) activity coefficient models were used to accurately fit the binary (PTx) data. Both models were able to fit the binary VLE data satisfactorily. However, the NRTL model was found to be slightly better than UNIQUAC model in fitting the VLE data for both systems. The results can be used for designing liquid-liquid extraction and distillation processes involving mixtures of THF, AA, and TCE., Engineering Thermodynamics

Recently, deep eutectic solvents (DES) have been considered as possible electrolytes for the electrochemical reduction of CO2 to value-added products such as formic and oxalic acids. The applicability of pure DES as electrolytes is hindered by high viscosities. Mixtures of DES with organic solvents can be a promising way of designing superior electrolytes by exploiting the advantages of each solvent type. In this study, densities, viscosities, diffusivities, and ionic conductivities of mixed solvents comprising DES (i.e., reline and ethaline), methanol, and propylene carbonate were computed using molecular simulations. To provide a quantitative assessment of the affinity and mass transport of CO2 and oxalic and formic acids in the mixed solvents, the solubilities and self-diffusivities of these solutes were also computed. Our results show that the addition of DES to the organic solvents enhances the solubilities of oxalic and formic acids, while the solubility of CO2 in the ethaline-containing mixtures are in the same order of magnitude with the respective pure organic components. A monotonic increase in the densities and viscosities of the mixed solvents is observed as the mole fraction of DES in the mixture increases, with the exception of the density of ethaline-propylene carbonate which shows the opposite behavior due to the high viscosity of the pure organic component. The self-diffusivities of all species in the mixtures significantly decrease as the mole fraction of DES approaches unity. Similarly, the self-diffusivities of the dissolved CO2 and the oxalic and formic acids also decrease by at least 1 order of magnitude as the composition of the mixture shifts from the pure organic component to pure DES. The computed ionic conductivities of all mixed solvents show a maximum value for mole fractions of DES in the range from 0.2 to 0.6 and decrease as more DES is added to the mixtures. Since for most mixtures studied here no prior experimental measurements ex, Engineering Thermodynamics

The electrochemical CO2 reduction reaction (CO2RR) is important for a sustainable future. Key insights into the reaction pathways have been obtained by density functional theory (DFT) analysis, but so far, DFT has been unable to give an overall understanding of selectivity trends without important caveats. We show that an unconsidered parameter in DFT models of electrocatalysts-the surface coverage of reacting species-is crucial for understanding the CO2RR selectivities for different surfaces. Surface coverage is a parameter that must be assumed in most DFT studies of CO2RR electrocatalysts, but so far, only the coverage of nonreacting adsorbates has been treated. Explicitly treating the surface coverage of reacting adsorbates allows for an investigation that can more closely mimic operating conditions. Furthermore, and of more immediate importance, the use of surface coverage-dependent adsorption energies allows for the extraction of ratios of adsorption energies of CO2RR intermediates (COOHads and HCOOads) that are shown to be predictive of selectivity and are not susceptible to systematic errors. This approach allows for categorization of the selectivity of several monometallic catalysts (Pt, Pd, Au, Ag, Zn, Cu, Rh, W, Pb, Sn, In, Cd, and Tl), even problematic ones such as Ag or Zn, and does so by only considering the adsorption energies of known intermediates. The selectivity of the further reduction of COOHads can now be explained by a preference for Tafel or Heyrovsky reactions, recontextualizing the nature of selectivity of some catalysts. In summary, this work resolves differences between DFT and experimental studies of the CO2RR and underlines the importance of surface coverage., Large Scale Energy Storage, Engineering Thermodynamics, Process and Energy

The Kirkwood-Buff theory is a cornerstone of the statistical mechanics of liquids and solutions. It relates volume integrals over the radial distribution function, so-called Kirkwood-Buff integrals (KBIs), to particle number fluctuations and thereby to various macroscopic thermodynamic quantities such as the isothermal compressibility and partial molar volumes. Recently, the field has seen a strong revival with breakthroughs in the numerical computation of KBIs and applications to complex systems such as bio-molecules. One of the main emergent results is the possibility to use the finite volume KBIs as a tool to access finite volume thermodynamic quantities. The purpose of this Perspective is to shed new light on the latest developments and discuss future avenues., Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Engineering Thermodynamics

Methane hydrate dissociation kinetics can be inhibited in NaCl solutions; however, this effect is reversed by promoting bubble formation that enhances dissociation. The negative and positive effects of inorganic salt injection on gas production from hydrate-bearing sediments are still controversial. Here, molecular dynamics simulations were performed to investigate the characteristics of NaCl solution invasion into hydrate-occupied nanopores and the effects on the confined hydrate dissociation kinetics. Two initial configurations comprising liquid and silica pore phases were studied with a low or high NaCl concentration, respectively. The results show that, under the simulation conditions, salt invasion decelerated hydrate dissociation within the silica pore as NaCl invasion into the pore is stepwise. Initially, few ions can diffuse into the pore phase, and gas nanobubbles form on the solid surface mainly via confinement and surface effects, independent of NaCl solution invasion. Subsequently, gradual salt diffusion immersed the residual hydrate in the salt solution and hindered hydrate decomposition until the dissociation finished. More ions could diffuse into the pore phase at the high NaCl concentrations with a low diffusion efficiency, leading to surface nanobubble growth toward the residual hydrate and somewhat accelerated hydrate dissociation. This severely hinders the escape of released methane from the pore. This study yields molecular-level insight into the origin of the negative effect of salt invasion on hydrate dissociation, which should be avoided during gas production from hydrate reservoirs with low permeabilities via salt injection combined with thermal stimulation., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Engineering Thermodynamics

In this work, adsorption properties of the UiO-66 metal–organic framework were investigated, with particular emphasis on the influence of structural defects. A series of UiO-66 samples were synthesized and characterized using a wide range of experimental techniques. Type I adsorption isotherms for low-temperature adsorption of N2 and Ar showed that micropore volume and specific surface area significantly increase with the number of defects. Adsorption of hexane isomers in UiO-66 was studied by means of quasi-equilibrated temperature-programmed desorption and adsorption (QE-TPDA) experimental and Monte Carlo simulation techniques. QE-TPDA profiles revealed that only defect-free UiO-66 exhibits distinct two adsorption states. This technique also yielded high-quality adsorption isobars that were successfully recreated using Grand-Canonical Monte Carlo molecular simulations, which, however, required refinement of the existing force fields. The calculations demonstrated the detailed mechanism of adsorption and separation of hexane isomers in the UiO-66 structure. The preferred tetrahedral cages provide suitable voids for bulky molecules, which is the reason for unusual “reverse” selectivity of UiO-66 towards di-branched alkanes. Interconnection of the tetrahedral cavities due to missing organic linkers greatly reduces the selectivity of the defected material., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Engineering Thermodynamics

There is a growing interest in the development of routes to produce formic acid from CO2, such as the electrochemical reduction of CO2 to formic acid. The solubility of CO2 in the electrolyte influences the production rate of formic acid. Here, the dependence of the CO2 solubility in aqueous HCOOH solutions with electrolytes on the composition and the NaCl concentration was studied by Continuous Fractional Component Monte Carlo simulations at 298.15 K and 1 bar. The chemical potentials of CO2, H2O, and HCOOH were obtained directly from single simulations, enabling the calculation of Henry coefficients and subsequently considering salting in or salting out effects. As the force fields for HCOOH and H2O may not be compatible due to the presence of strong hydrogen bonds, the Gibbs-Duhem integration test was used to test this compatibility. The combination of the OPLS/AA force field with a new set of parameters, in combination with the SPC/E force field for water, was selected. It was found that the solubility of CO2 decreases with increasing NaCl concentration in the solution and increases with the increase of HCOOH concentration. This continues up to a certain concentration of HCOOH in the solution, after which the CO2 solubility is high and the NaCl concentration has no significant effect., Engineering Thermodynamics

The thermophysical properties of aqueous electrolyte solutions are of interest for applications such as water electrolyzers and fuel cells. Molecular dynamics (MD) and continuous fractional component Monte Carlo (CFCMC) simulations are used to calculate densities, transport properties (i.e., self-diffusivities and dynamic viscosities), and solubilities of H2 and O2 in aqueous sodium and potassium hydroxide (NaOH and KOH) solutions for a wide electrolyte concentration range (0-8 mol/kg). Simulations are carried out for a temperature and pressure range of 298-353 K and 1-100 bar, respectively. The TIP4P/2005 water model is used in combination with a newly parametrized OH- force field for NaOH and KOH. The computed dynamic viscosities at 298 K for NaOH and KOH solutions are within 5% from the reported experimental data up to an electrolyte concentration of 6 mol/kg. For most of the thermodynamic conditions (especially at high concentrations, pressures, and temperatures) experimental data are largely lacking. We present an extensive collection of new data and engineering equations for H2 and O2 self-diffusivities and solubilities in NaOH and KOH solutions, which can be used for process design and optimization of efficient alkaline electrolyzers and fuel cells., Engineering Thermodynamics, Process and Energy, Team Poulumi Dey

We performed H-cell and flow cell experiments to study the electrochemical reduction of CO2 to oxalic acid (OA) on a lead (Pb) cathode in various nonaqueous solvents. The effects of anolyte, catholyte, supporting electrolyte, temperature, water content, and cathode potential on the Faraday efficiency (FE), current density (CD), and product concentration were investigated. We show that a high FE for OA can be achieved (up to 90%) at a cathode potential of -2.5 V vs Ag/AgCl but at relatively low CDs (10-20 mA/cm2). The FE of OA decreases significantly with increasing water content of the catholyte, which causes byproduct formation (e.g., formate, glycolic acid, and glyoxylic acid). A process design and techno-economic evaluation of the electrochemical conversion of CO2 to OA is presented. The results show that the electrochemical route for OA production can compete with the fossil-fuel based route for the base case scenario (CD of 100 mA/cm2, OA FE of 80%, cell voltage of 4 V, electrolyzer CAPEX of $20000/m2, electricity price of $30/MWh, and OA price of $1000/ton). A sensitivity analysis shows that the market price of OA has a huge influence on the economics. A market price of at least $700/ton is required to have a positive net present value and a payback time of less than 10 years. The performance and economics of the process can be further improved by increasing the CD and FE of OA by using gas diffusion electrodes and eliminating water from the cathode, lowering the cell voltage by increasing the conductivity of the electrolyte solutions, and developing better OA separation methods., Engineering Thermodynamics, Large Scale Energy Storage

In this manuscript, we show how to compute partial molar properties (e.g. partial molar volumes, energies, and enthalpies) of fluid mixtures from single Molecular Dynamics simulations in the microcanonical or canonical ensemble, using only relatively simple post-processing of trajectories. The method uses least squares linear regression of local fluctuations of particle numbers and energies in combination with the Small System Method, and is in principle valid for any number of components and for any type of intermolecular interactions. For multicomponent systems, only a single simulation is needed for a given composition of the mixture. Simulations of a binary WCA mixture are used to illustrate the method, and to investigate the effect of system size., Engineering Thermodynamics

A reactor model that deconvolutes thermodynamics of adsorption of hydrocarbon in the pores of zeolite Beta, obtained by Configurational-bias Monte Carlo simulations, from intrinsic, intraporous kinetics of hydroisomerization and hydrocracking reactions, provides a good quantitative description of all significant reactions in the kinetic network for interconversion and cracking of different heptane isomers. Activation enthalpies obtained for intraporous reactions follow the expected order according to the carbenium ion formalism: methyl shift< ethyl shift < isom(B) ∼ crack(B2) < crack(B1) < crack(C) ∼ crack(D) < crack(E) and apparently within each isomerization class, in terms of carbenium ions formally involved: sec → tert < sec → sec ∼ tert → tert < tert → sec. except for the ethyl shift reaction forming 3-ethylpentane. Cracking happens primarily through 2,4-dimethylpentane (type B2), regardless of the initial reactant. The model can be subsequently used to separate the effect of pore structure on selective adsorption and on intraporous reaction kinetics. Zeolite Beta will serve as a base case for a comparison of different zeolite structures., Process and Energy, Engineering Thermodynamics

Mass transport of different species plays a crucial role in electrochemical conversion of CO (Formula presented.) due to the solubility limit of CO (Formula presented.) in aqueous electrolytes. In this study, we investigate the transport of CO (Formula presented.) and other ionic species through the electrolyte and the membrane, and its impact on the scale-up process of HCOO (Formula presented.) /HCOOH formation. The mass transport of ions to the electrode and the membrane is modelled at constant current density. The mass transport limitations of CO (Formula presented.) on the formation of HCOO (Formula presented.) /HCOOH is investigated at different pressures ranges from 5–40 bar. The maximum achievable partial current density of formate/formic acid is increased with increasing CO (Formula presented.) pressure. We use an ion exchange membrane model to understand the ion transport behaviour for both the monopolar and bipolar membranes. The cation exchange (CEM) and anion exchange membrane (AEM) model show that ion transport is limited by the electrolyte salt concentrations. For 0.1 M KHCO (Formula presented.), the AEM reaches the limiting current density more quickly than the CEM. For the BPM model, ion transport across the diffusion layer on either side of the BPM is also included to understand the concentration polarization across the BPM. The model revealed that the polarization losses across the bipolar membrane depend on the pH of the electrolyte used for the CO (Formula presented.) reduction reaction (CO 2RR). The polarization loss on the anolyte side decreases with an increasing pH, while, on the cathode side, it increases with increasing catholyte pH. With this combined model for the electrode reactions and the membrane transport, we are able to account for the various factors influencing the polarization losses in the CO (Formula presented.) electrolyzer. To complete the analysis, we simulated the full cell polarization curve and fitted with t, Engineering Thermodynamics

Lithium (Li) dopants have garnered attention as a means to enhance hydrogen (H2) binding energies and capacities in 2D boron-based materials. However, it is unclear if and how these dopants affect H2 dissociation and chemisorption. Using density functional theory (DFT) and nudged elastic band (NEB) calculations, reaction pathways for H2 dissociation on borophene-hydride and striped-borophene are investigated in the presence of Li dopants. Our results indicate that tuning the Li-loading can increase the reversibility and the rate of dehydrogenation for both borophene-hydride and striped borophene. In particular, for Li-doped borophene-hydride, the heat of reaction for H2 release is reduced by more than 85% compared to the pristine structure (1.97 eV/H2). Our results signify that Li-doping can considerably change the H2 chemisorption properties of striped-borophene and borophene-hydride, and can lead to promising materials for H2 storage., Engineering Thermodynamics, Team Poulumi Dey

Cyclodextrins (CDs) are widely used in drug delivery, catalysis, food and separation processes. In this work, a comprehensive simulation study on the diffusion of the native α-, β- and γ-CDs in aqueous solutions is carried out using Molecular Dynamics simulations. The effect of the system size on the computed self-diffusivity is investigated and it is found that the required correction can be as much as 75% of the final value. The effect of the water force field is examined and it is shown that the q4md-CD/TIP4P/2005 force field combination predicts the experimentally measured self-diffusion coefficients of CDs very accurately. The self-diffusion coefficients of the three native CDs were also computed in aqueous-NaCl solutions using the Joung and Cheatham (JC) and the Madrid-2019 force fields. It is found that Na+ ions have higher affinity towards the CDs when the JC force field is used and for this reason the predicted diffusivity of CDs is lower compared to simulations using the Madrid-2019 force field. As a model system for drug delivery and waste-water treatment applications, the diffusion of the β-CD:Ibuprofen inclusion complex in water is studied. In agreement with experiments for similar components, it is shown that the inclusion complex and the free β-CD have almost equal self-diffusion coefficients. Our analysis revealed that this is most likely caused by the almost full inclusion of the ibuprofen in the cavity of the β-CD. Our findings show that Molecular Dynamics simulation can be used to provide reasonable diffusivity predictions, and to obtain molecular-level understanding useful for industrial applications of CDs., Engineering Thermodynamics

We present several new major features added to the Monte Carlo (MC) simulation code Brick-CFCMC for phase- and reaction equilibria calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one is thermodynamic integration for the computation of excess chemical potentials (μex). For this purpose, we implemented the computation of the ensemble average of the derivative of the potential energy with respect to the scaling factor for intermolecular interactions (⟨∂U∂λ⟩). Efficient bookkeeping is implemented so that the quantity ∂U∂λ is updated after every MC trial move with negligible computational cost. We demonstrate the accuracy and reliability of the calculation of μex for sodium chloride in water. Second, we implemented hybrid MC/MD translation and rotation trial moves to increase the efficiency of sampling of the configuration space. In these trial moves, short Molecular Dynamics (MD) trajectories are performed to collectively displace or rotate all molecules in the system. These trajectories are accepted or rejected based on the total energy drift. The efficiency of these trial moves can be tuned by changing the time step and the trajectory length. The new trial moves are demonstrated using MC simulations of a viscous fluid (deep eutectic solvent)., Engineering Thermodynamics

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau., Engineering Thermodynamics

Photoelectrochemical (PEC) splitting of water to make hydrogen is a promising clean-energy technology. The oxygen evolution reaction (OER) largely determines the energy efficiency in PEC water-splitting. Hematite, which is a cheap and sustainable semiconductor material with excellent chemical properties, a favourable band gap (2.1 eV) and composed of earth abundant elements is a suitable model photoanode material for studying OER. To understand the design of energy efficient anodes, it is highly desirable to have mechanistic insight into OER at an atomistic level which can be directly connected to experimentally measured quantities. We present a multiscale computational model of OER which connects the thermodynamics and kinetics of elementary charge transfer reactions in OER to kinetics of OER at laboratory length and time scales. We couple density functional theory (DFT) and DFT based molecular dynamics (DFT-MD) simulations with solvent effects at an atomistic level with kinetic Monte Carlo (kMC) simulations at a coarse-grained level in our multiscale model. The time and applied bias potential dependent surface coverage, which are experimentally not known, and the O2 evolution rate during OER at the hematite-water interface are calculated by the multiscale model. Furthermore, the multiscale model demonstrates the effect of explicitly modelling the interaction of water with the electrode surface via direct adsorption. This journal is, Accepted Author Manuscript, Engineering Thermodynamics

Direct electrochemical reduction of CO2 to C2 products such as ethylene is more efficient in alkaline media, but it suffers from parasitic loss of reactants due to (bi)carbonate formation. A two-step process where the CO2 is first electrochemically reduced to CO and subsequently converted to desired C2 products has the potential to overcome the limitations posed by direct CO2 electroreduction. In this study, we investigated the technical and economic feasibility of the direct and indirect CO2 conversion routes to C2 products. For the indirect route, CO2 to CO conversion in a high temperature solid oxide electrolysis cell (SOEC) or a low temperature electrolyzer has been considered. The product distribution, conversion, selectivities, current densities, and cell potentials are different for both CO2 conversion routes, which affects the downstream processing and the economics. A detailed process design and techno-economic analysis of both CO2 conversion pathways are presented, which includes CO2 capture, CO2 (and CO) conversion, CO2 (and CO) recycling, and product separation. Our economic analysis shows that both conversion routes are not profitable under the base case scenario, but the economics can be improved significantly by reducing the cell voltage, the capital cost of the electrolyzers, and the electricity price. For both routes, a cell voltage of 2.5 V, a capital cost of $10,000/m2, and an electricity price of <$20/MWh will yield a positive net present value and payback times of less than 15 years. Overall, the high temperature (SOEC-based) two-step conversion process has a greater potential for scale-up than the direct electrochemical conversion route. Strategies for integrating the electrochemical CO2/CO conversion process into the existing gas and oil infrastructure are outlined. Current barriers for industrialization of CO2 electrolyzers and possible solutions are discussed as well., Engineering Thermodynamics, Large Scale Energy Storage

Thermochemical heat-storage applications, based on the reversible endo-/exothermic hydration reaction of salts, are intensively investigated to search for compact heat-storage devices. To achieve a truly valuable storage system, progressively complex salts are investigated. For these salts, the equilibrium temperature and pressure conditions are not always easy to predict. However, these conditions are crucial for the design of thermochemical heat-storage systems. A biased grand-canonical Monte Carlo (GCMC) tool is developed, enabling the study of equilibrium conditions at the molecular level. The GCMC algorithm is combined with reactive force field molecular dynamics (ReaxFF), which allows bond formation within the simulation. The Weeks-Chandler-Andersen (WCA) potential is used to scan multiple trial positions for the GCMC algorithm at a small cost. The most promising trial positions can be selected for recomputation with the more expensive ReaxFF. The developed WCA-ReaxFF-GCMC tool was used to study the hydration of MgCl2·nH2O. The simulation results show a good agreement with experimental and thermodynamic equilibriums for multiple hydration levels. The hydration shows that water, present at the surface of crystalline salt, deforms the surface layers and promotes further hydration of these deformed layers. Additionally, the WCA-ReaxFF-GCMC algorithm can be used to study other, non-TCM-related, reactive sorption processes., Engineering Thermodynamics, Process and Energy

Hydrophobic deep eutectic solvents (DESs) have recently gained much attention as water-immiscible solvents for a wide range of applications. However, very few studies exist in which the hydrophobicity of these DESs is quantified. In this work, the interfacial properties of hydrophobic DESs with water were computed at various temperatures using molecular dynamics simulations. The considered DESs were tetrabutylammonium chloride-decanoic acid (TBAC-dec) with a molar ratio of 1:2, thymol-decanoic acid (Thy-dec) with a molar ratio of 1:2, and dl-menthol-decanoic acid (Men-dec) with a molar ratio of 2:1. The following properties were investigated in detail: interfacial tensions, water-in-DES solubilities (and salt-in-water solubilities for TBAC-dec/water), density profiles, and the number densities of hydrogen bonds. Different ionic charge scaling factors were used for TBAC-dec. Thy-dec and Men-dec showed a high level of hydrophobicity with negligible computed water-in-DES solubilities. For charge scaling factors of 0.7 and 1 for the thymol and decanoic acid components of Thy-dec, the computed interfacial tensions of the DESs are in the following order: TBAC-dec (ca. 4 mN m-1) < Thy-dec (20 mN m-1) < Men-dec (26 mN m-1). The two sets of charge scaling factors for Thy-dec did not lead to different density profiles but resulted in considerable differences in the DES/water interfacial tensions due to different numbers of decanoic acid-water hydrogen bonds at the interfaces. Large peaks were observed for the density profiles of (the hydroxyl oxygen of) decanoic acid at the interfaces of all DES/water mixtures, indicating a preferential alignment of the oxygen atoms of decanoic acid toward the aqueous phase., Engineering Thermodynamics

The circular economy demands waste utilization for the production of high-value products, and this requires the development of novel processing routes. In this study, rare earth (La and Ce) oxides were completely (>99%) recovered from polishing waste by a combined novel reductive acid leaching and alkali treatment process. About 70% of rare earths were dissolved during the first leaching step. The undissolved rare earth compounds are converted to oxides/hydroxides by alkali treatment and dissolved in the acid solution – the 2nd leaching step – for the complete recovery of rare earths. The recovered rare earth oxides were used for producing in-situ high-value Al-La-Ce alloys with fused salt electrolysis. Mechanical properties of our Al-La-Ce alloys are similar to the known high temperature Al-Ce alloys. This development of new alloys by our novel process helps in utilization of both overproduced primary La and Ce oxides as well as La and Ce recovered from polishing waste., Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., (OLD) MSE-1, Engineering Thermodynamics, Team Yongxiang Yang, Team Erik Offerman

With the emergence of hydrophobic deep eutectic solvents (DESs), the scope of applications of DESs has been expanded to include situations in which miscibility with water is undesirable. Whereas most studies have focused on the applications of hydrophobic DESs from a practical standpoint, few theoretical works exist that investigate the structural and thermodynamic properties at the nanoscale. In this study, Molecular Dynamics (MD) simulations have been performed to model DESs composed of tetraalkylammonium chloride hydrogen bond acceptor and decanoic acid hydrogen bond donor (HBD) at a molar ratio of 1:2, with three different cation chain lengths (4, 7, and 8). After fine-tuning force field parameters, densities, viscosities, self-diffusivities, and ionic conductivities of the DESs were computed over a wide temperature range. The liquid structure was examined using radial distribution functions (RDFs) and hydrogen bond analysis. The MD simulations reproduced the experimental density and viscosity data from the literature reasonably well and were used to predict diffusivities and ionic conductivities, for which experimental data are scarce or unavailable. It was found that although an increase in the cation chain length considerably affected the density and transport properties of the DESs (i.e., yielding smaller densities and slower dynamics), no significant influence was observed on the RDFs and the hydrogen bonds. The self-diffusivities showed the following order for the mobility of the various components: HBD > anion > cation. Strong hydrogen bonds between the hydroxyl and carbonyl groups of decanoic acid and between the hydroxyl group of decanoic acid and chloride were observed to dominate the intermolecular interactions., Engineering Thermodynamics

Two-dimensional (2D) boron-based materials are receiving much attention as H2 storage media due to the low atomic mass of boron and the stability of decorating alkali metals on the surface, which enhance interactions with H2. This work investigates the suitability of Li, Na, and K decorations on 2D honeycomb borophene oxide (B2O) for H2 storage, using dispersion corrected density functional theory (DFT-D2). A high theoretical gravimetric density of 8.3 wt % H2 is achieved for the Li-decorated B2O structure. At saturation, each Li binds to two H2 with an average binding energy of -0.24 eV/H2. Born-Oppenheimer molecular dynamics simulations at temperatures of 100, 300, and 500 K demonstrate the stability of the Li-decorated structure and the H2 desorption behavior at different temperatures. Our findings indicate that Li-decorated 2D B2O is a promising material for reversible H2 storage and recommend experimental investigation of 2D B2O as a potential H2 storage medium., Process and Energy, Engineering Thermodynamics, Team Poulumi Dey

Formic acid (FA) is an interesting hydrogen (H2) and carbon monoxide (CO) carrier that can be produced by the electrochemical reduction of carbon dioxide (CO2) using renewable energy. The separation of FA from water is challenging due to the strong (cross)association of the components and the presence of a high boiling azeotrope. For the separation of dilute FA solutions, liquid-liquid extraction is preferred over conventional distillation because distilling large amounts of water is very energy-intensive. In this study, we use 2-methyltetrahydrofuran (2-MTHF) to extract FA from the CO2 electrolysis process, which typically contains <20 wt % of FA. Vapor-liquid equilibrium (VLE) data of the binary system 2-MTHF-FA and liquid-liquid equilibrium (LLE) data of the ternary system 2-MTHF-FA-water are obtained. Continuous extraction and distillation experiments are performed to test the extraction power and recovery of 2-MTHF from the extract. The VLE and LLE data are used to design a hybrid extraction and distillation process to produce a commercial grade product (85 wt % of FA). A detailed economic analysis of this hybrid extraction-distillation process is presented and compared with the existing FA separation methods. It is shown that 2-MTHF is a cost-effective solvent for FA extraction from dilute streams (<20 wt % FA)., Engineering Thermodynamics, Large Scale Energy Storage