No Arabic abstract
We propose a mesoscopic model of binary fluid mixtures with tunable viscosity ratio based on a two-range pseudo-potential lattice Boltzmann method, for the simulation of soft flowing systems. In addition to the short range repulsive interaction between species in the classical single-range model, a competing mechanism between the short-range attractive and mid-range repulsive interactions is imposed within each species. Besides extending the range of attainable surface tension as compared with the single-range model, the proposed scheme is also shown to achieve a positive disjoining pressure, independently of the viscosity ratio. The latter property is crucial for many microfluidic applications involving a collection of disperse droplets with a different viscosity from the continuum phase. As a preliminary application, the relative effective viscosity of a pressure-driven emulsion in a planar channel is computed.
Soft particles at fluid interfaces play an important role in many aspects of our daily life, such as the food industry, paints and coatings, and medical applications. Analytical methods are not capable of describing the emergent effects of the complex dynamics of suspensions of many soft particles, whereas experiments typically either only capture bulk properties or require invasive methods. Computational methods are therefore a great tool to complement experimental work. However, an efficient and versatile numerical method is needed to model dense suspensions of many soft particles. In this article we propose a method to simulate soft particles in a multi-component fluid, both at and near fluid-fluid interfaces, based on the lattice Boltzmann method, and characterize the error stemming from the fluid-structure coupling for the particle equilibrium shape when adsorbed onto a fluid-fluid interface. Furthermore, we characterize the influence of the preferential contact angle of the particle surface and the particle softness on the vertical displacement of the center of mass relative to the fluid interface. Finally, we demonstrate the capability of our model by simulating a soft capsule adsorbing onto a fluid-fluid interface with a shear flow parallel to the interface, and the covering of a droplet suspended in another fluid by soft particles with different wettability.
We study a mesoscopic model for the flow of amorphous solids. The model is based on the key features identified at the microscopic level, namely peri- ods of elastic deformation interspersed with localised rearrangements of parti- cles that induce long-range elastic deformation. These long-range deformations are derived following a continuum mechanics approach, in the presence of solid boundaries, and are included in full in the model. Indeed, they mediate spatial cooperativity in the flow, whereby a localised rearrangement may lead a distant region to yield. In particular, we simulate a channel flow and find manifestations of spatial cooperativity that are consistent with published experimental obser- vations for concentrated emulsions in microchannels. Two categories of effects are distinguished. On the one hand, the coupling of regions subject to different shear rates, for instance,leads to finite shear rate fluctuations in the seemingly un- sheared plug in the centre of the channel. On the other hand, there is convinc- ing experimental evidence of a specific rheology near rough walls. We discuss diverse possible physical origins for this effect, and we suggest that it may be associated with the bumps of particles into surface asperities as they slide along the wall.
We implement the statistically sound G-JF thermostat for Langevin Dynamics simulations into the ESPREesSo molecular package for large-scale simulations of soft matter systems. The implemented integration method is tested against the integrator currently used by the molecular package in simulations of a fluid bilayer membrane. While the latter exhibits deviations in the sampling statistics that increase with the integration time step dt, the former reproduces near-correct configurational statistics for all dt within the stability range of the simulations. We conclude that, with very modest revisions to existing codes, one can significantly improve the performance of statistical sampling using Langevin thermostats.
We develop a modified Poisson-Nernst-Planck model which includes both the long-range Coulomb and short-range hard-sphere correlations in its free energy functional such that the model can accurately describe the ion transport in complex environment and under a nanoscale confinement. The Coulomb correlation including the dielectric polarization is treated by solving a generalized Debye-Huckel equation which is a Greens function equation with the correlation energy of a test ion described by the self Greens function. The hard-sphere correlation is modeled through the modified fundamental measure theory. The resulting model is available for problems beyond the mean-field theory such as problems with variable dielectric media, multivalent ions, and strong surface charge density. We solve the generalized Debye-Huckel equation by the Wentzel-Kramers-Brillouin approximation, and study the electrolytes between two parallel dielectric surfaces. In comparison to other modified models, the new model is shown more accurate in agreement with particle-based simulations and capturing the physical properties of ionic structures near interfaces.
Stochastic dynamical systems with continuous symmetries arise commonly in nature and often give rise to coherent spatio-temporal patterns. However, because of their random locations, these patterns are not well captured by current order reduction techniques and a large number of modes is typically necessary for an accurate solution. In this work, we introduce a new methodology for efficient order reduction of such systems by combining (i) the method of slices, a symmetry reduction tool, with (ii) any standard order reduction technique, resulting in efficient mixed symmetry-dimensionality reduction schemes. In particular, using the Dynamically Orthogonal (DO) equations in the second step, we obtain a novel nonlinear Symmetry-reduced Dynamically Orthogonal (SDO) scheme. We demonstrate the performance of the SDO scheme on stochastic solutions of the 1D Korteweg-de Vries and 2D Navier-Stokes equations.