 Award ID(s):
 2011401
 NSFPAR ID:
 10411333
 Date Published:
 Journal Name:
 Biomicrofluidics
 Volume:
 16
 Issue:
 5
 ISSN:
 19321058
 Page Range / eLocation ID:
 054109
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Nonlinear rheology of entangled wormlike micellar solutions predicted by a micelleslipspring modelWe examine linear and nonlinear shear and extensional rheological properties using a “micelleslipspring model” [T. Sato et al., J. Rheol. 64, 1045–1061 (2020)] that incorporates breakage and rejoining events into the slipspring model originally developed by Likhtman [Macromolecules 38, 6128–6139 (2005)] for unbreakable polymers. We here employ the Fraenkel potential for main chain springs and slipsprings to address the effect of finite extensibility. Moreover, to improve extensional properties under a strong extensional flow, stressinduced micelle breakage (SIMB) is incorporated into the micelleslipspring model. Thus, this model is the first model that includes the entanglement constraint, Rouse modes, finite extensibility, breakage and rejoining events, and stressinduced 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 startup 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 micelleslipspring 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 micelleslipspring model with SIMB show thickening followed by thinning, which is in qualitative agreement with experimental trends. Additionally, the extensional rheological properties of the micelleslipspring 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 slipspring model in shear and extensional flows and so provide an accurate, efficient model of micellar solution rheology.more » « less

Dissipative particle dynamics (DPD) simulations are performed on coarsegrained replicas of linear, monodisperse entangled polyethylene melts [Formula: see text] and [Formula: see text] undergoing both steadystate and transient planar elongational flow (PEF). The fidelity of the DPD simulations is verified by direct comparison of flow topological and rheological properties of a 334particle chain liquid against the unitedatom [Formula: see text] liquid, simulated using nonequilibrium molecular dynamics (NEMD). These DPD simulations demonstrate that a flowinduced coilstretch transition (CST) and its associated hysteresis caused by configurational microphase separation, as observed in previous NEMD simulations of PEF, can be replicated using a more computationally efficient coarsegrained system. Results indicate that the breadth of the CST hysteresis loop is enlarged for the longer molecule liquid relative to the shorter one. Furthermore, relaxation simulations reveal that reducing the applied flow Deborah number ([Formula: see text]) from a high value corresponding to a homogeneous phase of highly stretched molecules to a [Formula: see text] within the biphasic region results in a twostage relaxation process. There is a fast initial stratification of the kinetically trapped highly stretched chains into regions of highly extended and less extended chains, which displays similar behavior to a system undergoing a spinodal decomposition caused by spatial configurational free energy fluctuations. After a short induction period of apparently random duration, the less extended chain regions experience a stochastic nucleation event that induces configurational relaxation to domains composed of coiled molecules over a much longer time scale, leaving the more highly extended chains in surrounding sheetlike domains. The time scales of these two relaxation processes are of the same order of magnitude as the Rouse and disengagement times of the equilibrium liquids.

Abstract Extensional flow properties of polymer solutions in volatile solvents govern many industriallyrelevant coating processes, but existing instrumentation lacks the environment necessary to control evaporation. To mitigate evaporation during drippingontosubstrate (DoS) extensional rheology measurements, we developed a chamber to enclose the sample in an environment saturated with solvent vapor. We validated the evaporationcontrolled 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
for PEO in volatile solvents like dichloromethane and chloroform. PEO/chloroform solutions displayed an over 20fold increase in$$\lambda _{E}$$ ${\lambda}_{E}$ due to the formation of an evaporationinduced 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 semidilute PEO/chloroform solutions were measured with environmental control, where$$\lambda _{E}$$ ${\lambda}_{E}$ scaled with concentration by the exponent$$\lambda _{E}$$ ${\lambda}_{E}$ . These measurements validate the evaporationcontrolled DoS environment, and confirm that chloroform is a good solvent for PEO, with a Flory exponent of$$m=0.62$$ $m=0.62$ . 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.$$\nu =0.54$$ $\nu =0.54$ 
Monodisperse suspensions of Brownian colloidal spheres crystallize at high densities, and ordering under shear has been observed at densities below the crystallization threshold. We perform largescale simulations of a model suspension containing over [Formula: see text] particles to quantitatively study the ordering under shear and to investigate its link to the rheological properties of the suspension. We find that at high rates, for [Formula: see text], the shear flow induces an ordering transition that significantly decreases the measured viscosity. This ordering is analyzed in terms of the development of layering and planar order, and we determine that particles are packed into hexagonal crystal layers (with numerous defects) that slide past each other. By computing local [Formula: see text] and [Formula: see text] order parameters, we determine that the defects correspond to chains of particles in a squarelike lattice. We compute the individual particle contributions to the stress tensor and discover that the largest contributors to the shear stress are primarily located in these lower density, defect regions. The defect structure enables the formation of compressed chains of particles to resist the shear, but these chains are transient and shortlived. The inclusion of a contact friction force allows the stressbearing structures to grow into a systemspanning network, thereby disrupting the order and drastically increasing the suspension viscosity.more » « less

This work experimentally investigates the pressuredriven flow of a pure Newtonian fluid through threedimensional (3D) porous media models. The porous media model consists of square arrays of rods that also could be interpreted as a periodic tandem rod arrangement. We employed a timeresolved threedimensional particle tracking velocimetry (3D ShaketheBox) technique for a range of Reynolds numbers [Formula: see text] to observe flow structures and vortex formation between the rods in porous media structures with different porosities of [Formula: see text] which corresponds to the spacing ratio of [Formula: see text], where L is the distance between the centers of the rods, and D is the diameter of the rods. For all the examined cases, we further analyzed the effect of the Reynolds number and the spacing ratio on the instantaneous and averaged patterns of velocity, vorticity, and the other flow parameters after obtaining the twodimensional velocity fields using the binaveraging method. We observed both symmetrical and asymmetrical patterns of structure and recirculation regions between the rods depending on the Reynolds number and spacing ratio. Increasing the Reynolds number reduced the symmetrical patterns of flow structures with respect to the centerline of the gap region, while the spacing ratio was randomly affecting the symmetry degree. Vortex shedding was considerable for the two examined high Reynolds numbers of Re = 444 and Re = 890 behind the upstream rod as the porosity increased. The backward movement of the reattachment point has been observed by increasing the Reynolds number.more » « less