Multiscale Interface and Absorbing Boundary Condition for Non-Equilibrium Molecular Dynamics Simulation: Semi-Analytical Formulation and Machine Learning

Autor: Yu-Chuan Hsu, 許育銓
Rok vydání: 2018
Druh dokumentu: 學位論文 ; thesis
Popis: 106
Molecular dynamics (MD) simulation has been widely applied in studying materials behavior in the past decades. However, owning to its great demands on the modeling fidelity, MD is often limited by the length scale it can handle. Although it is much desired to develop multiscale modeling to overcome this limitation, the development of dynamic multiscale modeling is limited by the spurious wave-reflection problems on the interface between atoms and continuum. This study has developed a multiscale semi-analytical formulation for arbitrary lattices in MD. This method decomposes MD system into a coupled-domain of real and virtual domains, dissembles the contribution from the virtual domain in frequency, and derives the time-history kernel function (THKF) accounting for the interaction of virtual domain implicitly. The virtual atoms on the interface can further be controlled by this numerical method after pairing of degree-of-freedom of THKF and atoms on the interface. By doing so, the dynamics simulation without considering virtual atoms will be equivalent to the original full MD simulation. In a broad sense, multiscale interface (MI) from this study is a two-way implicit interface across real domains, which is not limited to the number of real domains and the way of coupling. However, in the special case of single real domain, multiscale interface reduces to the classical absorbing boundary condition (ABC). There is some key difference between these two applications. ABC totally ignores the information of the virtual domain, thus the time-cutoff of convolution can reduce the computational cost and maintain the bound of the non-reflecting system. On the other hand, due to the transition of the information between different real domains, the simulation of MI will result in spurious wave reflection if THKF are interfered by the time-cutoff. In addition, THKF analyzed from the MI has strong physical meanings of real space that contains the size effect, which can be further developed to resolve length scale limitation in MD. In this thesis, we also explore the possibility of applying machine learning (ML) in ABC based on the function form in THKF. The preliminary verification has been completed and ML for ABC can reduce the computational efforts substantially. This poses a good direction to develop ML for MI in the future.
Databáze: Networked Digital Library of Theses & Dissertations