Accelerating multiscale modelling of fluids with on-the-fly Gaussian process regression

• Published in: Microfluidics and nanofluidics Vol. 22; no. 12; pp. 139 - 12
• Main Authors: Stephenson, David, Kermode, James R., Lockerby, Duncan A.
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.12.2018; Springer Nature B.V
• Subjects: Accuracy; Analytical Chemistry; Biomedical Engineering and Bioengineering; Computational fluid dynamics; Computer applications; Computer simulation; Dynamics; Engineering; Engineering Fluid Dynamics; Fluids; Gaussian process; High aspect ratio; Modelling; Molecular dynamics; Nanotechnology and Microengineering; Noise threshold; Regression models; Research Paper; Simulation; Thermal noise; Thresholds; Uncertainty; Unsteady flow; Molecular dynamics; Hybrid methods; Multiscale modelling; Machine learning; Micro/nanofluidics
• ISSN: 1613-4982; 1613-4990
• DOI: 10.1007/s10404-018-2164-z

We present a scheme for accelerating hybrid continuum-atomistic models in multiscale fluidic systems by using Gaussian process regression as a surrogate model for computationally expensive molecular dynamics simulations. Using Gaussian process regression, we are able to accurately predict atomic-scale information purely by consideration of the macroscopic continuum-model inputs and outputs and judge on the fly whether the uncertainty of our prediction is at an acceptable level, else a new molecular simulation is performed to continually augment the database, which is never required to be complete. This provides a substantial improvement over the current generation of hybrid methods, which often require many similar atomistic simulations to be performed, discarding information after it is used once. We apply our hybrid scheme to nano-confined unsteady flow through a high-aspect-ratio converging–diverging channel, and make comparisons between the new scheme and full molecular dynamics simulations for a range of uncertainty thresholds and initial databases. For low thresholds, our hybrid solution is highly accurate—around that of thermal noise. As the uncertainty threshold is raised, the accuracy of our scheme decreases and the computational speed-up increases (relative to a full molecular simulation), enabling the compromise between accuracy and efficiency to be tuned. The speed-up of our hybrid solution ranges from an order of magnitude, with no initial database, to cases where an extensive initial database ensures no new MD simulations are required.