Effects of Bubbles on the Hydraulic Conductivity of Porous Materials Theoretical Results
|
|
- Johnathan Anthony
- 6 years ago
- Views:
Transcription
1 Transport in Porous Media 52: 51 65, Kluwer Academic Publishers. Printed in the Netherlands. 51 Effects of Bubbles on the Hydraulic Conductivity of Porous Materials Theoretical Results A. G. HUNT and MICHAEL MANGA Department of Earth and Planetary Science, University of California, Berkeley, CA 94720, U.S.A. (Received: 28 August 2001; in final form: 14 August 2002) Abstract. In a porous material, both the pressure drop across a bubble and its speed are nonlinear functions of the fluid velocity. Nonlinear dynamics of bubbles in turn affect the macroscopic hydraulic conductivity, and thus the fluid velocity. We treat a porous medium as a network of tubes and combine critical path analysis with pore-scale results to predict the effects of bubble dynamics on the macroscopic hydraulic conductivity and bubble density. Critical path analysis uses percolation theory to find the dominant (approximately) one-dimensional flow paths. We find that in steady state, along percolating pathways, bubble density decreases with increasing fluid velocity, and bubble density is thus smallest in the smallest (critical) tubes. We find that the hydraulic conductivity increases monotonically with increasing capillary number up to Ca 10 2, but may decrease for larger capillary numbers due to the relative decrease of bubble density in the critical pores. We also identify processes that can provide a positive feedback between bubble density and fluid flow along the critical paths. The feedback amplifies statistical fluctuations in the density of bubbles, producing fluctuations in the hydraulic conductivity. Key words: bubbles, critical path analysis, percolation theory, hydraulic conductivity. 1. Introduction The motion of bubbles in tubes is frequently studied for its value in understanding pore-scale multi-phase flow processes (Olbricht, 1996). While such model geometries are obviously gross oversimplifications of true pore geometries, the effects of bounding walls, viscosity ratio, volume fraction, bubble volume, and interfacial tension can all be considered. The model problem may thus provide physical insight into pore-scale hydrodynamics. The presence of a bubble in a pore affects the flow in that pore, and consequently global properties such as a macroscopic hydraulic conductivity. Of course the flow in a given pore is ultimately given in terms of the macroscopic hydraulic conductivity. When bubbles are present, the flow is not of simple Darcian form. In detail, because both the pressure drop across a bubble and the speed at which the bubble moves are nonlinear functions of the fluid flux, Darcy s law is not satisfied, even at the pore-scale. The goal of the present study is to use pore-scale results to
2 52 A. G. HUNT AND MICHAEL MANGA deduce properties of the porous material, in particular the macroscopic hydraulic conductivity and the density distribution of bubbles. One approach to calculate the linear (hydraulic) conductivity of disordered networks with a wide spread in local linear conductances is to use percolation theory (Kirkpatrick, 1973a, b; Seager and Pike, 1974; Pollak, 1985; Bernabe and Bruderer, 1998; Hunt, 2001). A general nonlinear treatment of transport using percolation theory is not universally agreed on (Pollak and Riess, 1976; Shklovskii, 1976; Aladashvili and Adamia, 2002). Nevertheless, as we argue here, percolation theory can still be applied when the local conductances are nonlinear so long as neither the form of the nonlinearity nor the choice of flow paths depend on the flow. Previous analytical and numerical studies (e.g. Ratulowski and Chang, 1989) of bubbles in porous materials have considered the effects of bubbles on the macroscopic hydraulic conductivity for ordered porous materials with a uniform (or nearly uniform) pore size. In the present work we apply percolation theory in order to derive analogous analytical results for a medium with a much wider range of tube sizes. An important aspect of our derivation is that the density of bubbles is no longer required to be uniform, rather it is determined dynamically. While our model is conceptually similar to droplet train models (e.g. Rossen, 1990; Foulser et al., 1991; Babchin and Yuan, 1997) we explicitly account for the nonlinear dynamics of the bubbles. The structure of the calculations depends on both the pore-scale dynamics and the application of percolation theory to find the macroscopic response. The particular application of percolation theory depends on the model chosen. Thus, in this paper we first describe the motion of a single bubble in a single tube (Section 2) and introduce a model for the pore space topology (Section 3). Then we introduce the relevant concepts of percolation theory, of which the most important are that certain (approximately one-dimensional, 1D) pathways dominate the flow, and that the properties of flow on these pathways are dominated by particular pores (Section 4). Given the assumption of a 1D pathway, it is possible to derive a relationship between the bubble density and tube size (Section 5). In Section 6 we use the bubble density, the local nonlinear dynamics, and the results from percolation theory to obtain an estimate of the effects of bubbles on the macroscopic hydraulic conductivity. Finally, in Section 7 we acknowledge that steady-state conditions do not always exist in real flows, and discuss the possible consequences for bubble dynamics. 2. Single Bubble Dynamics In order to develop a model for the motion of a bubbly fluid through a porous material, we represent a possibly complex pore geometry (Figure 1(a)) as a network of cylindrical tubes (Figure 1(b)). This geometric approximation is made in order to simplify the analysis. It also allows us to use well-established pore-scale results. While the geometry shown in Figure 1(b) may not resemble typical porous
3 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 53 Figure 1. (a) Schematic two-dimensional view of non-cylindrical bubbles in non-cylindrical pore space. (b) An idealized version of (1a) commensurate with cylindrical geometry, represented as a network of tubes. materials, for example, soils, we will argue in Section 4 that the macroscopic properties of the flow (in particular the hydraulic conductivity) are dominated by processes in the smaller pores through which fluid is forced to flow. For such pores our model problem may be a reasonable approximation that captures the essential physics of flow, at least for the case in which the gas phase exists as discrete bubbles. Bretherton (1961) found that an infinitely long bubble moves at a speed u which exceeds the fluid velocity v by an amount proportional to the bubble speed, with the proportionality constant defined to be w. The relationship between u and v is therefore u = v 1 w. (1) w depends on the capillary number, Ca. To leading order, w = 1.29(3Ca) 2/3, (2) wherecaisdefinedby Ca = µu σ. (3) Here µ is the dynamic viscosity of the suspending fluid, and σ is the interfacial tension between the bubble and the suspending fluid. Equation (2) is an excellent approximation up to Ca Although Bretherton s (1961) analysis considered infinitely long bubbles, Equation (2) applies for bubbles with volumes as small as 0.95(4π/3)r 3,wherer is the tube radius (Olbricht, 1996). Because Ca depends on u, Equation (1) represents an implicit equation for w. A solution for w in terms of the fluid velocity, v, accurate over the entire range of capillary numbers investigated, is given in Stark and Manga (2000). Here, however, we solve Equation (1) numerically whenever a value is required.
4 54 A. G. HUNT AND MICHAEL MANGA The pressure drop, P b, across a bubble within a tube of radius r is P b = 3.58 σ r (3Ca)2/3. (4) P b can be rewritten in terms of w and v by using Equation (3) for Ca, and Equation (1) for u, 7.45(σ/r)(µv/σ )2/3 P b = 7.45(σ/r)(Ca ) 2/3, (5) [1 w] 2/3 [1 w] 2/3 where Ca is defined based on the fluid speed v rather than the bubble speed u. The dynamics described by Equations (2) and (4) apply for surfactant-free discrete bubbles and if non-hydrodynamic effects can be ignored (reasonable for length scales greater than 10 nm). Many practical applications involving bubbles in porous materials use foams stabilized with surfactants, however (e.g. Hazlett and Furr, 2000). Although effects of surfactants on bubble motion have been studied (e.g. Ratulowski and Chang, 1990; He et al., 1991; Borhan and Mao, 1992), we nevertheless choose to neglect such effects in the interests of clarity and simplicity. One important final result is that in a chain of identical bubbles, Ratulowski and Chang (1989) noted that each bubble behaves independently. Thus, extension to treat a large number of bubbles is possible by multiplying the effects of one bubble by the number of bubbles, N. In the rest of this paper we consider the implications of bubble dynamics. In particular we will see that for steady-state flow along 1D paths, the conservation of bubbles, and the conservation of fluid, have some interesting consequences for the variation of bubble density. The relevance of 1D flow paths is understood within the concept of percolation theory. Percolation theory has a second important implication that the dynamics of bubbles in a critical tube (of a particular radius) dominate the hydraulic properties of the entire network. In order to introduce percolation theory, however, we must first briefly define the model pore space. 3. Pore Space Topology Some justification for fractal pore space descriptions exists (Turcotte, 1986; Sahimi, 1995). We thus consider a tube network with some fractal characteristics, and use a basic result for the porosity (Rieu and Sposito, 1991) as a test. The individual tubes are right circular cylinders with radius r, length Cr,andvolume,V = Cπr 3, with C large. The probability to measure a tube radius within dr of r is W(r)dr r 1 D. In detail, W(r)dr = 3 D rm 3 D Cπ r 1 D dr. (6) This power law is compatible with a (random) fractal (of dimensionality D) arrangement of tubes with radii r 0 r r m. r 0 and r m are the radii which bound the
5 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 55 validity of a fractal description. The power 1 D is used rather than D because, in contrast to Rieu and Sposito (1991), we are not modeling a discrete fractal, but a continuous distribution (Hunt, 2001). The proportionality constant, however, is chosen so as to produce the known result for the porosity, φ (Rieu and Sposito, 1991). The porosity is found from Equation (6) by integrating over all allowed r the product of VW(r), rm ( ) 3 D 3 D r 1 D Cπr 3 r0 dr = 1 = φ. (7) r 0 Cπr 3 D m Though the value of the permeability itself will depend on pore space topology and geometry (Katz and Thompson, 1986; Friedman and Seaton, 1998; Hunt, 2001), the structure of the result for the change in permeability due to bubbles will not depend sensitively on the choice of representation of the pore space. r m 4. Percolation Theory and Critical Path Analysis One of the most accurate means to calculate transport coefficients, including the electrical and hydraulic conductivities, in disordered media (Kirkpatrick, 1973a, b; Seager and Pike, 1974; Bernabe and Bruderer, 1998; Hunt, 2001), is to apply a version of percolation theory called critical path analysis (Ambegaokar et al., 1971; Shklovskii and Efros, 1971; Pollak, 1972; Katz and Thompson, 1986). Critical path analysis, explained here by an analogy to the electrical resistivity problem, exploits the following simple relationship. If two resistances differing by many orders of magnitude are configured in parallel, virtually all the current flows through the smaller resistor, and the effective resistance is equal to the smaller resistance value. If two such resistances are configured in series, virtually the entire potential drop occurs across the largest resistor, and the effective resistance is equal to the larger resistance value. By use of percolation theory it is possible to select that resistance value, out of a wide range of resistance values, that governs the current. Specifically, percolation theory can be used to find the largest resistance value, R c (called the critical resistance), on the path that has the smallest possible total resistance on a transit through the system. In keeping with the above simplified discussion, other paths through the system have a much higher resistance, and may be neglected. But the most important potential drops occur over the largest resistors on the critical path. The complete treatment is more complex (e.g. Friedman and Pollak, 1981; Hunt, 2001; Skal, 2001), and the hydraulic conductivity involves additional length scales. But in many applications, knowing the critical value of the hydraulic conductance, including, if necessary, how it changes with changing conditions, has proved sufficient (see review by Hunt, 2001), and will be assumed sufficient for the present application as well. The particular application of percolation theory depends on whether or not the pore space is represented as a regular network. On a regular network, with bonds describing resistances of varying strengths, the quantile of the distribution
6 56 A. G. HUNT AND MICHAEL MANGA of resistances defining R c is given by the bond percolation fraction, p c. When the pore space is treated as a fractal, it has proven advantageous to apply continuum percolation (Hunt, 2001; Hunt and Gee, 2001). The reason for this is that fractal geometries require pores at all length scales to have the same aspect ratio, and inclusion of, say, two orders of magnitude of pore sizes then requires a distribution of pore lengths over two orders of magnitude. Such a wide distribution of pore lengths is difficult to constrain on a regular lattice (with translationally invariant coordination number). When continuum percolation is applied, it is necessary to consider the total volume occupied by pores (Hunt, 2001). An example can be the network of tubes considered here. When a portion of the tube network, including tubes with the largest radii, is emplaced spatially at random (with no spatial correlations), a particular volume fraction is associated with the tubes. If enough tubes are included so that the set of tubes interconnects to infinite size, the critical value, α c,ofthe volume fraction has been reached, and the tube volumes are said to percolate. This critical volume fraction is not known a priori for arbitrary geometries and topologies of connection. Nevertheless a typical value of α c 0.15 can be used with some justification (Stauffer, 1975) although higher values are also possible (Balberg, 1985; Moldrup et al., 2001). What is the critical radius for percolation when the tube size distribution is consistent with a fractal, as described in the previous section? To find the critical radius, r c, it is necessary to integrate over the distribution of r-values from the critical value to the largest allowed in the distribution, and set that pore space volume equal to the critical volume fraction, α c, rm 3 D rm 3 D r 2 D dr = α c. (8) r c Solving Equation (8) for r c in terms of α c, we obtain r c = r m (1 α c ) 1/(3 D). (9) This is the smallest radius of any tube that the fluid must flow through to get from one side of a system to the other. It has been conjectured that for fractal soils, the fractal dimensionalities of the solid volume and the pore space are very similar (Rieu and Sposito, 1991). For soil particle-size distributions, typical values of D are 2.8 ± 0.1 (Wuet al., 1993; Bittelli et al., 1999) leading to values of r c which are between 2 and 5 times smaller than r m,ifα c = 0.15 is used. Because such fractal descriptions may be valid over orders of magnitude of pore sizes (Bittelli et al., 1991; Perrier et al., 1996; Hunt and Gee, 2001) r c can often be more than an order of magnitude larger than r 0. The value of the critical radius for percolation is not changed by considering nonlinearities in flow due to bubbles under the following conditions, both of which should apply for the problem being considered:
7 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS there is a unique relationship between the individual hydraulic conductance and the tube radius, that is, it is not possible to get two different conductivities with one tube radius, nor is it possible to get the same conductance from two different tube radii, 2. the tube conductance is an increasing function of tube radius. The first condition holds whenever the initial bubble distribution is constrained to be homogeneous in space. The reason is that it is, in this case, always possible to use transformation of variables in the integration to return to the same integral condition on r. Nevertheless, the condition is perhaps most easily understood by considering its violation through, for example, an inhomogeneous bubble distribution. If, at low Ca, two tubes with very slightly differing r values have greatly differing bubble contents, the tube with the larger r will have greater flow only if it has fewer bubbles. Under these circumstances, both the number of bubbles and the tube dimensions can be considered as independent random variables, and it is necessary to find the critical value of the hydraulic conductance g(r,φ b ).Such problems coupling different random variables tend, in 3D, to lead to larger values of a critical conductance, but smaller values of the smallest r value encountered along a percolation path. An analogous effect is considered in detail in Hunt (2001) where for the electrical conductivity both the radius and the ionic strength are treated as random variables. If the first condition mentioned above holds, but the latter condition is violated, it is necessary only to change the order of integration from the original smallest r required on the path to r m, to a final range from r 0 to the largest r required on the path. The calculated value of r c is useful to define a critical value of the hydraulic conductance, gc 0 = πρr4 c = πρr3 c 8µl c 8µC, (10) where the superscript 0 indicates that this is the value of g c in the absence of bubbles, and subscripts c indicate values for the critical pore. An added complication is the necessity, in steady state, to guarantee a time-independent flux of bubbles as well as a time-independent fluid flux. This problem is dealt with next. 5. Spatial Distribution of Bubbles We consider the case for which the bubble distribution has no random component. Associated with the bubbles is a fraction of the total pore space, which we call the bubble pore space, with bubble porosity, φ b, normalized to the porosity of the medium so that 0 φ b < 1. The number, N, of bubbles in a given tube of volume V depends on the bubble volume, V b (4π/3)r 3 b,aswellasonφ b, N = φ b V V b. (11)
8 58 A. G. HUNT AND MICHAEL MANGA Substituting the tube volume V = Cr(πr 2 ) for a tube of radius r and length Cr (constant aspect ratio, C) gives N = πcr3 φ b = 3 ( ) r 3 V b 4 C φ b. (12) r b The number of bubbles per unit length in such a tube is πr 2 φ b. (13) V b The actual length of a bubble is d = V b /πr 2, (14) so that the total length taken up by bubbles is Nd, Nd = πcr3 V b V b πr 2 φ b = Crφ b. (15) Now it is possible to write for the effective length (bubble-free portion) of the tube, l = Cr[1 φ b ]. (16) By using percolation theory to find the macroscopic hydraulic conductivity we can make the first order approximation that flow occurs along one-dimensional paths, and that off such paths there is comparatively little flow. As a consequence, it is possible to derive some simple relationships between bubble density and fluid flow velocities for tubes along the percolation path. Consider a tube of radius r on the percolation path. The mean distance between bubbles in such a tube is d s = V b. (17) πr 2 φ b In time t the mean number of bubbles escaping a tube of radius, r, is u t V b /(πr 2 ) φ b. (18) This number must be identical for all tubes on the percolation path (assuming, in a lowest order approximation, 1D flow along this path), and therefore uπr 2 φ b = v 1 w πr2 φ b = constant. (19) In addition, the flux of fluid volume, Q, must be identical everywhere along this path (under steady-state conditions). Using the result that v = Q/πr 2 yields φ b (1 w). (20)
9 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 59 Thus on the percolation path (with the largest fluid speeds, and hence, the largest bubble fluxes) bubble density in steady-state flow is not uniform, but depends on the tube radius, r,thatis,φ b = φ b (r). Under steady-state conditions the mean concentration of bubbles along the percolation paths is the same as the injection concentration, φ b. Rewrite Equation (20) as φ b (r) = φ 0 (1 w(r)) (21) with φ 0 a constant. The integral of φ b over all tubes found along the critical path for percolation is equal to the mean value, φ b, of the bubble porosity by assumption, so that it is possible to solve for φ 0, rm φ b r φ 0 = c W(r)dr rm r c (1 w(r))w(r) dr. (22) The upper limit of integration, r m, is the largest tube radius in the system, while the lower limit, r c, is the smallest tube encountered along the percolation path. Using Equation (22) for φ 0, we obtain φ b (r) = φ b[1 w(r)] r m r c W(r)dr rm (1 w)w(r) dr r c φ b(1 w(r)), (23) 1 w where w is the average value of w along the percolation path. For a system with a wide distribution of pore sizes, such as a fractal distribution with 3 D 0.2, w c /w>1sothat φ b (r c ) φ b (1 w(r c )). (24) 6. Application of Percolation Theory to Flow with Nonlinear Bubble Dynamics The fluid flux in a tube of radius r is given in terms of the pressure drop in the tube, less the pressure drops across the bubbles. Because the bubbles are identical, we need to multiply P b in Equation (5) by N from Equation (12), to find the total pressure drop across the bubbles, and subtract this value from the total pressure drop in the tube. The remaining pressure drop, which acts on the fluid, thus acts only on the length of tube which is bubble-free. So the length, l, which appears in Equation (10) must be replaced by l from Equation (16). Now it is possible to write a self-consistent equation for Q for any of the tubes on the percolation path, Q = πρr 4 8µCr[1 φ b (r)] [ 1 V V b φ b 1 P 7.54σ r [1 w(r)] 1/3 [1 w] 2/3 (Ca ) 2/3 ] P, (25)
10 60 A. G. HUNT AND MICHAEL MANGA where V is the volume of the tube under consideration, and φ b is given in Equation (23). Using the definition (Equation (10)) of the hydraulic conductance in the absence of bubbles, Equation (25) can be rewritten, Q = g 0 1 φ b (r) [ 1 V σ [1 w(r)] 1/3 φ V b (Ca ) 2/3 b P r [1 w] 2/3 ] P. (26) In order to characterize our results, we use the procedure of Stark and Manga (2000), and determine the effect of the presence of bubbles on the hydraulic conductivity by comparing the conductivity with bubbles to its value without bubbles. This procedure amounts to setting the values of P across the critical conductance (and hence, critical tube) in the presence and in the absence of bubbles equal, and solving for the ratio of the resulting discharges, Q. We thus obtain 8µC[1 (1 φ b (r c )]Q c πρr 3 c + V cφ b σ V b r c [ [ ] (Ca w(rc ) 1/3 ] c )2/3 (1 w) 2 = 8µCQ0 c. (27) πρrc 3 Equation (27) can be solved for the ratio of Q c /Q 0 c to yield { [ ] g c 1 w(rc ) = 1 φ gc 0 1 w b + πrc 3φ [ [ ] b (Ca w(rc ) 1/3 ]} 1 c V ) 1/3. (28) b (1 w) 2 Equation (28) provides an expression for effects of bubbles on the hydraulic conductivity, when combined with information about the pore geometry, in particular r c from Equation (9) and w from Equation (23). g c /gc 0 is equivalent to the liquid phase relative permeability; a relative permeability for the gas phase is not meaningful as the bubbles are not continuous. One last point needs to be addressed. It is tacitly assumed above that by not changing the value of the externally applied pressure, the value of P that appears in Equations (25) and (26) is not changed. This is justified within the basic approximation of critical path analysis that all (or nearly all) the pressure drop in the network occurs across critical valued resistances along the percolation path. Thus it is necessary only to find out how far apart such critical resistance values occur, a distance we call l c. Given a value for l c, it is possible to write ( ) lc P c = P, (29) z where z is the total width of the system (over which the pressure difference P is applied), and P c is the pressure difference applied across a critical tube (with r c as calculated in the absence of bubbles). Since, however, there has been no
11 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 61 Figure 2. Modification of the system hydraulic conductivity as a function of Ca referenced to the fluid flow velocity, and as a function of bubble density. change in r c due to the introduction of the bubbles (only in Q c ), it is reasonable to assume that l c does not change either. It is therefore unnecessary to know what P is because the same fraction of the macroscopic pressure drop is still recorded across the critical resistances. The previous argument implies that the change in conductance of the critical tube (in the form of a ratio) is likely to be the same as the result for the entire network. In reality, however, the presence of bubbles reduces the conductance of smaller tubes less (with larger Ca, see Figure 2) than that of larger tubes, thus changing P across the critical tubes slightly upon introduction of bubbles. The results of Equation (28) for various bubble concentrations as a function of Ca are shown in Figure 2. Our numerical solution uses the approximation, w(r c )/w 1. For Ca less than 10 2, the hydraulic conductivity is depressed by the presence of bubbles. The effect is most pronounced at small values of the capillary number and large bubble concentrations, a well-known result for foams (e.g. Ratulowski and Chang, 1989; Patzek, 2001). For example, for bubble concentrations φ b = 0.7andCa 10 7, the reduction is two orders of magnitude. At values of Ca > 10 2 bubbles enhance the hydraulic conductivity. This enhancement is strongest at the largest bubble concentrations. The form of the curves for g c /g 0 c shown in Figure 2 is similar to those obtained numerically by Stark and Manga (2000), especially for Ca < 10 3.
12 62 A. G. HUNT AND MICHAEL MANGA 7. Unsteady Flow and Flow Instabilities Up to this point, our analysis assumes that steady-state conditions are obeyed. But an interesting consequence of the present results is that they support the existence of flow fluctuations that are amplified by the nonlinear dynamics. To see this consider again Figure 2 in the regime of low Ca. Assume a drastically simplified model of the percolation path, with alternating tubes of radii r c and 5r c. The tubes with radius 5r c have 125 times the conductance, so the potential drops across them are only 1/125 times as large, in the mean, in order to produce the same flow, Q. The tubes with r = r c, on the other hand, have a 25 times larger value of a mean v, on account of their smaller r-values, and thus have larger values of the capillary number (both Ca and Ca ). Consider now the possibility that a fluctuation in bubble numbers occurs such that the number of bubbles in the tubes with r = r c drops. Their conductances are increased and the net flow increases as well, whether or not the bubble number in the larger tubes increases. Because, to first approximation, 125/126 of the pressure drop occurs across the critical tube, the larger tubes can increase their flow (and pressure drop) significantly without any significant change in the pressure drop across the critical tube. Thus v increases for both tubes, and by the same fractional increment. Such a proportional increase in v, however, causes a disproportionate increase in bubble velocity, u, because the critical tubes have a larger v, and the rate at which u increases with v increases with increasing u. Thus the bubble velocity is increased by a larger fraction in the critical tubes than in the larger tubes. Such a disproportionate increase tends to reduce the number of bubbles in the critical tubes even further, and this positive feedback leads to amplifications of flow fluctuations. By an analogous argument, an increase in bubble concentration in the critical tubes (at low Ca) leads to a further increase, with the result that these tubes can become blocked, a process we neglect in our analysis. In response, the fluid may find another path (with a larger value of the maximum resistance) leading to a fluctuation in flow paths, and hence the hydraulic conductivity. When bubbles are first injected into a porous medium, each blocking of a critical pore may be accompanied by the trapping of air in pores. If so, we expect to see a trend of diminishing water content and associated reduction in hydraulic conductivity along with the flow fluctuations, such as those described in Faybishenko (1995, 1999, 2000), which include both laboratory and field measurements. Stark and Manga (2000) simulated numerically nonlinear bubble flow in the same model considered here, but with less extreme geometry (ratio of maximum to minimum pore size is of a factor 3), and observed that, At low Ca, the presence of one or more bubbles on a high-flow pathway can significantly reduce the flow through that pathway. Other bubbles are thus more likely to leave that high-flow pathway, and the presence of nodes enables bubbles to avoid congested pathways." Indeed, in Figure 2(a) of Stark and Manga (2000) there are significant mean velocity fluctuations that are much larger than those due to variations in bubble
13 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 63 content alone (compare with the magnitude of the bubble density fluctuations in their Figure 2(b)). For cases in which the critical tubes have large Ca values (above 10 2 ) a reduction in bubble density in the critical tubes leads to a diminution of flow, because the flow here increases with bubble density. Consequently a greater reduction in u occurs in the critical tubes (with the larger Ca) and bubble concentrations increase in the critical tubes. Thus, when the critical tubes are in the large Ca region, fluctuations in bubble density are attenuated. This argument does not depend on the value of Ca in the larger tubes. 8. Conclusions We have treated nonlinear pore-scale dynamics of bubble and fluid flow in porous media using a network of tubes with some characteristics of a fractal. The calculations of the macroscopic system response utilized critical path analysis, defined in terms of percolation theory. In this analysis the predominant fluid flow and bubble transport occur along quasi-one-dimensional paths of interconnected tubes. The smallest tubes through which the bubbly fluid must flow dominate both the flow and bubble transport properties. Conservation of bubbles in steady state requires a systematic variation of bubble density with tube radius on these dominant paths: larger tubes have smaller capillary numbers and higher bubble densities. We find that K is reduced by the presence of bubbles when the critical tube for flow is in a low capillary number regime (Ca < 10 3 ), but is enhanced for larger values of Ca. In both cases, the effects increase with increasing bubble density. The dynamic redistribution of bubbles has little effect on the system-wide hydraulic conductivity, except at Ca > 10 3, where it reduces the enhancement of the conductivity resulting from the presence of bubbles. Although the preceding results should hold for steady-state flow, we find that the nonlinearities associated with bubble transport amplify flow fluctuations due to fluctuations in bubble density when Ca in the critical tube is less than At higher values of Ca,theenhancement of the hydraulic conductivity from the presence of bubbles suppresses flow variations due to fluctuations in bubble density. We must reemphasize the significant geometric simplifications in the model problem. Also, the critical path analysis is essentially one-dimensional. Moreover, approximating the dispersed gas phase as discrete bubbles will not describe the topology of the gas phase in the general case (Avraam and Payatakes, 1995), although it may describe situations such as foam mobilization (Rossen et al., 1994). Still, the theoretical results developed here provide a framework for interpreting experimental results and numerical simulations. Differences between measurements and our results should in principle provide insight into those processes not included in the model.
14 64 A. G. HUNT AND MICHAEL MANGA References Aladashvili, D. I. and Adamia, Z. A.: 2001, High-field hopping and hopping domains, Phil. Mag. B 81, Ambegaokar, V. N., Halperin, B. I. and Langer, J. S.: 1971, Hopping conductivity in disordered systems, Phys. Rev. B 4, Avraam, D. G. and Payatakes, A. C.: 1995, Flow regimes and relative permeabilities during steadystate two-phase flow in porous media, J. Fluid Mech. 293, Babchin, A. J. and Yuan, J.-Y.: 1997, On the capillary coupling between two phases in a droplet train model, Trans. Porous Med. 26, Balberg, I.: 1985, Universal percolation threshold limits in the continuum, Phys. Rev. B 31, Bernabe, Y. and Bruderer, C.: 1998, Effect of the variance of pore size distribution on the transport properties of heterogeneous networks, J. Geophys. Res. 103, 513. Bittelli, M., Campbell, G. S. and Flury, M.: 1999, Characterization of particle-size distribution in soils with a fragmentation model, Soil Sci. Soc. Am. J. 63, Borhan, A. and Mao, C. F.: 1992, Effects of surfactants on the motion of drops through circular tubes, Phys. Fluids A 4, Bretherton, F. P.: 1961, The motion of long bubbles in tubes, J. Fluid Mech. 10, Faybishenko, B. A.: 1995, Hydraulic behavior of quasi-saturated soils in the presence of entrapped air. Laboratory investigations, Water Resour. Res. 31, Faybishenko, B. A.: 1999, Comparison of laboratory and field methods for determining quasiunsaturated hydraulic conductivity of soils, in: Proceedings of the International Workshop on Characterization and Measurement of the Hydraulic Properties of Unsaturated Porous Media, pp Faybishenko, B. A.: 2000, Chaotic dynamics in flow through unsaturated fractured media, LBNL Report Foulser, R. W. S., Goodyear, S. G. and Sims, R. J.: 1991, New concepts in relative permeabilities at high capillary number for surfactant flooding, Trans. Porous Med. 6, Friedman, L. and Pollak, M.: 1981, The Hall effect due to hopping conduction in the localized states of amorphous semiconductors, J. Phys. Coll. 42, Friedman, S. P. and Seaton, N. A.: 1998, Critical path analysis of the relationship between permeability and electrical conductivity of three-dimensional pore networks, Water Resour. Res. 34, Hazlett, R. D. and Furr, M. J.: 2000, Percolation model for permeability reduction in porous media by continuous-gas foams, Ind. Eng. Chem. Res. 39, He, Z., Dagan, Z. and Maldarelli, C.: 1991, The influence of surfactant adsorption on the motion of a fluid sphere in a tube. Part 1. Uniform retardation controlled by sorption kinetics, J. Fluid Mech. 222, Hunt, A. G.: 2001, Applications of percolation theory to porous media with distributed local conductances, Adv. Water Res. 24, Hunt, A. G. and Gee, G. W.: 2002, Application of critical path analysis to fractal porous media: Comparison with examples from the Manford site, Adv. Water Res. 25, Katz, A. J. and Thompson, A. H.: 1986, Quantitative prediction of permeability in porous rock, Phys. Rev. B 34, Kirkpatrick, S.: 1973a, Percolation and conduction, Rev. Mod. Phys. 45, 574. Kirkpatrick, S.: 1973b, Hopping conduction: experiment versus theory, In: J. Stuke and W. Brenig (eds), Proceedings of the 5th Annual Conference on Amorphous and Liquid Semiconductors, Taylor & Francis, London, pp Moldrup, P., Olesen, T., Komatsu, T., Schjonning, P. and Roleston, D.: 2001, Tortuosity, diffusivity, and permeability in the soil liquid and gaseous phases, Soil Sci. Soc. Am. J. 65,
15 HYDRAULIC CONDUCTIVITY OF POROUS MATERIALS 65 Olbricht, W. L.: 1996, Pore-scale prototypes of multiphase flow in porous media, Ann. Rev. Fluid Mech. 28, Patzek, T. W.: 2001, Verification of a complete pore network simulator of drainage and imbibition, SPE Journal 6, Perrier, E., Rieu, M., Sposito, G. and DeMarsily, G.: 1996, Models of the water retention curve for soils with a fractal pore size distribution, Water Resour. Res. 32, Pollak, M.: 1972, A percolation treatment of dc hopping conduction, J. Non-Cryst. Solids 11, Pollak, M.: 1985, Disordered Semiconductors, Chap. 5b (Other transport properties), CRC Press, Boca Raton, FL. Pollak, M. and Riess, I.: 1976, A percolation treatment of high-field hopping transport, J. Phys. C 9, Ratulowski, J. and Chang, H.-C: 1989, Transport of gas bubbles in capillaries, Phys. Fluids A 1, Ratulowski, J. and Chang, H.-C.: 1990, Marangoni effects of trace impurities on the motion of long gas bubbles in capillaries, J. Fluid Mech. 210, Rieu, M. and Sposito, G.: 1991, Fractal fragmentation, soil porosity, and soil water properties. I. Theory, Soil Sci. Soc. Am. J. 55, Rossen, W. R.: 1990, Theory of mobilization pressure gradient of flowing foams in porous media, J. Colloid Interf. Sci. 136, Rossen, W. R., Shi, J. and Zeilinger, S. C.: 1994, Percolation modeling of foam generation in porous media, AIChE J. 40, Sahimi, M.: 1995, Flow and Transport in Porous Media and Fractured Rock: from Classical Methods to Modern Approaches, VCH, Weinheim, Federal Republic of Germany, Seager, C. H. and Pike, G. E.: 1974, Percolation and conductivity: A computer study. Phys. Rev. B 10, , Shklovskii, B. I.: 1976, Nonohmic hopping conduction, Sov. Phys. Semicond. 10, Shklovskii, B. I. and Efros, A. L.: 1971, Zh. Eksp. Teor. Fiz. 60, 867. Skal, A. S.: 2001, Superconductive percolation: transmission of order by two kinds of poles, zeros and critical fluctuations, Phil. Mag. B 81, Stark, J. and Manga, M.: 2000, The motion of long bubbles in a network of tubes, Trans. Porous Med. 40, Stauffer, D.: 1975, Introduction to Percolation Theory, Taylor & Francis, London. Turcotte, D. L.: 1986, Fractals and fragmentation, J. Geophys. Res. 91, Wu, Q., Borkovec, M. and Sticher, H.: 1993, On particle-size distributions in soils, Soil Sci. Am. J. 57,
Determination of permeability using electrical properties of reservoir rocks by the critical path analysis
World Essays Journal / 3 (2): 46-52, 2015 2015 Available online at www. worldessaysj.com Determination of permeability using electrical properties of reservoir rocks by the critical path analysis Ezatallah
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationExcess 1/f noise in systems with an exponentially wide spectrum of resistances and dual universality of the percolation-like noise exponent
Excess 1/f noise in systems with an exponentially wide spectrum of resistances and dual universality of the percolation-like noise exponent A. A. Snarski a) Ukrainian National Technical University KPI,
More informationResistance distribution in the hopping percolation model
Resistance distribution in the hopping percolation model Yakov M. Strelniker, Shlomo Havlin, Richard Berkovits, and Aviad Frydman Minerva Center, Jack and Pearl Resnick Institute of Advanced Technology,
More informationModels of Disordered Media and Predictions of Associated Hydraulic Conductivity
Wright State University CORE Scholar Browse all Theses and Dissertations Theses and Dissertations 2006 Models of Disordered Media and Predictions of Associated Hydraulic Conductivity L. Aaron Blank Jr.
More informationTHE EFFECT OF WATER SATURATION ON GAS SLIP FACTOR BY PORE SCALE NETWORK MODELING
SCA00-53 1/6 THE EFFECT OF WATER SATURATION ON GAS SLIP FACTOR BY PORE SCALE NETWORK MODELING Liu Qingjie, Liu Baohua, Li Xianbing, Yan Shouguo Research Institute of Petroleum Exploration and Development,
More informationPore-scale modeling extension of constitutive relationships in the range of residual saturations
WATER RESOURCES RESEARCH, VOL. 37, NO. 1, PAGES 165 170, JANUARY 2001 Pore-scale modeling extension of constitutive relationships in the range of residual saturations Rudolf J. Held and Michael A. Celia
More informationConductance and resistance jumps in finite-size random resistor networks
J. Phys. A: Math. Gen. 21 (1988) L23-L29. Printed in the UK LETTER TO THE EDITOR Conductance and resistance jumps in finite-size random resistor networks G G Batrounit, B Kahngt and S Rednert t Department
More informationFracture relative permeability revisited
Fracture relative permeability revisited NOROLLAH KASIRI and GHASEM BASHIRI, Iran University of Science and Technology Relative permeability is one of the most uncertain terms in multiphase flow through
More informationWater retention of prefractal porous media generated with the homogeneous and heterogeneous algorithms
WATER RESOURCES RESEARCH, VOL. 37, NO. 10, PAGES 2631 2636, OCTOBER 2001 Water retention of prefractal porous media generated with the homogeneous and heterogeneous algorithms Michael C. Sukop 1 Department
More informationPHYSICS OF FLUID SPREADING ON ROUGH SURFACES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Supp, Pages 85 92 c 2008 Institute for Scientific Computing and Information PHYSICS OF FLUID SPREADING ON ROUGH SURFACES K. M. HAY AND
More informationDiffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad
Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in
More informationModified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media
J. Phys. A: Math. Gen. 3 (998) 7227 7234. Printed in the UK PII: S0305-4470(98)93976-2 Modified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media Juris Robert Kalnin
More informationLIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE
Proceedings of the ASME/JSME 2011 8th Thermal Engineering Joint Conference AJTEC2011 March 13-17, 2011, Honolulu, Hawaii, USA AJTEC2011-44190 LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Youngbae
More informationDarcy's Law. Laboratory 2 HWR 531/431
Darcy's Law Laboratory HWR 531/431-1 Introduction In 1856, Henry Darcy, a French hydraulic engineer, published a report in which he described a series of experiments he had performed in an attempt to quantify
More informationThe Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media
Transport in Porous Media 44: 109 122, 2001. c 2001 Kluwer Academic Publishers. Printed in the Netherls. 109 The Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media RAMON G.
More informationI. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010
: an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS
More informationFluid Flow Fluid Flow and Permeability
and Permeability 215 Viscosity describes the shear stresses that develop in a flowing fluid. V z Stationary Fluid Velocity Profile x Shear stress in the fluid is proportional to the fluid velocity gradient.
More informationPORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Igor Bogdanov 1*, Sylvain Jardel 1, Anis Turki 1, Arjan Kamp 1 1 Open &
More informationPractical methodology for inclusion of uplift and pore pressures in analysis of concrete dams
Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty
More informationTransactions on Ecology and the Environment vol 7, 1995 WIT Press, ISSN
Non-Fickian tracer dispersion in a karst conduit embedded in a porous medium A.P. Belov*, T.C. Atkinson* "Department of Mathematics and Computing, University of Glamorgan, Pontypridd, Mid Glamorgan, CF3
More informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationPartial Saturation Fluid Substitution with Partial Saturation
Fluid Substitution with 261 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) 4 3.5 K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage 3 2.5 2 Fine-scale mixing 1 = S w + S o + S g K fluid
More informationKOZENY-CARMAN EQUATION REVISITED. Jack Dvorkin Abstract
KOZENY-CARMAN EQUATION REVISITED Jack Dvorkin -- 009 Abstract The Kozeny-Carman equation is often presented as permeability versus porosity, grain size, and tortuosity. When it is used to estimate permeability
More informationAdvanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur. Lecture 6
Advanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur Lecture 6 Good morning and welcome to the next lecture of this video course on Advanced Hydrology.
More informationPETROLEUM RESERVOIRS FLUID FLOW IN. ill OR 236 URBANA X Q ~ < o S z» 5 8. DIVISION OF THE ILLINOIS STATE GEOLOGICAL SURVEY JOHN C.
s 14.GS: OR 236 c. 1 ILLINOIS GEOLOGICAL SURVEY LIBRARY STATE OF ILLINOIS WILLIAM G. STRATTON, Governor DEPARTMENT OF REGISTRATION AND EDUCATION VERA M. BINKS, Director FLUID FLOW IN PETROLEUM RESERVOIRS
More informationSaturation-Dependence of Dispersion in Porous Media
Wright State University CORE Scholar Physics Faculty Publications Physics 12-2012 Saturation-Dependence of Dispersion in Porous Media B. Ghanbarian-Alavijeh Wright State University - Main Campus Thomas
More informationAbility of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium Situations
Research Article imedpub Journals http://www.imedpub.com Resources, Recycling and Waste Management Ability of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationTemperature ( o C)
Viscosity (Pa sec) Supplementary Information 10 8 10 6 10 4 10 2 150 200 250 300 Temperature ( o C) Supplementary Figure 1 Viscosity of fibre components (PC cladding blue; As 2 Se 5 red; CPE black) as
More informationSimulating Fluid-Fluid Interfacial Area
Simulating Fluid-Fluid Interfacial Area revealed by a pore-network model V. Joekar-Niasar S. M. Hassanizadeh Utrecht University, The Netherlands July 22, 2009 Outline 1 What s a Porous medium 2 Intro to
More informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
More informationSimplified thermodynamic model for equilibrium capillary pressure in a fractal porous medium
Simplified thermodynamic model for equilibrium capillary pressure in a fractal porous medium M. R. Deinert Theoretical and Applied Mechanics, Cornell University, Ithaca, New York 14850, USA J.-Y. Parlange
More informationCorrelation Between Resistivity Index, Capillary Pressure and Relative Permeability
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Correlation Between Resistivity Index, Capillary Pressure and Kewen Li Stanford Geothermal Program, Stanford University, Stanford,
More informationPROCEEDINGS THIRD WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING. December 14-15,1977
SGPTR258 PROCEEDINGS THIRD WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING December 1415,1977 *Conducted under Subcontract No. 16735 with Lawrence Berkeley Laboratory, Universityof California, sponsored by the
More informationFrom Solute Transport to Chemical Weathering
Wright State University CORE Scholar Physics Seminars Physics --203 From Solute Transport to Chemical Weathering Allen Hunt Wright State University - Main Campus, allen.hunt@wright.edu Thomas E. Skinner
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few
More informationAn approximate analytical solution for non-darcy flow toward a well in fractured media
WATER RESOURCES RESEARCH, VOL. 38, NO. 3, 1023, 10.1029/2001WR000713, 2002 An approximate analytical solution for non-arcy flow toward a well in fractured media Yu-Shu Wu Earth Sciences ivision, Lawrence
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationLiquid Jet Breakup at Low Weber Number: A Survey
International Journal of Engineering Research and Technology. ISSN 0974-3154 Volume 6, Number 6 (2013), pp. 727-732 International Research Publication House http://www.irphouse.com Liquid Jet Breakup at
More informationMeasurement of Liquid Film Thickness in Micro Square Channel
Measurement of Liquid Film Thickness in Micro Square Channel Youngbae Han and Naoki Shikazono Department of Mechanical Engineering, The University of Tokyo Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan
More informationTHE MISSING LINK BETWEEN PORE-SCALE ANCHORING AND PORE-SCALE PREDICTION
SCA22-25 /4 THE MISSING LINK BETWEEN PORE-SCALE ANCHORING AND PORE-SCALE PREDICTION Steven R. McDougall and Ken S. Sorbie Department of Petroleum Engineering, Heriot-Watt University ABSTRACT An earlier
More informationQuarterly Report for Contract DE-FG36-02ID14418 Stanford Geothermal Program October-December 2004
Quarterly Report for Contract DE-FG36-02ID14418 Stanford Geothermal Program October-December 2004 2 Table of Contents 1. THEORETICAL STUDY OF PHASE TRANSFORMATION EFFECTS ON STEAM-WATER RELATIVE PERMEABILITIES
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationMiao Boya and An Yu Department of Physics, Tsinghua University, Beijing , People s Republic of China
Localization in an acoustic cavitation cloud Miao Boya and An Yu Department of Physics, Tsinghua University, Beijing 100084, People s Republic of China Using a nonlinear sound wave equation for a bubbly
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationSatiated relative permeability of variable-aperture fractures
Satiated relative permeability of variable-aperture fractures Russell L. Detwiler Lawrence Livermore National Laboratory, University of California, 7000 East Avenue, Livermore, California 94551, USA Harihar
More informationMigrations of Fines in Porous Media
Migrations of Fines in Porous Media by Kartic C. Khilar Department of Chemical Engineering, Indian Institute of Technology, Bombay, India and H. Scott Fogler Department of Chemical Engineering, University
More informationMEASUREMENT OF CAPILLARY PRESSURE BY DIRECT VISUALIZATION OF A CENTRIFUGE EXPERIMENT
MEASUREMENT OF CAPILLARY PRESSURE BY DIRECT VISUALIZATION OF A CENTRIFUGE EXPERIMENT Osamah A. Al-Omair and Richard L. Christiansen Petroleum Engineering Department, Colorado School of Mines ABSTRACT A
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More information2nd International Conference Mechanics of Unsaturated Soils 7 th 9 th March 2007
2nd International Conference Determination of the Soil Water Retention Curve and the Unsaturated Hydraulic Conductivity from the Particle Size Distribution Alexander Scheuermann & Andreas Bieberstein Motivation
More informationA two-dimensional network simulator for two-phase flow in porous media
A two-dimensional network simulator for two-phase flow in porous media Eyvind Aker, Knut Jørgen Måløy Department of Physics University of Oslo, N-0316 Oslo, Norway Alex Hansen Department of Physics, Norwegian
More informationPaper E A pore network model for calculation of interfacial velocities
Paper E A pore network model for calculation of interfacial velocities Submitted to Advances in Water Resources, fall 21. A PORE NETWORK MODEL FOR CALCULATION OF INTERFACIAL VELOCITIES H.F. Nordhaug a,
More informationCORRELATION OF THE MULTIPHASE FLOW COEFFICIENTS OF POROUS MEDIA WITH WETTABILITY: A PORE NETWORK APPROACH
CMWRXVI CORRELATION OF THE MULTIPHASE FLOW COEFFICIENTS OF POROUS MEDIA WITH WETTABILITY: A PORE NETWORK APPROACH CHRISTOS D. TSAKIROGLOU Institute of Chemical Engineering and High Temperature Chemical
More informationarxiv: v1 [physics.geo-ph] 25 Mar 2010
Small-World Networks in Geophysics arxiv:03.4886v1 [physics.geo-ph] 5 Mar 0 Xin-She Yang Department of Applied Mathematics, University of Leeds, LEEDS LS 9JT, UK Abstract Many geophysical processes can
More informationInfluence of the Flow Direction on the Mass Transport of Vapors Through Membranes Consisting of Several Layers
Influence of the Flow Direction on the Mass Transport of Vapors Through Membranes Consisting of Several Layers Thomas Loimer 1,a and Petr Uchytil 2,b 1 Institute of Fluid Mechanics and Heat Transfer, Vienna
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More information2. Modeling of shrinkage during first drying period
2. Modeling of shrinkage during first drying period In this chapter we propose and develop a mathematical model of to describe nonuniform shrinkage of porous medium during drying starting with several
More informationPermeability and fluid transport
Permeability and fluid transport Thermal transport: Fluid transport: q = " k # $p with specific discharge (filter velocity) q [m s 1 ] pressure gradient p [N m 3 ] dynamic viscosity η [N s m 2 ] (intrinsic)
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationUse of Fractal Geometry for Determination of Pore Scale Rock Heterogeneity
Use of Fractal Geometry for Determination of Pore Scale Rock Heterogeneity Summary Dipak Mandal, DC Tewari, MS Rautela, TR Misra Institute of Reservoir Studies, ONGC, Chandkheda Campus, Ahmedabad Fractal
More informationDarcy s Law, Richards Equation, and Green-Ampt Equation
Darcy s Law, Richards Equation, and Green-Ampt Equation 1. Darcy s Law Fluid potential: in classic hydraulics, the fluid potential M is stated in terms of Bernoulli Equation (1.1) P, pressure, [F L!2 ]
More informationINJECTION, CONDUCTION AND PRODUCTION
Chapter III Injection, Conduction and Production Chapter III From a physical point of view strictly steady-state conditions of heterogeneous-fluid flow in oil-producing systems are virtually never encountered.
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationGravitational effects on the deformation of a droplet adhering to a horizontal solid surface in shear flow
PHYSICS OF FLUIDS 19, 122105 2007 Gravitational effects on the deformation of a droplet adhering to a horizontal solid surface in shear flow P. Dimitrakopoulos Department of Chemical and Biomolecular Engineering,
More informationA HYBRID SEMI-ANALYTICAL AND NUMERICAL METHOD FOR MODELING WELLBORE HEAT TRANSMISSION
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 5 SGP-TR-176 A HYBRID SEMI-ANALYTICAL AND NUMERICAL METHOD FOR MODELING
More informationLinear stability of radial displacements in porous media: Influence of velocity-induced dispersion and concentration-dependent diffusion
PHYSICS OF FLUIDS VOLUME 16, NUMBER 10 OCTOBER 2004 Linear stability of radial displacements in porous media: Influence of velocity-induced dispersion and concentration-dependent diffusion A. Riaz Department
More informationComparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media
Comparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media Charlie Matsubara *, Tim Kuo, Helen Wu PIM Consulting, No. 16, Gonyequ 1 st
More informationAn Adsorption Desorption-Controlled Surfactant on a Deforming Droplet
JOURNAL OF COLLOID AND INTERFACE SCIENCE 208, 68 80 (1998) ARTICLE NO. CS985816 An Adsorption Desorption-Controlled Surfactant on a Deforming Droplet Charles D. Eggleton* and Kathleen J. Stebe,1 *Department
More informationOnset of convection of a reacting fluid layer in a porous medium with temperature-dependent heat source
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2011, Science Huβ, http://www.scihub.org/ajsir ISSN: 2153-649X doi:10.5251/ajsir.2011.2.6.860.864 Onset of convection of a reacting fluid layer in
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
More informationSteady-state Two-phase Flow in Porous Media: Statistics and Transport Properties
Steady-state Two-phase Flow in Porous Media: Statistics and Transport Properties Ken Tore Tallakstad, 1 Henning Arendt Knudsen, 1 Thomas Ramstad, 2,3 Grunde Løvoll, 1 Knut Jørgen Måløy, 1 Renaud Toussaint,
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationSupporting Information
Supporting Information Retention and Release of Graphene Oxide in Structured Heterogeneous Porous Media under Saturated and Unsaturated Conditions Shunan Dong 1, Xiaoqing Shi 1, Bin Gao 3, Jianfeng Wu
More informationSunday September 28th. Time. 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn. CP: Contributed Presentation 25 mn
Sunday September 28th 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn CP: Contributed Presentation 25 mn Monday September 29th 07:30 am 08:30 am Registration 08:30 am 09:00
More informationPercolation for a model of statistically inhomogeneous random media
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 13 1 OCTOBER 1999 Percolation for a model of statistically inhomogeneous random media J. Quintanilla a) Department of Mathematics, University of North Texas,
More informationModeling of 1D Anomalous Diffusion In Fractured Nanoporous Media
LowPerm2015 Colorado School of Mines Low Permeability Media and Nanoporous Materials from Characterisation to Modelling: Can We Do It Better? IFPEN / Rueil-Malmaison - 9-11 June 2015 CSM Modeling of 1D
More informationEffect of Sorption/Curved Interface Thermodynamics on Pressure transient
PROCEEDINGS, Twentieth Workshop on Geothermal Rey~volr Englneerlng Stanford Unhrenlty, Stanfoni, Callfornla, January 2426 1995 SGP-m-150 Effect of Sorption/Curved Interface Thermodynamics on Pressure transient
More informationThe Impact of Pore-Scale Flow Regimes on Upscaling of Immiscible Two-Phase Flow in Geothermal Reservoirs
PROCEEDINGS 43rd Workshop on Geothermal Reservoir Engineering Stanford University Stanford California February 12-14 2018 SGP-TR-213 The Impact of Pore-Scale Flow Regimes on Upscaling of Immiscible Two-Phase
More informationSupporting Information. Technique for real-time measurements of endothelial permeability in a
Supporting Information Technique for real-time measurements of endothelial permeability in a microfluidic membrane chip using laser-induced fluorescence detection Edmond W.K. Young a,b,, Michael W.L. Watson
More information1. Water in Soils: Infiltration and Redistribution
Contents 1 Water in Soils: Infiltration and Redistribution 1 1a Material Properties of Soil..................... 2 1b Soil Water Flow........................... 4 i Incorporating K - θ and ψ - θ Relations
More informationA LABORATORY STUDY OF FOAM FOR EOR IN NATURALLY FRACTURED RESERVOIRS. William R. Rossen Bander. I. AlQuaimi
A LABORATORY STUDY OF FOAM FOR EOR IN NATURALLY FRACTURED RESERVOIRS William R. Rossen Bander. I. AlQuaimi Gravity Backround Gas-injection EOR can displace nearly all oil contacted, but sweep efficiency
More informationAnalysis of a drainage efficiency in stratified porous media
Analysis of a drainage efficiency in stratified porous media CHAKIB SELADJI Département de Génie Mécanique Université Abou Bekr BELKAID BP 119 Chétouane Tlemcen 13000 ALGERIA seladji@yahoo.fr Abstract:
More informationMohamed Daoud Claudine E.Williams Editors. Soft Matter Physics. With 177 Figures, 16 of them in colour
Mohamed Daoud Claudine E.Williams Editors Soft Matter Physics With 177 Figures, 16 of them in colour Contents 1. Droplets: CapiUarity and Wetting 1 By F. Brochard-Wyart (With 35 figures) 1.1 Introduction
More informationStochastic geometry and porous media
Stochastic geometry and transport in porous media Hans R. Künsch Seminar für Statistik, ETH Zürich February 15, 2007, Reisensburg Coauthors Thanks to the coauthors of this paper: P. Lehmann, A. Kaestner,
More informationSOLIDS AND LIQUIDS - Here's a brief review of the atomic picture or gases, liquids, and solids GASES
30 SOLIDS AND LIQUIDS - Here's a brief review of the atomic picture or gases, liquids, and solids GASES * Gas molecules are small compared to the space between them. * Gas molecules move in straight lines
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationSCA : A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP
SCA2003-53: A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP This paper was prepared for presentation at the International
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More information1 Modeling Immiscible Fluid Flow in Porous Media
Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.
More informationPHYSICAL REALITIES FOR IN DEPTH PROFILE MODIFICATION. RANDY SERIGHT, New Mexico Tech
PHYSICAL REALITIES FOR IN DEPTH PROFILE MODIFICATION RANDY SERIGHT, New Mexico Tech 1. Gel treatments (of any kind) are not polymer floods. 2. Crossflow makes gel placement challenging. 3. Adsorbed polymers,
More informationSTATISTICAL COMPARISON OF ELECTRICAL EFFICIENCY THEORY TO ARCHIE S EQUATIONS AND EFFECTIVE MEDIUM THEORY by Charles R.
STATISTICAL COMPARISON OF ELECTRICAL EFFICIENCY THEORY TO ARCHIE S EQUATIONS AND EFFECTIVE MEDIUM THEORY 21-1 by Charles R. Berg ABSTRACT A comparison is made between Archie, electrical efficiency theory,
More informationApplication of Green Function Method for Solution Advection Diffusion Equation of Nutrient Concentration in Groundwater
Annals of Pure and Applied Mathematics Vol. 7, No. 2, 2014, 1-5 ISSN: 2279-087X (P), 2279-0888(online) Published on 26 August 2014 www.researchmathsci.org Annals of Application of Green Function Method
More informationDRAINAGE AND IMBIBITION CAPILLARY PRESSURE CURVES OF CARBONATE RESERVOIR ROCKS BY DIGITAL ROCK PHYSICS
SCA2012-56 1/6 DRAINAGE AND IMBIBITION CAPILLARY PRESSURE CURVES OF CARBONATE RESERVOIR ROCKS BY DIGITAL ROCK PHYSICS Y. Mu 1, Q. Fang 1, C. Baldwin 1, J. Töelke 1, A. Grader 1, M. Dernaika 2, and Z. Kalam
More informationEssay 3. Mathematical Background of Heat Conduction
Essay 3 Mathematical Background of Heat Conduction 3.1 Fourier s Law of Heat Conduction In a groundbreaking treatise published in 1822, Théorie Analytique de la Chaleur, Jean Baptiste Joseph Fourier set
More informationREVIEW OF THE WINLAND R35 METHOD FOR NET PAY DEFINITION AND ITS APPLICATION IN LOW PERMEABILITY SANDS
REVIEW OF THE WINLAND R35 METHOD FOR NET PAY DEFINITION AND ITS APPLICATION IN LOW PERMEABILITY SANDS Mike Spearing, Tim Allen and Gavin McAulay (AEA Technology) INTRODUCTION The definition of net sand
More informationError in the Estimation of Effective Diffusion Coefficients from Sorption Measurements*
Error in the Estimation of Effective Diffusion Coefficients from Sorption Measurements* D. BOBOK and E. BESEDOVÁ Department of Chemical and Biochemical Engineering, Faculty of Chemical Technology, Slovak
More informationGeneralised Separable Solution of Double Phase Flow through Homogeneous Porous Medium in Vertical Downward Direction Due to Difference in Viscosity
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 932-9466 Vol. 8, Issue (June 203), pp. 305-37 Applications and Applied Mathematics: An International Journal (AAM) Generalised Separable Solution
More information