equation of state in fortran program

-

 Equation of state in fortran program

In thermodynamics, different gas laws, for example, Boyle's law, Charles law, have been derived by using the concept of ideal gas. Ideal gas is actually an ideal concept. No ideal gas exists in atmosphere. 

Equation of state is the thermodynamical equation based on concept of ideal gas. Real gases, that are present in the atmosphere, are studied by using this equation. 

For n mole of ideal gas, the ideal gas equation or equation of state is given as,

PV = nRT  ................. (1)

 Where, 

P = pressure

V = volume

T = absolute temperature

R = Universal gas constant.

At normal temperature and pressure, that is, normal conditions, 

P = `1.01 \times 10^(5)` Pascal

V = `22.4 \times 10^(-3)` `m^3`

T= `0^o` C = (0 + 273) K = 273 K

  Then equation (1) gives,

`R =  \frac(P \times V)(n \times T)`

`R = \frac(1.01 \times 10^5 \times 22.4 \times 10^(-3))(n \times 273)`

For one mole ideal gas, n =1,

`R = \frac{1.01 \times 10^5 \times 22.4 \times 10^{-3}}{1 \times 273}`

`\therefore` R = 8.31 `J mol^(-1) K^(-1)`

 This is the value of universal gas constant at normal temperature and pressure conditions.

Also learn, Why speed of light is taken in mass energy relation, `E =mc^2` ?

Fortran program on Equation of State

program eqofstate !PV= nRT
! This program helps you to calculate different variables in equation of state
!==================Technical Things===========
!=============TribhuvanUniversity=====================
!========Central Department of Physics=========
!=================Nepal=======================================
implicit none
real::p,t,v,n ! Variable declaration
real::R1 , R2  !Universal gas constant
integer::r !integer for defining condition of calculation of different variables
write(*,*) 'Euation of State, PV = nRT'
write(*,*) " " !for space
write(*,*) "What do you want to calculate?"
write(*,*) " "
write(*,*) "Note: Enter necessary values in SI units, to get results in SI units"
write(*,*) " "
write(*,*) "Enter 1 to calculate no. of mole(n), enter 2 to calculate pressure(P), "
write(*,*) " "
write(*,*) "Enter 3 to calculate volume(V), enter 4 to calculate temperature(T) or"
write(*,*) " "
write(*,*) "Enter 5 to calculate Universal gas constant(R) at normal conditions"
!Conditions for calculation of different variables of equation of state
read(*,*) r
if (r==1) then ! for calculation of no. of  moles
goto 100
elseif (r==2) then ! for calculation of presure
goto 200
elseif (r==3) then ! for calculation of volume
goto 300
elseif (r==4) then ! for calculation of temperature
goto 400
elseif (r==5) then ! for calculation of gas constant
goto 500
else
stop
end if
!=========================================
100 write(*,*) "To calculate n, enter value of p, v,t"
read(*,*) p,v,t
R1=8.31 ! value of universal gas constant at normal conditions.
n =(p*v)/(t*R1)
write(*,*) "value of n is", n, 'moles'
stop
!========================================
200 write(*,*) "To calculate P, enter value of n, V, T"
read(*,*) n, v, t
R1=8.31
p =(n*t*R1)/v
write(*,*) "value of P is", p, 'Pascal'
stop
!========================================
300 write(*,*) "To calculate V, enter value of n, P, T"
read(*,*) n, p, t
R1=8.31
v =(n*t*R1)/p
write(*,*) "value of V is", v, 'm^3'
stop
!========================================
400 write(*,*) "To calculate T, enter value of n, p, v"
read(*,*) n, p, v
R1=8.31
t =(p*v)/(n*R1)
write(*,*) "value of T is", t, 'K'
stop
!========================================
500 write(*,*) "To calculate R, enter value of n, v, t, p"
read(*,*) n, v, t, p
R2=(p*v)/(n*t)
write(*,*) "value of R is", R2, 'J/K.mol'
!========================================
end program

For downloading fortran file of the program,

Click the Link, Download File.

Follow this blog for more updates.

Comments

Post a Comment

These posts may also be useful for you

top 3 websites for downloading research papers for free

-
  This post provides you valid links of top 3 websites from where you can download research papers published on high ranked journals for free. Just copy the Title or DOI of the papers, then go to either of these websites and paste on search bar and then hit search button.  1. Sci -hub  2. Library genesis 3.Core Download the file to get the links of these websites. Get File Doing research works is not only learning, it's innovation and interpretation of field of interest.  Watch video 👇👇 Watch Now Don't forget to Follow this blog.  Follow this Blog  
Read More »

plot graph in gnuplot from csv and data file | knowledge of physics

-
 If you are using LINUX based OS, for example Ubuntu, you must be familiar with plotting data in gnuplot saved in a file with .dat extension. It is very simple technique.  To plot data in .dat format, you simply go to gnuplot and run following command, gnuplot> plot 'file_name.dat' using 1:2  This command gives you simple plot of your data saved in  file_name.dat. You can further modify this plot by adding title, axes names, legends and colors of plot and width etc. Here 1:2 means first column of data is along x-axis and second column along y-axis. But what about if you have data in a file of  .csv format ? Let's discuss. CSV means coma separated variables. This means each column of the data in the .csv file is separated by coma (,), and gnuplot does not accept this format directly just as you do for the .dat file.  For gnuplot to read data from .csv file, you should tell it that data you are applying are in coma separated format. So, write code as below: gnuplot > un
Read More »

projectile motion in fortran | relation of angle of projection and horizontal range

-
Image
 Have you ever tried to visualize projectile motion in fortran program? This post will help you to generate data and plot the graph of projectile motion.   Above figure shows you the angle of projection(°) versus horizontal range (meter) plot, where the speed of projection is 30 m/s. 👉 Click on link for NEB grade 12  physics numerical problem-solutions:-  LINK Let's start with some basic things: What is projectile and projectile motion ? Projectile:  A projectile is an object thrown to the sky with certain velocity against the earth's gravity. The motion that the projectile possesses during it's journey until it returns to the same surface, from where it was thrown, is called the projectile motion. In above figure, the projectile starts it's journey from point o and follows the curved path to reach another point on the same surface.  Suppose, `\theta` = angle of projection `\therefore` Angle of projection is the angle between the horizontal surface of earth and the
Read More »

NEB physics exam numerical problems with solutions for grade 12 students

-
Image
 This post provides you the numerical problems and their solutions of physics based on NEB(National Examination Board) syllabus grade 12. Solutions of Previously Asked Numerical Problems in NEB exams. Q.1. A coil consisting of 100 circular loops with radius 0.60 cm carries a current of 5A. At what distance from the center, along the axis, the magnetic field has magnitude as `\frac{1}{8}` great as it is at the center? Solution:- Given, number of circular loops = N = 100 radius of the loop, `r = 0.60cm = 6.0 \times 10^-3` m. current, I = 5 A Magnetic field at the center of loop = B (suppose) magnetic field at distance x from the center = `\frac{1}{8}B`   To find, x = ? The magnetic field at distance x from center of the coil is given as, `B = \frac{\mu_o I N r^2}{2(x^2  + r^2)^(\frac{3}{2})}`. By given in question, magnetic field at x = `\frac{1}{8} \times` magnetic field at x = 0.       or, `\frac{\mu_o I N r^2}{2(x^2 + r^2)^(\frac{3}{2})} =\frac{1}{8} \frac{\mu_o I N r^2}{2(0^2 +
Read More »

how to create GIF animation in gnuplot | animation using gnuplot software

-
Image
 GNUPLOT software is not only 2D and/or 3D plotting tool, but it also provides flexibility to create animation and GIF animation plots.  GIF means Graphics Interchange Format. It is a repeated display of images from same of different sides. It helps you to view the image clearly so that you don't need to pay much attention to understand it.  This post provides you tools and techniques to create GIF animation in gnuplot. Codes necessary for creating the GIF are given in links. Please check it as well. GIF can also be both 2D and 3D.  You must use 'plot' command to create 2D GIF and 'splot' command to create 3D GIF. 👉 To learn to create plots in gnuplot: Click Here   👉 To learn gnuplot basics for beginners: Click Here Example of 2D GIF: GNUPLOT code for creating above GIF, 👇 gnuplot > set title 'gif image of sine and cosine functions' gnuplot > set xlabel 'x' gnuplot > set ylabel 'sine and cosine' gnuplot > f1(x,t) = sin(x-t) gnu
Read More »

how to plot 3D and parametric graphs in gnuplot | three dimensional plot in gnuplot software

-
Image
Although being an open-source plotting software, gnuplot gives you almost all feasibilities to plot 2D, 3D, parametric plots and even the gif as well as animations.  Here you can learn to plot three dimensional graphs in gnuplot. Additionally, you also get the idea of parametric plots in gnuplot software. 1. 3D graphs in gnuplot To plot 3D graphs you just need to replace the command 'plot' for 2D graphs by 'splot' for 3D graphs.  Here, splot means surface plot.  Plotting graphs in 3D view allows you to observe it from all possible sides of the graph.  You rotate it upward, downward , left-side as well as to right-side.  Example: gnuplot >  splot sin(x) Then press ENTER key. You get the output as it looks as picture below: Fig. 1 3D graph of sin(x). 2. parametric plots Parametric plots can be created both in 2D and 3D.  For 2D graph, the dummy variable is t, whereas for 3D plot, they are u and v. Example: Image of parametric plot is shown below: Fig.2 Parametric plot
Read More »

NEB board exam maths solved problems with proper solutions - available for free | short answered mathematics problems solved group A

-
Image
 This post provides you the problems with solutions of NEB board exam of grade 11 mathematics 2076. Attempt All Questions. Group A  `5 \times 3 \times 2 = 30` Get full question paper from Here . 1. a. Find the truth value of the biconditional statements of p and q where p represents '2 is an even number' and q represents '4 is an even number'. Solution:- Biconditional statements imply if one statement is true then another is also true. if one is false, then another statement is also false. So, one statement supports the other in case of biconditional statement.  Here, p: 2 is an even number.          q: 4 is an even number. 4 is the square of 2 and 2 is the square root of 4. 2 is even, then 4 is also even. So, both p and q represent true value (T). p q p ⇔ q T T T Thus truth value of p and q is T, i.e. 2 is an even number implies 4 is also an even number.// b. If f
Read More »

NEB grade 11 and 12 maths exam MCQ solved problems | NEB mathematics MCQ solved problems

-
Image
 NEB  Mathematics grade 11 and 12 specification grid solutions Short Solution Problems(Group A) Rewrite the correct option in your answer sheet. 1. Which of the following is a statement? a) The fishes are beautiful. b) Study Mathematics c) x is the capital of country y. d) Water is essential for health. Correct option: (d) It because any sentence is considered as a statement if it is either true or false but not both. Among given options, (a), (b) and (c)  do not declare truth or falsity. But option (d) declares the truth that water is important for health. 2. The  value of: `\sqrt{-16} \times \sqrt{-25}` is,  (a) -20 (b) -20i  (c) 20i  (d)20 Correct option: (a) Here,  `\sqrt{-16} \times \sqrt{-25}` = `\sqrt{-1 \times 16} \times \sqrt{-1 \times 25}`  But `i^2 = -1`, i = complex number. = `\sqrt{i^2 \times 4^2} \times \sqrt{i^2 \times 5^2}` = ` i \times 4 \times i \times 5` = `i^2 \times 20` = -20. // Click here for NEB grade 11 maths short solution problems(Group A)   3.  If `\angle C
Read More »