A hybrid molecular-continuum simulation method for incompressible flows in micro/nanofluidic networks

Matthew K. Borg, Duncan A. Lockerby, Jason M. Reese

Research output: Contribution to journalArticle

22 Citations (Scopus)

Abstract

We present a hybrid molecular-continuum simulation method for modelling nano- and micro-flows in network-type systems. In these types of problem, a full molecular dynamics (MD) description of the macroscopic flow behaviour would be computationally intractable, or at least too expensive to be practical for engineering design purposes. Systems that exhibit multiscale traits, such as this, can instead be solved using a hybrid approach that distinguishes the problem into macroscopic and microscopic dynamics, modelled by their respective solvers. The technique presented in this study is an extension and addition to a hybrid method developed by Borg et al. (J Comput Phys 233:400-413, 2013) for high-aspect-ratio channel geometries, known as the internal-flow multiscale method (IMM). Computational savings are obtained by replacing long channels in the network, which are highly scale-separated, by much smaller, but representative, MD simulations, without a substantial loss of accuracy. On the other hand, junction components do not exhibit this length-scale separation, and so must be simulated in their entirety using MD. The current technique combines all network elements (junctions and channels) together in a coupled simulation using continuum conservation laws. For the case of steady, isothermal, incompressible, low-speed flows, we use the conservation of mass and momentum flux equations to derive a set of molecular-continuum constraints. An algorithm is presented here that computes at each iteration the new constraints on the pressure differences to be applied over individual MD micro-elements (channels and junctions), successively moving closer to macroscopic mass and momentum conservation. We show that hybrid simulations of some example network cases converge quickly, in only a few iterations, and compare very well to the corresponding full MD results, which are taken as the most accurate solutions. Major computational savings can be afforded by the IMM-type approximation in the channel components, but for steady-state solutions, even greater savings are possible. This is because the micro-elements are coupled to a steady-state continuum conservation expression, which greatly speeds up the relaxation of individual micro-components to steady conditions as compared to that of a full MD simulation. Unsteady problems with high temporal scale separation can also be simulated, but general transient problems are beyond the capabilities of the current technique.

LanguageEnglish
Pages541-557
Number of pages17
JournalMicrofluidics and Nanofluidics
Volume15
Issue number4
DOIs
Publication statusPublished - 1 Oct 2013

Fingerprint

Nanofluidics
incompressible flow
Incompressible flow
Molecular dynamics
molecular dynamics
continuums
Conservation
conservation
internal flow
simulation
iteration
Momentum
momentum
Computer simulation
conservation laws
high aspect ratio
low speed
Aspect ratio
engineering
Fluxes

Keywords

  • coupled solvers
  • hybrid methods
  • microfluidics
  • molecular dynamics
  • multiscale simulations
  • nanofluidics
  • scale separation

Cite this

Borg, Matthew K. ; Lockerby, Duncan A. ; Reese, Jason M. / A hybrid molecular-continuum simulation method for incompressible flows in micro/nanofluidic networks. In: Microfluidics and Nanofluidics. 2013 ; Vol. 15, No. 4. pp. 541-557.
@article{6565c15195f24660a5dc0aa50d89a5d8,
title = "A hybrid molecular-continuum simulation method for incompressible flows in micro/nanofluidic networks",
abstract = "We present a hybrid molecular-continuum simulation method for modelling nano- and micro-flows in network-type systems. In these types of problem, a full molecular dynamics (MD) description of the macroscopic flow behaviour would be computationally intractable, or at least too expensive to be practical for engineering design purposes. Systems that exhibit multiscale traits, such as this, can instead be solved using a hybrid approach that distinguishes the problem into macroscopic and microscopic dynamics, modelled by their respective solvers. The technique presented in this study is an extension and addition to a hybrid method developed by Borg et al. (J Comput Phys 233:400-413, 2013) for high-aspect-ratio channel geometries, known as the internal-flow multiscale method (IMM). Computational savings are obtained by replacing long channels in the network, which are highly scale-separated, by much smaller, but representative, MD simulations, without a substantial loss of accuracy. On the other hand, junction components do not exhibit this length-scale separation, and so must be simulated in their entirety using MD. The current technique combines all network elements (junctions and channels) together in a coupled simulation using continuum conservation laws. For the case of steady, isothermal, incompressible, low-speed flows, we use the conservation of mass and momentum flux equations to derive a set of molecular-continuum constraints. An algorithm is presented here that computes at each iteration the new constraints on the pressure differences to be applied over individual MD micro-elements (channels and junctions), successively moving closer to macroscopic mass and momentum conservation. We show that hybrid simulations of some example network cases converge quickly, in only a few iterations, and compare very well to the corresponding full MD results, which are taken as the most accurate solutions. Major computational savings can be afforded by the IMM-type approximation in the channel components, but for steady-state solutions, even greater savings are possible. This is because the micro-elements are coupled to a steady-state continuum conservation expression, which greatly speeds up the relaxation of individual micro-components to steady conditions as compared to that of a full MD simulation. Unsteady problems with high temporal scale separation can also be simulated, but general transient problems are beyond the capabilities of the current technique.",
keywords = "coupled solvers, hybrid methods, microfluidics, molecular dynamics, multiscale simulations, nanofluidics, scale separation",
author = "Borg, {Matthew K.} and Lockerby, {Duncan A.} and Reese, {Jason M.}",
year = "2013",
month = "10",
day = "1",
doi = "10.1007/s10404-013-1168-y",
language = "English",
volume = "15",
pages = "541--557",
journal = "Microfluidics and Nanofluidics",
issn = "1613-4982",
publisher = "Springer-Verlag",
number = "4",

}

A hybrid molecular-continuum simulation method for incompressible flows in micro/nanofluidic networks. / Borg, Matthew K.; Lockerby, Duncan A.; Reese, Jason M.

In: Microfluidics and Nanofluidics, Vol. 15, No. 4, 01.10.2013, p. 541-557.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A hybrid molecular-continuum simulation method for incompressible flows in micro/nanofluidic networks

AU - Borg, Matthew K.

AU - Lockerby, Duncan A.

AU - Reese, Jason M.

PY - 2013/10/1

Y1 - 2013/10/1

N2 - We present a hybrid molecular-continuum simulation method for modelling nano- and micro-flows in network-type systems. In these types of problem, a full molecular dynamics (MD) description of the macroscopic flow behaviour would be computationally intractable, or at least too expensive to be practical for engineering design purposes. Systems that exhibit multiscale traits, such as this, can instead be solved using a hybrid approach that distinguishes the problem into macroscopic and microscopic dynamics, modelled by their respective solvers. The technique presented in this study is an extension and addition to a hybrid method developed by Borg et al. (J Comput Phys 233:400-413, 2013) for high-aspect-ratio channel geometries, known as the internal-flow multiscale method (IMM). Computational savings are obtained by replacing long channels in the network, which are highly scale-separated, by much smaller, but representative, MD simulations, without a substantial loss of accuracy. On the other hand, junction components do not exhibit this length-scale separation, and so must be simulated in their entirety using MD. The current technique combines all network elements (junctions and channels) together in a coupled simulation using continuum conservation laws. For the case of steady, isothermal, incompressible, low-speed flows, we use the conservation of mass and momentum flux equations to derive a set of molecular-continuum constraints. An algorithm is presented here that computes at each iteration the new constraints on the pressure differences to be applied over individual MD micro-elements (channels and junctions), successively moving closer to macroscopic mass and momentum conservation. We show that hybrid simulations of some example network cases converge quickly, in only a few iterations, and compare very well to the corresponding full MD results, which are taken as the most accurate solutions. Major computational savings can be afforded by the IMM-type approximation in the channel components, but for steady-state solutions, even greater savings are possible. This is because the micro-elements are coupled to a steady-state continuum conservation expression, which greatly speeds up the relaxation of individual micro-components to steady conditions as compared to that of a full MD simulation. Unsteady problems with high temporal scale separation can also be simulated, but general transient problems are beyond the capabilities of the current technique.

AB - We present a hybrid molecular-continuum simulation method for modelling nano- and micro-flows in network-type systems. In these types of problem, a full molecular dynamics (MD) description of the macroscopic flow behaviour would be computationally intractable, or at least too expensive to be practical for engineering design purposes. Systems that exhibit multiscale traits, such as this, can instead be solved using a hybrid approach that distinguishes the problem into macroscopic and microscopic dynamics, modelled by their respective solvers. The technique presented in this study is an extension and addition to a hybrid method developed by Borg et al. (J Comput Phys 233:400-413, 2013) for high-aspect-ratio channel geometries, known as the internal-flow multiscale method (IMM). Computational savings are obtained by replacing long channels in the network, which are highly scale-separated, by much smaller, but representative, MD simulations, without a substantial loss of accuracy. On the other hand, junction components do not exhibit this length-scale separation, and so must be simulated in their entirety using MD. The current technique combines all network elements (junctions and channels) together in a coupled simulation using continuum conservation laws. For the case of steady, isothermal, incompressible, low-speed flows, we use the conservation of mass and momentum flux equations to derive a set of molecular-continuum constraints. An algorithm is presented here that computes at each iteration the new constraints on the pressure differences to be applied over individual MD micro-elements (channels and junctions), successively moving closer to macroscopic mass and momentum conservation. We show that hybrid simulations of some example network cases converge quickly, in only a few iterations, and compare very well to the corresponding full MD results, which are taken as the most accurate solutions. Major computational savings can be afforded by the IMM-type approximation in the channel components, but for steady-state solutions, even greater savings are possible. This is because the micro-elements are coupled to a steady-state continuum conservation expression, which greatly speeds up the relaxation of individual micro-components to steady conditions as compared to that of a full MD simulation. Unsteady problems with high temporal scale separation can also be simulated, but general transient problems are beyond the capabilities of the current technique.

KW - coupled solvers

KW - hybrid methods

KW - microfluidics

KW - molecular dynamics

KW - multiscale simulations

KW - nanofluidics

KW - scale separation

UR - http://www.scopus.com/inward/record.url?scp=84885322750&partnerID=8YFLogxK

U2 - 10.1007/s10404-013-1168-y

DO - 10.1007/s10404-013-1168-y

M3 - Article

VL - 15

SP - 541

EP - 557

JO - Microfluidics and Nanofluidics

T2 - Microfluidics and Nanofluidics

JF - Microfluidics and Nanofluidics

SN - 1613-4982

IS - 4

ER -