<!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>Analysis of dependence of the 2FSI subsystem parameters on the labyrinth seal diameter of the centrifugal compressor in the gas pumping unit</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>L N Butymova</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>V Y Modorskii</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Perm National Research Polytechnic University (PNRPU)</institution>
          ,
          <addr-line>Komsomolsky Av. 29, Perm, Russia, 614990</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2018</year>
      </pub-date>
      <fpage>409</fpage>
      <lpage>414</lpage>
      <abstract>
        <p>The system of centrifugal compressor GTU has several 2FSI subsystems: impeller, labyrinth seals (LS), dummies, etc. The article analyzes dependence of one subsystem - the labyrinth seal - on its diameter. The effect of this parameter on the wave processes occurring in labyrinth seals during vibrations of the rotor is studied. A solution is obtained for a nonstationary aeroelastic three-dimensional gas flow problem in deformable LS. The amplitude-frequency characteristics of the gas dynamic and vibration parameters of the LS CC of GTU are obtained. Time dependences of the gas-dynamical force and the work of gasdynamic forces in the LS zone at different rotor diameters are obtained.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Physical statement</title>
      <p>
        The computational model of the labyrinth seal was considered, which is a set of the LS gas-dynamic
gap and the LS steel ring in a 3D statement. The change in the gas dynamic force in time, calculated
through the pressure field acting directly on the rotor [
        <xref ref-type="bibr" rid="ref1 ref2 ref3 ref4 ref5 ref6">1-6</xref>
        ], was monitored. The rotor motion of the
was described by the harmonic law with given amplitude and frequency [
        <xref ref-type="bibr" rid="ref10 ref11 ref7 ref8 ref9">7-11</xref>
        ]. When assigning the
rotor oscillation frequency, the blade frequency and the rotational frequency were taken into account.
The rotor geometric characteristics were varied and their influence on the wave processes occurring in
the labyrinth seals during rotor vibrations was estimated. At the LS inlet and outlet, a differential
pressure was set. Particular attention was paid to evaluation of wave processes in the circumferential
direction [
        <xref ref-type="bibr" rid="ref12 ref13 ref14 ref15 ref16 ref17 ref18">12-18</xref>
        ]. As a parameter of the investigated dependences, the rotor diameter in the LS zone is
adopted.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. The solid staging</title>
      <p>The solid model (basic variant) is a circular gas-dynamic gap of LS with width 37 mm, height 500 μm
and diameter 422 mm. The solid-state LS ring is a complex stepped construction (Figure 1). The
construction material is steel, the pressure drop is not taken into account, the working medium is
perfect gas.</p>
      <p>The plan for carrying out the computational experiment on studying resonances in the 2FSI subsystem
“labyrinth seals” (LS) of GTU is shown in Table 1. The parameters vary with respect to the “basic
variant”.
(BV1) - basic variant</p>
      <p>The calculation lasts 8.5 hours when examining the perfect gas and absence of pressure drop.
Preliminary calculations have shown that the calculation time increases and reaches 48 hours for
variants corresponding to real gases and presence of pressure drop in the LS. The amount of
computational resources required is limited by the ANSYS capabilities for parallelizing the durability
part of 2FSI calculations. Acceleration of calculations of the structure stress-strain state in the 2FSI
statement was observed when only up to 8 cores (1 node) were used. With a further increase in the
number of calculation cores, the count rate did not increase.</p>
    </sec>
    <sec id="sec-3">
      <title>3. The mesh model</title>
      <p>The minimum size of the finite elements of the mesh models of the LS ring and the gas-dynamic gap
of the LS is 3.75 mm, the element type is hexahedral (Figure 2). For the feasibility of 2FSI
calculations in the construction of mesh models of the gas dynamic region and durability region, it was
ensured that the corresponding network elements coincide at the data transfer boundary from one type
of analysis to another.
a)
b)
Figure 2. The mesh model of the 2FSI subsystem of the LS CC GTU: (a) - mesh model of the
labyrinth seal ring; (b) - mesh model of the gas-dynamic region of the labyrinth seal</p>
    </sec>
    <sec id="sec-4">
      <title>4. Preparation and carrying out calculations on estimation of influence of LS geometrical characteristics on circumferential oscillations in 2FSI subsystem LS CC GTU</title>
      <p>In accordance with the plan, the influence of the LS geometric characteristics on circumferential
oscillations in the 2FSI subsystem of the LS GTU is investigated. Time dependences of pressure, gas
dynamic force, work of gas dynamic force and total displacements in the LS ring LU on the rotor
diameter are obtained (Figure 3).
d)</p>
      <p>Figure 3. Time dependences of the 2FSI parameters of the LS CC GTU subsystem on the rotor
diameter: (a) pressure in the LS gas-dynamic gap; (b) gas-dynamic force; (c) work of the gas dynamic
force; (d) total displacements in the LS ring.
e)</p>
      <p>Figure 4. Analysis of the time dependences of the LS subsystem 2FSI parameters CC GTU on the
rotor diameter: (a) - maximum amplitude of the braking pressure; (b) - maximum amplitude of the gas
dynamic force; (c) - maximum amplitude of the gas dynamic force; (d) - maximum amplitude of the
total displacements; (e) - total work of the gas dynamical force.</p>
      <p>When analyzing time dependences of the 2FSI subsystem of the LS CC GTU on the LS geometric
parameters (Figure 4), it was revealed:</p>
      <p>- with increase in the shaft diameter from 294 mm to 550 mm, the maximum amplitude of the gas
dynamic force is observed at the diameter of 550 mm and equal to 28.9 kN, the minimum value of the
gas dynamic force amplitude is observed at the diameter of 294 mm and equal to 26.6 kN;
- the maximum amplitude frequency of the gas dynamic force is constant and equal to 1.5 kHz for
the considered range of diameters and is not equal to the frequency of the driving force (0.95 kHz);
- with an increase in diameter from 294 mm to 550 mm, the maximum amplitude of the oscillation
of the braking pressure does not change;</p>
      <p>- with increase in diameter from 294 mm to 550 mm, the maximum amplitude of total
displacements does not change;</p>
      <p>- with increase in diameter from 294 mm to 550 mm (by 1.9 times), increase in the maximum
amplitude of the work done by gas dynamic force from 1.45 J to 2.14 J (1.5 times) is observed;
- in figure 4a, it is seen that at the diameter of 294 mm, the pressure oscillation amplitude is
maximal. Spectral analysis showed that the diameter equal to 294 mm corresponds to an oscillatory
process with a resonant frequency of 1800 Hz, which is multiple of the rotor oscillation frequency (2
times). It should be clarified that the rotor oscillation frequency was set equal to the blade frequency.</p>
      <p>The labyrinth seals with diameters equal to 422 mm and 550 mm correspond to oscillatory
processes with resonant frequencies of 1260 Hz and 720 Hz, respectively, which are not multiple of
the exciting frequency. In this case, the pressure oscillation amplitude at these diameters is reduced
(Figure 4a.)
5. Conclusions
1. The influence of the DN diameter on the character of the oscillatory process in the LS is revealed.</p>
      <p>2. If the frequencies of oscillations of the rotor and gas coincide (in the circumferential direction),
the maximum amplitude of the braking pressure in the LS is observed.</p>
      <p>3. The total work of gas-dynamic forces is positive in the considered range of diameters, which
indicates the possibility of “build-up” of the system (Figure 4e).</p>
      <p>4. To reduce the amplitude of the oscillatory process, it is necessary to "detune" from the
coincidence of the oscillation frequencies of the rotor and the gas dynamic volume of the LS by
changing the diameter of the LS at the design stage.</p>
    </sec>
    <sec id="sec-5">
      <title>Acknowledgments</title>
      <p>This study was supported by the Russian Science Foundation, project No 14-19-00877.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Butymova</surname>
            <given-names>L N</given-names>
          </string-name>
          ,
          <article-title>Modorskii V Ya and Petrov V Yu 2016 Numerical modeling of interaction in the dynamic system “gas-structure” with harmonic motion of the piston in the variable section pipe</article-title>
          <source>ICMAR AIP Conference Proceedings 1770 030103</source>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Butymova L N and Modorskii</surname>
            <given-names>V Ya</given-names>
          </string-name>
          <year>2016</year>
          <article-title>One-way FSI simulation of the phase and the geometric parameters of the model of compressor blades on the oscillating gas-dynamic processes pipe MATEC Web Conf</article-title>
          .
          <volume>75</volume>
          <fpage>04007</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>Mekhonoshina</surname>
            <given-names>E V</given-names>
          </string-name>
          and
          <string-name>
            <surname>Modorskii</surname>
            <given-names>V Ya</given-names>
          </string-name>
          <year>2016</year>
          <article-title>Impact of magnetic suspension stiffness on aeroelastic compressor rotor vibrations of gas pumping units</article-title>
          <source>ICMAR AIP Conference Proceedings 1770 030113</source>
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Makarov</surname>
            <given-names>A A</given-names>
          </string-name>
          and
          <string-name>
            <surname>Zaitsev N N 2015</surname>
          </string-name>
          <article-title>Engineering and theoretical problems of using labyrinth seals in high-speed rotary machines Bulletin of PNRPU</article-title>
          . Aerospace engineering
          <volume>42</volume>
          <fpage>61</fpage>
          -
          <lpage>81</lpage>
          (in Russian)
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Butymova L N and Modorskii</surname>
            <given-names>V Ya</given-names>
          </string-name>
          <year>2017</year>
          <article-title>Development and application of a unified algorithm for solving the interdisciplinary problem of modeling aeroelastic processes in the labyrinth seal of centrifugal compressors AIP</article-title>
          <source>Conference Proceedings 1893 030067</source>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Modorskii</surname>
            <given-names>V Y</given-names>
          </string-name>
          and
          <string-name>
            <surname>Shevelev N A 2016</surname>
          </string-name>
          <article-title>Research of aerohydrodynamic and aeroelastic processes on PNRPU HPC system</article-title>
          <source>ICMAR AIP Conference Proceedings 1770 020001</source>
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Akzholov</surname>
            <given-names>M Zh</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Alabuzhev P M</surname>
            , Aliev
            <given-names>A V</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Artistsheva L M and Boyko L G 1995</surname>
          </string-name>
          <article-title>Numerical investigation of actual problems of mechanical engineering and mechanics of continuous and granular media by the method of large particles. Investigation of actual problems of mechanics and machine building (Moscow: National Academy of Applied Sciences, International Association of Developers and Users of the Large-Particle Method</article-title>
          ) p
          <fpage>1658</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <surname>Mekhonoshina</surname>
            <given-names>E V</given-names>
          </string-name>
          and
          <string-name>
            <surname>Modorskii</surname>
            <given-names>V Y</given-names>
          </string-name>
          <year>2015</year>
          <article-title>On a phase-shift of waves at the middle interface</article-title>
          <source>Computer Optics</source>
          <volume>39</volume>
          <fpage>385</fpage>
          -
          <lpage>391</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Butymova L N and Modorskii</surname>
            <given-names>V Ya</given-names>
          </string-name>
          <year>2017</year>
          <article-title>Numerical modeling of the labyrinth seal taking into account vibrations of the gas transmittal unit rotor in aeroelastic formulation</article-title>
          <source>CEUR Workshop Proceedings</source>
          <volume>1902</volume>
          <fpage>10</fpage>
          -
          <lpage>17</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <surname>di Mare</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Imregun</surname>
            <given-names>M</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Green J S and Sayma A I 2010</surname>
          </string-name>
          <article-title>A numerical study of labyrinth seal flutter Journal of tribology-transactions of the ASME 132(2) 022201</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <surname>Zhang</surname>
            <given-names>H</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jia</surname>
            <given-names>XY</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Pan</surname>
            <given-names>XJ</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jiang</surname>
            <given-names>B</given-names>
          </string-name>
          and
          <string-name>
            <surname>Zheng</surname>
            <given-names>Q 2016</given-names>
          </string-name>
          <article-title>Interaction between rotor and annular seals: interlaced and straight-through labyrinth seals Journal of propulsion</article-title>
          and power
          <volume>32</volume>
          (
          <issue>6</issue>
          )
          <fpage>1483</fpage>
          -
          <lpage>1493</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <surname>Vannini</surname>
            <given-names>G</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bertoneri</surname>
            <given-names>M</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nielsen K K</surname>
            , Iudiciani
            <given-names>P</given-names>
          </string-name>
          and
          <string-name>
            <surname>Stronach</surname>
            <given-names>R 2016</given-names>
          </string-name>
          <article-title>Experimental results and computational fluid dynamics simulations of labyrinth and pocket damper seals for wet gas compression Journal of engineering for gas turbines and power-transactions of the ASME 138(5) 052501</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <surname>Dairien</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thouverez</surname>
            <given-names>F</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Blanc</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Helies</surname>
            <given-names>P</given-names>
          </string-name>
          and
          <article-title>Dehouve J Fluid-structure interaction in a labyrinth gas seal coupled to a flexible stator Topics in Modal Analysis</article-title>
          &amp;
          <source>Testing</source>
          <volume>10</volume>
          <fpage>1</fpage>
          -
          <lpage>9</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <surname>Zhang</surname>
            <given-names>E J</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jiao</surname>
            <given-names>Y H</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chen Z</surname>
            <given-names>B</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mo</surname>
            <given-names>W C</given-names>
          </string-name>
          and
          <article-title>Wang S Parametric influences on the nonlinear dynamic responses of a rotor-bearing-foundation-labyrinth seal system Proceedings of the ASME international mechanical engineering congress and exposition 4a UNSP V04AT05A032</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>A V</given-names>
            <surname>Srinivasan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R A</given-names>
            <surname>Arnoldi</surname>
          </string-name>
          and
          <string-name>
            <given-names>A J</given-names>
            <surname>Dennis 1984</surname>
          </string-name>
          <article-title>Aeroelastic instabilities in labyrinth air seal systems</article-title>
          <source>Mechanical Engineering</source>
          <volume>106</volume>
          (
          <issue>6</issue>
          )
          <fpage>91</fpage>
          -
          <lpage>96</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [16]
          <string-name>
            <surname>D R Abbott 1981</surname>
          </string-name>
          <article-title>Advances in labyrinth seal aeroelastic instability prediction and prevention Journal of Engineering for Power-Transactions of the ASME 103(2</article-title>
          )
          <fpage>308</fpage>
          -
          <lpage>312</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          [17]
          <string-name>
            <surname>Phibel</surname>
            <given-names>R</given-names>
          </string-name>
          <article-title>and di Mare L 2012 Comparison between a CFD code and a three control-volume model for labyrinth seal flutter predictions</article-title>
          <source>Proceedings of the ASME Turbo Expo</source>
          <volume>6</volume>
          <fpage>597</fpage>
          -
          <lpage>606</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          [18]
          <string-name>
            <surname>Phibel</surname>
            <given-names>R</given-names>
          </string-name>
          ,
          <string-name>
            <surname>di Mare</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Green J S and Imregun</surname>
            <given-names>M 2009</given-names>
          </string-name>
          <article-title>Numerical investigation of labyrinth seal aeroelastic stability</article-title>
          <source>Proceedings of the ASME Turbo Expo</source>
          <volume>6</volume>
          <fpage>573</fpage>
          -
          <lpage>583</lpage>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>