Investigation of shear thinning behavior and microstructures of MWCNT/epoxy and CNF/epoxy suspensions under steady shear conditions

This report investigates the steady-state viscosities of multiwall carbon nanotube (MWCNT)/epoxy and carbon nanofiber (CNF)/epoxy suspensions with varying filler concentrations under different shear rates at various temperatures. In situ observation of filler networks suggests the build-up of shear induced MWCNT and CNF agglomerates at low shear rates, which correlates with the measured shear thinning behavior. The agglomeration process in MWCNT/epoxy suspensions is enhanced at lower shear rates in the case of higher temperatures, whereas, at high shear rates, both nano-fillers show good dispersion. Shear thinning behavior is observed for both types of fillers, and shear thinning exponential parameters are evaluated as a function of filler content. The shear thinning exponent increases in conjunction with increase in filler content, but it is found to saturate at a specific value, independently of filler material. Finally, the micromechanical elasticity-based analogy model is applied to the prediction of steady state shear viscosity of suspensions at higher shear rates with the assumption of complete dispersion and alignment of individual nano-fillers in suspensions. The predicted viscosities and the experimental data at higher shear rates are compared. The results conclude that fairly good agreement can be seen for the cases of CNF/epoxy suspensions at lower temperatures, whereas MWCNT/epoxy suspensions and CNF/epoxy suspensions at higher temperatures show discrepancy between the prediction and the experimental data.     

Association under shear flow in aqueous solutions of pectin

Effects of oscillatory and steady shear flows on intermolecular associations in dilute and semidilute aqueous solutions of pectin in the absence and presence of the hydrogen bond breaking agent urea are reported. A weak oscillatory shear perturbation builds up, depending on polymer concentration, multichain aggregates or networks in the course of time and these association structures are mainly stabilized through hydrogen bonds. The association effect is more pronounced at higher concentrations, and the growth of intermolecular interactions is inhibited by the addition of urea. Steady shear measurements on the pectin-water solutions reveal shear thickening at low shear rates for all the concentrations, except the lowest one, and disruption of intermolecular junctions at high shear rates. In the presence of urea, no shear thickening is detected. The polymer concentration dependence of the viscosity at a low shear rate can be described by a power law η ∼ c^x, with x = 1.9 and 1.4 without and with urea, respectively. When a low constant shear rate is applied to pectin solutions and this monitoring shear rate is interrupted periodically by transitory high shear rates perturbations during a short time, prominent association structures evolve upon return to the monitoring shear rate. This effect is more evident at a lower polymer concentration, and in the presence of urea, the growth of the association complexes is damped. The shear-induced alignment and stretching of polymer chains and the formation of hydrogen-bonded structures are analyzed in the framework of a model, where cooperative zipping of stretched chains play an important role. Viscosity enhancement is found for a semidilute pectin-water solution in the presence of moderate levels of salt addition (NaCl), suggesting that partial screening of electrostatic interactions promotes growth of energetic cross-links.     

Structure and transport properties of polymer grafted nanoparticles

We perform molecular dynamics simulations on a bead-spring model of pure polymer grafted nanoparticles (PGNs) and of a blend of PGNs with a polymer melt to investigate the correlation between PGN design parameters (such as particle core concentration, polymer grafting density, and polymer length) and properties, such as microstructure, particle mobility, and viscous response. Constant strain-rate simulations were carried out to calculate viscosities and a constant-stress ensemble was used to calculate yield stresses. The PGN systems are found to have less structural order, lower viscosity, and faster diffusivity with increasing length of the grafted chains for a given core concentration or grafting density. Decreasing grafting density causes depletion effects associated with the chains leading to close contacts between some particle cores. All systems were found to shear thin, with the pure PGN systems shear thinning more than the blend; also, the pure systems exhibited a clear yielding behavior that was absent in the blend. Regarding the mechanism of shear thinning at the high shear rates examined, it was found that the shear-induced decrease of Brownian stresses and increase in chain alignment, both correlate with the reduction of viscosity in the system with the latter being more dominant. A coupling between Brownian stresses and chain alignment was also observed wherein the non-equilibrium particle distribution itself promotes chain alignment in the direction of shear.     

Rheology of hydrophobically modified polyacrylamide-co-poly(acrylic acid) on addition of surfactant and variation of solution pH

Constant shear and shear dependent viscosity measurements are reported in aqueous systems of co- and terpolymers of acrylamide (AM), N-n-alkylacrylamide (C10, C12, and C14 alkyl groups), and acrylic acid (AA) with added anionic surfactant sodium dodecyl sulfate (SDS). The results are presented as three-dimensional plots of viscosity vs surfactant concentration and pH at constant shear rate or viscosity vs shear rate and surfactant concentration at constant pH. For terpolymers incorporating AA, a strong viscosity maximum is observed at intermediate pH values (pH 4-6) where the AA groups are partially ionized and at SDS concentrations close to the critical micelle concentration. At high pH, all AA incorporating terpolymer solutions with SDS are strongly shear thinning, but at pH 3-4 the systems of terpolymers with SDS are strongly shear thickening at low shear, followed by a shear-thinning region at high shear. These results are explained in terms of surfactant-mediated network formation with polymer coil expansion and hydrogen bonding between partially ionized AA groups as additional factors.     

Shear thinning in dilute polymer solutions

We use bead-spring models for a polymer coupled to a solvent described by multiparticle collision dynamics to investigate shear thinning effects in dilute polymer solutions. First, we consider the polymer motion and configuration in a shear flow. For flexible polymer models we find a sharp increase in the polymer radius of gyration and the fluctuations in the radius of gyration at a Weissenberg number approximately 1. We then consider the polymer viscosity and the effect of solvent quality, excluded volume, hydrodynamic coupling between the beads, and finite extensibility of the polymer bonds. We conclude that the excluded volume effect is the major cause of shear thinning in polymer solutions. Comparing the behavior of semiflexible chains, we find that the fluctuations in the radius of gyration are suppressed when compared to the flexible case. The shear thinning is greater and, as the rigidity is increased, the viscosity measurements tend to those for a multibead rod.     

Alignment of lamellar diblock copolymer phases under shear: insight from dissipative particle dynamics simulations

Sheared self-assembled lamellar phases formed by symmetrical diblock copolymers are investigated through dissipative particle dynamics simulations. Our intent is to provide insight into the experimental observations that the lamellar phases adopt parallel alignment at low shear rates and perpendicular alignment at high shear rates and that it is possible to use shear to induce a transition from the parallel to perpendicular alignment. Simulations are initiated either from lamellar structures prepared under zero shear where lamellae are aligned into parallel, perpendicular, or transverse orientations with respect to the shear direction or from a disordered melt obtained by energy minimization of a random structure. We first consider the relative stability of the parallel and perpendicular phases by applying shear to lamellar structures initially aligned parallel and perpendicular to the shear direction, respectively. The perpendicular lamellar phase persists for all shear rates investigated, whereas the parallel lamellar phase is only stable at low shear rates, and it becomes unstable at high shear rates. At the high shear rates, the parallel lamellar phase first transforms into an unstable diagonal lamellar phase; and upon further increase of the shear rate, the parallel lamellar phase reorients into a perpendicular alignment. We further determine the preferential alignment of the lamellar phases at low shear rate by performing the simulations starting from either the initial transverse lamellar structure or the disordered melt. Since the low shear-rate simulations are plagued by the unstable diagonal lamellar phases, we vary the system size to achieve the natural spacing of the lamellae in the simulation box. In such cases, the unstable diagonal lamellar phases disappear and lamellar phases adopt the preferential alignment, either parallel or perpendicular. In agreement with the experimental observations, the simulations show that the lamellar phase preferentially adopts the parallel orientation at low shear rates and the perpendicular orientation at high shear rates. The simulations further reveal that the perpendicular lamellar phase has lower internal energy than the parallel lamellar phase, whereas the entropy production of the perpendicular lamellar phase is higher with respect to the parallel lamellar phase. Values of the internal energy and entropy production for the unstable diagonal lamellar phases lie between the corresponding values for the parallel and perpendicular lamellar phases. These simulation results suggest that the relative stability of the parallel and perpendicular lamellar phases at low shear rates is a result of the interplay between competing driving forces in the system: (a) the system's drive to adopt a structure with the lowest internal energy and (b) the system's drive to stay in a stationary nonequilibrium state with the lowest entropy production.     

Dynamical and structural properties of benzene in supercritical water

We have employed an anisotropic united atom model of benzene (R. O. Contreras, Ph.D. thesis, Universitat Rovira i Virgili 2002) that reproduces the quadrupolar moment of this molecule through the inclusion of seven point charges. We show that this kind of interaction is required to reproduce the solvation of these molecules in supercritical water. We have computed self-diffusion coefficient and Maxwell-Stefan coefficients as well as the shear viscosity for the mixture water-benzene at supercritical conditions. A strong density and composition dependence of these properties is observed. In addition, our simulations are in qualitative agreement with the experimental evidence that, at medium densities (0.6 g/cm(3) and 673 K), almost half of the benzene molecules have one hydrogen bond with water molecules. We also observe that these bonds are longer lived than the corresponding hydrogen bonds between water molecules. Similarly, we obtain an important reduction of the dielectric constant of the mixture with the increment of the amount of benzene molecules at medium and high densities.     

Examining the frictional forces between mixed hydrophobic-hydrophilic alkylsilane monolayers

Monolayers presenting methyl-terminated (hydrophobic) and hydroxyl-terminated (hydrophilic) surfaces on silica have been studied by molecular dynamics simulation and the effects of hydrogen bonding, chain length, and chain mixing on the frictional properties determined. The hydroxyl-terminated monolayers were found to show large adhesion zones as a result of strong interfacial interlayer hydrogen bonds; the interfacial sliding forces observed in the hydroxyl-terminated monolayers being one order of magnitude higher than the interfacial forces for the hydrophobic surfaces at the characteristic point of zero-load. Mixed hydroxyl- and methyl-terminated monolayers of equal length were found to exhibit intermediate shear stress values between those observed for pure monolayers, with the magnitude of the shear stress depending on the surface content of the hydroxyl-terminated chains. For mixed monolayers of unequal chain lengths, at high loads a maximum in the magnitude of the shear stress as a function of the length of the methyl-terminated chain was observed due to the creation of a buffer zone between the hydroxyl-terminated chains that produces strong hydrogen-bonding interactions. The effect of a constant normal load or constant separation simulation ensemble on the results has also been studied and in general found to have minimal influence on the observed behavior, although some differences are observed for the shear stress at intermediate normal loads due to the formation of stronger hydrogen bond networks at constant load compared to constant separation.     

Effects of electrolyte concentration and counterion valence on the microstructural flow regimes in dilute cetyltrimethylammonium tosylate micellar solutions

The shear thickening behavior and the transition to shear thinning are examined in dilute cetyltrimethylammonium tosylate (CTAT) micellar solutions as a function of surfactant concentration and ionic strength using electrolytes with different counterion valence. Newtonian behavior at low shear rates, followed by shear thickening and shear thinning at higher shear rates, are observed at low and intermediate surfactant and electrolyte concentrations. Shear thickening diminishes with increasing surfactant concentration and ionic strength. At higher surfactant or electrolyte concentration, only a Newtonian region followed by shear thinning is detected. A generalized flow diagram indicates two controlling regimes: one in which electrostatic screening dominates and induces micellar growth, and another, at higher electrolyte and surfactant concentrations, where chemical equilibrium among electrolyte and surfactant counterions controls the rheological behavior by modifying micellar breaking and reforming. Analysis of the shear thickening behavior reveals that not only a critical shear rate is required for shear thickening, but also a critical deformation, which appears to be unique for all systems examined, within experimental error. Moreover, a superposition of the critical shear rate for shear thickening with surfactant and electrolyte concentration is reported.     

Non-linear Viscoelastic Rheological Properties of PCC/PEG Suspensions

The shear thinning and shear thickening rheological properties of PCC/PEG suspension were investigated with the increase of oscillatory amplitude stress at different constant frequencies. The results show that the complex viscosity was initially independent of stress amplitude and obvious shear thinning occurred, then dramatic shear thickening took place after reaching the minimum viscosity. Typically, in a constant frequency of 5 rad/s, the elastic modulus, viscous modulus, and tanδ （δ is the out-of-phase angle） vs. the stress amplitude was investigated. It is found that the elastic modulus initially appeared to be independent of stress amplitude and then exhibited a rapid decrease, but the viscous modulus was independent of amplitude stress at lower amplitude stress. After reaching the minimum value the viscous modulus showed a rapid increase. On the other hand, tanδ increased from 0.6 to 92, which indicates that the transition from elastic to viscous had taken place and tanδ showed a steep increase when shear thickening occurred. Lissajous plots are shown for the dissipated energy vs. different maximum stress amplitude in the shear thinning and shear thickening regions. The relationship of dissipated energy vs. maximum stress amplitude was determined, which follows a power law. In the shear thinning region the exponent was 1.91, but it steeply increases to 3.97 in the shear thickening region.     

Flow-induced structure in a thermotropic liquid crystalline polymer as studied by SANS

Small-angle neutron scattering is utilized to determine the flow induced alignment of a model thermotropic liquid crystalline polymer (LCP) as a function of shear rate and temperature. The results demonstrate that the flow-induced structures in thermotropic liquid crystalline polymers have similarities and differences to those in lyotropic liquid crystalline polymer solutions. The shear rate dependence of the alignment shows that the flow-induced alignment correlates very well to the viscosity behavior of the LCP in the shear thinning regime, while temperature variation results in a change in the extent of alignment within the nematic phase. Relaxation results also demonstrate that the flow-induced alignment remains essentially unchanged for up to an hour after the shear field has been removed. Last, there exists a regime at low shear rate and low temperature where alignment of the LCP molecule perpendicular to the applied shear flow is stable. These results provide important experimental evidence of the molecular level changes that occur in a thermotropic liquid crystalline polymer during flow, which can be utilized to develop theoretical models and more efficiently process thermotropic polymers. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 3017–3023, 1998     

Simulating infrared spectra and hydrogen bonding in cellulose Iβ at elevated temperatures

We have modeled the transformation of cellulose Iβ to a high temperature (550 K) structure, which is considered to be the first step in cellulose pyrolysis. We have performed molecular dynamics simulations at constant pressure using the GROMOS 45a4 united atom forcefield. To test the forcefield, we computed the density, thermal expansion coefficient, total dipole moment, and dielectric constant of cellulose Iβ, finding broad agreement with experimental results. We computed infrared (IR) spectra of cellulose Iβ over the range 300-550 K as a probe of hydrogen bonding. Computed IR spectra were found to agree semi-quantitatively with experiment, especially in the O-H stretching region. We assigned O-H stretches using a novel synthesis of normal mode analysis and power spectrum methods. Simulated IR spectra at elevated temperatures suggest a structural transformation above 450 K, a result in agreement with experimental IR results. The low-temperature (300-400 K) structure of cellulose Iβ is dominated by intrachain hydrogen bonds, whereas in the high-temperature structure (450-550 K), many of these transform to longer, weaker interchain hydrogen bonds. A three-dimensional hydrogen bonding network emerges at high temperatures due to formation of new interchain hydrogen bonds, which may explain the stability of the cellulose structure at such high temperatures.     

Preferential solvation in urea solutions at different concentrations: properties from simulation studies

We performed molecular dynamics simulations of urea solutions at different concentrations with two urea models (OPLS and KBFF) to examine the structures responsible for the thermodynamic solution properties. Our simulation results showed that hydrogen-bonding properties such as the average number of hydrogen bonds and their lifetime distributions were nearly constant at all concentrations between infinite dilution and the solubility limit. This implies that the characterization of urea-water solutions in the molarity concentration scale as nearly ideal is a result of facile local hydrogen bonding rather than a global property. Thus, urea concentration does not influence the local propensity for hydrogen bonds, only how they are satisfied. By comparison, the KBFF model of urea donated fewer hydrogen bonds than OPLS. We found that the KBFF urea model in TIP3P water better reproduced the experimental density and diffusion constant data. Preferential solvation analysis showed that there were weak urea-urea and water-water associations in OPLS solution at short distances, but there were no strong associations. We divided urea molecules into large, medium, and small clusters to examine fluctuation properties and found that any particular urea molecule did not stay in the same cluster for a long time. We found neither persistent nor large clusters.     

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. II. Shear modulus, relaxation-elasticity connections, and rheology

We apply the quiescent and mechanically driven versions of nonlinear Langevin equation theory to study how particle softness influences the shear modulus, the connection between shear elasticity and activated relaxation, and nonlinear rheology of the repulsive Hertzian contact model of dense soft sphere fluids. Below the soft jamming threshold, the shear modulus follows a power law dependence on volume fraction over a narrow interval with an apparent exponent that grows with particle stiffness. To a first approximation, the elastic modulus and transient localization length are controlled by a single coupling constant determined by local fluid structure. In contrast to the behavior of hard spheres, an approximately linear relation between the shear modulus and activation barrier is predicted. This connection has recently been observed for microgel suspensions and provides a microscopic realization of the elastic shoving model. Yielding, shear and stress thinning of the alpha relaxation time and viscosity, and flow curves are also studied. Yield strains are relatively weakly dependent on volume fraction and particle stiffness. Shear thinning commences at values of the effective Peclet number far less than unity, a signature of stress-assisted activated relaxation when barriers are high. Apparent power law reduction of the viscosity with shear rate is predicted with a thinning exponent less than unity. In the vicinity of the soft jamming threshold, a power law flow curve occurs over an intermediate reduced shear rate range with an apparent exponent that decreases as fluid volume fraction and/or repulsion strength increase.     

Rheological properties of microfibrillar suspension of TEMPO-oxidized pulp

The rheology of cellulose microfibril suspensions from TEMPO-oxidized pulp was investigated. The suspension showed a pseudo-plastic and thixotropic behavior, slowly evolving with time under a given shear rate. The viscosity was proportional to the concentration up to the critical concentration of 0.23%. Above it, the viscosity followed a power law with exponents from 2 to 6 depending on the shear rate, and the system showed shear thinning behavior and behaved gel-like. Below this concentration, the system was more Newtonian. Birefringence measurement of 0.44% and 0.78% suspension showed that microfibrils alignment saturated at a small shear rate with a Herman’s orientation parameter below 0.65 probably due to the interconnection of microfibrils.     

Structure and rheological behavior of highly charged colloidal particles in a cylindrical pore I. Effect of pore size

In this work we performed nonequilibrium Brownian dynamics (NEBD) computer simulations of highly charged colloidal particles in diluted suspension under a parabolic flow in cylindrical pores. The influence of charged and neutral cylindrical pores on the structure and rheology of suspensions is analyzed. A shear-induced disorder-order-disorder-like transition was monitored for low shear rates and small pore diameters. We calculate the concentration profiles, axial distribution functions, and axial-angular pair correlation functions to determine the structural properties at steady state for a constant shear flow for different pore sizes and flow strengths. Similar behavior has been observed in a planar narrow channel in the case of charged interacting colloidal particles (M.A. Valdez, O. Manero, J. Colloid Interface Sci. 190 (1997) 81). The mobility of the particles in the radial direction decreases rapidly with the flow and becomes practically frozen. The flow exhibits non-Newtonian shear thinning behavior due to interparticle interactions and particle-wall interaction; the apparent viscosity is lower as the pore diameter decreases, giving rise to an apparent slip in the colloidal suspension. The calculated slip velocity was higher than that obtained in a rectangular slit under shear flow.     

A new lattice-based theory for hydrogen-bonding liquids in uniform electric fields

We propose a new lattice-based, mean-field theory for predicting alignment of molecular dipoles and hydrogen bonds in liquids subject to uniform electric fields. The theory is presently restricted to liquids whose molecules possess one (proton) donor and one acceptor sites each, and wherein the H-bond axis is collinear with the dipole moments of the bonded molecules. The final expressions for hydrogen bond stoichiometry and polarization are free of lattice parameters, are interpretable using simple phenomenological arguments, and reduce to known limiting forms. The theory is applied to understand the internal structure of hydrogen cyanide in the liquid state at different electric fields.     

Simulations of liquid crystals in Poiseuille flow

Lattice Boltzmann simulations are used to explore the behavior of liquid crystals subject to Poiseuille flow. In the nematic regime at low shear rates we find two possible steady-state configurations of the director field. The selected state depends on both the shear rate and the history of the sample. For both director configurations there is clear evidence of shear thinning, a decrease in the viscosity with increasing shear rate. Moreover, at very high shear rates or when the order parameter is large, the system transforms to a ‘log-rolling state’ with boundary layers that may exhibit oscillatory behavior.     

Molecular simulation of structure formation and rheological properties of mixtures of telechelic and monofunctional associating polymer

We study the phase behavior, the characteristic times, and the rheological properties under the steady shear flow of the mixtures consisting of telechelic and monofunctional associating polymers by a coarse‐grained molecular dynamics simulation. The mixtures form the transient networks, the closely packed spherical micelles, and the wormlike micelles. We confirm the molecular origins of the several characteristic times of the mixtures. The dependencies of the characteristic times on the composition ratio between telechelic and monofunctional associating polymers show good agreement with reported experimental results. Under the steady shear flow, the mixtures show the shear thinning induced by the change of the spatial configuration of the micelles. The telechelic associating polymers especially play an important role in connecting the micelles at the shear thinning regime and enhance the steady shear viscosity. Furthermore, at the wormlike micellar region, the mixtures show the second shear thinning initiated by the transformation of the association conformation of the telechelic associating polymers.     

Ab initio calculation of the adhesion and ideal shear strength of planar diamond interfaces with different atomic structure and hydrogen coverage

We propose a method to calculate the ideal shear strength τ of two surfaces in contact by ab initio calculations. This quantity and the work of adhesion γ are the interfacial parameters usually derived from tip-based friction force measurements. We consider diamond interfaces and quantitatively evaluate the effects of surface orientation and passivation. We find that in the case of fully passivated interfaces, γ is not affected by the orientation and the alignment of the surfaces in contact. On the contrary, τ does show a dependence on the atomic-scale roughness of the interface. The surface termination has a major impact on the tribological properties of diamond. The presence of dangling bonds, even at concentrations low enough to prevent the formation of interfacial C-C bonds, causes an increase in the resistance to sliding by 2 orders of magnitude with respect to the fully hydrogenated case. We discuss our findings in relation to experimental observations.