Molecular Dynamics Simulation in Arbitrary Geometries for Nanoscale Fluid Mechanics

Graham Bruce Macpherson

Research output: ThesisDoctoral Thesis

Abstract

Simulations of nanoscale systems where fluid mechanics plays an important role are required to help design and understand nano-devices and biological systems. A simulation method which hybridises molecular dynamics (MD) and continuum computational fluid dynamics (CFD) is demonstrated to be able to accurately represent the relevant physical phenomena and be computationally tractable. An MD code has been written to perform MD simulations in systems where the geometry is described by a mesh of unstructured arbitrary polyhedral cells that have been spatially decomposed into irregular portions for parallel processing. The MD code that has been developed may be used for simulations on its own, or may serve as the MD component of a hybrid method. The code has been implemented using OpenFOAM, an open source C++ CFD toolbox (www.openfoam.org). Two key enabling components are described in detail. 1) Parallel generation of initial configurations of molecules in arbitrary geometries. 2) Calculation of intermolecular pair forces, including between molecules that lie on mesh portions assigned to different, and possibly non-neighbouring processors. To calculate intermolecular forces, the spatial relationship of mesh cells is calculated once at the start of the simulation and only the molecules contained in cells that have part of their surface closer than a cut-off distance are required to interact. Interprocessor force calculations are carried out by creating local copies of molecules from other processors in a layer around the processor in question. The process of creating these copied molecules is described in detail. A case study of flow in a realistic nanoscale mixing channel, where the geometry is drawn and meshed using engineering CAD tools, is simulated to demonstrate the capabilities of the code for complex simulations.
LanguageEnglish
Awarding Institution
  • University Of Strathclyde
Publication statusPublished - 2008

Fingerprint

Fluid mechanics
Molecular dynamics
Molecules
Geometry
Computer simulation
Fluid dynamics
Biological systems
Computer aided design
Processing

Keywords

  • molecular dynamics
  • thesis
  • fluid mechanics
  • nanoscale fluid dynamics
  • nanoscale systems
  • computational fluid dynamics
  • mechanical engineering

Cite this

@phdthesis{45ac4b4eea9e4b0c91e9e10953d3eaed,
title = "Molecular Dynamics Simulation in Arbitrary Geometries for Nanoscale Fluid Mechanics",
abstract = "Simulations of nanoscale systems where fluid mechanics plays an important role are required to help design and understand nano-devices and biological systems. A simulation method which hybridises molecular dynamics (MD) and continuum computational fluid dynamics (CFD) is demonstrated to be able to accurately represent the relevant physical phenomena and be computationally tractable. An MD code has been written to perform MD simulations in systems where the geometry is described by a mesh of unstructured arbitrary polyhedral cells that have been spatially decomposed into irregular portions for parallel processing. The MD code that has been developed may be used for simulations on its own, or may serve as the MD component of a hybrid method. The code has been implemented using OpenFOAM, an open source C++ CFD toolbox (www.openfoam.org). Two key enabling components are described in detail. 1) Parallel generation of initial configurations of molecules in arbitrary geometries. 2) Calculation of intermolecular pair forces, including between molecules that lie on mesh portions assigned to different, and possibly non-neighbouring processors. To calculate intermolecular forces, the spatial relationship of mesh cells is calculated once at the start of the simulation and only the molecules contained in cells that have part of their surface closer than a cut-off distance are required to interact. Interprocessor force calculations are carried out by creating local copies of molecules from other processors in a layer around the processor in question. The process of creating these copied molecules is described in detail. A case study of flow in a realistic nanoscale mixing channel, where the geometry is drawn and meshed using engineering CAD tools, is simulated to demonstrate the capabilities of the code for complex simulations.",
keywords = "molecular dynamics, thesis, fluid mechanics, nanoscale fluid dynamics, nanoscale systems, computational fluid dynamics, mechanical engineering",
author = "Macpherson, {Graham Bruce}",
year = "2008",
language = "English",
school = "University Of Strathclyde",

}

Molecular Dynamics Simulation in Arbitrary Geometries for Nanoscale Fluid Mechanics. / Macpherson, Graham Bruce.

2008. 255 p.

Research output: ThesisDoctoral Thesis

TY - THES

T1 - Molecular Dynamics Simulation in Arbitrary Geometries for Nanoscale Fluid Mechanics

AU - Macpherson, Graham Bruce

PY - 2008

Y1 - 2008

N2 - Simulations of nanoscale systems where fluid mechanics plays an important role are required to help design and understand nano-devices and biological systems. A simulation method which hybridises molecular dynamics (MD) and continuum computational fluid dynamics (CFD) is demonstrated to be able to accurately represent the relevant physical phenomena and be computationally tractable. An MD code has been written to perform MD simulations in systems where the geometry is described by a mesh of unstructured arbitrary polyhedral cells that have been spatially decomposed into irregular portions for parallel processing. The MD code that has been developed may be used for simulations on its own, or may serve as the MD component of a hybrid method. The code has been implemented using OpenFOAM, an open source C++ CFD toolbox (www.openfoam.org). Two key enabling components are described in detail. 1) Parallel generation of initial configurations of molecules in arbitrary geometries. 2) Calculation of intermolecular pair forces, including between molecules that lie on mesh portions assigned to different, and possibly non-neighbouring processors. To calculate intermolecular forces, the spatial relationship of mesh cells is calculated once at the start of the simulation and only the molecules contained in cells that have part of their surface closer than a cut-off distance are required to interact. Interprocessor force calculations are carried out by creating local copies of molecules from other processors in a layer around the processor in question. The process of creating these copied molecules is described in detail. A case study of flow in a realistic nanoscale mixing channel, where the geometry is drawn and meshed using engineering CAD tools, is simulated to demonstrate the capabilities of the code for complex simulations.

AB - Simulations of nanoscale systems where fluid mechanics plays an important role are required to help design and understand nano-devices and biological systems. A simulation method which hybridises molecular dynamics (MD) and continuum computational fluid dynamics (CFD) is demonstrated to be able to accurately represent the relevant physical phenomena and be computationally tractable. An MD code has been written to perform MD simulations in systems where the geometry is described by a mesh of unstructured arbitrary polyhedral cells that have been spatially decomposed into irregular portions for parallel processing. The MD code that has been developed may be used for simulations on its own, or may serve as the MD component of a hybrid method. The code has been implemented using OpenFOAM, an open source C++ CFD toolbox (www.openfoam.org). Two key enabling components are described in detail. 1) Parallel generation of initial configurations of molecules in arbitrary geometries. 2) Calculation of intermolecular pair forces, including between molecules that lie on mesh portions assigned to different, and possibly non-neighbouring processors. To calculate intermolecular forces, the spatial relationship of mesh cells is calculated once at the start of the simulation and only the molecules contained in cells that have part of their surface closer than a cut-off distance are required to interact. Interprocessor force calculations are carried out by creating local copies of molecules from other processors in a layer around the processor in question. The process of creating these copied molecules is described in detail. A case study of flow in a realistic nanoscale mixing channel, where the geometry is drawn and meshed using engineering CAD tools, is simulated to demonstrate the capabilities of the code for complex simulations.

KW - molecular dynamics

KW - thesis

KW - fluid mechanics

KW - nanoscale fluid dynamics

KW - nanoscale systems

KW - computational fluid dynamics

KW - mechanical engineering

UR - http://www.mecheng.strath.ac.uk/multiscaleflows/files/GBM_thesis.pdf

M3 - Doctoral Thesis

ER -