No Arabic abstract
A lattice Boltzmann model for amphiphilic fluid dynamics is presented. It is a ternary model, in that it conserves mass separately for each chemical species present (water, oil, amphiphile), and it maintains an orientational degree of freedom for the amphiphilic species. Moreover, it models fluid interactions at the microscopic level by introducing self-consistent forces between the particles, rather than by positing a Landau free energy functional. This combination of characteristics fills an important need in the hierarchy of models currently available for amphiphilic fluid dynamics, enabling efficient computer simulation and furnishing new theoretical insight. Several computational results obtained from this model are presented and compared to existing lattice-gas model results. In particular, it is noted that lamellar structures, which are precluded by the Peierls instability in two-dimensional systems with kinetic fluctuations, are not observed in lattice-gas models, but are easily found in the corresponding lattice Boltzmann models. This points out a striking difference in the phenomenology accessible to each type of model.
We present an energy-conserving multiple-relaxation-time finite difference lattice Boltzmann model for compressible flows. This model is based on a 16-discrete-velocity model. The collision step is first calculated in the moment space and then mapped back to the velocity space. The moment space and corresponding transformation matrix are constructed according to the group representation theory. Equilibria of the nonconserved moments are chosen according to the need of recovering compressible Navier-Stokes equations through the Chapman-Enskog expansion. Numerical experiments showed that compressible flows with strong shocks can be well simulated by the present model. The used benchmark tests include (i) shock tubes, such as the Sod, Lax, Sjogreen, Colella explosion wave and collision of two strong shocks, (ii) regular and Mach shock reflections, and (iii) shock wave reaction on cylindrical bubble problems. The new model works for both low and high speeds compressible flows. It contains more physical information and has better numerical stability and accuracy than its single-relaxation-time version.
We study the phase behavior of ternary amphiphilic systems in the framework of a curvature model with non-vanishing spontaneous curvature. The amphiphilic monolayers can arrange in different ways to form micellar, hexagonal, lamellar and various bicontinuous cubic phases. For the latter case we consider both single structures (one monolayer) and double structures (two monolayers). Their interfaces are modeled by the triply periodic surfaces of constant mean curvature of the families G, D, P, C(P), I-WP and F-RD. The stability of the different bicontinuous cubic phases can be explained by the way in which their universal geometrical properties conspire with the concentration constraints. For vanishing saddle-splay modulus $bar kappa$, almost every phase considered has some region of stability in the Gibbs triangle. Although bicontinuous cubic phases are suppressed by sufficiently negative values of the saddle-splay modulus $bar kappa$, we find that they can exist for considerably lower values than obtained previously. The most stable bicontinuous cubic phases with decreasing $bar kappa < 0$ are the single and double gyroid structures since they combine favorable topological properties with extreme volume fractions.
The dynamics of dry active matter have implications for a diverse collection of biological phenomena spanning a range of length and time scales, such as animal flocking, cell tissue dynamics, and swarming of inserts and bacteria. Uniting these systems are a common set of symmetries and conservation laws, defining dry active fluids as a class of physical system. Many interesting behaviours have been observed at high densities, which remain difficult to simulate due to the computational demand. Here, we show how two-dimensional dry active fluids in a dense regime can be studied using a simple modification of the lattice Boltzmann method. We apply our method on a model that exhibits motility-induced phase separation, and an active model with contact inhibition of locomotion, which has relevance to collective cell migration. For the latter, we uncover multiple novel phase transitions: two first-order and one potentially critical. We further support our simulation results with an analytical treatment of the hydrodynamic equations obtained via a Chapman-Enskog coarse-graining procedure.
A rigorous free energy model for ternary fluid flows with density ratio up to of order $O(10^3)$ is presented and implemented using the entropic lattice Boltzmann scheme. The model is thermodynamically consistent and allows a broad range of surface tension ratios, covering both partial wetting states where Neumann triangles are formed, and full wetting states where complete encapsulation of one of fluid components is observed. We further demonstrate that we can capture the bouncing, adhesive and insertive regimes for the binary collisions between immiscible droplets suspended in air. Our approach opens up a vast range of multiphase flow applications involving one gas and several liquid components.
We generalize the force-level, microscopic, Nonlinear Langevin Equation (NLE) theory and its elastically collective generalization (ECNLE theory) of activated dynamics in bulk spherical particle liquids to address the influence of random particle pinning on structural relaxation. The simplest neutral confinement model is analyzed for hard spheres where there is no change of the equilibrium pair structure upon particle pinning. As the pinned fraction grows, cage scale dynamical constraints are intensified in a manner that increases with density. This results in the mobile particles becoming more transiently localized, with increases of the jump distance, cage scale barrier and NLE theory mean hopping time; subtle changes of the dynamic shear modulus are predicted. The results are contrasted with recent simulations. Similarities in relaxation behavior are identified in the dynamic precursor regime, including a roughly exponential, or weakly supra-exponential, growth of the alpha time with pinning fraction and a reduction of dynamic fragility. However, the increase of the alpha time with pinning predicted by the local NLE theory is too small, and severely so at very high volume fractions. The strong deviations are argued to be due to the longer range collective elasticity aspect of the problem which is expected to be modified by random pinning in a complex manner. A qualitative physical scenario is offered for how the three distinct aspects that quantify the elastic barrier may change with pinning. ECNLE theory calculations of the alpha time are then presented based on the simplest effective-medium-like treatment for how random pinning modifies the elastic barrier. The results appear to be consistent with most, but not all, trends seen in recent simulations. Key open problems are discussed with regards to both theory and simulation.