Esempio n. 1
0
class DuctTestCase(unittest.TestCase):

  

    def tearDown(self): 
        comp = None

    def test_duct(self): 
    	   
    	# initial test case based on JT9D 34000 .8 operation    
    	self.duct = set_as_top(Duct())
        self.fs = FlowStation()
        self.fs.W = 100.
        self.fs.setTotalTP( 1000, 100 );
        
        duct = self.duct         
        duct.Fl_I = self.fs 
        duct.dPqP = .1
        duct.run()

        TOL = .001
        
        # check the flow and mechanical values
        assert_rel_error( self, duct.Fl_O.W, 100., TOL )
        assert_rel_error( self, duct.Fl_O.Pt, 90., TOL )
        assert_rel_error( self, duct.Fl_O.Tt, 1000., TOL )
Esempio n. 2
0
    def setUp(self): 
        self.comp = set_as_top(Compressor())

        self.fs = FlowStation()
        self.fs.W = 1.08
        self.fs.setTotalTP(630.74523, 0.0271945)
        self.fs.Mach = .6
Esempio n. 3
0
class InletTestCase(unittest.TestCase):

  

    def tearDown(self): 
        comp = None

    def test_duct(self): 
    	   
    	# initial test case based on JT9D 34000 .8 operation    
    	self.inlet = set_as_top(Inlet())
        self.fs = FlowStation()
        self.fs.W = 100.
        self.fs.setTotalTP( 1000, 100 );
        
        inlet = self.inlet       
        inlet.Fl_I = self.fs 
        inlet.ram_recovery = .9
        inlet.run()

        TOL = .001
        
        # check the flow and mechanical values
        assert_rel_error( self, inlet.Fl_O.W, 100., TOL )
        assert_rel_error( self, inlet.Fl_O.Pt, 90., TOL )
        assert_rel_error( self, inlet.Fl_O.Tt, 1000., TOL )
Esempio n. 4
0
        def configure(self):

            hx = self.add('hx', HeatExchanger())
            driver = self.add('driver',BroydenSolver())
            driver.add_parameter('hx.T_hot_out',low=0.,high=1000.)
            driver.add_parameter('hx.T_cold_out',low=0.,high=1000.)
            driver.add_constraint('hx.residual_qmax=0')
            driver.add_constraint('hx.residual_e_balance=0')

            #hx.Wh = 0.49
            #hx.Cp_hot = 1.006
            #hx.T_hot_in = 791
            fs = FlowStation()
            fs.setTotalTP(1423.8, 0.302712118187) #R, psi
            fs.W = .49
            hx.Fl_I = fs
            hx.W_cold = 0.45
            hx.Cp_cold = 4.186
            hx.T_cold_in = 288.15
            effectiveness = 0.9765

            #initial guess
            avg = ( hx.Fl_I.Tt*.555555556 + hx.T_cold_in )/2.
            hx.T_cold_out = avg
            hx.T_hot_out = avg  

            driver.workflow.add(['hx'])
Esempio n. 5
0
class CompressorRlineTestCase(unittest.TestCase):

  

    def tearDown(self): 
        comp = None

    def test_compressor(self): 
    	   
    	# initial test case based on JT9D 34000 .8 operation    
    	self.comp = set_as_top(CompressorRline())
        self.fs = FlowStation()
        self.fs.W = 111.61
        self.fs.setDryAir()
        self.fs.setTotalTP( 707.81, 22.395 )
        comp = self.comp
        comp.PR = 5.609
        comp.eff = .8651
        comp.Wfrac1=.055
        comp.hfrac1=1
        comp.Pfrac1=1
        comp.Wfrac2=.035
        comp.hfrac2=1
        comp.Pfrac2=1
        comp.Nmech=8000

        comp.Fl_I = self.fs
        comp.run()

        TOL = .001
        
        # check the flow and mechanical values
        assert_rel_error( self,comp.Fl_O.W, 101.57, TOL )
        assert_rel_error( self,comp.Fl_O.Pt, 125.611, TOL )
        assert_rel_error( self,comp.Fl_O.Tt, 1209.83, TOL )
 
        assert_rel_error( self,comp.Fl_bld1.W, 6.1386, TOL )
        assert_rel_error( self,comp.Fl_bld1.Pt, 125.611, TOL )
        assert_rel_error( self,comp.Fl_bld1.Tt, 1209.83, TOL )
        
        assert_rel_error( self,comp.Fl_bld2.W, 3.9064, TOL )
        assert_rel_error( self,comp.Fl_bld2.Pt, 125.611, TOL )
        assert_rel_error( self,comp.Fl_bld2.Tt, 1209.83, TOL )
        
        
        assert_rel_error( self,comp.pwr, -19639.2, TOL )
        assert_rel_error( self,comp.trq, -1350.435, TOL )
        
        assert_rel_error( self,comp.eff, .8651, TOL )

        # adjust the bleed fraction to ensure the temperature is about half
        # way between input and output
        comp.hfrac1 = .5
        comp.run()
 
        assert_rel_error( self,comp.Fl_bld1.Tt, 961.654, TOL )
Esempio n. 6
0
class CompressorTestCase(unittest.TestCase):

    def setUp(self): 
        self.comp = set_as_top(Compressor())

        self.fs = FlowStation()
        self.fs.W = 1.08
        self.fs.setTotalTP(630.74523, 0.0271945)
        self.fs.Mach = .6

    def tearDown(self): 
        comp = None

    def test_compressor(self): 
        comp = self.comp

        comp.PR_des = 12.47
        comp.MNexit_des = .4
        comp.eff_des = .80
        comp.Fl_I = self.fs
        comp.design = True

        comp.run()

        TOL = .001
        assert_rel_error(self,comp.Fl_O.W, 1.08,TOL)
        assert_rel_error(self,comp.Fl_O.Pt, .33899, TOL)
        assert_rel_error(self,comp.Fl_O.Tt, 1424.01, TOL)
        assert_rel_error(self,comp.Fl_O.rhos, .000594, TOL)
        assert_rel_error(self,comp.Fl_O.Mach, .4 ,TOL)
        assert_rel_error(self,comp.Fl_O.area, 364.7, TOL)
        assert_rel_error(self,comp.pwr, 303.2, TOL)
        assert_rel_error(self,comp.eff_poly, .8545, TOL)

        #run off design
        comp.run()

        #values should remain unchanged in off-design at design condition
        assert_rel_error(self,comp.Fl_O.W, 1.08,TOL)
        assert_rel_error(self,comp.Fl_O.Pt, .33899, TOL)
        assert_rel_error(self,comp.Fl_O.Tt, 1424.01, TOL)
        assert_rel_error(self,comp.Fl_O.rhos, .000594, TOL)
        assert_rel_error(self,comp.Fl_O.Mach, .4 ,TOL)
        assert_rel_error(self,comp.Fl_O.area, 364.7, TOL)
        assert_rel_error(self,comp.pwr, 303.2, TOL)
        assert_rel_error(self,comp.eff_poly, .8545, TOL)

        #try changing something
        self.fs.W *= 1.1
        comp.Fl_I = self.fs
        comp.run()
        assert_rel_error(self, comp.PR, 13.52995, TOL)
Esempio n. 7
0
class NozzleConvergentCase(unittest.TestCase):

  

    def tearDown(self): 
        comp = None

    def test_nozzleconvergent(self): 
    	   
    	# initial test case based on JT9D 34000 .8 operation    
    	# choked
        self.nozz = set_as_top(NozzleConvergent())
        self.fs = FlowStation()
        self.fs.W = 562.61
        self.fs.setTotalTP( 527.55, 8.978 );
        
        nozz = self.nozz       
        nozz.Fl_I = self.fs 
        nozz.Cfg = .9962
        nozz.PsExh = 3.626
        nozz.run()

        TOL = .001
        
        # check the flow and mechanical values
        assert_rel_error( self, nozz.Fg, 21033.3, TOL )
   
   	# initial test case based on JT9D 0. 0. operation    
    	# unchoked
        self.nozz = set_as_top(NozzleConvergent())
        self.fs = FlowStation()
        self.fs.W = 1293.91
        self.fs.setTotalTP( 632.71, 23.195 );
        
        nozz = self.nozz       
        nozz.Fl_I = self.fs 
        nozz.Cfg = .9975
        nozz.PsExh = 14.696
        nozz.run()

        TOL = .001
        
        # check the flow and mechanical values
        assert_rel_error( self, nozz.Fg, 38786, TOL )
Esempio n. 8
0
    def configure(self):

        #Add Compressor Cycle Components
        tube = self.add('tube', FlowStartStatic())
        #tube.W = 1.521
        tube.Ps = 0.01436
        tube.Ts = 525.6

        inlet = self.add('inlet', Inlet())
        inlet.ram_recovery = 1.0
        #inlet.MNexit_des = .6

        comp1 = self.add('comp1', Compressor())
        comp1.PR_des = 12.47
        comp1.MNexit_des = .4
        comp1.eff_des = .80

        duct1 = self.add('duct1', Duct())
        duct1.Q_dot = 0  # no heat exchangers
        duct1.dPqP = .1  #no losses

        split = self.add('split', SplitterW())
        split.W1_des = .44
        split.MNexit1_des = 1.0
        split.MNexit2_des = 1.0

        nozzle = self.add('nozzle', Nozzle())
        nozzle.dPqP = 0  #no losses
        nozzle.Fl_ref = FlowStation()

        comp2 = self.add('comp2', Compressor())
        comp2.PR_des = 5.0
        comp2.MNexit_des = .4
        comp2.eff_des = .80

        duct2 = self.add('duct2', Duct())  #to bearings
        duct2.Q_dot = 0  #no heat exchangers
        duct2.dPqP = .1  #no losses

        perf = self.add('perf', Performance())

        #Inter Component Connections
        self.connect('tube.Fl_O', 'inlet.Fl_I')
        self.connect('inlet.Fl_O', 'comp1.Fl_I')
        self.connect('comp1.Fl_O', 'duct1.Fl_I')
        self.connect('duct1.Fl_O', 'split.Fl_I')
        self.connect('split.Fl_O2', 'nozzle.Fl_I')
        self.connect('tube.Fl_O', 'nozzle.Fl_ref')
        self.connect('split.Fl_O1', 'comp2.Fl_I')
        self.connect('comp2.Fl_O', 'duct2.Fl_I')
        self.connect('comp1.pwr', 'perf.C1_pwr')
        self.connect('comp2.pwr', 'perf.C2_pwr')
        self.connect('duct2.Fl_O.Ps', 'perf.Ps_bearing')
        self.connect('nozzle.Fg', 'perf.Fg')
        self.connect('inlet.F_ram', 'perf.F_ram')

        #Input variable pass_throughs to the assembly boundary
        #Compress -> Tube
        self.connect('W_in', 'tube.W')
        self.connect('Ts_tube', 'tube.Ts')
        self.connect('Ps_tube', 'tube.Ps')
        self.connect('Mach_pod_max', 'tube.Mach')
        #Compress -> Inlet
        self.connect('Mach_c1_in', 'inlet.MNexit_des')
        #Compress -> C1
        self.connect('c1_PR_des', 'comp1.PR_des')
        #Compress -> C2
        self.connect('c2_PR_des', 'comp2.PR_des')
        #Compress -> Splitter
        self.connect('W_bearing_in', 'split.W1_des')
        #Compress -> Perf
        self.connect('Ps_bearing', 'perf.Ps_bearing_target')

        #Output variable pass_throughs to the assembly boundary
        self.connect('tube.Fl_O.rhot', 'rho_air')  #promoted for aero calc
        self.connect('tube.Fl_O.area', 'area_inlet_in')
        self.connect('tube.Fl_O.Vflow', 'speed_max')
        self.connect('inlet.Fl_O.area', 'area_c1_in')
        self.connect('comp1.Fl_O.area', 'area_c1_out')
        self.connect('nozzle.Fl_O.area', 'nozzle_flow_area')
        self.connect('nozzle.Fl_O', 'nozzle_Fl_O')
        self.connect('duct2.Fl_O', 'bearing_Fl_O')
        self.connect('perf.F_net', 'F_net')
        self.connect('perf.pwr', 'pwr_req')
        self.connect('perf.Ps_bearing_residual', 'Ps_bearing_residual')

        #driver setup
        design = self.driver
        comp_list = [
            'tube', 'inlet', 'comp1', 'duct1', 'split', 'nozzle', 'comp2',
            'duct2', 'perf'
        ]

        design.workflow.add(comp_list)
        for comp_name in comp_list:  #need to put everything in design mode
            design.add_event('%s.design' % comp_name)