Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical simulations of emulsions in shear flows
185
0
0.0
(
0
)
Added by
Marco Edoardo Rosti
Publication date
2018
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present a modification of a recently developed volume of fluid method for multiphase problems, so that it can be used in conjunction with a fractional step-method and fast Poisson solver, and validate it with standard benchmark problems. We then consider emulsions of two-fluid systems and study their rheology in a plane Couette flow in the limit of vanishing inertia. We examine the dependency of the effective viscosity on the volume-fraction (from 10% to 30%) and the Capillary number (from 0.1 to 0.4) for the case of density and viscosity ratio 1. We show that the effective viscosity decreases with the deformation and the applied shear (shear-thinning) while exhibits a non-monotonic behavior with respect to the volume fraction. We report the appearance of a maximum in the effective viscosity curve and compare the results with those of suspensions of rigid and deformable particles and capsules. We show that the flow in the solvent is mostly a shear flow, while it is mostly rotational in the suspended phase; moreover this behavior tends to reverse as the volume fraction increases. Finally, we evaluate the contributions to the total shear stress of the viscous stresses in the two fluids and of the interfacial force between them.
rate research
Read More
Multiphase shear flows often show banded structures that affect the global behavior of complex fluids e.g. in microdevices. Here we investigate numerically the banding of emulsions, i.e. the formation of regions of high and low volume fraction, alternated in the vorticity direction and aligned with the flow (shear bands). These bands are associated with a decrease of the effective viscosity of the system. To understand the mechanism of banding experimentally observed we have performed interface resolved simulations of the two-fluid system. The experiments were perfomed starting with a random distribution of droplets which, under the applied shear, evolves in time resulting in a phase separation. To numerically reproduce this process, the banded structures are initialized in a narrow channel confined by two walls moving in opposite direction. We find that the initial banded distribution is stable when droplets are free to merge and unstable when coalescence is prevented. In this case, additionally, the effective viscosity of the system increases, resembling the rheological behavior of suspensions of deformable particles. Droplets coalescence, on the other hand, allows emulsions to reduce the total surface of the system and hence the energy dissipation associated to the deformation, which in turn reduces the effective viscosity.
We perform $3$D numerical simulations to investigate the sedimentation of a single sphere in the absence and presence of a simple cross shear flow in a yield stress fluid with weak inertia. In our simulations, the settling flow is considered to be the primary flow, whereas the linear cross shear flow is a secondary flow with amplitude $10%$ of the primary flow. To study the effects of elasticity and plasticity of the carrying fluid on the sphere drag as well as the flow dynamics, the fluid is modeled using the elastovisco-plastic (EVP) constitutive laws proposed by cite{saramito2009new}. The extra non-Newtonian stress tensor is fully coupled with the flow equation and the solid particle is represented by an immersed boundary (IB) method. Our results show that the fore-aft asymmetry in the velocity is less pronounced and the negative wake disappears when a linear cross shear flow is applied. We find that the drag on a sphere settling in a sheared yield stress fluid is reduced significantly as compared to an otherwise quiescent fluid. More importantly, the sphere drag in the presence of a secondary cross shear flow cannot be derived from the pure sedimentation drag law owing to the non-linear coupling between the simple shear flow and the uniform flow. Finally, we show that the drag on the sphere settling in a sheared yield-stress fluid is reduced at higher material elasticity mainly due to the form and viscous drag reduction.
We simulated two particle-based fluid models, namely multiparticle collision dynamics and dissipative particle dynamics, under shear using reverse nonequilibrium simulations (RNES). In cubic periodic simulation boxes, the expected shear flow profile for a Newtonian fluid developed, consistent with the fluid viscosities. However, unexpected secondary flows along the shear gradient formed when the simulation box was elongated in the flow direction. The standard shear flow profile was obtained when the simulation box was longer in the shear-gradient dimension than the flow dimension, while the secondary flows were always present when the flow dimension was at least 25% larger than the shear-gradient dimension. The secondary flows satisfy the boundary conditions imposed by the RNES and have a lower rate of viscous dissipation in the fluid than the corresponding unidirectional flows. This work highlights a previously unappreciated limitation of RNES for generating shear flow in simulation boxes that are elongated in the flow dimension, an important consideration when applying RNES to complex fluids like polymer solutions.
To investigate the formation mechanism of energy spectra of internal waves in the oceans, direct numerical simulations are performed. The simulations are based on the reduced dynamical equations of rotating stratified turbulence. In the reduced dynamical equations only wave modes are retained, and vortices and horizontally uniform vertical shears are excluded. Despite the simplifications, our simulations reproduce some key features of oceanic internal-wave spectra: accumulation of energy at near-inertial waves and realistic frequency and horizontal wavenumber dependencies. Furthermore, we provide evidence that formation of the energy spectra in the inertial subrange is dominated by scale-separated interactions with the near-inertial waves. These findings support oceanographers intuition that spectral energy density of internal waves is the result of predominantly wave-wave interactions.
A kinematic approach for the identification of flow instabilities is proposed. By defining a flow instability in the Lagrangian frame as the increased folding of lines of fluid particles, subtle perturbations and unstable growth thereof are detected early based solely on the curvature change of material lines over finite time. The material line curvature is objective, parametrization independent, and can be applied to flows of general complexity without knowledge of the base flow. An analytic connection between the growth of Eulerian velocity modes perturbing a general shear flow and the induced flow map and Lagrangian curvature change is derived. The approach is verified to capture instabilities promptly in a temporally developing jet flow, an unstable separated shear flow over a cambered airfoil, and in the onset of a wake instability behind a circular cylinder.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
