Numerical study of electric field effects on the deformation of two-dimensional liquid drops in simple shear flow at arbitrary Reynolds number

2009 ◽  
Vol 626 ◽  
pp. 367-393 ◽  
Author(s):  
STEFAN MÄHLMANN ◽  
DEMETRIOS T. PAPAGEORGIOU

The effect of an electric field on a periodic array of two-dimensional liquid drops suspended in simple shear flow is studied numerically. The shear is produced by moving the parallel walls of the channel containing the fluids at equal speeds but in opposite directions and an electric field is generated by imposing a constant voltage difference across the channel walls. The level set method is adapted to electrohydrodynamics problems that include a background flow in order to compute the effects of permittivity and conductivity differences between the two phases on the dynamics and drop configurations. The electric field introduces additional interfacial stresses at the drop interface and we perform extensive computations to assess the combined effects of electric fields, surface tension and inertia. Our computations for perfect dielectric systems indicate that the electric field increases the drop deformation to generate elongated drops at steady state, and at the same time alters the drop orientation by increasing alignment with the vertical, which is the direction of the underlying electric field. These phenomena are observed for a range of values of Reynolds and capillary numbers. Computations using the leaky dielectric model also indicate that for certain combinations of electric properties the drop can undergo enhanced alignment with the vertical or the horizontal, as compared to perfect dielectric systems. For cases of enhanced elongation and alignment with the vertical, the flow positions the droplets closer to the channel walls where they cause larger wall shear stresses. We also establish that a sufficiently strong electric field can be used to destabilize the flow in the sense that steady-state droplets that can exist in its absence for a set of physical parameters, become increasingly and indefinitely elongated until additional mechanisms can lead to rupture. It is suggested that electric fields can be used to enhance such phenomena.

Author(s):  
Rémi Vachon ◽  
Mohsen Bazargan ◽  
Christoph F Hieronymus ◽  
Erika Ronchin ◽  
Bjarne Almqvist

Summary Elongate inclusions immersed in a viscous fluid generally rotate at a rate that is different from the local angular velocity of the flow. Often, a net alignment of the inclusions develops, and the resulting shape preferred orientation (SPO) of the particle ensemble can then be used as a strain marker that allows reconstruction of the fluid’s velocity field. Much of the previous work on the dynamics of flow-induced particle rotations has focused on spatially homogeneous flows with large-scale tectonic deformations as the main application. Recently, the theory has been extended to spatially varying flows, such as magma with embedded crystals moving through a volcanic plumbing system. Additionally, an evolution equation has been introduced for the probability density function (PDF) of crystal orientations. Here, we apply this new theory to a number of simple, two-dimensional flow geometries commonly encountered in magmatic intrusions, such as flow from a dyke into a reservoir or from a reservoir into a dyke, flow inside an inflating or deflating reservoir, flow in a dyke with a sharp bend, and thermal convection in a magma chamber. The main purpose is to provide a guide for interpreting field observations and for setting up more complex flow models with embedded crystals. As a general rule, we find that a larger aspect ratio of the embedded crystals causes a more coherent alignment of the crystals, while it has only a minor effect on the geometry of the alignment pattern. Due to various perturbations in the crystal rotation equations that are expected in natural systems, we show that the time-periodic behavior found in idealized systems is probably short-lived in nature, and the crystal alignment is well described by the time-averaged solution. We also confirm some earlier findings. For example, near channel walls, fluid flow often follows the bounding surface and the resulting simple shear flow causes preferred crystal orientations that are approximately parallel to the boundary. Where pure shear deformation dominates, there is a tendency for crystals to orient themselves in the direction of the greatest tensile strain rate. Where flow impinges on a boundary, for example in an inflating magma chamber or as part of a thermal convection pattern, the stretching component of pure shear aligns with the boundary, and the crystals orient themselves in that direction. In the field, this local pattern may be difficult to distinguish from a boundary-parallel simple shear flow. Pure shear also dominates along the walls of a deflating magma chamber and in places where the flow turns away from the reservoir walls, but in these locations, the preferred crystal orientation is perpendicular to the wall. Overall, we find that our calculated patterns of crystal orientations agree well with results from analogue experiments where similar geometries are available.


1974 ◽  
Vol 65 (2) ◽  
pp. 365-400 ◽  
Author(s):  
B. P. Ho ◽  
L. G. Leal

The familiar Segré-Silberberg effect of inertia-induced lateral migration of a neutrally buoyant rigid sphere in a Newtonian fluid is studied theoretically for simple shear flow and for two-dimensional Poiseuille flow. It is shown that the spheres reach a stable lateral equilibrium position independent of the initial position of release. For simple shear flow, this position is midway between the walls, whereas for Poiseuille flow, it is 0·6 of the channel half-width from the centre-line. Particle trajectories are calculated in both cases and compared with available experimental data. Implications for the measurement of the rheological properties of a dilute suspension of spheres are discussed.


1998 ◽  
Vol 361 ◽  
pp. 117-143 ◽  
Author(s):  
S. RAMANUJAN ◽  
C. POZRIKIDIS

The deformation of a liquid capsule enclosed by an elastic membrane in an infinite simple shear flow is studied numerically at vanishing Reynolds numbers using a boundary-element method. The surface of the capsule is discretized into quadratic triangular elements that form an evolving unstructured grid. The elastic membrane tensions are expressed in terms of the surface deformation gradient, which is evaluated from the position of the grid points. Compared to an earlier formulation that uses global curvilinear coordinates, the triangular-element formulation suppresses numerical instabilities due to uneven discretization and thus enables the study of large deformations and the investigation of the effect of fluid viscosities. Computations are performed for capsules with spherical, spheroidal, and discoidal unstressed shapes over an extended range of the dimensionless shear rate and for a broad range of the ratio of the internal to surrounding fluid viscosities. Results for small deformations of spherical capsules are in quantitative agreement with the predictions of perturbation theories. Results for large deformations of spherical capsules and deformations of non-spherical capsules are in qualitative agreement with experimental observations of synthetic capsules and red blood cells. We find that initially spherical capsules deform into steady elongated shapes whose aspect ratios increase with the magnitude of the shear rate. A critical shear rate above which capsules exhibit continuous elongation is not observed for any value of the viscosity ratio. This behaviour contrasts with that of liquid drops with uniform surface tension and with that of axisymmetric capsules subject to a stagnation-point flow. When the shear rate is sufficiently high and the viscosity ratio is sufficiently low, liquid drops exhibit continuous elongation leading to breakup. Axisymmetric capsules deform into thinning needles at sufficiently high rates of elongation, independent of the fluid viscosities. In the case of capsules in shear flow, large elastic tensions develop at large deformations and prevent continued elongation, stressing the importance of the vorticity of the incident flow. The long-time behaviour of deformed capsules depends strongly on the unstressed shape. Oblate capsules exhibit unsteady motions including oscillation about a mean configuration at low viscosity ratios and continuous rotation accompanied by periodic deformation at high viscosity ratios. The viscosity ratio at which the transition from oscillations to tumbling occurs decreases with the sphericity of the unstressed shape. Results on the effective rheological properties of dilute suspensions confirm a non-Newtonian shear-thinning behaviour.


1996 ◽  
Vol 320 (-1) ◽  
pp. 395 ◽  
Author(s):  
Xiaofan Li ◽  
R. Charles ◽  
C. Pozrikidis

1994 ◽  
Vol 23 (2) ◽  
pp. 251-278 ◽  
Author(s):  
M.R. Kennedy ◽  
C. Pozrikidis ◽  
R. Skalak

2016 ◽  
Vol 808 ◽  
Author(s):  
Chenggong Li ◽  
Mao Ye ◽  
Zhongmin Liu

We investigate numerically the rotational behaviour of a circular porous particle suspended in a two-dimensional (2D) simple shear flow with fluid inertia at particle shear Reynolds number up to 108. We use the volume-averaged macroscopic momentum equation to formulate the flow field inside and outside the moving porous particle, which is solved by a modified single relaxation time lattice Boltzmann method. The effects of fluid inertia, confinement of the bounding walls, and permeability of the particle are studied. Our two-dimensional simulation results confirm that the permeability has little effect on the rotation of a porous particle in unbounded shear flow without fluid inertia (Masoud, Stone & Shelley, J. Fluid Mech., vol. 733, 2013, R6), but also suggest that the role of permeability cannot be neglected when the confinement effect is significant, or the fluid inertia is not negligible. The fluid inertia and the confined walls have similar effects on the rotation of a porous particle as that on a solid impermeable particle. The angular velocity decays with an increase in fluid inertia, and the confinement effect suppresses the angular velocity to a shear rate ratio below 0.5. A simple scaling argument based on the balance of torque exerted by fluid flows adjacent to the two bounding walls and that due to the flow recirculation can explain our results.


Sign in / Sign up

Export Citation Format

Share Document