Spinodal decomposition to a lamellar phase: effects of hydrodynamic flow

Results are presented for the kinetics of domain growth of a two-dimensional fluid quenched from a disordered to a lamellar phase. At early times when a Lifshitz-Slyozov mechanism is operative the growth process proceeds logarithmically in time to a frozen state with locked-in defects. However when hydrodynamic modes become important, or the fluid is subjected to shear, the frustration of the system is alleviated and the size and orientation of the lamellae attain their equilibrium values.Comment: 4 Revtex pages, 4 figures, to appear in Physical Review Letter

Switching dynamics in cholesteric blue phases

Blue phases are networks of disclination lines, which occur in cholesteric liquid crystals near the transition to the isotropic phase. They have recently been used for the new generation of fast switching liquid crystal displays. Here we study numerically the steady states and switching hydrodynamics of blue phase I (BPI) and blue phase II (BPII) cells subjected to an electric field. When the field is on, there are three regimes: for very weak fields (and strong anchoring at the boundaries) the blue phases are almost unaffected, for intermediate fields the disclinations twist (for BPI) and unzip (for BPII), whereas for very large voltages the network dissolves in the bulk of the cell. Interestingly, we find that a BPII cell can recover its original structure when the field is switched off, whereas a BPI cell is found to be trapped more easily into metastable configurations. The kinetic pathways followed during switching on and off entails dramatic reorganisation of the disclination networks. We also discuss the effect of changing the director field anchoring at the boundary planes and of varying the direction of the applied field.Comment: 17 pages, 11 figure

Monte Carlo study of the growth of striped domains

We analyze the dynamical scaling behavior in a two-dimensional spin model with competing interactions after a quench to a striped phase. We measure the growth exponents studying the scaling of the interfaces and the scaling of the shrinking time of a ball of one phase plunged into the sea of another phase. Our results confirm the predictions found in previous papers. The correlation functions measured in the direction parallel and transversal to the stripes are different as suggested by the existence of different interface energies between the ground states of the model. Our simulations show anisotropic features for the correlations both in the case of single-spin-flip and spin-exchange dynamics.Comment: 15 pages, ReVTe

Correlation functions by Cluster Variation Method for Ising model with NN, NNN and Plaquette interactions

We consider the procedure for calculating the pair correlation function in the context of the Cluster Variation Methods. As specific cases, we study the pair correlation function in the paramagnetic phase of the Ising model with nearest neighbors, next to the nearest neighbors and plaquette interactions in two and three dimensions. In presence of competing interactions, the so called disorder line separates in the paramagnetic phase a region where the correlation function has the usual exponential behavior from a region where the correlation has an oscillating exponentially damped behavior. In two dimensions, using the plaquette as the maximal cluster of the CVM approximation, we calculate the phase diagram and the disorder line for a case where a comparison is possible with results known in literature for the eight-vertex model. In three dimensions, in the CVM cube approximation, we calculate the phase diagram and the disorder line in some cases of particular interest. The relevance of our results for experimental systems like mixtures of oil, water and surfactant is also discussed.Comment: 31 pages, LaTeX file, 7 figure

Bistable defect structures in blue phase devices

Blue phases (BPs) are liquid crystals made up by networks of defects, or disclination lines. While existing phase diagrams show a striking variety of competing metastable topologies for these networks, very little is known as to how to kinetically reach a target structure, or how to switch from one to the other, which is of paramount importance for devices. We theoretically identify two confined blue phase I systems in which by applying an appropriate series of electric field it is possible to select one of two bistable defect patterns. Our results may be used to realise new generation and fast switching energy-saving bistable devices in ultrathin surface treated BPI wafers.Comment: 4 pages, 3 figures. Accepted for publication in Phys. Rev. Let

Corner transport upwind lattice Boltzmann model for bubble cavitation

Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann (LB) model that describes a two-dimensional ($2D$) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner transport upwind (CTU) numerical scheme on large square lattices (up to $6144 \times 6144$ nodes). The numerical viscosity and the regularization of the model are discussed for first and second order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows to recover the solution of the $2D$ Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient $D$ and the capillary number $Ca$ is found at small $Ca$ but with a different factor than in equilibrium liquids. A non-linear regime is observed for $Ca \gtrsim 0.2$.Comment: Accepted for publication in Phys. Rev.

Cavitation inception of a van der Waals fluid at a sack-wall obstacle

Cavitation in a liquid moving past a constraint is numerically investigated by means of a free-energy lattice Boltzmann simulation based on the van der Waals equation of state. The fluid is streamed past an obstacle and, depending on the pressure drop between inlet and outlet, vapor formation underneath the corner of the sack-wall is observed. The circumstances of cavitation formation are investigated and it is found that the local bulk pressure and mean stress are insufficient to explain the phenomenon. Results obtained in this study strongly suggest that the viscous stress, interfacial contributions to the local pressure, and the Laplace pressure are relevant to the opening of a vapor cavity. This can be described by a generalization of Joseph's criterion that includes these contributions. A macroscopic investigation measuring mass flow rate behavior and discharge coefficient was also performed. As theoretically predicted, mass flow rate increases linearly with the square root of the pressure drop. However, when cavitation occurs, the mass flow growth rate is reduced and eventually it collapses into a choked flow state. In the cavitating regime, as theoretically predicted and experimentally verified, the discharge coefficient grows with the Nurick cavitation number