Review of Wet Gas Compressor Performance
Review of Wet Gas Compressor Performance
Review of Wet Gas Compressor Performance
Performance
Zaenuddin Ar Rofiqi
Petroleum Engineering
Submission date: June 2014
Supervisor: Michael Golan, IPT
Zaenuddin Ar Rofiqi
Petroleum Engineering
Subsea wet gas compression without preliminary separation give an opportunity to develop
mature and remote gas/condensate field. Predicting wet gas compressor performance is
important to develop those fields. Currently there is no standard exist for wet gas
compression performance calculation (Hundseid et al., 2008). The current standards (ISO
5389:2005) exist only for dry gas by utilizing Schultz method.
In wet gas compression, fluid properties are changing as fluid being compressed. Schultz
method may not be accurate to be applied in wet gas compression since it uses average fluid
properties between inlet and outlet of the compressor.
The necessities of standard in wet gas compression encourage this project to conduct a
research in wet gas compression performance. This project will evaluate Huntington direct
integration method (Huntington, 1985) which originally for dry gas to be used in wet gas
compression calculation and compare the result with Schultz method.
In this project, visual basic script was created to calculate wet gas compressor performance
by using direct integration method. Fluid properties are calculated by HYSYS TM software. It
is also desired to know how HYSYS TM calculate wet gas compression. Compressor off
design operation effect on compressor performance is also taken into consideration when
calculating wet gas compressor performance.
The results show that direct integration method is much more suitable for wet gas
compression calculation compared to Schultz method. As wet gas being compressed, more
liquid will form in the mixture and it will affect fluid properties. The presence of heavier
fluid in the mixture and cooling effect is one of the factor that affecting compressor
performance. This factor will reduce specific polytopic head of compressor. Direct
integration method calculation is considering this effect.
Off design operation affect compressor performance calculation significantly. The deviation
between original performance and corrected performance become larger as actual and
reference condition become much more different. At large deviation, the corrected
compressor performances become uncertain.
i
In HYSYSTM software, interaction between liquid and vapor phase is taken into
consideration. The software does not calculate compressor performance by treating each
phase separately.
Further recommendations for this project are to compare calculated performance using direct
integration method with actual experiment and improvement of current compressor off design
correction. By comparing compressor performance calculation with real data, the accuracy of
direct integration method will be discovered. In gas condensate field, fluid properties will
change significantly with decreasing pressure. Therefore, actual and test reference condition
will be different over time and will affect gas condensate field performance. This condition
need attention since current off design correction accuracy will be reduced if the deviation is
large. Good prediction of wet gas compression performance is required in developing gas
condensate field
ii
ACKNOWLEDGEMENT
This project is submitted to the Norwegian University of Science and Technology (NTNU)
for requirement of master thesis in petroleum engineering master degree program.
I would like to express my deepest appreciation to Professor Michael Golan who gives me
the opportunity to be his student and take master thesis.
I would like to thanks to Jesus and Gilberto who help me in the master thesis. In this thesis, I
often stuck in a problem. And they guide me a lot to finish this project.
I would like to thanks also to Mayembe and Milan. Who helps me a lot for one year. They
help me with my semester project and master thesis
I would like to thanks also to Mohamad. He guides me to create visual basic in excel.
Finally, I would like to express a gratitude for my family, my wife, and my friend for
supporting me in this project.
iii
TABLE OF CONTENTS
ABSTRACT ................................................................................................................................ i
LIST OF FIGURE..................................................................................................................... vi
1. BACKGROUND ............................................................................................................. 10
2. THEORY ......................................................................................................................... 12
iv
3.1.1. Fluid properties .................................................................................................. 34
4. RESULT .......................................................................................................................... 41
5. DISCUSSION .................................................................................................................. 69
NOMENCLATURE ................................................................................................................ 76
REFERENCE ........................................................................................................................... 79
APPENDIX .............................................................................................................................. 81
v
LIST OF FIGURE
Figure 1.1 Effect of installing subsea wet gas compression to gas/condensate field
performance (Courtesy of Vinterst, 2014). ............................................................................ 11
Figure 2.3 Cross section of centrifugal compressor (Courtesy of petroleum extension service)
.................................................................................................................................................. 14
Figure 2.4 Work in graphical representation (Courtesy of Moran et al., 2011) ...................... 15
Figure 2.5 Compressor reference processes (Courtesy of ISO 5389, 2005) ............................ 16
Figure 2.6 Comparison between isentropic and actual compression (Courtesy of Moran,
2011) ........................................................................................................................................ 18
Figure 2.9 Effect of wet gas presence to specific polytropic head (Courtesy of Hundseid et al,
2008) ........................................................................................................................................ 26
Figure 2.10 Mass flow corrected specific polytropic head (Courtesy of Hundseid et al, 2008)
.................................................................................................................................................. 27
Figure 2.11 Wood s corrected specific polytropic head (Courtesy of Hundseid et al, 2008) . 27
Figure 2.12 Effect of fluid density to compressor curve (Courtesy of Hanlon, 2001) ............ 28
Figure 2.13 Basic procedure for off design correction by iterative method (Courtesy of Ma et
al, 2013) ................................................................................................................................... 30
Figure 3.1 Process flow diagram of mixing liquid and gas phase ........................................... 35
vi
Figure 4.1 Compression process path at various GVF for case A ........................................... 41
Figure 4.2 Gas volume fraction change at each pressure step i with total number of pressure
step = 10 for case A ................................................................................................................. 42
Figure 4.3 Polytropic exponent change at each pressure step i, case A GVF = 1 ................... 43
Figure 4.4 Polytropic exponent change at each pressure step i, case A GVF = 0.99 .............. 43
Figure 4.5 Polytropic exponent change at each pressure step i, case A GVF = 0.97 .............. 44
Figure 4.6 Polytropic exponent change at each pressure step i, case A GVF = 0.96 .............. 44
Figure 4.7 Calculated head and error vs total pressure step size, case A GVF = 1 ................. 45
Figure 4.8 Calculated head and error vs total pressure step size, case A GVF = 0.99 ............ 46
Figure 4.9 Calculated head and error vs total pressure step size, case A GVF = 0.97 ............ 46
Figure 4.10 Calculated head and error vs total pressure step size, case A GVF = 0.96 .......... 47
Figure 4.11 Calculated head and power vs total pressure step size, case A GVF = 1 ............. 48
Figure 4.12 Calculated head and power vs total pressure step size, case A GVF = 0.99 ........ 48
Figure 4.13 Calculated head and power vs total pressure step size, case A GVF = 0.97 ........ 49
Figure 4.14 Calculated head and power vs total pressure step size, case A GVF = 0.96 ........ 49
Figure 4.15 Percentage of head on each pressure step i, case A GVF = 1 .............................. 50
Figure 4.16 Percentage of head on each pressure step i, case A GVF = 0.99 ......................... 51
Figure 4.17 Percentage of head on each pressure step i, case A GVF = 0.97 ......................... 51
Figure 4.18 Percentage of head on each pressure step i, case A GVF = 0.96 ......................... 52
Figure 4.20 Original vs corrected compressor curve, case B GVF = 0.99 .............................. 53
Figure 4.21 Original vs corrected compressor curve, case B GVF = 0.97 .............................. 54
Figure 4.22 Original vs corrected compressor curve, case B GVF = 0.96 .............................. 54
Figure 4.23 Compression process path at various GVF for case B ......................................... 55
Figure 4.24 Gas volume fraction change at each pressure step i with total number of pressure
step = 10 for case B.................................................................................................................. 56
Figure 4.25 Polytropic exponent change at each pressure step i, case B GVF = 1.................. 57
Figure 4.26 Polytropic exponent change at each pressure step i, case B GVF = 0.99............. 57
vii
Figure 4.27 Polytropic exponent change at each pressure step i, case B GVF = 0.97............. 58
Figure 4.28 Polytropic exponent change at each pressure step i, case B GVF = 0.96............. 58
Figure 4.29 Calculated head and error vs total pressure step size, case B GVF = 1 ............... 59
Figure 4.30 Calculated head and error vs total pressure step size, case B GVF = 0.99 .......... 60
Figure 4.31 Calculated head and error vs total pressure step size, case B GVF = 0.97 .......... 60
Figure 4.32 Calculated head and error vs total pressure step size, case B GVF = 0.96 .......... 61
Figure 4.33 Calculated head and power vs total pressure step size, case B GVF = 1 ............. 62
Figure 4.34 Calculated head and power vs total pressure step size, case B GVF = 0.99 ........ 62
Figure 4.35 Calculated head and power vs total pressure step size, case B GVF = 0.97 ........ 63
Figure 4.36 Calculated head and power vs total pressure step size, case B GVF = 0.96 ........ 64
Figure 4.37 Percentage of head on each pressure step i, case B GVF = 1 ............................... 65
Figure 4.38 Percentage of head on each pressure step i, case B GVF = 0.99 .......................... 65
Figure 4.39 Percentage of head on each pressure step i, case B GVF = 0.97 .......................... 66
Figure 4.40 Percentage of head on each pressure step i, case B GVF = 0.96 .......................... 66
Figure 5.1 Compression process path case A with bubble point and dew point line .............. 70
Figure 5.2 Compression process path case B with bubble point and dew point line ............... 70
Figure 5.3 Polytropic exponent in a p-v diagram (Courtesy of Moran et al., 2010)................ 71
Figure 5.4 Cooling effect in wet gas compression for case A ................................................. 72
viii
LIST OF TABLES
Table 1 Separator gas and liquid compositions (Courtesy of McCain, 1990) ......................... 34
Table 10 Comparison of wet gas compression using separator and without separator (case
A) .............................................................................................................................................. 68
Table 11 Comparison of wet gas compression using separator and without separator (case
B) .............................................................................................................................................. 68
ix
1. BACKGROUND
As technology developed and giant gas/condensate field have been explored, major oil and
gas companies start to explore and develop mature and remote gas/condensate field (Knudsen
and Solvik, 2011). A big portion of remaining hydrocarbon reserve is located in remote deep
water area and harsh condition (Hjelmeland et al., 2011). Installing floating production
facilities on those areas is a big challenge. In case of mature gas field, the pressure has been
declined yielding lower production rate and increasing liquid content in the flow. Subsea wet
gas compression without preliminary separation process give a promising solution to those
problems by increasing the pressure and to avoid installing floating production facilities in
harsh condition.
Subsea wet gas compression gives many benefits such as increased recovery factor,
accelerated production, reduced CAPEX & OPEX, operation flexibility and HSE benefit
(Vinterst, 2014). Figure 1.1 shows the effect of installing subsea wet gas compressor on
field production. By installing compressor, production rate will increase and field life time
will be extended. In some case, CAPEX & OPEX is reduced because subsea wet gas
compressor eliminate the need of implementing expensive floating production facilities in
harsh offshore condition (Hjelmeland et al., 2011).
10
Figure 1.1 Effect of installing subsea wet gas compression to gas/condensate field performance
(Courtesy of Vinterst, 2014).
To give satisfying forecast of a gas condensate field project installed with subsea wet
gas compressor, a good compressor performance model for wet gas fluid should be used.
Currently no standard exist for wet gas compressor performance (Hundseid et al., 2008). For
dry gas, Schultz method (Schultz, 1962) gives satisfying result and already used in a standard
ISO 5389:2005. But for wet gas compression, Schultz method may not be accurate. It is
because in a wet gas fluid, the composition on each phase is changing from inlet condition
into outlet condition of the compressor as the pressure increase. Therefore, fluid properties
changes as gas compressed. Schultz made a simplification in his method by using average
gas properties at exit and outlet compressor. However, this assumption may have a significant
impact on the Schultz compressor performance that reduces accuracy of compressor
performance prediction (Hundseid et al, 2006).
The necessities of standard in wet gas compression encourage this project to conduct
a research in wet gas compression performance. This project will evaluate Huntington direct
integration method (Huntington, 1985) which originally for dry gas to be used in wet gas
compression calculation. Direct integration method predicted to be suitable for wet gas
compression since fluid properties change along compression is taken into account.
11
2. THEORY
Compressors are a device to pass work to the fluid in order to change the state of the fluid,
generally to increase pressure and/or elevation (Moran et al., 2011). Fluid goes into the
compressor at low pressure. Inside the compressor, energy converted from electricity into
pressure energy. This conversion process in general convert electricity into kinetic energy to
a part of the compressor. This kinetic energy then transferred into the fluid and after several
processes, it produces higher pressure at outlet compressor. In depth how the process
occurred, depend on type of compressor.
12
Choosing the type of compressor to be used is depend on several factors such as flow
rate, pressure, cost and gas characteristic (Forsthoffer, 2005). Figure 2.2 show the range of
pressure and flow rate range for each compressor type. Beside flow range and pressure
criteria, cost also becomes an important factor. In general, dynamic compressor (centrifugal
and axial compressor) usually become the first choices since it have low maintenance cost.
In oil and gas industries, the centrifugal compressor type is the most dominant type of
compressor that being used. It is because centrifugal compressor have high efficiency, can
handle high flow capacity and robust (Brenne et al, 2008; Guo et al, 2007).
13
Figure 2.3 Cross section of centrifugal compressor (Courtesy of petroleum extension service)
Equation 1:
Graphical representation of work described by Equation 1 can be seen in Figure 2.4. In the
figure, line A and B represent state change during compression process (for compression, the
arrow should be reversed). The initial and end point are the same for both lines. From
Equation 1, work done on compression process is area under curve pressure-volume which is
area under line A or B in Figure 2.4. It can be seen that work done in a compression process
depends on state change during entire compression process. It is not only depends on initial
and final state. Thus to find work in a compression process, the relation between pressure and
volume have to be known for the entire compression process.
14
Figure 2.4 Work in graphical representation (Courtesy of Moran et al., 2011)
15
Figure 2.5 Compressor reference processes (Courtesy of ISO 5389, 2005)
Equation 2:
By using Equation 2 and Equation 1, specific head for isothermal compression is (ISO
5389, 2005):
Equation 3:
( )
Mean value of z between initial and final condition should be used in Equation 3:
Equation 4:
16
2.2.1. Isentropic compression
Isentropic is an adiabatic reversible (no heat transfer, no friction loss) process (Moran et al.,
2011). Pressure-volume relationships for isentropic processes are:
Equation 5:
Using real gas equation, the relation between pressure and temperature become:
Equation 6:
( )
From Equation 5, Equation 6, and Equation 1 specific compression head for real gas
isentropic compression become (ISO 5389, 2005)
Equation 7:
[( ) ]
Mean value of k between initial and final condition in Equation 7 should use
Equation 8:
Figure 2.6 Show isentropic and actual compression process in a Mollier diagram
(Moran et al., 2011). It is shown that actual compression processes requires more energy
compared to isentropic process. The relation between actual compression and isentropic
compression related by isentropic efficiency:
Equation 9:
17
Figure 2.6 Comparison between isentropic and actual compression (Courtesy of Moran, 2011)
Equation 10:
In polytropic compression, the ratios of differential enthalpy to work for the entire
compression are constant and it is the same for all compression ratios (ISO 5389, 2005):
Equation 11:
[( ) ]
From Equation 10, relation between pressure and temperature for real gas are:
Equation 12:
( )
The polytropic exponent n can be defined from (NEL Fluid Rep, 1959):
18
Equation 13
( ) ( )
Equation 14
isothermal condition, it assumes constant temperature. This will cause inaccuracy since gas
temperature increase as pressure rise (except for compressor for cooler).
Even though polytropic reference process calculation already take into account effect
of different operating condition, it is still not good enough to predict compressor
performance. It is because Equation 11 was derived by assuming constant polytropic
exponent along the compression path. It is not true since polytropic exponent is a function of
polytropic efficiency and isentropic exponent, which is changing along compression path
(Hundseid et al., 2006).
To resolve the changing gas property along compression path, several authors propose
solution. There are three methods that widely used in industry. They are:
19
2.3.1. Schultz method
To settle the issue of changing gas property along the path, Schultz propose new solution by
introducing new variable (Schultz, 1962). The new variable will take consideration the effect
of changing gas property along the path. In his solution, he still use polytropic compression
path, but add f factor (Schultz correction factor) into Equation 11. Schultz equation is used in
ISO 5389:2005 standard (ISO 5389, 2005). The modified equation is:
Equation 15:
[( ) ]
Equation 16:
nv is given by
Equation 17:
( ) ( )
Equation 18:
( )
Equation 19:
( )
And
Equation 20:
20
Or
Equation 21:
Equation 22:
Equation 23:
Equation 24:
Use polytropic path equation (Equation 25) when deriving his method
Equation 25:
Equation 26:
( ) ( )
Huntington starts the derivation from Equation 25, Equation 26, and Maxwell
Equation 27 below to calculate polytropic efficiency (Equation 28).
21
Equation 27:
Equation 28:
( ) [( ) ] [ ( )]
The constant a, b, and c in Equation 26 are calculated from initial point, end point,
and intermediate point along the compression path. The intermediate pressure given by
Equation 29:
Equation 30:
Equation 31:
Equation 32:
[( ) ]
Equation 33:
( )
( )
Equation 34:
The error of this method can be further reduced by improve the midpoint temperature
estimate as follow
22
Equation 35:
[ ]
Equation 36:
( ) [( ) ] [ ( )]
{ }
( ) [( ) ] [ ( )]
Huntington found out that iteration of Tm once is enough to reduce variation of Tm to less
than 0.1 K.
In the same paper as section 2.3.2, Huntington calculates compression process by numerical
integration (Huntington, 1985). He divides the polytropic compression path into series of
small pressure steps. If the pressure step divided sufficiently small enough, compression work
at each step can be calculated by
Equation 37:
And
Equation 38:
Work/head calculation on each part calculated by trial and error. For the case where
initial pressure and temperature known, but only final pressure known an estimate of sub step
final temperature must be assumed. Head calculated by Equation 37 compared to head
calculated by Equation 38. The sub step final temperature is altered until Equation 37 equal to
Equation 38 (Huntington, 1985).
Total work calculated by this method is sum of each sub step work:
Equation 39:
23
2.4. Compressor curve map
Manufacturer of compressor give their compressor performance in a curve map. Figure 2.7
show typical compressor curve map of compressor. The abscissa is normalized volume flow
and the ordinate is head, discharge pressure, or pressure ratio. The solid lines represent
performance at several constant compressor speeds. The dashed curve in the middle is
constant efficiency curve. There is an optimum operating point as shown in the graph. The
operating point of the compressor must between surge limit and stonewall limit to avoid
damage or to avoid reduced efficiency (Hanlon, 2001).
Stonewall in Figure 2.7 is a condition which increasing capacity will give rapid
decrease in head as flow increased. This is occurred because fluid velocity approaching sound
velocity (M=1.0). When compressor operates at rate higher than stonewall limit, it will
reduce compressor efficiency and may damage the compressor (Hanlon, 2001).
Surge limit is peak head of compressor curve (Hanlon, 2001). If flow rate reduced
below surge limit, head will decrease. Surge will damage compressor and should be avoided.
When surge take place, nearly all of compressor component will reverse bended (Hanlon,
2001).
24
2.5. Effect of wet gas
Wet gas flow is defined by fluid with high gas volume fraction with liquid volume fraction up
to 5% (Hundseid et al., 2008). Wet gas fluid flow in gas condensate field usually flows as an
annular flow. Figure 2.8 show annular flow in horizontal pipe. Annular flow characterized by
very high gas velocity that moves in the center of pipe. Liquid will flow as a thin film
covering the pipe and as a droplet in the gas.
The presence of liquid alters compressor performance by momentum, heat, and mass
transfer between phases (Hundseid et al., 2008). In more detail, there are several phenomena
that will occur in wet gas compressor such as evaporate cooling, heat transfer, liquid
entrainment/deposition, and liquid film.
Two phenomena that will reduce temperature of the compressed fluid are evaporate
cooling and heat transfer (Hundseid et al., 2008). Evaporate cooling is a process of
evaporation of liquid into gas by heat transfer from gas phases into liquid phases. The phase
change event will result internal loss in the fluid and increase in entropy. The cooling effect
will take place even no phase change event occurred. This is because of high heat capacity of
the liquid phase.
Liquid entrainment is a situation that will reduce the total kinetic energy of the fluid
(Hundseid et al., 2008). The kinetic energy of the fluid will reduced if entrained liquid is
accelerated. And also, total kinetic energy is reduced if the droplets are move into the liquid
film.
Liquid film in the compression process will increase frictional loss of the system
(Hundseid et al., 2008). It can be explained as follow. The liquid film will cover the surface
of impeller. This will increase the surface roughness due to liquid wave and droplet impact.
Another effect is the liquid film act as a blockage to the flow area. Both effects will further
increase friction loss in the system.
25
Figure 2.9 show effect of gas volume fraction to compressor performance (Hundseid
et al., 2008). The data obtained from experiment data using closed loop system at
StatoilHydro test facility. From the figures it can be seen that by increasing liquid content in
the flow, the specific polytropic head is reduced.
Figure 2.9 Effect of wet gas presence to specific polytropic head (Courtesy of Hundseid et al,
2008)
In the paper, Hundseid et al describe that the reduced specific polytropic head is
caused by the effect of high liquid density (mass flow), high flow rate, and compressibility
effect (Hundseid et al., 2008). Their calculation is based on direct integration method in
section 2.3.3. From Equation 37 it is clear that by increasing density, polytropic head will be
reduced.
26
Figure 2.10 Mass flow corrected specific polytropic head (Courtesy of Hundseid et al, 2008)
Figure 2.11 Wood s corrected specific polytropic head (Courtesy of Hundseid et al, 2008)
Figure 2.12 show the effect of different fluid density into compressor curve (Hanlon,
2001). From the figure it can be seen that by increasing density, the head curve will increase.
27
Figure 2.12 Effect of fluid density to compressor curve (Courtesy of Hanlon, 2001)
28
2.7.1. Dimensional analysis method
The corrections of compressor performance are based on dimensional analysis. The following
equations are used to correct compressor parameter (Boyce, 2003):
29
Equation 45 (corrected pressure):
Ma et al (2013) propose a new procedure for off design compressor operation by iterative
method. Figure 2.13 show general correction procedure by iterative method. First task is to
transform operating inlet condition into test reference condition. After correcting operating
inlet condition into reference test condition, the discharge condition is calculated iteratively
based on corrected inlet condition. The final procedures are to recalculate polytropic head,
power, and efficiency based on corrected inlet condition.
Figure 2.13 Basic procedure for off design correction by iterative method (Courtesy of Ma et al,
2013)
30
In the paper, they claimed that this method eliminate some limitation in ASME-PTC-
10(ASME, 1998) standard. The limitations of the standard are it does not account the change
in pressure and temperature change in the inlet, thus it cannot be applied in general case.
Another limitation of the standard is it requires complex input data to do the correction.
Equation 47:
Equation 48:
Equation 49:
Equation 50:
Equation 51:
Equation 52:
Equation 53:
3. Calculate nc initial
31
Equation 54:
[ ]
4. Set nc = nc initial
Equation 55:
Equation 56:
( )
Equation 57:
( )
Equation 58:
( ) ( )
Equation 59:
9. Repeat step 5-8 100 with new nc calculated from Equation 59.
32
10. Calculate Schultz correction factor f by
Equation 60:
Where
Equation 61:
Equation 62:
12. Correct speed and mass flow using fan law (Godse, 2006).
Equation 63:
Equation 64:
Equation 65:
33
3. OBJECTIVE AND METHODOLOGY
Table 1 show fluid properties data used as input in this project. The fluid is a typical wet
gas/gas condensate hydrocarbon fluid (McCain, 1990). The table consists of separator gas
compositions and separator liquid compositions. Table 2 show the property of C 7+ component
Fluid compositions
Mole fraction Mole fraction
component
separator gas separator liquid
C1 0.837 0.0869
C2 0.096 0.0569
C3 0.046 0.0896
i-C4 0.006 0.0276
n-C4 0.009 0.0539
i-C5 0.003 0.0402
n-C5 0.002 0.0400
n-C6 0.001 0.0782
C7+ 0.000 0.5267
Mw total 19.68 83.84
34
Table 2 C7+ properties (Courtesy of McCain, 1990)
Figure 3.1 Process flow diagram of mixing liquid and gas phase
Cases 1 2 3 4
Gas volume fraction (at inlet
1 0.99 0.97 0.96
condition Case A and B)
35
3.1.2. Calculation cases and workflow
Figure 3.2 shows general flowchart for calculating compressor performance. In summary, the
first thing to do is to check whether compressor test design condition differ from operating
condition. If it is different, then the compressor curve map needs to be corrected by
dimensional analysis method or iterative method. If no correction needed, the calculation
goes directly into compressor performance calculation by Schultz method (Schultz, 1962) or
direct integration method (Huntington, 1985). The final outputs are compressor performance,
outlet condition, or head.
36
In this project, two cases of compressor calculation were created (Figure 3.2). The
first case (case A) is compressor performance calculation without compressor curve map, and
second case (Case B) is case with compressor curve map. On both cases, fluid properties
were varied according to Table 3 and produce fluid variation case in Figure 3.3.
Table 4 gives input parameter such as inlet condition, compressor efficiency, and
outlet pressure for case A. In this case, no compressor curve map was given and compressor
test reference condition is assumed to be the same with actual condition. Thus, no off design
correction needs to be done.
For case B, the inlet condition is given by Table 5. Table 6 and Figure 3.4 is the
compressor curve map tested at condition described in Table 7. When performing compressor
curve map, dry gas composition in Table 3 was used as test fluid. As can be seen from Table
5 and Table 7, the reference test condition and actual condition is different. Thus, off design
correction need to be done in this case.
37
Table 5 Case B inlet condition
Inlet condition
P1 500.0 psia
T1 563.0 R
molarflow 1500 kgmole/hr
38
Table 7 Compressor test reference
The calculation for off design correction is using iterative method described in section
2.7.2. The off design result from iterative method will be compared with dimensional analysis
method described in section 2.7.1.
For compressor performance calculation, both Schultz method and direct integration
method will be calculated. The calculation procedures have been described in section 2.3.1
and 2.3.3. The direct integration method needs to be programmed to facilitate large step size.
ExcelTM visual basic program for direct integration method is on Appendix. Result from
direct integration method will be compared by Schultz method calculated by HYSYSTM.
It is desired to know how HYSYSTM calculate wet gas compressor performance. One
hypothesis is that the wet gas compression can be represented by compression on of each
phase separately as seen in Figure 3.5. In the figure, the wet gas fluid will be separated before
being compressed. Vapor phase goes into compressor while liquid phase goes into pump.
After both fluids being compressed, the fluids will be mixed again. The result of this method
will be compared with normal compression without separator as seen in Figure 3.6. Schultz
method will be used to calculate both compressor performances.
39
Figure 3.6 Compressor performance calculation without separator
40
4. RESULT
As described in previous section, two cases was created (Case A and B) in this project. On
each case, the fluid gas volume fraction properties were varied according to Table 3. The
difference between case A and case B is the data input and off design correction. In case A,
no compressor curve map was given and no off design correction. For case B, the input data
are compressor curve and off design correction need to be done in case B. In both case, the
result of gas compression using prior separator and no separator will also be compared.
4.1. Case A
Figure 4.1 show compression path graph (pressure-molar volume diagram) for case A with
various gas volume fraction (GVF=1, GVF=0.98, GVF=0.97, GVF=0.96). Fluid compressed
from 700 psia to 1000 psia. The pressure-volume data was gathered using direct integration
method at 1, 10, and 60 total pressure step size. The direct integration method compression
path compared with Schultz method in the graph.
41
Figure 4.2 shows fluid properties (gas volume fraction) change on each pressure step
as the fluid compressed from initial pressure into final pressure for case A. Each point in the
graph calculated using direct integration method with total pressure step size equal to 10.
Four initial GVF (GVF=1, GVF= 0.99, GVF=0.97, and GVF=0.96) cases also shown in the
figure.
Figure 4.2 Gas volume fraction change at each pressure step i with total number of pressure
step = 10 for case A
Figure 4.3 through Figure 4.6 Show polytropic exponent change as the pressure step
increase toward final pressure for case A. The graph show direct integration calculation data
for various gas volume fractions (GVF) at total pressure step size of 1, 4, 10, 18, 30, and 60.
The figures also show Schultz polytropic exponent parameter as a line (Schultz polytropic
exponent is not a function of pressure step size, but it is shown in the figure for comparison
purpose).
42
Figure 4.3 Polytropic exponent change at each pressure step i, case A GVF = 1
Figure 4.4 Polytropic exponent change at each pressure step i, case A GVF = 0.99
43
Figure 4.5 Polytropic exponent change at each pressure step i, case A GVF = 0.97
Figure 4.6 Polytropic exponent change at each pressure step i, case A GVF = 0.96
44
Direct integration method calculation requires Equation 37 and Equation 38 equal.
The difference between those equation calculated by Equation 66 below. The calculated
errors and total specific head for case A shown in Figure 4.7 through Figure 4.10 below.
Equation 66:
Figure 4.7 Calculated head and error vs total pressure step size, case A GVF = 1
45
Figure 4.8 Calculated head and error vs total pressure step size, case A GVF = 0.99
Figure 4.9 Calculated head and error vs total pressure step size, case A GVF = 0.97
46
Figure 4.10 Calculated head and error vs total pressure step size, case A GVF = 0.96
Figure 4.11 through Figure 4.14 show the result of calculated total specific head and
power using Schultz and direct integration method with various gas volume fractions (GVF)
for case A. The abscissa represents total pressure step size in direct integration method. The
ordinate show total specific head and power required to compress the gas.
47
Figure 4.11 Calculated head and power vs total pressure step size, case A GVF = 1
Figure 4.12 Calculated head and power vs total pressure step size, case A GVF = 0.99
48
Figure 4.13 Calculated head and power vs total pressure step size, case A GVF = 0.97
Figure 4.14 Calculated head and power vs total pressure step size, case A GVF = 0.96
49
Figure 4.15 through Figure 4.18 show percentage of calculated specific head at each
pressure step relative to total specific head for case A. The percentages are calculated by
Equation 67 below
Equation 67:
Specific head on each pressure step are calculated by Equation 38. Data from various
gas volume fraction and several total pressure step size (6, 8, and 10) also shown in the
figures.
50
Figure 4.16 Percentage of head on each pressure step i, case A GVF = 0.99
Figure 4.17 Percentage of head on each pressure step i, case A GVF = 0.97
51
Figure 4.18 Percentage of head on each pressure step i, case A GVF = 0.96
4.2. Case B
In case B, compressor curve map is an input and actual condition are different from test
reference condition (see Table 5 and Table 7). The compressor curve map needs to be
corrected first before used as calculation. Figure 4.19 through Figure 4.22 show the original
compressor curve map and corrected compressor curve map for various GVF (GVF=1,
GVF=0.99, GVF=0.97, and GVF=0.96). Compressor curve map corrected using dimensional
analysis method and iterative method described in section 2.7.1 and 2.7.2.
52
Figure 4.19 Original vs corrected compressor curve, case B GVF = 1
53
Figure 4.21 Original vs corrected compressor curve, case B GVF = 0.97
54
Figure 4.23 shows compression process path calculated using corrected compressor
curve with direct integration method. Fluid compressed from initial pressure of 500 psia into
final pressure. In the figure, compression processes of four different fluid properties (gas
volume fraction) are shown. Data shown in the figure are for direct integration method with
total pressure step size of 1, 10, and 25. Schultz compression path also shown in the figure.
Figure 4.24 shows fluid properties (gas volume fraction) change on each pressure step
as the fluid compressed from initial pressure into final pressure for case B. Data in the graph
calculated by direct integration method with total pressure step size of 10. Four initial gas
volume fraction (GVF=1, GVF= 0.99, GVF=0.97, and GVF=0.96) cases also shown in the
figure.
55
Figure 4.24 Gas volume fraction change at each pressure step i with total number of pressure
step = 10 for case B
Figure 4.25 through Figure 4.28 show polytropic exponent on each pressure step size
for case B. The data show calculation result using direct integration method at total pressure
step size of 1, 4, 10, and 18. Four different fluid properties (GVF=1, GVF=0.99, GVF=0.97,
GVF=0.96) calculation result also shown in the graphs. In addition, Schultz polytropic
exponent also show in the figure.
56
Figure 4.25 Polytropic exponent change at each pressure step i, case B GVF = 1
Figure 4.26 Polytropic exponent change at each pressure step i, case B GVF = 0.99
57
Figure 4.27 Polytropic exponent change at each pressure step i, case B GVF = 0.97
Figure 4.28 Polytropic exponent change at each pressure step i, case B GVF = 0.96
58
Figure 4.29 through Figure 4.32 show calculated error using Equation 66 with various
initial gas volume fraction (GVF=1, GVF=0.99, GVF=0.97, GVF=0.96). The abscissa in the
figure is total pressure step size. The ordinate show calculated error and total specific head
calculated by direct integration method.
Figure 4.29 Calculated head and error vs total pressure step size, case B GVF = 1
59
Figure 4.30 Calculated head and error vs total pressure step size, case B GVF = 0.99
Figure 4.31 Calculated head and error vs total pressure step size, case B GVF = 0.97
60
Figure 4.32 Calculated head and error vs total pressure step size, case B GVF = 0.96
Figure 4.33 through Figure 4.34 shows calculated total specific head and power using
direct integration method and Schultz method for case B. Calculated data using four initial
gas volume fractions are also shown in the figures. The abscissa show total pressure step size
in direct integration method. The ordinate show total specific head and power required to
compress the fluid. Calculated compressor performance using Schultz method shown as a
line in the figure as comparison purpose only. Actually, Schultz method is not a function of
total pressure step size.
61
Figure 4.33 Calculated head and power vs total pressure step size, case B GVF = 1
Figure 4.34 Calculated head and power vs total pressure step size, case B GVF = 0.99
62
Figure 4.35 Calculated head and power vs total pressure step size, case B GVF = 0.97
63
Figure 4.36 Calculated head and power vs total pressure step size, case B GVF = 0.96
Figure 4.37 through Figure 4.40 show head percentage on each step calculated by
Equation 67 for various gas volume fraction (GVF=1, GVF=0.99, GVF=0.97, GVF=0.96).
Data shown in the figure are calculated from direct integration method with total pressure
step size of 6, 8, and 10.
64
Figure 4.37 Percentage of head on each pressure step i, case B GVF = 1
Figure 4.38 Percentage of head on each pressure step i, case B GVF = 0.99
65
Figure 4.39 Percentage of head on each pressure step i, case B GVF = 0.97
Figure 4.40 Percentage of head on each pressure step i, case B GVF = 0.96
66
Table 8 and Table 9 show summary of compressor performance for case A and case
B. Direct integration method data are shown using total pressure step size of 10. The head
and power difference between Schultz method and direct integration method calculated by
using the following equations Equation 68 and Equation 69
Equation 68:
Equation 69:
Total Case A
pressure
step size
Parameter Method
(For direct GVF=1 GVF=0.99 GVF=0.97 GVF=0.96
integration
method)
Direct
Total specific head (Kj/Kg) integration 10 43.3 35.3 30.0 26.1
Schultz - 43.2 35.2 30.1 26.1
Direct
Power (kW) integration 10 157.5 148.3 141.6 135.7
Schultz - 157.5 148.2 141.7 135.5
specific head
Schultz-direct - -
integration (%) -0.02% -0.19% 0.28% -0.20%
power Schultz-direct
- -
integration (%) -0.01% -0.10% 0.08% -0.13%
Total Case B
pressure
step size
Parameter Method
(For direct GVF=1 GVF=0.99 GVF=0.97 GVF=0.96
integration
method)
Direct
Total specific head (Kj/Kg)
integration 10 76.9 81.9 105.0 127.5
67
Schultz - 72.8 82.8 112.8 143.6
Direct
Power (kW) integration 10 902.4 1184.0 1920.7 2789.4
Schultz - 817.1 1081.9 1725.8 2587.3
specific head
Schultz-direct - -
integration (%) -5.24% 1.08% 7.46% 12.63%
power Schultz-direct
- -
integration (%) -9.45% -8.63% -10.15% -7.24%
Table 10 and Table 11 show result of compression of wet gas using separator and
without separator as in Figure 3.5 and Figure 3.6. For case A, the difference between wet gas
compression with separator and without separator are calculated by Equation 70 below
Equation 70:
Table 10 Comparison of wet gas compression using separator and without separator (case A)
CASE A
Table 11 Comparison of wet gas compression using separator and without separator (case B)
CASE B
68
5. DISCUSSION
Figure 4.19 through Figure 4.22 show that as gas volume fraction decrease, corrected head
for case B increased and the deviation between dimensional analysis methods with iterative
method increasing. As gas volume fraction decrease, molecular weight increase (more liquid
present in the mixture). Increasing molecular weight relative to test reference condition will
lead to increasing head, which is consistent with Figure 2.12. Correction method using
dimensional analysis and iterative method give different result as molecular weight increase.
This is because of fluid in test reference condition is using dry gas composition. As gas
molecular weight increase, the molecular weight difference between actual and test reference
condition become larger. And both methods may not be accurate for large deviation.
From pressure-molar volume graph in Figure 4.1 and Figure 4.23 (case A and B) it
can be seen that as gas volume fraction decrease, the compression path move to the left. This
is because the presence of hydrocarbon liquid which is heavier (high density) than gas at the
same pressure and temperature.
Figure 4.1 also show that there is no visible difference of compression path in case A
using different total pressure step size (1, 10, and 60) and Schultz method.
For case B in Figure 4.23, final compressor outlet pressure increases with decreasing
gas volume fraction. This is because of in case B, corrected curve head is higher at low gas
volume fraction. Thus, it will give higher compressor discharge pressure. The presence of
liquid also increase discharge pressure.
As can be seen in Figure 4.2 and Figure 4.24 (case A and B), gas volume fraction
(GVF) is decreasing with increasing pressure step size at initial gas volume fraction less than
1. This event can be described in Figure 5.1 and Figure 5.2. The figures show compression
process path in a pressure-volume diagram with bubble point and dew point line. It can be
seen that as the pressure increase, fluid states move closer to the bubble point line. Therefore,
an increase in liquid content is expected.
69
Figure 5.1 Compression process path case A with bubble point and dew point line
Figure 5.2 Compression process path case B with bubble point and dew point line
Figure 4.2 and Figure 4.24 (case A and B) also shows that the change in gas volume
fraction is bigger at lower GVF. At lower gas volume fraction, polytropic exponent values
become smaller (See Figure 4.3 through Figure 4.6 for case A and Figure 4.25 through Figure
70
4.28 for case B). Lower polytropic exponent implies higher molar volume change as pressure
change (See Figure 5.3). Accordingly, the fluid state moves faster into bubble point line
(lower gas volume fraction/higher liquid content).
Figure 5.3 Polytropic exponent in a p-v diagram (Courtesy of Moran et al., 2010)
Figure 4.3 through Figure 4.6(case A) and Figure 4.25 through Figure 4.28(case B)
show that polytopic exponent increase as pressure step increases. The increase of polytropic
exponent was caused by changing property of fluid as pressure and temperature increase. As
discussed before, more liquid will form at higher pressure. Properties of gas and liquid are
different, thus it will also contribute to the change of polytropic exponent by altering mixture
compressibility.
From Figure 4.7 through Figure 4.10(case A) and Figure 4.29 through Figure 4.32
(case B) show that there is an optimum number of total pressure step size in direct integration
method. From section 2.3.3 it is described that Equation 37 and Equation 38 must equal, the
difference between the equations was calculated from Equation 66. At low total number of
pressure step size, the error calculated by Equation 66 is large (up to 1.5%). It is because the
small pressure step size cannot take into account change in fluid properties. But at high
71
number of pressure step size, the calculated head become divergence because of numerical
error (for case A only). From the figures, the optimum number of total pressure step size is
equal to 10 (calculated error approaching 0%).
Table 8 show that calculated head decreased with increasing gas volume fraction for
case A. As described in section 2.5, factors that causing decrease in head with increasing
GVF is evaporate cooling and heat transfer. From Figure 4.2 it is clear that evaporate cooling
is occurred as gas compressed. The result of evaporate cooling is decreased temperature with
more liquid in fluid mixture as can be seen in Figure 5.4.
For case B in Table 9, as gas volume fraction decreased, calculated head is increasing.
This is because of in case B, compression ratio becomes higher at low gas volume fraction
(see Figure 4.23). Thus it will need more head to compress the fluid.
From Table 8, Table 9, Figure 4.11 through Figure 4.14, and Figure 4.33 through
Figure 4.36 show that Schultz method and direct integration head give different results.
Schultz method give lower head and power compared to direct integration method except for
GVF=0.97. Schultz methods give lower result because polytropic exponent in Schultz is
72
lower compared to entire direct integration methods polytropic exponent (See Figure 4.3
through Figure 4.6 for case A and Figure 4.25 through Figure 4.28 for case B). But
calculated head in Schultz method will be corrected by Schultz correction factor f.
Figure 4.15 through Figure 4.18 (case A) and Figure 4.37 through Figure 4.40 (case
B) show that as pressure step increase, calculated head will slightly increase on each pressure
step. It is caused by changing fluid properties as pressure and temperature increase. As
discussed before, as pressure step increase polytropic exponent will increase. From section
2.2 it is clear that increasing polytropic exponent will increase calculated head.
The results in Table 10 and Table 11 show that HYSYSTM does not calculate wet gas
compression where each phase compressed separately using compressor (for vapor) and
pump (for liquid) as seen in Figure 3.5. For case A, the difference between compression using
separator and without separator is up to -3%. For case B the difference is larger (up to 57%).
Therefore, HYSYSTM calculate wet gas compression by considering interaction between
liquid and vapor phase.
73
6. CONCLUSION AND RECOMMENDATION
Compressor performance running at off design condition will have large deviation
from original tested performance, therefor it needs correction. The correction can be done by
using iterative method or dimensional analysis method. Both methods give the same result if
the difference between test condition and actual condition is small. At large difference
between actual and reference condition, both method give significantly different result.
Accordingly, corrected compressor performance becomes uncertain.
The presence of liquid in compressor will reduce specific head of the compressor.
Specific head reduced because of liquid having higher density compared to gas. Another
factor that reduces compressor specific head is evaporate cooling and heat transfer effect. The
presence of liquid will cool off the fluid since it has high heat capacity.
Direct integration method is much more appropriate for application in wet gas
compression compared to Schultz. This is due to changing fluid properties as wet gas
compressed. The direct integration method calculation will take consideration of changing
fluid properties as pressure goes up. In Schultz method, fluid properties change does not
completely taken into account. It only uses average fluid properties and add Schultz
correction factor.
There is an optimum number of total pressure step size in calculation using direct
integration method. If the pressure step size is too low, compressor performance accuracy is
low because of changing fluid properties does not well represented in the calculation. If the
total pressure step size is too high, there will be numerical error. The optimum number of
total pressure step size is when head calculated by enthalpy difference and head calculated by
molar volume- pressure is nearly equal.
74
separation. The result between case with separator and without separator gives large
difference.
75
NOMENCLATURE
HP = Power, j/s
k = Isentropic exponent
MW = Molecular weight
n = Polytropic exponent
q = Heat flow, kJ
S = Entropy, j/K
T = Temperature, C [K]
V = Volume, m3
W = Work, J
76
X = Compressibility function
Y = Compressibility function
z = Compressibility factor
= Efficiency, % [fraction]
= Difference
Subscript
1 = Inlet condition
2 = Outlet condition
act = Actual
avg = Average
corr = Corrected
h = Huntington method
i = Pressure step i
m = Midpoint
p = Polytropic
P = Isobar
s = Isentropic
S = Schultz method
77
T = Isothermal
t = Total
v = Volume corrected
superscript
c = Corrected
78
REFERENCE
ASME 1998. Performance Test Code on Compressors and Exhausters, the Society.
BOYCE, M. P. 2003. Centrifugal Compressors: A Basic Guide, PennWell.
, ., , ., , . ., CH, . . , H. 200 .
Performance evaluation of a centrifugal compressor operating under wet gas
conditions. Proceedings Of The Thirty-Fourth Turbomachinery Symposium.
FORSTHOFFER, W. E. 2005. 1. Forsthoffer's Rotating Equipment Handbooks:
Fundamentals of Rotating Equipment, Elsevier Science.
GIAMPAOLO, T. 2010. Compressor Handbook: Principles and Practice, Fairmont Press.
GODSE, A. 2006. Compressors and Modern Process Applications, Wiley.
GUO, B., WILLIAM C. LYONS, P. D. P. E. & GHALAMBOR, A. 2011. Petroleum
Production Engineering, A Computer-Assisted Approach, Elsevier Science.
HANLON, P. 2001. Compressor Handbook, Mcgraw-hill.
HJELMELAND, M., OLSEN, A. B. & MARJOHAN, R. 2011. Advances in Subsea Wet Gas
Compression Technologies. International Petroleum Technology Conference.
HUNDSEID, ., BAKKEN, L. E. & HELDE, T. 2006. A REVISED COMPRESSOR
POLYTROPIC PERFORMANCE ANALYSIS. ASME Turbo Expo 2006: Power for
Land, Sea and Air.
HUNDSEID, O. Y. & BAKKEN, L. E. 2006. WET GAS PERFORMANCE ANALYSIS.
ASME Turbo Expo 2006: Power for Land, Sea and Air.
H S , . ., , . ., , . ., , . , . 2008.
Wet Gas Performance of a Single Stage Centrifugal Compressor Proceedings of
ASME Turbo Expo 2008: Power for Land, Sea and Air.
HUNTINGTON, R. A. 1985. Evaluation of Polytropic Calculation Methods for
Turbomachinery Performance. Journal of Engineering for Gas Turbines and Power,
107, 872-876.
ISO 2005. Turbocompressors Performance test code Standard ISO 5389:2005.
KNUDSEN, T. W. & SOLVIK, N. A. 2011. SWorld First Submerged Testing of Subsea Wet
Gas Compressor. Offshore Technology Conference.
MA, Y., FRETHEIM, H., PERSSON, E. & HAUGEN, T. 2013. AN ITERATIVEMETHOD
APPLIED TO CORRECT THE ACTUAL COMPRESSOR PERFORMANCE TO
THE EQUIVALENT PERFORMANCE UNDER THE SPECIFIED REFERENCE
CONDITIONS.
MCCAIN, W. D. 1990. The Properties of Petroleum Fluids, PennWell Books.
MORAN, M. J., SHAPIRO, H. N., BOETTNER, D. D. & BAILEY, M. 2010. Fundamentals
of Engineering Thermodynamics, John Wiley & Sons.
NEL 1959. NEL Fluids Report, DSIR, National Engineering Laboratory, Fluid Mechanics
Division.
SCHLUMBERGER. 2014. Annular Flow. [Accessed 20 May 2014]. Available from:
http://www.glossary.oilfield.slb.com/en/Terms/a/annular_flow.aspx
79
SCHULTZ, J. M. 1962. The Polytropic Analysis of Centrifugal Compressors. Journal of
Engineering for Gas Turbines and Power, 84, 69-82.
VINTERST, T. 2014. Subsea Gas Compression Now And In The Future. Energy Claims
Confrence 2013.
WOOD, A. B. 1941. A textbook of sound: being an account of the physics of vibrations with
special reference to recent theoretical and technical developments, The Macmillan
company.
80
APPENDIX
VISUAL BASIC SCRIPT FOR DIRECT INTEGRATION METHOD
Dim N As Integer
81
connect with H S S
hyApp.Visible = True
fileName = Worksheets("Integration").Range("Filename_int")
simCase.Visible = True
End If
End If
'Clear Content
Sheets("Integration").Select
Worksheets("Integration").Range("TopLeft_Calc_Int").Select
Range(Selection, Selection.End(xlDown)).Select
Range(Selection, Selection.End(xlToRight)).Select
Selection.ClearContents
rowcalc = Worksheets("Integration").Range("TopLeft_Calc_Int").Row
colcalc = Worksheets("Integration").Range("TopLeft_Calc_Int").Column
rowstep = Worksheets("Integration").Range("TopLeft_Step_Int").Row
colstep = Worksheets("Integration").Range("TopLeft_Step_Int").Column
82
'Topleft of result table
rowresult = Worksheets("Integration").Range("Result_Table").Row
colresult = Worksheets("Integration").Range("Result_Table").Column
Sheets("Integration").Select
Worksheets("Integration").Range("Topleft_Curve_Int_Ref").Select
rowcurve = Worksheets("Integration").Range("Topleft_Curve_Int_Ref").Row
colcurve = Worksheets("Integration").Range("Topleft_Curve_Int_Ref").Column
Sheets("Integration").Select
Worksheets("Integration").Range("Topleft_Comp_Int").Select
rowcomp = Worksheets("Integration").Range("Topleft_Comp_Int").Row
colcomp = Worksheets("Integration").Range("Topleft_Comp_Int").Column
Sheets("Integration").Select
Worksheets("Integration").Range("Topleft_Comp_Int_Ref").Select
rowcomp_ref = Worksheets("Integration").Range("Topleft_Comp_Int_Ref").Row
colcomp_ref = Worksheets("Integration").Range("Topleft_Comp_Int_Ref").Column
rowcalc_curve = Worksheets("Integration").Range("TopLeft_Calc_Curve_Int").Row
colcalc_curve = Worksheets("Integration").Range("TopLeft_Calc_Curve_Int").Column
83
'Count Column calculation table
Sheets("Integration").Select
Worksheets("Integration").Range("Header_Calc_Table").Select
number = Worksheets("Integration").Range("Number_step_Int").Value
row_curve_result_stepi = Worksheets("Integration").Range("TopLeft_Curve_Res_Int").Row
col_curve_result_stepi = Worksheets("Integration").Range("TopLeft_Curve_Res_Int").Column
'Curve_corr_Result_Int
g=1
h = col_curveRPM_result
84
'loop for different number of integration step (to the right)
For y = 1 To number
Worksheets("Integration").Range("molarflow_int").Value
=Worksheets("Integration").Range("molarflow_parameter_Int").Value
compositions = inlet.ComponentMolarFractionValue
Next item
inlet.ComponentMolarFraction.Values = compositions
compositions = inlet.ComponentMolarFractionValue
Next item
inlet.ComponentMolarFraction.Values = compositions
85
Worksheets("Integration").Range("MW_INT").Value = inlet.MolecularWeight.GetValue("")
For z = 1 To step
Sheets("Integration").Select
Worksheets("Integration").Range("TopLeft_Calc_Int").Select
Worksheets("Integration").Cells(rowcalc, k + 1).Value =
Worksheets("Integration").Range("Press_init_Int").Value
Worksheets("Integration").Cells(rowcalc, k + 2).Value =
Worksheets("Integration").Range("Temp_init_int").Value
86
Worksheets("Integration").Cells(rowcalc + (z), k + 1).Value =
Worksheets("Integration").Cells(rowcalc + (z - 1), k + 1).Value *
Worksheets("Integration").Cells(rowstep + 1, colstep + (y - 1)).Value
molarflow = Worksheets("Integration").Range("molarflow_int").Value
87
'get inlet molar entropy
efficiency = 100
compressor.CompPolytropicEff.SetValue efficiency
88
'get outlet isentropic molar enthalpy
compressor.CompPolytropicEff.SetValue efficiency
89
Worksheets("Integration").Cells(rowcalc + (z - 1), k + 14).Value =
compressor.CompAdiabaticEff.GetValue("")
'calculate dh isentropic
'calculate dh
'polytropic head
90
Worksheets("Integration").Cells(rowcalc + (z - 1), k + 19).Value =
Worksheets("Integration").Cells(rowcalc + (z - 1), k + 18).Value * 1000 /
Worksheets("Integration").Range("gravity").Value
'set outlet molar enthalpy (find temperature, known pressure and enthalpy)
91
'get molar density
'Calculate Power
'calculate average v
92
Worksheets("Integration").Cells(rowcalc + (z - 1), k + 20).Value = DeltaP * densityaverage /
1000
Next z
sum_head = 0
sum_power = 0
vdeltap = 0
i_sum = 1
For i_sum = 1 To z
Next i_sum
93
Worksheets("Integration").Cells(rowresult + 1, colresult + (y - 1)) = sum_power
Next y
End Sub
94