Solving Van Der Waals gas equation for volume of a real gas by the method of successive approximation.

 To solve van der Waals gas equation for volume of a real gas by the method of successive approximation.







Python Code:


# ------------------------------------------------
# Exp -4: To solve van der Waals gas equation for
# volume of a real gas by the method of
# successive approximation ( iteration method ).
# ------------------------------------------------
# V = R*T/(P + a/(V*V)) + b = g(V)
# range of V = [3 , 10]
# a = 1.5
# T = 300 K
# R = 8.3144 J/(K. mol )
# -----------------------------------------------

import numpy as np
import math

a = 1.5 # constant
T = 300 #K
R = 8.3144 # J/(K. mol )

# --------------------------
P = 40 # variable
b = 0; # vaiiable
# --------------------------

V = 10.0 # initial guass

for it in range (1 , 101 , 1) : # number of iteration

  V_new = R*T /( P + a /( V*V )) + b

  if ( abs ( V_new - V) < 0.008322) :
    break
  V = V_new

# ----------------------------------------
print (" root of VDW equation : ", V_new )
print (" number of iteration : ", it )

# ----------------------------------------

Comments

Post a Comment

Popular posts from this blog

Solving Kepler’s equation by Newton-Raphson method

Image
 To solve Kepler’s equation by Newton-Raphson method Python Program: # ------------------------------------------------------------- # Exp -3: To solve K e p l e r s equation by Newton - Raphson method # ------------------------------------------------------------- import  numpy  as  np import  math # ------------------------------------------------------------- def   kepler_eq  ( M  ,  E  ,  e  ,  err  ): # M [ rad ], E [ rad ], e = between 0 to 1, E inital assumption   print ( " --------------------------------------------------------------" )   print ( " value of M:" , M ,  ", value of e:" ,e ,  ", intial value of E:" ,E...
Read more

Velapanti Python Program

 print("Hello!!") print("Welcome to Physics Tech Talks") print("What is your name ? ") name=str(input()) print("Are u a male(1) or female(2)? ") print("Enter 1 if male and 2 if female: ") sex=int(input()) if sex==1:     print("Hlo Sir, this is soo nice talking to u !!!")     print("Mr.",name,"this program is for time pass !!! Do u know that?? If yes press 1 if no then 0:  ") else:      print("Hlo mam, this is soo nice talking with u !!!")     S=int(input("Are u married (Yes(1)or No(0)): "))     if S==1:         print("Mrs.",name,"this program is for time pass !!! Do u know that?? If yes press 1 if no then 0: ")     if S==0:         print("Miss",name,"this program is for time pass !!! Do u know that?? If yes press 1 if no then 0: ")     k=int(input()) if k==1:     print("I am glad to know that you are aware of your velapanti. Yo know that you a...
Read more

To find mean , standard deviation and frequency distribution of an actual data set (Python Program)

Image
  Python Program: # --------------------------------------------------------------------- # Exp -1: To find mean , standard deviation and frequency distribution # of an actual data set from any physics experiment # --------------------------------------------------------------------- import  numpy  as  np  # numpy related numerical python import  statistics import  matplotlib . pyplot  as  plt x = [ 102  ,  104  ,  107  ,  112  ,  108  ,  110  ,  103  ,  101  ,  100  ,  104  ,  103  ,  104  ,  101  ,  105 ]  print ( "-------------------------------------------------------" ) print ( "             Re...
Read more