<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Mathematical modeling of separation of watered oil-containing mixture</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>V.A. Zelenskiy</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>A.A. Sushin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>A.I. Shchodro</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Samara National Research University</institution>
          ,
          <addr-line>34 Moskovskoe Shosse, 443086, Samara</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2017</year>
      </pub-date>
      <fpage>16</fpage>
      <lpage>19</lpage>
      <abstract>
        <p>A mathematical description of the separation process of the watered oil-containing mixture has been proposed. The mathematical model is based on the determination of the speed of movement of oil globules in a constrained flow. The formula for determining the velocity of the globules is derived from theoretical provisions and experimental data. The key issue of simulation is determining the separation time. This characteristic is calculated as a time period from the moment of arrival of the portion of oil-containing mixture to the first chamber of the separator to the moment of formation of the continuous phase layer of the given thickness on the surface. The mathematical model enables us to determine the separation time taking into account geometric characteristics of the separation device at the high level of water content. Thus, a topical problem of improving the performance of the separation device without sacrificing the quality of the commercial oil is solved.</p>
      </abstract>
      <kwd-group>
        <kwd>separation device</kwd>
        <kwd>watered oil-containing mixture</kwd>
        <kwd>oil globules</kwd>
        <kwd>constrained flow</kwd>
        <kwd>separation time</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        Computer Modeling / V.A. Zelenskiy, A.A. Sushin, A.I. Shchodro
indirect measuring of density the mixture watercut value is obtained [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. This information is used to determine the separation
time [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. Three approaches to obtain this parameter can be highlighted.
      </p>
      <p>1) Separation time is a period of time from the moment the first chamber of the separator is filled (up to a given level, not
exceeding the overflow level) until full decomposition of the oil-containing mixture. This approach did not find practical
use as it does not take into account the real velocity of demulsification process.
2) Separation time is a period of time from the moment of arrival of the portion of oil-containing mixture to the first
chamber of oil and gas separator to the moment of its decomposition to the state determined by GOST P 51858-2002. This
approach is more progressive than the first one, but does not take into account the continuous nature of OTCS operation.
3) Separation time is a period of time from the moment of arrival of the proportion of oil-containing mixture to the first
chamber of the separator to the moment of formation of the continuous phase layer of the given thickness on the surface.</p>
      <p>This is the preferred method to be used in the framework of this research.</p>
    </sec>
    <sec id="sec-2">
      <title>3. Mathematical description of separation process</title>
      <p>
        There is no strict mathematical description of the stratification of the oil-containing mixture as Navier-Stokes equation is valid
just for laminar flow and some special cases [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. However, there are a large number of empirical and semi-empirical
relationships that determine the nature of the processes in the stratification of emulsions that can be used as a basis [
        <xref ref-type="bibr" rid="ref6 ref7">6, 7</xref>
        ]. Under
the conditions of strong watercut of oil mixture (70% ... 90%), oil is a dispersed phase and water is a dispersed medium. The oil
particle rising to the surface experiences the difference between the gravity force and the lifting Archimedes' force [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ].
F  g d 3p ,
      </p>
      <p>6
where Δρ is the difference between dispersed phase and dispersed medium particles' density, g is the gravitational acceleration,
d is the particle's diameter. Resistance force of the continuous medium:</p>
      <p>Fc   о d42 2о  с ,
where ξ0 is the coefficient of hydraulic resistance of the continuous medium to the movement of a single particle in it, ω0 is the
velocity of a single particle relative to the medium, ρc is the density of the continuous medium. Let us assume that the
temperature at all points of a separation device chamber is the same. Then there is no thermal convection. If the particle's
velocity in the medium is constant: ΔF = Fc the Reynolds criterion is determined by the following ratio:
Reo  0d c ,
 c
where μc is the dynamic viscosity of the continuous medium. Archimedes' criterion is as follows:</p>
      <p>Ar 
d 3 g  д   с
 c2  с
,
where vc is the kinematic viscosity of continuous medium, ρd is the density of the dispersed phase. Taking these criteria into
account, we can derive the following equality:
 од   о (1  ) n ,</p>
      <p>Reд  (1  ) Reо .
where ξod is the coefficient of hydraulic resistance for a continuous medium for a single particle under the constrained flow
conditions, φ is the volume fraction of the dispersed phase within the system. It would be useful to define the type of function f
(φ) for the small and large Reynolds criterion values. Experimental studies showed that the velocities of particles depositing are
connected by the following relation:
where ω0d is the depositing rate of the particle relative to the continuous medium in constrained flow conditions, ω0d is the rate
of free depositing of the particle, and n index is to be determined. Using the obtained parameter called the volume fraction of the
dispersed phase, we have:</p>
      <p>Computer Modeling / V.A. Zelenskiy, A.A. Sushin, A.I. Shchodro</p>
      <p>
        It is experimentally shown [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] that with small values of Reynolds criterion (Re &lt; 500) the hydraulic resistance coefficient of
the medium equals:
 d 
0.843 lg( / 0.065) Re
where θ is the coefficient of the particle shape, equal to the ratio of the surface area of the spherical particle to the surface area of
the real particle of the same volume. With small values of Reynolds coefficient, we can assume that:
.
      </p>
      <p>For large values of Reynolds coefficient Re the following expression is valid:
f ( )  (1  )n
f ( )  (1  )2n .</p>
      <p>
        It is experimentally shown [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] that the function f (φ), both in case of large and small values of Re varies from (1-φ)-4.65 to
(1φ)-4.78. Then we can take an average index and write:
      </p>
      <p>f ( )  (1  )4.72 .</p>
      <p>Taking into account the expressions derived, it is obtained that the ratio of the particle depositing rate relative to the
continuous medium under the constrained flow conditions to the particle's free depositing rate is equal to:
оd /0  (1  )4.72 .</p>
      <p>
        There exist empirical formulas which enable us to account for the influence of constrainity [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ]. For example, when φ &lt; 0.3,
the following formula is applied:
      </p>
      <p>оd /o  (1  )2101.82 .</p>
      <p>For φ&gt;0,3 the following formula is applied:
 оd / o  0.123 (1  )3 .</p>
      <p></p>
      <p>Table 1 shows comparative data of the calculation of the velocity of the oil globules rising for the known relations and
calculated according to the obtained formula. Numerical values are reduced to the velocity of oil globules freely rising to t he
surface, i.e. relative ones.</p>
      <p>
        Data for the watercut degree from 70% to 90% have not been found. The velocities of oil globules rising distributed by
fractions are of practical interest. The known relations are obtained for a case when the dispersed phase is water. In our case, the
dispersed phase is oil. In paper [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] it is assumed that the distribution of water drops in oil after filling the chamber of separ ator is
uniform. Therefore, the watercut B in any vertical section is the same. The relative velocity of the constrained surfacing of oil
globules of di diameter in this case equals:
where dmax is the maximum size of the globule. Generally, it is suggested to define the separation time through the velocity of
the oil globules rising and the geometric parameters of the oil and gas separator. The velocity of constrained rising of oil
globules is as follows:
(w0d / w0 )i   1  B 
1  B 1  (di / dmax )2 
4.72
18 c
      </p>
      <p>Computer Modeling / V.A. Zelenskiy, A.A. Sushin, A.I. Shchodro
where pd, pc is the density of dispersed and continuous medium, kg/m3; μс is the viscosity of the continuous medium, Pa s; d is
the diameter of the globule; g is the gravitational acceleration, m/s2. While rising, globules of different sizes are moving at
different speeds. It is proposed to use the following expression to describe the calculation of the constrained rising of globules:
4.72
(w0d / w0 )i  ( pd  pc )d 2 g  1  B  ,</p>
      <p>18 c 1  B 1  (di / dmax )2 </p>
      <p>The equation obtained enables us to calculate the spectrum of velocities of constrained rising of the oil globules, taking into
account the change of watercut of emulsion according to the height of the partition wall between the chambers of the separation
device. In Table 2 there are given the data for the watercut values of the studied range from 70% to 90%. In accordance with
table 2 and geometric characteristics of the separation device, it is possible to determine the separation time for oil-containing
mixture, taking into account the composition of the mixture and watercut value.
100
150
200</p>
      <p>The velocity of the oil globules rising, cm/s</p>
    </sec>
    <sec id="sec-3">
      <title>4. Conclusion</title>
    </sec>
    <sec id="sec-4">
      <title>Acknowledgements References</title>
      <p>The obtained relations enable us to mathematically describe the separation process through the velocity of the globules in the
case of a high watercut value for oil-containing mixture, which demonstrates the scientific novelty. Using the mathematical
model data and knowing the geometric characteristics of the separation device, separation time can be calculated. The precise
definition of separation time improves the performance of the separation device without sacrificing the quality of the
commercial oil, which is of great practical importance.</p>
      <p>The authors express their gratitude for LLC Coordination (Ufa), LLC "GIRS", "Neftestroy" Educational Centre (Samara) for
the fruitful cooperation and experimental data provided.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Zelenskiy</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shchodro</surname>
            <given-names>AI</given-names>
          </string-name>
          .
          <article-title>Increasing the efficiency of separation by controlling the differential density of the oil and gas mixture</article-title>
          .
          <source>Bulletin of Samara State Technical University. Series "Technical Sciences"</source>
          <year>2015</year>
          ;
          <volume>1</volume>
          (
          <issue>45</issue>
          ):
          <fpage>178</fpage>
          -
          <lpage>183</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Zelenskiy</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shchodro</surname>
            <given-names>AI</given-names>
          </string-name>
          .
          <article-title>The automated control system of the oil and gas separator providing control over the density of the oil-containing mixture</article-title>
          .
          <source>Bulletin of Samara State Technical University. Series "Technical Sciences"</source>
          <year>2016</year>
          ;
          <volume>1</volume>
          (
          <issue>49</issue>
          ):
          <fpage>15</fpage>
          -
          <lpage>23</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>Zelenskiy</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shchodro</surname>
            <given-names>AI</given-names>
          </string-name>
          .
          <article-title>Analysis of errors of measuring the density of the oil-containing mixture and their impact on the determination of the separation time</article-title>
          .
          <source>Bulletin of Samara Scientific Centre of the Russian Academy Sciences</source>
          <year>2016</year>
          ;
          <volume>18</volume>
          (
          <issue>3</issue>
          ):
          <fpage>896</fpage>
          -
          <lpage>901</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Zelenskiy</surname>
            <given-names>VA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shchodro</surname>
            <given-names>AI</given-names>
          </string-name>
          .
          <article-title>Method, mathematical model, and the algorithm of control over the oil separation process</article-title>
          .
          <source>Bulletin of Samara State Technical University. Series "Technical Sciences"</source>
          <year>2016</year>
          ;
          <volume>3</volume>
          (
          <issue>51</issue>
          ):
          <fpage>21</fpage>
          -
          <lpage>28</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Cochin</surname>
            <given-names>NE</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kibel</surname>
            <given-names>IA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rose</surname>
            <given-names>NV</given-names>
          </string-name>
          .
          <article-title>Theoretical hydromechanics, part 2</article-title>
          . Moscow: Fizmatlit,
          <year>1963</year>
          ; 727 p.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Ponomarev</surname>
            <given-names>SV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mishchenko</surname>
            <given-names>SV</given-names>
          </string-name>
          .
          <article-title>Methods and devices for measuring the effective thermal performance of the flows of technological liquids</article-title>
          .
          <source>Tambov: Publ. Tambov State Technical University</source>
          ,
          <year>1997</year>
          ; 249 p.
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Astarita</surname>
            <given-names>J</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Maruchchi</surname>
            <given-names>J.</given-names>
          </string-name>
          <article-title>The basics of hydromechanics of non-Newtonian fluids</article-title>
          .
          <source>Мoscow: "Mir"</source>
          ,
          <year>1978</year>
          ; 312 p.
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>