Reactive Grand-Canonical Monte Carlo Simulations for Modeling Hydration of MgCl2.

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.

Gibbs Ensemble Monte Carlo for Reactive Force Fields to Determine the Vapor-Liquid Equilibrium of CO2 and H2O.

Absorption and reactive properties of fluids in porous media are key to the design and improvement of numerous energy related applications. Molecular simulations of these systems require accurate force fields that capture the involved chemical reactions and have the ability to describe the vapor-liquid equilibrium (VLE). Two new reactive force fields (ReaxFF) for CO2 and H2O are developed, which are capable of not only modeling bond breaking and formation in reactive environments but also predicting their VLEs at saturation conditions. These new force fields include extra terms (ReaxFF-lg) to improve the long-range interactions between the molecules. For validation, we have developed a new Gibbs ensemble Monte Carlo (GEMC-ReaxFF) approach to predict the VLE. Computed VLE data show good agreement with National Institute of Standards and Technology reference data as well as existing nonreactive force fields. This validation proves the applicability of the GEMC-ReaxFF method to test new reactive force fields, and simultaneously it proves the applicability to extend newly developed ReaxFF force fields to other more complex reactive systems.

Industrial lignin from 2G biorefineries - Assessment of availability and pricing strategies.

With a view to boost practical implementation of lignin conversion technologies, this paper assesses the availability of industrial lignin and evaluates pricing strategies applicable to multi-product biorefineries. The biorefineries, producing either denatured ethanol or sugar hydrolysate as a main product, can yield 43% and 61% of lignin residue (LR) comprising 33% and 23% of lignin by mass, respectively, without sacrificing the output of the main product and before electricity import has become indispensable. Analysis of the pricing strategies reveals that LR must be treated as a low-value by-product, and its minimum selling price (MSP) is driven mainly by the prevailing electricity price. Under the biorefinery net zero energy balance, and taking into account the LR market price adequacy, as well as the main probabilistic conditions, the upper range for the MSP is calculated at $43-70 and $18-37 per ton for biorefineries producing ethanol and hydrolysate, respectively.

Water nucleation in helium, methane, and argon: A molecular dynamics study.

Nucleation of highly supersaturated water vapor in helium, methane, and argon carrier gases at 350 K was investigated using molecular dynamics simulations. Nucleation rates obtained from the mean first passage time (MFPT) method are typically one order of magnitude lower than those from the Yasuoka and Matsumoto method, which can be attributed to the overestimation of the critical cluster size in the MFPT method. It was found that faster nucleation will occur in carrier gases that have better thermalization properties such that latent heat is removed more efficiently. These thermalization properties are shown to be strongly dependent on the molecular mass and Lennard-Jones (LJ) parameters. By varying the molecular mass, for unaltered LJ parameters, it was found that a heavier carrier gas removes less heat although it has a higher collision rate with water than a lighter carrier. Thus, it was shown that a clear distinction between water vapor-carrier gas collisions and water cluster-carrier gas collisions is indispensable for understanding the effect of collision rates on thermalization. It was also found that higher concentration of carrier gas leads to higher nucleation rate. The nucleation rates increased by a factor of 1.3 for a doubled concentration and by almost a factor of two for a tripled concentration.

Homogeneous nucleation of water in argon. Nucleation rate computation from molecular simulations of TIP4P and TIP4P/2005 water model.

Molecular dynamics (MD) simulations were conducted to study nucleation of water at 350 K in argon using TIP4P and TIP4P/2005 water models. We found that the stability of any cluster, even if large, strongly depends on the energetic interactions with its vicinity, while the stable clusters change their composition almost entirely during nucleation. Using the threshold method, direct nucleation rates are obtained. Our nucleation rates are found to be 1.08×1027 cm-3 s-1 for TIP4P and 2.30×1027 cm-3 s-1 for TIP4P/2005. The latter model prescribes a faster dynamics than the former, with a nucleation rate two times larger due to its higher electrostatic charges. The non-equilibrium water densities derived from simulations and state-of-art equilibrium parameters from Vega and de Miguel [J. Chem. Phys. 126, 154707 (2007)] are used for the classical nucleation theory (CNT) prediction. The CNT overestimates our results for both water models, where TIP4P/2005 shows largest discrepancy. Our results complement earlier data at high nucleation rates and supersaturations in the Hale plot [Phys. Rev. A 33, 4156 (1986)], and are consistent with MD data on the SPC/E and the TIP4P/2005 model.

Cross-plane heat transfer through single-layer carbon structures.

Graphene-based nano-structures have been recently proposed to function as additives to improve the conductivity of thermally sluggish phase change materials (PCMs). Based on the existing research studies, the improvement is dependent not only on the matrix material, but also on the geometry of the carbon structure. To gain more insight into the nano-scale thermal transport problem, we launched the current pilot research using water as the matrix material, to represent the hydroxyl-group-rich sugar alcohols as PCMs. We have found that the heat conduction across a graphene layer to water is much faster than the heat conduction to the graphene layer itself. Also, the high graphene-water thermal contact resistance fails to acknowledge the fast thermal kinetics of the low frequency phonons. In the investigation of the geometry effect, the cross-plane heat transfer coefficient is found to decrease with decreasing CNT diameter except CNT(9,9).