CALIB053 - Mathcad 15 Equations
CALIB053 - Mathcad 15 Equations
CALIB053 - Mathcad 15 Equations
Flow Calculations to ISO 5167 Measurement of fluid flow by means of pressure differential devices
Note:- This file carries out the core calculation in strict SI units (Data is entered in 'Oilfield metric')
Calculation Title
OPTIONS
Wet Gas correction Method
MC 4
Gas CV
GasCV 2
1 = Mass, 2 = Volumetric
Condensate CV
CondCV 2
1 = Mass, 2 = Volumetric
Methanol Density
RhoCH3OH 2
1 = Keypad, 2 = Calculated
GDCoeff 1
1 = Fixed, 2 = Table
n 8
Temperature Input
TempInput 1
1 = Upstream, 2= Downstream
Texp 1
Water Density
RhoH2O 2
1 = Keypad, 2 = Calculated
CONFIGURATION PARAMETERS
System Parameters
tcal 20
Dcal.mm 229.92
ecs 0.000013
Thermal Expansion
dcal.mm 161.13
ess 0.000013
Thermal Expansion
Additional wet gas calculation factors
or
or
D 0.0001732
15w.in 1000
G 0.0000713
15m.in 633.6
K0 1489.067
15c 633.6
Pe 1.01325 10
g 9.8066
Pref 1.01325 10
tref 15
E 0.1307
F 1040
H 0.3344
J 760
K1 0
K2 0.0018684
A 0.59
Printed at 21/07/2016
B 0.86
C 0.436
INPUT PARAMETERS
Gas discharge coefficient
CG.in 0.99735
CL 1.2
Liquid Properties
Condensate to gas mass ratio
c 0.0110489919
w 0.0110489919
m 0.00526810719
CVg 70
CVc 70
w.in 1000
m.in 633.6
Gas Properties
Gas Isentropic exponent
1.30873
cp 0.012132
Temperature Exponent
KTe 0.2359
b 0.774972911
Process Conditions
Upstream Pressure (bar.g)
pop 25
tf.in 18
pmbar 721
f 20.6810271
DCLTable
Reynolds
Number at D
Discharge
Coefficient C
92972.00000
0.99751
8004183.95420
0.99617
12613500.71600
1.00018
20938522.39857
0.99581
21084325.29000
0.99523
25914898.81200
1.00007
30746226.64400
1.00187
44714447.81400
0.99931
48714447.81400
0.99800
10
50714447.81400
0.99850
CALIB053
Page 2 of 8
Printed at 21/07/2016
1.002
1
DischargeCoeff 0.998
0.996
0.994
0
10000000
20000000
30000000
40000000
50000000
Reynolds
Convert Dimensions to SI
Dcal
dcal
D cal.mm
0.22992
1000
d cal.mm
1000
0.16113
cp 10
Convert mbar to Pa
Beta Ratio
D cal
metre
Pf pf 1.01325 10
dcal
metre
pf p op 10
PaA
0.000012132
0.700808977035
A B C dp 8873.554713081878
(Pascal)
(Pascal)
(Pascal)
Pf 2601325
Upstream Temperature
Temperature Exponent
Te
if Texp = 2
0.2359
KTe otherwise
Tf
Te
P3
Tf.in P
if TempInput = 2
291.15
(K)
tf Tf 273.15 18
(degC)
Tf.in otherwise
Throat Diameter Temperature Correction
dop d cal 1 ess tf tcal 0.16112581062
(metre)
(metre)
Liquid to Gas Mass Ratio (from liquid component to gas mass ratios)
L c w m 0.02736609099
Gas Mass Fraction
CALIB053
Page 3 of 8
1
1 L
0.973362863316
Printed at 21/07/2016
c c g 0.010754678393
K0 K1 15c
2
K2 0.001840831181
t tf 15 3
(deg C)
15c
Fc 10
Cpl
870960
2
15c
4209.2 tf
2
15c
1.005274634893
0.000209878364
(/bar)
c C tl C pl 15c 633.418740831676
1 Pf Pe Fc 10
(kg/m3)
D tf2 E tf F
if RhoH2O = 2
w 1037.5912832
(kg/m3)
if RhoCH3OH = 2
m 754.0039012
(kg/m3)
15w 1038.00053
(kg/m3)
15m 755.0000425
(kg/m3)
w.in otherwise
Methanol Line Density
m
G tf2 H tf J
m.in otherwise
( 225 G 15 H J) if RhoCH3OH = 2
15m.in otherwise
L
c
c
w
w
780.128308508431
(kg/m3)
m
m
Pf dp
Pf dp
Pf
1
Pf
2
1
Pf dp
1 4 Pf dp
Pf
Pf
0.5
0.977494961548354
CALIB053
Page 4 of 8
Printed at 21/07/2016
Cd 1
for k 1 20
Re
Cd
Cd
qmutable
dop
2 dp f
4 qmutable
Dop
DCLTable1 2 if Re DCLTable1 1
1
1
2
1
DCLTable DCLTable 1 1 if Re max DCLTable
lookup max DCLTable
1
2
linterp submatrix( DCLTable 1 n 1 2) submatrix( DCLTable 1 n 1 2 ) Re otherwise
kk1
qmutable
Re Cd
qmutable
Re C d qmcalc
qmutable
Re C d
39.40889 17988978.67585
CG
Cd if GDCoeff = 2
3
0.99736
CG 0.99735
CG.in otherwise
Uncorrected Mixture Mass Flow Rate
qmu
qmu 39.408566841053
qmutable if GDCoeff = 2
d 2 C
op
G
4
1
1
(kg/s)
2 dp f otherwise
Re 17988978.675846118
17988829.242673423
otherwise
0.746 Frg
Define Sub-Routine
CALIB053
Page 5 of 8
DeLew
qgov.lew
Ylew
Printed at 21/07/2016
q mu
f
1 g Cd
f
if MC = 5
g
L
CL
1 g f
otherwise
g L
for run 1 15
lew qgov.lew
2
D
op
lew
f
Fr
g Dop L f
Klew
L frd f frd
f L
Flew
qmu
qm.lew
Flew
qgov.lew
Fr
frd
F
lew
Ylew
Fr
frd
F
lew
Ylew
q m.lew
f
lew
Klew
Fr
frd
F
lew
lew DeLew
Klew
Ylew
lew
Klew
4.949569680023
0.004455698537
0.590902925176
45.037281362282
1.019120579237
8.660123985453
Flew 1.019120579237104
Murdock (Dickinson/Jamieson) Correction Factor
M
1.5 if MC = 2
1.26
1.26 otherwise
1 g Cd
g CL
Fdj 1 M
f
L
1.004561112938649
Ydj
CALIB053
Ych
1 g Cd
g CL
1 g
g
f
L
f
L
0.003619931
0.004455698537
Page 6 of 8
L
Kch
f
Fchm
0.25
f
L
Printed at 21/07/2016
0.25
2.881776634396
Interim Results
Throat Tapping Pressure, Pa
P2 2529225
P3 2592451.445286918
tf 18
Throat Diameter, m
dop 0.16112581062
L 0.02736609099
0.977494961548354
Fwg
Fdj if MC = 1 MC = 2
Fwg 1.019120579237104
Fchm if MC = 3
Flew if MC = 4 MC = 5
0 otherwise
Dry Gas Gross Volume Flow Rate, m3/s
qvgG
qmL
qmu
Fwg f
qmu
Fwg
1.869790574301
1 g
1.058224559419
Beta Ratio
0.700808977035
ReD 17988829.24
Cd 0.9973582849846854
Results
Uncorrected Mixture Mass Flow Rate (kg/s)
qmu 39.40856684
qmG
qvsG
qeG
qmu
Fwg
qmG
b
38.66918954
49.897472530429
CVg qmG
CVg qvsG
if GasCV = 1
otherwise
qmc c qm 0.42725556199
qvsc
CALIB053
qec
qmc
15c
3492.823077130001
0.000674330117
CVc qmc
CVc qvsc
if CondCV = 1
0.047203108174
otherwise
qmhc q mG q mc 39.096445100439
Page 7 of 8
qmm m q m 0.203713435439
qvsm
q mm
0.000269819105
15m
qmw w qm 0.42725556199
qvsw
(Sm3/s)
(Sm3/s)
Printed at 21/07/2016
qmw
0.000411614011
15w
qmhcThr qmhc
qmGThr q mG
qmuThr q mu
3600
1000
3600
1000
3600
1000
140.747202361579
139.209082338416
141.870840627792
qeG.hr q eG
qec.hr qec
3600
6
12.574163077668
10
3600
0.000169931189
6
10
3600
qehc.hr qehc
12.574333008857
6
10
CALIB053
Page 8 of 8