Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows

Journal Article · Fri Jun 01 04:00:00 EDT 2018 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2018.02.039· OSTI ID:1565655

Zhao, Lifei; Li, Zhen; Caswell, Bruce; Ouyang, Jie; Karniadakis, George Em

We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number RE = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.

• Research Organization: Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Argonne National Lab. (ANL), Argonne, IL (United States); UT-Battelle LLC/ORNL, Oak Ridge, TN (Unted States)
• Sponsoring Organization: USDOE Office of Science
• DOE Contract Number: AC02-06CH11357; AC05-00OR22725
• OSTI ID: 1565655
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: C Vol. 363; ISSN 0021-9991
• Publisher: Elsevier
• Country of Publication: United States
• Language: English

Related Subjects

Computer Science
Physics