def Expander():
Ref = 'Nitrogen'
Tinlet = 300
pinlet = 2000
poutlet = 500
SE = ScrollExpander()
SE.set_scroll_geo(100e-6,2.3,0.005,0.005)
SE.set_disc_geo('2Arc', r2 = 0)
SE.geo.delta_suction_offset = 0.0e-3
SE.geo.phi_ie_offset = 0.0
SE.omega = 60*2*pi
SE.Tamb = 298.0
#Temporarily set the bearing dimensions
SE.mech = core.struct()
SE.mech.D_upper_bearing = 0.02540
SE.mech.L_upper_bearing = 0.031
SE.mech.c_upper_bearing = 10e-6 #??
SE.mech.D_crank_bearing = 0.0254
SE.mech.L_crank_bearing = 0.023
SE.mech.c_crank_bearing = 10e-6 #??
SE.mech.D_lower_bearing = 0.01895
SE.mech.L_lower_bearing = 0.023
SE.mech.c_lower_bearing = 10e-6 #??
SE.mech.thrust_ID = 0.05
SE.mech.thrust_friction_coefficient = 0.02 #Tuned for point A1
SE.mech.L_ratio_bearings = 5
SE.mech.mu_oil = 0.008
SE.mech.scroll_plate_thickness = 0.007
SE.mech.scroll_density = 8100
SE.mech.scroll_added_mass = 0
SE.mech.scroll_plate_diameter = 0.092 #[m]
m, zcm = SE.calculate_scroll_mass()
SE.mech.orbiting_scroll_mass = m
SE.mech.scroll_zcm__thrust_surface = zcm
SE.mech.detailed_analysis = False
SE.h_shell = 10
SE.A_shell = 0.05
SE.HTC = 0.5
SE.generator = Motor()
SE.generator.set_eta(1.0)
SE.generator.suction_fraction = 0.0
inletState = CPState(Ref,{'T':Tinlet,'P':pinlet})
T2s = SE.guess_outlet_temp(inletState,poutlet)
outletState = CPState(Ref,{'T':T2s,'P':poutlet})
SE.geo.delta_radial = 0e-6
SE.geo.delta_flank = 0e-6
mdot_guess = inletState.rho*SE.Vdisp/SE.Vratio*SE.omega/(2*pi)
SE.add_tube(Tube(key1 = 'inlet.1',
key2 = 'inlet.2',
L = 0.03,
ID = 0.02,
mdot = mdot_guess,
State1 = inletState.copy(),
fixed = 1,
TubeFcn = SE.TubeCode))
SE.add_tube(Tube(key1 = 'outlet.1',
key2 = 'outlet.2',
L = 0.03,
ID = 0.02,
mdot = mdot_guess,
State2 = outletState.copy(),
fixed = 2,
TubeFcn = SE.TubeCode))
SE.auto_add_CVs(inletState, outletState)
# SE.auto_add_leakage(flankFunc = SE.FlankLeakage,
# radialFunc = SE.RadialLeakage)
FP = FlowPath(key1='outlet.1',
key2='da',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.02**2/4
SE.add_flow(FP)
SE.add_flow(FlowPath(key1 = 'da',
key2 = 'd1',
MdotFcn = SE.DA_D,
MdotFcn_kwargs = dict(X_d = 0.5)
)
)
SE.add_flow(FlowPath(key1 = 'da',
key2 = 'd2',
MdotFcn = SE.DA_D,
MdotFcn_kwargs = dict(X_d = 0.5)
)
)
FP = FlowPath(key1='inlet.2',
key2='ss',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.007**2/4
SE.add_flow(FP)
FP = FlowPath(key1='inlet.2',
key2='sss',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.007**2/4
SE.add_flow(FP)
SE.add_flow(FlowPath(key1='s1',
key2='ss',
MdotFcn=SE.S_to_SS))
SE.add_flow(FlowPath(key1='s2',
key2='ss',
MdotFcn=SE.S_to_SS))
SE.connect_callbacks(step_callback = SE.step_callback,
lumps_energy_balance_callback = SE.lump_energy_balance_callback,
endcycle_callback= SE.endcycle_callback)
# beta = np.linspace(0,2*pi,1000)
# V_e11 = np.array([SE.V_e1(b,1)[0] for b in beta])
# V_e21 = np.array([SE.V_e2(b,1)[0] for b in beta])
# V_e12 = np.array([SE.V_e1(b,2)[0] for b in beta])
# V_e22 = np.array([SE.V_e2(b,2)[0] for b in beta])
# V_d1 = np.array([SE.V_d1(b)[0] for b in beta])
# V_d2 = np.array([SE.V_d2(b)[0] for b in beta])
# V_s1 = np.array([SE.V_s1(b)[0] for b in beta])
#
# plt.plot(beta,V_e11)
# plt.plot(beta,V_e12)
# plt.plot(beta,V_d1)
# plt.plot(beta,V_s1)
# plt.gca().axvline(2*pi-SE.theta_d)
# plt.show()
#try:
SE.eps_cycle = 0.003
SE.precond_solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
UseNR = False, #If True, use Newton-Raphson ND solver to determine the initial state
OneCycle = False,
plot_every_cycle= False,
hmin = 1e-8
)
def Compressor(Te=0,
DTsh=11.1,
Tc=20,
Tamb=25,
Nmot=3600,
f=None,
OneCycle=False,
Ref='R134a',
HDF5file='scroll_compressor.h5'):
ScrollComp = Scroll()
#Set the scroll wrap geometry
ScrollComp.set_scroll_geo(82.7e-6,
2.5,
0.004,
0.005,
phi_i0=0.0,
phi_os=0.3,
phi_is=3.142)
ScrollComp.set_disc_geo('2Arc', r2=0.0010)
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = Nmot / 60 * 2 * pi
ScrollComp.slip_ratio = 0.01
ScrollComp.Tamb = Tamb + 273.15 #[K]
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.detailed_analysis = True
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.013 #0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.005
ScrollComp.mech.rho_oil = 930
ScrollComp.mech.oldham_ring_radius = 0.06 #m
ScrollComp.mech.oldham_mass = 0.1 #kg
ScrollComp.mech.oldham_thickness = 0.008 #m
ScrollComp.mech.oldham_key_height = 0.006 #m
ScrollComp.mech.oldham_key_width = 0.006 #m
ScrollComp.mech.oldham_key_friction_coefficient = 0.01 #-
ScrollComp.mech.oldham_rotation_beta = 0 #rad
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 450.75e-3 * ((246.126e-3)**2 * pi / 4)
ScrollComp.HTC = 0.0
ScrollComp.HT_corr = 'Pereira-Deschamps' #'Jang-Jeong'
# Temperature Lumps
ScrollComp.OEB_type = 'multi-lump' #'single-lump'
ScrollComp.OEB_solver = 'MDNR'
ScrollComp.Rshell_oil = 190 #K/kW from Chen (2000) - PhD thesis
# Define motor efficiency
ScrollComp.motor.set_eta(0.95)
ScrollComp.motor.suction_fraction = 0.5
# Operating conditions
Tevap = Te + 273.15 #K
Tcond = Tc + 273.15 #K
Tin = Tevap + DTsh
DT_sc = 7 #K
temp = State.State(Ref, {'T': Tevap, 'Q': 1})
pe = temp.p
temp.update(dict(T=Tcond, Q=1))
pc = temp.p
inletState = State.State(Ref, {'T': Tin, 'P': pe})
ScrollComp.pressure_ratio = pc / pe
print('Pressure Ratio:', pc / pe, pc, pe)
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
print('Vdisp:', ScrollComp.Vdisp * 1e6, 'cm3/rev')
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
#Leakage flow
ScrollComp.geo.delta_flank = -9.615e-7 * (ScrollComp.pressure_ratio -
1.67) + 15e-6 #10e-6
ScrollComp.geo.delta_radial = 1.1e-6 * (ScrollComp.pressure_ratio -
1.67) + 1e-6 #10e-6
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.8)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.8)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.8)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.8)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='d1', key2='dd', MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(
FlowPath(key1='d2', key2='dd', MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
from time import clock
t1 = clock()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
ScrollComp.verbosity = 10
ScrollComp.precond_solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
x0=[330, 330, 350], #Guesses [Td,Tlump[0],Tlump[1]]
#hmin = 1e-3
eps_cycle=3e-3,
eps_energy_balance=0.1 #relaxed multi-lump convergence
)
print('time taken', clock() - t1)
#debug_plots(ScrollComp)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
return ScrollComp
def Compressor(Te = 273, Tc = 300, f = None,TTSE = False, OneCycle = False):
if TTSE:
CP.set_TTSESinglePhase_LUT_size("Propane", 500, 500)
CP.enable_TTSE_LUT('Propane')
global Injection
ScrollComp=Scroll()
#This runs if the module code is run directly
ScrollComp.set_scroll_geo(83e-6, 3.3, 0.005, 0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc',r2 = 0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = 3000/60*2*pi
ScrollComp.Tamb = 298.0
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 0.05
ScrollComp.HTC = 1.0
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Ref = 'Propane'
#Ref = 'REFPROP-MIX:R410A.mix'
Te = -20 + 273.15
Tc = 20 + 273.15
Tin = Te + 11.1
DT_sc = 7
pe = CP.PropsSI('P','T',Te,'Q',1.0,Ref)/1000.0
pc = CP.PropsSI('P','T',Tc,'Q',1.0,Ref)/1000.0
inletState = State.State(Ref,{'T':Tin,'P':pe})
T2s = ScrollComp.guess_outlet_temp(inletState,pc)
outletState = State.State(Ref,{'T':T2s,'P':pc})
mdot_guess = inletState.rho*ScrollComp.Vdisp*ScrollComp.omega/(2*pi)
ScrollComp.add_tube(Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc = ScrollComp.FlankLeakage,
radialFunc = ScrollComp.RadialLeakage)
FP = FlowPath(key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.01**2/4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs = dict(X_d = 0.7)
)
)
ScrollComp.add_flow(FlowPath(key1 = 'sa',
key2 = 's2',
MdotFcn = ScrollComp.SA_S2,
MdotFcn_kwargs = dict(X_d = 0.7)
)
)
ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
key2 = 'dd',
MdotFcn = ScrollComp.DISC_DD,
MdotFcn_kwargs = dict(X_d = 0.7)
)
)
ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
key2 = 'ddd',
MdotFcn = ScrollComp.DISC_DD,
MdotFcn_kwargs = dict(X_d = 0.7)
)
)
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
ScrollComp.add_flow(FlowPath(key1='d1',
key2='dd',
MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(FlowPath(key1='d2',
key2='dd',
MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(step_callback = ScrollComp.step_callback,
endcycle_callback = ScrollComp.endcycle_callback,
heat_transfer_callback = ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback = ScrollComp.lump_energy_balance_callback
)
from time import clock
t1=clock()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.precond_solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle = OneCycle,
plot_every_cycle= False,
#hmin = 1e-3
eps_cycle = 3e-3
)
except BaseException as E:
print(E)
raise
print 'time taken',clock()-t1
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
import CoolProp
h5.write_to_file(ScrollComp, 'CPgit_'+CoolProp.__gitrevision__+'.h5')
return ScrollComp
def Compressor(ScrollClass,
Te=273,
Tc=300,
f=None,
OneCycle=False,
Ref='R410A',
HDF5file='scroll_compressor.h5'):
ScrollComp = ScrollClass()
#This runs if the module code is run directly
ScrollComp.set_scroll_geo(83e-6, 3.3, 0.005,
0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = 3000 / 60 * 2 * pi
ScrollComp.Tamb = 298.0
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 0.05
ScrollComp.HTC = 0.0
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Te = -20 + 273.15
Tc = 20 + 273.15
Tin = Te + 11.1
DT_sc = 7
temp = State.State(Ref, {'T': Te, 'Q': 1})
pe = temp.p
temp.update(dict(T=Tc, Q=1))
pc = temp.p
inletState = State.State(Ref, {'T': Tin, 'P': pe})
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='d1', key2='dd', MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(
FlowPath(key1='d2', key2='dd', MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
t1 = timeit.default_timer()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.precond_solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
#hmin = 1e-3
eps_cycle=3e-3)
except BaseException as E:
print(E)
raise
print('time taken', timeit.default_timer() - t1)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
# debug_plots(ScrollComp, family='Scroll Compressor')
return ScrollComp
def Compressor(ScrollClass,
Te=253,
Tc=310,
f=None,
OneCycle=False,
Ref='R410A',
HDF5file='scroll_compressor.h5',
discharge_valve=True):
ScrollComp = ScrollClass()
ScrollComp.set_scroll_geo(83e-6, 2.2, 0.005,
0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
# print(ScrollComp.geo); quit()
ScrollComp.omega = 3000 / 60 * 2 * pi
ScrollComp.Tamb = 298.0
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 0.05
ScrollComp.HTC = 0.01
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Tin = Te + 11.1
temp = State.State(Ref, {'T': Te, 'Q': 1})
pe = temp.p
temp.update(dict(T=Tc, Q=1))
pc = temp.p
inletState = State.State(Ref, {'T': Tin, 'P': pe})
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.7)))
if not discharge_valve:
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
else:
E = 1.93e11 # Youngs Modulus, [Pa]
h_valve = 0.0006 # Valve thickness, [m]
d_discharge = ScrollComp.geo.ra_arc1 * 1.9 # Port diameter [m]
l_valve = 5 * d_discharge # Total length of valve, [m]
a_valve = l_valve / 1.5 # Distance from anchor to force, [m]
assert (a_valve < l_valve)
rho_valve = 8000 # Density of spring steel, [kg/m^3]
C_D = 1.17 # Drag coefficient [-]
d_valve = d_discharge * 1.5 # Valve Diameter [m]
x_stopper = 0.006 # Stopper location [m]
I = (d_valve * h_valve**3) / 12 # Moment of Inertia for valve,[m^4]
k_valve = (6 * E * I) / (a_valve**2 *
(3 * l_valve - a_valve)) # Valve stiffness
m_eff = (
1 / 3
) * rho_valve * l_valve * d_valve * h_valve # Effective mass of valve reeds
x_tr_discharge = 0.25 * (d_discharge**2 / d_valve)
# Construct the valve
ScrollComp.discharge_valve = ValveModel(
d_valve=d_valve,
d_port=d_discharge,
C_D=C_D,
m_eff=m_eff,
k_valve=k_valve,
rho_valve=rho_valve, # Not used directly
x_stopper=x_stopper,
x_tr=x_tr_discharge,
key_up=['ddd', 'dd'],
key_down='outlet.1')
# Inform the model about the valve
ScrollComp.add_valve(ScrollComp.discharge_valve)
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DischargeValve,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DischargeValve,
MdotFcn_kwargs=dict(X_d=0.7)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='d1', key2='dd', MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(
FlowPath(key1='d2', key2='dd', MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
t1 = timeit.default_timer()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
#hmin = 1e-3
eps_cycle=3e-3)
except BaseException as E:
print(E)
raise
print('time taken', timeit.default_timer() - t1)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
debug_plots(ScrollComp, family='Scroll Compressor')
return ScrollComp
def Expander():
Ref = 'Nitrogen'
Tinlet = 300
pinlet = 2000
poutlet = 500
SE = ScrollExpander()
SE.set_scroll_geo(100e-6,2.3,0.005,0.005)
SE.set_disc_geo('2Arc', r2 = 0)
SE.geo.delta_suction_offset = 0.0e-3
SE.geo.phi_ie_offset = 0.0
SE.omega = 60*2*pi
SE.Tamb = 298.0
#Temporarily set the bearing dimensions
SE.mech = core.struct()
SE.mech.D_upper_bearing = 0.02540
SE.mech.L_upper_bearing = 0.031
SE.mech.c_upper_bearing = 10e-6 #??
SE.mech.D_crank_bearing = 0.0254
SE.mech.L_crank_bearing = 0.023
SE.mech.c_crank_bearing = 10e-6 #??
SE.mech.D_lower_bearing = 0.01895
SE.mech.L_lower_bearing = 0.023
SE.mech.c_lower_bearing = 10e-6 #??
SE.mech.thrust_ID = 0.05
SE.mech.thrust_friction_coefficient = 0.02 #Tuned for point A1
SE.mech.L_ratio_bearings = 5
SE.mech.mu_oil = 0.008
SE.mech.scroll_plate_thickness = 0.007
SE.mech.scroll_density = 8100
SE.mech.scroll_added_mass = 0
SE.mech.scroll_plate_diameter = 0.092 #[m]
m, zcm = SE.calculate_scroll_mass()
SE.mech.orbiting_scroll_mass = m
SE.mech.scroll_zcm__thrust_surface = zcm
SE.mech.detailed_analysis = False
SE.h_shell = 10
SE.A_shell = 0.05
SE.HTC = 0.5
SE.generator = Motor()
SE.generator.set_eta(1.0)
SE.generator.suction_fraction = 0.0
inletState = CPState(Ref,{'T':Tinlet,'P':pinlet})
T2s = SE.guess_outlet_temp(inletState,poutlet)
outletState = CPState(Ref,{'T':T2s,'P':poutlet})
SE.geo.delta_radial = 0e-6
SE.geo.delta_flank = 0e-6
mdot_guess = inletState.rho*SE.Vdisp/SE.Vratio*SE.omega/(2*pi)
SE.add_tube(Tube(key1 = 'inlet.1',
key2 = 'inlet.2',
L = 0.03,
ID = 0.02,
mdot = mdot_guess,
State1 = inletState.copy(),
fixed = 1,
TubeFcn = SE.TubeCode))
SE.add_tube(Tube(key1 = 'outlet.1',
key2 = 'outlet.2',
L = 0.03,
ID = 0.02,
mdot = mdot_guess,
State2 = outletState.copy(),
fixed = 2,
TubeFcn = SE.TubeCode))
SE.auto_add_CVs(inletState, outletState)
# SE.auto_add_leakage(flankFunc = SE.FlankLeakage,
# radialFunc = SE.RadialLeakage)
FP = FlowPath(key1='outlet.1',
key2='da',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.02**2/4
SE.add_flow(FP)
SE.add_flow(FlowPath(key1 = 'da',
key2 = 'd1',
MdotFcn = SE.DA_D,
MdotFcn_kwargs = dict(X_d = 0.5)
)
)
SE.add_flow(FlowPath(key1 = 'da',
key2 = 'd2',
MdotFcn = SE.DA_D,
MdotFcn_kwargs = dict(X_d = 0.5)
)
)
FP = FlowPath(key1='inlet.2',
key2='ss',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.007**2/4
SE.add_flow(FP)
FP = FlowPath(key1='inlet.2',
key2='sss',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.007**2/4
SE.add_flow(FP)
SE.add_flow(FlowPath(key1='s1',
key2='ss',
MdotFcn=SE.S_to_SS))
SE.add_flow(FlowPath(key1='s2',
key2='ss',
MdotFcn=SE.S_to_SS))
SE.connect_callbacks(step_callback = SE.step_callback,
lumps_energy_balance_callback = SE.lump_energy_balance_callback,
endcycle_callback= SE.endcycle_callback)
# beta = np.linspace(0,2*pi,1000)
# V_e11 = np.array([SE.V_e1(b,1)[0] for b in beta])
# V_e21 = np.array([SE.V_e2(b,1)[0] for b in beta])
# V_e12 = np.array([SE.V_e1(b,2)[0] for b in beta])
# V_e22 = np.array([SE.V_e2(b,2)[0] for b in beta])
# V_d1 = np.array([SE.V_d1(b)[0] for b in beta])
# V_d2 = np.array([SE.V_d2(b)[0] for b in beta])
# V_s1 = np.array([SE.V_s1(b)[0] for b in beta])
#
# plt.plot(beta,V_e11)
# plt.plot(beta,V_e12)
# plt.plot(beta,V_d1)
# plt.plot(beta,V_s1)
# plt.gca().axvline(2*pi-SE.theta_d)
# plt.show()
#try:
SE.eps_cycle = 0.003
SE.precond_solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
UseNR = False, #If True, use Newton-Raphson ND solver to determine the initial state
OneCycle = False,
plot_every_cycle= False,
hmin = 1e-8
)
def Compressor(hs,
Vdisp,
Vratio,
Thickness,
OrbitingRadius,
phi_i0=0.0,
phi_os=0.3,
phi_is=pi,
phi_ie=21.76458,
phi_o0=-1.06276,
rb=0.004319,
f=None,
OneCycle=False,
Ref='R410A',
Qct=1,
do=0.123,
r_upper_bearing=0.02,
r_lower_bearing=0.02,
c_ub=1e-5):
#def Compressor(Te = 273, Tc = 300,Tamb = 25, f = None, OneCycle = False, Ref = 'R410A'):
Te = 280.37
Tc = 327.594
Tamb = 35
ScrollComp = Scroll()
#This runs if the module code is run directly
#ScrollComp.set_scroll_geo(1.53204e-5, 2.63, 0.00459, 0.004099,0.0,1.442,3.6)
print(Vdisp, Vratio, Thickness, OrbitingRadius, phi_os, phi_is)
ScrollComp.set_scroll_geo(Vdisp,
Vratio,
Thickness,
OrbitingRadius,
phi_i0=0.0,
phi_os=phi_os,
phi_is=phi_is,
phi_ie=phi_ie,
phi_o0=phi_o0,
hs=hs,
rb=rb)
#ScrollComp.set_scroll_geo(2.1723e-4, 2.63, 0.00459582, 0.00565,0.0,1.442,3.639) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 1.9e-5
ScrollComp.geo.delta_radial = 1.9e-5
#X_d=1
ScrollComp.geo.delta_suction_offset = 0.0e-3
#ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.Tamb = Tamb + 273.15 #[K]
ScrollComp.omega = 3600 / 60 * 2 * pi
f = 3600 / 60 # frequency for port size calculaton
#ScrollComp.Tamb = 298.0
#ScrollComp.journal_tune_factor=1.3
# if Qct<2.6:
# d_inlet=0.019
# elif Qct> 2.6 and Qct<7.6:
# d_inlet=0.022
# elif Qct> 7.6 and Qct<15:
# d_inlet=0.02857
# elif Qct> 15 and Qct<20:
# d_inlet=0.0412
# elif Qct> 20 and Qct<29:
# d_inlet=0.0539
# elif Qct> 29 and Qct<41:
# d_inlet=0.0457
# elif Qct> 41 and Qct<88:
# d_inlet=0.0666
# elif Qct> 88 and Qct<105:
# d_inlet=0.092
# else:
# d_inlet=0.092
#
#
if Qct < 4.5:
d_outlet = 0.0127
elif Qct > 4.5 and Qct < 13.5:
d_outlet = 0.020
elif Qct > 13.5 and Qct < 16.5:
d_outlet = 0.02857
elif Qct > 16.5 and Qct < 29:
d_outlet = 0.0349
elif Qct > 29 and Qct < 87:
d_outlet = 0.041
elif Qct > 87 and Qct < 105:
d_outlet = 0.0539
else:
d_outlet = 0.0539
# d_inlet=0.01
# d_outlet=0.003
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 2 * r_upper_bearing
ScrollComp.mech.L_upper_bearing = 2 * r_upper_bearing
ScrollComp.mech.c_upper_bearing = c_ub
ScrollComp.mech.D_crank_bearing = 2 * r_upper_bearing
ScrollComp.mech.L_crank_bearing = 2 * r_upper_bearing
ScrollComp.mech.c_crank_bearing = c_ub
ScrollComp.mech.D_lower_bearing = 2 * r_lower_bearing
ScrollComp.mech.L_lower_bearing = 2 * r_lower_bearing
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = do
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 5
ScrollComp.mech.mu_oil = 0.02
ScrollComp.mech.oldham_ring_radius = 0.06 #m
ScrollComp.mech.oldham_mass = 0.1 #kg
ScrollComp.mech.oldham_thickness = 0.008 #m
ScrollComp.mech.oldham_key_height = 0.006 #m
ScrollComp.mech.oldham_key_width = 0.006 #m
ScrollComp.mech.oldham_key_friction_coefficient = 0.01 #-
ScrollComp.mech.oldham_rotation_beta = 0 #rad
# ScrollComp.h_shell = 0.02
# ScrollComp.A_shell = 0.05
# ScrollComp.HTC = 3
ScrollComp.h_shell = 0.01
ScrollComp.A_shell = (2 * do + hs) * ((do)**2 * pi / 4)
#ScrollComp.HTC = 0.0
ScrollComp.HT_corr = 'Dittus-Boelter' #'Jang-Jeong'
# Temperature Lumps
ScrollComp.OEB_type = 'single-lump' #'single-lump'
ScrollComp.OEB_solver = 'MDNR'
ScrollComp.Rshell_oil = 190 #K/kW from Chen (2000) - PhD thesis
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
# Te = Te+ 273.15
# Tc = Tc+ 273.15
Tin = Te + 11.1
DT_sc = 8.3
temp = State.State(Ref, {'T': Te, 'Q': 1})
pe = temp.p
#pe=2950
temp.update(dict(T=Tc, Q=1))
pc = temp.p
PR = pc / pe
X_d = 0.7
#pc=4500
inletState = State.State(Ref, {'T': Tin, 'P': pe})
# defining the port area w.r.t displacement volume
V = 10
A = Vdisp * f / V
d_inlet = (4 * A / 3.14)**0.5
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=d_inlet,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=d_outlet,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
#
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * d_inlet**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=X_d)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=X_d)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=X_d)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=X_d)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='d1',
key2='dd',
MdotFcn=ScrollComp.D_to_DD,
MdotFcn_kwargs=dict(X_d=1)))
ScrollComp.add_flow(
FlowPath(key1='d2',
key2='dd',
MdotFcn=ScrollComp.D_to_DD,
MdotFcn_kwargs=dict(X_d=1)))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
from time import clock
t1 = clock()
ScrollComp.RK45_eps = 1e-4
ScrollComp.eps_cycle = 1e-3
try:
ScrollComp.precond_solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
x0=[320, 320], #Guesses [Td,Tlump[0],Tlump[1]]
hmin=1e-4,
eps_cycle=1e-3,
eps_energy_balance=0.1 #relaxed multi-lump convergence
)
except BaseException as E:
print(E)
raise
print('time taken', clock() - t1)
#debug_plots(ScrollComp)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
import CoolProp
Filepath1 = 'C:\Users\Mohsin\OneDrive - Oklahoma A and M System\Documents\Phd\Scroll comp scalling\cases_final\case 8/'
Filepath2 = 'E:\OneDrive - Oklahoma A and M System\Documents\Phd\Scroll comp scalling\cases_final\case 8/'
if os.path.isfile(Filepath1) == True:
Filepath = Filepath1
else:
Filepath = Filepath2
FileName = Filepath + 'Qc_t' + '_' + str(Qct) + '.h5'
h5.write_to_file(ScrollComp, FileName)
return ScrollComp
def Compressor(**kwargs):
recip = PURecip() #Instantiate the model
recip.Vdead = 0.5e-6 #Dead volume [m3]
recip.Vdisp = 8e-6 #Displacement/rev [m3]
recip.omega = 377 #Frequency, rad/sec (60Hz)
recip.d_discharge = 0.0059
#discharge port diameter [m]
recip.d_suction = recip.d_discharge
#suction diameter [m]
recip.A_discharge = pi * recip.d_discharge**2 / 4
recip.A_suction = pi * recip.d_suction**2 / 4
#These are parameters needed for the ambient heat transfer model
recip.h_shell = 0.010 #[kW/m2/K]
recip.A_shell = pi * 10 * 2 * (0.0254**2) #[m2]
recip.Tamb = 298 #[K]
#Parameters for the mechanical losses model (simplified)
recip.Wdot_parasitic = 0.01 #Parasitic losses [kW]
Ref = 'Air'
inletState = State.State(Ref, dict(T=298.15, P=101.325))
outletState = State.State(Ref, {'T': 400, 'P': inletState.p * 4})
mdot_guess = inletState.rho * recip.Vdisp * recip.omega / (2 * pi)
#First add the control volumes.
recip.add_CV(
ControlVolume(key='A',
initialState=outletState.copy(),
VdVFcn=recip.V_dV,
becomes='A'))
#Add the inlet tube
recip.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.03,
ID=0.01,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=recip.TubeCode))
#Add the outlet tube
recip.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.03,
ID=0.01,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=recip.TubeCode))
#Add the flow paths that link flow nodes together
recip.add_flow(FlowPath(key1='inlet.2', key2='A', MdotFcn=recip.Suction))
recip.add_flow(FlowPath(key1='outlet.1', key2='A',
MdotFcn=recip.Discharge))
recip.connect_callbacks(
endcycle_callback=recip.endcycle_callback, # Provided by PDSimCore
heat_transfer_callback=recip.heat_transfer_callback,
lumps_energy_balance_callback=recip.lump_energy_balance_callback)
t1 = timeit.default_timer()
recip.solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method=kwargs.get('solver_method', 'Euler'),
OneCycle=False,
UseNR=False)
print('time taken', timeit.default_timer() - t1, 's')
def Compressor(ScrollClass,
Te=273,
Tc=300,
f=None,
OneCycle=False,
Ref='R410A',
HDF5file='scroll_compressor.h5'):
ScrollComp = ScrollClass()
# This runs if the module code is run directly
ScrollComp.set_scroll_geo(83e-6, 3.3, 0.005,
0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = 3000 / 60 * 2 * pi
ScrollComp.Tamb = 298.0
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 0.05
ScrollComp.HTC = 0.0
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Te = -20 + 273.15
Tc = 20 + 273.15
Tin = Te + 11.1
temp = State.State(Ref, {'T': Te, 'Q': 1})
pe = temp.p
temp.update(dict(T=Tc, Q=1))
pc = temp.p
inletState = State.State(Ref, {'T': Tin, 'P': pe})
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
sim = ScrollComp
phi_list = [sim.geo.phi_fie - 2 * pi - 0.3]
involute_list = ['i']
D_list = [sim.geo.t] * len(phi_list)
offset_list = [d / 1.5 for d in D_list]
X_d_list = [1.0] * len(phi_list)
X_d_backflow_list = [0.0] * len(phi_list)
port_tube_list = [0]
Ltube_list = [0.001]
IDtube_list = [0.003]
# A superheated state
Tsat_inj = (Te + Tc) / 2 # [K]
temp = State.State(Ref, {'T': Tsat_inj, 'Q': 1})
temp = State.State(Ref, {'T': Tsat_inj + 10, 'P': temp.p})
Ttube_list = [temp.T]
RHOtube_list = [temp.rho]
sim.fixed_scroll_ports = []
for phi, involute, offset, D, parent, X_d, X_d_backflow in zip(
phi_list, involute_list, offset_list, D_list, port_tube_list,
X_d_list, X_d_backflow_list):
p = Port()
p.phi = phi
p.involute = involute
p.offset = offset
p.D = D
p.parent = parent
p.X_d = X_d
p.X_d_backflow = X_d_backflow
sim.fixed_scroll_ports.append(p)
more_keys = []
for i, (L, ID, T, D) in enumerate(
zip(Ltube_list, IDtube_list, Ttube_list, RHOtube_list)):
sim.add_tube(
Tube(key1='VITube' + str(i + 1) + '.1',
key2='VITube' + str(i + 1) + '.2',
L=L,
ID=ID,
mdot=mdot_guess * 0.1,
State1=State.State(inletState.Fluid, dict(T=T, D=D)),
fixed=1,
TubeFcn=sim.TubeCode))
more_keys.append('VITube' + str(i + 1) + '.1')
sim.additional_inlet_keys = more_keys
# # This block can be uncommented to plot the scroll wrap with the VI ports
# plotScrollSet(0, geo = sim.geo)
# for port in sim.fixed_scroll_ports:
# x, y = scroll_geo.coords_inv(port.phi, sim.geo, 0, port.involute)
# nx, ny = scroll_geo.coords_norm(port.phi, sim.geo, 0, port.involute)
# x0 = x - nx*port.offset
# y0 = y - ny*port.offset
# t = np.linspace(0, 2*pi)
# plt.plot(port.D/2.0*np.cos(t) + x0, port.D/2.0*np.sin(t) + y0,'b')
# plt.savefig('wrap_with_VI_ports.png',transparent=True)
# plt.close()
# quit()
# Calculate the areas between each VI port and every control volume
sim.calculate_port_areas()
for port in sim.fixed_scroll_ports:
for partner in port.area_dict:
# Create a spline interpolator object for the area between port and the partner chamber
A_interpolator = interp.splrep(port.theta,
port.area_dict[partner],
k=2,
s=0)
# Add the flow between the EVI control volume and the chamber through the port
sim.add_flow(
FlowPath(key1='VITube' + str(port.parent + 1) + '.2',
key2=partner,
MdotFcn=sim.INTERPOLATING_NOZZLE_FLOW,
MdotFcn_kwargs=dict(X_d=port.X_d,
X_d_backflow=port.X_d_backflow,
upstream_key='VITube' +
str(port.parent + 1) + '.2',
A_interpolator=A_interpolator)))
ScrollComp.add_flow(
FlowPath(key1='d1', key2='dd', MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(
FlowPath(key1='d2', key2='dd', MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
t1 = timeit.default_timer()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
#hmin = 1e-3
eps_cycle=3e-3)
except BaseException as E:
print(E)
raise
print('time taken', timeit.default_timer() - t1)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
return ScrollComp
def Compressor(ScrollClass, Te=273, Tc=300, f=None, OneCycle=False, Ref='R410A', HDF5file='scroll_compressor.h5'):
ScrollComp = ScrollClass()
# This runs if the module code is run directly
ScrollComp.set_scroll_geo(83e-6, 3.3, 0.005, 0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = 3000/60*2*pi
ScrollComp.Tamb = 298.0
#Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 #From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02
ScrollComp.A_shell = 0.05
ScrollComp.HTC = 0.0
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Te = -20 + 273.15
Tc = 20 + 273.15
Tin = Te + 11.1
DT_sc = 7
temp = State.State(Ref,{'T':Te,'Q':1})
pe = temp.p
temp.update(dict(T=Tc, Q=1))
pc = temp.p
inletState = State.State(Ref,{'T':Tin,'P':pe})
T2s = ScrollComp.guess_outlet_temp(inletState,pc)
outletState = State.State(Ref,{'T':T2s,'P':pc})
mdot_guess = inletState.rho*ScrollComp.Vdisp*ScrollComp.omega/(2*pi)
ScrollComp.add_tube(Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi*0.01**2/4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.7)
)
)
ScrollComp.add_flow(FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.7)
)
)
ScrollComp.add_flow(FlowPath(key1='outlet.1',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)
)
)
ScrollComp.add_flow(FlowPath(key1='outlet.1',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)
)
)
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
sim = ScrollComp
phi_list = [sim.geo.phi_fie-2*pi-0.3]
involute_list = ['i']
D_list = [sim.geo.t]*len(phi_list)
offset_list = [d/1.5 for d in D_list]
X_d_list = [1.0]*len(phi_list)
X_d_backflow_list = [0.0]*len(phi_list)
port_tube_list = [0]
Ltube_list = [0.001]
IDtube_list = [0.003]
# A superheated state
Tsat_inj = (Te+Tc)/2 # [K]
temp = State.State(Ref,{'T':Tsat_inj,'Q':1})
temp = State.State(Ref,{'T':Tsat_inj+10,'P':temp.p})
Ttube_list = [temp.T]
RHOtube_list = [temp.rho]
sim.fixed_scroll_ports = []
for phi,involute,offset,D,parent,X_d,X_d_backflow in zip(phi_list,involute_list,offset_list,D_list,port_tube_list,X_d_list,X_d_backflow_list):
p = Port()
p.phi = phi
p.involute = involute
p.offset = offset
p.D = D
p.parent = parent
p.X_d = X_d
p.X_d_backflow = X_d_backflow
sim.fixed_scroll_ports.append(p)
more_keys = []
for i, (L, ID, T, D) in enumerate(zip(Ltube_list, IDtube_list, Ttube_list, RHOtube_list)):
sim.add_tube(Tube(key1='VITube'+str(i+1)+'.1',
key2='VITube'+str(i+1)+'.2',
L=L,
ID=ID,
mdot=mdot_guess*0.1,
State1=State.State(inletState.Fluid,
dict(T=T, D=D)
),
fixed=1,
TubeFcn=sim.TubeCode
)
)
more_keys.append('VITube'+str(i+1)+'.1')
sim.additional_inlet_keys = more_keys
# # This block can be uncommented to plot the scroll wrap with the VI ports
# plotScrollSet(0, geo = sim.geo)
# for port in sim.fixed_scroll_ports:
# x, y = scroll_geo.coords_inv(port.phi, sim.geo, 0, port.involute)
# nx, ny = scroll_geo.coords_norm(port.phi, sim.geo, 0, port.involute)
# x0 = x - nx*port.offset
# y0 = y - ny*port.offset
# t = np.linspace(0, 2*pi)
# plt.plot(port.D/2.0*np.cos(t) + x0, port.D/2.0*np.sin(t) + y0,'b')
# plt.savefig('wrap_with_VI_ports.png',transparent=True)
# plt.close()
# quit()
# Calculate the areas between each VI port and every control volume
sim.calculate_port_areas()
for port in sim.fixed_scroll_ports:
for partner in port.area_dict:
# Create a spline interpolator object for the area between port and the partner chamber
A_interpolator = interp.splrep(port.theta, port.area_dict[partner], k=2, s=0)
# Add the flow between the EVI control volume and the chamber through the port
sim.add_flow(FlowPath(key1='VITube'+str(port.parent+1)+'.2',
key2=partner,
MdotFcn=sim.INTERPOLATING_NOZZLE_FLOW,
MdotFcn_kwargs=dict(X_d=port.X_d,
X_d_backflow=port.X_d_backflow,
upstream_key='VITube'+str(port.parent+1)+'.2',
A_interpolator=A_interpolator
)
)
)
ScrollComp.add_flow(FlowPath(key1='d1',
key2='dd',
MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(FlowPath(key1='d2',
key2='dd',
MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback
)
t1 = timeit.default_timer()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.precond_solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
#hmin = 1e-3
eps_cycle=3e-3
)
except BaseException as E:
print(E)
raise
print('time taken', timeit.default_timer()-t1)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
return ScrollComp
def Expander(**kwargs):
expander = PistonExpander() #Instantiate the class
Ref = 'Nitrogen'
inletState = State.State(Ref,dict(T = 298.15, P = 501.325))
outletState = State.State(Ref,dict(T = 200, P = inletState.p/10))
mdot_guess = inletState.rho*expander.Vdisp*expander.omega/(2*pi)
#First add the control volume.
expander.add_CV(ControlVolume(key='A',
initialState=inletState.copy(),
VdVFcn=expander.V_dV,)
)
#These are parameters needed for the ambient heat transfer model
expander.h_shell = 0.010 #[kW/m2/K]
expander.A_shell = pi*10*2*(0.0254**2) #[m2]
expander.Tamb = 298 #[K]
#Parameters for the mechanical losses model (simplified)
expander.Wdot_parasitic = 0.01 #Parasitic losses [kW]
#Add the inlet tube
expander.add_tube(Tube(key1 = 'inlet.1',
key2 = 'inlet.2',
L = 0.03,
ID = 0.01,
mdot = mdot_guess,
State1 = inletState.copy(),
fixed = 1,
TubeFcn = expander.TubeCode)
)
#Add the outlet tube
expander.add_tube(Tube(key1 = 'outlet.1',
key2 = 'outlet.2',
L = 0.03,
ID = 0.01,
mdot = mdot_guess,
State2 = outletState.copy(),
fixed = 2,
TubeFcn = expander.TubeCode)
)
#Add the flow paths that link flow nodes together
expander.add_flow(FlowPath(key1='inlet.2',key2='A',MdotFcn=expander.Suction))
expander.add_flow(FlowPath(key1='outlet.1',key2='A',MdotFcn=expander.Discharge))
t1 = timeit.default_timer()
expander.EulerN = 4000
expander.RK45_eps = 1e-10
expander.connect_callbacks(step_callback = expander.step_callback,
endcycle_callback=expander.endcycle_callback, # Provided by PDSimCore
heat_transfer_callback=expander.heat_transfer_callback,
lumps_energy_balance_callback = expander.lump_energy_balance_callback)
expander.solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method = kwargs.get('solver_method', 'Euler'),
OneCycle = False,
UseNR = True,
plot_every_cycle = False,
eps_cycle=3e-3,
eps_energy_balance=3e-3
)
print('time taken',timeit.default_timer()-t1,'s')
def Compressor():
recip=Recip()
recip.piston_stroke = 0.02 #Piston stroke, m
recip.piston_diameter = 0.02 #Piston diameter, m
recip.piston_length = 0.02 #Piston Length, m
recip.omega = 377 #Frequency, rad/sec (60Hz)
recip.crank_length = 0.01 #length of crank, m
recip.connecting_rod_length = 0.04 #length of connecting rod, m
recip.x_TDC = 0.003 #Distance to the TDC position of the piston from the valve plate
recip.d_discharge=0.0059; #discharge port diameter in meters
recip.d_suction=recip.d_discharge; #suction diameter in meters
#These are parameters needed for the ambient heat transfer model
recip.h_shell = 0.010 #[kW/m2/K]
recip.A_shell = pi*10*2*(0.0254**2) #[m2]
recip.Tamb = 298 #[K]
recip.mu_oil = 0.0086
recip.delta_gap = 10e-6
recip.eta_motor = 0.95
recip.shell_volume = 100e-6
#Calculate Vdisp
recip.pre_solve()
recip.Vdisp = recip.Vdisp()
Ref='R410A'
inletState=State.State(Ref,dict(T=289.15, D=33.1))
p_outlet = inletState.p*2.5
T2s = recip.guess_outlet_temp(inletState,p_outlet)
outletState=State.State(Ref,{'T':T2s,'P':p_outlet})
mdot_guess = inletState.rho*recip.Vdisp*recip.omega/(2*pi)
#First add the control volumes.
recip.add_CV( ControlVolume(key='A',
initialState=outletState.copy(),
VdVFcn=recip.V_dV,
becomes='A') )
recip.add_CV( ControlVolume(key='shell',
initialState=inletState.copy(),
VdVFcn=recip.V_shell,
becomes='shell') )
recip.add_tube( Tube(key1='inlet.1',key2='inlet.2',L=0.03,ID=0.02,
mdot=mdot_guess, State1=inletState.copy(),
fixed=1,TubeFcn=recip.TubeCode) )
recip.add_tube( Tube(key1='outlet.1',key2='outlet.2',L=0.03,ID=0.02,
mdot=mdot_guess, State2=outletState.copy(),
fixed=2,TubeFcn=recip.TubeCode) )
recip.add_flow(FlowPath(key1='shell', key2='inlet.2', MdotFcn=recip.Inlet))
recip.add_flow(FlowPath(key1='inlet.2', key2='A', MdotFcn=recip.Suction))
recip.add_flow(FlowPath(key1='outlet.1', key2='A', MdotFcn=recip.Discharge))
recip.add_flow(FlowPath(key1='shell', key2='A', MdotFcn=recip.PistonLeakage))
E = 1.93e11 #Youngs Modulus, [Pa]
h_valve = 0.0001532 #Valve thickness, [m]
l_valve = 0.018 #Total length of valve, [m]
a_valve = 0.0140 #Distance from anchor to force, [m]
rho_valve = 8000 #Density of spring steel, [kg/m^3]
C_D = 1.17 #Drag coefficient [-]
d_valve = 0.007 #Valve Diameter [m]
x_stopper = 0.0018 #Stopper location [m]
I=(d_valve*h_valve**3)/12 #Moment of Intertia for valve,[m^4]
k_valve=(6*E*I)/(a_valve**2*(3*l_valve-a_valve)) #Valve stiffness
m_eff=(1/3)*rho_valve*l_valve*d_valve*h_valve #Effective mass of valve reeds
x_tr_suction = 0.25*(recip.d_suction**2/d_valve)
x_tr_discharge = 0.25*(recip.d_discharge**2/d_valve)
#The suction valve parameters
recip.suction_valve=ValveModel(
d_valve=d_valve,
d_port=recip.d_suction,
C_D=C_D,
rho_valve=rho_valve,
x_stopper=x_stopper,
m_eff = m_eff,
k_valve = k_valve,
x_tr = x_tr_suction,
key_up='inlet.2',
key_down='A'
)
recip.add_valve(recip.suction_valve)
#The discharge valve parameters
recip.discharge_valve=ValveModel(
d_valve=d_valve,
d_port=recip.d_discharge,
C_D=C_D,
rho_valve=rho_valve,
x_stopper=x_stopper,
m_eff = m_eff,
k_valve = k_valve,
x_tr = x_tr_discharge,
key_up='A',
key_down='outlet.1'
)
recip.add_valve(recip.discharge_valve)
recip.connect_callbacks(endcycle_callback=recip.endcycle_callback,
heat_transfer_callback=recip.heat_transfer_callback,
lumps_energy_balance_callback = recip.lump_energy_balance_callback
)
t1 = timeit.default_timer()
recip.precond_solve(key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method = 'RK45',
OneCycle = False,
UseNR = True,
)
print('time taken', timeit.default_timer()-t1)
debug_plots(recip)
del recip.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(recip, 'recipsample.h5')
def Compressor(**kwargs):
recip = PURecip() #Instantiate the model
recip.Vdead = 8e-8 #Dead volume [m3]
recip.Vdisp = 8e-6 #Displacement/rev [m3]
recip.omega = 377 #Frequency, rad/sec (60Hz)
recip.d_discharge = 0.0059
#Discharge port diameter [m]
recip.d_suction = recip.d_discharge
#Suction diameter [m]
recip.A_discharge = pi * recip.d_discharge**2 / 4
recip.A_suction = pi * recip.d_suction**2 / 4
recip.RK45_eps = 1e-8
recip.h_shell = 0.010 #[kW/m2/K]
recip.A_shell = pi * 10 * 2 * (0.0254**2) #[m2]
recip.Tamb = 293 #[K]
recip.plot_names = 'recip_plot_buttons'
recip.Wdot_parasitic = 0 #Parasitic losses [kW]
Ref = 'R134a'
inletState = State.State(Ref, dict(T=293.15, D=23.75))
outletState = State.State(Ref, {'T': 400, 'P': inletState.p * 2.5})
mdot_guess = inletState.rho * recip.Vdisp * recip.omega / (2 * pi)
#First add the control volumes.
recip.add_CV(
ControlVolume(key='A',
initialState=outletState.copy(),
VdVFcn=recip.V_dV,
becomes='A'))
#Add the inlet tube
recip.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.03,
ID=0.01,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=recip.TubeCode))
#Add the outlet tube
recip.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.03,
ID=0.01,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=recip.TubeCode))
#Add the flow paths that link flow nodes together
recip.add_flow(FlowPath(key1='inlet.2', key2='A', MdotFcn=recip.Suction))
recip.add_flow(FlowPath(key1='outlet.1', key2='A',
MdotFcn=recip.Discharge))
E = 1.93e11 #Youngs Modulus, [Pa]
h_valve = 0.0001532 #Valve thickness, [m]
l_valve = 0.018 #Total length of valve, [m]
a_valve = 0.0140 #Distance from anchor to force, [m]
rho_valve = 8000 #Density of spring steel, [kg/m^3]
C_D = 1.17 #Drag coefficient [-]
d_valve = 0.0059 #Valve Diameter [m]
x_stopper = 0.0018 #Stopper location [m]
I = (d_valve * h_valve**3) / 12 #Moment of Intertia for valve,[m^4]
k_valve = (6 * E * I) / (a_valve**2 *
(3 * l_valve - a_valve)) #Valve stiffness
m_eff = (
1 / 3
) * rho_valve * l_valve * d_valve * h_valve #Effective mass of valve reeds
x_tr_suction = 0.25 * (recip.d_suction**2 / d_valve)
x_tr_discharge = 0.25 * (recip.d_discharge**2 / d_valve)
#The suction valve parameters
recip.suction_valve = ValveModel(d_valve=d_valve,
d_port=recip.d_suction,
C_D=C_D,
rho_valve=rho_valve,
x_stopper=x_stopper,
m_eff=m_eff,
k_valve=k_valve,
x_tr=x_tr_suction,
key_up='inlet.2',
key_down='A')
recip.add_valve(recip.suction_valve)
#The discharge valve parameters
recip.discharge_valve = ValveModel(d_valve=d_valve,
d_port=recip.d_discharge,
C_D=C_D,
rho_valve=rho_valve,
x_stopper=x_stopper,
m_eff=m_eff,
k_valve=k_valve,
x_tr=x_tr_discharge,
key_up='A',
key_down='outlet.1')
recip.add_valve(recip.discharge_valve)
recip.connect_callbacks(
endcycle_callback=recip.endcycle_callback, # Provided by PDSimCore
heat_transfer_callback=recip.heat_transfer_callback,
lumps_energy_balance_callback=recip.lump_energy_balance_callback)
t1 = clock()
recip.precond_solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=False,
UseNR=True,
)
print('time taken', clock() - t1)
debug_plots(recip)
del recip.FlowStorage
# from PDSim.plot.plots import debug_plots
# debug_plots(recip)
return recip
def Compressor(ScrollClass,
*,
Te=273,
Tc=300,
OneCycle=False,
Ref='R410A',
HDF5file='scroll_compressor_disc.h5',
Nports=4):
ScrollComp = ScrollClass()
#This runs if the module code is run directly
ScrollComp.set_scroll_geo(83e-6, 3.3, 0.005,
0.006) #Set the scroll wrap geometry
ScrollComp.set_disc_geo('2Arc', r2=0)
ScrollComp.geo.delta_flank = 10e-6
ScrollComp.geo.delta_radial = 10e-6
ScrollComp.geo.delta_suction_offset = 0.0e-3
ScrollComp.geo.phi_ie_offset = 0.0
ScrollComp.omega = 3000 / 60 * 2 * pi
ScrollComp.Tamb = 298.0
# Temporarily set the bearing dimensions
ScrollComp.mech = struct()
ScrollComp.mech.D_upper_bearing = 0.04
ScrollComp.mech.L_upper_bearing = 0.04
ScrollComp.mech.c_upper_bearing = 20e-6
ScrollComp.mech.D_crank_bearing = 0.04
ScrollComp.mech.L_crank_bearing = 0.04
ScrollComp.mech.c_crank_bearing = 20e-6
ScrollComp.mech.D_lower_bearing = 0.025
ScrollComp.mech.L_lower_bearing = 0.025
ScrollComp.mech.c_lower_bearing = 20e-6
ScrollComp.mech.thrust_ID = 0.05
ScrollComp.mech.thrust_friction_coefficient = 0.028 # From Chen thesis
ScrollComp.mech.orbiting_scroll_mass = 2.5
ScrollComp.mech.L_ratio_bearings = 3
ScrollComp.mech.mu_oil = 0.008
ScrollComp.h_shell = 0.02 # kW/m^2/K
ScrollComp.A_shell = 0.05 # m^2
ScrollComp.HTC = 0.050 # Heat transfer coefficient between scroll and gas in chambers, in kW/m^2/K
ScrollComp.motor = Motor()
ScrollComp.motor.set_eta(0.9)
ScrollComp.motor.suction_fraction = 1.0
Tin = Te + 11.1
temp = State.State(Ref, {'T': Te, 'Q': 1})
pe = temp.p
temp.update(dict(T=Tc, Q=1))
pc = temp.p
inletState = State.State(Ref, {'T': Tin, 'P': pe})
T2s = ScrollComp.guess_outlet_temp(inletState, pc)
outletState = State.State(Ref, {'T': T2s, 'P': pc})
mdot_guess = inletState.rho * ScrollComp.Vdisp * ScrollComp.omega / (2 *
pi)
ScrollComp.add_tube(
Tube(key1='inlet.1',
key2='inlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State1=inletState.copy(),
fixed=1,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.add_tube(
Tube(key1='outlet.1',
key2='outlet.2',
L=0.3,
ID=0.02,
mdot=mdot_guess,
State2=outletState.copy(),
fixed=2,
TubeFcn=ScrollComp.TubeCode))
ScrollComp.auto_add_CVs(inletState, outletState)
ScrollComp.add_CV(
ControlVolume(
key='discplenum',
initialState=outletState.copy(),
VdVFcn=lambda theta:
(1e-4, 0
) # Constant volume, outputs of function are volume, dV/dtheta
))
ScrollComp.auto_add_leakage(flankFunc=ScrollComp.FlankLeakage,
radialFunc=ScrollComp.RadialLeakage)
FP = FlowPath(
key1='inlet.2',
key2='sa',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
FP = FlowPath(
key1='outlet.1',
key2='discplenum',
MdotFcn=IsentropicNozzleWrapper(),
)
FP.A = pi * 0.01**2 / 4
ScrollComp.add_flow(FP)
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s1',
MdotFcn=ScrollComp.SA_S1,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='sa',
key2='s2',
MdotFcn=ScrollComp.SA_S2,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='discplenum',
key2='dd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
ScrollComp.add_flow(
FlowPath(key1='discplenum',
key2='ddd',
MdotFcn=ScrollComp.DISC_DD,
MdotFcn_kwargs=dict(X_d=0.7)))
# ScrollComp.add_flow(FlowPath(key1 = 'outlet.1',
# key2 = 'd1',
# MdotFcn = ScrollComp.DISC_D1,
# MdotFcn_kwargs = dict(X_d = 0.7)
# )
# )
#
# FP = FlowPath(key1='outlet.1',
# key2='dd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
#
# FP = FlowPath(key1='outlet.1',
# key2='ddd',
# MdotFcn=IsentropicNozzleWrapper(),
# )
# FP.A = pi*0.006**2/4
# ScrollComp.add_flow(FP)
# ***************
# Discharge Ports
# ***************
sim = ScrollComp
for involute in ['i', 'o']:
# Define the parameters for the ports
if involute == 'i':
phi_list = np.linspace(sim.geo.phi_fie - 2 * np.pi - 0.1,
sim.geo.phi_fis + np.pi / 2, Nports)
elif involute == 'o':
phi_list = np.linspace(sim.geo.phi_foe - 2 * np.pi - 0.1,
sim.geo.phi_fos + np.pi / 2, Nports)
else:
raise ValueError(involute)
involute_list = [involute] * len(phi_list)
D_list = [sim.geo.t] * len(
phi_list) # All ports are thickness of scroll in diameter
offset_list = [
d * 0.5 for d in D_list
] # Distance of translation of center of port normal to wall
X_d_list = [1.0] * len(phi_list) # Flow coefficient for forward flow
X_d_backflow_list = [0.0] * len(phi_list) # Backflow coefficient
# Add the ports to the model
sim.fixed_scroll_ports = []
for i in range(Nports):
p = Port()
p.phi = phi_list[i]
p.involute = involute_list[i]
p.offset = offset_list[i]
p.D = D_list[i]
p.parent = 'discplenum'
p.X_d = X_d_list[i]
p.X_d_backflow = X_d_backflow_list[i]
sim.fixed_scroll_ports.append(p)
# Calculate the areas between each port and every control volume
#
# Can turn on plotting (plot=True) to see the flow areas.
sim.calculate_port_areas(plot=False)
# Iterate over the ports, connecting them to discharge plenum
for port in sim.fixed_scroll_ports:
# Iterate over the chambers each port is connected to at least
# at some point
for partner in port.area_dict:
# Create a spline interpolator object for the area between port and the partner chamber
A_interpolator = interp.splrep(port.theta,
port.area_dict[partner],
k=2,
s=0)
# Add the flow path connecting the working chamber and discharge plenum
#
# It is a simplified nozzle model, allowing for flow out the nozzle,
# according to the flow area at the given crank angle, but no backflow.
# This represents an idealized case.
sim.add_flow(
FlowPath(key1='discplenum',
key2=partner,
MdotFcn=sim.INTERPOLATING_NOZZLE_FLOW,
MdotFcn_kwargs=dict(
X_d=port.X_d,
X_d_backflow=port.X_d_backflow,
upstream_key=partner,
A_interpolator=A_interpolator)))
ScrollComp.add_flow(
FlowPath(key1='d1', key2='dd', MdotFcn=ScrollComp.D_to_DD))
ScrollComp.add_flow(
FlowPath(key1='d2', key2='dd', MdotFcn=ScrollComp.D_to_DD))
#Connect the callbacks for the step, endcycle, heat transfer and lump energy balance
ScrollComp.connect_callbacks(
step_callback=ScrollComp.step_callback,
endcycle_callback=ScrollComp.endcycle_callback,
heat_transfer_callback=ScrollComp.heat_transfer_callback,
lumps_energy_balance_callback=ScrollComp.lump_energy_balance_callback)
t1 = timeit.default_timer()
ScrollComp.RK45_eps = 1e-6
ScrollComp.eps_cycle = 3e-3
try:
ScrollComp.solve(
key_inlet='inlet.1',
key_outlet='outlet.2',
solver_method='RK45',
OneCycle=OneCycle,
plot_every_cycle=False,
#hmin = 1e-3
eps_cycle=3e-3)
except BaseException as E:
print(E)
raise
print('time taken', timeit.default_timer() - t1)
del ScrollComp.FlowStorage
from PDSim.misc.hdf5 import HDF5Writer
h5 = HDF5Writer()
h5.write_to_file(ScrollComp, HDF5file)
if '--plot' in sys.argv:
debug_plots(ScrollComp, family='Scroll Compressor')
return ScrollComp