Manipulating thermal resistance at the solid–fluid interface through monolayer deposition

Unraveling the Regimes of Interfacial Thermal Conductance at a Solid/Liquid Interface

The interfacial thermal conductance at a solid/liquid interface (G) exhibits an exponential-to-linear crossover with increasing solid/liquid interaction strength, previously attributed to the relative strength of solid/liquid to liquid/liquid interactions. Instead, using a simple Lennard-Jones setup, our molecular simulations reveal that this crossover occurs due to the onset of solidification in the interfacial liquid at high solid/liquid interaction strengths. This solidification subsequently influences interfacial energy transport, leading to the crossover in G. We use the overlap between the spectrally decomposed heat fluxes of the interfacial solid and liquid to pinpoint when "solid-like energy transport" within the interfacial liquid emerges. We also propose a novel decomposition of G into (i) the conductance right at the solid/liquid interface and (ii) the conductance of the nanoscale interfacial liquid region. We demonstrate that the rise of solid-like energy transport within the interfacial liquid influences the relative magnitude of these conductances, which in turn dictates when the crossover occurs. Our results can aid engineers in optimizing G at realistic interfaces, critical to designing effective cooling solutions for electronics among other applications.

Probing the thermal resistance of solid-liquid interfaces in nanofluids with molecular dynamics

The significance of interfacial thermal resistance in the thermal conductivity of nanofluids is not well understood, in part because of the absence of measurements of this quantity. Here, we study the interfacial thermal resistance for metal-oil nanofluids of interest as heat transfer fluids for concentrating solar power, using density functional theory and molecular dynamics simulations. Insights on the role of chemical interactions in determining the interfacial thermal resistance are revealed. The results presented here showcase a general picture in which the stronger the chemical interactions between species at the interface, the lower the associated interfacial thermal resistance. The implications toward nanofluid design are discussed. We show that, for this important family of metal-oil nanofluids, the interfacial thermal resistance values are low enough so that it is possible to afford a reduction in particle size, minimizing stability and rheological issues while still offering enhancement in the effective thermal conductivity with respect to the base fluid.

Solid-like Behaviors Govern Evaporative Transport in Adsorbed Water Nanofilms

The thermophysical attributes of water molecules confined in a sub-nanometer thickness significantly differ from those in bulk liquid where their molecular behaviors start governing interfacial physics at the nanoscale. In this study, we elucidate nanothin film evaporation by employing a computational approach from a molecular perspective. As the liquid thickness decreases, the solid-like characteristics of adsorbed water nanofilms make the resistance at solid-liquid interfaces or Kapitza resistance significant. Kapitza resistances not only show a strong correlation with the surface wettability but also dominate the overall thermal resistance during evaporation rather than the resistance at evaporating liquid-vapor interfaces. Once the liquid thickness reaches the critical value of 0.5-0.6 nm, the evaporation kinetics is suppressed due to the excessive forces between the liquid and solid atoms. The understanding of molecular-level behaviors explains how a hydrophilic surface plays a role in determining evaporation rates from an atomistic perspective.

Molecular transportation phenomena of simple liquids through a nanoporous graphene membrane

The mechanism of pressure-driven transport of simple liquid through a nanoporous graphene membrane has been analyzed using nonequilibrium molecular dynamics simulation. In this study, we investigate liquid dynamics properties such as local density, pressure variation, and local viscosity depending on the flow region. With movement of the specular reflection wall at the end of the front and back reservoirs, a pressure difference occurs mainly due to the change in the relative distance between the liquid molecules in the corresponding reservoir. The interfacial pressure difference strongly depends on the intermolecular force of the graphene membrane governed by the layered structure of the simple liquid and the applied flow velocity. The local viscosity was calculated for a nanochannel of simple liquid sheared by graphene walls. The liquid velocity adjacent to the pore edge was considered as the slip velocity, which provides updates in the Sampson flow equation. We observed that the entrance interfacial pressure and higher local viscosity in the vicinity of the graphene membrane, which are associated with the optimized definition of the wall-liquid boundary near the pore edge, play a critical role in the permeation of simple liquids through the nanoporous graphene membrane.

Thermal Assessment of Laminar Flow Liquid Cooling Blocks for LED Circuit Boards Used in Automotive Headlight Assemblies

This research work presents a comparative thermal performance assessment of the laminar flow cooling blocks produced for automotive headlight assembly using a high power Light Emitting Diode (LED) chip. A three-dimensional numerical model with conjugate heat transfer in solid and fluid domains was used. Laminar flow was considered in the present analysis. The validation of the numerical model was realized by using the measured data from the test rig. It was observed that substantial temperature variations were occurred around the LED chip owing to volumetric heat generation. The cooling board with lower height performs better thermal performance but higher pressure drop for the same mass flow rates. The cooling board with the finned cover plate performs better thermal performance but results in an increased pressure drop for the same mass flow rates. Increasing the power of the LED results in higher temperature values for the same mass flow rates. The junction temperature is highly dependent on the mass flow rates and LED power. It can be controlled by means of the mass flow rate of the coolant fluid. New Nusselt number correlations are proposed for laminar flow mini-channel liquid cooling block applications.

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young's modulus, Poisson's ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young's modulus and Poisson's ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a "safe zone" regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

Prediction of Kapitza resistance at fluid-solid interfaces

Understanding the interfacial heat transfer and thermal resistance at an interface between two dissimilar materials is of great importance in the development of nanoscale systems. This paper introduces a new and reliable linear response method for calculating the interfacial thermal resistance or Kapitza resistance in fluid-solid interfaces with the use of equilibrium molecular dynamics (EMD) simulations. The theoretical predictions are validated against classical molecular dynamics (MD) simulations. MD simulations are carried out in a Lennard-Jones (L-J) system with fluid confined between two solid slabs. Different types of interfaces are tested by varying the fluid-solid interactions (wetting coefficient) at the interface. It is observed that the Kapitza length decreases monotonically with an increasing wetting coefficient as expected. The theory is further validated by simulating under different conditions such as channel width, density, and temperature. Our method allows us to directly determine the Kapitza length from EMD simulations by considering the temperature fluctuation and heat flux fluctuations at the interface. The predicted Kapitza length shows an excellent agreement with the results obtained from both EMD and non-equilibrium MD simulations.

Subatomic-Level Solid/Fluid Boundary of Lennard-Jones Atoms: A Molecular Dynamics Study of Metal-Inert Fluid Interface

At the molecular scale, the definition of solid/fluid boundary is ambiguous since its defining precision is comparable to the size of the electron orbitals. It is important to figure out the sub-atomic-level solid/fluid boundary as the definition of the solid/fluid interface is related to estimating various properties such as slip length, Kapitza resistance, confined volume, thermodynamic properties, and material properties. In this work, molecular dynamics (MD) simulations were conducted to show the effects of the solid/fluid boundary on estimating thermodynamic properties. Our results reveal that the different definitions of solid/fluid boundary can cause a considerable impact on quantitative analysis and even qualitative analysis of a nanoscale system. The solid/fluid boundary for Lennard-Jones atoms is determined within sub-atomic precision via heat transfer MD simulations and microscopic heat flux relation. The result shows that solid/fluid boundary is slightly shifted to the fluid regime as the temperature increase. We suggested a mathematical expression of solid/fluid boundary of LJ atom that is theoretically estimated by ignoring the thermal vibration. The results presented in this work are expected to improve the accuracy of analyzing nanoscale phenomena as well as the continuum-based models for nanoscale heat and mass transport.