Shear-stress-controlled dynamics of nematic complex fluids

Complex dynamics of a sheared nematic fluid

Nonlinearities in constitutive equations of extended objects in shear flow lead to novel phenomena, e.g. 'rheochaos' in solutions of wormlike micelles and 'elastic turbulence' in polymer solutions. Since both phenomena involve anisotropic objects, their contributions to the deviatoric stress are likely to be similar. However, these two fields have evolved rather independently and an attempt at connecting these fields is still lacking. We show that a minimal model in which the anisotropic nature of the constituting objects is taken into account by a nematic alignment tensor field reproduces several statistical features found in rheochaos and elastic turbulence. We numerically analyse the full non-linear hydrodynamic equations of a sheared nematic fluid under shear stress and strain rate controlled situations, incorporating spatial heterogeneity only in the gradient direction. For a certain range of imposed stress and strain rates, this extended dynamical system shows signatures of spatiotemporal chaos and transient shear banding. In the chaotic regime the power spectra of the order parameter stress, velocity fluctuations and the total injected power show power law behavior and the total injected power shows a non-gaussian, skewed probability distribution. These dynamical features bear resemblance to elastic turbulence phenomena observed in polymer solutions. The scaling behavior is independent of the choice of shear rate/stress controlled method.

Analyzing spinodal decomposition of an anisotropic fluid mixture

Spinodal decomposition leads to spontaneous fluctuations of the local concentration. In the early stage, the resulting pattern provides explicit information about the material properties of the mixture. In the case of two isotropic fluids, the static structure factor shows the characteristic ring shape. If one component is a liquid crystal, the pattern is typically anisotropic and the structure factor is more complex. Using numerical methods, we investigate how structure factors can be used to extract information about material properties like the diffusion constant or the isotropic and the anisotropic contributions to the interfacial tension. The method is based on momenta taken from structure factors in the early stage of the spinodal demixing. We perform phase field calculations for an isotropic and an anisotropic spinodal decomposition. A comparison of the extracted results with analytic values is made. The calculations show that linear modes dominate in the beginning of the growth process, while non-linear modes grow monotonously in the region of the k-space for which damping is predicted by the linearized theory. As long as non-linear modes are small enough, linearized theory can be applied to extract material properties from the structure factor.

The Connection between Biaxial Orientation and Shear Thinning for Quasi-Ideal Rods

The complete orientational ordering tensor of quasi-ideal colloidal rods is obtained as a function of shear rate by performing rheo-SANS (rheology with small angle neutron scattering) measurements on isotropic fd-virus suspensions in the two relevant scattering planes, the flow-gradient (1-2) and the flow-vorticity (1-3) plane. Microscopic ordering can be identified as the origin of the observed shear thinning. A qualitative description of the rheological response by Smoluchowski, as well as Doi⁻Edwards⁻Kuzuu theory is possible, as we obtain a master curve for different concentrations, scaling the shear rate with the apparent collective rotational diffusion coefficient. However, the observation suggests that the interdependence of ordering and shear thinning at small shear rates is stronger than predicted. The extracted zero-shear viscosity matches the concentration dependence of the self-diffusion of rods in semi-dilute solutions, while the director tilts close towards the flow direction already at very low shear rates. In contrast, we observe a smaller dependence on the shear rate in the overall ordering at high shear rates, as well as an ever-increasing biaxiality.

Binary mixtures of rod-like colloids under shear: microscopically-based equilibrium theory and order-parameter dynamics

This paper is concerned with the dynamics of a binary mixture of rod-like, repulsive colloidal particles driven out of equilibrium by means of a steady shear flow (Couette geometry). To this end we first derive, starting from a microscopic density functional in Parsons-Lee approximation, a mesoscopic free energy functional whose main variables are the orientational order parameter tensors. Based on this mesoscopic functional we then explore the stability of isotropic and nematic equilibrium phases in terms of composition and rod lengths. Second, by combining the equilibrium theory with the Doi-Hess approach for the order parameter dynamics under shear, we investigate the orientational dynamics of binary mixtures for a range of shear rates and coupling parameters. We find a variety of dynamical states, including synchronized oscillatory states of the two components, but also symmetry breaking behavior where the components display different in-plane oscillatory states.

Amplitude-phase coupling drives chimera states in globally coupled laser networks

For a globally coupled network of semiconductor lasers with delayed optical feedback, we demonstrate the existence of chimera states. The domains of coherence and incoherence that are typical for chimera states are found to exist for the amplitude, phase, and inversion of the coupled lasers. These chimera states defy several of the previously established existence criteria. While chimera states in phase oscillators generally demand nonlocal coupling, large system sizes, and specially prepared initial conditions, we find chimera states that are stable for global coupling in a network of only four coupled lasers for random initial conditions. The existence is linked to a regime of multistability between the synchronous steady state and asynchronous periodic solutions. We show that amplitude-phase coupling, a concept common in different fields, is necessary for the formation of the chimera states.

Manipulating shear-induced non-equilibrium transitions in colloidal films by feedback control

Using Brownian Dynamics (BD) simulations we investigate non-equilibrium transitions of sheared colloidal films under controlled shear stress σxz. In our approach the shear rate [small gamma, Greek, dot above] is a dynamical variable, which relaxes on a time scale τc such that the instantaneous, configuration-dependent stress σxz(t) approaches a pre-imposed value. Investigating the dynamics under this "feedback-control" scheme we find unique behavior in regions where the flow curve σxz([small gamma, Greek, dot above]) of the uncontrolled system is monotonic. However, in non-monotonic regions our method allows to select between dynamical states characterized by different in-plane structure and viscosities. Indeed, the final state strongly depends on τc relative to an intrinsic relaxation time of the uncontrolled system. The critical values of τc are estimated on the basis of a simple model.

Feedback control of flow alignment in sheared liquid crystals

Based on a continuum theory, we investigate the manipulation of the nonequilibrium behavior of a sheared liquid crystal via closed-loop feedback control. Our goal is to stabilize a specific dynamical state, that is, the stationary "flow alignment," under conditions where the uncontrolled system displays oscillatory director dynamics with in-plane symmetry. To this end we employ time-delayed feedback control (TDFC), where the equation of motion for the ith component q(i)(t) of the order parameter tensor is supplemented by a control term involving the difference q(i)(t)-q(i)(t-τ). In this diagonal scheme, τ is the delay time. We demonstrate that the TDFC method successfully stabilizes flow alignment for suitable values of the control strength K and τ; these values are determined by solving an exact eigenvalue equation. Moreover, our results show that only small values of K are needed when the system is sheared from an isotropic equilibrium state, contrary to the case where the equilibrium state is nematic.

Oscillatory motion of sheared nanorods beyond the nematic phase

We study the role of the control parameter triggering nematic order (temperature or concentration) on the dynamical behavior of a system of nanorods under shear. Our study is based on a set of mesoscopic equations of motion for the components of the tensorial orientational order parameter. We investigate these equations via a systematic bifurcation analysis based on a numerical continuation technique, focusing on spatially homogeneous states. Exploring a wide range of parameters we find, unexpectedly, that states with oscillatory motion can exist even under conditions where the equilibrium system is isotropic. These oscillatory states are characterized by a wagging motion of the paranematic director, and they occur if the tumbling parameter is sufficiently small. We also present full nonequilibrium phase diagrams in the plane spanned by the concentration and the shear rate.

Helicopter rotation and smectic-isotropic coexistence of strongly attractive rods

Hydrodynamic simulations of strongly attractive rodlike colloids are performed with and without shear flow. In the absence of flow, the isotropic-nematic coexistence becomes isotropic-smectic A, and the interfacial properties clearly vary with increasing attraction strength. In the presence of shear flow, a new collective rotation appears in which the director rotates in the vorticity-flow plane in a similar fashion to the movement of the rotor of a helicopter.

Nonlinear rheology of active particle suspensions: insights from an analytical approach

We consider active suspensions in the isotropic phase subjected to a shear flow. Using a set of extended hydrodynamic equations we derive a variety of analytical expressions for rheological quantities such as shear viscosity and normal stress differences. In agreement to full-blown numerical calculations and experiments we find a shear-thickening or -thinning behavior depending on whether the particles are contractile or extensile. Moreover, our analytical approach predicts that the normal stress differences can change their sign in contrast to passive suspensions.