skip to main content


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 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.  more » « less
Award ID(s):
1727870
NSF-PAR ID:
10332579
Author(s) / Creator(s):
;
Date Published:
Journal Name:
International Journal of Computational Methods
ISSN:
0219-8762
Page Range / eLocation ID:
2143003
Format(s):
Medium: X
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 creep 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. 
    more » « 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. 
    more » « less
  3. Three-dimensional printing (3DP) of functional materials is increasingly important for advanced applications requiring objects with complex or custom geometries or prints with gradients or zones with different properties. A common 3DP technique is direct ink writing (DIW), in which printable inks are comprised of a fluid matrix filled with solid particles, the latter of which can serve a dual purpose of rheology modifiers to enable extrusion and functional fillers for performance-related properties. Although the relationship between filler loading and viscosity has been described for many polymeric systems, a thorough description of the rheological properties of three-dimensional (3D) printable composites is needed to expedite the creation of new materials. In this manuscript, the relationship between filler loading and printability is studied using model paraffin/photopolymer composite inks containing between 0 and 73 vol. % paraffin microbeads. The liquid photopolymer resin is a Newtonian fluid, and incorporating paraffin microbeads increases the ink viscosity and imparts shear-thinning behavior, viscoelasticity, and thixotropy, as established by parallel plate rheometry experiments. Using Einstein and Batchelor's work on colloidal suspension rheology, models were developed to describe the thixotropic behavior of inks, having good agreement with experimental results. Each of these properties contributes to the printability of highly filled ([Formula: see text]43 vol. % paraffin) paraffin/photopolymer composite inks. Through this work, the ability to quantify the ideal rheological properties of a DIW ink and to selectively control and predict its rheological performance will facilitate the development of 3D printed materials with tunable functionalities, thus, advancing 3DP technology beyond current capabilities. 
    more » « less
  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 the 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. 
    more » « less
  5. 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.

     
    more » « less