Numerical Simulation of Two-Dimensional Drops Suspended in Simple Shear Flow at Nonzero Reynolds Numbers

2011 ◽  
Vol 133 (3) ◽  
Author(s):  
S. Mortazavi ◽  
Y. Afshar ◽  
H. Abbaspour

The motion of deformable drops suspended in a linear shear flow at nonzero Reynolds numbers is studied by numerical simulations in two dimensions. It is found that a deformable drop migrates toward the center of the channel in agreement with experimental findings at small Reynolds numbers. However, at relatively high Reynolds numbers (Re=80) and small deformation, the drop migrates to an equilibrium position off the centerline. Suspension of drops at a moderate areal fraction (φ=0.44) is studied by simulations of 36 drops. The flow is studied as a function of the Reynolds number and a shear thinning behavior is observed. The results for the normal stress difference show oscillations around a mean value at small Reynolds numbers, and it increases as the Reynolds number is raised. Simulations of drops at high areal fraction (φ=0.66) show that if the Capillary number is kept constant, the effective viscosity does not change in the range of considered Reynolds numbers (0.8–80). The normal stress difference is also a weak function of the Reynolds number. It is also found that similar to flows of granular materials, suspension of drops at finite Reynolds numbers shows the same trend for the density and fluctuation energy distribution across the channel.

2011 ◽  
Vol 266 ◽  
pp. 130-134
Author(s):  
Jin Yan Wang ◽  
Jing Bo Chen ◽  
Chang Yu Shen

The paper presents a numerical simulation for the isothermal flow-induced crystallization of polyethylene under a simple shear flow. The effect of flow on crystllization is considered through the simple mathematical relationship between the additional number of nuclei induced by shear treatment and the first normal stress difference. Leonov viscoelastic model and Avrami model are used to describe the normal stress difference and the crystallization kinetics, respectively. It is found that the short-term shear treatment has a large effect on the crystallization dynamics of polyethylene , but the effect of the intensity of the shear flow is not infinite ,which shows a saturation phenomenon, namely, the accelerated degree of crystallization tending to level off when the shear rate or shear time is large enough.


2015 ◽  
Vol 786 ◽  
pp. 110-127 ◽  
Author(s):  
D. Matsunaga ◽  
Y. Imai ◽  
T. Yamaguchi ◽  
T. Ishikawa

We present a numerical analysis of the rheology of a dense suspension of spherical capsules in simple shear flow in the Stokes flow regime. The behaviour of neo-Hookean capsules is simulated for a volume fraction up to${\it\phi}=0.4$by graphics processing unit computing based on the boundary element method with a multipole expansion. To describe the specific viscosity using a polynomial equation of the volume fraction, the coefficients of the equation are calculated by least-squares fitting. The results suggest that the effect of higher-order terms is much smaller for capsule suspensions than rigid sphere suspensions; for example,$O({\it\phi}^{3})$terms account for only 8 % of the specific viscosity even at${\it\phi}=0.4$for capillary numbers$Ca\geqslant 0.1$. We also investigate the relationship between the deformation and orientation of the capsules and the suspension rheology. When the volume fraction increases, the deformation of the capsules increases while the orientation angle of the capsules with respect to the flow direction decreases. Therefore, both the specific viscosity and the normal stress difference increase with volume fraction due to the increased deformation, whereas the decreased orientation angle suppresses the specific viscosity, but amplifies the normal stress difference.


Soft Matter ◽  
2020 ◽  
Vol 16 (43) ◽  
pp. 9864-9875
Author(s):  
Young Ki Lee ◽  
Kyu Hyun ◽  
Kyung Hyun Ahn

The first normal stress difference (N1) as well as shear stress of non-Brownian hard-sphere suspensions in small to large amplitude oscillatory shear flow is investigated.


1999 ◽  
Vol 17 (No. 1) ◽  
pp. 23-30 ◽  
Author(s):  
P. Novotna ◽  
M. Houska ◽  
V. Sopr ◽  
H. Valentova ◽  
P. Stern

The shear flow rheological properties of sugar solutions (70% w/w concentration) modified by different cellulose derivatives have been measured. Thickeners  were expected to cause the viscoelastic behaviour of the resulting sol ution. Therefore, the elastic rheological parameters were measured by oscillatory shear technique (phase angle, elastic modulus) and also the first normal stress difference N<sub>1</sub>. The increase of thickener concen tration caused a moderate increase of non-Newtonian behaviour in the shear flow. The sensory viscosity (ra nged between 0 and 100%) was evaluated by five different methods - as an effort for stirring with teaspoon, time for flowing down the spoon, slurping from spoon, compression between tongue and palate and swallowing. The influence of shear viscosity and first normal difference on sensory viscosity was tested. Correlation procedu re between change of sensory viscosity .tlSE and change of shear viscosity .tlJ.Iz showed that only for swallowing there is a statistically evident de­pendence. The correlation between change of sensory viscosity t.SE and first normal stress difference N<sub>1</sub> is not statistically   evident. For all the methods of sensory evaluation the dependence between these parameters is only weak and indirect (with increasing normal stress difference the sensory viscosity is decreasing).


Author(s):  
Orest Shardt ◽  
J. J. Derksen ◽  
Sushanta K. Mitra

When droplets collide in a shear flow, they may coalesce or remain separate after the collision. At low Reynolds numbers, droplets coalesce when the capillary number does not exceed a critical value. We present three-dimensional simulations of droplet coalescence in a simple shear flow. We use a free-energy lattice Boltzmann method (LBM) and study the collision outcome as a function of the Reynolds and capillary numbers. We study the Reynolds number range from 0.2 to 1.4 and capillary numbers between 0.1 and 0.5. We determine the critical capillary number for the simulations (0.19) and find that it is does not depend on the Reynolds number. The simulations are compared with experiments on collisions between confined droplets in shear flow. The critical capillary number in the simulations is about a factor of 25 higher than the experimental value.


2002 ◽  
Vol 455 ◽  
pp. 21-61 ◽  
Author(s):  
ALEXANDER Z. ZINCHENKO ◽  
ROBERT H. DAVIS

An efficient algorithm for hydrodynamical interaction of many deformable drops subject to shear flow at small Reynolds numbers with triply periodic boundaries is developed. The algorithm, at each time step, is a hybrid of boundary-integral and economical multipole techniques, and scales practically linearly with the number of drops N in the range N < 1000, for NΔ ∼ 103 boundary elements per drop. A new near-singularity subtraction in the double layer overcomes the divergence of velocity iterations at high drop volume fractions c and substantial viscosity ratio γ. Extensive long-time simulations for N = 100–200 and NΔ = 1000–2000 are performed up to c = 0.55 and drop-to-medium viscosity ratios up to λ = 5, to calculate the non-dimensional emulsion viscosity μ* = Σ12/(μeγ˙), and the first N1 = (Σ11−Σ22)/(μe[mid ]γ˙[mid ]) and second N2 = (Σ22−Σ33)/(μe[mid ]γ˙[mid ]) normal stress differences, where γ˙ is the shear rate, μe is the matrix viscosity, and Σij is the average stress tensor. For c = 0.45 and 0.5, μ* is a strong function of the capillary number Ca = μe[mid ]γ˙[mid ]a/σ (where a is the non-deformed drop radius, and σ is the interfacial tension) for Ca [Lt ] 1, so that most of the shear thinning occurs for nearly non-deformed drops. For c = 0.55 and λ = 1, however, the results suggest phase transition to a partially ordered state at Ca [les ] 0.05, and μ* becomes a weaker function of c and Ca; using λ = 3 delays phase transition to smaller Ca. A positive first normal stress difference, N1, is a strong function of Ca; the second normal stress difference, N2, is always negative and is a relatively weak function of Ca. It is found at c = 0.5 that small systems (N ∼ 10) fail to predict the correct behaviour of the viscosity and can give particularly large errors for N1, while larger systems N [ges ] O(102)show very good convergence. For N ∼ 102 and NΔ ∼ 103, the present algorithm is two orders of magnitude faster than a standard boundary-integral code, which has made the calculations feasible.


2010 ◽  
Vol 646 ◽  
pp. 255-296 ◽  
Author(s):  
R. VIVEK RAJA ◽  
GANESH SUBRAMANIAN ◽  
DONALD L. KOCH

The behaviour of an isolated nearly spherical drop in an ambient linear flow is examined analytically at small but finite Reynolds numbers, and thereby the first effects of inertia on the bulk stress in a dilute emulsion of neutrally buoyant drops are calculated. The Reynolds numbers, Re = a2ρ/μ and $\hat{\Rey} \,{=}\, \dot{\gamma}a^2\rho/\hat{\mu}$, are the relevant dimensionless measures of inertia in the continuous and disperse(drop) phases, respectively. Here, a is the drop radius, is the shear rate, ρ is the common density and and μ are, respectively, the viscosities of the drop and the suspending fluid. The assumption of nearly spherical drops implies the dominance of surface tension, and the analysis therefore corresponds to the limit of the capillary number(Ca) based on the viscosity of the suspending fluid being small but finite; in other words, Ca ≪ 1, where Ca = μa/T, T being the coefficient of interfacial tension. The bulk stress is determined to O(φRe) via two approaches. The first one is the familiar direct approach based on determining the force density associated with the disturbance velocity field on the surface of the drop; the latter is determined to O(Re) from a regular perturbation analysis. The second approach is based on a novel reciprocal theorem formulation and allows the calculation, to O(Re), of the drop stresslet, and hence the emulsion bulk stress, with knowledge of only the leading-order Stokes fields. The first approach is used to determine the bulk stress for linear flows without vortex stretching, while the reciprocal theorem approach allows one to generalize this result to any linear flow. For the case of simple shear flow, the inertial contributions to the bulk stress lead to normal stress differences(N1, N2) at O(φRe), where φ(≪1) is the volume fraction of the disperse phase. Inertia leads to negative and positive contributions, respectively, to N1 and N2 at O(φRe). The signs of the inertial contributions to the normal stress differences may be related to the O(ReCa) tilting of the drop towards the velocity gradient direction. These signs are, however, opposite to that of the normal stress differences in the creeping flow limit. The latter are O(φCa) and result from an O(Ca2) deformation of the drop acting to tilt it towards the flow axis. As a result, even a modest amount of inertia has a significant effect on the rheology of a dilute emulsion. In particular, both normal stress differences reverse sign at critical Reynolds numbers(Rec) of O(Ca) in the limit Ca ≪ 1. This criterion for the reversal in the signs of N1 and N2 is more conveniently expressed in terms of a critical Ohnesorge number(Oh) based on the viscosity of the suspending fluid, where Oh = μ/(ρaT)1/2. The critical Ohnesorge number for a sign reversal in N1 is found to be lower than that for N2, and the precise numerical value is a function of λ. In uniaxial extensional flow, the Trouton viscosity remains unaltered at O(φRe), the first effects of inertia now being restricted to O(φRe3/2). The analytical results for simple shear flow compare favourably with the recent numerical simulations of Li & Sarkar (J. Rheol., vol. 49, 2005, p. 1377).


1995 ◽  
Vol 283 ◽  
pp. 273-285 ◽  
Author(s):  
H. Nirschl ◽  
H. A. Dwyer ◽  
V. Denk

Three-dimensional solutions have been obtained for the steady simple shear flow over a spherical particle in the intermediate Reynolds number range 0.1 [les ] Re [les ] 100. The shear flow was generated by two walls which move at the same speed but in opposite directions, and the particle was located in the middle of the gap between the walls. The particle-wall interaction is treated by introducing a fully three-dimensional Chimera or overset grid scheme. The Chimera grid scheme allows each component of a flow to be accurately and efficiently treated. For low Reynolds numbers and without any wall influence we have verified the solution of Taylor (1932) for the shear around a rigid sphere. With increasing Reynolds numbers the angular velocity for zero moment for the sphere decreases with increasing Reynolds number. The influence of the wall has been quantified with the global particle surface characteristics such as net torque and Nusselt number. A detailed analysis of the influence of the wall distance and Reynolds number on the surface distributions of pressure, shear stress and heat transfer has also been carried out.


Sign in / Sign up

Export Citation Format

Share Document