Scientific Machine Learning for Modeling and Simulating Complex Fluids

The formulation of rheological constitutive equations -- models that relate internal stresses and deformations in complex fluids -- is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine learning constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data, and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these 'digital fluid twins' to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation -- a task that is not achievable using any previously developed data-driven rheological equation of state.Comment: 13 pages, 4 figure

The Medium Amplitude Response of Nonlinear Maxwell-Oldroyd Type Models in Simple Shear

A general framework for Maxwell-Oldroyd type differential constitutive models is examined, in which an unspecified nonlinear function of the stress and rate-of-deformation tensors is incorporated into the well-known corotational version of the Jeffreys model discussed by Oldroyd. For medium amplitude simple shear deformations, the recently developed mathematical framework of medium amplitude parallel superposition (MAPS) rheology reveals that this generalized nonlinear Maxwell model can produce only a limited number of distinct signatures, which combine linearly in a well-posed basis expansion for the third order complex viscosity. This basis expansion represents a library of MAPS signatures for distinct constitutive models that are contained within the generalized nonlinear Maxwell model. We describe a framework for quantitative model identification using this basis expansion, and discuss its limitations in distinguishing distinct nonlinear features of the underlying constitutive models from medium amplitude shear stress data. The leading order contributions to the normal stress differences are also considered, revealing that only the second normal stress difference provides distinct information about the weakly nonlinear response space of the model. After briefly considering the conditions for time-strain separability within the generalized nonlinear Maxwell model, we apply the basis expansion of the third order complex viscosity to derive the medium amplitude signatures of the model in specific shear deformation protocols. Finally, we use these signatures for estimation of model parameters from rheological data obtained by these different deformation protocols, revealing that three-tone oscillatory shear deformations produce data that is readily able to distinguish all features of the medium amplitude, simple shear response space of this generalized class of constitutive models.Comment: 26 pages, 11 figure

Maximum Likelihood Estimation for Single Particle, Passive Microrheology Data with Drift

Volume limitations and low yield thresholds of biological fluids have led to widespread use of passive microparticle rheology. The mean-squared-displacement (MSD) statistics of bead position time series (bead paths) are either applied directly to determine the creep compliance [Xu et al (1998)] or transformed to determine dynamic storage and loss moduli [Mason & Weitz (1995)]. A prevalent hurdle arises when there is a non-diffusive experimental drift in the data. Commensurate with the magnitude of drift relative to diffusive mobility, quantified by a P\'eclet number, the MSD statistics are distorted, and thus the path data must be "corrected" for drift. The standard approach is to estimate and subtract the drift from particle paths, and then calculate MSD statistics. We present an alternative, parametric approach using maximum likelihood estimation that simultaneously fits drift and diffusive model parameters from the path data; the MSD statistics (and consequently the compliance and dynamic moduli) then follow directly from the best-fit model. We illustrate and compare both methods on simulated path data over a range of P\'eclet numbers, where exact answers are known. We choose fractional Brownian motion as the numerical model because it affords tunable, sub-diffusive MSD statistics consistent with typical 30 second long, experimental observations of microbeads in several biological fluids. Finally, we apply and compare both methods on data from human bronchial epithelial cell culture mucus.Comment: 29 pages, 12 figure

Pioneering Tree Improvement in Oklahoma

The pioneering tree improvement work in Oklahoma started in 1965 when Clayton Posey moved from Auburn University to Oklahoma State University. Clayton was hired by Glen Durrell (Department Head) to fill a teaching/research position in the Department of Forestry. As a native Oklahoman, Clayton recognized the need to start some long-term studies with the economically important timber species in the state. Fortunately he had access to McIntire-Stennis funds which he used to initiate studies with loblolly pine (Pinus taeda) shortleaf pine (Pinus echinata) and eastern cottonwood (Populus deltoides). Tree selection started in 1966 and concurrently the Kiamichi Field Station was transferred to the Forestry Department from Horticulture. In typical Oklahoma fashion a strong spirit of cooperation emerged with Dierks Lumber Company (soon to be acquired by Weyerhaeuser), Herron Lumber Company, Oklahoma Forestry Division, and the Tiak District of the Ouachita National Forest all assisting with the program. The cooperative spirit was formalized in 1980 when the Oklahoma Forestry Division officially joined the Western Gulf Forest Tree Improvement Program.Papers and abstracts from the 27th Southern Forest Tree Improvement Conference held at Oklahoma State University in Stillwater, Oklahoma on June 24-27, 2003