Exact Analytical Solution for Large‐Amplitude Oscillatory Shear Flow

The Oldroyd 8-constant framework for continuum constitutive theory contains a rich diversity of popular special cases for polymeric liquids. In this paper, we focus on the normal stress difference responses to large-amplitude oscillatory shear (LAOS) flow. The nonlinearity of the polymeric liquids, triggered by LAOS, causes these responses at even multiples of the test frequency. We call responses at a frequency higher than twice the test frequency higher harmonics. The normal stress difference responses for the Oldroyd 8-constant framework has recently yielded to the exact analytical solution. However, in its closed form, Bessel functions appear 24 times, each within summations to infinity. In this paper, to simplify the exact solution, we expand it in a Taylor series. We truncate the series after its 17th power of the shear rate amplitude. Our main result reduces to the well-known expression for the special cases of the corotational Jeffreys and corotational Maxwell fluids. Whereas these special cases yielded to the Goddard integral expansion (GIE), the more general Oldroyd 8-constant framework does not. We use Ewoldt grids to show our main result to be highly accurate for the corotational Jeffreys and corotational Maxwell fluids. For these two special cases, our solutions agree closely with the exact solutions as long as Wi/De<3310. We compare our main result, for the special case of the Johnson–Segalman fluid, with measurements on dissolved polyisobutylene in the isobutylene oligomer. For this, we use the Spriggs relations to generalize our main result to multimode, which then agrees closely with the measurements.

Experiments in which an oscillatory shear flow is superimposed on a steady-state circular shear flow between a cone and a plate were performed on non-Newtonian solutions by means of aWeissenberg Rheogoniometer. The steady-state shear stress and in a first approximation also the normal stress difference arising from the steady shear flow appear not to be influenced by the superimposed oscillatory flow. On the other hand, the dynamic moduli as obtained from the oscillatory parts of shear stress and shear flow are highly dependent on the superimposed steady rate of shear. The absolute value of the complex shear modulus decreases and the phase difference between oscillatory shear stress and shear flow increases in all cases and for all frequencies if the superimposed shear rate is increased. Consequently, this phase difference can become equal to and even larger than π/2. Between the angular frequency ω0 at which the phase difference is π/2 and the steady shear rateq the relation ω0= 1/2,q was experimentally found to exist in most cases. These dynamic results cannot be described by the current theories of viscoelasticity. The large and fast deformations imposed on the material should explicitly be taken into account.

Corrigenda: “Large-amplitude oscillatory shear flow from the corotational Maxwell model” [Journal of Non-Newtonian Fluid Mechanics 166 (2011) 1081–1099

We present a numerical analysis of the rheology of a suspension of red blood cells (RBCs) for different volume fractions in a wall-bounded, effectively inertialess, small amplitude oscillatory shear (SAOS) flow for a wide range of applied frequencies. The RBCs are modeled as biconcave capsules, whose membrane is an isotropic and hyperelastic material following the Skalak constitutive law. The frequency-dependent viscoelasticity in the bulk suspension is quantified by the complex viscosity, defined by the amplitude of the particle shear stress and the phase difference between the stress and shear. SAOS flow basically impedes the deformation of individual RBCs as well as the magnitude of fluid-membrane interactions, resulting in a lower specific viscosity and first and second normal stress differences than in steady shear flow. Although it is known that the RBC deformation alone is sufficient to give rise to shear-thinning, our results show that the complex viscosity weakly depends on the frequency-modulated deformations or orientations of individual RBCs but rather depends on combinations of the frequency-dependent amplitude and phase difference. The effect of the viscosity ratio between the cytoplasm and plasma and of the capillary number is also assessed.

Analytical solutions for either the shear stress or the normal stress differences in large-amplitude oscillatory shear flow, both for continuum or molecular models, often take the form of the first few terms of a power series in the shear rate amplitude. Here, we explore improving the accuracy of these truncated series by replacing them with ratios of polynomials. Specifically, we examine replacing the truncated series solution for the corotational Maxwell model with its Pade approximants for the shear stress response and for the normal stress differences. We find these Pade approximants to agree closely with the corresponding exact solution, and we learn that with the right approximants, one can nearly eliminate the inaccuracies of the truncated expansions.

We investigated the nonlinear and nonsymmetric responses of viscoelastic fluids under large amplitude oscillatory squeeze (LAOSQ) flow. The nonlinear and nonsymmetric stress response is a unique feature of oscillatory squeeze (OSQ) flow under larger deformation, but has rarely been investigated. The goal of this study is to systematically characterize the responses of viscoelastic fluids at larger deformation under oscillatory squeeze flow and to provide a platform for the analysis of nonsymmetric stress signals. We report a framework for the analysis of nonlinear-and-nonsymmetric stress signals at larger strain amplitude under oscillatory squeeze flow. The storage and loss modulus showed strain thinning behavior as the strain amplitude increases in both oscillatory shear and oscillatory squeeze flow. However the normal stress under LAOSQ was found to be nonsymmetric in both magnitude and shape unlike the shear stress under oscillatory shear flow. In addition the energy dissipation was found to be larger in extension than in compression. This study is expected to provide a platform to understand the nonlinear and nonsymmetric characteristics of complex fluids under LAOSQ flow.

The rheological behavior of blends of poly(methyl methacrylate) (PMMA) and poly(acrylonitrile-stat-styrene)-graft-polybutadiene (ABS) was investigated using a cone-and-plate rheometer. The rheological properties measured were shear stress (σ12), viscosity (η), and first normal stress difference (N1) as functions of shear rate (\documentclass{article}\pagestyle{empty}\begin{document}$ \dot \gamma$ \end{document}) in steady shearing flow, and storage modulus (G′) and loss modulus (G″) as functions of frequency (ω) in oscillatory shearing flow. It has been found that the rheological behavior of blends of ABS and PMMA was very similar to that of blends of poly(styrene-stat-acrylonitrile) (SAN) and PMMA, in that N1 in logarithmic plots of N1 versus σ12, and G′ in logarithmic plots of G′ versus G″, vary regularly with blend composition. This has led us to conclude that the rubber particles that are grafted on an SAN resinous matrix in ABS resin plays only a minor role in influencing the compatibility of ABS/PMMA blends, and that the SAN chains attached to the surface of rubber particles, and the SAN matrix phase, play a major role in compatibilizing ABS resin with PMMA.

We examine the simplest relevant molecular model for large-amplitude shear (LAOS) flow of a polymeric liquid: the suspension of rigid dumbbells in a Newtonian solvent. We find explicit analytical expressions for the shear rate amplitude and frequency dependences of the first and third harmonics of the alternating shear stress response. We include a detailed comparison of these predictions with the corresponding results for the simplest relevant continuum model: the corotational Maxwell model. We find that the responses of both models are qualitatively similar. The rigid dumbbell model relies entirely on the dumbbell orientation to explain the viscoelastic response of the polymeric liquid, including the higher harmonics in large-amplitude oscillatory shear flow. Our analysis employs the general method of Bird and Armstrong ["Time-dependent flows of dilute solutions of rodlike macromolecules," J. Chem. Phys. 56, 3680 (1972)] for analyzing the behavior of the rigid dumbbell model in any unsteady shear flow. We derive the first three terms of the deviation of the orientational distribution function from the equilibrium state. Then, after getting the "paren functions," we use these for evaluating the shear stress for LAOS flow. We find the shapes of the shear stress versus shear rate loops predicted to be reasonable.

We demonstrate that erythrocyte deformations, specifically of a type as occur in splenic flow (Zhu et al., 2017), and of the type that promote vesiculation can be caused by simple, yet tailored, oscillatory shear flow. We show that such oscillatory shear flow provides an ideal environment to explore a wide variety of metabolic and biochemical effects that promote erythrocyte vesiculation. Deformation details, typical of splenic flow, such as in-folding and implications for membrane/skeleton interaction are demonstrated and quantitatively analyzed. We introduce a theoretical, essentially analytical, vesiculation model that directly couples to our more complex numerical, multilevel, model that clearly delineates various fundamental elements, i.e., sub-processes, that are involved and mediate the vesiculation process. This analytical model highlights particulary important vesiculation precursors such as areas of membrane/skeleton disruptions that trigger the vesiculation process. We demonstrate, using flow cytometry, that the deformations we experimentally induce on cells, and numerically simulate, do not induce lethal forms of cell damage but do induce vesiculation as theoretically forecasted. This, we demonstrate, provides a direct link to cell membrane/skeletal damage such as is associated with metabolic and aging damage. An additional noteworthy feature of this approach is the avoidance of artificial devices, e.g., micro-fluidic chambers, in which deformations and their time scales are often unrepresentative of physiological processes such as splenic flow.

Results are reported for the dynamic moduli,G′ andG″, measured mechanically, and the dynamic third normal stress difference, measured optically, of a series bidisperse linear polymer melts under oscillatory shear. Nearly monodisperse hydrogenated polyisoprenes of molecular weights 53000 and 370000 were used to prepare blends with a volume fraction of long polymer,Φ L, of 0.10, 0.20, 0.30, 0.50, and 0.75. The results demonstrate the applicability of birefringence measurements to solve the longstanding problem of measuring the third normal stress difference in oscillatory flow. The relationship between the third normal stress difference and the shear stress observed for these entangled polymer melts is in agreement with a widely predicted constitutive relationship: the relationship between the first normal stress difference and the shear stress is that of a simple fluid, and the second normal stress difference is proportional to the first. These results demonstrate the potential use of 1,3-birefringence to measure the third normal stress difference in oscillatory flow. Further, the general constitutive equation supported by the present results may be used to determine the dynamic moduli from the measured third normal stress difference in small amplitude oscillatory shear. Directions for future research, including the use of birefringence measurements to determineN 2/N 1 in oscillatory shear, are described.

Exact static solutions to piezoelectric smart beams including peel stresses: I: Theoretical formulation

This paper has introduced a pseudo-potential in bond-fluctuation model to simulate oscillatory shear flow of multiple self-avoiding chains in three dimensions following our previous work under simple shear flow. The oscillatory flow field was reasonably reproduced by lattice Monte Carlo simulation using this pseudo-potential neglecting hydrodynamic interaction. By sampling the configuration distribution functions, the macroscopic viscoelasticity of semi-concentrated polymer solution was determined. Both Newtonian and non-Newtonian regimes were studied. The complex modulus and dynamic viscosity exhibit a reasonable power relation with oscillatory frequency, which is consistent with present theories and experiments. Consequently, lattice Monte Carlo simulation has been extended to model free-draining self-avoiding multi chains subject to oscillatory shear flow and to investigate associated viscoelasticity on the molecular level.

The Al project—an overview

Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.

Dynamics of rigid dumbbells in confined geometries Part II. Time-dependent shear flow