Streaming-potential phenomena in the thin-Debye-layer limit. Part 2. Moderate Péclet numbers

2012 ◽  
Vol 704 ◽  
pp. 109-136 ◽  
Author(s):  
Ory Schnitzer ◽  
Itzchak Frankel ◽  
Ehud Yariv

AbstractMacroscale description of streaming-potential phenomena in the thin-double-layer limit, and in particular the associated electro-viscous forces, has been a matter of long-standing controversy. In part 1 of this work (Yariv, Schnitzer & Frankel, J. Fluid Mech., vol. 685, 2011, pp. 306–334) we identified that the product of the Hartmann ($\lambda $) and Péclet ($\mathit{Pe}$) numbers is $O({\delta }^{\ensuremath{-} 2} )$, $\delta $ being the dimensionless Debye thickness. This scaling relationship defines a one-family class of limit processes appropriate to the consistent analysis of this singular problem. In that earlier contribution we focused on the generic problems associated with moderate $\lambda $ and large $\mathit{Pe}$, where the streaming-potential magnitude is comparable to the thermal voltage. Here we consider the companion generic limit of moderate Péclet numbers and large Hartmann numbers, deriving the appropriate macroscale model wherein the Debye-layer physics is represented by effective boundary conditions. Since the induced electric field is asymptotically smaller, calculation of these conditions requires higher asymptotic orders in analysing the Debye-scale transport. Nonetheless, the leading-order electro-viscous forces are of the same $O({\delta }^{2} )$ relative magnitude as those previously obtained in the large-$\mathit{Pe}$ limit. The structure of these forces is different, however, first because the small Maxwell stresses do not contribute at leading order, and second because salt polarization results in a dominant diffuso-osmotic slip. Since the salt distribution is governed by an advection–diffusion equation, this slip gives rise to electro-viscous forces which are nonlinear in the driving flow. The resulting scheme is illustrated by the calculation of the electro-viscous excess drag in the prototype problem of a translating sphere.

2011 ◽  
Vol 685 ◽  
pp. 306-334 ◽  
Author(s):  
Ehud Yariv ◽  
Ory Schnitzer ◽  
Itzchak Frankel

AbstractElectrokinetic streaming-potential phenomena are driven by imposed relative motion between liquid electrolytes and charged solids. Owing to non-uniform convective ‘surface’ current within the Debye layer Ohmic currents from the electro-neutral bulk are required to ensure charge conservation thereby inducing a bulk electric field. This, in turn, results in electro-viscous drag enhancement. The appropriate modelling of these phenomena in the limit of thin Debye layers $\delta \ensuremath{\rightarrow} 0$ ($\delta $ denoting the dimensionless Debye thickness) has been a matter of ongoing controversy apparently settled by Cox’s seminal analysis (J. Fluid Mech., vol. 338, 1997, p. 1). This analysis predicts electro-viscous forces that scale as ${\delta }^{4} $ resulting from the perturbation of the original Stokes flow with the Maxwell-stress contribution only appearing at higher orders. Using scaling analysis we clarify the distinction between the normalizations pertinent to field- and motion-driven electrokinetic phenomena, respectively. In the latter class we demonstrate that the product of the Hartmann & Péclet numbers is $O({\delta }^{\ensuremath{-} 2} )$ contrary to Cox (1997) where both parameters are assumed $O(1)$. We focus on the case where motion-induced fields are comparable to the thermal scale and accordingly present a singular-perturbation analysis for the limit where the Hartmann number is $O(1)$ and the Péclet number is $O({\delta }^{\ensuremath{-} 2} )$. Electric-current matching between the Debye layer and the electro-neutral bulk provides an inhomogeneous Neumann condition governing the electric field in the latter. This field, in turn, results in a velocity perturbation generated by a Smoluchowski-type slip condition. Owing to the dominant convection, the present analysis yields an asymptotic structure considerably simpler than that of Cox (1997): the electro-viscous effect now already appears at $O({\delta }^{2} )$ and is contributed by both Maxwell and viscous stresses. The present paradigm is illustrated for the prototypic problem of a sphere sedimenting in an unbounded fluid domain with the resulting drag correction differing from that calculated by Cox (1997). Independently of current matching, salt-flux matching between the Debye layer and the bulk domain needs also to be satisfied. This subtle point has apparently gone unnoticed in the literature, perhaps because it is trivially satisfied in field-driven problems. In the present limit this requirement seems incompatible with the uniform salt distribution in the convection-dominated bulk domain. This paradox is resolved by identifying the dual singularity associated with the limit $\delta \ensuremath{\rightarrow} 0$ in motion-driven problems resulting in a diffusive layer of $O({\delta }^{2/ 3} )$ thickness beyond the familiar $O(\delta )$-wide Debye layer.


2015 ◽  
Vol 786 ◽  
pp. 84-109 ◽  
Author(s):  
Ory Schnitzer ◽  
Ehud Yariv

We employ the moderate-Péclet-number macroscale model developed in part 2 of this sequence (Schnitzer et al., J. Fluid Mech., vol. 704, 2012, pp. 109–136) towards the calculation of electroviscous forces on charged solid particles engendered by an imposed relative motion between these particles and the electrolyte solution in which they are suspended. In particular, we are interested in the kinematic irreversibility of these forces, stemming from the diffusio-osmotic slip which accompanies the salt-concentration polarisation induced by that imposed motion. We illustrate the electroviscous irreversibility using two prototypic problems, one involving side-by-side sedimentation of two spherical particles, and the other involving a force-free spherical particle suspended in the vicinity of a planar wall and exposed to a simple shear flow. We focus on the pertinent limit of near-contact configurations, where use of lubrication approximations provides closed-form expressions for the leading-order lateral repulsion. In this approximation scheme, the need to solve the advection–diffusion equation governing the salt-concentration polarisation is circumvented.


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 215-229 ◽  
Author(s):  
Lichi Deng ◽  
Michael J. King

Summary Spontaneous and forced imbibition are recognized as important recovery mechanisms in naturally fractured reservoirs because the capillary force controls the movement of the fluid between the matrix and the fracture. For unconventional reservoirs, imbibition is also important because the capillary pressure is more dominant in these tighter formations, and a theoretical understanding of the flow mechanism for the imbibition process will benefit the understanding of important multiphase-flow phenomena such as waterblocking. In this paper, a new semianalytic method is presented to examine the interaction between spontaneous and forced imbibition and to quantitatively represent the transient imbibition process. The methodology solves the partial-differential equation (PDE) of unsteady-state immiscible, incompressible flow with arbitrary saturation-dependent functions using the normalized water flux concept, which is identical to the fractional-flow terminology used in the traditional Buckley-Leverett analysis. The result gives a universal inherent relationship between time, normalized water flux, saturation profile, and the ratio between cocurrent and total flux. The current analysis also develops a novel stability envelope outside of which the flow becomes unstable caused by strong capillary forces, and the characteristic dimensionless parameter shown in the envelope is derived from the intrinsic properties of the rock and fluid system, and it can describe the relative magnitude of capillary and viscous forces at the continuum scale. This dimensionless parameter is consistently applicable in both capillary-dominated and viscous-dominated flow conditions.


1976 ◽  
Vol 78 (3) ◽  
pp. 489-512 ◽  
Author(s):  
Jorg Imberger ◽  
Rory Thompson ◽  
Chris Fandry

The time development of two-dimensional fluid motion induced by a line sink in a rectangular, density stratified reservoir with a free surface is given. It is shown that the initiation of such a sink gives birth to a spectrum of internal expanding shear fronts with a progressively decreasing vertical wavelength. These fronts move out from the sink and travel towards the far wall, where they are reflected. This process ceases once the front with a vertical wavelength equal to the steady withdrawal-layer thickness has reached the end wall. The fronts so introduced continue to move back and forth, expanding to standing waves if the viscosity of the fluid is small enough. The evolution and nature of the withdrawal layer are shown to depend critically on the relative magnitude of the convective inertia and viscous forces, the number of reflexions from the rear wall and the Prandtl number.


2018 ◽  
Vol 841 ◽  
pp. 203-227 ◽  
Author(s):  
A. L. Sánchez ◽  
C. Martínez-Bazán ◽  
C. Gutiérrez-Montes ◽  
E. Criado-Hidalgo ◽  
G. Pawlak ◽  
...  

Radionuclide scanning images published inNatureby Di Chiro in 1964 showed a downward migration along the spinal canal of particle tracers injected in the brain ventricles while also showing an upward flow of tracers injected in the lumbar region of the canal. These observations, since then corroborated by many radiological measurements, have been the basis for the hypothesis that there must be an active circulation mechanism associated with the transport of cerebrospinal fluid (CSF) deep down into the spinal canal and subsequently returning a portion back to the cranial vault. However, to date, there has been no physical explanation for the mechanism responsible for the establishment of such a bulk recirculating motion. To investigate the origin and characteristics of this recirculating flow, we have analyzed the motion of the CSF in the subarachnoid space of the spinal canal. Our analysis accounts for the slender geometry of the spinal canal, the small compliance of the dura membrane enclosing the CSF in the canal, and the fact that the CSF is confined to a thin annular subarachnoid space surrounding the spinal cord. We apply this general formulation to study the characteristics of the flow generated in a simplified model of the spinal canal consisting of a slender compliant cylindrical pipe with a coaxial cylindrical inclusion, closed at its distal end, and subjected to small periodic pressure pulsations at its open entrance. We show that the balance between the local acceleration and viscous forces produces a leading-order flow consisting of pure oscillatory motion with axial velocities on the order of a few centimetres per second and amplitudes monotonically decreasing along the length of the canal. We then demonstrate that the nonlinear term associated with the convective acceleration contributes to a second-order correction consisting of a steady streaming that generates a bulk recirculating motion of the CSF along the length of the canal with characteristic velocities two orders of magnitude smaller than the leading-order oscillatory flow. The results of the analysis of this idealized geometry of the spinal canal are shown to be in good agreement not only with experimental measurements in anin-vitromodel but also with radiological measurements conducted in human adults.


2003 ◽  
Vol 19 (3) ◽  
pp. 164-174 ◽  
Author(s):  
Stephen N. Haynes ◽  
Andrew E. Williams

Summary: We review the rationale for behavioral clinical case formulations and emphasize the role of the functional analysis in the design of individualized treatments. Standardized treatments may not be optimally effective for clients who have multiple behavior problems. These problems can affect each other in complex ways and each behavior problem can be influenced by multiple, interacting causal variables. The mechanisms of action of standardized treatments may not always address the most important causal variables for a client's behavior problems. The functional analysis integrates judgments about the client's behavior problems, important causal variables, and functional relations among variables. The functional analysis aids treatment decisions by helping the clinician estimate the relative magnitude of effect of each causal variable on the client's behavior problems, so that the most effective treatments can be selected. The parameters of, and issues associated with, a functional analysis and Functional Analytic Clinical Case Models (FACCM) are illustrated with a clinical case. The task of selecting the best treatment for a client is complicated because treatments differ in their level of specificity and have unequally weighted mechanisms of action. Further, a treatment's mechanism of action is often unknown.


1999 ◽  
Vol 09 (PR10) ◽  
pp. Pr10-129-Pr10-132 ◽  
Author(s):  
J. P. McCarten ◽  
T. C. Jones ◽  
X. Wu ◽  
J. H. Miller ◽  
I. Pirtle ◽  
...  

TAPPI Journal ◽  
2016 ◽  
Vol 15 (5) ◽  
pp. 331-335 ◽  
Author(s):  
LEBO XU ◽  
JEREMY MYERS ◽  
PETER HART

Retention of cationic dispersed rosin size was studied via turbidity measurements on stock filtrate with different alum and dispersed rosin size dosages. Stock charge characteristics were analyzed using both an analysis of charge demand determined via a streaming current detector and an evaluation of zeta potential of the fibers by streaming potential measurement. The results indicated that an optimum amount of alum existed such that good sizing retention was maintained throughout a wide range of dispersed rosin size dosages. However, when an excessive amount of alum was used and fines and colloidal particles were transitioned from anionic to cationic, the cationic size retention was reduced. Laboratory results were confirmed with a paper machine trial. All data suggested that a stock charge study was necessary to identify optimal alum dosage for a cationic dispersed rosin sizing program.


Sign in / Sign up

Export Citation Format

Share Document