/contrib/drea_usecase/DREA_HSRnoise.py

https://github.com/thearn/OpenMDAO-Framework · Python · 131 lines · 94 code · 23 blank · 14 comment · 1 complexity · c644f21680057c67c7a21539e3d6b72a MD5 · raw file 

    

  1. """

    2. 

  2.     DREA_HSRnoise.py - Code to run a DREA followed by HSRnoise for M-E off-design analysis

    3. 

  3. """

    4. 

  4. 

    5. 

  5. from openmdao.main.api import Assembly, Slot, Case, Instance

    6. 

  6. from openmdao.main.datatypes.api import Float, Enum

    7. 

  7. from openmdao.lib.casehandlers.api import ListCaseIterator, DBCaseRecorder

    8. 

  8. from openmdao.lib.drivers.api import CaseIteratorDriver

    9. 

  9. #from openmdao.lib.drivers.api import DOEdriver

    10. 

  10. #from openmdao.lib.doegenerators.api import FullFactorial

    11. 

  11. from DREA import DREA

    12. 

  12. from hsrnoise import HSRNOISE

    13. 

  13. from DREAprep import DREAprep

    14. 

  14. from geometry import Geometry

    15. 

  15. from ACDgen import ACDgen

    16. 

  16. 

    17. 

  17. class DREA_HSRnoise(Assembly):

    18. 

  18.     """Assembly to execute on DREA followed by HSRnoise."""

    19. 

  19.     

    20. 

  20.     geo = Instance(Geometry, iotype='in')

    21. 

  21.     alt = Float(0.0, iotype='in', units='ft', desc='Altitude')

    22. 

  22.     point = Enum(1, [1,2,3],  iotype='in', desc='Certification observer point')

    23. 

  23.     

    24. 

  24.     def __init__(self):

    25. 

  25.         """Creates an Assembly to run DREA and HSRnoise."""

    26. 

  26.         super(DREA_HSRnoise, self).__init__()

    27. 

  27.         

    28. 

  28.         FO1 = Case(inputs=[('point', 1),('dreaprep.Mach',0.28),('alt',2000.0),('dreaprep.PC',100.0),('hsrnoise.phi', 0.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    29. 

  29.         FO2 = Case(inputs=[('point', 1),('dreaprep.Mach',0.28),('alt',2000.0),('dreaprep.PC', 65.0),('hsrnoise.phi', 0.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    30. 

  30.         App = Case(inputs=[('point', 2),('dreaprep.Mach',0.20),('alt', 394.0),('dreaprep.PC', 30.0),('hsrnoise.phi', 0.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    31. 

  31.         SL1 = Case(inputs=[('point', 3),('dreaprep.Mach',0.25),('alt',1000.0),('dreaprep.PC',100.0),('hsrnoise.phi', 0.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    32. 

  32.         SL2 = Case(inputs=[('point', 3),('dreaprep.Mach',0.25),('alt',1000.0),('dreaprep.PC',100.0),('hsrnoise.phi',30.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    33. 

  33.         SL3 = Case(inputs=[('point', 3),('dreaprep.Mach',0.25),('alt',1000.0),('dreaprep.PC',100.0),('hsrnoise.phi',60.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    34. 

  34.         SL4 = Case(inputs=[('point', 3),('dreaprep.Mach',0.25),('alt',1000.0),('dreaprep.PC',100.0),('hsrnoise.phi',90.0)], outputs=[('drea.CFG',0.0),('hsrnoise.thetas',0.0),('hsrnoise.Freq',0.0),('hsrnoise.SPL_corr',0),('hsrnoise.OASPL30',0.0),('hsrnoise.OASPL60',0.0),('hsrnoise.OASPL90',0.0),('hsrnoise.OASPL120',0.0),('hsrnoise.OASPL150',0.0)])

    35. 

  35.         

    36. 

  36.         cases = ListCaseIterator([FO1,FO2,App,SL1,SL2,SL3,SL4])

    37. 

  37.         db_recorder = DBCaseRecorder()

    38. 

  38.         

    39. 

  39.         self.add('geo', Geometry())

    40. 

  40.         

    41. 

  41.         self.add('dreaprep', DREAprep())

    42. 

  42.         self.add('drea', DREA())

    43. 

  43.         self.add('hsrnoise', HSRNOISE())

    44. 

  44.         self.add('ACDgen', ACDgen())

    45. 

  45.         

    46. 

  46.         self.add('analysis',CaseIteratorDriver())

    47. 

  47.         self.analysis.iterator = cases

    48. 

  48.         self.analysis.recorders = [db_recorder]

    49. 

  49.         self.ACDgen.case_data = db_recorder.get_iterator()

    50. 

  50.         

    51. 

  51.         # Set up the workflows

    52. 

  52.         #---------------------------

    53. 

  53.         #self.analysis.workflow.add(['dreaprep', 'drea', 'hsrnoise'])

    54. 

  54.         #self.driver.workflow.add(['analysis','ACDgen'])

    55. 

  55.         

    56. 

  56.         self.driver.workflow.add(['dreaprep', 'drea', 'hsrnoise'])

    57. 

  57.                 

    58. 

  58.         # Connections

    59. 

  59.         #---------------------------

    60. 

  60.         self.connect('geo',['drea.geo_in','hsrnoise.geo_in'])

    61. 

  61.         self.connect('alt',['dreaprep.alt','hsrnoise.ALTEVO'])

    62. 

  62.         self.connect('dreaprep.flow_out','drea.flow_in')

    63. 

  63.         self.connect('drea.flow_out','hsrnoise.flow_in')

    64. 

  64.         self.connect('drea.CFG','hsrnoise.CFG')

    65. 

  65.         

    66. 

  66. if __name__ == "__main__":

    67. 

  67.     asy = DREA_HSRnoise()

    68. 

  68. 

    69. 

  69.     import time

    70. 

  70.     tt = time.time()

    71. 

  71.     asy.drea.mode = 'Subsonic'

    72. 

  72.     asy.dreaprep.dTs = 18.3

    73. 

  73.     asy.dreaprep.Mach = 0.0

    74. 

  74.     asy.alt = 0.0

    75. 

  75.     asy.hsrnoise.phi = 90.0

    76. 

  76.     asy.dreaprep.PC = 80

    77. 

  77.     asy.dreaprep.deltaPt = 0.0

    78. 

  78.     asy.dreaprep.deltaTt = 0.0

    79. 

  79.     asy.dreaprep.deltaM = 0.0

    80. 

  80. 

    81. 

  81.     asy.geo.Apri = 12.0

    82. 

  82.     asy.geo.AsAp = 1.75

    83. 

  83.     asy.geo.AR = 3.0 

    84. 

  84.     asy.geo.AeAt = .95

    85. 

  85.     asy.geo.length = 8.0 

    86. 

  86.     asy.geo.ChuteAngles = 10

    87. 

  87.     asy.geo.calc_geom(asy.geo.length,asy.geo.Apri,asy.geo.AsAp,asy.geo.AR,asy.geo.AeAt,asy.geo.LhMh,asy.geo.LhWave)

    88. 

  88. 

    89. 

  89.     

    90. 

  90.     asy.run()

    91. 

  91.     print '\n\n'

    92. 

  92.     print 'alt             : ',asy.dreaprep.alt

    93. 

  93.     print 'Mach            : ',asy.dreaprep.Mach

    94. 

  94.     print 'HMIC            : ',asy.hsrnoise.HMIC

    95. 

  95.     print 'SL              : ',asy.hsrnoise.SL

    96. 

  96.     print '\n'

    97. 

  97.     print 'INPUTS            DREA       HSRNoise'

    98. 

  98.     print 'Length          : ',asy.drea.geo_in.length, '     ', asy.hsrnoise.geo_in.length

    99. 

  99.     print 'Width           : ',asy.drea.geo_in.width, '     ', asy.hsrnoise.geo_in.width

    100. 

  100.     print 'Primary Area    : ',asy.drea.geo_in.Apri, '     ', asy.hsrnoise.geo_in.Apri

    101. 

  101.     print 'Secondary Area  : ',asy.drea.geo_in.Asec, '     ', asy.hsrnoise.geo_in.Asec

    102. 

  102.     print 'Exit Area       : ',asy.drea.geo_in.Aexit, '     ', asy.hsrnoise.geo_in.Aexit

    103. 

  103.     print 'Primary Pressure: ',asy.drea.flow_in.pri.Pt,  '     ', asy.hsrnoise.flow_in.pri.Pt

    104. 

  104.     print 'Primary Temp    : ',asy.drea.flow_in.pri.Tt,  '     ', asy.hsrnoise.flow_in.pri.Tt

    105. 

  105.     print 'Primary Mach    : ',asy.drea.flow_in.pri.Mach,  '     ', asy.hsrnoise.flow_in.pri.Mach

    106. 

  106.     print '\n'

    107. 

  107.     print 'RESULTS'

    108. 

  108.     print 'Gross Thrust    : ',asy.drea.GrossThrust

    109. 

  109.     print 'Exit Velocity   : ',asy.drea.ExitVelocity

    110. 

  110.     print 'Exit Mach       : ',asy.drea.ExitMach

    111. 

  111.     print 'CFG:            : ',asy.drea.CFG

    112. 

  112.     print 'Pumping Ratio   : ',asy.drea.PumpingRatio

    113. 

  113.     print 'Mass Flow (Pri) : ',asy.drea.flow_out.pri.W

    114. 

  114.     print 'Mass Flow (Sec) : ',asy.drea.flow_out.sec.W

    115. 

  115.     print 'Secondary Mach  : ',asy.drea.flow_out.sec.Mach

    116. 

  116.     print 'NPR             : ',asy.drea.NPR

    117. 

  117.     print 'Degree of mixing: ',asy.drea.DegreeOfMixing

    118. 

  118.     print 'OASPL30         : ',asy.hsrnoise.OASPL30

    119. 

  119.     print 'OASPL60         : ',asy.hsrnoise.OASPL60

    120. 

  120.     print 'OASPL90         : ',asy.hsrnoise.OASPL90

    121. 

  121.     print 'OASPL120        : ',asy.hsrnoise.OASPL120

    122. 

  122.     print 'OASPL150        : ',asy.hsrnoise.OASPL150

    123. 

  123.     #print 'test: ', asy.hsrnoise.SPL_corr

    124. 

  124.     print "\n"

    125. 

  125.     print "Elapsed time: ", time.time()-tt, "seconds"
    126.