Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

1-2P-Flash Separator - Soka

Download as xls, pdf, or txt
Download as xls, pdf, or txt
You are on page 1of 27

COSTUMER DATE ISSUED ISSUED FOR PROJECT NAME

PT. DOOSAN INDONESIA 14-Jan-02 PROPOSAL MEDIUM PRESSURE


CLIENT CALCULATED BY CHECKED BY WELL TEST SEPARATOR
PT MEDCO E&P IND. AST ASY (OPTIONAL)

Gas Capacity Constrain Formulation


Vessel is 25% full of liquid. Gas Velocity, Vg. A is in ft2, D in ft, d in inches, Q in ft3/s
Q
V g=
Ag

( ) ( ) d2
2
75 π 2 75 π d
A g= 100 D = 100 =244
4 4 144

Qg is in MMscfd
scf day hr 14 . 7 TZ
Q=Q g X 106 X X X X
MMscf 24 hr 3 , 600 s P 520
TZ
=0. 327 Q
P g

Vg =
( 0 . 327
TZ
P
Q )
x
g
244
d2
TZQ g
V g= 80
2
Pd
Set the residence time of the gas equal to the time required for
the droplet to fall to gas-liquid interface:
Leff D d
t g= t d= =
V 2 V t 24 V t
g
Leff
t g=
80
( TZQ g
Pd 2 )
[( ]
1/ 2
Remembering that V t = 0. 0119
ρl CD)
ρ l− ρ g d m

[( ]
1/ 2

t d=
d
( 24 ) ( 0 . 0119 )
ρg
ρ l− ρ g ) CD
dm

Derived and Calculated


EQ-Derivation Page 1 By Ahmad Santoso
[( ]
1/ 2

t d=
d
( 24 ) ( 0 . 0119 )
ρg
ρ l− ρ g ) CD
dm

Setting tg = td

[( ]
1/ 2

Leff
=
d
ρg CD
)
ρ l− ρ g d m

( )
80 TZQ g (24 ) (0 . 0119 )
Pd 2

So for a vessel 25% full of liquid, and seperation of 100-micron liquid drops
from the gas, the following equations apply;

[( ]
1/ 2

Leff d = 280
TZQ g
P
ρg
ρl − ρ g ) CD
dm

Leff = Effective length of vessel where seperation occures, ft


d = Vessel internal diameter, in
T = Operating Temperature, 8R
Z = Gas Copressibility
Qg = Gas Flow Rate, MMscfd
P = Operating Pressure, psia
rl = Density of liquid, lb/ft3
rg = Density of gas, lb/ft3
CD = Drag coefficient
dm = Liquid drop do to be seperated, micron

LIQUID RETENTION TIME CONSTRAIN FORMULATION

t is in s, Vol in ft3, Q in ft3/s, D in ft, d in inches, Leff in ft

Vol
t=
Q

( πD 2 Leff
)
πd 2 Leff 2
Vol=
25 25
= 100 ( 4 )( 144 ) =
1.36E-03 d Leff
100
4

( Vol )o = 1.36E-03 d 2 L
eff ( )
Ao
Al
( Vol )w = 1.36E-03 d 2 L
eff ( )
Aw
Al

Derived and Calculated


EQ-Derivation Page 2 By Ahmad Santoso
Qo and Qw are in bpd

3
ft day 24 hr
Q=Qo X 5 . 61 X X =6 . 49 x 10−5 Qo
barrel 24 hr 3 ,600 s
−5
Q=6 . 49 x 10 Q w
2
21
d Leff
t=
Ql
Liquid Capacity Formula for Two Phase Sizing:

t r Ql
d 2 Leff = 0.350

Ao, Aw and A, are cross-sectional areas of oil, water and liquid

21
( ) Ao
Al
t Q
= 2o o
d L eff
21
( )
Aw
Al
t Qw
= w
d 2 L eff

(tr)o and (tr)w are in minute

( )
21 A o
60 A
l
t Q
= 2o o
d L eff
21
60 ( )
Aw
Al
t w Qw
= 2
d L eff

60(
21 A o + A w
A l
=
)t o Qo + tw Qw
d2 L eff

Liquid Capacity Formula for Three Phase Sizing:

d Leff =2.856 [( t r )o Q o +( t r )w Q w ]
2

Derived and Calculated


EQ-Derivation Page 3 By Ahmad Santoso
RESULTS - HORIZONTAL VESSEL DESIGN
----------------------------------
Gas Flow Rate Converted to Macf/d … 3.695
Compressibility Factor … 0.989
Fluid Density, lb/ft3:
Gas ..... 0.21
Oil ..... 51.02
Water ... 61.71
2 - Phase 3 - Phase
Vessel Vessel
Critical Gas Velocity, ft/sec ..... 6.64 6.86
Gas Flow Area, ft2 ................ 7.57 8.51
Liquid Flow Area, ft2 ............. 14.06 15.81
Gas Capacity, Mscf/d .............. 14.69 17.06
MINIMUM I.D., inches .............. 62.97 66.77
MINIMUM GAS FLOW LENGTH, feet ..... 20.80 22.20
Length to Diameter Ratio .......... 3.96 3.99
Outside Diameter, inches .......... 66.00 72.00
Min. Shell Thickness, inches ...... 0.495 0.528
Ellipsoidal Head Thickness, inches. 0.493 0.526
Estimated Weight Empty, k-lbs...... 8.95 11.16
Water Filled, k-lbs ............. 38.54 48.77
Fluid Calculations

Viscosity (Cp)

Formula
m = 10x - 1 m = #REF! cp
By Beggs and Robinson (the data set from which the above relationship
was obtained included a range of between 168 and 588 API and 708 F to
2958 F.
m = oil viscosity, cp
T = oil temperature, 8F T = 150 8F
x = y(T)-1.163 x = #REF!
y = 10 z y = #REF!
z = 3.0324-0.02023G z = #REF!
G = oil gravity, 8API G = #REF! 8API

Gas Capacity Constraint in Horizontal Separator Sizing

Three Phase Test Separator


Gas Capacity Constraint: T = 150 8F = 610 8R

)[( ]
1 /2 Z = 0.95
d Leff = 280 ( TZQ g
P
ρg
) CD
ρl − ρg d m
Qg
P
=
=
#REF! MMscfd
64.6959 Psia
rl = 48.80 lb/ft3
= #REF! rg = #REF! lb/ft3
CD = 4.99
dm = 500 Micron

Two Phase Production Separator


Gas Capacity Constraint: T = 100 8F = 560 8R

)[( ]
1 /2 Z = 0.95
d Leff = 280
( TZQ g
P
ρg
) CD
ρl − ρg d m
Qg
P
=
=
1 MMscfd
49.7 Psia
rl = 48.80 lb/ft3
= 37.57451533475 rg = 0.15 lb/ft3
CD = 4.99
dm = 100 Micron

Fluid Properties Page 5 Calculated By Ahmad Santoso


Page 6

Horizontal 2-Phase LP Flash Separator Sizing For Soka Station PT. MEPI

GIVEN:
rl = 48.80 lb/ft3

rg = 0.15 lb/ft3

dm = 100 micron
m = 0.0104 cp

P = 49.7 Psia

T = 560 8R

From Iteration
Vt = 0.950

Re = 6.845

CD = 4.9927

Qo = 500 bopd

Oil = 45.4 8API


DS.G. = 0.2200

(S.G.)o = 0.8

(S.G.)w = 1.02

Sg = 0.7

(tr)o = 3 min

(tr)w = 3 min

Qw = 500 bwpd

Qg = 35 MMscfd

MW = 20.3 Z = 0.95
(From GPSA)
Gas Capacity Constraint

)[( ]
1 /2

d Leff = 280 (
TZQ g
P
ρg
)
CD
ρl − ρg d m
= 37.574515

Diameter vs. Length for Gas Capacity Restraint


d Leff 12Lss/d Remark
(in) (ft) For Slenderness Ratio

30 1.25 0.50 Denied


32 1.17 0.44 Denied
34 1.11 0.39 Denied
36 1.04 0.35 Denied
38 0.99 0.31 Denied
40 0.94 0.28 Denied
42 0.89 0.26 Denied
44 0.85 0.23 Denied
46 0.82 0.21 Denied
48 0.78 0.20 Denied
50 0.75 0.18 Denied
52 0.72 0.17 Denied

The accepted size is a minimum requirement

file:///conversion/tmp/activity_task_scratch/735781123.xls Calculated By Ahmad Santoso


Page 7

Liquid Capacity Constraint

t r Ql
d 2 Leff = 0.350
= 8567.415

Horizontal Separator Diameter vs. Length


d Liquid Gas Lss
d2 12Lss/d
(in) Leff Leff (ft) Remark
26 12.67 1.45 676 14.84 6.8 Denied
28 10.93 1.34 784 13.26 5.7 Denied
30 9.52 1.25 900 12.02 4.8 Denied
32 8.37 1.17 1024 11.03 4.1 Denied
34 7.41 1.11 1156 10.24 3.6 Accepted
36 6.61 1.04 1296 9.61 3.2 Accepted
38 5.93 0.99 1444 9.10 2.9 Accepted
40 5.35 0.94 1600 8.69 2.6 Denied
42 4.86 0.89 1764 8.36 2.4 Denied
44 4.43 0.85 1936 8.09 2.2 Denied
46 4.05 0.82 2116 7.88 2.1 Denied
48 3.72 0.78 2304 7.72 1.9 Denied
50 3.43 0.75 2500 7.59 1.8 Denied
52 3.17 0.72 2704 7.50 1.7 Denied
Slenderness ratio (12Lss/d) Choices in the range of 2.75 to 4

Gas Capacity Constraint


14 18 22
44
Vessel Dia, in

36

28

20

Vessel Length, ft
Note : (1) Liquid Capacity Restrain Does Not Govern

The reasonable sizes to choose from above graph that does not violete Liquid
capacity restraint is 66 in x 16' - 6" Seam to Seam

Liquid Level

Normal Liquid Level :


(from grafic on next page)
- Percentage of total Volume is 75%
- Percentage of total Diameter 100 H/D = 72
- Normal Liquid Level is H= 72 D
(used to be normal oil level = 72 x 38

file:///conversion/tmp/activity_task_scratch/735781123.xls Calculated By Ahmad Santoso


Page 8

in three phase separator) = 27.36 Inch

High High Liquid Level (HHLL) = 31 Inch


High Liquid Level (HLL) = 29 Inch
Normal Liquid Level (NLL) = 27 Inch
Low Liquid Level (LLL) = 7 Inch
Low Low Liquid Level (LLLL) = -13 Inch

75%

file:///conversion/tmp/activity_task_scratch/735781123.xls Calculated By Ahmad Santoso


Page 9

DEMISTER SERVICE AREA SIZING CALCULATION

For general design, the following equation can be used for determining the mesh area.

V m =K
√ ρ l− ρ g
ρg

Where
Design
Vm = velocityphase
Gas
rg = density
Liquid
Mesh phase
rl = density
capacity
K = factor.

Typical K values are in the 0.15-0.25 m/s range for vane type and 0.1 m/s for mesh type.
For horizontal vessels, the required demister area (Ad) is given by:

TZQ g
0 .327
P
A d=
Kd
√ ρ g − ρl
ρg

For Vane Type Mist: For Mesh Type Mist:


= 10.471 ft2 = 20.941 ft2
Dd = 3.65 ft Dd = 5.16 ft
= 43.82 in = 61.96 in
Or Ld = 3.00 in
Wd = 14.6 in

Calculating Vane Efficiency

Several authors [2,3,4] use the theory of inertial impaction to provide a meaningful way to estimate the
droplet removal efficiency of a vane mist eliminator. These theories assume the gravity effect to be
negligible (valid for small droplets) and consider the gas-phase drag force on the droplets to be of
prime importance. Use these assumptions, together with Newton’s 1st law to formulate the following
efficiency equation for vanes:

h(Dp) = 1 - e - [rp*Dp2*Vv*N*Q / (515.7*mg*S*Cos2 Q)] = 0.99998


99.998%
Where h(Dp) = Fractional removal efficiency of droplets of diameter Dp.
rp = Droplet density (g/cm3) = 0.78072 (g/cm3)
Dp = Droplet diameter (cm) = 0.001 (cm)
Vv = Superficial vapor velocity (cm/sec) = 28.943801033 (cm/sec)
N = Number of successive blades = 3
Q = Angle of inclination of the vane blades to the vapor flow direction (deg) = 45 (deg)
µg = Vapor viscosity (g/cm-sec) = 0.000104 (g/cm-sec)
S = Spacing between parallel plates (cm) = 25.4 (cm)

file:///conversion/tmp/activity_task_scratch/735781123.xls Calculated By Ahmad Santoso


NOZZLE CALCULATION

Given:

Gas Flow Rate : Qg = 35 MMscfd


Oil Flow Rate : Qo = 250.00 BPD
Water Flow Rate : Qw = 250.00 BPD
Gas S.G. : S = 0.7
Liquid S.G. : S.G. = 0.91
Gas Desity : rg = 0.70 lb/ft3
Liquid Density : rl = 48.80 lb/ft3
Operating Pressure : P = 150.00 Psia
Operating Temperature : T = 560 8R
Gas Compressibility : Z = 0.95
Gas/Liquid Ratio : R = 70,000.00 ft3/bbl

Inlet Nozzle :
(Two-Phase Line)

[ ]
1

d=
(
11. 9+
ZRT
Q
16 . 7 P l ) 2

1 , 000 V
Criteria
Vmax = 60 ft/s or Ve
Vmin = 10 to 50 ft/s

IDmin = 86.25012
[V]1/2

V ID
10 27.27
50 12.20
60 11.13

ID V
4-inch Sch. 40 12.000 51.660 Denied
6-inch Sch. 40 14.000 37.955 Accepted
8-inch Sch. 40 16.000 29.059 Accepted

Use 12-inch Sch. 40

FWKO-PTEN Page 10 Calculated by Ahmad Santoso


Gas Outlet Nozzle :
(Gas-Line)

[ ]
1
60 Q g TZ 2
d=
VP

Criteria
Vmax = 60 ft/s or Ve
Vmin = 10 to 50 ft/s

IDmin = 86.3018
[V]1/2

V ID
10 27.29
50 12.20
60 11.14

ID V
3-inch Sch. 40 12.000 51.722 Denied
4-inch Sch. 40 14.000 38.000 Accepted
6-inch Sch. 40 16.000 29.094 Accepted

Use 4-inch Sch. 40

Liquid Outlet Nozzle :


(Liquid-Line)

[ ]
1
Ql 2
d=0 . 012
V
Criteria
Vmax = 15 ft/s
Vmin = 3 ft/s

IDmin = 2.45
[V]1/2

V ID
3 1.41
15 0.63

ID V
4-inch Sch. 40 1.000 6.000 Accepted

FWKO-PTEN Page 11 Calculated by Ahmad Santoso


6-inch Sch. 40 3.000 0.667 Denied
8-inch Sch. 40 5.000 0.240 Denied

Use 6-inch Sch. 40

FWKO-PTEN Page 12 Calculated by Ahmad Santoso


PT ASDAR WAHYU BERKAH -CONSULTING ENGINEERS- (62)(021) 801-6734
Pressure Drop Calculations
Inputs
Nominal Pipe Diameter 1 inches Inside Diameter 1.049
Schedule 40
Flow Rate 68.17099567 GPM Click this Button
Viscosity 0.6564 cps to see a list of
Density 6.55 lbs/gal acceptable pipe diameters
Absolute Roughness 0.00015 ft and pipe schedules

Reynolds Number Delta P Per 100ft. of pipe Delta P Per 100ft. of pipe
245,479 90.94 psi 90.45 psi
267.13 feet 265.69 feet
Linear Velocity Friction Factor (Darcy) Friction Factor (Darcy)
25.28 ft/sec 0.0235 0.0233
Churchill Equation Chen Equation
(1977) (1979)
Gas Liquid
62.33766234 5.8333333

INDEX 5 INDEX2
Nominal Pipe 5 10 20 30 40 60
0.125 no such Sch. 0.307 no such Sc no such Sch. 0.269 no such Sch.
0.25 no such Sch. 0.41 no such Sc no such Sch. 0.364 no such Sch.
0.375 no such Sch. 0.545 no such Sc no such Sch. 0.493 no such Sch.
0.5 0.71 0.674 no such Sc no such Sch. 0.622 no such Sch.
0.75 0.92 0.884 no such Sc no such Sch. 0.824 no such Sch.
1 1.185 1.097 no such Sc no such Sch. 1.049 no such Sch.
1.25 0.065 0.109 no such Sc no such Sch. 1.38 no such Sch.
1.5 1.77 1.682 no such Sc no such Sch. 1.61 no such Sch.
2 2.245 2.157 no such Sc no such Sch. 2.067 no such Sch.
2.5 2.709 2.635 no such Sc no such Sch. 2.469 no such Sch.
3 3.334 3.26 no such Sc no such Sch. 3.068 no such Sch.
3.5 3.834 3.76 no such Sc no such Sch. 3.548 no such Sch.
4 4.334 4.26 no such Sc no such Sch. 4.026 no such Sch.
5 5.345 5.295 no such Sc no such Sch. 5.047 no such Sch.
6 6.407 6.357 no such Sc no such Sch. 6.065 no such Sch.
8 8.407 8.329 8.125 8.071 7.981 7.813
10 10.482 10.42 10.25 10.136 10.02 9.75
12 12.438 12.39 12.25 12.09 11.938 11.626
14 13.624 13.5 13.375 13.25 13.124 12.814
Acceptable Pipe Schedules:
5 80 5 S XS TYPE THESE IN
10 100 10 S XSS AS A STRING I.E.
20 120 40 S '10S
30 140 80 S
40 160
60

Acceptable Nominal Pipe Diameters in inches:

1/8 1-1/4 4
1/4 1-1/2 5 To return to the input sheet
3/8 2 6 click the input button
1/2 2-1/2 8 below.
3/4 3 10
1 3-1/2 12

{goto}a200~ \m

{BREAKOFF} \0
/WGPE
{GOTO}preday~
{Wait @now+@time(0,0,8)}
{HOME}
{down}
{down}

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
*
*
* PIPEFLOW VERSION 2.0
* 1993 & 1995
*
* R.C.COSTELLO & ASSOC., INC.
*
* -CONSULTING ENGINEERS-
*
* (please wait)
*
* For best results set:
* Zoom Control to 84% BY ROCKY C. COSTELLO
* and have two (2)
* Toolbars Showing November 15 , 1995
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *

FRICTION FACTOR EQUATIONS

For information on the original Colebrook-White Equation


Click Here Pressure Drop
Equation

For information on the Churchill Equation


Click Here References
on equations

For information on the Chen Equation


Click Here Return to input page
Colebrook-White Equation

k
1 D 1.2613
4. log k=Roughness in feet
3.7065
f Re. f
D=inside diameter in feet

Re= Reynolds number


This is the original formula for calculating friction factor. It is implicit in friction
factor f , therefore it is not used. Friction factor (f) must be solved by trial and error.
In addition in the laminar flow regime, one must use f=64/Re, for Re < 2,100
See references for further detail.

Churchill Equation

16
0.9 16
7 k 37530
A2 2.457. ln 0.27. A3
Re D Re

1 k=Roughness in feet
12
12
8 1
f 8. D=inside diameter in feet
Re 3
This is a very good formula and in most cases agrees with the Chen equation.
2
A 2 use
In very low Reynolds numbers A this
3 (i.e. transition and laminar
Re=ranges.)
Reynolds number
This is a full regime equation.

Chen Equation

1.1098
k
0.8981
D 7.149
A4
2.8257 Re

2
1 .4
f k=Roughness in feet
k
D 5.0452 . D=inside diameter in feet
4. log log A 4
3.7065 Re
Re= Reynolds number
2
1 .4
f k=Roughness in feet
k
D 5.0452 . D=inside diameter in feet
4. log log A 4
3.7065 Re
Considered the most accurate. Re= Reynolds number
Use only in turbulent range.

Instructions

1. Input a nominal pipe diameter and a pipe schedule. The program automatically
looks up the correct I.D. in a table. For the (S) schedules type an apostrophe before
typing the schedule. For example '10s. All other schedules type in 40 or 80 etc.

2. Input flowrate of the fluid in GPM, density of the fluid in lbs/gallon and viscosity of the
fluid in centipoise.

3. For pipe roughness use .00015 ft.for commercial steel and stainless steel.

4. Program then calculates pressure drop in psi per 100 feet of pipe.

5. A word of caution. This program calculates pressure drop due to frictional losses
only. For difference in elevation use:
Delta P in psig=(H2-H1)*0.4331* specific gravity
Add this to the frictional loss.

Pressure Drop Equation

Length of pipe in feet


Darcy friction factor density in lbs./cubic foot

2
0.000216 . f. L.  . Q flow rate in GPM
P
5
d
pressure in psig
pipe ID in inches

REFERENCES
1. Churchill, S.W. , Friction Factor Equation Spans All Fluid Flow Regimes, Chemical Engineering
Nov. 7, 1977, page 91.

2. Chen, N. H. , An Explicit Equation for Friction Factor in Pipe, Ind. Eng. Chem. Fund., 18, 296, 1979.

3. Gregory G.A. and Forgarasi M., Alternate to Standard Friction Factor Equation, Oil and Gas
Journal, April 1, 1985, page 120

4. Flow of Fluids, By Crane Technical Paper 410, 1979 edition


021) 801-6734

inches 0.0266446 Nominal Pipe Diameter in meters D

0.00171592555 K/D

0.00004572015 Absolute Roughness in meters K

(7/RE)^0.9
8.1203724E-05
A 1.827823584E+20
B 8.9091961592E-14

C 0.00038621707696

80 100 120 140 160 5s


0.215 no such Sch. no such Sch. no such Sch. N/A
0.302 no such Sch. no such Sch. no such Sch. N/A
0.423 no such Sch. no such Sch. no such Sch. N/A
0.546 no such Sch. no such Sch. no such Sch. 0.466 0.71
0.742 no such Sch. no such Sch. no such Sch. 0.612 0.92
0.957 no such Sch. no such Sch. no such Sch. 0.815 1.185
1.278 no such Sch. no such Sch. no such Sch. 1.16 1.53
1.5 no such Sch. no such Sch. no such Sch. 1.338 1.77
1.939 no such Sch. no such Sch. no such Sch. 1.687 2.245
2.323 no such Sch. no such Sch. no such Sch. 2.125 2.709
2.9 no such Sch. no such Sch. no such Sch. 2.624 3.334
3.364 no such Sch. no such Sch. no such Sch. no such Sch. 3.834
3.826 no such Sch. 3.624 no such Sch. 3.428 4.334
4.813 no such Sch. 4.563 no such Sch. 4.313 5.345
5.761 no such Sch. 5.501 no such Sch. 5.187 6.407
7.625 7.437 7.187 7.001 6.813 8.407
9.562 9.312 9.062 8.75 8.5 10.482
11.374 11.062 10.75 10.5 10.126 12.438
12.5 12.125 11.814 11.5 11.188 13.624
* * * *
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
* * * *

Pressure Drop

on equations

Return to input page


16
37530
Re

ness in feet

diameter in feet

nolds number

hness in feet

e diameter in feet

ynolds number
hness in feet

e diameter in feet

ynolds number

./cubic foot

w rate in GPM

hes
(0.27*I10)/I5
0.0004632999 7.5156359616

10s 40s 80s xs xxs


0.307 0.269 0.215 0.215 N/A
0.41 0.364 0.302 0.302 N/A
0.545 0.493 0.423 0.423 N/A
0.674 0.622 0.546 0.546 0.252
0.884 0.824 0.742 0.742 0.434
1.097 1.049 0.957 0.957 0.599
1.442 1.38 1.278 1.278 0.896
1.682 1.61 1.5 1.5 1.1
2.157 2.067 1.939 1.939 1.503
2.635 2.469 2.323 2.323 1.771
3.26 3.068 2.9 2.9 2.3
3.76 3.548 3.364 3.364 N/A
4.26 4.026 3.826 3.826 3.152
5.295 5.047 4.813 4.813 4.063
6.357 6.065 5.761 5.761 4.897
8.329 7.981 7.625 7.625 6.875
10.42 10.02 9.75 9.75 8.75
12.39 12 11.75 11.75 10.75
13.5 13.124 12.5 13
Iteration to Solve the Values Of Vt (Terminal Velocity), Re (Renold's Number) and CD (Drag Coefficient)

Assume CD = 300

(1). Start with

[ ]
1/ 2
( ρ1− ρg ) d m
V t =0. 0204
ρg
(2). Calculate Re
ρo d m V
Re=0 . 0049
μ
(3). From Re, Calculate CD using
24 3
CD= + 1/ 2 +0 .34
Re Re
(4). Recalculate Vt using,

[( ]
1/ 2

V t = 0. 0199
(5). Go to step 2 and iterate
ρ1− ρg d m
ρg CD )
No. of Iteration 1 2 3 4 5 6 7
Vt = 0.105 0.356 0.624 0.797 0.883 0.922 0.938

Re = 0.757 2.567 4.499 5.744 6.368 6.645 6.763

CD = 35.4954 11.5607 7.0895 5.7696 5.2978 5.1154 5.0425

No. of Iteration 8 9 10 11 12 13 14
Vt = 0.945 0.948 0.949 0.949 0.949 0.950 0.950

Re = 6.811 6.831 6.840 6.843 6.844 6.845 6.845

CD = 5.0130 5.0009 4.9960 4.9940 4.9932 4.9929 4.9928

No. of Iteration 15 16 17 18 19 20 21
Vt = 0.950 0.950 0.950 0.950 0.950 0.950 0.950

Re = 6.845 6.845 6.845 6.845 6.845 6.845 6.845

CD = 4.9927 4.9927 4.9927 4.9927 4.9927 4.9927 4.9927

No. of Iteration 22 23 24 25 26 27 28
Vt = 0.950 0.950 0.950 0.950 0.950 0.950 0.950

Re = 6.845 6.845 6.845 6.845 6.845 6.845 6.845

CD = 4.9927 4.9927 4.9927 4.9927 4.9927 4.9927 4.9927

No. of Iteration 29 30
Vt = 0.950 0.950 OK

Re = 6.845 6.845 OK

CD = 4.9927 4.9927 OK

Iteration of Fluid Properties Page 27 Calculated By Ahmad Santoso

You might also like