Python steam_ph Examples

Example #1

# these are the pressures we're interested in here
pp = [freesteam.psat_T(50 + 273.15), 3e6, 12e6, 165e5, 300e5]

print("low p =", pp[0])

hh = arange(50., 4500., 100) * 1e3

hh1 = arange(50., 4500., 20) * 1e3

x = []
y = []
u = []
v = []
for p in pp:
    plot(hh1 / 1e3, [freesteam.steam_ph(p, h).T for h in hh1], alpha=0.8)
    if 0:
        for h in hh:
            try:
                S = freesteam.steam_ph(p, h)
                x += [S.h / 1.e3]
                y += [S.T]
                dy = freesteam.freesteam_deriv(S, 'T', 'h', 'p')
                dx = 0.0005
                m = math.sqrt(dx**2 + dy**2)
                u += [dx / m]
                v += [dy / m]
            except:
                pass

plot([0, 4500], [550 + 273.15, 550 + 273.15], 'b--')

Example #2

File: Thdiagram.py Project: zcarwile/conda_build_demo

# these are the pressures we're interested in here
pp = [freesteam.psat_T(50+273.15), 3e6, 12e6,  165e5, 300e5]

print "low p =",pp[0]

hh = arange(50.,4500.,100)*1e3

hh1 = arange(50.,4500.,20)*1e3


x = []
y = []
u = []
v = []
for p in pp:
	plot(hh1/1e3,[freesteam.steam_ph(p,h).T for h in hh1],alpha=0.8)
	if 0:
		for h in hh:
			try:
				S = freesteam.steam_ph(p,h)
				x += [S.h/1.e3]
				y += [S.T]
				dy = freesteam.freesteam_deriv(S,'T','h','p')
				dx = 0.0005
				m = math.sqrt(dx**2 + dy**2)
				u += [dx/m]
				v += [dy/m]
			except:
				pass

plot([0,4500],[550+273.15,550+273.15],'b--')

Example #3

File: pTdiagram.py Project: qingfengxia/freesteam

	for T in TT:
		S = freesteam.steam_pT(p,T)
		#print "p = %f, T = %f" % (p,T)
		r = S.region
		#print "p = %f MPa, T = %f K, region[%d,%d] = %d" % (p/1e6,T,x,y,r)
		im[x,y] = float(r) / 4.
		y += 1
	x += 1

imshow(im,extent=[Tmin,Tmax,pmin/1e6,pmax/1e6],origin='lower',aspect='auto',interpolation='nearest',alpha=0.6)

# LINES OF CONSTANT ENTHALPY

hh = arange(100,4500,200)*1e3
for h in hh:
	TT2 = [freesteam.steam_ph(p,h).T for p in pp]
	plot(TT2,pp/1e6,'g-')

plot(TT0,psat,'b-')
axis([Tmin,Tmax,pmin/1e6,pmax/1e6])
#quiver(x,y,u,v,alpha=0.6)

# LINES OF CONSTANT ENTROPY

ss = arange(0,12,1./3)*1e3
for s in ss:
	TT2 = [freesteam.steam_ps(p,s).T for p in pp]
	plot(TT2,pp/1e6,'r:')

plot(TT0,psat,'b-')
axis([Tmin,Tmax,pmin/1e6,pmax/1e6])

Example #4

#!/usr/bin/env python

import freesteam

S = freesteam.steam_ph(100e5, 300)

print("TESTING RESULTS")

print("region =", S.region)
print("h =", S.h)
print("v =", S.v)
print("p =", S.p)
print("s =", S.s)
print("mu=", S.mu)

print("TESTING PV")

p = 500e5
v = 1. / 401.
print("(p,v) = (%f, %f)" % (p, v))
print("bounds errors?", freesteam.bounds_pv(p, v, 1))
print("region?", freesteam.region_pv(p, v))

S = freesteam.steam_pv(p, v)

print("region =", S.region)
print("h =", S.h)
print("v =", S.v)
print("p =", S.p)
print("s =", S.s)
print("mu=", S.mu)

Example #5

File: test.py Project: zcarwile/conda_build_demo

		sys.stderr.write("can do this with 'export PATH=..:$PATH', assuming.\n")
		sys.stderr.write("you are in the directory containing test.py.\n")
		sys.exit(1)

else:
	if not os.environ.get(LDVAR) or up not in os.environ[LDVAR]:
		os.environ[LDVAR] = up
		script = os.path.join(sys.path[0],"test.py")					
		os.execvp("python",[script] + sys.argv)

#-------------------------------------------------------------------------------
# now for the tests...

import freesteam

S = freesteam.steam_ph(100e5,300)

print "TESTING RESULTS"

print "region =",S.region
print "h =",S.h
print "v =",S.v
print "p =",S.p
print "s =",S.s
print "mu=",S.mu

print "TESTING PV"

p = 500e5
v = 1./401.
print "(p,v) = (%f, %f)" % (p, v)