Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m(3)) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe. (C) 2010 Elsevier Inc. All rights reserved.

Thick dyke emplacement and internal flow: A structural and magnetic fabric study of the deep-seated dolerite dyke of Foum Zguid (southern Morocco)

Knowledge on forced magma injection and magma flow in dykes is crucial for the understanding of how magmas migrate through the crust to the Earth's surface. Because many questions still persist, we used the long, thick, and deep-seated Foum Zguid dyke (Morocco) to investigate dyke emplacement and internal flow by means of magnetic methods, structural analysis, petrography, and scanning electron microscopy. We also investigated how the host rocks accommodated the intrusion. Regarding internal flow: 1. Important variations of the rock magnetic properties and magnetic fabric occur with distance from dyke wall; 2. anisotropy of anhysteretic remanent magnetization reveals that anisotropy of magnetic susceptibility (AMS) results mainly from the superposition of subfabrics with distinct coercivities and that the imbrication between magnetic foliation and dyke plane is more reliable to deduce flow than the orientation of the AMS maximum principal axis; and 3. a dominant upward flow near the margins can be inferred. The magnetic fabric closest to the dyke wall likely records magma flow best due to fast cooling, whereas in the core the magnetic properties have been affected by high-temperature exsolution and metasomatic effects due to slow cooling. Regarding dyke emplacement, this study shows that the thick forceful intrusion induced deformation by homogeneous flattening and/or folding of the host sedimentary strata. Dewatering related to heat, as recorded by thick quartz veins bordering the dyke in some localities, may have also helped accommodating dyke intrusion. The spatial arrangement of quartz veins and their geometrical relationship with the dyke indicate a preintrusive to synintrusive sinistral component of strike slip.

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Slip flows of Newtonian and viscoelastic fluids in a 4:1 contraction

This work presents a numerical study of the 4:1 planar contraction flow of a viscoelastic fluid described by the simplified Phan-Thien–Tanner model under the influence of slip boundary conditions at the channel walls. The linear Navier slip law was considered with the dimensionless slip coefficient varying in the range ½0; 4500. The simulations were carried out for a small constant Reynolds number of 0.04 and Deborah numbers (De) varying between 0 and 5. Convergence could not be achieved for higher values of the Deborah number, especially for large values of the slip coefficient, due to the large stress gradients near the singularity of the reentrant corner. Increasing the slip coefficient leads to the formation of two vortices, a corner and a lip vortex. The lip vortex grows with increasing slip until it absorbs the corner vortex, creating a single large vortex that continues to increase in size and intensity. In the range De = 3–5 no lip vortex was formed. The flow is characterized in detail for De ¼ 1 as function of the slip coefficient, while for the remaining De only the main features are shown for specific values of the slip coefficient.

Why are slip lengths so large in carbon nanotubes?

The enhanced flow in carbon nanotubes is explained using a mathematical model that includes a depletion layer with reduced viscosity near the wall. In the limit of large tubes the model predicts no noticeable enhancement. For smaller tubes the model predicts enhancement that increases as the radius decreases. An analogy between the reduced viscosity and slip-length models shows that the term slip-length is misleading and that on surfaces which are smooth at the nanoscale it may be thought of as a length-scale associated with the size of the depletion region and viscosity ratio. The model therefore provides a physical interpretation of the classical Navier slip condition and explains why `slip-lengths' may be greater than the tube radius.

Reorganization of a deeply incised drainage: role of deformation, sedimentation and groundwater flow

Deeply incised drainage networks are thought to be robust and not easily modified, and are commonly used as passive markers of horizontal strain. Yet, reorganizations (rearrangements) appear in the geologic record. We provide field evidence of the reorganization of a Miocene drainage network in response to strike-slip and vertical displacements in Guatemala. The drainage was deeply incised into a 50-km-wide orogen located along the North America-Caribbean plate boundary. It rearranged twice, first during the Late Miocene in response to transpressional uplift along the Polochic fault, and again in the Quaternary in response to transtensional uplift along secondary faults. The pattern of reorganization resembles that produced by the tectonic defeat of rivers that cross growing tectonic structures. Compilation of remote sensing data, field mapping, sediment provenance study, grain-size analysis and Ar(40)/Ar(39) dating from paleovalleys and their fill reveals that the classic mechanisms of river diversion, such as river avulsion over bedrock, or capture driven by surface runoff, are not sufficient to produce the observed diversions. The sites of diversion coincide spatially with limestone belts and reactivated fault zones, suggesting that solution-triggered or deformation-triggered permeability have helped breaching of interfluves. The diversions are also related temporally and spatially to the accumulation of sediment fills in the valleys, upstream of the rising structures. We infer that the breaching of the interfluves was achieved by headward erosion along tributaries fed by groundwater flow tracking from the valleys soon to be captured. Fault zones and limestone belts provided the pathways, and the aquifers occupying the valley fills provided the head pressure that enhanced groundwater circulation. The defeat of rivers crossing the rising structures results essentially from the tectonically enhanced activation of groundwater flow between catchments.

Termination of major ductile strike-slip shear and differential cooling along the Insubric line (Central Alps): U-Pb, Rb-Sr and 40Ar/39Ar ages of cross-cutting pegmatites

To constrain the age of strike-slip shear, related granitic magmatism, and cooling along the Insubric line, 29 size fractions of monazite and xenotime were dated by the U-Pb method, and a series of 25 Rb-Sr and Ar-40/Ar-39 ages were measured on different size fractions of muscovite and biotite. The three pegmatitic intrusions analyzed truncate high-grade metamorphic mylonite gneisses of the Simplon shear zone, a major Alpine structure produced in association with dextral strike-slip movements along the southern edge of the European plate, after collision with its Adriatic indenter. Pegmatites and aplites were produced between 29 and 25 Ma in direct relation to right-lateral shear along the Insubric line, by melting of continental crust having Sr-87/Sr-86 between 0.7199 and 0.7244 at the time of melting. High-temperature dextral strike-slip shear was active at 29.2 +/- 0.2 (2 sigma) Ma, and it terminated before 26.4 +/- 0.1 Ma. During dike injection, temperatures in the country rocks of the Isorno-Orselina and Monte Rosa structural units did not exceed approximate to 500 degrees C, leading to fast initial cooling, followed by slower cooling to approximate to 350 degrees C within several million years. In one case, initial cooling to approximate to 500 degrees C was significantly delayed by about 4 m.y., with final cooling to approximate to 300 degrees C at 20-19 Ma in all units. For the period between 29 and 19 Ma, cooling of the three sample localities was non-uniform in space and time, with significant variations on the kilometre scale. These differences are most likely due to strongly varying heat flow, and/or heterogeneous distribution of unroofing rates within the continuously deforming Insubric line. If entirely ascribed to differences in unroofing, corresponding rates would vary between 0.5 and 2.5 mm/y, for a thermal gradient of 30 degrees/km.

Separation of Gas Bubbles from Circulation Waters and Pulp Suspensions in Open Channel Flow

The objective of this thesis was to study the removal of gases from paper mill circulation waters experimentally and to provide data for CFD modeling. Flow and bubble size measurements were carried out in a laboratory scale open gas separation channel. Particle Image Velocimetry (PIV) technique was used to measure the gas and liquid flow fields, while bubble size measurements were conducted using digital imaging technique with back light illumination. Samples of paper machine waters as well as a model solution were used for the experiments. The PIV results show that the gas bubbles near the feed position have the tendency to escape from the circulation channel at a faster rate than those bubbles which are further away from the feed position. This was due to an increased rate of bubble coalescence as a result of the relatively larger bubbles near the feed position. Moreover, a close similarity between the measured slip velocities of the paper mill waters and that of literature values was obtained. It was found that due to dilution of paper mill waters, the observed average bubble size was considerably large as compared to the average bubble sizes in real industrial pulp suspension and circulation waters. Among the studied solutions, the model solution has the highest average drag coefficient value due to its relatively high viscosity. The results were compared to a 2D steady sate CFD simulation model. A standard Euler-Euler k-ε turbulence model was used in the simulations. The channel free surface was modeled as a degassing boundary. From the drag models used in the simulations, the Grace drag model gave velocity fields closest to the experimental values. In general, the results obtained from experiments and CFD simulations are in good qualitative agreement.

Genesis and Evolution of Brittle Structures In Southwestern Finland and Western South Africa. Insights into Fault Reactivation, Fluid Flow and Structural Maturity in Precambrian Cratons

The bedrock of old crystalline cratons is characteristically saturated with brittle structures formed during successive superimposed episodes of deformation and under varying stress regimes. As a result, the crust effectively deforms through the reactivation of pre-existing structures rather than by through the activation, or generation, of new ones, and is said to be in a state of 'structural maturity'. By combining data from Olkiluoto Island, southwestern Finland, which has been investigated as the potential site of a deep geological repository for high-level nuclear waste, with observations from southern Sweden, it can be concluded that the southern part of the Svecofennian shield had already attained structural maturity during the Mesoproterozoic era. This indicates that the phase of activation of the crust, i.e. the time interval during which new fractures were generated, was brief in comparison to the subsequent reactivation phase. Structural maturity of the bedrock was also attained relatively rapidly in Namaqualand, western South Africa, after the formation of first brittle structures during Neoproterozoic time. Subsequent brittle deformation in Namaqualand was controlled by the reactivation of pre-existing strike-slip faults.In such settings, seismic events are likely to occur through reactivation of pre-existing zones that are favourably oriented with respect to prevailing stresses. In Namaqualand, this is shown for present day seismicity by slip tendency analysis, and at Olkiluoto, for a Neoproterozoic earthquake reactivating a Mesoproterozoic fault. By combining detailed field observations with the results of paleostress inversions and relative and absolute time constraints, seven distinctm superimposed paleostress regimes have been recognized in the Olkiluoto region. From oldest to youngest these are: (1) NW-SE to NNW-SSE transpression, which prevailed soon after 1.75 Ga, when the crust had sufficiently cooled down to allow brittle deformation to occur. During this phase conjugate NNW-SSE and NE-SW striking strike-slip faults were active simultaneous with reactivation of SE-dipping low-angle shear zones and foliation planes. This was followed by (2) N-S to NE-SW transpression, which caused partial reactivation of structures formed in the first event; (3) NW-SE extension during the Gothian orogeny and at the time of rapakivi magmatism and intrusion of diabase dikes; (4) NE-SW transtension that occurred between 1.60 and 1.30 Ga and which also formed the NW-SE-trending Satakunta graben located some 20 km north of Olkiluoto. Greisen-type veins also formed during this phase. (5) NE-SW compression that postdates both the formation of the 1.56 Ga rapakivi granites and 1.27 Ga olivine diabases of the region; (6) E-W transpression during the early stages of the Mesoproterozoic Sveconorwegian orogeny and which also predated (7) almost coaxial E-W extension attributed to the collapse of the Sveconorwegian orogeny. The kinematic analysis of fracture systems in crystalline bedrock also provides a robust framework for evaluating fluid-rock interaction in the brittle regime; this is essential in assessment of bedrock integrity for numerous geo-engineering applications, including groundwater management, transient or permanent CO2 storage and site investigations for permanent waste disposal. Investigations at Olkiluoto revealed that fluid flow along fractures is coupled with low normal tractions due to in-situ stresses and thus deviates from the generally accepted critically stressed fracture concept, where fluid flow is concentrated on fractures on the verge of failure. The difference is linked to the shallow conditions of Olkiluoto - due to the low differential stresses inherent at shallow depths, fracture activation and fluid flow is controlled by dilation due to low normal tractions. At deeper settings, however, fluid flow is controlled by fracture criticality caused by large differential stress, which drives shear deformation instead of dilation.

The effect of flow field & geometry on the dynamic contact angle

A number of recent experiments suggest that, at a given wetting speed, the dynamic contact angle formed by an advancing liquid-gas interface with a solid substrate depends on the flow field and geometry near the moving contact line. In the present work, this effect is investigated in the framework of an earlier developed theory that was based on the fact that dynamic wetting is, by its very name, a process of formation of a new liquid-solid interface (newly “wetted” solid surface) and hence should be considered not as a singular problem but as a particular case from a general class of flows with forming or/and disappearing interfaces. The results demonstrate that, in the flow configuration of curtain coating, where a liquid sheet (“curtain”) impinges onto a moving solid substrate, the actual dynamic contact angle indeed depends not only on the wetting speed and material constants of the contacting media, as in the so-called slip models, but also on the inlet velocity of the curtain, its height, and the angle between the falling curtain and the solid surface. In other words, for the same wetting speed the dynamic contact angle can be varied by manipulating the flow field and geometry near the moving contact line. The obtained results have important experimental implications: given that the dynamic contact angle is determined by the values of the surface tensions at the contact line and hence depends on the distributions of the surface parameters along the interfaces, which can be influenced by the flow field, one can use the overall flow conditions and the contact angle as a macroscopic multiparametric signal-response pair that probes the dynamics of the liquid-solid interface. This approach would allow one to investigate experimentally such properties of the interface as, for example, its equation of state and the rheological properties involved in the interface’s response to an external torque, and would help to measure its parameters, such as the coefficient of sliding friction, the surface-tension relaxation time, and so on.

Liquid flow over a substrate structured by seeded nanoparticles

Recent experiments have demonstrated that nanoparticles which sparsely distributed over a solid substrate can substantially change the flow conditions at the solid surface in the presence of slip. Inspired by these observations, the flow past tiny particles seeded on a solid substrate is investigated theoretically in the framework of an interface formation model. It has been shown, that even a single seeded nanoparticle can reduce significantly the measurable tangential component of hydrodynamic velocity at the substrate and affect the amount of the observed apparent slippage of the liquid. The effect from the particle manifests in a form of a long relaxation tail defined by the characteristic time of the interface formation process. A comparison with experiments has demonstrated a good agreement between theoretically predicted and experimentally observed values of the relaxation tail length scale.

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

Computations of the flow past a still sphere at moderate reynolds numbers using an immersed boundary method

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Drag reduction phenomenon in viscous oil-water dispersed pipe flow: Experimental investigation and phenomenological modeling

An experimental study on drag-reduction phenomenon in dispersed oil-water flow has been performed in a 26-mm-i.d. Twelve meter long horizontal glass pipe. The flow was characterized using a novel wire-mesh sensor based on capacitance measurements and high-speed video recording. New two-phase pressure gradient, volume fraction, and phase distribution data have been used in the analysis. Drag reduction and slip ratio were detected at oil volume fractions between 10 and 45% and high mixture Reynolds numbers, and with water as the dominant phase. Phase-fraction distribution diagrams and cross-sectional imaging of the flow suggested the presence of a higher amount of water near to the pipe wall. Based on that, a phenomenology for explaining drag reduction in dispersed flow in a flow situation where slip ratio is significant is proposed. A simple phenomenological model is developed and the agreement between model predictions and data, including data from the literature, is encouraging. (c) 2011 American Institute of Chemical Engineers AIChE J, 2012