Page 128 - Percolation Models for Transport in Porous Media With
P. 128

6.3  MERCURY ELECTRIC POROMETRY                                      121


            It easily follows  both from  the physical  meaning of the  percolation  model of
         transfer phenomena in a heterogeneous medium and from  the immediate analysis
         of the relation  (6.28)  that the latter expression, as well  as  (6.30),  is valid within
         the interval 0 < Ti  < Tc·
            In measuring a  11  ( r  i) some intervals r  i < r < r  i +  ~i  can be obtained where elec-
         tric  conductivity  of  the  medium  does  not  change,  i.e.,   da 11 (ri)fdri
         =  0.  If ~i do  not  contain  rc,  then  in  these  intervals  A(ri)  =  <f>(ri)  =  0,  and
         from  (6.30),  f(ri)  = 0  when  ri  < r  < ri + ~i·  This  means  that  there  are  no
         capillaries with  Ti  E ~i in  the medium.  If, however,  some  ~i contains an rc  to-
         gether with its neighborhood, then after substituting A(ri) =  <f>(ri)  =  0 in  (6.30),
         we also obtain  f(ri) =  0.  In this case,  the fact  that the integral of f(r)  is  in  the
         denominator of (6.29)  formally  causes  a  zero-over-zero indeterminacy.  However
         the equality da 11 fdri  = f(ri)  = 0  near  the initially introduced  rc  simply  means
         that there are no capillaries with  Ti  close  to rc.  Therefore it suffices  to decrease
         rc  up  to the closest  ri's, for  which  da 11 /dri  f. 0,  and take it as the new  Tc·  This
         changes neither the meaning nor the contents of all formulas, but allows to get rid
         of the formal indeterminacy.  These cases can be encountered only in those media
         which  possess an  f(r)  function  with  two  global  maximums,  i.e.,  in  those having
         two different types of porosities (e.g., porous and capillary or block and inter  block)
         described by the same f(r).
            Taking account of the new definition of rc, estimate A(ri) and <f>(ri).  It follows
         from  (6.27)  and  (6.29), noting that 0 ~ f(r)  ~  Mo  < oo,  that

                                                                            (6.31)




            Here Nc  is the average value of f(r) in the first segment r1  ~  r  ~  rc  (rc- r1  =
         t5').
            Since the contraction mappings principle is valid for the Volterra integral equa-
         tion of the second type for  every finite  A(ri)  and <f>(ri)  [77,  78],  the relationships
         (6.31)  speak in  favor  of the existence  and  the  uniqueness  of the solution  to the
         equation  (6.30)  and justify  the  use  of the  method  of successive  approximations
         with an arbitrary initial function  j( 0 )(r) for  solving this equation.
            In  this case,  once  we  have a  normalized function  j(n)(r)  as  the n-th approx-
         imation,  to obtain the next approximation according to (6.30),  we  first  calculate
         its non-normalized value
                          j<n+l)(ri) = A(n)(ri) J                          (6.32)
                                            rc
                                              f(n)(r) ~~ + <f>(n)(ri)

                                            r;
   123   124   125   126   127   128   129   130   131   132   133