.. examples

Examples
========

State Point
-----------

Create a listing of properties at a given T and D::
    
    from engcoolprop.ec_fluid import EC_Fluid
    
    # Print state point
    ec = EC_Fluid(symbol="N2", T=530.0,P=100.0 ) # T=degR, P=psia
    
    ec.setProps(T=500., D=0.1) # Set T and D

    ec.printProps() # Print state point at given T,P

Ouput::
    
    State Point for fluid Nitrogen (N2)
    T =     500  degR (Tc= 227.146 , Tnbp= 139.239 Ttriple= 113.672 )
    P = 19.1444  psia (Pc= 492.519 )
    D =     0.1  lbm/cu ft (Dc= 19.5587 )
    E = 88.3793  BTU/lbm
    H = 123.806  BTU/lbm
    S = 1.59617  BTU/lbm degR
    Cv=0.177465  BTU/lbm degR
    Cp=0.248873  BTU/lbm degR
    g = 1.40238  Cp/Cv (-)
    A =    1115  ft/sec
    V = 1.13254  viscosity [1.0E5 * lb/ft-sec]
    C =0.0140818  thermal conductivity [BTU/ft-hr-R]
    MW= 28.0135  lbm/lbmmole
    Q =      -1  Vapor Quality (mass fraction gas)
    Z =0.999722  (-)    
    

Critical Properties
-------------------

Create a listing of critical properties::

    from engcoolprop.ec_fluid import EC_Fluid

    ec = EC_Fluid(symbol="O2", T=530.0,P=100.0 ) # T=degR, P=psia

    ec.printCriticalProps() # Print critical properties

Output::

    Critical Properties for fluid Oxygen (O2)
    Tc= 278.246 degR
    Tnbp= 162.338 degR
    Pc = 731.425 psia
    Dc = 27.2276 lbm/cu ft
    MW= 31.9988  lbm/lbmmole
    Zc =0.287943  (-)
    Hc = 13.9393  BTU/lbm
    Ec = 8.96846  BTU/lbm
    Sc = 1.00335  BTU/lbm degR
    Ac = 502.502  ft/sec
    Vc =     nan  viscosity [1.0E5 * lb/ft-sec]
    Cc =0.0221309  thermal conductivity [BTU/ft-hr-R]


Table of Values
---------------

Create a table of properties::

    from engcoolprop.ec_fluid import EC_Fluid

    ec = EC_Fluid(symbol="AR")
    print('Properties for: %s (%s)'%(ec.name, ec.symbol) )
    for T in range(460, 570, 10):
        ec.setProps(T=T, P=14.7)
        ec.printTPD()

Output::

    Properties for: Argon (AR)
    AR T= 460.0 P=  14.7 D=0.1191 E= 34.21 H= 57.05 S=0.906 Q=-1.00
    AR T= 470.0 P=  14.7 D=0.1166 E= 34.96 H= 58.30 S=0.908 Q=-1.00
    AR T= 480.0 P=  14.7 D=0.1141 E= 35.71 H= 59.55 S=0.911 Q=-1.00
    AR T= 490.0 P=  14.7 D=0.1118 E= 36.46 H= 60.79 S=0.914 Q=-1.00
    AR T= 500.0 P=  14.7 D=0.1095 E= 37.21 H= 62.04 S=0.916 Q=-1.00
    AR T= 510.0 P=  14.7 D=0.1074 E= 37.95 H= 63.29 S=0.919 Q=-1.00
    AR T= 520.0 P=  14.7 D=0.1053 E= 38.70 H= 64.53 S=0.921 Q=-1.00
    AR T= 530.0 P=  14.7 D=0.1033 E= 39.45 H= 65.78 S=0.923 Q=-1.00
    AR T= 540.0 P=  14.7 D=0.1014 E= 40.20 H= 67.02 S=0.926 Q=-1.00
    AR T= 550.0 P=  14.7 D=0.0995 E= 40.94 H= 68.27 S=0.928 Q=-1.00
    AR T= 560.0 P=  14.7 D=0.0978 E= 41.69 H= 69.52 S=0.930 Q=-1.00

TD Plot
-------

Create a TD plot using CoolProp plotting routine::

    from CoolProp.Plots.Common import BasePlot, PropertyDict, BaseDimension
    from engcoolprop.ec_fluid import EC_Fluid

    symbol = 'N2'

    class ENGunits(PropertyDict):
        def __init__(self):
            self._D = BaseDimension(add_SI=0.0, mul_SI=1./16.01843417, off_SI=0.0, label='Density',       symbol=u'd', unit=u'lbm/ft3')
            self._H = BaseDimension(add_SI=0.0, mul_SI=0.000429923, off_SI=0.0, label='Specific Enthalpy',symbol=u'h', unit=u'BTU/lbm')
            self._P = BaseDimension(add_SI=0.0, mul_SI=1.0/6894.76, off_SI=0.0, label='Pressure',         symbol=u'p', unit=u'psia')
            self._S = BaseDimension(add_SI=0.0, mul_SI=0.0002388461,off_SI=0.0, label='Specific Entropy', symbol=u's', unit=u'BTU/lbm/R')
            self._T = BaseDimension(add_SI=0.0, mul_SI=1.8,         off_SI=0.0, label='Temperature',      symbol=u'T', unit=u'R')
            self._U = BaseDimension(add_SI=0.0, mul_SI=0.000429923, off_SI=0.0, label='Specific Internal Energy', symbol=u'u', unit=u'BTU/lbm')
            self._Q = BaseDimension(add_SI=0.0, mul_SI=1.0,         off_SI=0.0, label='Vapour Quality',   symbol=u'x', unit=u'')
            
    BasePlot.UNIT_SYSTEMS['ENG'] = ENGunits() # ['EUR','KSI','SI']

    ec = EC_Fluid(symbol)
    ec.printCriticalProps()

    # Plot Types: 'PS', 'PT', 'HS', 'TS', 'PD', 'TD', 'PH'
    # tp_limits:  'NONE', 'DEF','ACHP','ORC'
    plot = PropertyPlot(symbol, 'TD', unit_system='ENG', tp_limits='DEF')
    plot.calc_isolines(CP.iHmass, num=11)
    plot.calc_isolines(CP.iQ, num=11)
    plot.title('%s (%s)'%(ec.name, ec.symbol))
    plot.grid()
    plot.show()

.. image:: _static/N2_TD_plot.jpg
    :width: 60%

And printed Output::

    Critical Properties for fluid Nitrogen (N2)
    Tc= 227.146 degR
    Tnbp= 139.239 degR
    Pc = 492.519 psia
    Dc = 19.5587 lbm/cu ft
    MW= 28.0135  lbm/lbmmole
    Zc =0.289458  (-)
    Hc = 12.5678  BTU/lbm
    Ec = 7.90791  BTU/lbm
    Sc = 1.00671  BTU/lbm degR
    Ac = 444.365  ft/sec
    Vc = 1.22997  viscosity [1.0E5 * lb/ft-sec]
    Cc = 886.461  thermal conductivity [BTU/ft-hr-R]