Continuum modeling and simulation of fluids at nanoscalePopadić, Aleksandar (Avtor)
Praprotnik, Matej (Mentor)
fullerenescarbon nanotubespolymer meltscomputational fluid dynamicsfinite volume methodpartial slipmolecular dynamicsMonte CarloNanofluidics, a scientific field exploring nanoconfined fluids, has seen increasing interest in recent years. This interest is in large part due to the fascinating phenomena that occur at such small scales, such as ultra-efficient fluid transport through tubes with diameters of the order of nanometers. The development of computer science and technology in recent years has also allowed for efficient studies of matter at nanoscale. In this thesis, we study fluids at nanoscale with the aid of computer simulations. The thesis is roughly divided into two parts. In the first part, we study fluid flow past and through objects with sizes of the order of nanometers, and in the second part, we study conservation equations for polymer chains and the associated coupling of density and directional fluctuations.
A widely used approach for studying nanoconfined fluids is molecular dynamics. It is, however, computationally expensive, often mandating the use of supercomputers. Hence, in this thesis we resort to the use of a continuum approach. Various recent studies have suggested that Navier-Stokes equations could be valid down to nanometer scale. However, at such small scales, where the ratio between the surface and bulk volume is not small anymore, the physics at the surface becomes increasingly important. The effective no-slip boundary condition, which is the standard boundary condition at macroscale, is no longer valid. At nanoscale, a slip between the fluid and the surface is observed. We thus use the Navier boundary condition, which takes the slip into account by assuming a linear dependence of shear stress at the surface on the relative velocity between the fluid and the solid at the interface. The boundary condition is parametrized by the slip length, which for straight walls represents the depth to which the fluid velocity profile must be extrapolated to vanish.
We perform computational fluid dynamics simulations, subject to the Navier boundary condition, past spherical molecules and through carbon nanotubes. Water flow through carbon nanotubes is subjected to significant slip. This results in energy efficient flow through the carbon nanotubes. Due to this energy efficiency, energy dissipation at the nanotube entrance and exit becomes significant. We thus explore the energy dissipation in the vicinity of carbon nanotube ends. We observe a nonmonotonic dependence of energy dissipation on the slip length. We successfully model and explain the origin of this dependence and accurately predict the optimal slip length at which the energy dissipation is at its minimum.
In light of recent reports of phonon modes of carbon nanotubes and subsequent diffusion enhancements, we explore the influence of the oscillating carbon nanotubes on the water flow. To examine the effect of carbon nanotube oscillations, we employ fluctuating hydrodynamics, in which a random fluctuating stress tensor representing thermal fluctuations is introduced to the Navier-Stokes equations. We solve the equations and derive the diffusion coefficient for the center of mass of water in an oscillating carbon nanotube. We show that in the continuum description, the phonon modes of the carbon nanotubes do not contribute to the diffusion.
In the second part of this thesis, we study the conservation equations for polymer chain melts. In polymer chain melts, the defects in orientational order are closely related to density. This connection is described by the continuity equation, which is polar in its nature. This presents a problem for nematic polymers as the direction cannot be uniquely defined. In case of long polymer chains with abundant chain folding, the polar order disappears while the nematic order is conserved. By performing Monte Carlo simulations, we show that the continuity equation can be applied if we define a “recovered polar order” by introducing chain cuts at chain fold positions. In the last part, we examine the tensorial conservation equation, which is based on the quadrupolar order. In case of polymer chain folding, the quadrupolar order, in contrast to the polar order, is not lost, thus making the tensorial conservation more appropriate compared to the vectorial conservation equation. We perform Monte Carlo simulations of polymer chains in an istoropic phase and show first evidence of presence of the tensorial constraint.20182018-06-01 07:45:16Doktorsko delo/naloga101420VisID: 87306COBISS_ID: 3207780sl