No Arabic abstract
We present a lattice Boltzmann algorithm based on an underlying free energy that allows the simulation of the dynamics of a multicomponent system with an arbitrary number of components. The thermodynamic properties, such as the chemical potential of each component and the pressure of the overall system, are incorporated in the model. We derived a symmetrical convection diffusion equation for each component as well as the Navier Stokes equation and continuity equation for the overall system. The algorithm was verified through simulations of binary and ternary systems. The equilibrium concentrations of components of binary and ternary systems simulated with our algorithm agree well with theoretical expectations.
We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles velocity space. Different velocity sets, i.e., higher- and lower-order lattices, are consistently and efficiently coupled, allowing us to use the higher-order lattice only when and where needed. This includes regions of either high Mach number or high Knudsen number. The coupling procedure of different lattices consists of either projection of the moments of the higher-order lattice onto the lower-order lattice or lifting of the lower-order lattice to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme can be formulated both in a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is first validated through a convected athermal vortex and also studied in a jet flow setup. The performance of the proposed scheme is further investigated through the shock structure problem and a high Knudsen Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.
We study liquid-vapor phase separation under shear via the Shan-Chen lattice Boltzmann model. Besides the rheological characteristics, we analyze the Kelvin-Helmholtz(K-H) instability resulting from the tangential velocity difference of the fluids on two sides of the interface. We discuss also the growth behavior of droplets. The domains being close to the walls are lamellar-ordered, where the hydrodynamic effects dominate. The patterns in the bulk of the system are nearly isotropic, where the domain growth results mainly from the diffusion mechanism. Both the interfacial tension and the K-H instability make the liquid-bands near the walls tend to rupture. When the shear rate increases, the inequivalence of evaporation in the upstream and coagulation in the downstream of the flow as well as the role of surface tension makes the droplets elongate obliquely. Stronger convection makes easier the transferring of material particles so that droplets become larger.
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.
It is shown that the Shan-Chen (SC) model for non-ideal lattice fluids can be made compliant with a pseudo free-energy principle by simple addition of a gradient force, whose expression is uniquely specified in terms of the fluid density. This additional term is numerically shown to provide fairly negligible effects on the system evolution during phase-separation. To the best of our knowledge, these important properties of the SC model were not noted before. The approach developed in the present work is based on a continuum analysis: further extensions, more in line with a discrete lattice theory (X. Shan, {it Phys Rev E}, {bf 77} 066702 (2008)) can be envisaged for the future.
We use computer simulations to investigate the stability of a two-component polymer brush de-mixing on a curved template into phases of different morphological properties. It has been previously shown via molecular dynamics simulations that immiscible chains having different length and anchored to a cylindrical template will phase separate into striped phases of different widths oriented perpendicularly to the cylindrical axis. We calculate free energy differences for a variety of stripe widths, and extract simple relationships between the sizes of the two polymers, N_1 and N_2, and the free energy dependence on the stripe width. We explain these relationships using simple physical arguments based upon previous theoretical work on the free energy of polymer brushes.