Free energy studies of freezing in slit pores: an order-parameter approach using Monte Carlo simulation

The Nanoscale Wetting Parameter and Its Role in Interfacial Phenomena: Phase Transitions in Nanopores

Through analysis of the statistical mechanical equations for a thin adsorbed film (gas, liquid, or solid) on a solid substrate or confined within a pore, it is possible to express the equilibrium thermodynamic properties of the film as a function of just two dimensionless parameters: a nanoscale wetting parameter, αw, and pore width, H*. The wetting parameter, αw, is defined in terms of molecular parameters for the adsorbed film and substrate and so is applicable at the nanoscale and for films of any phase. The main assumptions in the treatment are that (a) the substrate structure is not significantly affected by the adsorbed layer and (b) the diameter of the adsorbate molecules is not very small compared to the spacing of atoms in the solid substrate. We show that different surface geometries of the substrate (e.g., slit, cylindrical, and spherical pores) and various models of wall heterogeneity can be accounted for through a well-defined correction to the wetting parameter; no new dimensionless variables are introduced. Experimental measurements are reported for contact angles for various liquids on several planar substrates and are shown to be closely correlated with the nanoscale wetting parameter. We apply this approach to phase separation in nanopores of various geometries. Molecular simulation results for the phase diagram in confinement, obtained by the flat histogram Monte Carlo method, are reported and are shown to be closely similar to experimental results for capillary condensation, melting, and the triple point. The value of the wetting parameter, αw, is shown to determine the qualitative behavior (e.g., increase vs decrease in the melting temperature, capillary condensation vs evaporation), whereas the pore width determines the magnitude of the confinement effect. The triple point temperature and pressure for the confined phase are always lower than those for the bulk phase for all cases studied.

Effects of Pore Structure and Pore Size of Expanded Graphite on the Properties of Paraffin Wax/Expanded Graphite Composite Phase Change Materials

In order to study the effects of the pore structure and pore size of porous carrier material on the phase transformation behavior of composite phase change materials (CPCMs), the pore structure and pore sizes of three different sizes of expanded graphites (EGs) (50, 80, 100 mesh) were studied using N2 adsorption-desorption isotherms and scanning electron microscopy. Then, the thermal characterization of CPCMs prepared with paraffin wax and EGs were tested using differential scanning calorimetry and a thermal conductivity tester. The results showed that EG-50 had a more web-like pore structure, and thus, higher adsorption capacity for paraffin wax. The addition of EG could reduce the supercooling degree of CPCMs and improve the thermal conductivity of CPCMs. CPCM with EG-50 had better performance due to its large specific surface area and low interfacial resistance. Compared with paraffin wax, the phase transition temperature (ΔT) of CPCMs increased slightly and the latent heat of CPCMS decreased to varying degrees. As the pore size of EG decreased, its constraint on PCM increased, but ΔT of CPCMs decreased, which was due to the combination of phase transformation behavior of different components in paraffin wax, which violates the conventional change law. It could be seen that the phase change behavior of CPCMs was related not only to the pore size of EG but also to the composition of PCM.

Study of Capric-Palmitic Acid/Clay Minerals as Form-Stable Composite Phase-Change Materials for Thermal Energy Storage

As an important method to effectively improve energy efficiency, the study of thermal energy storage is particularly important. In this study, six types of clay mineral-based form-stable phase-change materials (FSPCMs) were prepared by the vacuum adsorption method. The adsorption capacity of vermiculite and diatomite was satisfactory, and sepiolite showed the worst adsorption capacity. Clay minerals can delay the thermal decomposition rate of capric-palmitic acid (CA-PA), and the specific surface area and pore capacity of clay minerals all decrease dramatically after the clay minerals have absorbed CA-PA. FSPCMs exhibited a higher heat storage and release efficiency and reflects a certain temperature control performance. In addition, only physical adsorption between the CA-PA and the clay material occurred, and no chemical reaction occurred. Finally, FSPCMs still have high latent heat of phase transition, and they can be used for low-temperature thermal energy storage.

Bio-Based Phase Change Materials Incorporated in Lignocellulose Matrix for Energy Storage in Buildings—A Review

Due to growing consciousness regarding the environmental impact of fossil-based and non-sustainable materials in construction and building applications, there have been an increasing interest in bio-based and degradable materials in this industry. Due to their excellent chemical and thermo-physical properties for thermal energy storage, bio-based phase change materials (BPCMs) have started to attract attention worldwide for low to medium temperature applications. The ready availability, renewability, and low carbon footprint of BPCMs make them suitable for a large spectrum of applications. Up to now, most of the BPCMs have been incorporated into inorganic matrices with only a few attempts to set the BPCMs into bio-matrices. The current paper is the first comprehensive review on BPCMs incorporation in wood and wood-based materials, as renewable and sustainable materials in buildings, to enhance the thermal mass in the environmentally-friendly buildings. In the paper, the aspects of choosing BPCMs, bio-based matrices, phase change mechanisms and their combination, interpretation of life cycle analyses, and the eventual challenges of using these materials are presented and discussed.

Synthesis of monolithic shape-stabilized phase change materials with high mechanical stability via a porogen-assisted in situ sol–gel process

The confinement of phase change materials (PCMs) in construction materials has recently solved leakage, supercooling and low thermal conductivity problems in the industrial use of PCMs as monolithic thermal energy storage materials. To produce shape-stabilized PCMs (ss-PCMs) as crack-free monoliths, less than 15–30% v/v pure or encapsulated PCMs can be used in construction materials. Therefore, the heat storage capacity of these monolithic ss-PCM boards is comparatively low. In this study, we synthesized a novel class of monolithic ss-PCM boards with high compressive strength of 0.7 MPa at 30 °C (1.2 MPa at 10 °C), high PCM loadings of 86 wt%, and latent heats in the range of 100 J g⁻¹via a porogen-assisted in situ sol–gel process. We confined butyl stearate (BS) as PCM in a core-shell-like silica matrix via stabilized silica sol as silica source, sodium dodecyl sulfate as surfactant and poly(vinyl alcohol) as co-polymer. The ss-PCMs obtained are hydrophobic, thermally stable up to 320 °C and perform 6000 state transitions from solid to liquid and vice versa, without losing melting or freezing enthalpies. We analyzed the silica structure in the ss-PCMs to understand in detail the reasons for the high mechanical stability. The silica structure in ss-PCMs consists of spherical meso- and macropores up to 10 000 nm filled with PCM, formed mostly by BS droplets in water as templates during gelation. With an increasing BS amount in the synthesis of ss-PCMs, the total nanopore volume filled with PCM in ss-PCMs increases, resulting in higher compressive strengths up to 500% and thermal conductivities up to 60%.

Nanoconfinement of PCMs in core-shell-like silica structures via an inexpensive porogen-assisted sol–gel process to produce shape-stabilized PCMs as monoliths with high mechanical stability and high loading capacity.

Applications of Sustainable Polymer-Based Phase Change Materials in Mortars Composed by Different Binders

Eco-sustainable, low toxic and low flammable poly-ethylene glycol (PEG) was forced into flakes of the porous Lecce stone (LS), collected as stone cutting wastes, employing a very simple cheap method, to produce a “form-stable” phase change material (PCM). The experimental PCM was included in mortars based on different binders (hydraulic lime, gypsum and cement) in two compositions. The main thermal and mechanical characteristics of the produced mortars were evaluated in order to assess the effects due to the incorporation of the PEG-based PCM. The mortars containing the PEG-based PCM were found to be suitable as thermal energy storage systems, still displaying the characteristics melting and crystallization peaks of PEG polymer, even if the related enthalpies measured on the mortars were appreciably reduced respect to pure PEG. The general reduction in mechanical properties (in flexural and compressive mode) measured on all the mortars, brought about by the presence of PEG-based PCM, was overcome by producing mortars possessing a greater amount of binder. The proposed LS/PEG composite can be considered, therefore, as a promising PCM system for the different mortars analyzed, provided that an optimal composition is identified for each binder.

Evolution of methane density during melting in nanopores

Phase properties of gases adsorbed in small nanopores are mainly determined by the pore size and shape as well as the structural heterogeneity of the adsorbate. Here we analyze the evolution of the melting mechanism that occurs in pores <3 nm in size. Melting in slit-shaped graphene pores is compared with melting in SURMOF channel pores with square cross-sections. We show how the melting transformation is related to the adsorption mechanism. We use a graphical representation of the evolution of molecular density as a function of temperature in the nanopores.

Macro-Encapsulated PCM Cylinder Module Based on Paraffin and Float Stones

Organic phase change material (PCM) with macro-encapsulation is attractive in energy storage applications as it has relatively low cost. This study focuses on using PET plastic pipes to encapsulate paraffin and using low-cost float stones to increase the thermal conductivity of PCM modules as they have a special structure of high porosity. Float stones were immersed in the liquid PCM and an ultrasonic welding method used to prevent leakage of the PET plastic pipes. Scanning electron microscopy (SEM) was used to discover the appearance of the composite PCM. The thermal performance of the PCM cylinder module was analyzed through experimental tests of a constant-temperature water bath and numerical simulations. The result indicates that this PCM Ccylinder module is superior in thermal energy storage compared with the reference module even though fewer PCM was contained and the latent heat loss is considerable. The pipe diameter is an important parameter when using this kind of PCM cylinder module in water tanks.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcyclotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium with a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is computed for a specific confinement using the two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self-diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late α-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q4, was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

Diverting the phase transition behaviour of adipic acid via mesoporous silica confinement

Phase transition regulation was achieved by imbibing adipic acid into mesoporous silica chambers, on the basis of pore size-dependent thermal behaviour.

Examination of the phase transition behavior of nano-confined fluids by statistical temperature molecular dynamics

The fluid-solid phase transition behavior of nano-confined Lennard-Jones fluids as a function of temperature and degree of nanoconfinement has been studied via statistical temperature molecular dynamics (STMD). The STMD method allows the direct calculation of the density of states and thus the heat capacity with high efficiency. The fluids are simulated between parallel solid surfaces with varying pore sizes, wall-fluid interaction energies, and registry of the walls. The fluid-solid phase transition behavior has been characterized through determination of the heat capacity. The results show that for pores of ideal-spacing, the order-disorder transition temperature (T(ODT)) is reduced as the pore size increases until values consistent with that seen in a bulk system. Also, as the interaction between the wall and fluid is reduced, T(ODT) is reduced due to weak constraints from the wall. However, for non-ideal spacing pores, quite different behavior is obtained, e.g., generally T(ODT) are largely reduced, and T(ODT) is decreased as the wall constraint becomes larger. For unaligned walls (i.e., whose lattices are not in registry), the fluid-solid transition is also detected as T is reduced, indicating non-ideality in orientation of the walls does not impact the formation of a solid, but results in a slight change in T(ODT) compared to the perfectly aligned systems. The STMD method is demonstrated to be a robust way for probing the phase transitions of nanoconfined fluids systematically, enabling the future examination of the phase transition behavior of more complex fluids.

Liquid-solid and solid-solid phase transition of monolayer water: high-density rhombic monolayer ice

Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.

Development of Composite PCMs by Incorporation of Paraffin into Various Building Materials

In this research, we focused on the development of composite phase-change materials (CPCMs) by incorporation of a paraffin through vacuum impregnation in widely used building materials (Kaolin and ground granulated blast-furnace slag (GGBS)). The composite PCMs were characterized using environmental scanning electron microscopy (ESEM), Fourier transform infrared spectroscopy (FT-IR), differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) techniques. Moreover, thermal performance of cement paste composite PCM panels was evaluated using a self-designed heating system. Test results showed that the maximum percentage of paraffin retained by Kaolin and GGBS was found to be 18% and 9%, respectively. FT-IR results show that CPCMs are chemically compatible. The phase-change temperatures of CPCMs were in the human comfort zone, and they possessed considerable latent-heat storage capacity. TGA results showed that CPCMs are thermally stable, and they did not show any sign of degradation below 150 °C. From thermal cycling tests, it was revealed that the CPCMs are thermally reliable. Thermal performance tests showed that in comparison to the control room model, the room models prepared with CPCMs reduced both the temperature fluctuations and maximum indoor center temperature. Therefore, the prepared CPCMs have some potential in reducing peak loads in buildings when applied to building facade.

Effect of confinement on the solid-liquid coexistence of Lennard-Jones fluid

The solid-liquid coexistence of a Lennard-Jones fluid confined in slit pores of variable pore size, H, is studied using molecular dynamics simulations. Three-stage pseudo-supercritical transformation path of Grochola [J. Chem. Phys. 120(5), 2122 (2004)] and multiple histogram reweighting are employed for the confined system, for various pore sizes ranging from 20 to 5 molecular diameters, to compute the solid-liquid coexistence. The Gibbs free energy difference is evaluated using thermodynamic integration method by connecting solid-liquid phases under confinement via one or more intermediate states without any first order phase transition among them. Thermodynamic melting temperature is found to oscillate with wall separation, which is in agreement with the behavior seen for kinetic melting temperature evaluated in an earlier study. However, thermodynamic melting temperature for almost all wall separations is higher than the bulk case, which is contrary to the behavior seen for the kinetic melting temperature. The oscillation founds to decay at around H = 12, and beyond that pore size dependency of the shift in melting point is well represented by the Gibbs-Thompson equation.

Effect of pressure on the freezing of pure fluids and mixtures confined in nanopores

Monte Carlo simulations combined with the parallel tempering technique are used to study the freezing of Ar, CH(4), and their mixtures in a slit graphite nanopore. For all systems, the solid/liquid coexistence line is located at higher temperature than that for the bulk phase, as expected for fluids for which the wall/fluid interaction is stronger than the fluid/fluid interaction. In the case of the mixtures, the phase diagram for the confined system is of the same type as that for the bulk (azeotropic). It is also found that the freezing temperatures for the confined fluids and mixture are much more affected by pressure than those for the bulk phase. By calculating the isothermal compressibility of the confined fluids and determining the slope of the solid/liquid coexistence line (T,P) from the Clapeyron equation, we show that such a strong effect of pressure is not related to reduced compressibility within the pores. On the other hand, the pressure dependence of the in-pore freezing temperature is correctly described in the frame of the model proposed by Miyahara et al. [ Miyahara , M. ; Kanda , H. ; Shibao , M. ; Higashitani , K. J. Chem. Phys. 2000 , 112 , 9909. ], which is based on the pressure difference between the bulk and confined phases (capillary effect). In this model, a change in the in-pore freezing temperature with pressure is explained by a drastic change in the in-pore pressure, which varies very sharply with the bulk external pressure. We present an extended version of this model to confined systems for which an increase in the freezing temperature is observed.

Freezing in the bulk controlled by prefreezing at a surface

We use Monte Carlo simulations of the Lennard-Jones model to study the nucleation of a crystal phase at a flat surface. Our motivation is the observation that crystal phases almost always nucleate at a surface. We find that a surface phase transition (prefreezing) can control nucleation of the bulk crystal. This finding should be general and so surface phase behavior should be considered whenever nucleation of bulk phases at surfaces is considered. Also, nucleation of the bulk crystal transforms smoothly into the nucleation of a surface crystal layer as the bulk transition is crossed.

Structure of liquid and glassy methanol confined in cylindrical pores

We present a neutron scattering analysis of the density and the static structure factor of confined methanol at various temperatures. Confinement is performed in the cylindrical pores of MCM-41 silicates with pore diameters D=24 and 35 A. A change of the thermal expansivity of confined methanol at low temperature is the signature of a glass transition, which occurs at higher temperature for the smallest pore. This is evidence of a surface induced slowing down of the dynamics of the fluid. The structure factor presents a systematic evolution with the pore diameter, which has been analyzed in terms of excluded volume effects and fluid-matrix cross correlation. Conversely to the case of Van der Waals fluids, it shows that stronger fluid-matrix correlations must be invoked most probably in relation with the H-bonding character of both methanol and silicate surface.

Solvation force, structure and thermodynamics of fluids confined in geometrically rough pores

The effect of periodic surface roughness on the behavior of confined soft sphere fluids is investigated using grand canonical Monte Carlo simulations. Rough pores are constructed by taking the prototypical slit-shaped pore and introducing unidirectional sinusoidal undulations on one wall. For the above geometry our study reveals that the solvation force response can be phase shifted in a controlled manner by varying the amplitude of roughness. At a fixed amplitude of roughness, a, the solvation force for pores with structured walls was relatively insensitive to the wavelength of the undulation, lambda for 2.3<lambda/sigma(ff)<7, where sigma(ff) is the Lennard-Jones diameter of the confined fluid. This was not the case for smooth walled pores, where the solvation force response was found to be sensitive to the wavelength, for lambda/sigma(ff)<7.0 and amplitudes of roughness, a/sigma(ff)>/=0.5. The predictions of the superposition approximation, where the solvation force response for the rough pores is deduced from the solvation force response of the slit-shaped pores, was in excellent agreement with simulation results for the structured pores and for lambda/sigma(ff)>/=7 in the case of smooth walled pores. Grand potential computations illustrate that interactions between the walls of the pore can alter the pore width corresponding to the thermodynamically stable state, with wall-wall interactions playing an important role at smaller pore widths and higher amplitudes of roughness.

Triple point of Lennard-Jones fluid in slit nanopore: Solidification of critical condensate

We report the results of a molecular dynamics simulation that looked for the triple point of Lennard-Jones fluid in slit-shaped nanopores. The simulation method employed for this purpose is able to maintain vapor-liquid coexistence in a nanopore at a specific equilibrium bulk-phase pressure. The triple point is the freezing point of the critical condensate. The triple-point temperature could be higher or lower than the bulk triple point, depending on the pore size. This is thought to be due to two opposing factors: the elevating effect of the pore-wall potential energy, and the depressing effect of the capillary condensate's tensile condition. Because of the cancellation, the deviation of the triple-point temperature from the bulk triple-point temperature was not considered significant. The pressure of the triple point, however, was significantly different from that of the bulk triple point. A simple model to describe the triple point is developed and shown to agree well with the results of the simulation. The importance of the two factors in nanoscale pores, which cannot be described by the classic Gibbs-Thomson equation, is emphasized.