Advances in Colloid and Interface Science 82 Ž1999. 127᎐187 Foam in porous media: thermodynamic and hydrodynamic peculiarities Konstantin G. Kornev a , Alexander V. Neimark b,U , Aleksey N. Rozhkov a a Institute for Problems in Mechanics, Russian Academy of Sciences, Prospect Vernadskogo 101, Moscow 117526, Russia b TRI r Princeton, 601 Prospect A¨ e., Princeton, NJ 08542-0625, USA Abstract Thermodynamic and hydrodynamic properties of foams in porous media are examined from a unified point of view. We show that interactions between foam films Žlamellae. and wetting films covering the pore walls play an important role in treating experimental data and constructing a general theory of foam residence and motion through porous media. Mechanisms of in situ bubble generation, foam patterning, and rheological peculiarities of foams in pores are discussed in detail. In particular, we clarify the difference between foam lamellae and liquid lenses, focusing on intermolecular forces in thin foam and wetting films. A consistent description of conditions of mechanical equilibrium of curved lamellae, including dynamic effects, is presented for the first time. This microlevel approach enables us to describe the dependence of the capillary pressure in the Plateau border on the current state of the pair ‘wetting film᎐foam lamella’. We review a theory of foam patterning under a load. Two driving forces are invoked to explain specific interactions between the solid skeleton and foams. The binding forces caused by bubble compressibility and the pinning forces due to capillarity determine a specific ordering of lamellae in porous media. The microscopic bubble train model predicts asymptotic expressions for the start-up-yield pressure drop. We consider key problems that underlay the understanding of physical mechanisms of anomalous foam resistance. Different micromechanical models of foam friction are U Corresponding author. Tel.: q1-609-430-4818; fax: q1-609-683-7149; e-mail: aneimark@triprinceton.org 0001-8686r99r$ - see front matter ᮊ 1999 Elsevier Science B.V. All rights reserved. PII: S 0 0 0 1 - 8 6 8 6 Ž 9 9 . 0 0 0 1 3 - 5 128 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 thoroughly discussed. Bretherton’s Ž1961. theory of the forced, steady fluid᎐fluid displacement is reviewed in application to bubble transport through pore channels. The origins of disagreement of the theory and experiment are discussed. The Bretherton theory is augmented based on a new sailboat model, which accounts for thermodynamic coupling of foam lamellae and wetting films. Special attention is paid to studies of stick-slip motion of lamellae and bubbles in pores of varying diameter. Finally, we discuss macroscale models and analyze topical problems of foam behavior in porous media, including reservoirs, granular, and fibrous materials. ᮊ 1999 Elsevier Science B.V. All rights reserved. Keywords: Foam; Lamella; Wetting films; Porous media; Foam patterning and flow Contents 1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2. Thermodynamic peculiarities of foam distribution in confining systems . . . . . . . 2.1. Bulk foam and foam films: capillary and disjoining pressures . . . . . . . . . . 2.2. Foam patterning in porous media . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.3. Wetting films and menisci in pores: limiting capillary pressure . . . . . . . . . 2.4. Mechanisms of bubble generation in porous media: gas path closure . . . . . 2.5. Transformation of lens into lamella . . . . . . . . . . . . . . . . . . . . . . . . . . 2.6. Lamellae in convergent Ždivergent. pores . . . . . . . . . . . . . . . . . . . . . . . 2.7. Chain of lamellae: correlation length . . . . . . . . . . . . . . . . . . . . . . . . . 2.8. Superstructures: analysis of the phase portrait . . . . . . . . . . . . . . . . . . . . 2.9. Start-up pressure drop . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3. Foam transport in smooth capillaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.1. The Bretherton mechanism of gas mobility reduction: motion of individual bubbles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.2. Modifications of the Bretherton theory: bubble train motion . . . . . . . . . . 3.3. Sailboat model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.3.1. Continuity condition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.3.2. Force balance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.3.3. Lamella thickness and dynamic tension . . . . . . . . . . . . . . . . . . . . 3.3.4. Pressure continuity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.3.5. Mass balance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.4. Sailboat model: traveling wave solution . . . . . . . . . . . . . . . . . . . . . . . . 4. Flow trough bamboo-like capillaries and porous media . . . . . . . . . . . . . . . . . 4.1. Lamella in bamboo-like capillaries: stick-slip motion . . . . . . . . . . . . . . . 4.2. Stick-slip motion of bubble trains . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4.3. Weak foams: flow of ‘solutions’ of bubble chains through porous media . . . 4.4. Problems in describing strong foams . . . . . . . . . . . . . . . . . . . . . . . . . . 4.5. Foams in fiber systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5. Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129 . 129 . 129 . 134 . 136 . 138 . 143 . 146 . 148 . 151 . 153 . 157 . . . . . . . . . . . . . . . . . . 158 . 160 . 163 . 164 . 164 . 165 . 166 . 167 . 167 . 169 . 169 . 173 . 175 . 179 . 180 . 181 . 182 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 129 1. Introduction Foams display a variety of behaviors that distinguish them into the special subject of interface science. The cell structure of foams, exhibiting foam irregularity, instability, chaos and the like, can be observed daily. The specific properties of foams have attracted the attention of researches for centuries. At the present time, a substantial collection of books and reviews w1᎐12x, cover a wide range of fundamental questions. Especially when confined in porous media, foams have a unique hierarchical structure and rheological properties. Novel applications of foams in porous media have become an intriguing subject of fundamental, systematic studies. In particular, because of their extremely efficient blocking action, foams assist groundwaterrsoil remediation, selective blockage and extraction of oil or other liquids from porous materials. Industrial foam technologies include dyeing and finishing of textile fabrics, paper coating, resin-impregnation of fibrous mats and fabrics. However, while technologically effective, foams in porous media remain enigmatic. In this review, we discuss selected problems associated with the thermodynamics and hydrodynamics of foams in porous media. Earlier reviews w13᎐17x are useful guides to this branch of studies. We focus on the physical mechanisms governing foam flow and residence in confined geometries. Mechanistic models of foam evolution in bulk and in porous media, which are based on the population balance method w18,19x, dynamic scaling w20᎐22x, and also percolation models of foam flow in pore networks Žsee the latest review of Rossen w17x and references therein . are beyond the scope of this survey. 2. Thermodynamic peculiarities of foam distribution in confining systems 2.1. Bulk foam and foam films: capillary and disjoining pressures Prior to analyzing the specifics of foam behavior in porous media, let us briefly review the properties of bulk foam. Bulk foam is often encountered in everyday life Žsoap foam. and nature Žsea froth., and finds numerous technological applications including: fire fighting, surface and in situ cleaning, separations, etc. Traditional and potential applications of foam in engineering processes are reviewed by several authors w3,8,9,11x. Bulk foams are dispersions of gas bubbles in a liquid matrix. We distinguish between two types of foams. The first is the so-called quasi-spherical foam, a system of nearly spherical bubbles dispersed in a continuos liquid phase. This foam is similar to an emulsion. The second type, the focus of this review, is the so-called polyhedral foam, a system of coarsened bubbles separated by thin liquid films, usually called lamellae. The bubbles in polyhedral foams have a distorted polyhedral shape and may differ in their size and number of neighbors, which equals the number of polyhedral faces formed by the lamellae. The lamellae separating non-equal bubbles may be curved due to pressure differences between the bubbles. Real foams are unstable, and evolve because a decrease in the lamella area leads 130 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 to a decrease in the interfacial free energy. The physical factors controlling foam evolution, patterning and stability are gravity-driven liquid drainage, interbubble gas diffusion caused by the pressure difference across the lamellae, and lamella thinning and rupture. To guarantee a long lifetime, the lamellae are stabilized by surface active molecules or surfactants, such as sodium lauryl sulfate and the like Žwhich reduce the interfacial tension and establish an additional repulsion between the lamella surfaces ., and additives, e.g. hydroxylethylcellulose, to increase the liquid viscosity. Even in the ideal situation of an initially perfectly ordered foam, comprised of uniform bubbles with a uniform gas pressure, finite-size perturbations in bubble shape would lead to an irreversible growth of the larger bubbles at the expense of the smaller. This effect, known as coarsening w3x is unavoidable, yet can be slowed by using surfactants and additives. Although foam, as a non-equilibrium thermodynamic system, evolves in time, the instantaneous foam configuration may be regarded as a quasi-equilibrium configuration for the current distribution of pressures. The characteristic time to establish a local mechanical equilibrium between the bubbles and the lamellae is much smaller than the characteristic time for pressure changes due to mass transfer. Thus, the current distribution of pressures in the bubbles determines foam geometry. Although global foam structure is difficult to quantify, mechanical equilibrium necessarily minimizes the surface free energy with respect to infinitesimal changes in the local configuration of bubbles and lamellae, determining the local geometry, or microgeometry of foams. Foam microgeometry was studied in detail by JosefAntoine-Ferdinand Plateau over the course of a quarter century Ž1843᎐1869.; for a contemporary description see Nitsohe w23x. The Plateau law, one of the most elegant physical laws of nature, can be formulated as follows: in ‘dry’ foams with a small liquid content, the lamellae between the bubbles form a polyhedral cell structure with certain universal features. The faces Žlamellae. always meet three at a time at equal angles of 120Њ and form the edges named Plateau borders. The edges ŽPlateau borders. always meet each other in groups of four at the equal tetrahedral or Maraldi angle of 109Њ28Љ16Ј, whose cosine equals y1r3. The four-way junctions of Plateau borders are usually referred to as vertices. The Plateau borders constitute a network of channels for liquid transport through the foam. The vertices can be regarded as the nodes of this network, and the Plateau borders as the bonds. These specific features of the foam microgeometry determine to a great extent the rheological, mechanical and transport properties of foams. The Plateau law expresses the general condition of mechanical stability, provided that the lamellae are considered as hypothetical membranes with an inherent uniform surface tension. The validity of the Plateau law has been confirmed for various cellular systems including foams, biological tissues, polycrystals, and magnetic domains w6,10,23᎐26x. The Plateau law is not limited to cells with flat interfaces. but is valid for any foam, with curved lamellae between bubbles of different pressure. For curved interfaces, the inclination angles are measured between the planes tangent to the lamellae. Bolton and Weaire w27,28x have proven K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 131 the so-called decoration lemma and have shown that the Plateau law should be valid even in the case of thick Plateau borders, provided that the angles are measured using the perpendicular cross-section of the Plateau border between the extensions of the curved lamellae meeting in the center of the border. In reality, the films between the bubbles tend to squeeze out liquid. Consider a single Plateau border in ideal uniform foam, assuming that all the lamellae are flat ŽFig. 1.. The curvature of the Plateau border walls causes a pressure difference across the walls, equal to the Laplacian capillary pressure, Pc s Pg y Pw s ž 1 r1 q 1 r2 / , Ž1. where is the surface tension, ri are the principal radii of curvature of a meniscus, Pg is the gas pressure and Pw is the hydrostatic pressure in the liquid. By virtue of the Laplacian capillary pressure Pc , the hydrostatic pressure within the edge will be reduced with respect to the gas pressure. On the other hand, we would expect that the hydrostatic pressure within the flat lamellae would be equal to the gas pressure. At first sight, even in perfectly ordered foam with uniform gas pressure, a pressure difference exists between the Plateau borders and the lamellae, and the liquid should be squeezed out of the lamellae into the Plateau borders. However, in stabilized foams, the local mechanical equilibrium in the liquid phase is ensured by the additional, disjoining, pressure acting in the lamellae, which counterbalances the difference between the gas and hydrostatic pressures. The disjoining pressure, introduced by Derjaguin w29x, plays a key role in thin film stability. Here, it should be remembered that the foam cannot be formed from an ordinary pure liquid, i.e. without a so-called foaming agent. The ‘mother fluid’ for foams, bubbles and lamellae under consideration, is practically restricted to surfac- Fig. 1. Section through the Plateau border. Pl , pressure in lamellae; Pg ,, gas pressure; Pp b , pressure in the Plateau border. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 132 tant solutions. Typically, the surfactant molecules have two dissimilar parts ᎏ one hydrophilic, tending to be surrounded by water; the other, hydrophobic, tending to be oozed out by the attractive forces between the water molecules. The surface active molecules concentrate at the lamella surfaces so that their hydrophilic heads are in water, while their hydrophobic tails reside outside ŽFig. 2.. Drainage should occur until the lamella outer surfaces begin to ‘feel’ one another. For ionic surfactants, the electrostatic repulsion between the lamella surfaces prevents lamella bursting. Some other additives, either non-ionic or ionic Žorganic compounds, polymers, proteins, etc.., can also protect lamellae against bursting. The chemical and physical mechanisms of lamella stabilization have been intensively discussed w11,12,30᎐34x. The concept of disjoining pressure, which effectively accounts for the interactions between the lamella surfaces, allows us to quantitatively formulate the conditions of film equilibrium. Thermodynamic and experimental foundations of the existence of disjoining pressure in thin films were established almost 60 years ago in the pioneering papers of Derjaguin and Kussakov w35x and Derjaguin w36x. Disjoining pressure can stabilize perfectly ordered foam against small perturbations. Indeed, in thermodynamic equilibrium, the hydrostatic pressure within the lamellae equals that in the Plateau borders. Therefore, the capillary pressure is counterbalanced by the disjoining pressure, ⌸ Ž h., which is a function of the lamellae thickness, h: Pc s ⌸ Ž h . . Ž2. Disjoining pressure can be either positive Ždisjoining. or negative Žconjoining.. The magnitude and the sign of this excess pressure depends upon the physico᎐chemical properties of the film and its environment. In particular, in wetting films, the disjoining pressure is affected by the properties of both pairs: the ‘surrounding fluid-film’ and the ‘film-substrate’. The three characteristic types of disjoining pressure isotherms ŽFig. 3. correspond to stable Ž1., metastable Ž2., and unstable Ž3. films. The second type, with two descending branches, is also common for foam films ŽFig. 4.. The first descending branch of the disjoining pressure isotherm corresponds to the Newton black films of the smallest thickness. The Fig. 2. Schematic illustration of the surfactant distribution over the lamellar surfaces. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 133 Fig. 3. A schematic illustration of typical disjoining pressure isotherms: Ž1. stable films; Ž2. metastable films; Ž3. unstable films Žarbitrary units.. second descending branch corresponds to common black films w3,8,9,30x. The disjoining pressure isotherms vary with the type and concentration of the surfactant and electrolyte, or other foaming agents. The disjoining pressure isotherm of foam films intersects the h-axis when the film thickness exceeds a certain critical value ŽFig. 4., i.e. equilibrium lamellae exist only within a narrow range of the film thickness. Above a critical capillary pressure, the lifetimes of the lamellae and corresponding bulk foam become exceedingly short. Kruglyakov and Exerowa w8,9x analyzed the lifetime of bulk foam and individual foam lamellae with respect to the difference between the gas and liquid-phase pressures, i.e. in comparison to the capillary pressure. For a wide spectrum of experimental data, the thermodynamic properties of individual lamellae influence the lifetime of bulk foam. ŽFor foams in porous media, a similar analysis has been performed by Aronson et al. w40x.. The capillary pressure should correlate with the disjoining pressure in the lamella. Capillary pressure forces the lamella surfaces to come closer to one another. More stable foam results from foaming agents, which are able to form the Newton black films. Table 1 shows experimental data by Khristov et al. w41x, for single films and bulk foam made from 0.001 M sodium dodecyl sulfate with a sodium chloride additive. The critical micelle concentration of the solution corresponds to 0.18 M NaCl w32x. If the initial lamella is a thick film, liquid depletion increases the capillary suction pressure on the lamella. The film thins to a critical thickness near the primary maximum and spontaneously jumps to the Newton film ŽFig. 4.. For sufficiently high capillary suction pressure, even along the stable branches of ⌸ Ž h., macroscopic disturbances may initiate rupture. Despite the long history of the transition from common black films to Newton black films, a standard explanation is lacking w2᎐4,8,9,30x. 134 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 4. Experimental disjoining pressure isotherm at room temperature for sodium dodecyl sulfate ŽSDS. in brine ŽNaCI. w37᎐39x. Thus, foam stability is highly sensitive to the thermodynamics of the lamellae, which also affects the topological structure of real foams. In particular, for Newton black films, the deviations of vertex angles from 120Њ can reach several degrees w42x. 2.2. Foam patterning in porous media Foam patterning in confining systems, such as porous media, has drawn attention only recently. Apparently, Fried w43x was the first to show that foam, because of its unique structure, reduces gas flow in porous media. Foam within porous media possesses a unique structure and rheological properties. In particular, in both preand in situ-generated foams, pore sizes impose constraints on the foam texture. If the characteristic bubble size is much smaller than the characteristic pore size, the Table 1 Capillary pressures at which single films and bulk foams created from solutions of 0.001 M sodium dodecyl sulfate break w41x NaCl concentration Žmol ly1 . Film type 0.001 0.1 0.4 Ordinary thin Common black Newtonian black Critical capillary pressure ŽkPa. Single film Bulk foam G 20 G 100 G 120 5 50 100 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 135 foam confined in pores does not differ from bulk foam. In the opposite case, the foam in a porous medium is a network of thin liquid films Žlamellae. spanning the pore channels. The transition from the bulk phase to the ‘pore-confined’ regime is very complex and, as a rule, is discussed for ideal homogeneous capillaries or bead packs w44,45x. However, a more reasonable conceptual model of natural porous media is a network of interconnected capillaries of different sizes, which may contain constrictions and enlargements. Recent experiments w46x show that the motion of pregenerated foam through channels with alternating narrow and wide regions demonstrates characteristic patterns. Just after the beginning of injection in a channel comprised of alternating narrow and wide parts, as depicted in Fig. 5, the lamella moves freely along the narrow part. If the wide part is free of foam, the lamella, reaching the junction, stops, blocked by the surface energy required to flow through the enlargement of the channel. The gas pressure stretches the lamella forming a bubble. The bubble either disappears, and the bubble swelling restarts when the next lamella reaches this point; or, after touching the wall of the enlargement, the lamella jumps to the free zone and stops in a new, downstream position. Advancing bubbles jump sequentially from the narrow to the wide parts of the channels and stop. When the wide part of the channel is filled entirely with foam, the first lamella passes freely through the next narrow part, and the process repeats in the next enlargement. The manner of foam generation makes the distance between the drifting lamellae of the same order as the length of the narrow regions. In doing so, the specific film channels of approximately the size of the narrow part of the original channel arise. In a quasi-steady state, when the channel is completely filled, the lamellae drift along the narrow zones freely while their flow through the effective channels in the enlargements looks like stick-slip motion. Therewith, advancing lamellae change Fig. 5. The sequence of frames depicts the displacements of the lamellae in the enlargement of the channel w46x. Flow is from the bottom upwards. The white and black arrows indicate the positions of the pertinent lamellae at different times. The channel consisted of 10 cubic cell elements, 8 mm3 , and 10 narrow elements of 8-mm length, and 2.5 mm2 cross-section. 136 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 the topology of bulk foam in the enlargements, and their motion is jumpy, due to the tendency of lamellae to minimize their surface energy. No rupture or birth of moving lamellae occur. When lamellae pass each other, the bulk foam in the enlargements deforms, but its structure does not change. Foam flowing through the channel of an abruptly broadened cross-section, forms an effective channel, through which lamellae and gas are transported. The foam lamellae are astonishingly stable. Such observations w46᎐48x favor a crucial role for the solid skeleton in foam patterning in porous media. In the presence of solid surfaces, the interface between the two fluids takes up configurations whose equilibrium and stability depend upon the contact areas and the interfacial free energies between the various phases. The solid skeleton imposes a structure, balancing the surface and capillary forces. The wetting films covering the pore walls compete with the lamellae. The energetics of lamellae and wetting films dominate foam behavior in porous media, especially in strong foams, with unbranched lamellae causing strong blocking of flows Žsee Fig. 6. w17,49,50x. 2.3. Wetting films and menisci in pores: limiting capillary pressure As an introduction to the distribution of phases in porous media, we briefly discuss an ideal capillary system of two immiscible fluids, gas and liquid, where the liquid wets the solid matrix completely, as for surfactant solutions and other liquids with a low surface tension. In general, the solid matrix surface has a varying curvature, and Žexcluding external constraints such as gravity. if the pressure in the gas phase is constant, then, by the Laplacian capillary pressure acting on curved liquid᎐gas interfaces and the disjoining pressure acting in the wetting films, sets up a pressure gradient causing the liquid to flow. This flow tends to equalize the pressures and to establish constant curvature conditions over the menisci of the liquid᎐gas interface. The Fig. 6. Sketch of gas and liquid distribution in a strong foam in a granular medium. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 137 wetting phase can form a bulk phase bounded by menisci and wetting films coating the varying curvature pore walls. The wetting films have a varying curvature and thickness to balance the capillary and disjoining pressures. How do the menisci coexist with the films? A detailed thermodynamic theory of equilibrium between wetting films and menisci of the bulk fluid in pores of a different geometry w51᎐59x, can be summarized as follows. In heterogeneous pore systems, the narrowest pore regions are completely filled with the wetting fluid held there by capillary forces. Thin films coat the walls of unfilled pores, interconnecting menisci of capillary held fluid. Consequently, fluid in pores is interlinked, and the local capillary and disjoining pressures are dependent throughout the entire global pore network due to this interlinking. Fluid flows rapidly in response to local changes in capillary pressure. Drainage may cause thinning of wetting films and their transformations into ultra-thin films in the case of type 2 disjoining pressure isotherms ŽFig. 3.. In Fig. 3, the first stable branch of curve 2 Žthinnest films., corresponds to so-called ␣-films of. thickness of a few molecular diameters, which behave as matrix-bonded films w54x. The second stable branch corresponds to the ‘thick’ -films. A metastable -film can undergo a spontaneous transition to a more stable state of ␣-film, sometimes, via ‘rupture’ of a thick -film accompanied by the formation of a number of droplets, connected by ␣-films to the menisci in the pore corners w56,60x ŽFig. 7.. As reported by Derjaguin and others, the probability of the formation of a critical nucleus of ␣-film depends upon the capillary pressure of the menisci, and upon the film surface area w56᎐58x. In foam coarsening, the transitions from metastable -films to ␣-films are important also. Victorina et al. w61x and Churaev et al. w62x experimented with gas bubbles immersed in liquid-filled capillaries. A water-filled capillary with a gas bubble inside was placed in water-filled container with controlled external pressure. The bubble volumerlength depends on the applied pressure in the liquid phase. When the liquid pressure was varied periodically, the authors observed hysteresis: after the first stage of the pressure decrease and corresponding bubble increasing in size, the initial bubble length could not be restored upon exceeding a certain critical length. The bubble extends irreversibly as the water pressure is Fig. 7. Illustration of a transition from metastable -films to ␣-films in a wedge w60x. 138 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 decreased to a certain critical pressure drop. The menisci of such an extended bubble forms a finite contact angle despite their semi-spherical shape in the initial state. This effect has been explained by Derjaguin and co-authors as a transition of -films, coating the capillary wall inside, into ␣-films Thus, the transition of -films to ␣-films can drastically change the global structure of the liquid-filled capillary network w63x, and the conditions for generation, flow and residence of foam. On the macroscale of a porous medium, the two-phase system is described by the capillary pressure or the Leverett function w54,64᎐66x, which depends upon the saturation of one of the phases. This function links the pressures in coexisting phases and implicitly expresses the pressure difference, or the capillary pressure, via the local degree of saturation. In other words, the Leverett function serves as an equation of state for the two-phase system in the porous medium and determines foam patterning and rheology at a given degree of saturation. The ‘macroscopic’ capillary pressure reflects the main peculiarities of the pore network, with a critical capillary pressure above which no flux of the wetting fluid occurs in displacement experiments w64,65x. The associated saturation of the wetting phase is called residual saturation. In foams, however, the limiting capillary pressure also plays an important role w67x. The limiting capillary pressure is the pressure at which foam lamellae rapidly coalesce. Khatib et al. w67x experimentally studied the capillary pressure and the effluent bubble texture of foam generated before injection into sand or glass bead packs as the gas fractional flow is increased, while the gas flow rate stayed constant. Recall that the fraction of gas in the stream depends on the moving-gas saturation, estimated from the measured residence time of flowing gas w64,68x. The increase of gas fractional flow is associated with the growth of gas saturation, and correspondingly, with the increase in capillary suction. With increasing gas fractional flow at a constant gas flow rate, the capillary pressure approaches the limiting capillary pressure. At this pressure, foam coalescence causes the capillary pressure to drop abruptly and the foam texture to coarsen. So far, no model completely explains this phenomenon Žreviewed in w17,50,69x.. Foam demonstrates a characteristic of self-organized criticality w70᎐72x: coalescence and displacement of a coarse-textured foam usually maintains the capillary pressure near the limiting value w67,73x. 2.4. Mechanisms of bubble generation in porous media: gas path closure Foam may be forced Žinjected. into a porous medium, or it may be generated in situ. Garrett w34x has recently reviewed the injection of pregenerated foam into a porous medium. In-situ generated foams in the critical regime of limiting capillary pressure exhibit specific features. The three main mechanisms of in-situ foam generation are w48x: snap-off, lamella division and leave-behind ŽFig. 8.. Although all of these foam creation processes occur simultaneously, most authors believe that snap-off dominates blockage capability Že.g. w17,50x.. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 139 Fig. 8. Mechanisms of bubble generation in porous media. Ža. ‘Leave-behind’: arrow indicates a lamella formed by gas fingering; Žb. snap-off. when gas emerges from a wet constriction into a water-filled pore, the interface forces allow a leading portion of the gas to separate into a bubble; Žc. lamellar division: arrows indicate the directions of movement w73x. Roof w74x proposed that bubbles or drops are generated on certain sites. Consider gas injection into a pore filled with a wetting liquid, with the pore modeled as the inside wall of a torus Ž‘a bagel hole’., as analyzed by Mohanty w53x. This model reflects the main features of pore constrictions and allows some general conclusions about snap-off. We assume that the flow is slow enough to ensure that the liquid interfaces approximate their equilibrium shapes w53x. Three steps are key in bubble snap-off. First, the bubble moves into the constriction, depositing a film of wetting liquid ŽFig. 9.. We consider the wetting film as a layer of uniform thickness in the AT region. Then, the total curvature of the gas finger, 2 H, is approximately that of the pore wall in region AT, i.e. 2 HŽ z. s y 1 rm q 1 RH . Ž3. In Figs. 9 and 10, rt is the pore radius at the narrowest position Žpore throat., point T; z is the axial distance from midpoint of the throat, R H and rm are the principal radii of curvature of the pore wall Ž rm is the meridional radius of the torus.. In region BC, the meniscus looks like a spherical cap of radius r h . its K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 140 Fig. 9. Dynamics of collar formation in a toroidal pore. G, gas; W, water w53x. curvature becomes 2 Hc s 2 rh . Ž4. The curvature must continuously vary from point A to point B. Therefore, the curvature of the gas᎐fluid interface has a minimum at point A. Whether the curvature in region BC is greater or less than that at point T depends on the extent of the gas protrusion, i.e. the value of radius r b . The gas pressure is uniform, so, the pressure in the wetting phase Pw can be estimated using the Laplace equation of capillarity, Eq. Ž1., written in the form: Fig. 10. Radii of curvature of the surface of revolution described by the curve LM in Fig. 9. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Pw s Pg y 2 H , 141 Ž5. Because the pressure Pw achieves a maximum at point A ŽFig. 11., it drives the wetting fluid away from point A. Hence, the film pinches down at point A. At the same time, the wetting fluid, which resides between points A and B, will migrate toward the throat, point T, its interface taking the shape of a constant curvature surface. The accumulation of the wetting fluid at the throat is usually referred to as a wetting collar w48,53x. The flow rate limits the volume of the growing collar through hydraulic constraints. Therefore, we select as possible shapes the collar and the lens, which are stable to volume conserving perturbations. Fig. 12 demonstrates the characteristic features of the transition ‘stable collar to stable lens’ w75x. If the volume of the wetting film collected at the throat is small, it can be accommodated in a wetting collar around a stable neck meniscus, points A-B-C. But if the volume is larger than the maximum stable collar volume, point C, then the collar snaps off, creating a lens of wetting fluid at the throat. The schematic illustrates that the hydrodynamics of gas fingering is a secondary factor in lens formation because capillary forces dominate. Hence, the problem of the bridging transition has, mainly, a thermodynamic, rather than as hydrodynamic, origin, resembling capillary condensation in pores w54,57,58x. Recall that if a solid surface is placed in contact with a vapor of wetting liquid, an adsorption film is formed. The film thickness depends on the vapor pressure and surface curvature. In narrow pores, as the vapor pressure increases, the films on the opposite walls of the pore coalesce and form a lens as the film thickness exceeds a specific value for a given pore geometry. Capillary condensation has been studied for different pore geometries. In particular, Everett and Haynes w76x comprehensively explained the example of a cylindrical capillary, spherical particles have been studied in detail by Neimark w77x and Neimark and Rabinovich w78,79x. Disjoining pressure effects have been included Že.g. w54,56x.. Turning to the problem of foam distribution, the most important conclusion is Fig. 11. The pressure Žarbitrary units. in the wetting film at the end of the second stage w53x. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 142 Fig. 12. The typical capillary pressure Pc Žarbitrary units. as a function of the liquid saturation S Žfraction of liquid phase within the pore, arbitrary normalization. for a pore shape which is a surface of revolution w75x. that the pore size, r, where foam lamellae may occur, can be estimated using the condition of thermodynamic stability for a film coating the pore walls. For a cylindrical capillary, following Derjaguin’s supposition of the additive action of the capillary and surface forces in curved wetting films, the capillary pressure in the bulk liquid, in equilibrium with the film in the cylindrical capillary, equals: Pc Ž h . s ryh q ⌸ w Ž h. , Ž6. here ⌸ w is the disjoining pressure of the wetting film and r is the radius of the capillary. The film stability condition is defined by the inequality ѨPcrѨh F 0 w56x, namely: Ž r y h.2 q ⌸Xw Ž h . F 0. So the upper limit of thickness for a thermodynamically stable film, h c , is defined K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 143 by, ⌸Xw Ž h c . s y Ž r y hc .2 , and the corresponding capillary pressure by Eq. Ž6. at h s h c . Eq. Ž6. assumes that the disjoining pressure in the curved film can be approximated by the disjoining pressure isotherm for flat films, which is credible when h < r, and the radius of the capillary is greater than the characteristic length of action of molecular forces. The estimate is a required condition for lamellae to be generated by capillary instability. A similar thermodynamic approach describes the process of bubble generation in mesoporous catalysts due to liquid decomposition w80x. However, in the nanometer scale pores fluid᎐solid interactions should be taken into account explicitly w81᎐85x. Modifying the conditions for capillary condensation, density functional theory gives a better understanding of capillary condensation phenomena in nanoconfinements Že.g. w86,87x.. Thus, whether or not the lens generation takes place at a potentially generating site depends on the capillary pressure, which is always limited from below. At least two pressures are related to snap-off phenomenon Ždynamic capillary condensation.: the pressure required for snap-off, Psn , and the pressure Psh required to displace wetting liquid from the pore throat. As a rule. Psh ) Psn . Thus, for snap-off, the capillary pressure must first rise and then fall ŽFig. 11.. The ratio PsnrPsh f O Ž1. is a function of the pore geometry and the substrate wettability w17,53,74x, e.g. PsnrPsh s 1r2 for a wide cylindrical capillary w74x. 2.5. Transformation of lens into lamella When a surfactant is present, the configuration at point D ŽFig. 12. with two interfaces touching one other, may be stable. A stable flat lamella may form as the result of the transition of the ‘thick’ film Žwith bulk solution properties. to a new thermodynamic state. The characteristic feature of such a transition is the appearance of a contact angle between the film and the bulk liquid w4,88,89x. Consider a simplified physical model of lamella formation from a lens spanning the cylindrical capillary of radius R ŽFig. 13. w93x. This is the limiting case of the toroidal pore when rm ª ϱ and R s rt s R H . Under the capillary suction pressure, which here equals Pg y Pw s liquid gH, the lens squeezes and forms a lamella of radius R 0 Žsee notes in Fig. 13.. In sufficiently wide pores Ž10᎐100 m wider., the characteristic size of the meniscus is greater than the characteristic size of transition zones between the thin liquid film and the bulk liquid w90᎐92x. Therefore, in the description of the transition from a lens to a lamella, classic thermodynamics, in which disjoining pressure effects are neglected, seems suitable. Then contact conditions may be modeled via the contact angles of the meniscus with the substrate, R , and with the lamella, 2 0. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 144 Fig. 13. Idealized picture of lamella formation within a tube w93x. The equilibrium shape of the meniscus can be obtained by integrating the Laplace equation 2 H s 1 d r dr Ž r sin . s Pg y Pw s Pc , tan s dz dr . Ž7. Eq. Ž7. has the integral r sin s Pc 2 r 2 q B, B s const Ž8. Making use of the boundary conditions d z r d r s tan 0 , r s R 0 , and d zrd r s tan R , r s R we find the relation between the capillary suction pressure and the lamella radius R 0 as: RPc 2 s cos R y ␣ sin 0 1 y ␣2 , Ž9. where ␣ s R 0rR. This equation demonstrates hysteresis in lamella formation. Indeed, as soon as the capillary suction pressure reaches a value Pg y Pw s P0 s 2 cos RrR, the lens immediately generates a lamella of radius r 0U . If the pressure rises, the lamella radius increases monotonically. If the capillary suction pressure decreases, the lamella shrinks to a critical radius r 0min , then suddenly disappears. The critical radius r 0min minimizes Eq. Ž9.. Zorin et al. w93x performed a detailed experimental analysis of lamella formation by interferometry in transmitted light w94x. Water solutions of sodium dodecyl sulfate ŽSDS. and sodium chloride ŽNaCl. were used as liquids. The authors used capillaries of radius R s 120 m drilled within a finegrained filter, R f 0. Fig. 14 demonstrates the characteristic features of the transition. Line 1 corresponds to a solution of 0.32 M ŽNaCl. q 0.05% ŽSDS., s 31.8 mN my1 , and line II a solution K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 145 Fig. 14. The capillary pressure as a function of the lamella radius. Dashed lines represent the theoretical curves obtained by ignoring the action of the line tension w93x. of 0.45 M ŽNaCl. q 0.05% ŽSDS., s 31.25 mN my1 . The contact angle 0 and the line tension of lamella were also measured. The latter is defined as ␥ s 2 cos0 y R0 , where ␥ is the lamella tension. For the first solution, the line tension was k I s 3 = 10y9 N, for the second k II s y1.5 = 10y9 N. For both solutions within the whole range of the measured lamella radii, the contact angle 2 0 was close to r10. In this model, the lamella contact angle 2 0 and the lamella tension are assumed to be constant. Thermodynamics of lamella tension and contact angle, which are indeed functions of lamella thickness and disjoining pressure, was elaborated by Kralchevsky et al. w95x. Because of the overlapping of the interfacial layers, in thin films like foam lamellae, the film tension becomes a function of the disjoining pressure and film thickness w4,57,58x. When the film thickness exceeds a few nanometers, the film tension ␥ usually differs only slightly from twice the value of the surface tension, , of the interface between the bulk solution from which the film is made and the surrounding phase w90x. Nevertheless, even this small deviation causes a finite contact angle between the meniscus and the lamella. We 146 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 15. Ža. Scheme of lamellae distribution in a bamboo-like channel. Dashed lines are the initial positions of lamellae, and bold lines represent the lamellae under a load. Žb. Parameters needed for calculation of the bubble volume variation and capillary force. present the relation between the lamella tension and its disjoining pressure without proof Že.g. w4x.: ␥ s 2 f q ⌸ Ž h iiq1 . h iiq1 s 2 q ϱ Hh ⌸ Ž h . d h q ⌸ Ž h iiq1 . h iiq1 . Ž 10 . iiq1 Note, that his relation is valid for flat lamellae and may serve as a reasonable approximation for curved lamellae. 2.6. Lamellae in con¨ ergent (di¨ ergent) pores The thermodynamic stability of curved lamellae, spanning convergent Ždivergent. pores, differ from those of flat lamellae. For simplicity, consider a pore channel as a cylindrical capillary with a radius that varies periodically with the distance ŽFig. 15.. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 147 To study how the lamella thickness depends upon the pressure drop, acting across a curved lamella, one must generalize the Laplace formula for films with finite thickness. Using statistical average, Buff w96x showed how the Gibbs thermodynamic surface tension w97x can be determined Žreviewed in w98᎐102x.. A sophisticated thermodynamic analysis of curved thin films has been done by Ivanov and Kralchevsky w103x. We, however, use a somewhat simplified model of the thin film, in which the disjoining pressure acts as an additional load distributed over each lamella surface w35,36x.. The model assumes that the contributions of capillary and surface forces are additive. Thus, assuming that the inequality Pi - Piq 1 , holds, the balance of forces at the interface from the side of i-th bubble gives the following expression w42,57,58x Pi q 2 R iiq1 s ⌸ Ž h iiq1 . q Pw , Ž 11. where h iiql is the film thickness, ⌸ Ž h. is the disjoining pressure isotherm of a flat lamella, R ii q 1 is the mean radius of curvature of the lamella, provided that the film thickness is much less than the radius of curvature and is the surface tension of bulk liquid. Note that the use of the film surface tension, f , defined for flat films through Eq. Ž10., in Eq. Ž11. instead of the bulk value, , leads to inconsistent corrections. For the interface between the liquid film and bubble i q 1, the symmetric equation holds, Piq1 y 2 R iiq1 s ⌸ Ž h iiq1 . q Pw . Ž 12 . Eqs. Ž11. and Ž12. imply that the ordinary Laplace equation for a thick film between two phases gives the mean radius of film curvature: Ry1 iiq1 s Ž Piq1 y Pi . 4 Ž 13 . and the isotherm of the disjoining pressure gives the film thickness as: ⌸ Ž h i iq1 . s 12 Ž ⌬ Piq1 q ⌬ Pi . , Ž 14. where ⌬ Pi s Pi y Pw , is the capillary overpressure in the i-th bubble. First order corrections for a finite film thickness lead to the substitution of the bulk surface tension by the film tension ␥ as defined by Eq. Ž10.. However, these corrections seem to be insignificant for practical estimates. Recall that the capillary overpressures in Eq. Ž14. can be expressed via the disjoining pressure isotherms prescribed for ‘substrate-wetting film-gas’ conditions, i.e. via the background capillary pressure Pc s Pi y Pw . Thus, the equilibrium properties of curved lamellae also depend upon environmental conditions. Eq. Ž13. has many solutions, i.e. lamellae might be placed at different positions K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 148 along the pore. However, the mechanical equilibrium conditions at the contact line impose a selection criterion for lamella positions within an axially symmetric pore. Discussing mechanical equilibrium, we can ignore the Plateau borders and assume that the lamella intersects the pore walls at a 90Њ angle, as observed experimentally w48x. Eq. Ž13. reflects minimization of film surface areas Žsee Fig. 15.: R iiq1 s rrsin␣ Ž r . N xsa i , tan␣ s dr dx . Ž 15. where x is the longitudinal coordinate, r s r Ž x . is the pore profile, and a is the lamella chord coordinate. Laplace Eq. Ž13. defines the equilibrium locations of lamella for a given pressure differential; variations of the pressure differential lead to varying equilibrium lamella positions. However, not all the equilibrium configurations are stable. For stability, a small shift in the lamella position must give rise to a net restoring force, implying w104x d dx ž Piq1 y Pi y 4 R iiq1 / - 0. Ž 16 . If the pressure differential does not depend on the lamellar position, the stability condition depends only on the shape of the channel through Eqs. Ž15. and Ž16.. Then Eq. Ž16. shows that stable equilibrium lamellae can exist only over convex, d␣rd x - 0, segments of the pore channel. If, however, the pore has both convergent and divergent parts, then the equilibrium positions correspond to all local extremum points of the channel radius. Minimum corresponds to stable and maximum to unstable equilibrium positions; because the larger the lamella perimeter is, the larger its surface energy. 2.7. Chain of lamellae: correlation length Capillary forces tend to fix the lamellae at the pore throats, while elastic forces, caused by gas compressibility, shift the lamellae into new equilibrium positions, resulting in the equilibrium states observed in experiments. The situation resembles ‘commensurate᎐incommensurate’ phase transitions in solid state physics w105,106x and related problems, in which binding and pinning forces compete w107᎐112x.. Specific features can be revealed by the simple one-dimensional bubble train model w113᎐115x. Imagine a wavy channel as a rigid capillary with radius, r s r 0 q ␦cos 2 x ž / , Ž 17 . where the x-axis coincides with the capillary axis of symmetry. r 0 , , and ␦ are characteristic scales of the porous medium. Assume, for simplicity, that the pore K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 149 aspect ratio is low, Ž ␦rr 0 < 1. This assumption, however, can be altered without changing most of the physical conclusions. Consider the train of N lamellae, whose chord centers lie at x s a1 ,a2,...,a N . and assume that in the initial undeformed state, the foam is perfectly ordered, with the lamellae connected by links Žbubbles, or train carriages. of ‘carriage length’ K with integer K. Under a load, the i-th bubble is deformed, to a length a i y a iy 1 with an accuracy of O Ž ␦rr0 ., w113᎐115x. The resulting elastic force acting upon the i-th lamella can be written with the same accuracy as: f e s r 02 Ž Piq1 y Pi . Ž 18 . Making use of the ideal gas law, yields Piq1Ž a iq1 y a i . s Pi Ž a i y a iy1 . s . . . s Pg K , Ž 19. where Pg is the initial gas pressure in an individual bubble. We introduce a new function, i the displacement of the i-th lamella from its initial position at the throat. Then a i s iK q i q r2 and Eq. Ž19. takes the form: iq1 y i s Pg K Piq1 y K, Ž 20 . Note that the capillary force has a form similar to Eq. Ž18., but where the Laplacian pressure drops, Eq. Ž13. replaces the pressure differential. Then the balance of capillary and elastic forces gives, with accuracy O Ž ␦rr 0 ., the equation Piq1 y Pi s y 8␦ r0 sin 2 i ž / . Ž 21 . The boundary conditions for the 0th and the last, N y 1. lamellae are, P0 s Pg Pext , PN s Pg P , Ž 22 . where Pext and P are the dimensionless external pressures applied to the train. Thus all lamellae in the loaded train obey Eqs. Ž20. ᎐ Ž22.. The corresponding Cauchy problem for the Ulam map wEqs. Ž20. and Ž21.x in a different context w109x, applies to foam patterning and to foam elasticity. To estimate the limits of validity of the continuum mechanics, consider the linearized version of Eqs. Ž20. ᎐ Ž22., i.e. the asymptotic case of i ª 0, P, Pext ª 1. For the semi-infinite train N ª ϱ the solution can be written in the form w116x, i s aexp Ž ywi . , K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 150 where a is a constant defined by the boundary conditions, and W satisfies the following transcendental equation: cosh w s 1 q 2 2 K Ž 23 . Here, the parameter s 4␦ Pg r 0 , Ž 24 . measures the strength of capillary forces with respect to elastic forces. In regimes where cooperative properties of foam play a crucial role, K < 1, and the screening length Žcorrelation length. Ncor f 1rw of a train is estimated as: Ncor f 1 2 K ' K < 1. , Ž 25. In the opposite case, when either gas compressibility or carriage length are large, pinning forces dominate binding ones, so: Ncor f 1 ln Ž 2 q 4 K . 2 , K 4 1 42 . Ž 26 . Explicit barrier K s 1r4 2 is the stability criterion for two bubbles w115x. Thus. in the former case, the load is distributed over a large number of lamellae. In the latter case, the first few lamellae bear most of the load. Therefore, the continuum model wEqs. Ž20. ᎐ Ž22.x applies to foam patterning when K < 1; and discrete effects should be taken into account for K 4 1r4 2 . The above mentioned inequality assigns a specific meaning to the terms of compressible and incompressible and determines the range of parameters over which a theory of foam elasticity can be constructed. Restriction Ž25. allows simplified physical treatment w117x. Consider two characteristic values of pressure perturbation associated with foam in porous media on the microlevel of pores. The first characteristic pressure drop is the Laplacian capillary barrier ␦ Pc f 4rr 0 imposed by the inherent structure of the pore matrix. This pressure drop is required to push a single lamella through a pore constriction. To estimate the second pressure variation caused by the change of bubble volume, we imagine displacing a single lamella, with other lamellae remaining at their initial positions. For an ideal gas, the pressure perturbation within the deformed bubble is ␦ Pg f Pg ␦VrV f Pg r K, where V is the initial cell volume. This pressure variation does not depend on the external pressure drop, and is an inherent characteristic of the foam. We expect that whenever Pgr␦ Pc ) 1 Ž4K r Pg r 0 - 1., the foam will move as a whole, because the external pressure drop will be redistributed over all foam cells. In the opposite case, lamellae withstand variations of pressure of order Pg , so the ability of each lamella to sustain a local pressure variation of order Pg determines foam patterning. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 151 2.8. Superstructures: analysis of the phase portrait Phase portraits clarify the possible equilibrium states of trains with different interparticle interactions w106,111᎐115x. In the continuum limit, going over from a discrete variable to a continuous variable s s 2 Ki , and introducing dimensionless variables and functions 2 i r ª , Pi rPg ª , we rewrite Eqs. Ž20. and Ž21. as: y dp d ds ds s s 1 p K sin. Ž 27 . y 1, Ž 28 . Eqs. Ž27. and Ž28. has first integral, ln p y p s E q K cos, Ž 29. where E is a constant. The translational symmetry ª q 2 n, n s 0, "ϱ, allows us to consider Eqs. Ž27. and Ž28. over the whole range, g Žyϱ,ϱ.. The phase portrait ŽFig. 16. shows two kinds of singular points: hyperbolic and elliptic Žsee also w106x.. Points p s 1, s 2 n, n s 0 " ϱ are hyperbolic, representing configurations of the bubble train where lamellae are attached to the narrowest parts of the channel; so that the surface energy of lamellae is minimal. Points p s 1, s Ž2 n q 1., n s 0, "ϱ are elliptic, with cycles describing states in which lamellae avoid wide parts of the channel because the surface energy is maximal. The hyperbolic points are connected by a separatrix, which has two branches within each period of the phase portrait. One exits the left point. passes above the line Fig. 16. Phase portrait of the system wEqs. Ž27. and Ž28.x. The separatrix is denoted as ‘S’. Arrows show the direction of the increasing arclength s. 152 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 p s 1 and enters the right hyperbolic point. The other branch goes in the opposite direction under the line p s 1. The separatrix is described by, ln p y p s K Ž cos y 1 . , Ž 30 . which follows from Eq. Ž29.. At the separatrix, E s EUU , EUU s y1 y K . Ž 31 . The separatrix subdivides the phase plane Ž p, . into three parts. Parts I and II represent infinite slippage Ž E - EUU ., while part III corresponds to bounded motion Ž E ) EUU .. All integral curves in domain III pass through points at which pN s 1, so they describe bubble trains with free boundary lamellae. The integral curves of domain I correspond to a bubble train with a blocking lamella at N y 1 s 0, pN ) 1. Domain II represents extended bubble trains with fixed boundary lamellae N y 1 s 0, pN - 1 Each of the singular points p s 1, s 2 n, creates four branches of the separatrix corresponding to the configurations of the chain. Therefore, if some pattern incorporates a single or a few singular points, the respective texture of the train will be complex. If we associate the current number of lamellae with time, then the finite motion in domain III can be treated as a motion with a finite period T, i.e. the number of bubbles accumulated by a single cycle is finite. Near the separatrix, T tends to infinity, i.e., the solution describes an unconfined bubble train. The solution also demonstrates the effect of the ‘irreversibility’ of bubble train displacement. After the action of a critical pressure drop some of the lamellae never recover their former positions ŽFig. 17.. The solution looks like a solitary ‘domain wall’ w106᎐111,113x. s 4tany1 exp y wŽ x y x 0 .x q O ž( / K , xs ( K s. Ž 32 . For a given pressure drop, all the lamellae behind the ‘wall’ at s ª yϱ shift by the period s 2 and stay there after unloading. Lamellae at s ª ϱ, ahead of the wall, keep their undeformed state s 0. The domain wall matches the second zone Žthe undisturbed zone. where the lamellae are pinned at the equilibrium state with the first zone, where lamellae are displaced by one period ŽFig. 17.. The wall thickness can be estimated as O Ž 6Kr6 . w113x. The continuous model wEqs. Ž27. and Ž28.x, is valid whenever collective effects dominate and lamellar displacements and pressure vary weakly over the train. In the opposite case, the discrete map wEqs. Ž20. and Ž21.x resembling the Ularn model w112,118x. has stochastic solutions w114,115x. Discreteness results in more complex foam patterning; in particular, glass-like ordering of lamellae. The ability of K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 153 Fig. 17. Ža. Sketch of a domain wall and the pressure distribution Pi s PirPg over the chain; 1, 2 ᎏ chain under load; 3 ᎏ after unloading. Žb. The displacements of lamellae i s ir; 1, 2 ᎏ chain under load; 3 ᎏ after unloading. N s 80, s 10y3 , K s 1. variations of the initial data Ny 1 to abruptly alter the chosen branch of the separatrix and the sensitivity to the initial data of the distribution of the points along the trajectory manifest through chaotic ordering. While the domains may occupy varying regions, the associated energy remains almost the same w114,115x. Sensitivity to initial conditions results from the finite correlation length. Fig. 18 shows the phase portrait associated with glass-like patterning. 2.9. Start-up pressure drop The previous analysis shows that, in porous media, foam drastically changes the rheological behavior of the gas phase. In a foam, gas flows as if it were an homogeneous fluid, with a start-up yield pressure drop, required to initiate flow through porous media w14,17,45,50,119x. 154 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 18. The typical phase portrait for Eqs. Ž20. ᎐ Ž22., s 0.016, K s 10, N s 2000, and six starting points. The dashed line represents the separatrix of the continuous model w114,115x. The explanation of the yield pressure drop usually assumes that foam is incompressible w120x, leading to an enormous capillary barrier. Considering capillarity solely, we estimate the pressure required to displace a lamella from the pore as approximately ⌬ Pind f 4rr 0 . We can also estimate a similar barrier for a strong foam, assuming that before deformation foam lamellae are attached to each throat so that the distance between them is of the order of the pore radius, r f r 0 . If we neglect the gas compressibility, we need to overcome a huge capillary barrier ⌬ Ptotal f 4rr = number of lamellae, to change the pressure by an amount 4rr 0 for each lamella. In a strong foam of size L, the number of lamellae A Lrr so ⌬ P f 4 L r r 2. Experimental values are significantly lower than this estimate w121,122x. Thus, the capillarity alone cannot explain the observed critical pressure drop. However, as shown above, more important is the ratio of pinning Žcapillary. forces to binding Želastic . forces, expressed as the parameter wEq. Ž24.x. In the continuum model, the pressure variation is finite for all reasonable solutions Ždomain III in Fig. 16.. The upper boundary for the maximum dimensionless pressure drop, G, is the maximum of the separatrix, i.e. the solution of algebraic equation, ln Ž 1 q G . y G q 2 K s 0. Ž 33 . The bubble train is unable to sustain a pressure drop above G, if one of the train ends is free. The maximum external pressure rises as the bubble train acquires a K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 155 glass-like structure w114,115x. Determining the train ground state is difficult, because it depends on the pre-history of foam generation and motion. Another significant experimental parameter is the ratio between the correlation length Ncor and the size of the sample, L, which determines finite-size effects. Even for small K, when the continuum model holds, Eqs. Ž22., Ž27. and Ž28. show w115x that for Ncor 4 L, the critical pressure drop depends linearly on the number of lamellae in the train. In the dimensional form, the critical pressure drop G is: Gs 2 ␦ Kr 0 N Ž 34 . The linearity of GŽ N . shows the independent contribution of each lamella to the total pressure drop. Similar behavior might be expected for the limiting case of absolutely compressible bubbles, for which the start-up yield pressure drop is completely determined by the overall Laplacian barrier. For an absolutely incompressible gas, the critical pressure drop is also linear. For long samples, when Ncor < L, collective effects dominate and cause the screening effect: ahead of the domain wall lamellae stay undisturbed Žsee Fig. 17.. The critical pressure drop corresponds to that required to form a domain wall in an infinitely long train, as in Eq. Ž32. w113,115x: ( Gs4 Pg ␦ r0 K . Ž 35 . The two limiting regimes wEqs. Ž34. and Ž35.x demonstrate the importance of the size ratio. For short foams, the critical pressure drop depends on the sample length, regardless of foam texture and physical parameters. In Fig. 19, the experimental data reported by Falls et al. w121x, their Table 2, are fitted to Eq. Ž34. Ž G A r by 3, r b is the effective radius of a bubble., where the number of lamellae per unit length is inversely proportional to the volume of gas per bubble w44x. In the same picture, the expected result for a long train is presented as a dashed line r2 . Ž G A ry3 . Since the number of lamellae in the train was approximately eight, b Eq. Ž34. adequately describes the experiment. The correlation length covers the entire bubble train. Since glass-like ordering improves foam screening, determining the regimes of loading which result in a coarse-grained foam superstructure is desirable. This superstructure might be imagined as a random array of blocks with an internal crystalline-like order, separated by domain walls. In other words, the domain walls serve as apparent lamellae, and the blocks play the role of gas bubbles. Another approach to determining the critical pressure drop takes the critical pressure drop as that required to keep lamellae moving by overcoming the capillary and viscous forces, which resist their advance w117,121,123᎐127x. Rossen w124᎐127x attributes a crucial role in the appearance of the critical pressure gradient to the sharp edges, or cusps, within a pore channel modeled as a tube with a periodically varying radius. For a piecewise linear distribution of pore radii, the channel is a 156 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 19. The critical pressure drop as a function of the effective bubble radius. series of frustums of cones. If the bubbles were incompressible, some lamella positions would be prohibited by geometric restrictions, and the lamella motion would become a sequence of alternating equilibrium positions. For certain shapes of the channel, e.g. with edges, or cusps, the bubble has to jump a certain distance in order to conserve volume during the motion. Jumps occur when the bubble volume becomes a non-monotonic function of the lamella position. To estimate the critical pressure gradient, we calculate the volume-averaged value of Laplacian pressure drop per bubble w121,123᎐127x Žprovided that bubble volumes are random and uncorrelated., which represents the net work required to push a lamella through a pore. The driving pressure is the time-averaged pressure drop per bubble. ŽFor steady motion, the relation between the time and volume occupied by the bubble is linear.. The critical pressure gradient in smooth inhomogeneous capillaries, e.g. sinusoidal, is zero, while in most experiments on homogeneous bead packs w121᎐123x., in which, active channels should be smooth, it is not. Thus, both the experiment and theory cast doubt on the existing ‘dynamic approach’ to calculating the critical pressure drop. The critical pressure drop is hysteretic: to re-initiate flow, the ‘yield’ pressure drop has to increase above the initial level by 10᎐20% above that required to keep lamellae moving w128x. At least two different explanations are possible. Treating the foam as an elastic body, the critical pressure drop occurs when a dislocation deepens w110x. Equations of foam hydrodynamics, a solution of which determines an internal structure of the K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 157 foam, lead to a different critical pressure drop by considering the limit of infinitely slow flows. 3. Foam transport in smooth capillaries Even above the critical pressure drop, when the foam flows, the gas flux is several orders of magnitude less than for pure gas w13,14,17,50,123,129x. Bretherton w130x explained the non-linear friction of a bubble as the result of competition between capillary and viscous forces. Later, Hirasaki and Lawson w44x and followers w121,131,50x extended the Bretherton approach to bubble train motion through smooth capillaries and bead packs. They implicitly assumed that foam motion resembles that of a single-phase system, with the joint motion of the gas, the liquid slugs separating the isolated bubbles, and the wetting film coating the pore walls, treated as pseudohomogeneous Darcy flow Žrelating the applied pressure gradient to the resulting flow rate.. Within this approach, foam does not move relative to lamellae: all bubbles move in unison. However, the average velocities of different phases can be distinguished. In particular, the steady velocity of the gas bubbles is determined by the average velocity of the faster flowing fluid near the center of lamellae Žor liquid lenses.. Then, the volume swept out by a long bubble when moving with speed U equals the average speed Uliquid of the liquid, multiplied by the cross-sectional area of the tube. This balance estimate reveals that the average speed of the wetting fluid is less than the bubble velocity, U, by an amount WU, precisely, in the ratio Ž1 y W . s Ž1 y hϱ rR . 2 , where hϱ is the thickness of the uniform film between the front and rear menisci and R is the pore radius ŽFig. 20.. Therefore, the smaller the film thickness, the smaller the difference in the average speeds of different phases. The average speed of the gas in bubbles slightly exceeds the average speed of the wetting fluid. Because of the large viscosity of the wetting fluid with respect to the gas viscosity, the resulting increase in apparent gas viscosity is expected to be appreciable. Another scheme for foam flow in a capillary model is the wave-like motion of lamellae w113,132x. For each lamella in the train, the physical cause for lamella motion resembles that for a sailboat. Thus, the model is called the ‘sailboat model’. The lamella, exerting a tensile force, draws a meniscus Žthe Plateau border. so that a lubrication flow occurs within the meniscus and the wetting film. As for Bretherton, the main flow patterns concentrate within the wetting film. Since the lamellar thickness, radius of curvature, and tension, all self-consistently depend on the changes of the wetting film, collective motion of lamellae cannot be predicted in advance. Stationary motion requires restrictions on the pressure distribution within the train. Thus, the two models for foam friction appear to differ: the first deals with bubbles, while the second treats each lamella independently. The first considers the creep of the bubble train as a whole. The second treats motion as a wave-like displacement of lamellae in caravans. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 158 Fig. 20. Schematic of an air bubble traveling in the positive x-direction in a circular capillary. 3.1. The Bretherton mechanism of gas mobility reduction: motion of indi¨ idual bubbles Bretherton’s w130x approach, which aimed at a simple description of bubble motion, obscures the physics of foam motion. It applies to many similar problems, e.g. the motion of drops w55,129,133᎐137x and some problems of dynamic wetting Žthe so-called Landatu᎐Levich problem w138᎐142x. In the literature on soap films, this problem is called the Frankel problem w2x. It diminishes the unique dynamic coexistence of the wetting film and foam lamellae. For bubble motion in unison with the surrounding fluid, Bretherton invoked the lubrication approximation. Assuming that the gas viscosity is zero and the pressure within a bubble is constant, the volumetric flow rate is: qsy h 3 ѨPw 3 Ѩ x . Ž 36 . Here, is the viscosity of wetting fluid, h is the wetting film thickness, and x is the spatial coordinate measured along the pore. At each point x, the pressure Pw across the film is assumed to be constant which depends on the film thickness as: Pw s Pg y 2 H. Ž 37 . Here, the effects of the disjoining pressure, which are important in sufficiently thin films, are ignored. Thus, the flow within the film is controlled by a ‘dynamic’ capillary pressure and tends to establish a mechanical equilibrium between capillary and viscous forces, e.g. if the driving force is switched off, the bubble assumes a shape of constant mean curvature, as predicted by Eqs. Ž36. and Ž37.. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 159 The mass balance equation, written in the form: Ѩh q Ѩt Ѩq Ѩx s 0, Ž 38 . completes the model. Eqs. Ž36. ᎐ Ž38. plus boundary conditions specify the physical situation. Bretherton considered steady motion of the bubbles and reformulated Eqs. Ž36. ᎐ Ž38. in a bubble fixed reference frame, s s x y Ut, where U denotes the speed of the bubble. He subdivided the steady profile of the wetting film at each bubble end into three regions ŽFig. 20. within which different mechanisms determined the shape of the interface. Where flow is slow, capillarity dominates the shape of the meniscus Žregion 3., so this region is designated as part of the equilibrium meniscus. Behind lies a uniform film Žregion 1., reached asymptotically. Between is the transition region 2, within which the film profile is governed by interplay of viscous and capillary forces. Assume that within region 2, the slope of the gas᎐liquid interface is small and the film thickness is much less than the pore radius. The interface shape in a bubble-fixed reference frame is d3 h ds 3 s 3U h y hϱ h3 . Ž 39. Here hϱ is the ‘asymptotic’ film thickness in the middle part of the bubble. Transforming, h ª hrhϱ and s ª sŽ3U r .1r3rhϱ , scales Eq. Ž39. so that the capillary number Ca s 3Ur Žwhich measures the ratio of viscous forces to capillary forces. and the ‘asymptotic’ film thickness do not appear explicitly: d3 h ds 3 s hy1 h3 . Ž 40 . Bretherton numerically matched the thin film solution to the thick film solution at the ends of the bubble and assumed that the menisci have the curvatures equal to their static value Ža modern numerical solution of autonomous third-order ordinary differential equations, like Eq. Ž40., is reviewed by Tuck and Schwartz w143x., i.e. he assumed, that equilibrium menisci exist and look like hemispherical caps. He did not assume that a film with uniform thickness hϱ , coexists with the menisci; the ‘asymptotic’ film thickness, hϱ tends to zero with the capillary number, i.e. the equilibrium bubble consists of two hemispherical menisci and a segment of the dry tube, regardless of the wettability of the tube’s surface! If the fluid wets the substrate completely, the Bretherton solution is not the asymptotic limit Ca ª 0, but an intermediate asymptote w144,145x. Nevertheless, the velocity of an individual bubble does slightly exceed the average speed of the fluid, Uliquid , by an amount, Ws U y Uliquid U sA 3U ž / 2r3 , A f 2.29, Ca s 3U ª 0. Ž 41. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 160 Bretherton also pointed out that the dynamic pressure drop across such a bubble, besides the Laplacian pressure drop: ⌬Ps B 3U ž / 2r3 R , B f 3.58 Ž 42 . Experimental results by Fairbrother and Stubbs w146x and Taylor w148x on an air᎐water system, flowing within circular tubes, indicated that, in the range of 10y4 - Ca - 10y2 , the wetting film thickness was proportional to Ca1r2 , in contradiction to the Bretherton theory. ŽThe Marchessault ᎐Mason data also differ from the Bretherton predictions w147x.. Bretherton also experimented using aniline and benzene as the liquids. Theory and experiment disagreed for capillary numbers smaller than 0y4 . Schwartz et al. w149x in 1986 Žsee also w150x. repeated Bretherton’s experiments, using water as the liquid phase ŽFig. 21.. The main parameters and Ca, were measured independently with high accuracy was expressed via the decrease in the length of the liquid slug AS and the distance of travel D as W s ⌬ SrD. The capillary number was found by measuring the A-to-B travel time of the slug. The experimental findings are presented in Fig. 21d.. As in the Bretherton experiments, the discrepancy between theory and measurements increases as the speed is reduced, where the asymptotic theory should work best. Explanations for the discrepancy between the theory and the experiment proposed in the literature include wall roughness w150x, instability of the meniscus w130x, intermolecular forces w55x, and adsorbed impurities w55,130,151x. The action of impurities and the wall roughness is difficult to control. However, Ratulowski and Chang w151x disproved Chen’s assertion that surface roughness causes the film thickness observed by Schwartz-Princen-Kiss. Teletzke w55x challenged the significance of intermolecular forces in the experiments. No simple hydrodynamics explains the Bretherton experiments, only an empirical treatment of the adsorption kinetics. 3.2. Modifications of the Bretherton theory: bubble train motion Since the hydrodynamic picture of the flow within the Plateau border influences the resulting pressure drop very little, the Bretherton model applies to the case of touching bubbles. Hirasaki and Lawson w44x assumed that the bubble train moves like a piston so that the distance between the lamellae remains the same. The Plateau border is a region of constant curvature, i.e. the Plateau border is similar to the liquid ahead of and behind an isolated bubble, except that the radius of curvature of its meniscus can be less than the pore radius. The characteristic length of hydrodynamic perturbations becomes of order the size of the Plateau border, much smaller than the length of the bubble. Within uniform regions, the wetting film does not move but resides at the corresponding gas pressure Pi Žsee Eq. Ž37... Because of the difference in the gas pressures within adjacent bubbles, the resulting pressure drop acts on the liquid, which flows from one bubble to another, ‘bearing’ the lamella. As a result, the total dynamic pressure drop per lamella K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 161 Fig. 21. Ža. Endless and Žb. finite air bubbles pushing a liquid slug. Žc. Layout of experiment. Žd. The variation in W s 2 hϱ rR with capillary number: the solid line corresponds to Bretherton’s theory. The dashed line shows the empirical correlation. W s Ca.1r 2 w149x. becomes proportional to Ž Ca. 2r3 rrc , where rc is the radius of curvature of the interface at the Plateau border. Then, applying the Hagen᎐Poiseuille formula to the bubble train, the pressure drop can be compared with that required to maintain steady flow of the wetting fluid within the same pore at the same average speed. The apparent viscosity of foam relative to the liquid, or the ratio of the respective 162 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 pressure drops, is of order NCay1r3 where N is the number of lamellae in the train. For typical reservoir conditions, Ca f 10y6 to y10y3 hence the apparent viscosity of foam is a factor of 10᎐100 N larger than for pure liquid. Ratulowski and Chang w131x, using the same Eq. Ž40., extended the Hirasaki᎐Lawson results by accounting for the difference in the radii of curvature for the forward and backward menisci. Their numerical study, verified by matched asymptotic expansions, showed that at low capillary numbers the pressure drop across each bubble in the train is identical to that of an isolated bubble wEq. Ž42.x; where the radius of the Plateau border, rc , is used instead of R .. That is, the bubbles in the train do not ‘feel’ one another when the train creeps as a whole through the pore. The result is more apparent than real, because the model employed neglects many physical processes at the contact line. The hydrodynamic picture of the flow within the Plateau border is quite complex. However, Park and Homsy w152x showed theoretically that the integral characteristics of the flow, such as the total pressure drop and the ‘asymptotic’ film thickness, are appropriately predicted by the Bretherton theory. The next contribution to these integral characteristics would be of order of O Ž Ca. w152x. On the other hand, in the range of capillary numbers appropriate for the Bretherton theory, Ca ) 10y4 , Hirasaki and Lawson w44x experimentally showed that the pressure drop per bubble exceeds the Bretherton prediction by about one order of magnitude. They attributed the discrepancy to the effect of the surface tension gradient. Recall that the surfactant is adsorbed along the fluid interfaces, where it lowers the interfacial tension. Convection reduces the surfactant concentration at the front of the bubble and increases its surface concentration near the stagnation zone where the flow converges. Adsorption᎐desorption and the bulk diffusion diminish gradients in surface concentration, and, if either of these fluxes are low, a non-uniform distribution of the adsorbed surfactant causes a gradient in the interfacial tension. The resulting Marangoni stress w139x is directed from regions of lower surface tension toward regions of higher tension and causes interfacial convection. Herbolzheimer w153x and Chang and Ratulowski w151,154x demonstrated theoretically that the upper estimate for the film thickness and the pressure drop resulting from the Marangoni effect is 4 2r3 times the Bretherton expressions ŽHirasaki and Lawson w44x took these constants as fitting parameters.. Studies on the behavior of soluble surfactants at various convective, diffusive and sorptive kinetic timescales w151,155x and the bubble length w156x, indicate that the pressure drop for long bubbles are always bounded by a maximum predicted by the Herbolzheimer᎐ Chang᎐Ratulowski model. Thus, this effect cannot explain the experimental results of Hirasaki and Lawson. In biological applications, another type of ‘dynamic’ screening by foam could play an important role, the flexibility of channel walls w157x. A moving bubble creates convergent and divergent regions. which lead to a critical pressure drop. The start-up pressure drop must be distinguished from that required to keep the foam moving, as in Section 2.9. Gaver et al. w158x in the spirit of the Bretherton theory, demonstrates that the critical ‘dynamic’ pressure drop K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 163 required to keep the bubble moving steadily corresponds to the lower boundary for admissible pressure drops and serves as a solvability condition for the Bretherton traveling-wave solution. Gaver et al. w158x termed this critical pressure drop, the start-up, yield᎐pressure drop. Thus, the Bretherton theory includes a number of different physical effects, but, because of its asymptotic character, does not elucidate the interaction between the lamella itself and the flowing liquid within the Plateau border. It implies that each lamella in the train moves steadily, which is only one possible scenario of lamellar motion. Other schemes, e.g. motion of lamellae with acceleration, or with formation of a ring vortex in the Plateau borders w45x are possible. However, to describe these flows, a rigorous mathematical formulation of the problem is needed. 3.3. Sailboat model As discussed in the previous sections, the line of contact of the lamella and the wetting film plays a role of a three-phase contact line. During motion, all geometrical parameters, such as the lamella thickness, the wetting film thickness, the inclination angles and the curvatures of the menisci at the Plateau border, are mutually consistent. Recall that at thermodynamic equilibrium, contact conditions, e.g. ‘free film᎐wetting film’, can be found by applying ordinary thermodynamic rules. However, hydrodynamics also plays an important role, since the lamella exerts a time-dependent tensile force on the wetting film surface. Therefore, the determination of the shape of the moving wetting film is a free boundary problem. We have to find two functions ᎏ the thickness of the wetting film behind and ahead of the contact line, and the time-dependent coordinates of the contact line, by accounting for the coexistence conditions of the free and wetting film. To approximate the Bretherton approach, we use the lubrication approximation wEq. Ž36.x. However, in order to guarantee the existence of the equilibrium wetting film at the pore wall, we modify Eq. Ž37. by introducing the disjoining pressure ⌸ w Ž h. via the equation, Pw s Pg y 2 H y ⌸ w Ž h . . Ž 43 . Recall that 2 H is the curvature of the interface. For a flat substrate Ži.e. in the asymptotic limit hrR ª 0, where R is the pore radius.: 2H s d2 h d x2 1q dh 2 3r2 ž ž // dx . Ž 44 . Thus, Eqs. Ž36., Ž38., Ž43. and Ž44. describe two film profiles h i and h iq1 associated with the i-th and Ž i q 1.-th bubbles of the train. These functions have to match at the three-phase contact line. The free boundary, i.e. at the three-phase contact line, has five boundary conditions, because the curvature in Eqs. Ž43. and Ž44. is expressed through the second derivatives of the fields, and Eqs. Ž36. and Ž38. give two additional derivatives. 164 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 3.3.1. Continuity condition The first boundary condition is the continuity condition: h i < x ª y0 s h iq1 < x ª q0 Ž 45. where is the coordinate of the contact line. 3.3.2. Force balance Local mechanical equilibrium gives the second boundary condition, with restrictions dictated by the geometry. Denote i as the inclination angle between the lamella and the i-th air-wetting fluid interface ŽFig. 22.. Then the angle between the interfaces i and i q 1 will be: s 2 y ␣ i y ␣ iq1 . Ž 46 . Since the profile of the i-th film is increasing, while h iql decreases with x, the angle can be expressed as; s  iq1 y  i , tan i s Ѩh i Ѩx . Ž 47 . xs Thus, Eqs. Ž46. and Ž47. connect the inclination angles with the first spatial derivative, i.e. the slope of the film profile at the contact line. On the other hand, the inclination angles depend upon the lamella tension, as follows from the conditions of mechanical equilibrium: ␥d s i cos Ž y ␣ i . q iq1cos Ž y ␣ iq1 . y F5 , Ž 48 . FH qi sin Ž y ␣ i . s iq1 sin Ž y ␣ iq1 . , Ž 49 . where ␥d is the dynamic tension of the lamella and i is the surface tension of the wetting film in the i-th bubble. A specific friction of the contact line causes the extra-force F Ž F5 , FH ). The latter is unknown beforehand and must be found. Generally speaking, the surface tensions in adjacent bubbles can be distinguished by virtue of different ‘dynamic’ concentrations of the surfactant. But if we neglect Fig. 22. Schematic of a lamella sliding over a wetting film. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 165 this concentration effect, the surface tensions in different bubbles become equal to one another. As shown in Section 2.6, the equilibrium lamella tension is a function of the film thickness. However, the ‘dynamic’ lamella tension differs from its thermodynamic value because of flow occurring within the lamella. 3.3.3. Lamella thickness and dynamic tension To demonstrate the effect of dynamic tension, we assume for simplicity that during the motion, the lamella bulges uniformly like a spherical cap. We also assume that the lamella does not change its volume, but fluid redistributes over a new area by an elongation flow. The first assumption permits a description of the flow in spherical coordinates. The second assumption imposes a diagonal stress tensor with non-zero components, in addition to the hydrostatic pressure, r r , , . Together all tensor components are r-independent functions of the film thickness. Following Section 6 instead of Eqs. Ž10. and Ž12., we write the balance of normal stresses at both outer lamella surfaces as: y Pi q ž ž 2 R i iq1 y Piq1 y / q ⌸ Ž h iiq1 . s yPw q r r , 2 R iiq1 / q ⌸ Ž h iiq1 . s yPw q r r . Ž 50 . Ž 51 . Eqs. Ž50. and Ž51. give the mean radius of the lamella from the ordinary Laplace equation such as Eq. Ž13., with the lamella tension and the pressure difference in Eq. Ž13. treated as dynamic values, because they differ from the equilibrium values by an amount lost due to lamella motion. Eq. Ž13. describes the motion of individual lamella, provided that the radius R iiq1 is expressed through the lamella chord radius, R y h i Ž ., and the apparent contact angle as: sin R y hi Ž . s 1 2 Ž Piq1 y Pi . Ž 52 . Eqs. Ž50. and Ž51. relate the lamella thickness, the disjoining pressure and normal stress as: ⌸ Ž h i iq1 . s 12 Ž ⌬ Piq1 q ⌬ Pi . q r r , ⌬ Pi s Pi y Pw Ž ,t . . Ž 53 . At this step, the rheological model for the lamella has to be specified. To date, no conventional rheological model for foam lamellae exists, because the structure of the foam films remains enigmatic. A smectic-like structure w37,159x, a cubic lattice of ordered micelles w160x, a fluid with a specific exponential correlation function w170x, a bilayer of surfactants with aqueous core w161,162x and a gel-like structure w163᎐165x are only a few of the possible structures, each requiring a K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 166 specific approach to the rheology. Specific surface properties of the foam films also modify rheological relations w166᎐169x, so designating the form of the normal stress is difficult. Nevertheless, many rheological laws permit the expression of the normal stress, r r , in terms of the extensional strains rates, e r r and e . Due to the definition, they are e r r s d ln h iiq1rdt, e s d lnŽ R y h i Ž ..rdt. Then the normal stress, r r , can be expressed through the time derivative of the lamella thickness and chord radius. For instance, assuming that the lamella deforms as a Newtonian fluid, we have the linear relation between the respective component of the stress tensor and the extensional strain rate, i.e. r r s l d ln h iiq1rdt, where l is the lamella viscosity. Thus, the lamella thickness implicitly depends upon the dynamic pressure drop, the magnitude of the pressure within the lamella and the height of the crest. However, the pressure within the lamella must be equal to the pressure in the wetting fluid. Eq. Ž53. relates the parameters of the wetting film and the lamella thickness. Therefore, giving the functional form of the disjoining pressure, the lamella thickness can be found at any instant of time as a function of the input parameters by solving Eq. Ž53.. Finally, the dynamic tension of the lamella can be found from the balance of tangential forces assuming that the film surface tension is approximated by Eq. Ž10. for flat films: ␥d s 2 q ϱ Hh ⌸ Ž h . d h y h iiq1Ž Pw Ž . y Piq1 . q h i iq1 . Ž 54 . iiq1 Here we choose h iiq1 Piq1 as a ‘reference or background tension’. Substituting the difference Ž Pw Ž . y Piq 1 . from Eq. Ž51. into Eq. Ž54. and omitting the term of order h r R, we arrive at the formula for dynamic tension: ␥d s ␥ q h iiq1Ž y r r . . Ž 55 . The above analysis demonstrates how the lamella rheology influences its dynamic tension, e.g. for a Newtonian fluid, the theta-component of the stress tensor is a linear function of the respective component of the extensional strain rate, e . Due to fluid incompressibility, e is a function of the two others, i.e. e s ye r r y e , and can be expressed through the film thickness and the lamella radius. All the lamella parameters, such as thickness, apparent contact angle and tension, are self-consistently connected with the wetting film parameters. 3.3.4. Pressure continuity The lubrication approximation implies that the pressure in the wetting fluid at the contact line varies continuously. Making use of Eq. Ž43., the desired boundary condition is written as: Pi y 2 H N x ª y0 s Piq1 y 2 H N x ª q0 . Ž 56. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 167 3.3.5. Mass balance The last boundary condition is the balance of local mass fluxes at the moving boundary: if the lamella is passing the distance d for the time interval dt, the mass change, on the one hand, is Ž h i Ž . y minŽ h i , h iql . N N x N ªϱ .d and, on the other hand, this change is the flow influx qdt. Hence the boundary condition is, d Ž h i Ž . y min Ž h i ,h iq1 . < x <ªϱ . dt s q < x ª y0, q < x ª y0 s q < x ª q0 . Ž 57. Here h i N x ªyϱ is the thickness of the film at infinity, an unknown parameter, so Eqs. Ž36., Ž38., Ž43. and Ž57. are subject to additional boundary conditions at the ends of the wetting films. At the end where the film remains at rest, the film thickness has a well-defined value. At the same time, the thickness of the residual film at the opposite end is unknown and has to be found. We assume that all the spatial derivatives at both ends vanish, i.e. the wetting films tend to be flat at the ends. Also, we have to presume some law for fluid depletion in the system ‘lamella q wetting film’ 3.4. Sailboat model: tra¨ eling wa¨ e solution The model formulated in the previous section constitutes the basis for a self-consistent theory of lamella residence and motion. We demonstrate its characteristic features by analyzing the traveling-wave solution of form, h s hŽ s . , s s x y Ut, where is the velocity of the lamella propagation over the wetting film. This solution can describe the steady motion of an individual lamella over the prewetted pore, as well as the steady motion of the i-th lamella of a bubble train, when the latter is moving as a whole. The first integration of the governing equations with respect to S, leads to: y 3U h 3m Ž h m y h" . s d ds Pw Ž h m . , m s i ,i q 1, Ž 58 . where ‘q’ is the film thickness at left infinity and ‘y’ at right infinity. The integration of Eq. Ž58. relates between the pressure drop and the lamellar speed y3U ž 0 Hyϱ h y hq h3 ds q ϱ H0 h y hy h3 d s s yPi q Piq1 q ⌸ w Ž hy . y ⌸ w Ž hq . . / Ž 59 . On the other hand, the mass balance, Eq. Ž57., reveals that the residual film thickness and the film thickness at right infinity are the same, hqs h y. At low speed, the displacement can be classified as quasi-static, and the film profile is governed by the capillary and the disjoining pressures. The residual film thickness K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 168 hq and the apparent contact angle have essentially the same values as if the meniscus was motionless, hqs h yf hϱ , where hϱ . is the equilibrium film thickness at infinity. As a result, to leading order in U we have the linear ŽNewtonian. friction law: 3app U hϱ s ⌬ P, Ž 60 . where ⌬ P is the dynamic pressure drop for a steadily moving lamella. The apparent viscosity, app , is expressed through the thermodynamic parameters of the equilibrium film as: ϱ app s 2 H0 Hy1 H 3 d X,H s h hϱ , Xs x hϱ . Ž 61 . In particular, for the Derjaguin᎐Landau approximation to the disjoining pressure w171x, ⌸ w s A r h 2 , A s 2.36 = 10y12 J my2 , the equilibrium shape of the Plateau border can be described analytically w91x. Therefore, the apparent viscosity can be completely investigated as a function of the physical input parameters. For ordinary black films, the contact angle 2 0 , is very small w8,9x, even vanishing. Then the principal term in the asymptotic expansion of the apparent viscosity with respect to  DL s Arhϱ < 1 , has the form w172x: app A 1 ' DL qO ž ' / . DL Ž 62. The thicker the equilibrium film, the greater the apparent viscosity of the foam. Since the equilibrium film thickness is indissolubly connected with the size of the Plateau borders, the growth of the film thickness implies an increase in the size of the Plateau borders. Therefore, the bigger the border, the greater is the effective frictional area. The same tendency might be expected for other forms of the disjoining pressure isotherms because the transition zone makes the main contribution to the apparent viscosity of a lamella. However, the overall foam resistance increases when the foam becomes dryer. For example, substituting Eq. Ž62. into Eq. Ž60., we arrive at the following scaling relation ⌬ P A Ur6hϱ . A linear relation between the pressure drop and the lamella velocity has been observed, but not explained by Nutt and Burley w45x. The two limiting regimes of the lamella friction are Newtonian friction and non-Newtonian friction Žthe power-law Bretherton model.. The Bretherton regime selects a unique film thickness, over which the lamella can only propagate steadily, with no other traveling-wave solutions w131x. As the displacement speed increases, the apparent contact angle between the substrate and bubble interface decreases to zero, when a visible film is deposited, regardless of the wettability of the surface w55,173x. A similar change of regimes of motion might be expected for foam K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 169 lamellae as well, but remains an open question. A detailed experimental analysis of foam motion in smooth capillaries at low capillary numbers, Ca - 10y4 , is highly desirable. 4. Flow trough bamboo-like capillaries and porous media 4.1. Lamella in bamboo-like capillaries: stick-slip motion In the previous sections, we discussed resistance to creeping during steady lamella motion with a constant speed. Creeping occurs only in smooth capillaries. Any wave-like irregularity of the channel produces an additional Laplacian force, which is non-monotonic with respect to spatial coordinates, resulting in an irregular motion of a lamella. After the lamella leaves the pore throat, it bulges forward and resists movement. In the pore enlargement, the lamella bulges backward and accelerates. Visual experiments w46,48x reveal that, when pushed at a constant pressure drop, often the lamellae exhibit stick-slip motion w174,175x. The stick-slip motion can be qualitatively described by using the following model: m d2 i dt 2 s r 02⌬ P , Ž 63. where m is the effective mass of the lamella, t is the dimensional time, and ⌬ P is the pressure drop across the lamella, from the Laplacian pressure drop, Eq. Ž13., the dynamic pressure drop, Eqs. Ž42. and Ž60., and the external pressure drop. Consider a quasi-static motion so that all deviations from Newtonian friction are small. As in Sections 2.6᎐2.9, we assume that the lamella intersects the pore walls at 90Њ, in dimensionless variables, Eq. Ž63. takes the form: 2 d2 i dt 2 q d i dt s ⌬ p y 2 sin i , Ž 64 . where 2 s 2 mhϱ2 Pgr92app r 02 and the dimensionless time is t s 2 hϱ Pg tr3app . ⌬ p is the external pressure drop scaled by Pg . Eq. Ž64. has a very broad range of interpretations, depending on the choice of parameters w111,112,176,177x. In lamella motion, is usually small, so the first term in the left hand side of Eq. Ž64. dominates the second only for short times. At long times, the inertial term may be dropped and Eq. Ž64. takes the form: d i dt s ⌬ p y 2 sin i , ⌬ p s Const. Ž 65 . K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 170 This equation has an analytic solution w178x Žsee also w111x., Ž t . s 2tany1 ž ⌬ p 2 y 4 2 2 ⌬p 2 1r2 t / ž / tan T q 2 ⌬p , Ž 66 . where the period T is Ts 2 '⌬ p 2 Ž 67. y 4 2 2 The characteristic time for translation over a period of the channel depends upon the applied pressure drop. Thus, the periodic motion represented by Eq. Ž66. looks like stick-slip motion. In the vicinity of a threshold, i.e. when the applied pressure drop is of the same order as the Laplacian barrier, ⌬ p f 2 , the lamella spends most time in a slow creep. After overcoming the maximum pinning force Ži.e. after reaching the point s r2., the lamella suddenly jumps over an enlargement and the motion repeats itself. As the pressure drop increases, the jump time becomes negligibly small, and the lamella moves almost steadily. The lamella velocity, averaged over the period T, becomes a non-linear function of the external pressure drop: ²˙ : s 1 T H T 0 d d dt s 2 T s '⌬ p 2 y 4 2 2 , Ž 68 . with the expected threshold ⌬ p f . Eq. Ž68. displays an increase of the apparent viscosity of the lamella with increasing rate of shear, as determined by the lamella speed ŽFigs. 23 and 24., ‘shear-thickening’ w179x. The averaging which leads to Eq. Ž68. somewhat disguises its cause: near the vicinity of the critical pressure drop, capillary forces drive the lamella almost entirely. In this critical regime, a ‘bare’ external pressure drop is needed to compel the lamella to shift from its equilibrium position. Once shifted, the lamella drifts almost autonomously. Only in a high-speed regime does the apparent viscosity of the lamella tend to its bare value. Therefore, this ‘shear thickening’ behavior of a lamella is caused solely by irregularities of the channel. If we increase the pressure drop further and further, we inevitably arrive at the range of validity of the Bretherton theory. In this flow regime, ‘shear thickening’ behavior can change to the ‘shear-thinning’ Ži.e. the reduction of viscosity with increasing rate of shear in an averaged steady flow. w113x. One possible explanation of the experimentally observed ‘shear thinning’ ŽFig. 25. is the non-linearity of the friction law as expected from the Bretherton theory w121x, though the theory, which deals with a uniform capillary and which diminishes the effect of gas compressibility, underestimates the foam viscosity under all experimental conditions ŽFig. 26., and incorrectly predicts dependence of the gas phase velocity. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 171 Fig. 23. Transition from a stick-slip motion to a sliding motion wEqs. Ž66. and Ž67., s 10y3 .. To resolve the contradictions, Falls et al. w121x modified the Hirasaki᎐Lawson theory to account for the contribution of pore constrictions to the apparent viscosity. They approximated the Laplacian contribution to the dynamic pressure gradient as ٌP f n l 4rrt , where n l is the number of lamellae per unit length and Fig. 24. ‘Shear-thickening’ behavior of an individual lamella. The increase of viscosity with the increasing rate of shear in an averaged steady shear flow wsee Eq. Ž68., s 10y3 x. 172 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 25. Typical rheological curves. Foams from 0.25% Hostapur SAS carbon number fractions. The length of the pack, 100 cm; the permeability, 8 darcy; the pressure at the exciting end, 7 bar. ⌬ Pt , denotes the threshold pressure drop; vertical axis, the gas flow rate Žml miny1 .; horizontal axis, the pressure droprthe length of the pack Žbar my1 . w122x. Fig. 26. Apparent foam viscosity in glass bead packs as a function of gas-phase velocity. Physical conditions: the capillary pressure Pc f 1900 dynes cm -2 ; the ordinary gas-phase relative permeability, k r g f 0.056. Ža. The controlled-bubble-size regime, r b r rcap f 3.8, where r b is the effective radius of a foam bubble and rcap is an equivalent capillary radius. Žb. The pack-generated bubble-size regime, r brrcap f 0.92. The parameter X f s S g frS g is the ratio of the moving-gas saturation S g f , estimated from the measured residence time of flowing gas, and the gas saturation, S g . For the pack-generated bubble-size regime, X f s 1; s , is the Hirasaki᎐Lawson viscosity w121x. rt is the radius of the pore throat. As a result, the Laplacian contribution augments the Hirasaki᎐Lawson apparent viscosity defined through the Hagen᎐Poiseuille law as, con f 4 n lr¨ g rt , where ¨ g is the gas-phase velocity. Thus, the apparent viscosity must vary with the y1 power of the gas velocity at low speeds. Falls and K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 173 coworkers studied two regimes of foam generation: so-called controlled-bubble-size regime, where the bubble size was unchanged in situ, and the liquid was transported within the lamellae; and the so-called pack-generated-bubble-size regime, where one bead pack was used to generate the foam for a second. For both flow regimes the apparent viscosity measured at low speeds varies as y1 power of the gas-phase velocity. The authors fit the Hirasaki᎐Lawson theory to the experimental data even at higher speeds. At higher rates, the viscosity varies as the y1r3 power of the velocity in the first flow regime, and changes to the y2r3 power of the velocity in the second, so the surface tension gradient governs the foam resistance in the second flow regime w44x Žsee also Section 3.2.. Thus, the key assumption in treating the experimental data is the contribution of the pore constrictions to lamella resistance, which determines the changes in apparent viscosity as a function of the gas velocity. The results apply only for capillary numbers 10y4 - Ca < 1 and for short bubble trains, when the length of the sample is of the order of magnitude of the correlation length Ncor of the train. ŽThe Falls᎐Musters᎐Ratulowski experiments fall into this range, see Section 2.9.. Otherwise, the chain elasticity cannot be ignored and modifies the theory. 4.2. Stick-slip motion of bubble trains Consider waves caused by oscillations of gas pressure within bubbles w113,132x, which apply to acoustic flows, for which the high frequency modes are important, and transport of bubble trains. We first consider the propagation of acoustic waves through a perfectly ordered one-dimensional foam, i.e. through a bubble chain. We first specify the dispersion relations w116x. For the linearized version of Eqs. Ž64. and Ž20., the dispersion relation has the form Žneglecting the lamella friction.: 2 s 20 sin 2 ž Kk 2 q 2 2 K , / 20 s 2 Pg liquid hK , Ž 69 . where is the frequency, k is the wave number, and h is the lamella thickness. Eq. Ž69. has two critical frequencies, below and above which acoustic waves will not propagate through the chain. The lower critical boundary corresponds to the natural frequency of the individual lamella, reached when the wavelength tends to infinity Žor k ª 0., i.e. when the bubble train oscillates as a whole, and when all the lamellae vibrate at the throats of pore channels in unison. Another critical frequency selects the range of wave numbers allowed by the theory Žwhen the wavelength becomes comparable with the distance between adjacent lamellae, the theory fails to describe the triple substrate ᎐gas᎐lamella interactions .. The theory does not require pore channels with small pore aspect ratios, but is valid for all channels with perfectly ordered foam. Because screening originates from the Laplacian blocking action of the lamellae, we expect that any perturbation will be suppressed if its frequency lies outside the allowed range. Our analysis applies only to waves of small amplitude, and cannot be extended directly to the non-linear case. For the stick-slip motion of lamellae, 174 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 Fig. 27. Wave motion of foam. The head and tail lamellae move in unison, while the middle lamella moves in antiphase. Numerical experiment based on Eq. Ž64. and Eq. Ž20., s 0 w180x. we should focus on the wave-propagating solutions of the non-linear model, Eqs. Ž64. and Ž20., which dynamics occurs when we excite just one lamella, with all others initially at rest. In its general features, the stick-slip motion of the lamellae resembles the Frenkel᎐Kontorova mechanism for the flow of dislocations in crystals w107x. The excited lamella initiates a wave propagating through the bubble chain in a ‘falling dominos’ type of motion towards the other end of the chain, which is then reflected, moving in the opposite direction towards the initially excited lamella ŽFig. 27.. In such a motion, practically no energy is lost during the flow, because the surface energy of the lamellae is compensated almost entirely by the elastic energy of the bubbles. Focusing on the short time dynamics within an interval, we find a solution like the Frenkel᎐Kontorova kink w107,181,182x. Just after the external pressure drop has overcome the threshold G wEqs. Ž33. and Ž35.x, the domain wall, which separates the stagnation zones ahead of and behind the front of the wave, deepens and runs forward. Therefore, in soliton dynamics, the competition between the inertial, elastic, and pinning forces is important, and the viscous forces do not crucially change the character of the soliton propagation over short periods. If a bubble train moves for a long time, the inertial forces are negligibly small with respect to the three major forces ᎏ the elastic, viscous, and capillary forces. These forces together drive the displacement waves. If the external pressure drop is maintained for a long time, the number of the domain walls increases proportio- K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 175 nally to the number of channel periods over which the front domain wall has skipped. Therefore, the longer the bubble train, the greater the number of domain walls, and, consequently, the greater the bubble train resistance. To illustrate collective effects in bubble train friction, we use a simplified model of lamellar interactions. A linear elastic spring represents the elastic interactions of the lamellae. Because each individual domain wall plays a role in the chain friction, we focus on the motion of a solitary wave. Then, in the continuum limit, Eq. Ž65. is w113x: Ѩ Ѩt s Ѩ2 Ѩs 2 y K sin q ⌿. Ž 70 . Here ⌿ is a dimensionless external force which represents the total action of the external pressure gradient on each lamella in the chain. In the vicinity of the threshold, the profile of a domain wall can be adequately approximated by Eq. Ž32., provided that the condition G f 6 < 1 holds. Then, rewriting Eq. Ž70. in a wave-fixed reference frame, the velocity of the domain wall is found as the solvability condition to Eq. Ž70. w113,183x, y¨ ϱ Hyϱ d ž / dx 2 dxs⌿ ϱ d Hyϱ d x d x, Ž 71. where x s 6Ž s y ¨ t . and the function in Eq. Ž71. is expressed by Eq. Ž32.. Substituting Eq. Ž32. into Eq. Ž71., we have: ¨A ⌿ ' . Ž 72 . Thus, Eq. Ž72. may be treated as an equation for the motion of effective individual lamellae along an active channel. The friction coefficient of such a macro-lamella is affected by the binding and pinning energies of the chain. Braiman᎐Family᎐Hentschel w184x and Musin’s w180x numerical experiments indicate that the friction coefficient grows with the number of lamellae in the chain Žthe number of moving domain walls. and the scales similarly to Eq. Ž72... 4.3. Weak foams: flow of ‘solutions’ of bubble chains through porous media As shown in Section 4.2, during motion, a bubble train is subdivided into the regions within which the lamellae are pinned at the narrowest parts of the channel so that the most of the dimensionless displacements, i , remains almost constant. The size of each transition zone connecting adjacent regions is much greater than the pore size. Therefore, on the scale of the sample, the train behaves as a one-dimensional coarse-grained foam in which the transition zones serve as macrolamellae. Such lamellae interact between themselves via a recommended elastic potential, which incorporates all the microlevel elastic and pinning interac- 176 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 tions and displays the collective properties of the train w181,182x. The distance between neighboring effective lamellae crucially depends upon the history of loading. For instance, if a piston is applied to the tail end of the chain, and such a load is maintained long enough to shift the head lamella of the chain, then the number of domain walls will be a well-defined constant, characteristic of the bubble train, the pore channel and the applied pressure drop. The resulting structure of the chain on the scale of the sample will be periodic, i.e. macrolamellae will form a perfect lattice. However, if the load continuously redistributes over the train Žthe situation most likely in a natural porous medium., then the chain generally acquires an irregular structure: in the coarse-grained foam, the cell size varies. Designating the bubble train as an object allows a hydrodynamic theory for a weak foam, which as a very small fraction of trapped bubbles, so that the gas and bubble trains flow together. Gas mobility is reduced due to enhanced friction of the bubble trains. The hydrodynamics of polymer solutions is similar w185x; the action of bubble trains on a free gas, flowing through a porous medium, resembles the action of polymer molecules on a flowing solvent. Because the bubble train resistance is accumulated almost entirely within the domain walls, and, at the same time, each domain wall is much longer than a pore size, the hydrodynamics of a weak foam resembles that of an ensemble of bubble trains ‘dissolved’ in a Darcian fluid, provided that the free gas flow obeys Darcy’s law Ži.e. the relation between the pressure gradient and the gas velocity is linear. In a non-uniform flow, where each bubble train feels an external pressure gradient imposed by the free gas, any train will be stretched or contracted Žbecause each macrolamella in the train acquires a velocity imposed by kinematic conditions.. The inherent elasticity of the train tends to restore its equilibrium length, moving the surrounding gas and, consequently, other trains. Therefore, trains which are stretched only slightly flow to regions in which chains are more extended, creating a net restoring force. A schematic is presented in Fig. 28, where a dumbbell plays the role of a foam filled macrobubble or bubble train. This scenario for foam flow can be described using Eq. Ž72.. Then, phenomenologically expressing ⌿ in Eq. Ž72. as the pressure gradient, the flux of foam can be written as w113,132x.: Fig. 28. Sketch of the creation of a restoring force in a gradient flow. Trains, which are more extended, move the surrounding gas, pulling other trains. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 l J s yc H0 d s - uŽ s,t . ) иٌP. 177 Ž 73. Here c is the concentration of the bubble trains, l is the length of the train, and  is a phenomenological constant. In Eq. Ž73., uŽ s . a vector which joins two adjacent macrolamellae, and arclength s is measured along the train. The angle brackets denote an average over the orientations of the macrobubbles which link the macrolamellae. It is convenient to rewrite Eq. Ž73. in the following form: J s J0 y c l H0 d sS и ٌP , J0 s yc  3 Ž 74 . ٌP , S s ²u Ž s,t . u Ž s,t . y 13 I: , Ii j s ␦ i j , S и u s Si j u j , S:u s S i j u i u j , where S is the tensor of the order parameter. The first term on the right hand side of Eq. Ž74. may be included in the ordinary Darcy’s law associated with the flow of free gas. Thus, Eq. Ž74. conveniently treats foam motion as a flow of a ‘solution’ of bubble chains w185x. The detailed derivation of the constitutive equations has been presented by Kornev and Kurdyumov w132x. They contain a generalized Darcy’s law for gas flow in the presence of a foam: v sy k f ٌP y L,L s J y J0 , Ž 75 . and a kinetic equation which describes an evolution of the ‘blocking force’, L, due to the action of velocity gradients. In a one-dimensional case, e.g. for foam motion in a porous tube, the model has the following form: Ѩ¨ Ѩx ѨL Ѩt s 0, q¨ ¨ s yk f Ž 76 . ѨL Ѩx ѨP Ѩx sy y L. L , Ž 77. Ž 78 . Here ¨ is the gas velocity, a relaxation time, and k f an effective seepage coefficient. Eqs. Ž76. ᎐ Ž78. may correspond to the lamella ‘break-and-reform’ mechanism of foam motion w186x. Eq. Ž78. balances the forces of the moving foam. The viscous force on the left hand side of the equation balances the pressure gradient and blocking force, which arises from blockage of gas channels by ‘valve’ K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 178 lamellae. In the first approximation, the blocking force is directly proportional to the number of ‘valve’ lamellae. We thus can express the lamella density through the function L and, ignoring the generation rate of lamellae, treat Eq. Ž77. as a kinetic equation for the lamella population balance in a moving foam. The parameter should thus be treated as the lifetime of a ‘valve’ lamella, and represented as: 1 s 1 t q 1 h , where t is the ‘thermodynamic’ lifetime of a lamella, and h is its hydrodynamic traveling time. So, each ‘valve’ lamella blocks a gas path until it rips due to inherent thermodynamic or hydrodynamic instability. Within the bubble-train-flow-mechanism, the relaxation time can be obtained by analyzing a random walk of the bubble train as a whole. We hypothesize that once formed, the active channels cannot be destroyed, so the random walk resembles reptation of a polymer molecule through obstacles w132,185,187x.. A bubble train can change its active channel only by moving through a network of active channels in a worm-like manner. In fact, only the ends move with the rest of the macrolamellae effectively trapped within the currently existing active channels. Eqs. Ž76. ᎐ Ž78. reproduce the main features of the flow of a foam through a porous medium ŽFigs. 29 and 30.. The one-dimensional flow is governed by the parameters and b s L0 k f rH, where L0 is a boundary value of the blocking force L, and H is the length of the sample. These parameters can be extracted from experimental data using a piecewise-linear approximation ŽFig. 29.: kf kf ⌬P H ⌬P H s Ž1 q b . ¨ , s¨ q H b, ¨ ª 0, Ž 79 . ¨ ª ϱ, Ž 80. where ⌬ P is the applied pressure drop. The point at which the straight lines intersect has the coordinates ⌬ P k frH 2 s 1 q b Ž 1 y ey1 . f 1 q b. ¨ rH s 1, Ž 81 . The coordinates of the intersection point can be used as fitting parameters. An additional parameter is the breakthrough time T U s tU , which the front of the foam first reaches the exit of the sample. The dimensionless parameter tU can be found from the transcendental equation 1s ⌬ P k f H2 ½ tU 1qb q b 1qb и 1 y exp Ž y Ž 1 q b . tU . 1qb 5 . Ž 82 . K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 179 Fig. 29. The dimensionless velocity V s ¨ rH as a function of the dimensionless pressure drop ⌬P s ⌬ P k frH 2 . 1. y b s 1, 2. y b s 3, 3. y b s 5, 4. y b s 10. Eqs. Ž79. ᎐ Ž82. are closed, so that b, and k f can be found for each experiment, and plotted as functions of the water saturation, properties of the foaming agent, etc. Despite the small number of physical constants, the theory cannot directly apply to two- or three-dimensional flows Žbecause the parameter b depends upon the boundary value of the order parameter S.. A similar problem occurs in polymer viscoelasticity w188x and remains unresolved. Fig. 30. The dimensionless breakthrough time as a function of the dimensionless pressure drop ⌬P s ⌬ P k frH 2 . 1. y b s 1, 2. y b s 3, 3. y b s 5, 4. y b s 10. 180 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 4.4. Problems in describing strong foams The main feature of foam theology is ‘pseudoplasticity’, which is very sensitive both to the foam texture inside the porous medium and to the fraction of pores blocked by lamellae. The physics of blocking by weak and strong foams differ. For a weak foam, even in the absence of a start-up pressure gradient, enhanced friction of bubble trains results in apparent pseudoplasticity. For a strong foam, when the bubble trains flow through a small fraction of the pore channels, blocking originates both from trapped lamellae which do not participate in the movement, leading to a permeability orders of magnitude lower than for single phase flow, and from an enhanced friction of the bubble trains. Thus, the total gas mobility, i.e. the coefficient in the generalized Darcy’s law, is the product of permeability divided by apparent viscosity. If the liquid also participates in the movement so that its saturation alters the fraction of trapped pores, then, instead of the absolute permeability, the relative permeability to gas has to be considered. In addition to non-linearity of the gas mobility, the start-up pressure gradient required to deepen some bubble trains and to create a network of active channels w189᎐192x plays its key role to the hydrodynamics of strong foams. Thus, strong foams are plastic. The bubble trains and networks of active channels correspond to the dislocations and their network of dislocations in solid state physics. The trapped gas saturation of strong foams in porous media can be as high as 80% w17x. Therefore, the trapped foam forms an elastic field around the network of active channels. For short times, the trapped foam consists of motionless lamellae. However, diffusion affects long-term behavior w193,194x. The Nabarro᎐Herring᎐Lifshitz mechanism w195᎐198x for diffusion-induced plasticity may be important to the motion of trapped foam. Under a pressure gradient, gas diffuses through the lamellae so the gas pressure within the bubbles changes with time. The difference between the Laplacian and gas pressures causes the trapped lamellae to creep. Besides flowing foam’s maintenance of capillary pressure near the limiting value w67x ŽSection 2.3., the reduction of gas mobility suggests that, during flow, the network of active channels remains nearly in the same state, as if at the percolation threshold w120x. Scaling estimates of the gas mobility reduction by foams using ordinary percolation theory w117,120,104x, assumed the critical behavior of the network of active channels ad-hoc. Though most researches believe that transport in strong foams obeys percolation laws, the hydrodynamics remains unclear and puzzling. 4.5. Foams in fiber systems Foam applications in industrial fiber technologies include dyeing, printing, mercerizing, and finishing of textile fabrics, paper coating, and resin-impregnation of fibrous mats and fabrics w199᎐204x. Using foams instead of bulk liquids to deliver small amounts of liquid solutions to fiber surfaces, leads to substantial energy savings because of the small amount of residual solvent to be removed at K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 181 drying. Bubbles also occur in fabricating fiber-reinforced composites, during liquid resin impregnation, damaging uniform polymer distribution w205x. Studies try to minimize air entrapment and bubble interactions in a fiber network w206,207x. Interactions of bubbles with fibers and fiber networks have certain specifics compared with capillaries and pore networks in solid-wall materials. The major difference is that in contract with granular materials, where the sizes of pores and grains are commensurable, the typical diameter of fibers, which constitute a skeleton of a pore structure, is commonly smaller than the typical diameter of voidsrpores between the fibers. Therefore, identifying single pores, their shapes and dimensions is difficult. The most rational way to introduce the pore dimensions is to consider a model system of ‘effective’ pores, in which some characteristic processes would agree with experiments. For example, effective pore sizes are estimated from experiments with capillary equilibrium of immiscible fluids Žcommonly, wetting liquid and gas. and with steady or quasi-steady forced flow of a non-wetting fluid Žcommonly, gas. through a fibrous sample. An advanced technique of liquid porosimetry has been developed by Miller and Tyomkin w208x for determining pore volume distributions in fiber systems and other materials. The method is based on the consecutive, quasi-equilibrium wetting fluid᎐gas displacement under precisely controlled pressure in an isothermal environment. The effective radius of a pore, where the liquid and gas phases coexist at a given pressure, is related to the mean radius of curvature of the equilibrium meniscus between the phases through the Laplace equation. During gas᎐liquid displacement in fibrous materials, bubble formation occurs mostly due to the snap-off and leave-behind mechanisms. However, the fiber structure causes some peculiarities. Mostly, we deal with strong foams with the bubble size commensurate to the pores, as follows from the experiments of Gido et al. w209x, who examined the flow of foams through fibrous mats by characterizing the bubble sizes before and after injection. The size of the output bubbles exiting the fiber system was independent of the bubble size of the input bulk foam, but correlated with the pore size distributions measured by liquid porosimetry w208x. This result can be interpreted assuming that the bubbles during foam flow take into two configurations: immobilized bubbles strongly pinned by intersecting noncollinear fibers and unpinned bubble trains sliding along the active channels confined by the immobilized bubbles. The mobile bubbles are comparable in size to the pore constrictions in order to pass through them without essential deformation. The movement of lamellae in the active channels within a fibrous structure should be like the lamella stick-slip motion within porous solids. However, foam flow through a fiber mat is more difficult to formalize. We also have to account for wetting behavior of films on fiber surfaces. The lamellae Žbubble films. should coexist with wetting films and Plateau borders at the intersections of lamellae andror lamellae and wetting films. The equilibrium configurations are determined by capillary and disjoining pressures. Fiber surfaces are commonly convex, so the capillary pressure acting on the liquid᎐gas interface of the wetting film covering the fiber is positive and tends to squeeze liquid out of film into the regions of fiber crossings, leading to a reduced mobility of wetting 182 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 films on fibers compared with capillaries. Liquid spreading and drop residence on fibers were considered first by Carroll w210,211x, who accounted for capillary forces only, and then in great detail by French researchers w212᎐216x. Brochard emphasized the central role of long-range intermolecular forces in residence, stability, and spreading of films and drops on fiber surfaces w212,213x. Meglio w214x experimentally observed wetting films stabilized by van der Waals forces, and proved that mass transfer between the drops residing on a fiber occurs through these films. Similar effects should affect bubbles on fibers and in fiber networks. A theory of foams in fiber systems requires studies of the equilibrium shape of a bubble residing on a fiber, the transition zone between the bubble lamella and the wetting film coating the fiber, the slippage of a bubble along a fiber, bubble pinning at a fiber crossing, motion of a train of bubbles transfixed by a fiber Ž‘bubbles on a spit’., etc. 5. Conclusions Foams in porous media have unique structural, rheological, and transport properties determined by interactions between foam films Žlamellae., wetting films, and solid surfaces. The interfacial interactions in confined geometry impose specific restrictions on foam patterning and flow compared with bulk foams. Despite numerous applications of foams in industrial processes, foam behavior in porous media remains enigmatic, and a better understanding of physical mechanisms governing lamellae distribution and motion in pores is appealing. Thermodynamic and hydrodynamic properties of thin films in pores are determined by interplay of capillary ŽLaplacian pressure. and surface Ždisjoining pressure. forces. We reviewed recent experimental studies and theoretical models from the unified point of view, considering foam lamellae and wetting films as an interacting system. A consistent description of conditions of mechanical equilibrium of curved lamellae, including dynamic effects, is presented for the first time. Theoretical models of foam patterning under a load and foam friction show that the binding forces caused by bubble compressibility and the pinning forces due to capillarity determine specific ordering of lamellae in porous media and anomalous foam resistance. A new microscopic bubble train model is proven to predict asymptotic expressions for the start-up-yield pressure drop. The Bretherton theory of the forced fluid᎐fluid displacement as applied to bubble transport through pore channels is augmented based on a new sailboat model, which accounts for thermodynamic coupling of foam lamellae and wetting films. Stick-slip motion of lamellae and bubbles in pores of varying diameter is described in terms of the bubble train model. However, there exits a gap between the macroscopic models currently available for practical engineering applications and the pore level understanding of the physical mechanisms involved. Much work is required to turn from the microlevel to the pore network level and to construct macroscale models of foams in porous media, such as reservoirs, granular, and fibrous materials. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 183 References w1x C.V. Boys, Soap Bubbles and the Forces which Mould Them. Society for Promoting Christian Knowledge, E. and J.B. Young, London, 1890. Reprinted in Doubleday Anchor Books, New York, 1959. w2x K.J. Mysels, K. Shinoda, S. Frankel, Soap Films, Studies of their Thinning and a Bibliography, Pergamon Press, New York, 1959. w3x J.J. Bikerman, Foams, Springer-Verlag, New York, 1973. w4x I.B. Ivanov, Thin Liquid Films: Fundamentals and Applications, Marcel Dekker, New York, 1988. w5x A.M. Kraynik, Ann. Rev. Fluid Mech. 20 Ž1988. 325. w6x A. Wilson, Foams: Physics, Chemistry and Structure, Springer-Verlag, New York, 1989. w7x J. Glazier, Ph.D. Thesis, University of Chicago, 1989. w8x P.M. Kruglyakov, D.R. Exerowa, Foam and Foam Films, Khimia, Moscow, 1990. w9x D.R. Exerowa, P.M. Kruglyakov, Foam and Foam Films, Elsevier, Amsterdam, 1998. w10x J. Stavans, Rep. Prog. Phys. 56 Ž1993. 733. w11x R.K. Prud’homme, in: S.A. Khan ŽEd.., Foams: Fundamentals and Applications, Marcel Dekker, New York, 1995. w12x R.J. Pugh, Adv. Colloid Interf. Sci. 64 Ž1996. 67. w13x S.H. Raza, Soc. Petr. Eng. J. 10 Ž1970. 328. w14x S.S. Marsden, Foams in Porous Media-SUPRI TR-49, US DOE, 1986. w15x J.P. Heller, M.S. Kuntamukkula, Ind. Eng. Chem. Res. 26 Ž1987. 318. w16x L.L. Schramm, Foams: Fundamentals and Applications in the Petroleum Industry, Advances in Chemistry Series, 1994, p. 242. w17x W.R. Rossen, Rep. Prog. Phys. 56 Ž1993. 413. w18x T.W. Patzek, in: D.H. Smith ŽEd.., Surfactant Based Mobility Control Progress in Miscible-Flood Enhanced Oil Recovery, ACS Symposium Series 373, Washington DC, Chemical Society, 1988, p. 326. w19x A.H. Falls, G.J. Hirasaki, T.W. Patzek, P.A. Gauglitz, D.D. Miller, T. Ratulowski, SPE Res. Eng. 3 Ž1988. 884. w20x W.W. Mullins, J. Appl. Phys. 59 Ž1985. 1341. w21x M. Marder, Phys. Rev. A 36 Ž1987. 858. w22x H. Flyvbjerg, Phys. Rev. E 47 Ž1994. 4037. w23x J.C.C. Nitsohe, Am. Math. Mon. 81 Ž1974. 945. w24x N. Rivier, in: D. Bideau, A. Hansen ŽEds.., Disorder and Granular Media, Elsevier, New York, 1993. w25x D. Weaire, N. Rivier, Contemp. Phys. 25 Ž1984. 59. w26x J.A. Glazier, D. Weaire, J. Phys. Condens. Matter 4 Ž1992. 1867. w27x F. Bolton, D. Weaire, Phys. Rev. Lett. 65 Ž1990. 3440. w28x F. Bolton, D. Weaire, Philos. Mag. B 63 Ž1991. 795. w29x B.V. Derjaguin, Kolloid Zeits. 69 Ž1934. 155. w30x A.D. Scheludko, Adv. Colloid Interf. Sci. 1 Ž1967. 391. w31x H. Wennerstrom, ¨ B. Lindman, Phys. Rep. 52 Ž1979. 1. w32x F.D. Gunstone, J.L. Harwood, F.B. Padley, The Lipid Handbook, Chapman and Hall, London, NewYork, 1986. w33x E.H. Lucassen-Reynders, Anionic Surfactants ᎏ Physical Chemistry of Surfactant Action, Surfactant Science Series 11, Marcel Dekker, New York, 1981. w34x P.R. Garrett, Chem. Eng. Sci. 48 Ž1993. 367. w35x B.V. Derjaguin, M.M. Kussakov, Acta Physicochim. URSS 10 Ž1939. 25. w36x B.V. Derjaguin, Zh. Fiz. Khim. 16 Ž1940. 137. w37x V. Bergeron, C.J. Radke, Langmuir 8 Ž1992. 3020. w38x K.J. Mysels, M.N. Jones, Discuss. Faraday Soc. 42 Ž1966. 42. w39x D. Exerowa, T. Kolarov, K.H.R. Khristov, Colloids Surf. 22 Ž1987. 171. 184 K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 w40x A.S. Aronson, V. Bergeron, M.E. Fagan, C.J. Radke, Colloids Surf. A: Physicochem. Eng. Asp. 83 Ž1994. 109. w41x K.I. Khristov, D.R. Exerowa, P.M. Kruglyakov, Kolloidn. Zh. USSR 43 Ž1981. 101. w42x A.V. Neimark, M. Vignes-Adler, Phys. Rev. E 51 Ž1995. 788. w43x A.N. Fried, The Foam-Drive Process for Increasing the Recovery of Oil, Report US Bureau of Mines R.I.5866, 1961. w44x G.J. Hirasaki, J. Lawson, Soc. Petr. Eng. J. 25 Ž1985. 176. w45x C.W. Nutt, R.W. Burley, Foams: Physics, Chemistry and Structure, Springer-Verlag, New York, 1989, p. 105. w46x A.V. Bazilevsky, K. Kornev, A. Rozhkov, Rheology and Fluid Mechanics of Nonlinear Materials, vol. 217, ASME AMD 1996, pp. 101᎐104. w47x O.S. Owete, V.E. Brigham, SPE Res. Eng. 2 Ž1987. 315. w48x K.T. Chambers, C.J. Radke, in: N. Morrow ŽEd.., Interfacial Phenomena in Petroleum Recovery, Marcel Dekker, New York, 1990, p. 191. w49x R.A. Ettinger, C.J. Radke, SPE. Res. Eng. 7 Ž1992. 83. w50x A.R. Kovscek, C.J. Radke, Foams in Porous Media-SUPI-TR, US DOE, 1986, p. 113. w51x H.M. Princen, in: E. Matijevic ŽEd.., Surface and Colloid Science, vol 2, Wiley, New York, 1969, p. 1. w52x C. Hun, L.E. Scriven, J. Colloid Interf. Sci. 30 Ž1969. 323. w53x K.K. Mohanty, Ph.D. Thesis, University of Minnesota, 1981. w54x L.I. Kheifets, A.V. Neimark, Multiphase Processes in Porous Media, Khimia, Moscow, 1982. w55x G.F. Teletzke, Ph.D. Thesis, University of Minnesota, 1983. w56x B.V. Derjaguin, N.V. Churaev, Wetting Films, Nauka, Moscow, 1984. w57x B.V. Derjaguin, N.V. Churaev, V.M. Muller, Surface Forces, Nauka, Moscow, 1985. w58x B.V. Derjaguin, N.V. Churaev, V.M. Muller, Surface Forces, Consultants Bureau, New York, 1987. w59x G.J. Hirasaki, in: N.R. Morrow ŽEd.., Interfacial Phenomena in Petroleum Recovery, Marcel Dekker, New York, 1990, p. 23. w60x B. Derjaguin, N.V. Churaev, J.Colloid Interf. Sci. 54 Ž1976. 176. w61x M.M. Viktorina, B.V. Derjaguin, I.G. Ershova, N.V. Churaev, Dokl. Akad. Nauk SSSR 200 Ž1971. 1306. w62x N.V. Chureav, B.V. Derjaguin, I.G. Ershova, Res Surface Forces, 4, Consultants Bureau, New York, 1975, p. 165. w63x M. Sahimi, Rev. Mod. Phys. 65 Ž1993. 1393. w64x J. Bear, Dynamics of Fluids in Porous Media, Elsevier, New York, 1972. w65x F.A.L. Dullien, Fluid Transport and Pore Structure, Academic Press, New York, 1979. w66x L.W. Lake, Enhanced Oil Recovery, Prentice-Hall, New Jersey, 1989. w67x Z.I. Khatib, G.J. Hirasaki, A.H. Falls, SPE Res. Eng. 3 Ž1988. 919. w68x G.I. Barenblatt, V.M. Entov, V.M. Ryzhik, Theory of Fluid Flows Through Natural Rocks, Kluwer Academic Publishers, London, 1990. w69x G. Singh, G.J. Hirasaki, C.A. Miller, AIChE J. 43 Ž1997. 3241. w70x P. Bak, C. Tang, K. Wiesenfeld, Phys. Rev. A 38 Ž1988. 364. w71x P. Bak, K. Chen, Sci. Am. 264 Ž1991. 26. w72x P. Bak, How Nature Works. The Science of Self-Organized Criticality, Springer, Berlin, 1996. w73x J.E. Hanssen, in: S.M. Skjaeveland, J. Kleppe ŽEds.., Recent Advances in IOR Methods, SPOR Monograph, Norwegian Petroleum Directorate, Stavanger, 1992, p. 277. w74x Roof. Soc. Petr. Eng. J. 10 Ž1970. 85. w75x R.S. Schechter, A.C. Lam, M.S. Willis, in: G. Pe’tre’, A. Sanfeld ŽEds.., Capillarity Today, Springer-Verlag, Berlin, 1990. w76x D.H. Everett, J.M. Haynes, J. Colloid Interf. Sci. 38 Ž1972. 125. w77x A.V. Neimark, Kolloidn Zh. USSR 47 Ž1985. 531. w78x A.V. Neimark, A.B. Rabinovich, Kolloidn Zh. USSR 47 Ž1985. 1103. w79x A.V. Neimark, A.B. Rabinovich, Kolloidn Zh. USSR 52 Ž1990. 870. w80x A.V. Neimark, L.I. Kheifets, Kinet. Catal. 23 Ž1982. 1510. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 w81x w82x w83x w84x w85x w86x w87x w88x w89x w90x w91x w92x w93x w94x w95x w96x w97x w98x w99x w100x w101x w102x w103x w104x w105x w106x w107x w108x w109x w110x w111x w112x w113x w114x w115x w116x w117x w118x w119x w120x w121x w122x w123x 185 M.W. Cole, W.F. Saam, Phys. Rev. Lett. 32 Ž1974. 985. V.M. Nabutovsky, V.P. Belosludov, A.M. Korotkikh, Kolloidn Zh. USSR 41 Ž1979. 722. V.M. Nabutovsky, V.P. Belosludov, A.M. Korotkikh, Zh. Eksp. Teor. Fiz. 77 Ž1979. 700. M.P. Gelfand, R. Lipowsky, Phys. Rev. B 36 Ž1987. 8725. W.R. Osborn, J.M. Yeomans, Phys. Rev. E 51 Ž1995. 2053. A.V. Neimark, Langmuir 11 Ž1995. 4183. P.I. Ravikovitch, S.C. O’Domhnaill, A.V. Neimark, F. Schuth, ¨ K.K. Unger, Langmuir 11 Ž1995. 4765. B.V. Derjaguin, G.A. Martinov, Yu.V. Gutop, Kolloidn. Zh. USSR 27 Ž1965. 357. J.A. de Feijter, The Form-Drive Process for Increasing the Recovery of Oil, Report US Bureau of Mines R.I.5866, 1961, p. 1. J.A. de Feijter, A. Vrij, J. Electroanal. Chem. Interf. Electrochem. 37 Ž1972. 9. A.V. Neimark, L.I. Kheifetz, Kolloidn. Zh. USSR 43 Ž1981. 500. M. Kagan, W.V. Pinczewski, J. Colloid Interf. Sci. 180 Ž1996. 293. Z.M. Zorin, D. Platikanov, N. Rangelova, A. Scheludko, in: B.V. Derjaguin ŽEd.., Surface Forces and Liquid Boundary Layers, Nauka, Moscow, 1983, p. 200. Z.M. Zorin, Kolloidn. Zh. USSR 39 Ž1977. 1158. P.A. Kralchevsky, K. Dimitrov, I.B. Ivanov, Rep. Prog. Phys. 56 Ž1993. 1. F.P. Buff, J. Chem. Phys. 25 Ž1956. 146. G.W. Gibbs, The Collected Works, Longmans-Green, New York, 1928. A.I. Rusanov, Phase Equilibria and Surface Phenomena, Khimia, Leningrad, 1967. J.C. Melrose, Ind. Eng. Chem. 60 Ž1968. 53. F.C. Goodrich, in: E. Matijevic ŽEd.., Surface and Colloid Science, vol 1, Wiley, New York, 1969, p. 1. J.S. Rowlinson, B. Widom, Molecular Theory of Capillarity, Claredon Press, Oxford, 1982. V.S. Markin, M.M. Kozlov, S.L. Leikin, J. Chem. Soc. Faraday Trans. 84 Ž1988. 1149. I.B. Ivanov, P.A. Kralchevsky, Thin Liquid Films: Fundamentals and Applications, Marcel Dekker, New York, 1988, p. 49. V.M. Entov, R.M. Musin, ‘Preprint IPM RAN’, Moscow, No 560, 1996. V.L. Pokrovsky, A.L. Talapov, Zh. Eksp. Teor. Fiz. 78 Ž1980. 269. P. Bak, Rep. Prog. Phys. 45 Ž1982. 587. I. Frenkel’, T.A. Kontorova, Zh. Eksp. Teor. Fiz. 8 Ž1938. 1340. B.D. Josephson, Adv. Phys. 14 Ž1965. 419. I.O. Kulik, Zh. Eksp. Teor. Fiz. 51 Ž1966. 1952. A. Seeger, P. Schiller, in: W.P. Mason, R.N. Thurston ŽEds.., Physical Acoustics, vol IIIA, Academic, New York, 1966, p. 361. A. Barone, G. Paterno, Physics and Applications of the Josephson Effect, Wiley, 1982. A.J. Lichtenberg, M.A. Lieberman, Regular and Chaotic Dynamics, 2nd ed, Springer-Verlag, New York, 1992. K.G. Kornev, JETP 80 Ž1995. 1049. K. Kornev, V. Mourzenko, R. Dautov, Proceedings of the 9th International Conference on Computational Modeling and Computing in Physics, Sept. 16᎐21, 1996, D5, 11-97-112, Dubna, 1997, p 181. R. Dautov, K. Kornev, V. Mourzenko, Phys. Rev. E 56 Ž1997. 6929. C. Kittel, Introduction to Solid State Physics, Wiley, New York, 1956. P.G. de Gennes, Revue L’Institut Francais Petrole 47 Ž1992. 249. M.A. Lieberman, A.L. Lichtenberg, Phys. Rev A 5 Ž1972. 1852. G.R. Assar, R.W. Burley, in: N.P. Cheremisinoff ŽEd.., Encyclopedia of Fluid Mechanics, vol 3, Gulf Coleman, Houston, 1986, p. 26. W.R. Rossen, P.A. Gauglitz, AIChE J. 36 Ž1990. 1176. A.H. Falls, J.J. Musters, J. Ratulowski, SPE Res. Eng. 4 Ž1989. 55. J.E. Hanssen, M. Dalland, Foams in Porous Media-SUPRI -TR49, US DOE, 1986, p. 319. R.W. Flumerfelt, J. Prieditis, in: D.H. Smith ŽEd.., Surfactant-based mobility control: progress in miscible-flood enhanced oil recovery, ACS, Symposium Series 373, Washington, D.C, 1988, p. 295. 186 w124x w125x w126x w127x w128x w129x w130x w131x w132x w133x w134x w135x w136x w137x w138x w139x w140x w141x w142x w143x w144x w145x w146x w147x w148x w149x w150x w151x w152x w153x w154x w155x w156x w157x w158x w159x w160x w161x w162x w163x w164x w165x w166x w167x w168x w169x w170x w171x w172x K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 W.R. Rossen, J. Colloid Interf. Sci. 136 Ž1990. 1. W.R. Rossen, J. Colloid Interf. Sci. 136 Ž1990. 17. W.R. Rossen, J. Colloid Interf. Sci. 136 Ž1990. 38. W.R. Rossen, J. Colloid Interf. Sci. 139 Ž1990. 457. R.A. Albrecht, S.S. Marsden, Soc. Petr. Eng. J. 10 Ž1970. 51. R.J. Treinen, W.E. Brigham, L.M. Castanier, Apparent Viscosity Measurements of Surfactant Foam in Porous Media-Supri TR-48, US DOE, 1985. F.P. Bretherton, J. Fluid Mech. 10 Ž1961. 166. J. Ratulowski, H.C. Chang, Phys. Fluids A 1 Ž1989. 1642. K.G. Kornev, V.N. Kurdyumov, JETP 79 Ž1994. 252. H.L. Goldsmith, S.G. Mason, J.Colloid Sci. 18 Ž1963. 237. G. Hestroni, S. Haber, E. Wacholder, J. Fluid Mech. 41 Ž1970. 689. B.P. Ho, L.G. Leal, J. Fluid Mech. 71 Ž1975. 361. M.J. Martinez, K.S. Udell, J. Fluid Mech. 210 Ž1990. 565. C. Pozrikidis, J. Fluid Mech. 237 Ž1992. 627. L. Landau, B. Levich, Acta Physico-Chim. USSR 17 Ž1942. 42. B. Levich, Physicochemical Hydrodynamics, Prentice Hall, 1962. P.G. deGennes, Rev. Mod. Phys. 57 Ž1985. 827. L. Leger, J.F. Joanny, Rep. Prog. Phys. 55 Ž1992. 431. S.F. Kistler, in: J.C. Berg ŽEd.., Wettability, Surfactant Science Series 49, Marcel Dekker, New York, 1993, p. 311. E.O. Tuck, L.W. Schwartz, SIAM Rev. 32 Ž1990. 453. G.I. Barenblatt, Similarity, Self-Similarity and Intermediate Asymptotics: Theory and Applications to Geophysical Hydrodynamics, Gidrometeoizdat, Leningrad, 1978. G.I. Barenblatt, Ya.B. Zel’dovich, Ann. Rev. Fluid. Mech. 4 Ž1972. 285. F. Fairbrother, A.E. Stubbs, J.Chem. Soc. 1 Ž1935. 527. R.N. Marchessault, S.G. Mason, Ind. Chem. Eng. 52 Ž1960. 79. G.I. Taylor, J. Fluid Mech. 10 Ž1961. 161. L.W. Schwartz, H.M. Princen, A.D. Kiss, J. Fluid Mech. 172 Ž1986. 259. J.D. Chen, J. Colloid Interf. Sci. 109 Ž1986. 34. J. Ratulowski, H.C. Chang, J. Fluid Mech. 210 Ž1990. 303. C.W. Park, G.M. Homsy, J. Fluid Mech. 139 Ž1984. 291. E. Herbolzheimer, AIChE Annual Meeting, New York, 1987, Paper 68. H.C. Chang, J. Ratulowski, AIChE Annual Meeting, New York, 1987, Paper 681. K.L. Stebe, D. Barthes-Biesel, J. Fluid Mech. 286 Ž1995. 25. C.W. Park, Phys. Fluids A 4 Ž1992. 2335. J.B. Grotberg, Ann.Rev. Fluid Mech. 26 Ž1994. 529. D.P. Gaver III, D. Halpern, O.E. Jensen, I.B. Grotberg, J. Fluid Mech. 319 Ž1996. 25. J.W. Keuskamp, J. Lyklema, in: K.L. Mittal ŽEd.., Adsorption at Interfaces, ACS Symposium Series 8, American Chemical Society, Washington, DC, 1975, p. 191. A.D. Nikolov, P.A. Kralchevsky, I.B. Ivanov, D.T. Wasan, J. Colloid Interf. Sci. 133 Ž1989. 13. O. Belorgey, J.J. Benattar, Phys. Rev. Lett. 66 Ž1991. 313. S.E. Friberg, Langmuir 8 Ž1992. 1889. B. Cabane, R. Duplessix, J. Phys. 48 Ž1987. 651. S. Lionti-Addad, J.M. Meglio, Langmuir 8 Ž1992. 324. R. Bruinsma, J.M. Meglio, D. Quere, S. Cohen-Addad, Langmuir 8 Ž1992. 3161. J. Lucassen, Phys. Rep. 52 Ž1979. 217. D.A. Edwards, H. Brenner, D.T. Wasan, Interficial Transport Processes and Rheology, Buttenworth-Heinernann, Boston, 1991. J.B. Fournier, Phys. Rev. Lett. 75 Ž1995. 854. J.C. Earnshaw, D.J. Sharpe, J. Chem. Soc. Faraday Trans. 92 Ž1996. 611. O. Krichevsky, J. Stavans, Phys. Rev. Lett. 74 Ž1995. 2752. B.V. Derjaguin, L.D. Landau, Acta Physicochim. URSS 14 Ž1941. 633. K.G. Kornev, G.A. Shugai, Phys. Rev. E 58 Ž1998. 7606. K.G. Korne¨ et al. r Ad¨ . Colloid Interface Sci. 82 (1999) 127᎐187 w173x w174x w175x w176x w177x w178x w179x w180x w181x w182x w183x w184x w185x w186x w187x w188x w189x w190x w191x w192x w193x w194x w195x w196x w197x w198x w199x w200x w201x w202x w203x w204x w205x w206x w207x w208x w209x w210x w211x w212x w213x w214x w215x w216x 187 G.F. Teletzke, H.T. Davis, L.E. Scriven, Rev. Phys. Appl. 23 Ž1988. 989. E. Rabinowicz, Friction and Wear of Materials, Wiley, New York, 1965. F.P. Bowden, D.T. Tabor, Friction and Lubrication of Solids, Claredon Press, Oxford, 1986. L.D. Landau, E.M. Lifshitz, Mechanics, Nauka, Moscow, 1965. G.M. Zaslavsky, R.Z. Sagdeev, An Introduction to Nonlinear Physics: from Pendulum to Turbulence and Chaos, Nauka, Moscow, 1988. L.G. Aslamasov, A.I. Larkin, Zh. Eksp. Teor. Fiz. Pis’ma 9 Ž1969. 150. H.A. Barnes, J.F. Hutton, K. Walters, An Introduction to Rheology, 2nd ed, Elsevier, 1993. R. Musin, Ph.D. Thesis, Moscow, 1997. K. Lonngren, A. Scott, Solitons in Action, Academic Press, NY, 1978. R.K. Dodd, J.C. Eilbeck, J.D. Gibbon, H.C. Morris, Solitons and Nonlinear Wave Equations, Academic Press, London, 1984. J.S. Langer, in: G. Grinstein, G. Mazenko ŽEds.., Directions in Condensed Matter Physics, World Scientific, Singapore, 1986, p. 165. Y. Braiman, F. Family, H.G.E. Hentschel, Phys. Rev. E. 53 Ž1996. 3005. M. Doi, S.F. Edwards, The Theory of Polymer Dynamic, Claredon Press, Oxford, 1986. L.W. Holm, Soc. Petr Eng. J. 8 Ž1968. 359. P.G. deGennes, J. Chem. Phys. 55 Ž1971. 572. R.B. Bird, R.C. Armstrong, O. Hassager, C.F. Curtis, Dynamics of Polymeric Liquids, Wiley, New York, 1977. A.H. Cotrell, Dislocations and Plastic Flow in Crystals, Oxford University Press, London, 1953. W.T. Read, Dislocations in Crystals, McGraw Hill, New York, 1953. J.P. Hirth, J. Lothe, Theory of Dislocations, Wiley, New York, 1968. T. Suzuki, H. Yoshinaga, S. Takeuchi, Dynamics of Dislocations and Plasticity, Mir, Moscow, 1989. A.H. Falls, J.B. Lawson, G.J. Hirasaki, JPT 95 Ž1988.. D. Cohen, T.W. Patzek, C.J. Radke, J. Colloid Interf. Sci. 179 Ž1996. 357. R.N. Nabarro, Report of the Conference of the Strength of Solids, Phys. Soc. London, London, 1948, p. 75. C.J. Herring, J. Appl. Phys. 21 Ž1950. 5. I.M. Lifshitz, Zh. Eksp. Teor. Fiz. 44 Ž1963. 1349. I.M. Lifshitz, The Collected works of Il’ya Lifshitz, vol 1, Nauka, Moscow, 1987. T.F. Cooke, D.E. Hirt, in: R.K. Prud’homme, S.A. Khan ŽEds.., Foams, Marcel Dekker, NY, 1995, p. 339. C.C. Namboodri, M.W. Duke, Textile Res. J. 49 Ž1979. 156. R.C. Gregorian, Text. Chem. Color. 19 Ž1987. 13. D.E. Hirt, R.K. Prud’homme, L. Rebenfeld, AIChE J. 34 Ž1988. 326. D.E. Hirt, R.K. Prud’homme, L. Rebenfeld, Textile Res. J. 61 Ž1991. 47. E.L. Wright, Textile Res. J. 51 Ž1981. 251. N.C. Judd, W.W. Wright, SAMPE 14 Ž1978. 10. A.D. Mahale, R.K. Prud’homme, L. Rebenfeld, Polymer Eng. Sci. 32 Ž1992. 319. A.D. Mahale, R.K. Prud’homme, L. Rebenfeld 4 Ž1993. 199. B. Miller, I. Tyomkin, J. Colloid Interf. Sci. 162 Ž1994. 163. S.P. Gido, D.E. Hirt, S.M. Montgomery, R.K. Prud’homme, L. Rebenfeld, J. Dispersion Sci. Technol. 10 Ž1989. 785. B.J. Carrol, J. Colloid Interf. Sci. 57 Ž1976. 488. B.J. Carroll, Langmuir 2 Ž1986. 248. F. Brochard, J. Chem. Phys. 84 Ž1986. 4664. F. Brochard-Wyart, C.R. Acad. Sc. Ser. II 303 Ž1986. 1077. J.-M. d’Meglio, C.R. Acad. Sci. Paris Ser. II 303 Ž1996. 437. F. Brochard-Wyart, J.-M. diMeglio, D. Quere., C.R. Acad. Sci. Ser. II 304 Ž1987. 553. D. Quere, J.-M. diMeglio, F. Brochard-Wyart, Rev. Phys. Appl. 23 Ž1988. 1023.

© Copyright 2017 ExploreDoc