Molecular Dynamic Study of Local Interfacial Thermal Resistance of Solid-Liquid and Solid-Solid Interfaces: Water and Nanotextured Surface

Yoshitaka Ueki,Satoshi Matsuo,Masahiko Shibahara

DOI: https://doi.org/10.48550/arXiv.2011.03184

2021-04-13

Abstract:Degradation in performances of air conditioners and refrigerators is caused by a frost formation and adhesion on the surface. In the present study, by means of the classical molecular dynamics simulation, we investigate how and how much the nanotextured surface characteristics, such as surface wettability and geometry, influenced the interfacial thermal resistance (ITR) between the solid wall and the water/ice. The ITR of the interfacial region was comparable in both the water and the ice states. As the nanostructure gaps became narrower, the ITR of the interfacial region decreased. The local ITR had a weak negative correlation with the local H2O molecule density regardless of the phase of the H2O molecules. The local ITR decreased as the local density increased. The greater amount of the thermal energy was transferred through the material interface by means of the intermolecular interaction when more the H2O molecules were located in the proximity area, which was closer to the Pt solid wall than the first adsorption layer peak. When the H2O molecules were in the crystal form on the solid wall, the proximity molecules decreased, and then the local ITR significantly increased.

Mesoscale and Nanoscale Physics,Applied Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **To explore the influence of nanotextured surface properties (such as wettability and geometry) on the solid - liquid and solid - solid interfacial thermal resistance (ITR), especially in the states of water and ice**. Specifically, the research aims to understand how these surface properties affect the heat transfer performance in the interface region through classical molecular dynamics simulations.

### Problem Background:
1. **Performance degradation of air conditioners and refrigerators**: Due to the formation and adhesion of frost on the surface, the performance of air conditioners and refrigerators declines.
2. **Defrosting requirement**: When frost forms, the heat - transfer surface must be heated to melt the frost, which increases energy consumption.
3. **Anti - icing requirement**: In order to improve cooling performance and efficiency, it is necessary to suppress the formation of frost and the adhesion of ice.

### Research Objectives:
1. **Understand the influence of nanotextured surfaces**: Study how the geometry and wettability of nanotextured surfaces affect the solid - liquid and solid - solid interfacial thermal resistance.
2. **Explore the variation mechanism of interfacial thermal resistance**: Analyze the variation law of interfacial thermal resistance under different conditions (such as different nanostructure gap widths, different water molecule densities, etc.).
3. **Optimize the design of heat - transfer surfaces**: Provide a theoretical basis for designing more efficient heat - transfer surfaces, thereby reducing the formation of frost and the adhesion of ice and improving device performance.

### Main Findings:
1. **Influence of nanostructure gaps**: As the nanostructure gaps become narrower, the interfacial thermal resistance decreases. This is because more water molecules are close to the solid wall, enhancing the intermolecular interaction and promoting heat transfer.
2. **Relationship between local thermal resistance and density**: The local thermal resistance is negatively correlated with the local water molecule density, that is, the larger the local density, the smaller the thermal resistance.
3. **Influence of crystalline and amorphous states**: When water molecules are in the crystalline state, the number of molecules close to the solid wall decreases, resulting in a significant increase in local thermal resistance.

### Formula Representation:
The formulas involved in the paper include the calculation of intermolecular potential energy, for example:

- Potential energy function of the mW model:
\[
\phi_{ij} = \phi(r_{ij}) + \sum_{k \neq i, j} \phi_{ijk}(r_{ij}, r_{ik}, \theta_{ijk})
\]
where,
\[
\phi(r_{ij}) = \epsilon \left[ \left( \frac{\sigma}{r_{ij}} \right)^{12} - \left( \frac{\sigma}{r_{ij}} \right)^6 \right] e^{-\left(\frac{r_{ij} - a}{\lambda}\right)^2}
\]
\[
\phi_{ijk}(r_{ij}, r_{ik}, \theta_{ijk}) = \gamma \left( \frac{\sigma}{r_{ij}} \right)^3 \left( \frac{\sigma}{r_{ik}} \right)^3 \left[ \cos(\theta_{ijk}) - \cos(\theta_0) \right]^2 e^{-\left(\frac{r_{ij} - a}{\lambda}\right)^2} e^{-\left(\frac{r_{ik} - a}{\lambda}\right)^2}
\]

- Lennard - Jones potential energy function:
\[
\phi_{ij} = 4 \epsilon \left[ \left( \frac{\sigma}{r_{ij}} \right)^{12} - \left( \frac{\sigma}{r_{ij}} \right)^6 \right]
\]
For the interaction between Pt and H₂O molecules, the potential energy function is multiplied by a parameter \(\alpha\):
\[
\phi_{ij} = 4 \alpha \epsilon \left[ \left( \frac{\sigma}{r_{ij}} \right)^{12} - \left( \frac{\sigma}{r_{ij}} \right)^6 \right]