Ordering of the lamellar phase under a shear flow

Complex-tensor theory of simple smectics

Matter self-assembling into layers generates unique properties, including structures of stacked surfaces, directed transport, and compact area maximization that can be highly functionalized in biology and technology. Smectics represent the paradigm of such lamellar materials - they are a state between fluids and solids, characterized by both orientational and partial positional ordering in one layering direction, making them notoriously difficult to model, particularly in confining geometries. We propose a complex tensor order parameter to describe the local degree of lamellar ordering, layer displacement and orientation of the layers for simple, lamellar smectics. The theory accounts for both dislocations and disclinations, by regularizing singularities within defect cores and so remaining continuous everywhere. The ability to describe disclinations and dislocation allows this theory to simulate arrested configurations and inclusion-induced local ordering. This tensorial theory for simple smectics considerably simplifies numerics, facilitating studies on the mesoscopic structure of topologically complex systems.

Long-Range Phase Order in Two Dimensions under Shear Flow

We theoretically and numerically investigate a two-dimensional O(2) model where an order parameter is convected by shear flow. We show that a long-range phase order emerges in two dimensions as a result of anomalous suppression of phase fluctuations by the shear flow. Furthermore, we use the finite-size scaling theory to demonstrate that a phase transition to the long-range ordered state from the disordered state is second order. At a transition point far from equilibrium, the critical exponents turn out to be close to the mean-field value for equilibrium systems.

Shear-induced parallel and transverse alignments of cylinders in thin films of diblock copolymers

Coarse-grained Langevin dynamics simulations were performed to investigate the alignment behavior of monolayer films of cylinder-forming diblock copolymers under steady shear, a structure of significant importance for many technical applications such as nanopatterning. The influences of shear conditions, the interactions involved in the films, and the initial morphology of the cylinder-forming phase were examined. Our results showed that above a critical shear rate, the cylinders can align either along the shearing direction or transverse (log-rolling) to the shearing direction depending on the relative strength between the interchain attraction in the cylinders (εAA) and the surface attraction of the confining walls with the film (εBW). To understand the underlying mechanism, the microscopic properties of the films under shear were systematically investigated. It was found that at low εAA/εBW, the majority blocks of the diblock polymer that are adsorbed on the confining walls prefer to move synchronously with the walls, inducing the cylinder-forming blocks to align along the flow direction. When εAA/εBW is above a threshold value, a strong attraction between the cylinder-forming blocks restrains their movement during shear, leading to the log-rolling motions of the cylinders. To predict the threshold εAA/εBW, we developed an approach based on equilibrium thermodynamics data and found good agreement with our shear simulations.

Simulation of phase separation with temperature-dependent viscosity using lattice Boltzmann method

This paper presents an exploration of the phase separation behavior and pattern formation in a binary fluid with temperature-dependent viscosity via a coupled lattice Boltzmann method (LBM). By introducing a viscosity-temperature relation into the LBM, the coupling effects of the viscosity-temperature coefficient [Formula: see text] , initial viscosity [Formula: see text] and thermal diffusion coefficient [Formula: see text] , on the phase separation were successfully described. The calculated results indicated that an increase in initial viscosity and viscosity-temperature coefficient, or a decrease in the thermal diffusion coefficient, can lead to the orientation of isotropic growth fronts over a wide range of viscosity. The results showed that droplet-type phase structures and lamellar phase structures with domain orientation parallel or perpendicular to the walls can be obtained in equilibrium by controlling the initial viscosity, thermal diffusivity, and the viscosity-temperature coefficient. Furthermore, the dataset was rearranged for growth kinetics of domain growth and thermal diffusion fronts in a plot by the spherically averaged structure factor and the ratio of separated and continuous phases. The analysis revealed two different temporal regimes: spinodal decomposition and domain growth stages, which further quantified the coupled effects of temperature and viscosity on the evolution of temperature-dependent phase separation. These numerical results provide guidance for setting optimum temperature ranges to obtain expected phase separation structures for systems with temperature-dependent viscosity.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L(3/2) for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects.

Thermal and hydrodynamic effects in the ordering of lamellar fluids

Phase separation in a complex fluid with lamellar order has been studied in the case of cold thermal fronts propagating diffusively from external walls. The velocity hydrodynamic modes are taken into account by coupling the convection-diffusion equation for the order parameter to a generalized Navier-Stokes equation. The dynamical equations are simulated by implementing a hybrid method based on a lattice Boltzmann algorithm coupled to finite difference schemes. Simulations show that the ordering process occurs with morphologies depending on the speed of the thermal fronts or, equivalently, on the value of the thermal conductivity ξ. At large values of ξ, as in instantaneous quenching, the system is frozen in entangled configurations at high viscosity while it consists of grains with well-ordered lamellae at low viscosity. By decreasing the value of ξ, a regime with very ordered lamellae parallel to the thermal fronts is found. At very low values of ξ the preferred orientation is perpendicular to the walls in d=2, while perpendicular order is lost moving far from the walls in d=3.

Shear alignment of a disordered lamellar mesophase

The shear alignment of an initially disordered lamellar phase is examined using lattice Boltzmann simulations of a mesoscopic model based on a free-energy functional for the concentration modulation. For a small shear cell of width 8λ, the qualitative features of the alignment process are strongly dependent on the Schmidt number Sc=ν/D (ratio of kinematic viscosity and mass diffusion coefficient). Here, λ is the wavelength of the concentration modulation. At low Schmidt number, it is found that there is a significant initial increase in the viscosity, coinciding with the alignment of layers along the extensional axis, followed by a decrease at long times due to the alignment along the flow direction. At high Schmidt number, alignment takes place due to the breakage and reformation of layers because diffusion is slow compared to shear deformation; this results in faster alignment. The system size has a strong effect on the alignment process; perfect alignment takes place for a small systems of width 8λ and 16λ, while a larger system of width 32λ does not align completely even at long times. In the larger system, there appears to be a dynamical steady state in which the layers are not perfectly aligned--where there is a balance between the annealing of defects due to shear and the creation due to an instability of the aligned lamellar phase under shear. We observe two types of defect creation mechanisms: the buckling instability under dilation, which was reported earlier, as well as a second mechanism due to layer compression.

Mesoscale description of an asymmetric lamellar phase

The relationship between the parameters in a description based on a mesoscale free energy functional for the concentration field and the macroscopic properties, such as the bending and compression moduli and the permeation constant, are examined for an asymmetric lamellar phase where the mass fractions of the hydrophobic and hydrophilic parts are not equal. The difference in the mass fractions is incorporated using a cubic term in the free energy functional, in addition to the usual quadratic and quartic terms in the Landau-Ginsburg formulation. The relationship between the coefficient of the cubic term and the difference in the mass fractions of the hydrophilic and hydrophobic parts is obtained. For a lamellar phase, it is important to ensure that the surface tension is zero due to symmetry considerations. The relationship between the parameters in the free energy functional for zero surface tension is derived. When the interface between the hydrophilic and hydrophobic parts is diffuse, it is found that the bending and compression moduli, scaled by the parameters in the free energy functional, do increase as the asymmetry in the bilayer increases. When the interface between the hydrophilic and hydrophobic parts is sharp, the scaled bending and compression moduli show no dependence on the asymmetry in the bilayer. The ratio of the permeation constant in between the water and bilayer in a molecular description and the Onsager coefficient in the mesoscale description is O(1) for both sharp and diffuse interfaces and it increases as the difference in the mass fractions is increased.

Multiscale modeling of lamellar mesophases

The mesoscale simulation of a lamellar mesophase based on a free energy functional is examined with the objective of determining the relationship between the parameters in the model and molecular parameters. Attention is restricted to a symmetric lamellar phase with equal volumes of hydrophilic and hydrophobic components. Apart from the lamellar spacing, there are two parameters in the free energy functional. One of the parameters, r, determines the sharpness of the interface, and it is shown how this parameter can be obtained from the interface profile in a molecular simulation. The other parameter, A, provides an energy scale. Analytical expressions are derived to relate these parameters to r and A to the bending and compression moduli and the permeation constant in the macroscopic equation to the Onsager coefficient in the concentration diffusion equation. The linear hydrodynamic response predicted by the theory is verified by carrying out a mesoscale simulation using the lattice-Boltzmann technique and verifying that the analytical predictions are in agreement with simulation results. A macroscale model based on the layer thickness field and the layer normal field is proposed, and the relationship between the parameters in the macroscale model from the parameters in the mesoscale free energy functional is obtained.

Structural transitions and arrest of domain growth in sheared binary immiscible fluids and microemulsions

We investigate spinodal decomposition and structuring effects in binary immiscible and ternary amphiphilic fluid mixtures under shear by means of three-dimensional lattice Boltzmann simulations. We show that the growth of individual fluid domains can be arrested by adding surfactant to the system, thus forming a bicontinuous microemulsion. We demonstrate that the maximum domain size and the time of arrest depend linearly on the concentration of amphiphile molecules. In addition, we find that for a well-defined threshold value of amphiphile concentration, the maximum domain size and time of complete arrest do not change. For systems under constant and oscillatory shear we analyze domain growth rates in directions parallel and perpendicular to the applied shear. We find a structural transition from a sponge to a lamellar phase by applying a constant shear and the occurrence of tubular structures under oscillatory shear. The size of the resulting lamellae and tubes depends strongly on the amphiphile concentration, shear rate, and shear frequency.

Lamellar order, microphase structures, and glassy phase in a field theoretic model for charged colloids

In this paper we present a detailed analytical study of the phase diagram and of the structural properties of a field theoretic model with a short-range attraction and a competing long-range screened repulsion. We provide a full derivation and expanded discussion and digression on results previously reported briefly in M. Tarzia and A. Coniglio, Phys. Rev. Lett. 96, 075702 (2006). The model contains the essential features of the effective interaction potential among charged colloids in polymeric solutions. We employ the self-consistent Hartree approximation and a replica approach, and we show that varying the parameters of the repulsive potential and the temperature yields a phase coexistence, a lamellar and a glassy phase. Our results suggest that the cluster phase observed in charged colloids might be the signature of an underlying equilibrium lamellar phase, hidden on experimental time scales, and emphasize that the formation of microphase structures may play a prominent role in the process of colloidal gelation.

Pattern formation and glassy phase in the phi4 theory with a screened electrostatic repulsion

We study analytically the structural properties of a system with a short-range attraction and a competing long-range screened repulsion. This model contains the essential features of the effective interaction potential among charged colloids in polymeric solutions and provides novel insights on the equilibrium phase diagram of these systems. Within the self-consistent Hartree approximation and by using a replica approach, we show that varying the parameters of the repulsive potential and the temperature yields a phase coexistence, a lamellar, and a glassy phase. Our results strongly suggest that the cluster phase observed in charged colloids might be the signature of an underlying equilibrium lamellar phase, hidden on experimental time scales.

Large-scale grid-enabled lattice Boltzmann simulations of complex fluid flow in porous media and under shear

Well-designed lattice Boltzmann codes exploit the essentially embarrassingly parallel features of the algorithm and so can be run with considerable efficiency on modern supercomputers. Such scalable codes permit us to simulate the behaviour of increasingly large quantities of complex condensed matter systems. In the present paper, we present some preliminary results on the large-scale three-dimensional lattice Boltzmann simulation of binary immiscible fluid flows through a porous medium, derived from digitized X-ray micro-tomographic data of Bentheimer sandstone, and from the study of the same fluids under shear. Simulations on such scales can benefit considerably from the use of computational steering, and we describe our implementation of steering within the lattice Boltzmann code, called LB3D, making use of the RealityGrid steering library. Our large-scale simulations benefit from the new concept of capability computing, designed to prioritize the execution of big jobs on major supercomputing resources. The advent of persistent computational grids promises to provide an optimal environment in which to deploy these mesoscale simulation methods, which can exploit the distributed nature of computer, visualization and storage resources to reach scientific results rapidly; we discuss our work on the grid-enablement of lattice Boltzmann methods in this context.

Phase-separating binary fluids under oscillatory shear

We apply the lattice Boltzmann methods to study the segregation of binary fluid mixtures under oscillatory shear flow in two dimensions. The algorithm allows to simulate systems whose dynamics is described by the Navier-Stokes and the convection-diffusion equations. The interplay between several time scales produces a rich and complex phenomenology. We investigate the effects of different oscillation frequencies and viscosities on the morphology of the phase separating domains. We find that at high frequencies the evolution is almost isotropic with growth exponents 2/3 and 1/3 in the inertial (low viscosity) and diffusive (high viscosity) regimes, respectively. When the period of the applied shear flow becomes of the same order of the relaxation time T(R) of the shear velocity profile, anisotropic effects are clearly observable. In correspondence with nonlinear patterns for the velocity profiles, we find configurations where lamellar order close to the walls coexists with isotropic domains in the middle of the system. For particular values of frequency and viscosity it can also happen that the convective effects induced by the oscillations cause an interruption or a slowing of the segregation process, as found in some experiments. Finally, at very low frequencies, the morphology of domains is characterized by lamellar order everywhere in the system resembling what happens in the case with steady shear.