In [1]:
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
In [2]:
import libr3
In [3]:
# Pressure in [Pa]
bar = 1e5
P = 0.1*bar
gen = libr3.GeneratorLiBr(P,1,60,0.5,0.75)
In [4]:
gen
Out[4]:
In [5]:
gen.q(91.9642)
Out[5]:
In [6]:
T = np.linspace(gen.T_in,gen.T_out)
plt.plot(gen.q(T),T,'-o')
plt.xlabel('q (W)')
plt.ylabel('T ($^\circ$C)')
Out[6]:
In [7]:
from HRHX_integral_model import *
useHotT = True
cold = streamExample1(0)
hot = streamExample1(1,-1)
ci = counterflow_integrator(cold, hot, useHotT)
Qmax = ci.Qmax
print "Qmax=",Qmax
UA=ci.calcUA(Qmax/2)
print "UA=",UA
plotFlow(ci)
plotNN(ci)
c2 = streamExample2(-5.,100.,1.,1.,10.)
h2 = streamExample1(120.,-1.5)
ci2 = counterflow_integrator(c2,h2, useHotT)
print ci2.Qmax
plotFlow(ci2)
plotNN(ci2)
plt.show()
c3 = streamExample2(-5.,100.,1.,1.,10.)
h3 = streamExample2(15.,100.,-1.,1.,10.)
ci3 = counterflow_integrator(c3,h3, useHotT)
print ci3.Qmax
plotFlow(ci3)
plotNN(ci3)
plt.show()
QQ=np.arange(0,35)
Tc=c3.T(QQ)
Th=h3.T(QQ)
TT=np.arange(95,110)
Qc=c3.q(TT)
Qh=h3.q(TT)
plt.figure();plt.plot(QQ,Tc,QQ,Th);plt.plot(Qc,TT,'.',Qh,TT,'.')
Out[7]:
In [8]:
# Pressure in [Pa]
bar = 1e5
P = 0.1*bar
m_in = 0.2 #kg/s
T_in = 60 # deg C
x_in,x_out = 0.5,0.69
gen2 = libr3.GeneratorLiBr(P,m_in,T_in,x_in,x_out)
print gen2
print gen2.Tmax
In [9]:
#%%prun
T = np.linspace(gen2.T_in,gen2.T_out)
Qc = gen2.q(T)
QQ = np.linspace(0,2e5)
Tc = gen2.T(QQ)
plt.figure()
plt.plot(QQ,Tc)
plt.plot(Qc,T,'.')
plt.xlabel('q (W)')
plt.ylabel('T ($^\circ$C)')
Out[9]:
In [10]:
#%%prun
genhot = streamExample1(130,-10000,1)
#cigen = counterflow_integrator(gen2,genhot, False)
""" 1446681 function calls (1445561 primitive calls) in 25.872 seconds
Ordered by: internal time
ncalls tottime percall cumtime percall filename:lineno(function)
46222 9.878 0.000 10.462 0.000 libr_props.py:53(thetaFun)
27880 9.728 0.000 9.930 0.000 {CoolProp.CoolProp.PropsSI}
5576 2.255 0.000 9.536 0.002 libr_props.py:187(massSpecificEnthalp"""
cigen = counterflow_integrator(gen2,genhot, True)
93629 function calls (93497 primitive calls) in 1.483 seconds
Ordered by: internal time
ncalls tottime percall cumtime percall filename:lineno(function)
2982 0.576 0.000 0.616 0.000 libr_props.py:53(thetaFun)
1720 0.569 0.000 0.569 0.000 {CoolProp.CoolProp.PropsSI}
344 0.074 0.000 0.493 0.001 libr_props.py:187(massSpecificEnthalpy)
In [14]:
# Warning: this calls inverse function, T(q), and is very slow.
plotFlow(cigen)
plt.legend(loc='best')
Out[14]:
In [16]:
# Pressure in [Pa]
bar = 1e5
P = 0.1*bar
m_in = 0.2 #kg/s
T_in = 60 # deg C
x_in,x_out = 0.5,0.69
gen3 = libr3.GeneratorLiBrInterpolated(P,m_in,T_in,x_in,x_out)
print gen3
print gen3.q
print gen3.q(115), gen3.q(116)
In [15]:
#%%prun
T = np.linspace(gen3.T_in,gen3.Tmax+10)
Qc = gen3.q(T)
QQ = np.linspace(0,gen3._q(gen3.Tmax)*1.2)
Tc = gen3.T(QQ)
plt.figure()
plt.plot(QQ,Tc)
plt.plot(Qc,T,'.')
plt.xlabel('q (W)')
plt.ylabel('T ($^\circ$C)')
print gen2.Tmax
In [13]:
import tabulate
print gen3.Tmax, gen3._q(gen3.Tmax), gen3.q(gen3.Tmax)
print tabulate.tabulate(zip(T,Qc,Tc,QQ))
In [14]:
%%prun
#cigen = counterflow_integrator(gen2,genhot, False)
cigen = counterflow_integrator(gen3,genhot, True)
In [20]:
%
In [ ]: