Colloidal particles or nanoparticles, with equal affinity for two fluids, are known to adsorb irreversibly to the fluid-fluid interface. We present large-scale computer simulations of the demixing of a binary solvent containing such particles. The newly formed interface sequesters the colloidal particles; as the interface coarsens, the particles are forced into close contact by interfacial tension. Coarsening is markedly curtailed, and the jammed colloidal layer seemingly enters a glassy state, creating a multiply connected, solidlike film in three dimensions. The resulting gel contains percolating domains of both fluids, with possible uses as, for example, a microreaction medium.

We establish the complete phase diagram of self-propelled hard disks in two spatial dimensions from the analysis of the equation of state and the statistics of local order parameters. The equilibrium melting scenario is maintained at small activities, with coexistence between active liquid and hexatic order, followed by a proper hexatic phase, and a further transition to an active solid. As activity increases, the emergence of hexatic and solid order is shifted towards higher densities. Above a critical activity and for a certain range of packing fractions, the system undergoes motility-induced phase separation and demixes into low and high density phases; the latter can be either disordered (liquid) or ordered (hexatic or solid) depending on the activity.

We present a generic mechanism by which reproducing microorganisms, with a diffusivity that depends on the local population density, can form stable patterns. For instance, it is known that a decrease of bacterial motility with density can promote separation into bulk phases of two coexisting densities; this is opposed by the logistic law for birth and death that allows only a single uniform density to be stable. The result of this contest is an arrested nonequilibrium phase separation in which dense droplets or rings become separated by less dense regions, with a characteristic steady-state length scale. Cell division predominates in the dilute regions and cell death in the dense ones, with a continuous flux between these sustained by the diffusivity gradient. We formulate a mathematical model of this in a case involving run-and-tumble bacteria and make connections with a wider class of mechanisms for density-dependent motility. No chemotaxis is assumed in the model, yet it predicts the formation of patterns strikingly similar to some of those believed to result from chemotactic behavior.

When a large set of discrete bodies passes through a bottleneck, the flow may become intermittent due to the development of clogs that obstruct the constriction. Clogging is observed, for instance, in colloidal suspensions, granular materials and crowd swarming, where consequences may be dramatic. Despite its ubiquity, a general framework embracing research in such a wide variety of scenarios is still lacking. We show that in systems of very different nature and scale -including sheep herds, pedestrian crowds, assemblies of grains, and colloids- the probability distribution of time lapses between the passages of consecutive bodies exhibits a power-law tail with an exponent that depends on the system condition. Consequently, we identify the transition to clogging in terms of the divergence of the average time lapse. Such a unified description allows us to put forward a qualitative clogging state diagram whose most conspicuous feature is the presence of a length scale qualitatively related to the presence of a finite size orifice. This approach helps to understand paradoxical phenomena, such as the faster-is-slower effect predicted for pedestrians evacuating a room and might become a starting point for researchers working in a wide variety of situations where clogging represents a hindrance.

The late-stage demixing following spinodal decomposition of a three-dimensional symmetric binary fluid mixture is studied numerically, using a thermodynamically consistent lattice Boltzmann method. We combine results from simulations with different numerical parameters to obtain an unprecedented range of length and time scales when expressed in reduced physical units. (These are the length and time units derived from fluid density, viscosity, and interfacial tension.) Using eight large (256³) runs, the resulting composite graph of reduced domain size l against reduced time t covers 1 [lsim ] l [lsim ] 10⁵, 10 [lsim ] t [lsim ] 10⁸. Our data are consistent with the dynamical scaling hypothesis that l(t) is a universal scaling curve. We give the first detailed statistical analysis of fluid motion, rather than just domain evolution, in simulations of this kind, and introduce scaling plots for several quantities derived from the fluid velocity and velocity gradient fields. Using the conventional definition of Reynolds number for this problem, Re_ϕ = ldl/dt, we attain values approaching 350. At Re_ϕ [gsim ] 100 (which requires t [gsim ] 10⁶) we find clear evidence of Furukawa's inertial scaling (l ∼ t^2/3), although the crossover from the viscous regime (l ∼ t) is both broad and late (10² [lsim ] t [lsim ] 10⁶). Though it cannot be ruled out, we find no indication that Re_ϕ is self-limiting (l ∼ t^1/2) at late times, as recently proposed by Grant & Elder. Detailed study of the velocity fields confirms that, for our most inertial runs, the RMS ratio of nonlinear to viscous terms in the Navier–Stokes equation, R₂, is of order 10, with the fluid mixture showing incipient turbulent characteristics. However, we cannot go far enough into the inertial regime to obtain a clear length separation of domain size, Taylor microscale, and Kolmogorov scale, as would be needed to test a recent ‘extended’ scaling theory of Kendon (in which R₂ is self-limiting but Re_ϕ not). Obtaining our results has required careful steering of several numerical control parameters so as to maintain adequate algorithmic stability, efficiency and isotropy, while eliminating unwanted residual diffusion. (We argue that the latter affects some studies in the literature which report l ∼ t^2/3 for t [lsim ] 10⁴.) We analyse the various sources of error and find them just within acceptable levels (a few percent each) in most of our datasets. To bring these under significantly better control, or to go much further into the inertial regime, would require much larger computational resources and/or a breakthrough in algorithm design.

We introduce a dissipative particle dynamics scheme for the dynamics of nonideal fluids. Given a free-energy density that determines the thermodynamics of the system, we derive consistent conservative forces. The use of these effective, density dependent forces reduces the local structure as compared to previously proposed models. This is an important feature in mesoscopic modeling, since it ensures a realistic length and time scale separation in coarse-grained models. We consider in detail the behavior of a van der Waals fluid and a binary mixture with a miscibility gap. We discuss the physical implications of having a single length scale characterizing the interaction range, in particular for the interfacial properties.