We present the details of a lattice Boltzmann approach to phase separation in nonideal one- and two-component fluids. The collision rules are chosen such that the equilibrium state corresponds to an input free energy and the bulk flow is governed by the continuity, Navier-Stokes, and, for the binary fluid, a convection-diffusion equation. Numerical results are compared to simple analytic predictions to confirm that the equilibrium state is indeed thermodynamically consistent and that the kinetics of the approach to equilibrium lie within the expected universality classes. The approach is compared to other lattice Boltzmann simulations of nonideal systems. © 1996 The American Physical Society.

We report hybrid lattice Boltzmann (HLB) simulations of the hydrodynamics of an active nematic liquid crystal sandwiched between confining walls with various anchoring conditions. We confirm the existence of a transition between a passive phase and an active phase, in which there is spontaneous flow in the steady state. This transition is attained for sufficiently “extensile” rods, in the case of flow-aligning liquid crystals, and for sufficiently “contractile” ones for flow-tumbling materials. In a quasi-one-dimensional geometry, deep in the active phase of flow-aligning materials, our simulations give evidence of hysteresis and history-dependent steady states, as well as of spontaneous banded flow. Flow-tumbling materials, in contrast, rearrange themselves so that only the two boundary layers flow in steady state. Two-dimensional simulations, with periodic boundary conditions, show additional instabilities, with the spontaneous flow appearing as patterns made up of “convection rolls.” These results demonstrate a remarkable richness (including dependence on anchoring conditions) in the steady-state phase behavior of active materials, even in the absence of external forcing; they have no counterpart for passive nematics. Our HLB methodology, which combines lattice Boltzmann for momentum transport with a finite difference scheme for the order parameter dynamics, offers a robust and efficient method for probing the complex hydrodynamic behavior of active nematics.

We describe a lattice Boltzmann algorithm to simulate liquid crystal hydrodynamics. The equations of motion are written in terms of a tensor order parameter. This allows both the isotropic and the nematic phases to be considered. Backflow effects and the hydrodynamics of topological defects are naturally included in the simulations, as are non-Newtonian flow properties such as shear thinning and shear banding.

Recent experiments showed that the linear double-stranded DNA in bacteriophage capsids is both highly knotted and neatly structured. What is the physical basis of this organization? Here we show evidence from stochastic simulation techniques that suggests that a key element is the tendency of contacting DNA strands to order, as in cholesteric liquid crystals. This interaction favors their preferential juxtaposition at a small twist angle, thus promoting an approximately nematic (and apolar) local order. The ordering effect dramatically impacts the geometry and topology of DNA inside phages. Accounting for this local potential allows us to reproduce the main experimental data on DNA organization in phages, including the cryo-EM observations and detailed features of the spectrum of DNA knots formed inside viral capsids. The DNA knots we observe are strongly delocalized and, intriguingly, this is shown not to interfere with genome ejection out of the phage.