Gil Shohet , Jacob Price , Jeffrey Haack , Mathieu Marciante , Michael S. Murillo
{"title":"Heterogeneous multiscale method for high energy-density matter: Connecting kinetic theory and molecular dynamics","authors":"Gil Shohet , Jacob Price , Jeffrey Haack , Mathieu Marciante , Michael S. Murillo","doi":"10.1016/j.jcpx.2020.100070","DOIUrl":null,"url":null,"abstract":"<div><p>We have developed a concurrent heterogeneous multiscale method (HMM) framework with a microscale molecular dynamics (MD) model and a macroscale kinetic Vlasov-BGK model. The kinetic model is formulated such that BGK collision times are the closure data obtained from MD. Using the H-theorem, we develop the mathematical link between the MD and the kinetic model. We examine three relaxation processes, energy, momentum, and bump-on-tail, using full microscale MD simulations as a reference solution. We find that solutions computed with the HMM framework offer a significant computational reduction (<span><math><mn>14</mn><mo>×</mo><mo>−</mo><mn>100</mn><mo>×</mo></math></span>) compared with computing a full MD solution, with significant improvements in accuracy compared with a kinetic model using analytical collision times.</p></div>","PeriodicalId":37045,"journal":{"name":"Journal of Computational Physics: X","volume":"8 ","pages":"Article 100070"},"PeriodicalIF":0.0000,"publicationDate":"2020-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1016/j.jcpx.2020.100070","citationCount":"3","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics: X","FirstCategoryId":"1085","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S2590055220300226","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2020/8/17 0:00:00","PubModel":"Epub","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 3
Abstract
We have developed a concurrent heterogeneous multiscale method (HMM) framework with a microscale molecular dynamics (MD) model and a macroscale kinetic Vlasov-BGK model. The kinetic model is formulated such that BGK collision times are the closure data obtained from MD. Using the H-theorem, we develop the mathematical link between the MD and the kinetic model. We examine three relaxation processes, energy, momentum, and bump-on-tail, using full microscale MD simulations as a reference solution. We find that solutions computed with the HMM framework offer a significant computational reduction () compared with computing a full MD solution, with significant improvements in accuracy compared with a kinetic model using analytical collision times.