Nonlinear Problem
Nonlinear Problem
RT a
p
v b v(v b) T
where R = the universal gas constant [= 0.518 kJ/(kg K)], T = absolute temperature (K), p =
absolute pressure (kPa), and v = the volume of a kg of gas (m3/kg). The parameters a and b are
calculated by
R 2TC 2.5 TC
a 0.427 b 0.0866R
pC pC
where pc = 4600 kPa and Tc = 191 K. As a chemical engineer, you are asked to determine the
amount of methane fuel that can be held in a 3-m3 tank at a temperature of 40 C with a pressure
of 65,000 kPa.
Use a root-locating method of your choice to calculate v and then determine the mass of methane
GIVEN:
R 0.518kJ / (kgK)
p 65000kPa
pC 4600kPa
TC 191K
vL 0.0026
vU 0.003
REQUIREMENTS:
A. the specific volume v and the mass using the exact solution of finding the roots
B. the specific volume v and the mass using the bisection method with 3 iterations with the
C. the specific volume v and the mass using the newton-raphson method with 3 iterations with the
D. the specific volume v and the mass using the secant method with 3 iterations with the absolute
E. the specific volume v and the mass using the regula-falsi method with 3 iterations with the
RT a
f(v) p
v b v(v b) T
TC 191
And b 0.0866R 0.0866(0.518) 0.001862615
pC 4600
(0.518)(233) 12.5577832
f(v) 65000
v 0.001862615 v(v 0.001862615) 233
120.694 0.822687731
f(v) 65000
v 0.001862615 v(v 0.001862615)
A. To determine the specific volume v and the mass using the exact solution of finding the roots,
0.002807378m3/kg
-0.000475273m3/kg
-0.000475273m3/kg
V
m where V = 3m3 and v = 0.002807378m3/kg
v
3
m
0.002807378
m 1068.612776762kg
B. To determine the specific volume v and the mass using the bisection method with 3 iterations
To check,
f(vU ) * f(v L ) 0
f(0.003) * f(0.0026) 0
1528.082313056 * 27774.156684279 0
424391400.311505500 0
Then the root is bounded and the approximated specific volume is,
vU v L 0.003 0.0026
vM
2 2
vM 0.0028m3 / kg
V 3
mM
v M 0.0028
mM 1071.428571429kg
To check,
f(vM ) * f(v L ) 0
f(0.0028) * f(0.0026) 0
740.554048453 * 27774.156684279 0
20568264.174918566 0
Then the root is not bounded and the approximated specific volume is,
vU v M 0.003 0.0028
vM1
2 2
vM1 0.0029m3 / kg
V 3
mM1
v M1 0.0029
mM1 1034.482758621kg
|A | 3.448275862%
| present past |
|A | * 100%
present
| 1034.482758621 1071.428571429 |
|A | * 100%
1034.482758621
|A | 3.571428571%
To check,
f(vM1 ) * f(vM ) 0
f(0.0029) * f(0.0028) 0
8820.592728173 * 740.554048453 0
6087793.225534353 0
Then the root is bounded and the approximated specific volume is,
v M1 v M 0.0029 0.0028
v M 2
2 2
vM2 0.00285m3 / kg
V 3
mM2
v M2 0.00285
mM2 1052.631578947kg
A 1.754385965%
A 1.724137931%
C. To determine the specific volume v and the mass using the newton-raphson method with 3
vI 0.003m3 / kg
f vI 15280.082313056
vI1 vI 0.003
f ' vI 62901414.370263815
vI1 0.002757079m3 / kg
V
mI
vI
3
mI
0.003
mI 1000.000000000kg
V 3
mI1
vI1 0.002757079
mI1 1088.108102815kg
A 8.810810281%
A 8.097366667%
vI1 0.002757079m3 / kg
f vI1 5343.353111776
vI2 vI1 0.002757079
f ' vI1 113446040.096620220
vI2 0.002804179m3 / kg
V 3
mI2
vI2 0.002804179
mI2 1069.831697331kg
A 1.679636000%
A 1.708343988%
vI2 0.002804179m3 / kg
f v I 2 319.343532363
v I 3 v I 2 0.002804179
f ' v I 2 100250902.714098590
vI3 0.002807364m3 / kg
V 3
mI3
vI3 0.002807364
mI3 1068.617937198kg
A 0.113451622%
A 0.113582235%
D. To determine the specific volume v and the mass using the secant method with 3 iterations with
vI 0.003m3 / kg
vI1 0.0026m3 / kg
f vI [vI vI1 ]
vI1 vI
f vI f(vI1 )
f vI f 0.003 1528.082313056
vI1 0.002858039m3 / kg
V 3
mI1
vI1 0.002858039
mI1 1049.670856945kg
f vI1 [vI1 vI ]
vI2 vI1
f vI1 f(vI )
f vI f 0.0026 27774.156684279
1528.082313056[0.002858039 0.0026]
vI2 0.002858039
1528.082313056 27774.156684279
vI2 0.002820503m3 / kg
V 3
mI2
vI2 0.002820503
mI2 1063.640089963kg
A 1.330828636%
A 1.313342093%
1283.432429450[0.002820503 0.002858039]
vI3 0.002820503
1283.432429450 (4727.975225643)
vI3 0.002806517m3 / kg
V 3
mI3
vI3 0.002806517
mI3 1068.940561746kg
A 0.498335374%
A 0.495862162%
E. To determine the specific volume v and the mass using the regula-falsi method with 3 iterations
To check,
f(vU ) * f(v L ) 0
f(0.003) * f(0.0026) 0
1528.082313056 * 27774.156684279 0
424391400.311505500 0
Then the root is bounded and the approximated specific volume is,
f v U v L v U
vR vU
f vL f vU
( 1528.082313056)(0.0026 0.003)
vR 0.003
27774.156684279 (1528.082313056)
vR 0.002858039m3 / kg
mR 1049.670856945kg
To check,
f(vR ) * f(v L ) 0
f 0.002858039 * f (0.0026) 0
4727.975225643 * 27774.156684279 0
424391400.311505500 0
Then the root is bounded and the approximated specific volume is,
f v R v L v R
v R 1 v R
f vL f vR
( 4727.975225643)(0.0026 0.002858039)
vR 1 0.002858039
27774.156684279 (4727.975225643)
vR1 0.002820503m3 / kg
V 3
mR 1
v R 1 0.002820503
mR1 1063.640089963kg
A 1.330826452%
For the absolute relative approximate error in percent in mass,
A 1.313342093%
To check,
f(vR 1 ) * f(vL ) 0
f 0.002820503 * f (0.0026) 0
1528.082313056 * 27774.156684279 0
-42441197.589292817 0
Then the root is bounded and the approximated specific volume is,
f v R 1 v L vR 1
v R 2 vR 1
f v L f v R 1
( 1283.432429450)(0.0026 0.002820503)
vR 2 0.002820503
27774.156684279 (1283.432429450)
vR2 0.002810764m3 / kg
V 3
mR 2
vR 2 0.002810764
mR2 1067.325582912kg
A 0.345301659%
SUMMARY:
NONLINEAR APPROXIMATION
APPROXIMATE
METHOD REQUIREMENT ITERATION |A| (%)
VALUE
SPECIFIC 1 0.0028 ---
VOLUME 2 0.0029 3.448275862
BISECTION (m3/kg) 3 0.00285 1.754385965
METHOD 1 1071.428571429 ---
MASS
2 1034.482758621 3.571428571
(kg)
3 1052.631578947 1.724137931
SPECIFIC 1 0.002757079 8.810810281
VOLUME 2 0.002804179 1.679639000
NEWTON- (m3/kg) 3 0.002807364 0.113451622
RAPHSON
1 1088.108102815 8.097366667
METHOD MASS
2 1069.831697331 1.708343988
(kg)
3 1068.617937198 0.113582235
SPECIFIC 1 0.002858039 ---
VOLUME 2 0.002820503 1.330828636
SECANT (m3/kg) 3 0.002806517 0.498335374
METHOD 1 1049.670856945 ---
MASS
2 1063.640089963 1.313342093
(kg)
3 1068.940561746 0.495862162
SPECIFIC 1 0.002858039 ---
VOLUME 2 0.002820503 1.330826452
REGULA-FALSI (m3/kg) 3 0.002810764 0.346489424
METHOD 1 1049.670856945 ---
MASS
2 1063.640089963 1.313342093
(kg)
3 1067.325582912 0.345301659