Morphologies and flow patterns in quenching of lamellar systems with shear

Defect interactions in a two-dimensional sheared lamellar mesophase

The interaction between two edge dislocations in a sheared lyotropic lamellar liquid-crystalline medium is examined. The model is a mesoscale hydrodynamic model based on a free-energy functional that is minimised for a sinusoidal concentration modulation coupled with concentration and momentum equations. The defect dynamics are analysed as a function of the system size and the Ericksen number, the ratio of the shear stress and the characteristic free-energy density for deformation. Three different regimes are observed as the Ericksen number is increased when the edge dislocations are sheared towards each other, such that there is compression of layers between the defects: (a) defect motion that reduces the cross-stream separation till there is a steady spacing with plug flow between the defects, (b) defect attraction and cancellation resulting in a well-aligned state, and (c) defect creation due to a compressional instability between the defects resulting in an increase in the disorder. When the edge dislocations are sheared away from each other, such that there is extension of the layers between the defects, two distinct regimes are observed as the Ericksen number is increased: (a) bending of layers and a plug flow between the defects at their initial separation, and (b) buckling of the layers leading to creation of more defects and a dynamical steady state between defect creation and cancellation. These regimes are found to be robust for different values of the system size, from 32 to 128 layers, and for different values of the dimensionless groups that characterise the ratio of mass/momentum convection and diffusion.

The relationship between structure and rheology in a three-dimensional sheared lamellar mesophase

The evolution of a lamellar mesophase from an initially disordered state under shear is examined using simulations of a mesoscale model based on a concentration field ψ that distinguishes the hydrophilic and hydrophobic components. The Landau-Ginzburg free-energy functional is augmented by a term that is minimised for sinusoidal modulations in the concentration field with wavelength λ = (2π/k), and the dynamical equations are the model H equations. The structure and rheology are determined by the relative magnitudes of the diffusion time for coarsening, (λ²/D) and the inverse of the strain rate ^-1, and the Ericksen number, which is the ratio of the shear stress and the layer stiffness. When the diffusion time is small compared with the inverse of the strain rate, there is a local formation of misaligned layers, which are deformed by the imposed flow. There is near-perfect ordering with isolated defects at low values of the Ericksen number, but the defects result in a significant increase in viscosity due to the high layer stiffness. At high values of the Ericksen number, the concentration field is deformed by the mean shear before layers form via diffusion. Cylindrical structures aligned along the flow direction form after about 8-10 strain units, and these evolve into layers with disorder through diffusion perpendicular to the flow. The layers are not perfectly ordered, even after hundreds of strain units, due to the creation and destruction of defects via shear. The excess viscosity is low because the layer stiffness is small compared with the applied shear at a high Ericksen number. This study provides guidance on how the material parameters and imposed flow can be tailored to achieve the desired rheological behaviour.

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.

Coupled lattice Boltzmann-large eddy simulation model for three-dimensional multiphase flows at large density ratio and high Reynolds number

A coupled lattice Boltzmann-large eddy simulation model is developed for modeling three-dimensional multiphase flows at large density ratios and high Reynolds numbers. In the framework of the lattice Boltzmann method, the model is proposed based on the standard Smagorinsky subgrid-scale approach, and a reconstructed multiple-relaxation-time collision operator is adopted. The conservative Allen-Cahn equation and Navier-Stokes equations are solved through the lattice Boltzmann discretization scheme for the interface tracking and velocity field evolution, respectively. Relevant benchmark cases are carried out to validate the performance of this model in simulating multiphase flows at a large density ratio and a high Reynolds number, including a stationary droplet, the process of spinodal decomposition, the Rayleigh-Taylor instability, the phenomenon of a droplet splashing on a thin liquid film, and the liquid jet breakup process. The maximum values of density ratio and Re number are 1000 and 10 240, respectively. The capability and reliability of the proposed model have been demonstrated by the good agreement between simulation results and the analytical solutions or the previously available results.

Multiparticle collision dynamics for fluid interfaces with near-contact interactions

We present an extension of the multiparticle collision dynamics method for flows with complex interfaces, including supramolecular near-contact interactions mimicking the effect of surfactants. The new method is demonstrated for the case of (i) short range repulsion of droplets in close contact, (ii) arrested phase separation, and (iii) different pattern formation during spinodal decomposition of binary mixtures.

Mesoscale modelling of miscible and immiscible multicomponent fluids

A mesoscopic simulation method based on the integration of dissipative particle dynamics (DPD), smoothed particle hydrodynamics (SPH) and computational thermodynamics (CT) has been developed. The kinetic behaviours of miscible and immiscible fluids were investigated. The interaction force between multicomponent mesoscopic particles is derived from the system free energy. The diffusivity of the components in non-ideal solution is determined by the chemical potential. The proposed method provides convincing predictions to the effects of convection, diffusion and microscopic interaction on the non-equilibrium evolution of engineering fluids, and demonstrates a potential to simulate more complicated phenomena in materials processing.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (ργL(2)/μ), the Schmidt number (μ/ρD), the Ericksen number (μγ/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts μ(r), and the ratio of the system size and layer spacing (L/λ). Here, ρ and μ are the fluid density and average viscosity, γ is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, μ(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/λ=32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with "grain boundaries," which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Polar-coordinate lattice Boltzmann modeling of compressible flows

We present a polar coordinate lattice Boltzmann kinetic model for compressible flows. A method to recover the continuum distribution function from the discrete distribution function is indicated. Within the model, a hybrid scheme being similar to, but different from, the operator splitting is proposed. The temporal evolution is calculated analytically, and the convection term is solved via a modified Warming-Beam (MWB) scheme. Within the MWB scheme a suitable switch function is introduced. The current model works not only for subsonic flows but also for supersonic flows. It is validated and verified via the following well-known benchmark tests: (i) the rotational flow, (ii) the stable shock tube problem, (iii) the Richtmyer-Meshkov (RM) instability, and (iv) the Kelvin-Helmholtz instability. As an original application, we studied the nonequilibrium characteristics of the system around three kinds of interfaces, the shock wave, the rarefaction wave, and the material interface, for two specific cases. In one of the two cases, the material interface is initially perturbed, and consequently the RM instability occurs. It is found that the macroscopic effects due to deviating from thermodynamic equilibrium around the material interface differ significantly from those around the mechanical interfaces. The initial perturbation at the material interface enhances the coupling of molecular motions in different degrees of freedom. The amplitude of deviation from thermodynamic equilibrium around the shock wave is much higher than those around the rarefaction wave and material interface. By comparing each component of the high-order moments and its value in equilibrium, we can draw qualitatively the main behavior of the actual distribution function. These results deepen our understanding of the mechanical and material interfaces from a more fundamental level, which is indicative for constructing macroscopic models and other kinds of kinetic models.

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.

High-order thermal lattice Boltzmann models derived by means of Gauss quadrature in the spherical coordinate system

We use the spherical coordinate system in the momentum space and an appropriate discretization procedure to derive a hierarchy of lattice Boltzmann (LB) models with variable temperature. The separation of the integrals in the momentum space into angular and radial parts allows us to compute the moments of the equilibrium distribution function by means of Gauss-Legendre and Gauss-Laguerre quadratures, as well as to find the elements of the discrete momentum set for each LB model in the hierarchy. The capability of the high-order models in this hierarchy to capture specific effects in microfluidics is investigated through a computer simulation of Couette flow by using the Shakhov collision term to get the right value of the Prandtl number.

Phase separation in thermal systems: a lattice Boltzmann study and morphological characterization

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ(0)[(T(c)-T)/(T(c)-T(0))](3/2), where T(c) is the critical temperature and σ(0) is the interfacial stress at a reference temperature T(0), which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

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.

Lattice Boltzmann study on Kelvin-Helmholtz instability: roles of velocity and density gradients

A two-dimensional lattice Boltzmann model with 19 discrete velocities for compressible fluids is proposed. The fifth-order weighted essentially nonoscillatory (5th-WENO) finite difference scheme is employed to calculate the convection term of the lattice Boltzmann equation. The validity of the model is verified by comparing simulation results of the Sod shock tube with its corresponding analytical solutions [G. A. Sod, J. Comput. Phys. 27, 1 (1978).]. The velocity and density gradient effects on the Kelvin-Helmholtz instability (KHI) are investigated using the proposed model. Sharp density contours are obtained in our simulations. It is found that the linear growth rate γ for the KHI decreases by increasing the width of velocity transition layer D(v) but increases by increasing the width of density transition layer D(ρ). After the initial transient period and before the vortex has been well formed, the linear growth rates γ(v) and γ(ρ), vary with D(v) and D(ρ) approximately in the following way, lnγ(v)=a-bD(v) and γ(ρ)=c+elnD(ρ)(D(ρ)<D(ρ)(E)), where a, b, c, and e are fitting parameters and D(ρ)(E) is the effective interaction width of the density transition layer. When D(ρ)>D(ρ)(E) the linear growth rate γ(ρ) does not vary significantly any more. One can use the hybrid effects of velocity and density transition layers to stabilize the KHI. Our numerical simulation results are in general agreement with the analytical results [L. F. Wang et al., Phys. Plasma 17, 042103 (2010)].

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.