BoHung Kim

Transport Phenomena of Water in Molecular Fluidic Channels

In molecular-level fluidic transport, where the discrete characteristics of a molecular system are not negligible (in contrast to a continuum description), the response of the molecular water system might still be similar to the continuum description if the time and ensemble averages satisfy the ergodic hypothesis and the scale of the average is enough to recover the classical thermodynamic properties. However, even in such cases, the continuum description breaks down on the material interfaces. In short, molecular-level liquid flows exhibit substantially different physics from classical fluid transport theories because of (i) the interface/surface force field, (ii) thermal/velocity slip, (iii) the discreteness of fluid molecules at the interface and (iv) local viscosity. Therefore, in this study, we present the result of our investigations using molecular dynamics (MD) simulations with continuum-based energy equations and check the validity and limitations of the continuum hypothesis. Our study shows that when the continuum description is subjected to the proper treatment of the interface effects via modified boundary conditions, the so-called continuum-based modified-analytical solutions, they can adequately predict nanoscale fluid transport phenomena. The findings in this work have broad effects in overcoming current limitations in modeling/predicting the fluid behaviors of molecular fluidic devices.

Pressure dependence of Kapitza resistance at gold/water and silicon/water interfaces

We conducted non-equilibrium molecular dynamics simulations to investigate Kapitza length at solid/liquid interfaces under the effects of bulk liquid pressures. Gold and silicon were utilized as hydrophilic and hydrophobic solid walls with different wetting surface behaviors, while the number of confined liquid water molecules was adjusted to obtain different pressures inside the channels. The interactions of solid/liquid couples were reparameterized accurately by measuring the water contact angle of solid substrates. In this paper, we present a thorough analysis of the structure, normal stress, and temperature distribution of liquid water to elucidate thermal energy transport across interfaces. Our results demonstrate excellent agreement between the pressures of liquid water in nano-channels and published thermodynamics data. The pressures measured as normal stress components were characterized using a long cut-off distance reinforced by a long-range van der Waals tail correction term. To clarify the effects of bulk liquid pressures on water structure at hydrophilic and hydrophobic solid surfaces, we defined solid/liquid interface spacing as the distance between the surface and the peak value of the first water density layer. Near the gold surface, we found that interface spacing and peak value of first water density layer were constant and did not depend on bulk liquid pressure; near the silicon surface, those values depended directly upon bulk liquid. Our results reveal that the pressure dependence of Kapitza length strongly depends on the wettability of the solid surface. In the case of the hydrophilic gold surface, Kapitza length was stable despite increasing bulk liquid pressure, while it varied significantly at the hydrophobic silicon surface.

Near-surface viscosity effects on capillary rise of water in nanotubes.

In this paper, we present an approach for predicting nanoscale capillary imbibitions using the Lucas-Washburn (LW) theory. Molecular dynamics (MD) simulations were employed to investigate the effects of surface forces on the viscosity of liquid water. This provides an update to the modified LW equation that considered only a nanoscale slip length. An initial water nanodroplet study was performed to properly elucidate the wetting behavior of copper and gold surfaces. Intermolecular interaction strengths between water and corresponding solid surfaces were determined by matching the contact angle values obtained by experimental measurements. The migration of liquid water into copper and gold capillaries was measured by MD simulations and was found to differ from the modified LW equation. We found that the liquid layering in the vicinity of the solid surface induces a higher density and viscosity, leading to a slower MD uptake of fluid into the capillaries than was theoretically predicted. The near-surface viscosity for the nanoscale-confined water was defined and calculated for the thin film of water that was sheared between the two solid surfaces, as the ratio of water shear stress to the applied shear rate. Considering the effects of both the interface viscosity and slip length of the fluid, we successfully predicted the MD-measured fluid rise in the nanotubes.

Atomic density effects on temperature characteristics and thermal transport at grain boundaries through a proper bin size selection

This study focuses on the proper characterization of temperature profiles across grain boundaries (GBs) in order to calculate the correct interfacial thermal resistance (ITR) and reveal the influence of GB geometries onto thermal transport. The solid-solid interfaces resulting from the orientation difference between the (001), (011), and (111) copper surfaces were investigated. Temperature discontinuities were observed at the boundary of grains due to the phonon mismatch, phonon backscattering, and atomic forces between dissimilar structures at the GBs. We observed that the temperature decreases gradually in the GB area rather than a sharp drop at the interface. As a result, three distinct temperature gradients developed at the GB which were different than the one observed in the bulk solid. This behavior extends a couple molecular diameters into both sides of the interface where we defined a thickness at GB based on the measured temperature profiles for characterization. Results showed dependence on the selection of the bin size used to average the temperature data from the molecular dynamics system. The bin size on the order of the crystal layer spacing was found to present an accurate temperature profile through the GB. We further calculated the GB thickness of various cases by using potential energy (PE) distributions which showed agreement with direct measurements from the temperature profile and validated the proper binning. The variation of grain crystal orientation developed different molecular densities which were characterized by the average atomic surface density (ASD) definition. Our results revealed that the ASD is the primary factor affecting the structural disorders and heat transfer at the solid-solid interfaces. Using a system in which the planes are highly close-packed can enhance the probability of interactions and the degree of overlap between vibrational density of states (VDOS) of atoms forming at interfaces, leading to a reduced ITR. Thus, an accurate understanding of thermal characteristics at the GB can be formulated by selecting a proper bin size.

Wetting of chemically heterogeneous striped surfaces: Molecular dynamics simulations

Using molecular dynamics simulations, we thoroughly investigated the wetting behaviors of a chemically heterogeneous striped substrate patterned with two different wetting materials, face-centered cubic gold and face-centered cubic silver. We analyzed the density distributions, normal stress distributions, surface tensions, and contact angles of a water droplet placed on the substrates at different heterogeneities. We found that the density and stress profile of the water droplet near the substrate-water interface were noticeably affected by altering the gold and silver contents in the substrate. Specifically, a greater portion of gold (more wetting) or smaller portion of silver (less wetting) in the substrate composition induced higher densities and higher normal stresses in the vicinity of the substrate surface. Also, it was observed that the surface tensions at liquid-vapor interface and solid-vapor interface were not largely impacted by the change of the substrate composition while the solid-liquid surface tension decreased exponentially with increasing fraction of gold. Most importantly, we found that contact angle of a nanometer-sized water droplet resting on the chemically heterogeneous striped substrate does not show linear dependence on corresponding surface fractions like that predicted by Cassie-Baxter model at the macro-scale. Consequently, we proposed a method for successfully predicting the contact angle by including the critical effects of the substrate heterogeneity on both surface tensions and line tension at the three-phase contact line of the water droplet and the chemically striped substrate.

Physical origins of temperature continuity at an interface between a crystal and its melt

We justify and discuss the physical origins for the assumption of temperature continuity at crystal/melt interfaces by performing atomistic simulations. We additionally answer why the crystal/melt interfaces differ from the typical solid/liquid interfaces, which usually exhibit dissimilarities and a resulting temperature drop. We present results for pure silver modeled using the embedded-atom method and Lennard-Jones potential function and contrast the results with each other. We find that the temperature continuity at an interface between a crystal and its melt originates from the perfect vibrational coupling, which is caused by the interfacial structural diffusivity. This study provides fundamental insights into the heat transfer for cases of extremely large heat flux and thermal gradients occurring during rapid melting and solidification. The findings additionally determine the role of rough surfaces in manipulating the thermal conductance in nanodevices.