Effect of the adsorption component of the disjoining pressure on foam film drainage

Similar documents
Introduction to Interfacial Segregation. Xiaozhe Zhang 10/02/2015

Lecture 2 Thermodynamics and stability October 30, 2009

Harmonic oscillator approximation

Two-Layered Model of Blood Flow through Composite Stenosed Artery

Chapter 6 The Effect of the GPS Systematic Errors on Deformation Parameters

You must not circulate this book in any other binding or cover and you must impose this same condition on any acquirer.

Additional File 1 - Detailed explanation of the expression level CPD

Specification -- Assumptions of the Simple Classical Linear Regression Model (CLRM) 1. Introduction

Chapter 11. Supplemental Text Material. The method of steepest ascent can be derived as follows. Suppose that we have fit a firstorder

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Small signal analysis

The influence of Stern layer conductance on the. dielectrophoretic behaviour of latex nanospheres

Problem #1. Known: All required parameters. Schematic: Find: Depth of freezing as function of time. Strategy:

8 Waves in Uniform Magnetized Media

Not at Steady State! Yes! Only if reactions occur! Yes! Ideal Gas, change in temperature or pressure. Yes! Class 15. Is the following possible?

No! Yes! Only if reactions occur! Yes! Ideal Gas, change in temperature or pressure. Survey Results. Class 15. Is the following possible?

Team. Outline. Statistics and Art: Sampling, Response Error, Mixed Models, Missing Data, and Inference

Improvements on Waring s Problem

The multivariate Gaussian probability density function for random vector X (X 1,,X ) T. diagonal term of, denoted

Module 5. Cables and Arches. Version 2 CE IIT, Kharagpur

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Supplementary Notes for Chapter 9 Mixture Thermodynamics

CHAPTER X PHASE-CHANGE PROBLEMS

On the SO 2 Problem in Thermal Power Plants. 2.Two-steps chemical absorption modeling

Supporting Information. Hydroxyl Radical Production by H 2 O 2 -Mediated. Conditions

Improvements on Waring s Problem

A novel protocol for linearization of the Poisson-Boltzmann equation

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

Solution Thermodynamics

Assignment 4. Adsorption Isotherms

2.3 Least-Square regressions

CHAPTER 9 LINEAR MOMENTUM, IMPULSE AND COLLISIONS

Scattering of two identical particles in the center-of. of-mass frame. (b)

NUMERICAL DIFFERENTIATION

Gouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer.

Plastic Analysis and Design of Steel Plate Shear Walls

Appendix II Summary of Important Equations

MULTIPLE REGRESSION ANALYSIS For the Case of Two Regressors

Kinetic-Energy Density-Functional Theory on a Lattice

Weierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN

HYDRODYNAMIC LIMIT FOR A GAS WITH CHEMICAL REACTIONS

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

KEY POINTS FOR NUMERICAL SIMULATION OF INCLINATION OF BUILDINGS ON LIQUEFIABLE SOIL LAYERS

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

SOLUTION MANUAL ENGLISH UNIT PROBLEMS CHAPTER 9 SONNTAG BORGNAKKE VAN WYLEN. FUNDAMENTALS of. Thermodynamics. Sixth Edition

Two Approaches to Proving. Goldbach s Conjecture

Extended Prigogine Theorem: Method for Universal Characterization of Complex System Evolution

728. Mechanical and electrical elements in reduction of vibrations

A NUMERICAL MODELING OF MAGNETIC FIELD PERTURBATED BY THE PRESENCE OF SCHIP S HULL

Physics 120. Exam #1. April 15, 2011

MODELLING OF TRANSIENT HEAT TRANSPORT IN TWO-LAYERED CRYSTALLINE SOLID FILMS USING THE INTERVAL LATTICE BOLTZMANN METHOD

STAT 511 FINAL EXAM NAME Spring 2001

Electrical double layer: revisit based on boundary conditions

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Statistical Properties of the OLS Coefficient Estimators. 1. Introduction

AP Statistics Ch 3 Examining Relationships

1. The number of significant figures in the number is a. 4 b. 5 c. 6 d. 7

MATHEMATICAL AND COMPUTER HOMOGENIZATION MODELS FOR BULK MIXTURE COMPOSITE MATERIALS WITH IMPERFECT INTERFACES

NAME and Section No. it is found that 0.6 mol of O

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

4.2 Chemical Driving Force

Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

Wind - Induced Vibration Control of Long - Span Bridges by Multiple Tuned Mass Dampers

Electrical Circuits II (ECE233b)

Information Acquisition in Global Games of Regime Change (Online Appendix)

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

and Statistical Mechanics Material Properties

PHYS 100 Worked Examples Week 05: Newton s 2 nd Law

Supplementary information: Efficient mass transport by optical advection

Chemical Engineering Department University of Washington

Physics 111. CQ1: springs. con t. Aristocrat at a fixed angle. Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468.

Kernel Methods and SVMs Extension

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

ESCI 341 Atmospheric Thermodynamics Lesson 13 Phase Changes Dr. DeCaria

The Price of Anarchy in a Network Pricing Game

Entropy generation in a chemical reaction

Pythagorean triples. Leen Noordzij.

Phys 402: Raman Scattering. Spring Introduction: Brillouin and Raman spectroscopy. Raman scattering: how does it look like?

Name: SID: Discussion Session:

Physics 114 Exam 2 Fall 2014 Solutions. Name:

Quick Visit to Bernoulli Land

Method Of Fundamental Solutions For Modeling Electromagnetic Wave Scattering Problems

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Chapter 13: Multiple Regression

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Chapter 8. Potential Energy and Conservation of Energy

Numerical Methods for Solving Turbulent Flows by Using Parallel Technologies

A Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients

I wish to publish my paper on The International Journal of Thermophysics. A Practical Method to Calculate Partial Properties from Equation of State

Formal solvers of the RT equation

3) Thermodynamic equation to characterize interfaces

Foundations of Arithmetic

Supporting Information

Physics 240: Worksheet 30 Name:

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

Canonical transformations

DUE: WEDS FEB 21ST 2018

This appendix presents the derivations and proofs omitted from the main text.

Transcription:

Collod J. 75 (2013) 176-180 [arxv 1209.4745] Effect of the adorpton component of the djonng preure on foam flm dranage Stoyan I. Karakahev 1, Anh V. Nguyen 2 and Roumen Tekov 1 1 Department of Phycal Chemtry, Unverty of Sofa, 1164 Sofa, Bulgara 2 Department of Chemcal Engneerng, Unverty of Queenland, Brbane, QLD 4072, Autrala The preent work tryng to explan a dcrepancy between expermental obervaton of the dranage of foam flm from aqueou oluton of odum dodecyl ulfate (SDS) and the theoretcal DLVO-accomplhed Reynold model. It hown that, due to overlap of the flm adorpton layer, an adorpton component of the djonng preure mportant. The pre-exponental factor of the adorpton component wa obtaned by fttng the expermental dranage curve. It correpond to a lght repulon, whch reduce not only the thnnng velocty a oberved expermentally but correct alo the flm equlbrum thckne. Karakahev et al. 1 have nvetgated the dranage of foam flm of dlute aqueou oluton of odum dodecyl ulfate wthn the concentraton range 1-100 µm. They tred to decrbe the knetc of foam flm thnnng by the Reynold lubrcaton approxmaton accountng for the Marangon effect, urface hear vcoty and DLVO force. Sgnfcant dcrepancy between the theoretcal predcton and the expermental reult wa oberved. The detaled analy howed that the devaton of the theory from the expermental data orgnate from the nteracton between the flm urface. Therefore, t wa concluded that the clacal DLVO theory only not uffcent to match the expermental data. It wa uggeted that the dcrepancy between theory and experment due to a neglected varaton of the adorpton component of the urface tenon durng the flm dranage. Large number of lterature confrm the applcablty of the DLVO theory to foam flm. However, number of paper 2-6 report devaton of th theory from expermental data. Th dcrepancy pronounced motly n thn flm between hydrophobc urface. To olve the problem ome author 7-11 ntroduced n the theory an addtonal non-dlvo force, the o-called hydrophobc force, whch can be attractve or repulve. 12-14 There number of attempt n the lterature 15-19 to explan the nature of the hydrophobc nteracton but tll no full agreement of the opnon reached. The clacal DLVO theory doe not account alo for other nteracton n the thn lqud flm. For ntance, the nteracton between the overlappng dffuve adorpton layer hould contrbute to the overall nteracton between the flm urface and th contrbuton hould ncreae wth decreang flm thckne. The dea of the adorpton nteracton between the flm urface orgnate from the work of Ah, Everett and Radke 20 and t further developed by the Ruan chool of collod chemtry. The dperon nteracton between the olute and the flm urface accounted for 21-23 and t reult n a correcton n the van der Wall component of the djonng preure. Th

addtonal adorpton term n the total nteracton between the flm urface could be mportant but t ha been overlooked n a large volume of lterature caung dverty of the opnon regardng the hydrophobc nteracton. The reaon for th that the reearcher cted above have decrbed the urfactant dtrbuton only a a reult of nteracton wth the urface but neglected the nteracton between olute molecule. Of coure, the latter are not mportant n dlute oluton far away from the urface, but when the adorpton condered the concentraton near a urface tremendouly ncreaed. Tekov and Schulze 17 uggeted frt a clear thermodynamc nterpretaton of the adorpton term n the total djonng preure. They called t hydrophobc force, nce the orgn of the adorpton the urface hydrophobcty and the urfactant ablty to reduce t. The am of th paper to employ th approach for explanaton of our expermental data. 24 The good agreement wll certanly draw attenton on the mportance of the adorpton djonng preure. Accordng to the thn lqud flm thermodynamc any change of the flm free energy F at contant temperature gven by df pdv da dn (1) where the extenve flm parameter are volume V, flm area A and number of mole { n } of the flm component. The relatonhp between the ntenve parameter preure p, flm tenon and electrochemcal potental { } gven by the Gbb-Duhem equaton Vdp Ad nd 0 (2) It known that the thn lqud flm are anotropc tructure 25 and ther preure tenor poee two dtnct component, the normal and tangental one. At equlbrum the normal component of the preure tenor equal to the ga preure outde, whle the tangental component equal to the preure n the mencu adjacent to the flm. The preure p the normal component of the preure tenor. The flm tenon cont n two addtve, 26 where h V / A the flm thckne, 2 h (3) The purely nterfacal part twce the flm urface tenon whle the bulk part accounted by the djonng preure. Introducng Eq. (3) n Eq. (1) the latter change to df padh ( p ) hda 2dA dn (4)

It obvou now that the normal and tangental component of the flm preure tenor are not equal and the djonng preure ther dfference. Ung Eq. (3) one can derve an alternatve form of Eq. (2) Vdp Ad(2 h) nd 0 (5) After Gbb the flm can be dealzed by fllng t wth the bulk lqud from the mencu. Hence, ubtractng from Eq. (5) the Gbb-Duhem relaton dpl cd for the lqud n the mencu, where { c } are the concentraton of the chemcal component there, and keepng n mnd that p p L, one yeld an mportant nterfacal Gbb-Duhem relaton 27 d d dh /2 (6) where { ( n cv)/ 2 A} are the component adorpton. Eq. (6) provde traghtforward an mportant defnton of the djonng preure a the thckne dervatve of the flm urface tenon 2( ) h (7) a well a the followng Maxwell relaton for the djonng preure ( ) h 2( ) h (8) The latter already hnt the mportant effect of adorpton on the djonng preure. 28,29 Snce the urfactant could be charged pece the flm urface tenon depend on electrotatc a well. It can be plt nto uperpoton of water, electrotatc and adorpton component, W EL, whch are ndependent f the urface potental doe not depend on the flm thckne. Thu, durng the flm dranage the adorpton and urface charge denty, repectvely, can vary but the electrotatc component wll not be affected by. At contant temperature the water component depend only on the flm thckne, whle the urfactant component depend on the adorpton. Subttutng th preentaton n Eq. (7) the djonng preure plt alo nto three dtnct component EL

VW EL 2 ( )( ) h (9) where VW 2( W / h) and EL 2( EL / h) are the well-known van der Waal and electrotatc component. Indeed, at low urface potental the electrotatc component of the urface tenon equal to, where the recprocal Debye length, 2 EL 0 tanh( h/ 2)/ 2 and the correpondng electrotatc djonng preure acqure 2 2 2 EL 0 / 2coh ( h/ 2) t clacal form. 30 Let u conder now the lat adorpton component of the djonng preure n Eq. (9). To calculate t thckne dependence of adorpton one can employ the Maxwell relaton (8). Introducng the followng defnton hx X( h) X( ) for a dfference between the value of a property X of the equlbrum flm wth thckne h and nfnty, repectvely, one can wrte p. Note, that changng the flm thckne only t tangental preure component h L change, whle the normal one p p ( h ) reman contant. Thu, the Maxwell relaton (8) can be conecutvely modfed to G L pl 2( ) ( ) h h( ) h hc h (10) Knowng the adorpton otherm c ( ) at contant urface potental one able to ntegrate th equaton to obtan the thckne dependence of adorpton. If the change of the concentraton and adorpton, repectvely, are mall one can employ the followng lnear relatonhp h a hc, where a ( / c ) h the adorpton length on a ngle flat lqud/ga nterface. The latter, repreentng the thckne of the adorpton layer, depend on the adorpton equlbrum contant and. Solvng now the lnearzed dfferental equaton (10) yeld h h 0 exp( ) (11) 2a where { 0 } the dfference between the adorpton n a urfactant blayer ( h 0 ) and on a ngle flat urface ( h ). Subttutng now th expreon nto the defnton of the adorpton djonng preure from Eq. (9) lead to h 0 ( ) exp( ) (12) a 2a

Note that dependng on the gn of { 0 }, the adorpton djonng preure can be ether potve or negatve. It could be alo zero f no change n the adorpton n a blayer and on a ngle flat urface take place. Th probably the mot wdepread cae, whch explan why the adorpton djonng preure tll not well tuded. The preent thermodynamc theory cannot gve any value of { 0 } but jut aumng them decrbe the thckne dependence of ther effect. The dranage of thn lqud flm depend ubtantally on the moblty of flm urface. 25 Our preent etmate how, however, that wthn the pecfed SDS concentraton range the Marangon effect alway trong enough to block the tangental flow on the flm nterface. Hence, the dranage velocty can be well approxmated by the clacal Stefan-Reynold equaton dh dt 3 2 h ( p ) 2 3R (13) where p the capllary preure, R the flm radu, and the lqud vcoty. The djonng preure n Eq. (13) of crucal mportance for the modelng of the dranage. How t wa hown above, a uperpoton on the van der Waal, electrotatc and adorpton component. To determne the effect of the adorpton djonng preure correctly relable expreon for the DLVO component are requred. The van der Waal djonng preure between the flm urface can be etmated from the expreon 3 /6. Snce VW A h the flm thckne h alway larger than 150 nm and the Hamaker contant about 21 A 2 10 J, the van der Waal djonng preure neglgble for the preent ytem. At contant urface potental, the electrotatc djonng preure, calculated by the exact numercal oluton of the non-lnear Poon-Boltzmann equaton, em-analytcally decrbed a 31 2 1 2 32 tanh ( / 4)[(1 coh ) nh ( / 4)exp( )] (14) EL RTc y h f y f h where f 2coh(0.332 y 0.779) for y 7, y F / RT a dmenonle urface potental on a ngle flat ar/oluton nterface. In the cae of a ngle urfactant Eq. (12) reduce to 0 exp( h/ 2 a)/ a (15) where the dfference of the adorpton component of the urface tenon on the blayer and on a ngle flat nterface. Here a the urfactant adorpton length. Snce could be ether potve or negatve dependng on the nteracton between the two monolayer

Flm thckne [nm] Flm thckne [nm] of the blayer, could be alo repulve or attractve, repectvely. In order to compare the above theory wth the expermental data, Eq. (13) wa numercally ntegrated ung a fourth-tep Runge-Kutta method. A macro wa wrtten n the VBA (Vual Bac for Applcaton) programmng language avalable from Mcrooft Excel. The meaured value for the zeta potental were adopted for the urface potental n EL.The adorpton length a are determned by the model of Kralchevky et al. 32 The computed flm thckne v. tme wa compared wth the expermental one for each of the SDS concentraton n order to obtan the bet ft of the free parameter (ee Fg. 1). 380 340 1 M SDS 300 260 220 180 140 100 0 1 2 3 4 5 6 7 Tme [] 310 10 M SDS 270 230 190 150 110 0 2 4 6 8 10 Tme []

Equlbrum flm thckne [nm] Flm thckne [nm] 300 260 100 M SDS 220 180 140 100 0 2 4 6 8 Tme [] 200 160 120 80 1 10 100 1000 Concentraton [M] Fgure 1. Example of expermental data and theoretcal ft of the flm thckne v. tme for ome SDS concentraton. The lat plot preent the equlbrum flm thckne v. SDS concentraton. Thee data along wth the capllary preure are preented n Table 1. The juxtapoton between the expermental and theoretcal flm thckne v. tme preented n the complementary materal attached to the paper. Table 1 how that the adorpton djonng preure potve nce 0 0, whch correpond to addtonal (non-dlvo) repulon between the flm urface. Th repulon orgnate from the overlap of the adorpton layer of the two flm urface. The latter ndcate that the adorpton component of the flm urface tenon ncreae wth the decreae of the flm thckne. Hence, the urfactant adorpton for th partcular

cae (SDS) dmnhe durng the flm thnnng. In general, the adorpton djonng preure hould dappear at zero SDS concentraton. A expected ncreae wth ncreang urfactant concentraton. Snce the adorpton length a reduce wth ncreang of c, the adorpton djonng preure become horter ranged and tronger. Table 1. Capllary preure p, zeta potental, adorpton length a and fttng parameter v. concentraton of SDS c [µm] p [Pa] [mv] a [nm] [µn/m] 1 72.6-63.0 275 7 5 72.5-56.1 275 7 8 72.5-52.7 275 7 10 72.4-49.2 275 7 50 72.4-48.6 274 7 100 71.8-48.0 274 10 500 70.0-52.8 249 13 200 EL Djonng preure (Pa) 100 0 p VW -100 0 100 200 300 400 Flm thckne (nm) Fgure 2. Electrotatc EL, van der Waal VW, adorpton and total djonng preure v. the flm thckne for 10 µm SDS oluton. A general queton here how large the contrbuton of the adorpton djonng preure nto the total djonng preure. Another queton f th non-dlvo force long or hort

ranged. One can fnd anwer of thee queton n Fg. 2, whch how the electrotatc van der Waal EL, VW, adorpton and total djonng preure v. the flm thckne. The traght dahed lne n the fgure repreent the uckng capllary preure p. A een the adorpton nteracton between the flm urface decayng upon the flm thckne much weaker than the electrotatc and van der Waal nteracton. The contrbuton of the adorpton nteracton to the total nteracton between the flm urface gnfcantly mall. However, above a gven flm thckne (ca. 300 nm for the partcular cae of 10 M SDS) the adorpton nteracton preval over the electrotatc one. Depte beng long-ranged, the adorpton djonng preure cannot become equal to the capllary preure at any flm thckne due to t mall value. Th mean that n abence of electrotatc djonng preure an equlbrum flm cannot be formed and the flm wll thn untl rupture. The abence of adorpton nteracton, however, wll reflect n a gnfcantly maller equlbrum thckne of the flm (95 nm ntead of 160 nm) and fater flm dranage. The preent paper prove the extence of adorpton non-dlvo djonng preure between the foam flm urface. It orgnate by the overlap between the adorpton layer and can be attractve, repulve or vanhng. The adorpton djonng preure related to the properte of the adorpton layer. It part of the hydrophobc nteracton between the flm urface. 17 If the urfactant adorpton dmnhe upon the decreae of the flm thckne the adorpton nteracton repulve and vce vera. We menton here a well that uch flm develop treamng potental upon ther dranage. 33 Th theory valdated by experment on knetc of thnnng of foam flm from SDS wthn the concentraton range 1-500 µm. Ft upon the parameter for each one of the concentraton performed. A expected ncreae wth ncreae on the urfactant concentraton. Thu defned the adorpton nteracton doe not dffer from th one defned by Tekov and Schulze 17 and Wang and Yoon. 12,13 A more detaled tudy for the effect of the adorpton otherm on the adorpton component of the djonng preure can be found n Ref. 34 1. Karakahev, S.I.; Manev, E.D.; Nguyen, A.V., Collod Surf. A 2008, 319, 34 2. Exerowa, D.; Kolarov, T.; Khrtov, K., Collod Surf. 1987, 22, 171 3. Buchavzov, N.; Stubenrauch, C., Langmur 2007, 23, 5315 4. Mhra, N.C.; Muruganathan, R.M.; Müller, H.-J.; Krutev, R., Collod Surf. A 2005, 256, 77 5. Bowen, W.R.; Doneva, T.A.; Stoton, J.A.G., Collod Surf. A 2002, 201, 73 6. Krutev, R.; Müller, H. J., Langmur 1999, 15, 2134 7. Iraelachvl, J.N.; Pahley, R.M.; Perez, E.; Tandon, R.K., Collod Surf. 1981, 2, 287 8. Iraelachvl, J.; Pahley, R., Nature (London) 1982, 300, 341 9. Rabnovch, Y.I.; Derjagun, B.V., Collod Surf. 1988, 30, 243 10. Fa, K.; Nguyen, A.V.; Mller, J.D., J. Phy. Chem. B 2005, 109, 13112

11. Erkon, J.C.; Ljunggren, S.; Claeon, P.M., J. Chem. Soc., Faraday Tran. 2 1989, 85, 163 12. Wang, L.; Yoon, R.-H., Langmur 2004, 20, 11457 13. Wang, L.; Yoon, R.-H., Collod Surf. A 2005, 263, 267 14. Karakahev, S.I.; Phan, C.; Nguyen, A.V., J. Collod Interface Sc. 2005, 291, 489 15. Schalchl, A.; Sentenac, D.; Benattar, J.J.; Bergeron, V., J. Chem. Soc., Faraday Tran. 1996, 92, 2317 16. Evan, D.R.; Crag, V.S.J.; Senden, T.J., Phyca A 2004, 339, 101 17. Tekov, R.; Schulze, H.J., Langmur 1997, 13, 5674 18. Yoon, R.-H.; Akoy, B.S., J. Collod Interface Sc. 1999, 211, 1 19. Bau, S.; Nandakumar, K.; Malyah, J.H., J. Collod Interface Sc. 1996, 182, 82 20. Ah, S.G.; Everett, D.H.; Radke, C., J. Chem. Soc., Faraday Tran. 1973, 69, 1256 21. Derjagun, B.V.; Churaev, N.V., Kollodn. Zh. 1975, 37, 1075 22. Derjagun, B.V., Collod Polym. Sc. 1980, 258, 433 23. Exerowa, D.; Churaev, N.V.; Kolarov, T.; Epova, N.E.; Panchev, N.; Zorn, Z.M., Adv. Collod Interface Sc. 2003, 104, 1 24. Karakahev S.I.; Nguyen A.V., Collod Surf. A 2007, 293, 229 25. Ivanov, I.B. (Ed.), Thn Lqud Flm, Marcel Dekker, New York, 1988 26. Ruanov, A.I., Phae Equlbra and Surface Phenomena, Khmya, Lenngrad, 1967 27. Tohev, B.V.; Ivanov, I.B., Collod Polym. Sc. 1975, 253, 558 28. Derjagun, B.V.; Starov, V.M.; Churaev, N.V., Collod J. (USSR) 1976, 38, 411 29. Valeff, C.S.; Tohev, E.T.; Ivanov, I.B., n Surface Force and Boundary Lqud Layer, B.V. Derjagun (Ed.), Nauka, Mocow, 1983, p. 168 30. Verwey, E.J.W.; Overbeek, J.T.G., Theory of the Stablty of Lyophobc Collod, Elever, Amterdam, 1948; p. 218 31. Nguyen, A.V.; Evan, G.M.; Jameon, G.J., J. Collod Interface Sc. 2000, 230, 205 32. Kralchevky, P.A.; Danov, K.D.; Broze, G.; Mehreteab, A., Langmur 1999, 15, 2351 33. Karakahev, S.I.; Tekov, R., Langmur 2011, 27, 265 34. Tekov, R., Ann. Unv. Sofa, Fac. Chem. 2011, 102/103, 273