コード例 #1
0
 def test_add_unit(self):
     try:
         units.add_unit('ft', '20*m')
     except KeyError, err:
         self.assertEqual(
             str(err),
             "'Unit ft already defined with different factor or powers'")
コード例 #2
0
def add_units():
    """ Add units we need if not already defined. """

    required = (
        ('lbf',  '4.4482216152605*N', 'Pounds force'),
        ('lbm',  '1.*lb',             'Pounds mass'),
        ('rpm',  '6.*deg/s',          'Revolutions per minute.'),
        ('slug', '14.5939*kg',        'Slug'),
    )
    for name, unit, comment in required:
        try:
            units.PhysicalQuantity(0., name)
        except ValueError:
            units.add_unit(name, unit, comment)
コード例 #3
0
 def test_add_unit(self):
     try:
         units.add_unit('ft','20*m')
     except KeyError,err: 
         self.assertEqual(str(err),"'Unit ft already defined with different factor or powers'")
コード例 #4
0
class HSRNOISE(ExternalCode):
    """OpenMDAO component wrapper for HSRNOISE."""

    add_unit('dB', 'kg**0', 'decibel')

    # Variables from MEflows and Geometry variable trees
    # -------------------------
    flow_in = VarTree(MEflows(), iotype='in')
    geo_in = VarTree(Geometry(), iotype='in')

    LinFrac = Float(
        0.9,
        iotype='in',
        desc='Fraction of ejector length covered by acoustic liner')
    phi = Float(0.0,
                iotype='in',
                units='deg',
                desc='Roll angle at which noise is estimated')

    # Variables GEOM namelist
    # -----------------------------
    HMIC = Float(4.0,
                 iotype='in',
                 units='ft',
                 desc='Height of microphone above ground level')
    SL = Float(0.0, iotype='in', units='ft', desc='Sideline distance')

    # Variables FLIPATH namelist
    # -----------------------------
    ALTEVO = Float(700.0,
                   iotype='in',
                   units='ft',
                   desc='Aircraft altitude overhead')
    FPA = Float(0.0,
                iotype='in',
                units='deg',
                desc='Flight path angle, positive for climb')
    PAE = Float(
        0.0,
        iotype='in',
        units='deg',
        desc=
        'Angle between engine axis and horizontal direction, positive for inlet up'
    )
    VAIR = Float(0.0, iotype='in', units='ft/s', desc='Aircraft velocity')

    # Variables SOURCE namelist
    # -----------------------------
    JETMETHOD = Enum(
        1, [1, 2],
        iotype='in',
        desc='Jet method flag; 1=Boeing JN8C4 Method, 2=Stone Method')
    NENG = Int(1, iotype='in', desc='Number of engines')

    # Variables for JET1IN namelist
    # -----------------------------
    #APRI = Float(8.65620, iotype='in', units='ft**2', desc='Primary nozzle exit area')
    #ASEC = Float(16.832, iotype='in', units='ft**2', desc='Secondary nozzle exit area')
    ATHP = Float(8.5466,
                 iotype='in',
                 units='ft**2',
                 desc='Throat area of the primary nozzle')
    CFG = Float(0.95, iotype='in', desc='Mixer-ejector thrust coefficient')
    DELMIX = Float(0, iotype='in', units='dB', desc='Mixed jet suppression')
    DELPE = Float(0,
                  iotype='in',
                  units='dB',
                  desc='Primary external jet suppression')
    DELPI = Float(0,
                  iotype='in',
                  units='dB',
                  desc='Primary internal jet suppression')
    #DIVANG = Float(0, iotype='in', units='deg', desc='Spoke divergence angle. Outer tangent line relative to ejector flow surface at nozzle downstream of mixer exit')
    #EJASP = Float(1.5, iotype='in', desc='Mixer-ejector aspect ratio major/minor axis')
    EJD = Float(
        5.274,
        iotype='in',
        units='ft',
        desc='Equivalent inner diameter of the ejector at the mixer exit')
    #EJL = Float(13.3, iotype='in', units='ft', desc='Ejector length')
    #EJLIN = Float(9.0, iotype='in', units='ft', desc='Length of ejector lining')
    FLIN = Float(2000,
                 iotype='in',
                 units='Hz',
                 desc='Liner design center frequency')
    #PEN = Float(0.925, iotype='in', desc='Mixer spoke penetration to nozzle full height ratio')
    PLUGD = Float(0,
                  iotype='in',
                  units='ft',
                  desc='Plug diameter at nozzle exit station')
    PSI0 = Float(
        0,
        iotype='in',
        units='deg',
        desc=
        'Rectangular exit plane major axis roll angle relative to wing plane')
    SPOKE = Float(18.0,
                  iotype='in',
                  desc='Number of spokes or lobes of the mixer')
    #TPRI = Float(1553.0, iotype='in', units='degR', desc='Primary jet total temperature')
    #TSEC = Float(529.0, iotype='in', units='degR', desc='Secondary jet total temperature')
    #VPRI = Float(2378.0, iotype='in', units='ft/s', desc='Primary jet fully expanded velocity')
    #VSEC = Float(449.0, iotype='in', units='ft/s', desc='Secondary jet velocity')
    #WPRI = Float(811.3, iotype='in', units='lbm/s', desc='Primary jet mass flow rate')
    #WSEC = Float(566.9, iotype='in', units='lbm/s', desc='Secondary jet mass flow rate')
    #XMAR = Float(0.97, iotype='in', desc='Ejector exit area to mixing plane area ratio')

    # Variables for JET2IN namelist
    # -----------------------------
    APT = Float(8.5466,
                iotype='in',
                units='ft**2',
                desc='Throat area of the primary nozzle')
    AS = Float(16.832,
               iotype='in',
               units='ft**2',
               desc='Secondary nozzle flow area')
    CER = Float(iotype='in', desc='Core expansion ratio')
    DHP = Float(
        0.5,
        iotype='in',
        units='ft',
        desc='Hydrolic diameter of primary flow(DHP = 4*Area/Perimeter)')
    DL = Float(0.7, iotype='dB/ft', desc='Peak perforate suppression at FPK')
    DM = Float(0,
               iotype='in',
               desc='Degree of mixing (calculated if not supplied)')
    FPK = Float(1995,
                iotype='in',
                units='Hz',
                desc='Peak frequency of broad-band suppression')
    GAMMAC = Float(
        0.,
        iotype='in',
        desc=
        'Specific heat ratio of primary stream. Internally calculated if set to 0'
    )
    HEX = Float(4.0, iotype='in', units='ft', desc='Ejector height at nozzle')
    IEX = Enum(
        0, [0, 1],
        iotype='in',
        desc=
        'Mixed jet parameter flag for TEX and VEX; 0=calculate internally, 1=use input values'
    )
    ISUPPR = Enum(
        0, [0, 1],
        iotype='in',
        desc=
        'Suppression flag; 0=negate suppression of internal noise sources, 1=suppress internal noise sources'
    )
    LBE = Float(10.0,
                iotype='in',
                units='ft',
                desc='Length to end of bulk treatment from primary exit')
    LBS = Float(0.0,
                iotype='in',
                units='ft',
                desc='Length to start of bulk treatment from primary exit')
    LE = Float(10.0, iotype='in', units='ft', desc='Ejector length')
    LPE = Float(
        0.0,
        iotype='in',
        units='ft',
        desc=
        'Length to end of perforate treatment from primary exit. For perforate over bultk treatment, add bulk length to perforate'
    )
    LPS = Float(
        0.0,
        iotype='in',
        units='ft',
        desc='Length to start of perforate treatment from primary exit')
    MMC = Float(0.0,
                iotype='in',
                desc='Molecular mass (weight) of primary stream')
    MPD = Float(1.2, iotype='in', desc='Design Mach number of primary nozzle')
    PC = Float(iotype='in',
               units='lbf/ft**2',
               desc='Total pressure of primary stream upstream of nozzle')
    PEN = Float(0.925,
                iotype='in',
                desc='Mixer spoke penetration to nozzle full height ratio')
    SAR = Float(iotype='in',
                desc='Ratio of total mixing area to primary nozzle throat')
    SUPPK = Float(1.6628,
                  iotype='in',
                  units='dB/ft',
                  desc='Peak bulk suppression')
    TC = Float(iotype='in',
               units='degR',
               desc='Total temperature of primary stream upstream of nozzle')
    TEX = Float(iotype='in', units='degR', desc='Mixed exit total temperature')
    VEX = Float(iotype='in', units='ft/s', desc='Mixed exit velocity')
    WEX = Float(6.0, units='ft', desc='Ejector width at nozzle exit')
    WSWP = Float(iotype='in', desc='Pumping ratio')

    # Output variables
    # -----------------------------
    TotalEPNL = Float(iotype='out',
                      units='dB',
                      desc='EPNL for the total aircraft')
    TotalMaxPNLT = Float(iotype='out',
                         units='dB',
                         desc='Maximum PNLT for the total aircraft')
    JetEPNL = Float(iotype='out', units='dB', desc='EPNL from the jet only')
    OASPL = Array(zeros((0, 17)),
                  iotype='out',
                  desc='OASPL values for total jet noise')
    OASPL30 = Float(iotype='out',
                    units='dB',
                    desc='OASPL from the jet only at emission angle 30')
    OASPL60 = Float(iotype='out',
                    units='dB',
                    desc='OASPL from the jet only at emission angle 60')
    OASPL90 = Float(iotype='out',
                    units='dB',
                    desc='OASPL from the jet only at emission angle 90')
    OASPL120 = Float(iotype='out',
                     units='dB',
                     desc='OASPL from the jet only at emission angle 120')
    OASPL150 = Float(iotype='out',
                     units='dB',
                     desc='OASPL from the jet only at emission angle 150')
    thetas = Array(zeros((0, 17)),
                   dtype=numpy_float,
                   iotype='out',
                   desc='Yaw angles')
    corr = Array(zeros((0, 17)),
                 dtype=numpy_float,
                 iotype='out',
                 desc='Correction for spherical spreading')
    SPL = Array(zeros((24, 17)),
                dtype=numpy_float,
                iotype='out',
                desc='SPL values produced by HSRNoise')
    SPL_corr = Array(
        zeros((24, 17)),
        dtype=numpy_float,
        iotype='out',
        desc='SPL values after being corrected back to a 1 ft lossless arc')
    Freq = Array(zeros((24, 0)),
                 dtype=numpy_float,
                 iotype='out',
                 desc='Frequencies')

    def __init__(self):
        super(HSRNOISE, self).__init__()
        self.command = ['hsrnoise', 'test.input', 'test.output']

        self.add('geo_in', Geometry())
        self.add('flow_in', MEflows())

    def execute(self):
        """ Executes our file-wrapped component. """

        self.ATHP = self.geo_in.Apri / (
            ((self.flow_in.gamma + 1) / 2)**((-self.flow_in.gamma - 1) /
                                             (2 * (self.flow_in.gamma - 1))) *
            (1 + (self.flow_in.gamma - 1) / 2 * self.flow_in.pri.Mach**2)**(
                (self.flow_in.gamma + 1) /
                (2 * (self.flow_in.gamma - 1))) / self.flow_in.pri.Mach)
        self.EJD = 2 * (self.geo_in.Aexit / pi)**0.5
        self.SPOKE = self.geo_in.Num_Lobes
        self.HMIC = self.ALTEVO - 10 * cos(radians(self.phi))
        self.SL = 10 * sin(radians(self.phi))

        #Prepare the input file for HSRNOISE
        self.generate_input()

        #Execute the component
        #super(HSRNOISE, self).execute()

        #Parse HSRNOISE output file
        self.parse_output()

        #Delete .plt files produced by HSRnoise
        for filename in glob('*.plt'):
            os.remove(filename)

        #Post process the data produced by HSRNoise for use in ANOPP
        self.dist = sqrt((self.ALTEVO - self.HMIC)**2 + self.SL**2) / sin(
            radians(self.thetas))
        self.corr = 20 * log10(self.dist)
        self.SPL_corr = self.SPL + self.corr

    def generate_input(self):
        """Creates the HSRNOISE input file."""

        self.EJLIN = self.LinFrac * self.geo_in.length

        parser = InputFileGenerator()
        parser.set_template_file('test_temp.input')
        parser.set_generated_file('test.input')

        # Set Geometry and Flight Conditions
        # --------------------------------------
        parser.mark_anchor("$GEOM")
        parser.transfer_var(self.HMIC, 1, 3)
        parser.transfer_var(self.SL, 2, 3)
        parser.mark_anchor("$FLIPATH")
        parser.transfer_var(self.ALTEVO, 1, 3)

        parser.mark_anchor("JETMETHOD")
        parser.transfer_var(self.JETMETHOD, 0, 3)

        # Replace JET1IN variables in test.input
        # --------------------------------------
        parser.mark_anchor("$JET1IN")
        parser.transfer_var(self.geo_in.Apri, 1, 3)  #APRI
        parser.transfer_var(self.geo_in.Asec, 2, 3)  #ASEC
        parser.transfer_var(self.ATHP, 3, 3)
        parser.transfer_var(self.CFG, 4, 3)
        parser.transfer_var(self.DELMIX, 5, 3)
        parser.transfer_var(self.DELPE, 6, 3)
        parser.transfer_var(self.DELPI, 7, 3)
        parser.transfer_var(self.geo_in.ChuteAngles, 8, 3)  #DIVANG
        parser.transfer_var(self.geo_in.AR, 9, 3)  #EJASP
        parser.transfer_var(self.EJD, 10, 3)
        parser.transfer_var(self.geo_in.length, 11, 3)  # EJL
        parser.transfer_var(self.EJLIN, 12, 3)
        parser.transfer_var(self.FLIN, 13, 3)
        parser.transfer_var(self.geo_in.LhMh, 14, 3)  #PEN
        parser.transfer_var(self.PLUGD, 15, 3)
        parser.transfer_var(self.PSI0, 16, 3)
        parser.transfer_var(self.SPOKE, 17, 3)
        parser.transfer_var(self.flow_in.pri.Tt, 18, 3)  #TPRI
        parser.transfer_var(self.flow_in.sec.Tt, 19, 3)  #TSEC
        parser.transfer_var(self.flow_in.pri.Vel, 20, 3)  #VPRI
        parser.transfer_var(self.flow_in.sec.Vel, 21, 3)  #VSEC
        parser.transfer_var(self.flow_in.pri.W, 22, 3)  #WPRI
        parser.transfer_var(self.flow_in.sec.W, 23, 3)  #WSEC
        parser.transfer_var(self.geo_in.AeAt, 24, 3)  #XMAR

        # Replace JET2IN variables in test.input
        # --------------------------------------
        parser.mark_anchor("$JET2IN")
        parser.transfer_var(self.APT, 1, 3)
        parser.transfer_var(self.AS, 2, 3)
        parser.transfer_var(self.CER, 3, 3)
        parser.transfer_var(self.DHP, 4, 3)
        parser.transfer_var(self.DL, 5, 3)
        parser.transfer_var(self.DM, 6, 3)
        parser.transfer_var(self.FPK, 7, 3)
        parser.transfer_var(self.GAMMAC, 8, 3)
        parser.transfer_var(self.HEX, 9, 3)
        parser.transfer_var(self.IEX, 10, 3)
        parser.transfer_var(self.ISUPPR, 11, 3)
        parser.transfer_var(self.LBE, 12, 3)
        parser.transfer_var(self.LBS, 13, 3)
        parser.transfer_var(self.LE, 14, 3)
        parser.transfer_var(self.LPE, 15, 3)
        parser.transfer_var(self.LPS, 16, 3)
        parser.transfer_var(self.MMC, 17, 3)
        parser.transfer_var(self.MPD, 18, 3)
        parser.transfer_var(self.PC, 19, 3)
        parser.transfer_var(self.PEN, 20, 3)
        parser.transfer_var(self.SAR, 21, 3)
        parser.transfer_var(self.SUPPK, 22, 3)
        parser.transfer_var(self.TC, 23, 3)
        parser.transfer_var(self.TEX, 24, 3)
        parser.transfer_var(self.VEX, 25, 3)
        parser.transfer_var(self.WEX, 26, 3)
        parser.transfer_var(self.WSWP, 27, 3)

        parser.generate()

    def parse_output(self):
        """Parses the HSRNOISE output file and extracts data."""

        outfile = FileParser()
        outfile.set_file('test.output')

        outfile.mark_anchor("JN8C4 JET NOISE MODULE")
        outfile.mark_anchor("TOTAL")
        self.thetas = outfile.transfer_array(4, 2, 4, 18)
        self.Freq = outfile.transfer_2Darray(7, 1, 30, 1)
        self.SPL = outfile.transfer_2Darray(7, 2, 30, 18)
        outfile.mark_anchor("DBA")
        self.OASPL = outfile.transfer_array(-1, 2, -1, 18)
        self.OASPL30 = outfile.transfer_var(-1, 4)
        self.OASPL60 = outfile.transfer_var(-1, 7)
        self.OASPL90 = outfile.transfer_var(-1, 10)
        self.OASPL120 = outfile.transfer_var(-1, 13)
        self.OASPL150 = outfile.transfer_var(-1, 16)

        outfile.mark_anchor("EPNL SUMMARY")
        self.TotalEPNL = outfile.transfer_var(9, 2)
        self.TotalMaxPNLT = outfile.transfer_var(9, 5)
        self.JetEPNL = outfile.transfer_var(8, 2)

    def load_model(self):
        """Loads an existing HSRNOISE input file."""