### 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.

Language | English |
---|---|

Pages | 541-557 |

Number of pages | 17 |

Journal | Microfluidics and Nanofluidics |

Volume | 15 |

Issue number | 4 |

DOIs | |

Publication status | Published - 1 Oct 2013 |

### Fingerprint

### Keywords

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

### Cite this

*Microfluidics and Nanofluidics*,

*15*(4), 541-557. https://doi.org/10.1007/s10404-013-1168-y

}

*Microfluidics and Nanofluidics*, vol. 15, no. 4, pp. 541-557. https://doi.org/10.1007/s10404-013-1168-y

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

Research output: Contribution to journal › Article

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 -