(20 intermediate revisions by the same user not shown)
Line 36: Line 36:
 
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image1.png|600px]]
 
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image1.png|600px]]
 
|-
 
|-
|style="text-align:center;font-size:85%;"| 1-housing case, 2-cylinder, 3-rotor, 4-rotor bearing, 5-crankshaft, 6-sliding vane
+
|style="text-align:center;font-size:85%;padding-bottom:5px;"| 1-housing case, 2-cylinder, 3-rotor, 4-rotor bearing, 5-crankshaft, 6-sliding vane
 
|- style="text-align: center; font-size: 75%;"
 
|- style="text-align: center; font-size: 75%;"
 
| colspan="1" style="padding-bottom:10px;"| '''Figure 1'''. Structure of the proposed ERMP
 
| colspan="1" style="padding-bottom:10px;"| '''Figure 1'''. Structure of the proposed ERMP
Line 50: Line 50:
 
| colspan="1" style="padding-bottom:10px;"| '''Figure 2'''. Schematic diagram of the working process
 
| colspan="1" style="padding-bottom:10px;"| '''Figure 2'''. Schematic diagram of the working process
 
|}
 
|}
 
  
 
==3. Kinematic model for the sliding vane and rotor==
 
==3. Kinematic model for the sliding vane and rotor==
 
 
===3.1 Velocity and acceleration of the sliding vane===
 
===3.1 Velocity and acceleration of the sliding vane===
  
The rotor driven by the crankshaft, which rotates at a uniform speed, exhibits a composite planar motion. Simultaneously, the sliding vane driven by the rotor reciprocates along the vane slot. In order to analyze the kinematic characteristics of the sliding vane and the rotor, a fixed coordinate system and a rotational coordinate system have been established. Shown as in Figure 3, the fixed coordinate system ''XOY'' locates the origin point at the centrum of pump’s cylinder and Y-axis crossing through the centerline of sliding vane. The rotational coordinate system ''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>'' locates the origin point at the centrum of the rotor and ''Y<sub>1</sub>''-axis crossing the center of round end of the sliding vane. The established rotational coordinate system ''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>'' rotates synchronously with the rotor based on the radius of the eccentric distance of the crankshaft.
+
The rotor driven by the crankshaft, which rotates at a uniform speed, exhibits a composite planar motion. Simultaneously, the sliding vane driven by the rotor reciprocates along the vane slot. In order to analyze the kinematic characteristics of the sliding vane and the rotor, a fixed coordinate system and a rotational coordinate system have been established. Shown as in Figure 3, the fixed coordinate system <math>X O Y</math> locates the origin point at the centrum of pump’s cylinder and <math>Y</math>-axis crossing through the centerline of sliding vane. The rotational coordinate system <math>X_1 O_1 Y_1</math> locates the origin point at the centrum of the rotor and <math>Y_1</math>-axis crossing the center of round end of the sliding vane. The established rotational coordinate system <math>X_1 O_1 Y_1</math> rotates synchronously with the rotor based on the radius of the eccentric distance of the crankshaft.
  
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
Line 73: Line 71:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| style="text-align: center;" | <math display="inline">l_s=e(1-cos\theta )+\frac{e^2}{4r}(1-cos2\theta )</math>
+
| style="text-align: center;" | <math>l_s=e(1-cos\theta )+\frac{e^2}{4r}(1-cos2\theta )</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
 
|}
 
|}
  
where, ''e'' is the eccentric distance of the crankshaft and ''r'' is the radius of the rotor. Shown as in Figure 2, the minimum extension length of the sliding vane is zero when the crankshaft angle, ''θ'', equals to zero or ''2π'' while its maximum value is ''2e'' when the crankshaft angle equals to ''π''. By differentiating the equation (1), the velocity of the sliding vane, ''v<sub>s</sub>'', and acceleration of the sliding vane, ''a<sub>s</sub>'', can be expressed as,
+
where <math>e</math> is the eccentric distance of the crankshaft and <math>r</math> is the radius of the rotor. Shown as in Figure 2, the minimum extension length of the sliding vane is zero when the crankshaft angle, <math>\theta</math>, equals to zero or <math>2\pi</math> while its maximum value is <math>2e</math> when the crankshaft angle equals to <math>\pi</math>. By differentiating the equation (1), the velocity of the sliding vane, <math>v_s</math>, and acceleration of the sliding vane, <math>a_s</math>, can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 100: Line 98:
 
|}
 
|}
  
where, ''ω'' is the angular speed of the crankshaft which can be obtained as ''ω=/30, where n ''is the rotation rate of the pump. Shown as in Figure 3, the round end of the sliding vane that hinged into the rotor rotates reversely along with the rotor. Its relative angular velocity, ''ω<sub>r</sub>'', can be expressed as,
+
where <math>\omega</math> is the angular speed of the crankshaft which can be obtained as <math>\omega =n\pi /30</math>, where <math>n</math> is the rotation rate of the pump. Shown as in Figure 3, the round end of the sliding vane that hinged into the rotor rotates reversely along with the rotor. Its relative angular velocity, <math>\omega_r</math>, can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 112: Line 110:
 
|}
 
|}
  
where, ''d'' is the distance from ''O<sub>1</sub>'' to point ''A'' , and ''β'' is the intersection angle between the center line of the sliding vane and line ''AO<sub>1</sub>''. The value of ''β'' is vary with crankshaft angle ''θ'' and it can be obtained as,
+
where <math>d</math> is the distance from <math>O_1</math> to point <math>A</math>, and <math>\beta</math> is the intersection angle between the center line of the sliding vane and line <math>AO_1</math>. The value of <math>\beta</math> is vary with crankshaft angle <math>\theta</math> and it can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 126: Line 124:
 
===3.2 Velocity and acceleration of the rotor===
 
===3.2 Velocity and acceleration of the rotor===
  
The rotor of the pump is restrained by both the sliding vane and crankshaft. The movements of the pump include the rotation around the center of the cylinder ''O'' and rotation around its own center ''O<sub>1</sub>''. As the sliding vane and the rotor hinged at the round end, the absolute velocity of the rotor at the round end equals to the velocity of the sliding vane. Based on the rotation of the entrained movement of the rotor, the direction of the entrained velocity, ''v<sub>e</sub>'', is vertical to ''OO<sub>1</sub>'' and can be obtained as,
+
The rotor of the pump is restrained by both the sliding vane and crankshaft. The movements of the pump include the rotation around the center of the cylinder <math>O</math> and rotation around its own center <math>O_1</math>. As the sliding vane and the rotor hinged at the round end, the absolute velocity of the rotor at the round end equals to the velocity of the sliding vane. Based on the rotation of the entrained movement of the rotor, the direction of the entrained velocity, <math>v_e</math>, is vertical to <math>OO_1</math> and can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 139: Line 137:
  
  
The movement at point ''A'' is rotating around the rotational coordinate system ''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>''. Then, the relative velocity at point ''A'' is vertical to ''O<sub>1</sub>A'' and its value can be obtained as,
+
The movement at point <math>A</math> is rotating around the rotational coordinate system <math>X_1 O_1 Y_1</math>. Then, the relative velocity at point <math>A</math> is vertical to <math>O_1 A</math> and its value can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 151: Line 149:
 
|}
 
|}
  
where, ''r<sub>0</sub>'' is the radius of the round end of the sliding vane, ''w<sub>p</sub>'' is the rotor’s angular velocity around its own center. According to the velocity vector at point ''A'', the relative velocity can be expressed as,
+
where <math>r_0</math> is the radius of the round end of the sliding vane, <math>w_p</math> is the rotor’s angular velocity around its own center. According to the velocity vector at point <math>A</math>, the relative velocity can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 163: Line 161:
 
|}
 
|}
  
Then, the angular velocity,'' w<sub>p</sub>'', can be expressed as,
+
Then, the angular velocity <math>w_p</math> can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 175: Line 173:
 
|}
 
|}
  
The coordinate value of point on the head surface of the rotor can be expressed under the rotational coordinate system ''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>'' as,
+
The coordinate value of point on the head surface of the rotor can be expressed under the rotational coordinate system <math>X_1 O_1 Y_1</math> as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 190: Line 188:
 
|}
 
|}
  
where, ''l<sub>B</sub>'' is the distance between point on the head surface and the center of the rotor, ''α'' is the intersection angle between the coordinate axis ''x<sub>1</sub>'' and line ''BO<sub>1</sub>''.
+
where <math>l_B</math> is the distance between point on the head surface and the center of the rotor, <math>\alpha</math> is the intersection angle between the coordinate axis <math>x_1</math> and line <math>BO_1</math>.
  
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
Line 200: Line 198:
  
  
Based on the geometry relationship between the fixed coordinate system ''XOY'' and the rotational coordinate system ''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>'', the coordinate value of point on the head surface of the rotor can be expressed under the fixed coordinate system ''XOY'' as,
+
Based on the geometry relationship between the fixed coordinate system ''XOY'' and the rotational coordinate system <math>X_1 O_1 Y_1</math>, the coordinate value of point on the head surface of the rotor can be expressed under the fixed coordinate system <math>X O Y</math> as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 207: Line 205:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| style="text-align: center;" | <math display="inline">\left\{ \begin{array}{c}
+
| style="text-align: center;" | <math display="inline">\left\{ \begin{array}{l}
 
x=l_Bcos\left(\alpha +\beta \right)+esin\theta \\
 
x=l_Bcos\left(\alpha +\beta \right)+esin\theta \\
 
y=l_Bsin\left(\alpha +\beta \right)+ecos\theta  
 
y=l_Bsin\left(\alpha +\beta \right)+ecos\theta  
Line 222: Line 220:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| style="text-align: center;" | <math display="inline">\left\{ \begin{array}{c}
+
| style="text-align: center;" | <math display="inline">\left\{ \begin{array}{l}
 
x=\omega \left(ecos\theta -l_Bsin(\alpha +\beta )\frac{ecos\theta }{\sqrt{d^2-e^2{sin}^2\theta }}\right)\\
 
x=\omega \left(ecos\theta -l_Bsin(\alpha +\beta )\frac{ecos\theta }{\sqrt{d^2-e^2{sin}^2\theta }}\right)\\
 
y=\omega \left(-esin\theta -l_Bcos(\alpha +\beta )\frac{ecos\theta }{\sqrt{d^2-e^2{sin}^2\theta }}\right)
 
y=\omega \left(-esin\theta -l_Bcos(\alpha +\beta )\frac{ecos\theta }{\sqrt{d^2-e^2{sin}^2\theta }}\right)
Line 239: Line 237:
 
===4.1 Force model for sliding vane===
 
===4.1 Force model for sliding vane===
  
The forces acting on the sliding vane are defined to be positive when their directions are the same as the positive direction of the ''X''-axis or ''Y''-axis of the established coordinate system ''XOY'', shown as in Figure 5. The forces acting on the sliding vane include:
+
The forces acting on the sliding vane are defined to be positive when their directions are the same as the positive direction of the <math>X</math>-axis or <math>Y</math>-axis of the established coordinate system <math>XOY</math>, shown as in Figure 5. The forces acting on the sliding vane include:
  
* The normal force (''F<sub>n</sub>'') and tangential force (''F<sub>t</sub>'') at the round end of the sliding vane.
+
* The normal force (<math>F_n</math>) and tangential force (<math>F_t</math>) at the round end of the sliding vane.
  
* The binding forces (''F<sub>R1</sub>'' and ''F<sub>R2</sub>'') caused by the groove that the sliding vane inserted.
+
* The binding forces (<math>F_{R1}</math> and <math>F_{R2}</math>) caused by the groove that the sliding vane inserted.
  
* The friction forces (''F<sub>f1</sub>'' and ''F<sub>f2</sub>'') caused by the movement of the sliding vane.
+
* The friction forces (<math>F_{f1}</math> and <math>F_{f2}</math>) caused by the movement of the sliding vane.
  
* The inertia force (''F<sub>Iv</sub>'') of the sliding vane.
+
* The inertia force (<math>F_{Iv}</math>) of the sliding vane.
  
* The force (''F<sub>pb</sub>'') caused by the pressure difference between the intake chamber and the discharge chamber.
+
* The force (<math>F_{pb}</math>) caused by the pressure difference between the intake chamber and the discharge chamber.
  
  
Line 273: Line 271:
 
|}
 
|}
  
where ''u<sub>s</sub>'' is the friction coefficient between the sliding vane and groove. The pressure difference induced force acing on the sliding vane can be obtained as,
+
where <math>u_{s}</math> is the friction coefficient between the sliding vane and groove. The pressure difference induced force acing on the sliding vane can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 285: Line 283:
 
|}
 
|}
  
where ''H'' is the thickness of the sliding vane, ''l'' is the extension length of the sliding vane, ''p<sub>c</sub>'' is the discharge pressure of the pump, ''p<sub>d</sub>'' is the intake  pressure of the pump.
+
where <math>H</math> is the thickness of the sliding vane, <math>l</math> is the extension length of the sliding vane, <math>p_{c}</math> is the discharge pressure of the pump, <math>p_{d}</math> is the intake  pressure of the pump.
  
 
Based on the rotation of the rotor, the tangential force at the round end of the sliding vane can be obtained as,
 
Based on the rotation of the rotor, the tangential force at the round end of the sliding vane can be obtained as,
Line 299: Line 297:
 
|}
 
|}
  
where ''r'' is the radius of the rotor, ''r<sub>0</sub>'' is the radius of the round end of the sliding vane, ''M<sub>c</sub>'' is the friction torque between the rotor and crankshaft, ''M<sub>JZ</sub>'' is the inertia moment acting on the rotor.
+
where <math>r</math> is the radius of the rotor, <math>r_{0}</math> is the radius of the round end of the sliding vane, <math>M_{c}</math> is the friction torque between the rotor and crankshaft, <math>M_{JZ}</math> is the inertia moment acting on the rotor.
  
 
According to the conditions of force equilibrium and moment equilibrium on the sliding vane with a fixed axis, the equilibrium equations of the sliding vane can be expressed as,
 
According to the conditions of force equilibrium and moment equilibrium on the sliding vane with a fixed axis, the equilibrium equations of the sliding vane can be expressed as,
Line 317: Line 315:
 
|}
 
|}
  
where ''l<sub>0</sub>'' is the length of the sliding vane.
+
where <math>l_{0}</math> is the length of the sliding vane.
  
 
===4.2 Force model for rotor===
 
===4.2 Force model for rotor===
  
The force acting on the rotor include the normal force and tangential force induced by the sliding vane, force (''F<sub>pg</sub>'') induced by the pressure difference of the multiphase flow, inertia force (''F<sub>Ip</sub>'') of the rotor.
+
The force acting on the rotor include the normal force and tangential force induced by the sliding vane, force (<math>F_{pg}</math>) induced by the pressure difference of the multiphase flow, inertia force (<math>F_{Ip}</math>) of the rotor.
  
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
Line 343: Line 341:
 
|}
 
|}
  
The resistance moment (''M<sub>pg</sub>'') caused by the force mentioned above can be expressed as,
+
The resistance moment (<math>M_{pg}</math>) caused by the force mentioned above can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 367: Line 365:
 
|}
 
|}
  
where ''m'' is the weight of the rotor.
+
where <math>m</math> is the weight of the rotor.
  
In this research, the friction torque (''M<sub>c</sub>'') between the rotor and crankshaft is considered. Based on the lubrication theory, the friction torque between the rotor and crankshaft can be obtained as,
+
In this research, the friction torque (<math>M_{c}</math>) between the rotor and crankshaft is considered. Based on the lubrication theory, the friction torque between the rotor and crankshaft can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 381: Line 379:
 
|}
 
|}
  
where ''R<sub>e</sub>'' is the radius of the crankshaft, ''δ<sub>e</sub>'' is the interstice between the rotor and crankshaft, ''l<sub>e</sub>'' is the length of the crankshaft, ''μ'' is the frictional coefficient.
+
where <math>R_{e}</math> is the radius of the crankshaft, <math>\delta_{e}</math> is the interstice between the rotor and crankshaft, <math>l_{e}</math> is the length of the crankshaft, <math>\mu</math> is the frictional coefficient.
  
The resistance moment (''M<sub>n</sub>'') caused by the normal force at round end of the sliding vane can be expressed as,
+
The resistance moment (<math>M_{n}</math>) caused by the normal force at round end of the sliding vane can be expressed as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 407: Line 405:
 
|}
 
|}
  
where ''F<sub>x</sub>'' and ''F<sub>y</sub>'' are the radial components on the X-axis and Y-axis, and can be obtained as,
+
where <math>F_{x}</math> and <math>F_{y}</math> are the radial components on the <math>X</math>-axis and <math>Y</math>-axis, and can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 422: Line 420:
 
|}
 
|}
  
The resistance moment (''M<sub>z</sub>'') acting on the rotor can be obtained as,
+
The resistance moment (<math>M_{z}</math>) acting on the rotor can be obtained as,
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 448: Line 446:
 
==5. Simulation results and discussions==
 
==5. Simulation results and discussions==
  
In this investigation, the ERMP porotype with rated flow of 6 m<sup>3</sup>/h is designed for multiphase transportation, where its main dimensions are given in Table 1. Based on the established model, kinematic characteristics of rotor’s components and forces acting on them can be simulated.
+
In this investigation, the ERMP porotype with rated flow of 6<math>m^3/h</math> is designed for multiphase transportation, where its main dimensions are given in Table 1. Based on the established model, kinematic characteristics of rotor’s components and forces acting on them can be simulated.
  
 
<div class="center" style="font-size: 75%;">'''Table 1'''. Operating specifications and main dimensions of the ERMP porotype
 
<div class="center" style="font-size: 75%;">'''Table 1'''. Operating specifications and main dimensions of the ERMP porotype
Line 458: Line 456:
 
|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|'''Main dimensions'''
 
|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|'''Main dimensions'''
 
|-
 
|-
|  style="border: 1pt solid black;vertical-align: top;"|Rated flow
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Rated flow
|  style="border: 1pt solid black;vertical-align: top;"|6m<sup>3</sup>/h
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|6m<math>m^3/h</math>
|  style="border: 1pt solid black;vertical-align: top;"|Length of sliding vane
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Length of sliding vane
|  style="border: 1pt solid black;vertical-align: top;"|136mm
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|136mm
 
|-
 
|-
|  style="border: 1pt solid black;vertical-align: top;"|Rated speed
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Rated speed
|  style="border: 1pt solid black;vertical-align: top;"|200/min
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|200/min
|  style="border: 1pt solid black;vertical-align: top;"|Width of sliding vane
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Width of sliding vane
|  style="border: 1pt solid black;vertical-align: top;"|20mm
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|20mm
 
|-
 
|-
|  style="border: 1pt solid black;vertical-align: top;"|Working fluid
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Working fluid
|  style="border: 1pt solid black;vertical-align: top;"|Gas-oil multiphase flow
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|Gas-oil multiphase flow
|  style="border: 1pt solid black;vertical-align: top;"|Cylinder radius
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Cylinder radius
|  style="border: 1pt solid black;vertical-align: top;"|90mm
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|90mm
 
|-
 
|-
|  style="border: 1pt solid black;vertical-align: top;"|Intake pressure
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Intake pressure
|  style="border: 1pt solid black;vertical-align: top;"|0.1~0.3MPa
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|0.1<math>\sim</math>0.3MPa
|  style="border: 1pt solid black;vertical-align: top;"|Rotor radius
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Rotor radius
|  style="border: 1pt solid black;vertical-align: top;"|82mm
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|82mm
 
|-
 
|-
|  style="border: 1pt solid black;vertical-align: top;"|Discharge pressure
+
|  style="border: 1pt solid black;vertical-align: top;padding-left:10px;"|Discharge pressure
|  style="border: 1pt solid black;vertical-align: top;"|0.3~3.0MPa
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|0.3<math>\sim</math>3.0MPa
|  style="border: 1pt solid black;;vertical-align: top;"|Radius of vane joint
+
|  style="border: 1pt solid black;;vertical-align: top;padding-left:10px;"|Radius of vane joint
|  style="border: 1pt solid black;vertical-align: top;"|7mm
+
|  style="border: 1pt solid black;vertical-align: top;text-align: center;"|7mm
 
|}
 
|}
  
Line 487: Line 485:
 
===5.1 Simulation results for kinematics of sliding vane===
 
===5.1 Simulation results for kinematics of sliding vane===
  
In this study, the eccentric distance of the crankshaft is set to be 0.007m, 0.008m and 0.009m to evaluate the kinematic characteristics of sliding vane. Figure 7 shows the simulation results of kinematics for sliding vane. From figure 7, it can be noted that the movement, velocity, acceleration of the sliding vane and the angular velocity of the round end of sliding vane increase with the eccentric distance. Shown as in Figure 7 (a), the sliding vane reaches its maximum values when the crankshaft rotates to the ''π'' position. The angular velocity of the round end of the sliding vane has the same variation of the movement of sliding vane, shown as in Figure 7 (d). The velocity of the sliding vane vary with degree of crankshaft in sine form, shown as in Figure 7 (b). On the country, the acceleration of the sliding vane vary with degree of crankshaft in cosine form, shown as in Figure 7 (c).
+
In this study, the eccentric distance of the crankshaft is set to be 0.007m, 0.008m and 0.009m to evaluate the kinematic characteristics of sliding vane. Figure 7 shows the simulation results of kinematics for sliding vane. From figure 7, it can be noted that the movement, velocity, acceleration of the sliding vane and the angular velocity of the round end of sliding vane increase with the eccentric distance. Shown as in Figure 7 (a), the sliding vane reaches its maximum values when the crankshaft rotates to the <math>\pi</math> position. The angular velocity of the round end of the sliding vane has the same variation of the movement of sliding vane, shown as in Figure 7 (d). The velocity of the sliding vane vary with degree of crankshaft in sine form, shown as in Figure 7 (b). On the country, the acceleration of the sliding vane vary with degree of crankshaft in cosine form, shown as in Figure 7 (c).
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
[[Image:Draft_Wang_912933259-image32.png|600px]] </div>
+
|-
 
+
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image32.png|650px]]
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|- style="text-align: center; font-size: 75%;"
Figure 7 Kinematic analysis for sliding vane </div>
+
| colspan="1" style="padding-bottom:10px;"| '''Figure 7'''. Kinematic analysis for sliding vane
 +
|}
  
 
===5.2 Simulation results for kinematics of rotor===
 
===5.2 Simulation results for kinematics of rotor===
  
Based on the established fixed coordinate system (''XOY'') and rotational coordinate system (''X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>''), movements of points on the end surface of rotor can be obtained. Shown as in Figure 8, the motion path of the points are ellipses in different size, except the motion path of the rotor’s center which represents as a circle.
+
Based on the established fixed coordinate system (<math>XOY</math>) and rotational coordinate system (<math>X_1 O_1 Y_1</math>), movements of points on the end surface of rotor can be obtained. Shown as in Figure 8, the motion path of the points are ellipses in different size, except the motion path of the rotor’s center which represents as a circle.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
[[Image:Draft_Wang_912933259-image33.png|294px]] </div>
+
|-
 +
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image33.png|294px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;"| '''Figure 8'''. Movements of points on end surface of rotor
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Figure 8 Movements of points on end surface of rotor</div>
 
  
In this study, four points on the contour of rotor (C<sub>2</sub>, C<sub>3</sub>, C<sub>4</sub> and C<sub>5</sub>) and four points on the middle position of end surface of rotor (B<sub>2</sub>, B<sub>3</sub>, B<sub>4</sub> and B<sub>5</sub>) are selected to simulate velocity of rotor on the end surface, shown as in Figure 8. Figure 9 shows the velocity variations of each selected points on the end surface of rotor. The velocities of each selected points vary with the rotation of the crankshaft in a same trend. The points located at /2 position of the rotor, such as points C<sub>5</sub> and B<sub>5</sub>, have larger velocities compared with the points in other area. Besides, the more the points approach to the center of the rotor, the smaller velocities they have. Then, the area located at opposite of the sliding vane (/2 position of the rotor) and located at contour part of the rotor, such as area of point C<sub>5</sub>, have wear problem during application of the pump and attention need to be paid.
+
In this study, four points on the contour of rotor (<math>C_2</math>, <math>C_3</math>, <math>C_4</math> and <math>C_5</math>) and four points on the middle position of end surface of rotor (<math>B_2</math>, <math>B_3</math>, <math>B_4</math> and <math>B_5</math>) are selected to simulate velocity of rotor on the end surface, shown as in Figure 8. Figure 9 shows the velocity variations of each selected points on the end surface of rotor. The velocities of each selected points vary with the rotation of the crankshaft in a same trend. The points located at <math>3\pi /2</math> position of the rotor, such as points <math>C_5</math> and <math>B_5</math>, have larger velocities compared with the points in other area. Besides, the more the points approach to the center of the rotor, the smaller velocities they have. Then, the area located at opposite of the sliding vane (<math>3\pi /2</math> position of the rotor) and located at contour part of the rotor, such as area of point <math>C_5</math>, have wear problem during application of the pump and attention need to be paid.
  
[[Image:Draft_Wang_912933259-image34.png|600px]]
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 +
|-
 +
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image34.png|600px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;"| '''Figure 9'''. Velocity variations of rotor under different position
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Figure 9 Velocity variations of rotor under different position</div>
 
  
 
Besides, different value of eccentric distance are set to evaluate the effects of eccentric distance on the kinematic behavior of the rotor, shown as in Figure 10. It can noted in Figure 10, the rotor rotates at a faster velocity, a larger angular velocity and angular acceleration when selecting a larger value of the eccentric distance.
 
Besides, different value of eccentric distance are set to evaluate the effects of eccentric distance on the kinematic behavior of the rotor, shown as in Figure 10. It can noted in Figure 10, the rotor rotates at a faster velocity, a larger angular velocity and angular acceleration when selecting a larger value of the eccentric distance.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
[[Image:Draft_Wang_912933259-image35.png|600px]] </div>
+
|-
 
+
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image35.png|600px]]
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|- style="text-align: center; font-size: 75%;"
Figure10 Variation of kinematic behavior of rotor</div>
+
| colspan="1" style="padding-bottom:10px;"| '''Figure 10'''. Variation of kinematic behavior of rotor</
 +
|}
  
 
===5.3 Simulation results for forces on sliding vane===
 
===5.3 Simulation results for forces on sliding vane===
  
Figure 11 shows the forces acting on the sliding vane. It can be noted that the binding force between the groove (''F<sub>R1</sub>'') neutralized most of the pressure difference induced force (''F<sub>pg</sub>''). The tangential force (''F<sub>t</sub>'') acting on the sliding vane is much less than the pressure difference induced force for the reason that it mainly rest on the rotational inertia of the rotor. Besides, the normal force (''F<sub>n</sub>''), which mainly affected by the friction forces (''F<sub>f1</sub>'' and ''F<sub>f2</sub>'') caused by the movement of the sliding vane, increase with rotation of the crankshaft at the first period and decrease at the last period. Moreover, compared with other forces acting on the sliding vane, the binding force (''F<sub>R1</sub>''), pressure difference induced force (''Fpg''), the normal force (''F<sub>n</sub>'') and friction forces (''F<sub>f1</sub>'') are much larger and need to be carefully analyzed.
+
Figure 11 shows the forces acting on the sliding vane. It can be noted that the binding force between the groove (<math>F_{R1}</math>) neutralized most of the pressure difference induced force (<math>F_{pg}</math>). The tangential force (<math>F_{t}</math>) acting on the sliding vane is much less than the pressure difference induced force for the reason that it mainly rest on the rotational inertia of the rotor. Besides, the normal force (<math>F_{n}</math>), which mainly affected by the friction forces (<math>F_{f1}</math> and <math>F_{f2}</math>) caused by the movement of the sliding vane, increase with rotation of the crankshaft at the first period and decrease at the last period. Moreover, compared with other forces acting on the sliding vane, the binding force (<math>F_{R1}</math>), pressure difference induced force (<math>F_{pg}</math>), the normal force (<math>F_{n}</math>) and friction forces (<math>F_{f1}</math>) are much larger and need to be carefully analyzed.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
  [[Image:Draft_Wang_912933259-image36.png|378px]] </div>
+
|-
 +
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image36.png|378px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;"| '''Figure 11'''. Forces acting on the sliding vane
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Figure 11 Forces acting on the sliding vane</div>
 
  
 
Figure 12 shows the frictional force and radial force acting on the rotor. Compared with the frictional force acting on the sliding vane, the frictional force acting on the rotor is much smaller, shown as in Figure 12 (a). The radial force acting on the rotor is also much larger than the frictional force acting on the rotor, shown as in Figure 12 (b). The radial force increase with rotation of the crankshaft until reach its maximum value at π position.
 
Figure 12 shows the frictional force and radial force acting on the rotor. Compared with the frictional force acting on the sliding vane, the frictional force acting on the rotor is much smaller, shown as in Figure 12 (a). The radial force acting on the rotor is also much larger than the frictional force acting on the rotor, shown as in Figure 12 (b). The radial force increase with rotation of the crankshaft until reach its maximum value at π position.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
[[Image:Draft_Wang_912933259-image37.png|600px]] </div>
+
|-
 +
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image37.png|600px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;"| '''Figure 12'''. Frictional force and radial force acting on the rotor
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Figure 12 Frictional force and radial force acting on the rotor</div>
 
  
Figure 13 shows the resistance moment acting on the rotor. It can be noted that the total resistance moment (''M<sub>z</sub>'') is mainly composed by the resistance moment (''M<sub>pg</sub>'') caused by the pressure difference induced force, for the reason that the values of these two resistances are almost the same, shown as in Figure 13 (a). Meanwhile, the values of the resistance moment caused by the normal force at round end of the sliding vane (''M<sub>n</sub>'') and the friction moment between the rotor and crankshaft (''M<sub>c</sub>'') are much smaller than that of the resistance moment (''M<sub>pg</sub>'') caused by the pressure difference induced force, shown as in Figure 13 (b). Then, the resistance moment ''M<sub>n </sub>and M<sub>c </sub>''acting on the rotor can be neglected during the analysis of the pump for the resistance moment is mainly contributed by the'' ''resistance moment of ''M<sub>pg</sub>'' which need carefully analyzed.
+
Figure 13 shows the resistance moment acting on the rotor. It can be noted that the total resistance moment (<math>M_{z}</math>) is mainly composed by the resistance moment (<math>M_{pg}</math>) caused by the pressure difference induced force, for the reason that the values of these two resistances are almost the same, shown as in Figure 13 (a). Meanwhile, the values of the resistance moment caused by the normal force at round end of the sliding vane (<math>M_{n}</math>) and the friction moment between the rotor and crankshaft (<math>M_{c}</math>) are much smaller than that of the resistance moment (<math>M_{pg}</math>) caused by the pressure difference induced force, shown as in Figure 13 (b). Then, the resistance moment <math>M_{n}</math> and <math>M_{c}</math> acting on the rotor can be neglected during the analysis of the pump for the resistance moment is mainly contributed by the resistance moment of <math>M_{pg}</math> which need carefully analyzed.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
[[Image:Draft_Wang_912933259-image38.png|600px]] </div>
+
|-
 
+
|style="padding:10px;"| [[Image:Draft_Wang_912933259-image38.png|600px]]
<div id="_GoBack" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|- style="text-align: center; font-size: 75%;"
Figure 13 Resistance moment acting on the rotor</div>
+
| colspan="1" style="padding-bottom:10px;"| '''Figure 13'''. Resistance moment acting on the rotor
 +
|}
  
 
==6. Conclusion==
 
==6. Conclusion==
  
The structural characteristics and working process of the proposed ERMP are introduced. The mathematic models of the kinematic and dynamic behaviors of key components of the pump are established. Based on the work done in this research, the following results and conclusions were obtained.
+
The structural characteristics and working process of the proposed ERMP are introduced. The mathematic models of the kinematic and dynamic behaviors of key components of the pump are established. Based on the work done in this research, the following results and conclusions were obtained:
 
+
(a) The velocity and acceleration of the sliding vane vary with degree of crankshaft in sine and cosine form, respectively. The sliding vane reaches its maximum value when the crankshaft rotates to the 'π' position.
+
 
+
(b) The areas on the surface of rotor located opposite the sliding vane and apart from the center have larger velocities and wear problem during application of the pump.
+
 
+
(c) The binding force, pressure difference induced force, the normal force and friction forces are much larger than other forces and need to be carefully analyzed.
+
 
+
(d) The resistance moment is mainly contributed by the resistance moment of ''M<sub>pg</sub>'' induced by the pressure difference of the multiphase flow.
+
  
 +
# The velocity and acceleration of the sliding vane vary with degree of crankshaft in sine and cosine form, respectively. The sliding vane reaches its maximum value when the crankshaft rotates to the <math>\pi</math> position.
 +
#The areas on the surface of rotor located opposite the sliding vane and apart from the center have larger velocities and wear problem during application of the pump.
 +
#The binding force, pressure difference induced force, the normal force and friction forces are much larger than other forces and need to be carefully analyzed.
 +
#The resistance moment is mainly contributed by the resistance moment of <math>M_{pg}</math> induced by the pressure difference of the multiphase flow.
  
 
==Acknowledgements==
 
==Acknowledgements==
  
 
The authors acknowledge the support of grant No. 2017QHZ010 from the scientific research-starting project of SWPU China.
 
The authors acknowledge the support of grant No. 2017QHZ010 from the scientific research-starting project of SWPU China.
 
  
 
==References==
 
==References==
 +
<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
  
[1] Porto F., Larson A. 1996. Multiphase Pump Field Trials Demonstrate Practical Applications for the Technology. SPE Prod. Facil., 12(3): 159–164.
+
[1] Porto F., Larson A. Multiphase pump field trials demonstrate practical applications for the technology. SPE Prod. Facil., 12(3):159–164, 1996.
  
[2] Chynoweth E. 1986. Twin Screw Pump is Developed for Multiphase Applications. Process Eng. London, 67(12): 27.
+
[2] Chynoweth E. Twin screw pump is developed for multiphase applications. Process Eng. London, 67(12):27, 1986.
  
[3] Neumann W. 1991. Efficient Multiphase Pump Station for Onshore Application and Prospects for Offshore Application. Proceedings of the Eighth International Pump Users Symposium, Texas A&M Univ., College Station, TX: 43–48.
+
[3] Neumann W. Efficient multiphase pump station for onshore application and prospects for offshore application. Proceedings of the Eighth International Pump Users Symposium, Texas A&M Univ., College Station, TX, 43–48, 1991.
  
[4] Beran W. T. 1994. On the Threshold: Subsea Multiphase Pumping. Proceedings of the 1994 SPE Annual Technical Conference and Exhibition, SPE, New Orleans: 569–581.
+
[4] Beran W.T. On the threshold: subsea multiphase pumping. Proceedings of the 1994 SPE Annual Technical Conference and Exhibition, SPE, New Orleans, 569–581, 1994.
  
[5] Salis D. J., Cordner M., Birnov M. 1998. Multiphase Pumping Comes of Age. World Pumps, 384: 53–54.
+
[5] Salis D.J., Cordner M., Birnov M. Multiphase pumping comes of age. World Pumps, 384:53–54, 1998.
  
[6] Leporcher E., Delaytermoz A., Renault J. F., Gerbier A., Burger O. 2001. Deployment of Multiphase Pumps on a North Sea Field. Proceedings of the 2001 SPE Annual Technical Conference and Exhibition, SPE, New Orleans: 1755–1766.
+
[6] Leporcher E., Delaytermoz A., Renault J.F., Gerbier A., Burger O. Deployment of multiphase pumps on a north sea field. Proceedings of the 2001 SPE Annual Technical Conference and Exhibition, SPE, New Orleans, 1755–1766, 2001.  
  
[7] Dogru A. H., Hamoud A. A., Barlow S. G. 2004. Multiphase Pump Recovers More Oil in a Mature Carbonate Reservoir. J. Petrol. Technol., 56(2): 64–67.
+
[7] Dogru A.H., Hamoud A.A., Barlow S.G. Multiphase pump recovers more oil in a mature carbonate reservoir. J. Petrol. Technol., 56(2):64–67, 2004.
  
[8] Bornemann P. 2009. Multiphase Pumps – the Path to Success. World Pumps, 511: 18–20.
+
[8] Bornemann P. Multiphase pumps – the path to success. World Pumps, 511:18–20, 2009.
  
[9] Cao S., Peng G., Yu Z. 2005. Hydrodynamic Design of Rotodynamic Pump Impeller for Multiphase Pumping by Combined Approach of Inverse Design and CFD Analysis. ASME J. Fluids Eng., 127(2): 330–338.
+
[9] Cao S., Peng G., Yu Z. Hydrodynamic design of rotodynamic pump impeller for multiphase pumping by combined approach of inverse design and cfd analysis. ASME J. Fluids Eng., 127(2):330–338, 2005.
  
[10] Cao F., Peng X., Xing Z., Shu P. 2001. Thermodynamic Performance Simulation of a Twin-Screw Multiphase Pump. Proc. Inst. Mech. Eng., Part E, 215(2): 157–163.
+
[10] Cao F., Peng X., Xing Z., Shu P. Thermodynamic performance simulation of a twin-screw multiphase pump. Proc. Inst. Mech. Eng., Part E, 215(2):157–163, 2001.
  
[11] Chynoweth E. Twin screw pump is developed for multiphase applications. Process Eng (London) 1986; 67(12): 27–31.
+
[11] Chynoweth E. Twin screw pump is developed for multiphase applications. Process Eng. London, 67(12):27–31, 1986.
  
[12] Dal P. D., Larson L. A. 1996. Multiphase pump field trials demonstrate practical applications for the technology. SPE Prod Facilities. 12(3): 159–164.
+
[12] Dal P.D., Larson L.A. Multiphase pump field trials demonstrate practical applications for the technology. SPE Prod. Facilities. 12(3):159–164, 1996.
  
[13] Ryazantsev V.M., Plyasov V.V., Shumakov N. D. 2004. Twin-screw pumps for transferring petroleum-water-gas multiphase mixture or diesel fuel. Chem Petrol Eng, 40(8): 407–410.
+
[13] Ryazantsev V.M., Plyasov V.V., Shumakov N.D. Twin-screw pumps for transferring petroleum-water-gas multiphase mixture or diesel fuel. Chem. Petrol. Eng., 40(8):407–410, 2004.
  
[14] Nakashima C. Y., Oliveira S. 2006. Heat transfer in a twin-screw multiphase pump: thermal modeling and one application in the petroleum industry. Energy, 31(15): 3415–3425.
+
[14] Nakashima C. Y., Oliveira S. Heat transfer in a twin-screw multiphase pump: thermal modeling and one application in the petroleum industry. Energy, 31(15):3415–3425, 2006.
  
[15] Muller L. D., Rohlfing G., Spelter H. 2009. Twin screw pumps help improving oil recovery in mature fields and transfer heavy crude oil over long distances. Oil Gas Eur Mag, 35(3): 127–131.
+
[15] Muller L.D., Rohlfing G., Spelter H. Twin screw pumps help improving oil recovery in mature fields and transfer heavy crude oil over long distances. Oil Gas Eur. Mag., 35(3):127–131, 2009.
 +
 
 +
</div>

Revision as of 13:37, 27 March 2020

Abstract

In order to simplify the manufacture of multiphase pump and improve the operating flexibility for gas void fraction (GVF) of the multiphase fluids that the pump transported, an eccentric rotary multiphase pump (ERMP) is presented. In this study, the structural characteristics and working principle of the ERMP are presented first. Then, the kinematic and force models are established for the key components- sliding vane and rotor. The velocity, acceleration, and force equations with shaft rotation angle are derived for each component. Based on the established models, simulations are performed for an ERMP prototype. The simulated results show that the areas opposite the sliding vane and apart from the center of the rotor have larger velocities and wear problem. Moreover, the binding force, pressure difference induced force and the normal force exercise a negative effect on the friction at the sliding vane sides and rotor. Lower shaft speed and smaller eccentric distance of the crankshaft are helpful to reduce this effect. The findings confirm that the proposed ERMP is suitable for multiphase transportation and has a higher mechanical efficiency for its advanced structure and working principle.

Key words: multiphase pump, ERMP, sliding vane, rotor, kinematic and force behavior

1. Introduction

Multiphase fluid transportation, which enables the transfer of multiphase fluid in a single pipeline, is considered to be an efficient and economical technology in the petrochemical industry. Due to the elimination of conventional transportation facilities such as bulky separators, fluid pumps, gas compressors and extra pipelines, the cost of the multiphase transportation technology is only about 70% of the conventional transportation technology [1].

Moreover, the multiphase transportation technology is especially helpful in the exploration of the petroleum resources at remote areas for this technology provides a simpler and more compact system. More and more petrochemical companies worldwide are attracted to the use of multiphase pumps for boosting multiphase flow and many investigations on multiphase pumps have been carried out [2–8]. According to the boosting principles of developed pumps, the multiphase pumps can be classified as the rotodynamic pump and positive displacement pump. Rotodynamic pump, such as the helicon-axial pump, transform the mechanical input power into kinetic energy of multiphase fluid by displacing it along suitable designed channels within the pumps and then converted into working pressure [9]. The rotodynamic pump can provide a large flow but is not suitable for handling multiphase fluids that gas void fraction (GVF) is more than 50%. Positive displacement pump, such as the twin-screw pump, works on the principle of enclosing a defined volume on the suction side and moving it to the discharge end [10]. The positive displacement pump can handle the multiphase fluids with a large difference in the inlet GVFs. Several multiphase boosting applications demonstrated that positive displacement pump could operate for multiphase duties successfully [11-15]. However, the rotors of the positive displacement pump are geometrically complex and require higher machining precision that results in a higher manufacture cost.

In order to simplify the manufacture of multiphase pump and improve its operating flexibility for GVF, a novel positive displacement multiphase pump names eccentric rotary multiphase pump (ERMP) is presented. The major components of the proposed ERMP are all geometry simple and easy to manufacture. Due to the synchronically movement of the rotor and the cylinder of the pump, the frictional loss between the major components are largely minimized and the mechanical efficiency of the proposed pump is higher than that of the conventional rotary mechanisms.

In this investigation, the structural characteristics and working process of the proposed ERMP are particularly introduced. Its kinematic behaviors as well as the dynamic behaviors are theoretically formulated and discussed. In addition, the effects of the eccentric distance of the crankshaft on the kinematic and dynamic behaviors of key components of the pump are theoretically analyzed.

2. Structural characteristics and working principle

Shown as in Figure 1, the proposed ERMP consists of a rotor, a sliding vane, a housing case, and a crankshaft. The rotor is eccentrically assembly in the cylinder that is assembly with the housing case. Due to the outside wall of the rotor and the inside wall of the cylinder are tangent to each other, a crescent-shaped working chamber is formed by the inside wall of the cylinder and the outside wall of the rotor. The upper end the sliding vane is inserted in the groove of the housing case while the lower end is assembly with the rotor by a joint. In such case, the formed crescent-shaped chamber is divided into an intake chamber and a discharge chamber by the sliding vane. The intake port and discharge port of the pump are connected to the intake chamber and discharge chamber, respectively.

Draft Wang 912933259-image1.png
1-housing case, 2-cylinder, 3-rotor, 4-rotor bearing, 5-crankshaft, 6-sliding vane
Figure 1. Structure of the proposed ERMP


During the operation of the pump, the rotor is driven by the crankshaft, which is rotate clockwise around its geometrical axis, to swing in the cylinder and change the volume of the crescent-shaped working chamber. At the same time, the sliding vane reciprocate straightly in the vane slot. Due to the compound motion of the vane, the volumes of the divided intake chamber and discharge chamber vary gradually, and thus the intake -discharge process is resulted, shown as in Figure 2. The volume of the intake chamber increase with the crankshaft angle of θ and result in a vacuum in this area during the intake process. Meanwhile, the volume of the discharge chamber decrease with the crankshaft angle of θ during the discharge process. In such case, the pump inhales the multiphase fluid continuously into the intake chamber through the intake port of the pump during the intake process and discharge the multiphase fluid out of the pump through the discharge port during the discharge process.

Draft Wang 912933259-image2.png
Figure 2. Schematic diagram of the working process

3. Kinematic model for the sliding vane and rotor

3.1 Velocity and acceleration of the sliding vane

The rotor driven by the crankshaft, which rotates at a uniform speed, exhibits a composite planar motion. Simultaneously, the sliding vane driven by the rotor reciprocates along the vane slot. In order to analyze the kinematic characteristics of the sliding vane and the rotor, a fixed coordinate system and a rotational coordinate system have been established. Shown as in Figure 3, the fixed coordinate system locates the origin point at the centrum of pump’s cylinder and -axis crossing through the centerline of sliding vane. The rotational coordinate system locates the origin point at the centrum of the rotor and -axis crossing the center of round end of the sliding vane. The established rotational coordinate system rotates synchronously with the rotor based on the radius of the eccentric distance of the crankshaft.

Draft Wang 912933259-image3.png
Figure 3. Coordinate systems of the analysis model


The extension length of the sliding vane, ls, can be expressed as,

(1)

where is the eccentric distance of the crankshaft and is the radius of the rotor. Shown as in Figure 2, the minimum extension length of the sliding vane is zero when the crankshaft angle, , equals to zero or while its maximum value is when the crankshaft angle equals to . By differentiating the equation (1), the velocity of the sliding vane, , and acceleration of the sliding vane, , can be expressed as,

(2)
(3)

where is the angular speed of the crankshaft which can be obtained as , where is the rotation rate of the pump. Shown as in Figure 3, the round end of the sliding vane that hinged into the rotor rotates reversely along with the rotor. Its relative angular velocity, , can be expressed as,

(4)

where is the distance from to point , and is the intersection angle between the center line of the sliding vane and line . The value of is vary with crankshaft angle and it can be obtained as,

(5)

3.2 Velocity and acceleration of the rotor

The rotor of the pump is restrained by both the sliding vane and crankshaft. The movements of the pump include the rotation around the center of the cylinder and rotation around its own center . As the sliding vane and the rotor hinged at the round end, the absolute velocity of the rotor at the round end equals to the velocity of the sliding vane. Based on the rotation of the entrained movement of the rotor, the direction of the entrained velocity, , is vertical to and can be obtained as,

(6)


The movement at point is rotating around the rotational coordinate system . Then, the relative velocity at point is vertical to and its value can be obtained as,

(7)

where is the radius of the round end of the sliding vane, is the rotor’s angular velocity around its own center. According to the velocity vector at point , the relative velocity can be expressed as,

(8)

Then, the angular velocity can be expressed as,

(9)

The coordinate value of point on the head surface of the rotor can be expressed under the rotational coordinate system as,

(10)

where is the distance between point on the head surface and the center of the rotor, is the intersection angle between the coordinate axis and line .

Draft Wang 912933259-image14.png
Figure 4. Coordinate systems of the analysis model


Based on the geometry relationship between the fixed coordinate system XOY and the rotational coordinate system , the coordinate value of point on the head surface of the rotor can be expressed under the fixed coordinate system as,

(11)

Then, the velocity of the head surface on the rotor can be obtained by differentiating the above equation as,

(12)

4. Force model for sliding vane and rotor

The force model established in this research include not only the principle forces but also the often-neglected pressure difference induced force between intake chamber and discharge chamber. In these models, we assume the following:

  1. The pressure pulsation of the intake and discharge processes is ignored, which means that the intake and discharge pressures are regarded as constants.
  2. The effect of oil film thickness change on friction coefficient is neglected.

4.1 Force model for sliding vane

The forces acting on the sliding vane are defined to be positive when their directions are the same as the positive direction of the -axis or -axis of the established coordinate system , shown as in Figure 5. The forces acting on the sliding vane include:

  • The normal force () and tangential force () at the round end of the sliding vane.
  • The binding forces ( and ) caused by the groove that the sliding vane inserted.
  • The friction forces ( and ) caused by the movement of the sliding vane.
  • The inertia force () of the sliding vane.
  • The force () caused by the pressure difference between the intake chamber and the discharge chamber.


Draft Wang 912933259-image17.png
Figure 5. Forces acting on the sliding vane


The friction forces acting on the sliding vane can be obtained as,

(13)

where is the friction coefficient between the sliding vane and groove. The pressure difference induced force acing on the sliding vane can be obtained as,

(14)

where is the thickness of the sliding vane, is the extension length of the sliding vane, is the discharge pressure of the pump, is the intake pressure of the pump.

Based on the rotation of the rotor, the tangential force at the round end of the sliding vane can be obtained as,

(15)

where is the radius of the rotor, is the radius of the round end of the sliding vane, is the friction torque between the rotor and crankshaft, is the inertia moment acting on the rotor.

According to the conditions of force equilibrium and moment equilibrium on the sliding vane with a fixed axis, the equilibrium equations of the sliding vane can be expressed as,

(16)

where is the length of the sliding vane.

4.2 Force model for rotor

The force acting on the rotor include the normal force and tangential force induced by the sliding vane, force () induced by the pressure difference of the multiphase flow, inertia force () of the rotor.

Draft Wang 912933259-image24.png
Figure 6. Forces acting on the rotor


Different pressure exist between the intake chamber and discharge chamber of the pump during its application. The pressure difference induced force would loaded on the rotor during the application of the pump, shown as in Figure 6. The value of the force can be obtained as,

(17)

The resistance moment () caused by the force mentioned above can be expressed as,

(18)

The inertia force of the rotor can be obtained as,

(19)

where is the weight of the rotor.

In this research, the friction torque () between the rotor and crankshaft is considered. Based on the lubrication theory, the friction torque between the rotor and crankshaft can be obtained as,

(20)

where is the radius of the crankshaft, is the interstice between the rotor and crankshaft, is the length of the crankshaft, is the frictional coefficient.

The resistance moment () caused by the normal force at round end of the sliding vane can be expressed as,

(21)

Based on the research mentioned above, the radial force acting on the rotor can be obtained as,

(22)

where and are the radial components on the -axis and -axis, and can be obtained as,

(23)

The resistance moment () acting on the rotor can be obtained as,

(24)

Then, the capacity of the multiphase pump can be obtained as,

(25)

5. Simulation results and discussions

In this investigation, the ERMP porotype with rated flow of 6 is designed for multiphase transportation, where its main dimensions are given in Table 1. Based on the established model, kinematic characteristics of rotor’s components and forces acting on them can be simulated.

Table 1. Operating specifications and main dimensions of the ERMP porotype
Operating specifications Main dimensions
Rated flow 6m Length of sliding vane 136mm
Rated speed 200/min Width of sliding vane 20mm
Working fluid Gas-oil multiphase flow Cylinder radius 90mm
Intake pressure 0.10.3MPa Rotor radius 82mm
Discharge pressure 0.33.0MPa Radius of vane joint 7mm


5.1 Simulation results for kinematics of sliding vane

In this study, the eccentric distance of the crankshaft is set to be 0.007m, 0.008m and 0.009m to evaluate the kinematic characteristics of sliding vane. Figure 7 shows the simulation results of kinematics for sliding vane. From figure 7, it can be noted that the movement, velocity, acceleration of the sliding vane and the angular velocity of the round end of sliding vane increase with the eccentric distance. Shown as in Figure 7 (a), the sliding vane reaches its maximum values when the crankshaft rotates to the position. The angular velocity of the round end of the sliding vane has the same variation of the movement of sliding vane, shown as in Figure 7 (d). The velocity of the sliding vane vary with degree of crankshaft in sine form, shown as in Figure 7 (b). On the country, the acceleration of the sliding vane vary with degree of crankshaft in cosine form, shown as in Figure 7 (c).

Draft Wang 912933259-image32.png
Figure 7. Kinematic analysis for sliding vane

5.2 Simulation results for kinematics of rotor

Based on the established fixed coordinate system () and rotational coordinate system (), movements of points on the end surface of rotor can be obtained. Shown as in Figure 8, the motion path of the points are ellipses in different size, except the motion path of the rotor’s center which represents as a circle.

Draft Wang 912933259-image33.png
Figure 8. Movements of points on end surface of rotor


In this study, four points on the contour of rotor (, , and ) and four points on the middle position of end surface of rotor (, , and ) are selected to simulate velocity of rotor on the end surface, shown as in Figure 8. Figure 9 shows the velocity variations of each selected points on the end surface of rotor. The velocities of each selected points vary with the rotation of the crankshaft in a same trend. The points located at position of the rotor, such as points and , have larger velocities compared with the points in other area. Besides, the more the points approach to the center of the rotor, the smaller velocities they have. Then, the area located at opposite of the sliding vane ( position of the rotor) and located at contour part of the rotor, such as area of point , have wear problem during application of the pump and attention need to be paid.

Draft Wang 912933259-image34.png
Figure 9. Velocity variations of rotor under different position


Besides, different value of eccentric distance are set to evaluate the effects of eccentric distance on the kinematic behavior of the rotor, shown as in Figure 10. It can noted in Figure 10, the rotor rotates at a faster velocity, a larger angular velocity and angular acceleration when selecting a larger value of the eccentric distance.

Draft Wang 912933259-image35.png
Figure 10. Variation of kinematic behavior of rotor</

5.3 Simulation results for forces on sliding vane

Figure 11 shows the forces acting on the sliding vane. It can be noted that the binding force between the groove () neutralized most of the pressure difference induced force (). The tangential force () acting on the sliding vane is much less than the pressure difference induced force for the reason that it mainly rest on the rotational inertia of the rotor. Besides, the normal force (), which mainly affected by the friction forces ( and ) caused by the movement of the sliding vane, increase with rotation of the crankshaft at the first period and decrease at the last period. Moreover, compared with other forces acting on the sliding vane, the binding force (), pressure difference induced force (), the normal force () and friction forces () are much larger and need to be carefully analyzed.

Draft Wang 912933259-image36.png
Figure 11. Forces acting on the sliding vane


Figure 12 shows the frictional force and radial force acting on the rotor. Compared with the frictional force acting on the sliding vane, the frictional force acting on the rotor is much smaller, shown as in Figure 12 (a). The radial force acting on the rotor is also much larger than the frictional force acting on the rotor, shown as in Figure 12 (b). The radial force increase with rotation of the crankshaft until reach its maximum value at π position.

Draft Wang 912933259-image37.png
Figure 12. Frictional force and radial force acting on the rotor


Figure 13 shows the resistance moment acting on the rotor. It can be noted that the total resistance moment () is mainly composed by the resistance moment () caused by the pressure difference induced force, for the reason that the values of these two resistances are almost the same, shown as in Figure 13 (a). Meanwhile, the values of the resistance moment caused by the normal force at round end of the sliding vane () and the friction moment between the rotor and crankshaft () are much smaller than that of the resistance moment () caused by the pressure difference induced force, shown as in Figure 13 (b). Then, the resistance moment and acting on the rotor can be neglected during the analysis of the pump for the resistance moment is mainly contributed by the resistance moment of which need carefully analyzed.

Draft Wang 912933259-image38.png
Figure 13. Resistance moment acting on the rotor

6. Conclusion

The structural characteristics and working process of the proposed ERMP are introduced. The mathematic models of the kinematic and dynamic behaviors of key components of the pump are established. Based on the work done in this research, the following results and conclusions were obtained:

  1. The velocity and acceleration of the sliding vane vary with degree of crankshaft in sine and cosine form, respectively. The sliding vane reaches its maximum value when the crankshaft rotates to the position.
  2. The areas on the surface of rotor located opposite the sliding vane and apart from the center have larger velocities and wear problem during application of the pump.
  3. The binding force, pressure difference induced force, the normal force and friction forces are much larger than other forces and need to be carefully analyzed.
  4. The resistance moment is mainly contributed by the resistance moment of induced by the pressure difference of the multiphase flow.

Acknowledgements

The authors acknowledge the support of grant No. 2017QHZ010 from the scientific research-starting project of SWPU China.

References

[1] Porto F., Larson A. Multiphase pump field trials demonstrate practical applications for the technology. SPE Prod. Facil., 12(3):159–164, 1996.

[2] Chynoweth E. Twin screw pump is developed for multiphase applications. Process Eng. London, 67(12):27, 1986.

[3] Neumann W. Efficient multiphase pump station for onshore application and prospects for offshore application. Proceedings of the Eighth International Pump Users Symposium, Texas A&M Univ., College Station, TX, 43–48, 1991.

[4] Beran W.T. On the threshold: subsea multiphase pumping. Proceedings of the 1994 SPE Annual Technical Conference and Exhibition, SPE, New Orleans, 569–581, 1994.

[5] Salis D.J., Cordner M., Birnov M. Multiphase pumping comes of age. World Pumps, 384:53–54, 1998.

[6] Leporcher E., Delaytermoz A., Renault J.F., Gerbier A., Burger O. Deployment of multiphase pumps on a north sea field. Proceedings of the 2001 SPE Annual Technical Conference and Exhibition, SPE, New Orleans, 1755–1766, 2001.

[7] Dogru A.H., Hamoud A.A., Barlow S.G. Multiphase pump recovers more oil in a mature carbonate reservoir. J. Petrol. Technol., 56(2):64–67, 2004.

[8] Bornemann P. Multiphase pumps – the path to success. World Pumps, 511:18–20, 2009.

[9] Cao S., Peng G., Yu Z. Hydrodynamic design of rotodynamic pump impeller for multiphase pumping by combined approach of inverse design and cfd analysis. ASME J. Fluids Eng., 127(2):330–338, 2005.

[10] Cao F., Peng X., Xing Z., Shu P. Thermodynamic performance simulation of a twin-screw multiphase pump. Proc. Inst. Mech. Eng., Part E, 215(2):157–163, 2001.

[11] Chynoweth E. Twin screw pump is developed for multiphase applications. Process Eng. London, 67(12):27–31, 1986.

[12] Dal P.D., Larson L.A. Multiphase pump field trials demonstrate practical applications for the technology. SPE Prod. Facilities. 12(3):159–164, 1996.

[13] Ryazantsev V.M., Plyasov V.V., Shumakov N.D. Twin-screw pumps for transferring petroleum-water-gas multiphase mixture or diesel fuel. Chem. Petrol. Eng., 40(8):407–410, 2004.

[14] Nakashima C. Y., Oliveira S. Heat transfer in a twin-screw multiphase pump: thermal modeling and one application in the petroleum industry. Energy, 31(15):3415–3425, 2006.

[15] Muller L.D., Rohlfing G., Spelter H. Twin screw pumps help improving oil recovery in mature fields and transfer heavy crude oil over long distances. Oil Gas Eur. Mag., 35(3):127–131, 2009.

Back to Top

Document information

Published on 27/03/20
Accepted on 24/03/20
Submitted on 10/03/20

Volume 36, Issue 1, 2020
DOI: 10.23967/j.rimni.2020.03.004
Licence: CC BY-NC-SA license

Document Score

0

Views 160
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?