|Louisiana State University
|Free Energy; CFD; MD; Navier-Stokes; fluid flow; Multiscale Simulation
|Full text PDF:
A multiscale framework for simulation of fluid flow is developed by using a hybrid atomistic-continuum approach to couple molecular dynamics (MD) to the incompressible Navier-Stokes equations. The framework incorporates the accuracy of atomistic methods to resolve the detailed physics at the nano- to micro-scale in regions with large gradients, and the efficiency of the continuum models to describe the rest of the fluid. A state exchange scheme is employed where the value of conserved variables are reciprocated between the atomistic and continuum descriptions. A polyatomic constrained dynamics is introduced to facilitate construction of the required atomistic information from the corresponding continuum domain. The developed framework is demonstrated by studying flow field evolution in shear driven Couette flows, ranging from simple LennardJones liquids to complex ones such as water and alkanes. The molecular details of interfaces at solidfluid and fluidfluid under shear are investigated and compare them with results from analytical predictions.As a prelude to developing the ability to handle large nanoparticles within the multiscale framework, the interaction of nanoparticles in water with planar polymeric walls is investigated by means of molecular dynamics and free energy calculations. A mean-field approximation is employed to represent an apolar fullerene-based nanoparticle as a spherical colloid which enables continuous variation of the solute size. This allows for a systematic size-dependent study of the solute potentials of mean force (PMFs) and solvation forces as a function of solutewall distance. It is shown that the solute has a strong tendencyfor deposition onto the PMMA surface. However, for solutes with diameters greater than 3 nm a free energy barrier starts to hinder the solutesurface deposition. The transition from barrierless PMFs at small solute sizes to PMFs with barriers at larger solute sizes is attributed to the loss of hydrogen bonding between water molecules in proximity to the solute which intensifies the role of the comparatively weak, attractive van der Waals interactions between the water molecules and the solute, and the enlarged solute surface which makes depletion of the last layers of water trapped between the solute and the wall extremely expensive.Advisors/Committee Members: Nikitopoulos, Dimitris E. (chair), Moldovan, Dorel (committee member), Kohlmeyer, Axel (committee member), Ramanujam, Jagannathan (committee member), Novak, Brian (committee member), Chance, Donald M. (committee member).