Numec

Download as rtf, pdf, or txt
Download as rtf, pdf, or txt
You are on page 1of 14

MEC500 Numerical Methods With Applications | 2017

FACULTY OF MECHANICAL ENGINEERING

Program : Bachelor of Mechanical Engineering (Hons.) EM220


Course : Numerical Methods With Applications :
Course code :MEC500
Lecturer : Mardziah Che Murad
Group : EMD6M9A

TECHNICAL REPORT (INDIVIDUAL ASSIGNMENT)

NO NAME STUDENT ID SIGNATURE

1 Dinie bin Abdullah zamawi 2014445493

Report Submission : 23 May 2017 Lecturers Validation : _________________


(Date) (Signature)
MEC500 Numerical Methods With Applications | 2017

TABLE OF CONTENTS

PAGE

1.0 Introduction 1

2.0 Objective(s) 1

3.0 Procedures 2

4.0 Sample Calculation. 3

5.0 Result.. 6

6.0 Discussion 7

7.0 Conclusion.. 7

8.0 Reference 8

9.0 Appendix. 8
MEC500 Numerical Methods With Applications | 2017

1.0 INTRODUCTION

The ideal gas law is given by

(1)

Where p is the absolute pressure, V is the volume, n is the number of moles, R is the universal gas
constant and T is the absolute temperature. Although this equation is widely by engineer and scientist,
it is accurate over only a limited range of pressure and temperature.

Furthermore, Eqs (1) is more appropriate for some gases that for others. An alternative equation of
state for gasses is given by:

(2)

Known as the van der Waals equation, where v=V/n is the molal volume and a and b are empirical
constants that depend on the particular gas.

2.0 OBJECTIVE(S)

1. To analyze and solve the problem using numerical technique(s).


2. To apply engineering problem solving method.
3. To use computer computational tools (MATLAB/Excel/etc.) to verify the result.
1
MEC500 Numerical Methods With Applications | 2017

3.0 PROCEDURES

1. v0 was calculated for temperature 300 K by using ideal gas law.


2. The equation of function, f(v) was calculated by using equation of state for gases according
to van der Waals equation.
3. Then, the equation of function, f(v) was differentiated to get f (v).
4. Next, the value of molal volume, v was determined by Newton-Raphson iterative method.
5. The value in (4) then was substituted into (2) and (3).
6. Then, the calculations were repeated until third iteration.
7. Next, find the approximate error for each iteration.
8. Step (1) to (7) were repeated using temperature of 400 K and 500 K.
9. All calculations were made using Microsoft Excel.
2
MEC500 Numerical Methods With Applications | 2017

4.0 CALCULATION

Given
P=2.5atm
T=300K,400K,500K
R=0.082054
a=14.09
b=0.0994

Find
Molar volume of acetone.

Theory
pV = nRT -ideal gas equation (1)
a
( p + 2 )(v - b) = RT -Van der waals equation (2)
v

-Newton Raphson method Equation (3)


Solution

a) T=300K

b) T=400K

c) T=500K

Rearrange equation (2)

Substitute value,

------(4)

Differentiate equation,

-------(5)

Note that :
Use result from (1) as initial guess ( and stop when . The true
value of V=9.247621.

(Sample calculation T=300K)


First iteration ;
For T=300K

i vi Vi+1 T (%)
0.720
0 9.84648 9.314488
1 9.314488 9.248603 0.008

For T = 300 K, molal volume of acetone is 9.248603 L/mol

For T=400K

i vi vi+1 T (%)
0.220
0 13.12864 12.722785
1 12.727785 12.694812 0.000

For T = 400 K, molal volume of acetone is 12.694812 L/mol

For T=500K

i vi vi+1 T (%)
0.090
0 16.4108 16.074818
1 16.074818 16.060137 0.001

For T = 500 K, molal volume of acetone is 16.060137 L/mol


MEC500 Numerical Methods With Applications | 2017

5
MEC500 Numerical Methods With Applications | 2017

5.0 RESULT

6
MEC500 Numerical Methods With Applications | 2017

6.0 DISCUSSION

The ideal gas law assumed that the value of molecular size and intermolecular force are approaching
zero. The equation was modified into van der Waals equations that consider the characteristic of
individual gases. The constant and are introduced to represents the intermolecular forces and
molecular size. Thus, the volume of acetone that had been calculated is more precise by using the van
der Walls equation.
Newton-Raphson method was used in this analysis. The reasons are this method has fastest way
to calculate the solution with less iteration. Unlike other methods, for example, graphical method
seems like easy but it require graph paper and hand skill to draw the best curve to get the value of
roots. It also not had formula to perform this method. Then, bisection method needs two initial guesses.
NR method just needs one guess to find the answer. So, NR method is the best method to solve the
problem given.

7.0 CONCLUSION

In conclusion, I was able to analyze and solve the problem using numerical technique, which
is Newton-Raphson iterative method. Next, I was able to use computer computational tool to
verify the result, which is Microsoft Excel. Molal volume of acetone, (v) for temperature of 300
K, 400 K and 500 K are 9.247621 L/mol, 12.694972 L/mol and 16.070247 L/mol.

8.0`REFERENCE

1. 1. Law, S., & Lielmezs, J. (1985). Modified Van der Waals equation of state for
saturated vapourliquid equilibrium. Thermochimica Acta, 84, 71-81. doi:10.1016/0040-
6031(85)85375-2
2. Van der Waals equation. (2017, April 17). Retrieved May 21, 2017, from
https://en.wikipedia.org/wiki/Van_der_Waals_equation
3. (n.d.). Retrieved May 21, 2017, from http://hyperphysics.phy-
astr.gsu.edu/hbase/Kinetic/waal.html

7
MEC500 Numerical Methods With Applications | 2017

9.0 APPENDIX
For T=300K

Figure 1 : NR Method with T= 300K


For T= 400K

Figure 2 : NR Method with T= 400K

For T=500K

Figure 3 : NR Method with T= 500K

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy