Amir Azim Bin Mohd Israpil - 2017466174 - EH2434

Ø = 0.2 u = 1 cp Pᴏ = 1 atm
k = 1000 md ∆T = 0.2 sec (initial pressure anywhere)
N = 10 c = 10⁻⁴ atm⁻ᴵ PL = 2 atm
T max = 0.2 sec L = 100 cm PR = 1 atm

1) All grid blocks at time intervals of 0.00025 s

*Transient flow because of time dependent flow

o o o o o o o o o o

←-------------- 100 cm ----------------- →

1st block until 10th block

*If steady state

PL= 2
PR= 1
X (cm) P
XL 0 2
X1 5 2
X2 15 1.9
X3 25 1.8
X4 35 1.7
X5 45 1.6
X6 55 1.5
X7 65 1.4
X8 75 1.3
X9 85 1.2
X10 95 1.1
XR 100 1
 For transient flow

sin
n X X(cm) PL PR Pie K (D) Ø c t (s) tgh
(n.pie.x/L)
1 XL 0 2 1 3.142 1 0.2 0.001 0.00025 0 -0.001234021
1 X1 5 2 1 3.142 1 0.2 0.001 0.00025 0.0027 -0.001234021
2 X2 15 2 1 3.142 1 0.2 0.001 0.00025 0.0164 -0.004936082
3 X3 25 2 1 3.142 1 0.2 0.001 0.00025 0.0411 -0.011106185
4 X4 35 2 1 3.142 1 0.2 0.001 0.00025 0.0767 -0.019744328
5 X5 45 2 1 3.142 1 0.2 0.001 0.00025 0.1231 -0.030850513
6 X6 55 2 1 3.142 1 0.2 0.001 0.00025 0.1800 -0.044424738
7 X7 65 2 1 3.142 1 0.2 0.001 0.00025 0.2469 -0.060467005
8 X8 75 2 1 3.142 1 0.2 0.001 0.00025 0.3231 -0.078977312
9 X9 85 2 1 3.142 1 0.2 0.001 0.00025 0.4073 -0.099955661
10 X10 95 2 1 3.142 1 0.2 0.001 0.00025 0.4977 -0.12340205
10 XR 100 2 1 3.142 1 0.2 0.001 0.00025 0.5213 -0.12340205

exp n Σ P(x,t)
0 1 1 1 1.36346
-3.3319E-06 1 1 2 0.67693
-8.0952E-05 2 0.5 2.5 0.25866
-0.00045646 3 0.33 2.83 -0.05140
-0.00151439 4 0.25 3.08 -0.31053
-0.0037977 5 0.20 3.28 -0.53784
-0.00799645 6 0.17 3.45 -0.74605
-0.0149293 7 0.14 3.59 -0.93517
-0.02551757 8 0.12 3.71 -1.11155
-0.04071194 9 0.11 3.82 -1.28157
-0.0614172 10 0.09 3.91 -1.43886
-0.06432949 10 0.09 4 -1.54615

2.00000
1.50000
1.00000
Pressure, P

0.50000
0.00000
-0.50000 0 20 40 60 80 100 120

-1.00000
-1.50000
-2.00000
Distance, X
2) Explicit formulation , ∆t= 0.00025, 0.0005, 0.00075,0.00125

𝑃𝑖+1 − 2𝑃𝑖 + 𝑃𝑖−1 ∅𝑢𝑐 𝑃𝑖 𝑡+∆𝑡 − 𝑃𝑖

= ×
∆𝑥 2 𝑘 ∆𝑡
i = 1,2,3,4,5,6,7,8,9,10

∆t= 0.00025 0.0005 0.00075 0.00125

X (cm) P(x,t) ∆x Pi^(t+∆t) Pi^(t+∆t) Pi^(t+∆t) Pi^(t+∆t)
XL 0 1.363463
5 0 0 0 0
X1 5 0.676926 -13.013 -26.70293 -40.39285 -67.77271
10 0 0 0 0
X2 15 0.258657 1.611299 2.96394 4.316582 7.021865
10 0 0 0 0
X3 25 -0.0514 0.585137 1.221674 1.858211 3.131286
10 0 0 0 0
X4 35 -0.310535 0.087301 0.485137 0.882973 1.678644
10 0 0 0 0
X5 45 -0.537842 -0.299141 -0.060439 0.178262 0.655665
10 0 0 0 0
X6 55 -0.746053 -0.507352 -0.268651 -0.029949 0.447454
10 0 0 0 0
X7 65 -0.935169 -0.776034 -0.6169 -0.457766 -0.139497
10 0 0 0 0
X8 75 -1.111553 -1.031986 -0.952419 -0.872852 -0.713717
10 0 0 0 0
X9 85 -1.281572 -1.122438 -0.963304 -0.804169 -0.485901
10 0 0 0 0
X10 95 -1.438861 -53.43003 -105.4212 -157.4124 -261.3947
5 0 0 0 0
XR 100 -1.546149

50

0
0 20 40 60 80 100
-50
0.00025
-100 0.0005

-150 0.00075
0.00125
-200

-250

-300
3) plots pressure vs time

Explicit

Block 5
0.8

0.6

0.4
Pressure

0.2
Block 5
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012 0.0014
-0.2

-0.4
time

Block 9
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012 0.0014
-0.2

-0.4
Pressure

-0.6
Block 9
-0.8

-1

-1.2
time

Analytical

Analytical
2.00000
1.50000
1.00000
0.50000
Pressure

0.00000
Analytical
-0.50000 0 2 4 6 8 10 12 14

-1.00000
-1.50000
-2.00000
time
4) re-run part 1,2 with 5 blocks and 10 blocks
5 BLOCK

If steady
X (cm) P
XL 0 2
X1 10 1.9
X2 30 1.7
X3 50 1.5
X4 70 1.3
X5 90 1.1
XR 100 1

For transient flow

n X X(cm) PL PR Pie K (D)

1 XL 0 2 1 3.142 1
1 X1 10 2 1 3.142 1
2 X2 30 2 1 3.142 1
3 X3 50 2 1 3.142 1
4 X4 70 2 1 3.142 1
5 X5 90 2 1 3.142 1
5 XR 100 2 1 3.142 1

sin
n Ø c t (s) tgh
(n.pie.x/L)
1 0.2 0.001 0.00025 0 -0.001234 0 1.0000
1 0.2 0.001 0.00025 0.0055 -0.001234 -6.79E-06 1.0000
2 0.2 0.001 0.00025 0.0329 -0.004936 -0.000162 0.9998
3 0.2 0.001 0.00025 0.0822 -0.011106 -0.000913 0.9991
4 0.2 0.001 0.00025 0.1529 -0.019744 -0.003019 0.9970
5 0.2 0.001 0.00025 0.2443 -0.030851 -0.007537 0.9925
5 0.2 0.001 0.00025 0.2708 -0.030851 -0.008354 0.9917

n Σ P(x,t)
1 1 1 1 P vs X
1 1 2 0 2
2 0.4999 2.4999 -0.4999
3 0.3330 2.8329 -0.8329 1
4 0.2493 3.0822 -1.0822 0
5 0.1985 3.2807 -1.2807 0 50 100 150
-1
5 0.1983 3.479 -1.479
-2