First-Principles Study of Chemical Mixtures of CaCl2 and MgCl2 Hydrates for Optimized Seasonal Heat Storage.

Chloride-based salt hydrates form a promising class of thermochemical materials (TCMs), having high storage capacity and fast kinetics. In the charging cycles of these hydrates however hydrolysis might appear along with dehydration. The HCl produced during the hydrolysis degrades and corrodes the storage system. Our GGA-DFT results show that the enthalpy charge during proton formation (an important step in hydrolysis) is much higher for CaCl2·2H2O (33.75 kcal/mol) than for MgCl2·2H2O (19.55 kcal/mol). This is a strong indicator that hydrolysis can be minimized by appropriate chemical mixing of CaCl2 and Mg Cl2 hydrates, which is also confirmed by recent experimental studies. GGA-DFT calculations were performed to obtain and analyze the optimized structures, charge distributions, bonding indicators and harmonic frequencies of various chemical mixtures hydrates and compared them to their elementary salts hydrates. We have further assessed the equilibrium products concentration of dehydration/hydrolysis of the chemical mixtures under a wide range of operating conditions. We observed that chemical mixing leads to an increase of the onset hydrolysis temperature with a maximum value of 79 K, thus increasing the resistance against hydrolysis with respect to the elementary salt hydrates. We also found that the chemical mixing of CaCl2 and MgCl2 hydrates widens the operating dehydration temperature range by a maximum value of 182 K (CaMg2Cl6·2H2O) and lowers the binding enthalpy with respect to the physical mixture by ≈65 kcal/mol for TCM based heat storage systems.

Gas-surface interactions using accommodation coefficients for a dilute and a dense gas in a micro- or nanochannel: heat flux predictions using combined molecular dynamics and Monte Carlo techniques.

The influence of gas-surface interactions of a dilute gas confined between two parallel walls on the heat flux predictions is investigated using a combined Monte Carlo (MC) and molecular dynamics (MD) approach. The accommodation coefficients are computed from the temperature of incident and reflected molecules in molecular dynamics and used as effective coefficients in Maxwell-like boundary conditions in Monte Carlo simulations. Hydrophobic and hydrophilic wall interactions are studied, and the effect of the gas-surface interaction potential on the heat flux and other characteristic parameters like density and temperature is shown. The heat flux dependence on the accommodation coefficient is shown for different fluid-wall mass ratios. We find that the accommodation coefficient is increasing considerably when the mass ratio is decreased. An effective map of the heat flux depending on the accommodation coefficient is given and we show that MC heat flux predictions using Maxwell boundary conditions based on the accommodation coefficient give good results when compared to pure molecular dynamics heat predictions. The accommodation coefficients computed for a dilute gas for different gas-wall interaction parameters and mass ratios are transferred to compute the heat flux predictions for a dense gas. Comparison of the heat fluxes derived using explicit MD, MC with Maxwell-like boundary conditions based on the accommodation coefficients, and pure Maxwell boundary conditions are discussed. A map of the heat flux dependence on the accommodation coefficients for a dense gas, and the effective accommodation coefficients for different gas-wall interactions are given. In the end, this approach is applied to study the gas-surface interactions of argon and xenon molecules on a platinum surface. The derived accommodation coefficients are compared with values of experimental results.

A DFT based equilibrium study on the hydrolysis and the dehydration reactions of MgCl2 hydrates.

Magnesium chloride hydrates are characterized as promising energy storage materials in the built-environment. During the dehydration of these materials, there are chances for the release of harmful HCl gas, which can potentially damage the material as well as the equipment. Hydrolysis reactions in magnesium chloride hydrates are subject of study for industrial applications. However, the information about the possibility of hydrolysis reaction, and its preference over dehydration in energy storage systems is still ambiguous at the operating conditions in a seasonal heat storage system. A density functional theory level study is performed to determine molecular structures, charges, and harmonic frequencies in order to identify the formation of HCl at the operating temperatures in an energy storage system. The preference of hydrolysis over dehydration is quantified by applying thermodynamic equilibrium principles by calculating Gibbs free energies of the hydrated magnesium chloride molecules. The molecular structures of the hydrates (n = 0, 1, 2, 4, and 6) of MgCl2 are investigated to understand the stability and symmetry of these molecules. The structures are found to be noncomplex with almost no meta-stable isomers, which may be related to the faster kinetics observed in the hydration of chlorides compared to sulfates. Also, the frequency spectra of these molecules are calculated, which in turn are used to calculate the changes in Gibbs free energy of dehydration and hydrolysis reactions. From these calculations, it is found that the probability for hydrolysis to occur is larger for lower hydrates. Hydrolysis occurring from the hexa-, tetra-, and di-hydrate is only possible when the temperature is increased too fast to a very high value. In the case of the mono-hydrate, hydrolysis may become favorable at high water vapor pressure and at low HCl pressure.

Parameterization of a reactive force field using a Monte Carlo algorithm.

Parameterization of a molecular dynamics force field is essential in realistically modeling the physicochemical processes involved in a molecular system. This step is often challenging when the equations involved in describing the force field are complicated as well as when the parameters are mostly empirical. ReaxFF is one such reactive force field which uses hundreds of parameters to describe the interactions between atoms. The optimization of the parameters in ReaxFF is done such that the properties predicted by ReaxFF matches with a set of quantum chemical or experimental data. Usually, the optimization of the parameters is done by an inefficient single-parameter parabolic-search algorithm. In this study, we use a robust metropolis Monte-Carlo algorithm with simulated annealing to search for the optimum parameters for the ReaxFF force field in a high-dimensional parameter space. The optimization is done against a set of quantum chemical data for MgSO4 hydrates. The optimized force field reproduced the chemical structures, the equations of state, and the water binding curves of MgSO4 hydrates. The transferability test of the ReaxFF force field shows the extend of transferability for a particular molecular system. This study points out that the ReaxFF force field is not indefinitely transferable.

Molecular dynamics study of the influence of wall-gas interactions on heat flow in nanochannels.

Especially at the nanometer scale interfaces play an important role. The effect of the wettability on the solid-liquid interface has already been studied with molecular dynamics. In this paper we study the dependence of wetting on the solid-gas interface for different density gases and investigate the influence of wetting on the heat transport properties over such an interface using molecular dynamics. Subsequently we show how the flow profile of a gas flowing along a surface also depends on this wettability. These simulations show that wettability increases the conductivity of a solid to a stationary gas and decreases the flow velocity near the interface for a gas flow. These two effects influence the cooling of a solid achieved by a cold gas flowing along its surface in opposite ways. However, we show that a higher wettability has a positive net effect on the cooling, explaining experimental results that showed an increased heat cooling effect of hydrophilic over hydrophobic microchannels.

Hybrid method coupling molecular dynamics and Monte Carlo simulations to study the properties of gases in microchannels and nanochannels.

We combine molecular dynamics (MD) and Monte Carlo (MC) simulations to study the properties of gas molecules confined between two hard walls of a microchannel or nanochannel. The coupling between MD and MC simulations is introduced by performing MD near the boundaries for accuracy and MC in the bulk because of the low computational cost. We characterize the influence of different densities and molecule sizes on the equilibrium properties of the gas in the microchannel. The effect of the particle size on the simulation results is very small in the case of a dilute gas and increases with the density. The hybrid MD-MC simulation method is validated by comparing the results for density and temperature profiles with those of pure MD and pure MC simulations. These results compare well for pure MD and pure MC, as well as hybrid MD-MC, both in the bulk and near the boundaries, when hard-sphere interactions are used. When Lennard-Jones potentials are used to accurately model the interactions between the gas and wall molecules instead, the results of pure MD simulations differ significantly from the pure MC simulations near the boundaries, but the results of the hybrid method compare well with the pure MD results near the wall, and with the pure MC and pure MD results in the middle of the channel. The hybrid method also very accurately simulates the interface between the MD and MC simulation domains. Comparisons between MD, MC, and hybrid MD-MC computational costs are outlined. The speedup when using 50% of the domain for MD simulations and 50% for MC simulations is very small compared to pure MD simulations times, but this speedup increases drastically for more realistic situations where the region near the wall is small compared to the bulk region.

Infinitely fast diffusion in single-file systems.

We have used dynamic Monte Carlo(DMC) methods and analytical techniques to analyze single-file systems for which diffusion is infinitely fast. We have simplified the master equation removing the fast reactions, and we have introduced a DMC algorithm for infinitely fast diffusion. The DMC method for fast diffusion give similar results as the standard DMC with high diffusion rates. We have investigated the influence of characteristic parameters, such as pipe length, adsorption, desorption, and conversion rate constants on the steady-state properties of single-file systems with a reaction, looking at cases when all the sites are reactive and when only some of them are reactive. We find that the effect of fast diffusion on single-file properties of the system is absent even when diffusion is infinitely fast. Diffusion is not important in these systems. Smaller systems are less reactive and the occupancy profiles for infinitely long systems show an exponential behavior.

Analytical results for the reactivity of a single-file system.

We derive analytical expressions for the reactivity of a single-file system with fast diffusion and particles entering and leaving the system at one end. If the conversion reaction is fast, then the reactivity depends only very weakly on the system size, and the conversion is about 100%. If the reaction is slow, then the reactivity becomes proportional to the system size, the loading, and the reaction rate constant. If the system size increases the reactivity goes to the geometric mean of the reaction rate constant and the rate of particles entering and leaving the system. For large systems, the number of unconverted particles decreases exponentially with distance from the open end.

Transient behavior in single-file systems.

We have used Monte Carlo methods and analytical techniques to investigate the influence of the characteristics, such as pipe length, diffusion, adsorption, desorption, and reaction rates on the transient properties of single-file systems. The transient or the relaxation regime is the period in which the system is evolving to equilibrium. We have studied the system when all the sites are reactive and also when only some of them are reactive. Comparisons between mean-field predictions, cluster approximation predictions, and Monte Carlo simulations for the relaxation time of the system are shown. We outline the cases where the mean-field analysis gives good results compared to dynamic Monte Carlo results. For some specific cases we can analytically derive the relaxation time. Occupancy profiles for different distributions of the sites both for the mean field and simulations are compared. Different results for slow and fast reaction systems and different distributions of reactive sites are discussed.