skip to main content

This content will become publicly available on February 14, 2023

Title: Surrogate Modeling with Gaussian Processes for an Inverse Problem in Polymer Dynamics
When rheological models of polymer blends are used for inverse modeling, they can characterize polymer mixtures from rheological observations. This requires repeated evaluation of potentially expensive rheological models. We explored surrogate models based on Gaussian processes (GP-SM) as a cheaper alternative for describing the rheology of polydisperse binary blends. We used the time-dependent diffusion double reptation (TDD-DR) model as the true model; it takes a 5-dimensional input vector specifying the binary blend as input and yields a function called the relaxation spectrum as output. We used the TDD-DR model to generate training data of different sizes [Formula: see text], via Latin hypercube sampling. The optimal values of the GP-SM hyper-parameters, assuming a separable covariance kernel, were obtained by maximum likelihood estimation. The GP-SM interpolates the training data by design and offers reasonable predictions of relaxation spectra with uncertainty estimates. In general, the accuracy of GP-SMs improves as the size of the training data [Formula: see text] increases, as does the cost for training and prediction. The optimal hyper-parameters were found to be relatively insensitive to [Formula: see text]. Finally, we considered the inverse problem of inferring the structure of the polymer blend from a synthetic dataset generated using the true more » model. Surprisingly, the solution to the inverse problem obtained using GP-SMs and TDD-DR was qualitatively similar. GP-SMs can be several orders of magnitude cheaper than expensive rheological models, which provides a proof-of-concept validation for using GP-SMs for inverse problems in polymer rheology. « less
Award ID(s):
Publication Date:
Journal Name:
International Journal of Computational Methods
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Fungi such as Candida albicans exist in biofilm phenotypes, which present as viscoelastic materials; however, a method to measure linear viscoelastic moduli, yield stress, and yield strain is lacking. Characterization methods for fungal materials have been limited to techniques specific to particular industries. Here, we present a method to measure the shear stress, strain amplitude, and creep of C. albicans BWP17 biofilms. Our method includes features tailored to the analysis of fungi including an in vitro growth protocol attuned to the slow growth rates of C. albicans biofilms and a resultant cultured biofilm that has sufficient integrity to be transferred to the rheometer tooling without disrupting its structure. The method's performance is demonstrated by showing that results are insensitive to gap, evaporative sealant, length of experiment, and specimen radius. Multiscale imaging of the fungal biofilm showed complex entanglement networks at the hundred-micrometer scale. For a wild-type strain cultivated for 14 days, using small-amplitude oscillatory rheology, we found that the elastic (G′) and viscous (G″) moduli were nearly independent of frequency over the range 0.1–10 s −1 , with magnitudes of [Formula: see text] and [Formula: see text], respectively. The yield stress was approximately [Formula: see text]. We modeled the linear creepmore »response of the fungal biofilm and found that C. albicans has a characteristic relaxation time of [Formula: see text] and a viscosity of [Formula: see text]. We applied this method to probe the effects of altered chitin deposition in the C. albicans cell wall. Differences between the biofilm's phenotypic cell shape and rheological properties in mutants with altered chitin synthase activity were resolved. Discovering how genotypic, phenotypic, and environmental factors impact the material properties of these microbial communities can have implications for understanding fungal biofilm growth and aid in the development of remediation strategies.« less
  2. It is computationally expensive to predict reliability using physical models at the design stage if many random input variables exist. This work introduces a dimension reduction technique based on generalized sliced inverse regression (GSIR) to mitigate the curse of dimensionality. The proposed high dimensional reliability method enables active learning to integrate GSIR, Gaussian Process (GP) modeling, and Importance Sampling (IS), resulting in an accurate reliability prediction at a reduced computational cost. The new method consists of three core steps, 1) identification of the importance sampling region, 2) dimension reduction by GSIR to produce a sufficient predictor, and 3) construction of a GP model for the true response with respect to the sufficient predictor in the reduced-dimension space. High accuracy and efficiency are achieved with active learning that is iteratively executed with the above three steps by adding new training points one by one in the region with a high chance of failure.
  3. Abstract

    Extensional flow properties of polymer solutions in volatile solvents govern many industrially-relevant coating processes, but existing instrumentation lacks the environment necessary to control evaporation. To mitigate evaporation during dripping-onto-substrate (DoS) extensional rheology measurements, we developed a chamber to enclose the sample in an environment saturated with solvent vapor. We validated the evaporation-controlled DoS device by measuring a model high molecular weight polyethylene oxide (PEO) in various organic solvents both inside and outside of the chamber. Evaporation substantially increased the extensional relaxation time$$\lambda _{E}$$λEfor PEO in volatile solvents like dichloromethane and chloroform. PEO/chloroform solutions displayed an over 20-fold increase in$$\lambda _{E}$$λEdue to the formation of an evaporation-induced surface film; evaporation studies confirmed surface features and skin formation reminiscent of buckling instabilities commonly observed in drying polymer solutions. Finally, the relaxation times of semi-dilute PEO/chloroform solutions were measured with environmental control, where$$\lambda _{E}$$λEscaled with concentration by the exponent$$m=0.62$$m=0.62. These measurements validate the evaporation-controlled DoS environment, and confirm that chloroform is a good solvent for PEO, with a Flory exponent of$$\nu =0.54$$ν=0.54. Our results are the first to control evaporation during DoS extensional rheology, and provide guidelines establishing when environmental control is necessary to obtain accurate rheological parameters.

  4. We examine linear and nonlinear shear and extensional rheological properties using a “micelle-slip-spring model” [T. Sato et al., J. Rheol. 64, 1045–1061 (2020)] that incorporates breakage and rejoining events into the slip-spring model originally developed by Likhtman [Macromolecules 38, 6128–6139 (2005)] for unbreakable polymers. We here employ the Fraenkel potential for main chain springs and slip-springs to address the effect of finite extensibility. Moreover, to improve extensional properties under a strong extensional flow, stress-induced micelle breakage (SIMB) is incorporated into the micelle-slip-spring model. Thus, this model is the first model that includes the entanglement constraint, Rouse modes, finite extensibility, breakage and rejoining events, and stress-induced micelle breakage. Computational expense currently limits the model to micellar solutions with moderate numbers of entanglements ([Formula: see text]), but for such solutions, nearly quantitative agreement is attained for the start-up of the shearing flow. The model in the extensional flow cannot yet be tested owing to the lack of data for this entanglement level. The transient and steady shear properties predicted by the micelle-slip-spring model for a moderate shear rate region without significant chain stretch are fit well by the Giesekus model but not by the Phan–Thien/Tanner (PTT) model, which is consistent with themore »ability of the Giesekus model to match experimental shear data. The extensional viscosities obtained by the micelle-slip-spring model with SIMB show thickening followed by thinning, which is in qualitative agreement with experimental trends. Additionally, the extensional rheological properties of the micelle-slip-spring model with or without SIMB are poorly predicted by both the Giesekus and the PTT models using a single nonlinear parameter. Thus, future work should seek a constitutive model able to capture the behavior of the slip-spring model in shear and extensional flows and so provide an accurate, efficient model of micellar solution rheology.« less
  5. We consider the problem of estimating the discrete clustering structures under the sub-Gaussian mixture model. Our main results establish a hidden integrality property of a semidefinite programming (SDP) relaxation for this problem: while the optimal solution to the SDP is not integer-valued in general, its estimation error can be upper bounded by that of an idealized integer program. The error of the integer program, and hence that of the SDP, are further shown to decay exponentially in the signal-to-noise ratio. In addition, we show that the SDP relaxation is robust under the semirandom setting in which an adversary can modify the data generated from the mixture model. In particular, we generalize the hidden integrality property to the semirandom model and thereby show that SDP achieves the optimal error bound in this setting. These results together highlight the “global-to-local” mechanism that drives the performance of the SDP relaxation. To the best of our knowledge, our result is the first exponentially decaying error bound for convex relaxations of mixture models. A corollary of our results shows that in certain regimes, the SDP solutions are in fact integral and exact. More generally, our results establish sufficient conditions for the SDP to correctly recovermore »the cluster memberships of [Formula: see text] fraction of the points for any [Formula: see text]. As a special case, we show that under the [Formula: see text]-dimensional stochastic ball model, SDP achieves nontrivial (sometimes exact) recovery when the center separation is as small as [Formula: see text], which improves upon previous exact recovery results that require constant separation.« less