
-²lZ˜  ã               @   s‰  d  Z  d d l Z d d l j Z d Z d d d d „ Z d d d d	 d
 d „ Z e j	 d d d d d ƒ\ Z
 Z d d „  e
 Dƒ Z d d „  e
 Dƒ Z d Z d Z e j d e e f ƒ e j e
 e d d d d d d d d ƒe j e
 e d d d d d d d d ƒe j e
 e d d  d d! d d d" d# d d$ ƒe j d% d& d' ƒe j d( d& d' ƒe j d) d* ƒ e j d+ ƒ e j ƒ  d S),a¨  
@author:       Ole Havik Bjoerkedal
@copyright:    Ole Havik Bjoerkedal
@contact:      ole.h.bjorkedal@ntnu.no
@organization: Department of Chemical Engineering, NTNU, Norway
@since:        2018-01-19
@purpose:      Calculating and visualizing equations of state for gases
               Ideal Gas and van der Waal
               Support material for Exercise 3 in TKP4120
@requires:     Python 3
               matplotlib
é    NgTã¥›Ä  @g      r@g      ð?c             C   s   | t  | |  } | S)zå
    @parameters:
    -------------
    V = Volume [m3]
    T = Temperature [K], default value is 298 K (25 C)
    n = Amount of gas [mol], default value is 1 mol

    @returns:
    -------------
    Ideal pressure [Pa]     
    )ÚR)ÚVÚTÚnÚP© r   úE/Users/vegardjervell/Documents/2semester/Prossess/py_oving/Ex3_vdw.pyÚIdealGasPressure   s    r	   g.ÿ!ýöuñ?gT?Ì¹o?c             C   s,   |  | } t  | | | | | d } | S)aA  
    @parameters:
    -------------
    V = Volume [m3]
    T = Temperature [K], default value is 298 K (25 C)
    n = Amount of gas [mol], default value is 1 mol
    a = van der Waal parameter [J m3 mol-2]
    b = van der Waal parameter [m3 mol-1]
    @returns:
    -------------
    van der Waal pressure [Pa]     
    é   )r   )r   r   r   ÚaÚbZVmr   r   r   r   ÚVanDerWaalPressure!   s    
r   zpropaneData_2.txtZcommentsú#ZunpackTc             C   s   g  |  ] } t  | ƒ ‘ q Sr   )r	   )Ú.0r   r   r   r   ú
<listcomp>B   s   	 r   c             C   s   g  |  ] } t  | ƒ ‘ q Sr   )r   )r   r   r   r   r   r   C   s   	 g      4@ZfigsizeZcolorr   Z	linestyleú-Zlwg      à?Zlabelz	Ideal GasÚrzvan der Waal GasÚgÚNoneZmarkerÚxzReal GaszVolume [m3]Zfontsizeé   zPressure [Pa]Úloczupper rightzpressure.pdf)Ú__doc__ZnumpyZnpZmatplotlib.pyplotZpyplotZpltr   r	   r   ZloadtxtZvolumeZrealPressureZidealPressureZvdwPressureÚwidthZheightZfigureZsemilogxZxlabelZylabelZlegendZsavefigZshowr   r   r   r   Ú<module>   s&   !((.