Level 1#

What is level 1#

This part is a tutorial type introduction. We assume users have read Aviary User Interface. In this doc page we discuss Level 1’s interface, but also other details, such as inputs and outputs.

If you have not yet, make sure to install Aviary.

To run Aviary in the Level 1 command-line interface (CLI), type this in a terminal:

!aviary
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
usage: aviary [-h] [--version]  ...

aviary Command Line Tools

options:
  -h, --help            show this help message and exit
  --version             show version and exit

Tools:
  
    convert_aero_table  Converts FLOPS- or GASP-formatted aero data files into
                        Aviary csv format.
    convert_engine      Converts FLOPS- or GASP-formatted engine decks into
                        Aviary csv format. FLOPS decks are changed from
                        column-delimited to csv format with added headers.
                        GASP decks are reorganized into column based csv. T4
                        is recovered through calculation. Data points whose T4
                        exceeds T4max are removed.
    convert_prop_table  Converts GASP-formatted propeller map file into Aviary
                        csv format.
    dashboard           Run the Dashboard tool
    draw_mission        Allows users to draw a mission profile for use in
                        Aviary.
    fortran_to_aviary   Converts legacy Fortran input decks to Aviary csv
                        based decks
    hangar              Allows users that pip installed Aviary to download
                        models from the Aviary hangar
    plot_drag_polar     Plot a Drag Polar Graph using a provided polar data
                        csv input
    run_mission         Runs Aviary using a provided input deck

We’ll now discuss the different parts of the CLI, starting with how to run an Aviary problem.

First level 1 run#

We are ready to run some models. By default, Aviary runs the optimizer for up to 50 iterations. In order to reduce the run time, we will limit the maximum number of iterations to 1. It will be faster for these examples, but you will not get optimal solutions. Issue the following command:

!aviary run_mission models/test_aircraft/aircraft_for_bench_GwGm.csv --optimizer IPOPT --max_iter 1
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
Loaded default phase_info for 2dof equations of motion
User has specified Design.NUM_PASSENGERS but CrewPayload.NUM_PASSENGERS has been left blank or set to zero.
Assuming they are equal to maintain backwards compatibility with GASP and FLOPS output files.
If you intended to have no passengers on this flight, please set Aircraft.CrewPayload.TOTAL_PAYLOAD_MASS to zero in aviary_values.
/home/runner/work/Aviary/Aviary/aviary/utils/preprocessors.py:375: UserWarning: Mount location for engines of type <turbofan_23k_1> not specified. Wing-mounted engines are assumed.
  warnings.warn(
The following variables have been overridden:
  'aircraft:fuselage:wetted_area  4000  ft**2
  'aircraft:horizontal_tail:form_factor  1.25  unitless
  'aircraft:nacelle:form_factor  [1.5]  unitless
  'aircraft:strut:fuselage_interference_factor  0  unitless
  'aircraft:vertical_tail:form_factor  1.25  unitless
  'aircraft:wing:form_factor  1.25  unitless
--- Constraint Report [traj] ---
    --- groundroll ---
        None
    --- rotation ---
        [final]   0.0000e+00 == normal_force [lbf]
    --- ascent ---
        [final]   5.0000e+02 == altitude [ft]
        [path]    0.0000e+00 <= load_factor <= 1.1000e+00  [unitless]
        [path]    0.0000e+00 <= fuselage_pitch <= 1.5000e+01  [deg]
    --- accel ---
        [final]   2.5000e+02 == EAS [kn]
    --- climb1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- climb2 ---
        [final]   3.7500e+04 == altitude [ft]
        [final]   1.0000e-01 <= altitude_rate  [ft/min]
        [final]   8.0000e-01 == mach [unitless]
    --- cruise ---
        None
    --- desc1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- desc2 ---
        [final]   1.0000e+03 == altitude [ft]
'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.ascent' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.


'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.climb1' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.


'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.desc1' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.
Full total jacobian for problem 'aircraft_for_bench_GwGm' was computed 3 times, taking 1.8758711499999663 seconds.
Total jacobian shape: (238, 212) 
Jacobian shape: (238, 212)  (8.73% nonzero)
FWD solves: 51   REV solves: 0
Total colors vs. total size: 51 vs 212  (75.94% improvement)

Sparsity computed using tolerance: 1e-25
Time to compute sparsity:   1.8759 sec
Time to compute coloring:   0.3298 sec
Memory to compute coloring:   1.7500 MB
Coloring created on: 2024-12-17 21:05:22
This is Ipopt version 3.14.17, running with linear solver MUMPS 5.7.3.

Number of nonzeros in equality constraint Jacobian...:     7490
Number of nonzeros in inequality constraint Jacobian.:     1592
Number of nonzeros in Lagrangian Hessian.............:        0
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/total_jac.py:1728: DerivativesWarning:Constraints or objectives [('traj.phases.climb2->final_boundary_constraint->mach', inds=[(11, 0)])] cannot be impacted by the design variables of the problem.
Total number of variables............................:      212
                     variables with only lower bounds:        2
                variables with lower and upper bounds:      210
                     variables with only upper bounds:        0
Total number of equality constraints.................:      204
Total number of inequality constraints...............:       33
        inequality constraints with only lower bounds:        1
   inequality constraints with lower and upper bounds:       32
        inequality constraints with only upper bounds:        0

iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls
   0  5.2731333e+00 1.74e+03 1.65e+01  -5.0 0.00e+00    -  0.00e+00 0.00e+00   0
Number of Iterations....: 1

                                   (scaled)                 (unscaled)
Objective...............:   4.4308077759954099e+00    4.4308077759954099e+00
Dual infeasibility......:   1.6268337948231213e+01    1.6268337948231213e+01
Constraint violation....:   5.3436691787363870e+01    3.0963174509844049e+02
Variable bound violation:   0.0000000000000000e+00    0.0000000000000000e+00
Complementarity.........:   1.0000016687732632e+05    1.0000016687732632e+05
Overall NLP error.......:   1.0000016687732632e+05    1.0000016687732632e+05


Number of objective function evaluations             = 2
Number of objective gradient evaluations             = 2
Number of equality constraint evaluations            = 2
Number of inequality constraint evaluations          = 2
Number of equality constraint Jacobian evaluations   = 2
Number of inequality constraint Jacobian evaluations = 2
Number of Lagrangian Hessian evaluations             = 0
Total seconds in IPOPT                               = 2.439

EXIT: Maximum Number of Iterations Exceeded.
minimal_print is not available for this solution


Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
    Objective Function: _objfunc

    Solution: 
--------------------------------------------------------------------------------
    Total Time:                    2.4406
       User Objective Time :       0.7056
       User Sensitivity Time :     1.6981
       Interface Time :            0.0235
       Opt Solver Time:            0.0134
    Calls to Objective Function :       3
    Calls to Sens Function :            3


   Objectives
      Index  Name                               Value
          0  mission:objectives:fuel     4.430808E+00

   Variables (c - continuous, i - integer, d - discrete)
      Index  Name                                      Type      Lower Bound            Value      Upper Bound     Status
          0  mission:design:gross_mass_0                  c     5.714286E-05     9.943768E-01     1.000000E+30           
          1  mission:summary:gross_mass_0                 c     5.714286E-05     9.943769E-01     1.000000E+30           
          2  mission:takeoff:ascent_t_initial_0           c     0.000000E+00     5.650347E-01     3.333333E+00           
          3  mission:takeoff:ascent_duration_0            c     1.000000E-01     5.388178E-01     1.000000E+02           
          4  tau_gear_0                                   c     1.000000E-02     2.003476E-01     1.000000E+00           
          5  tau_flaps_0                                  c     1.000000E-02     9.687929E-01     1.000000E+00           
          6  traj.accel.t_duration_0                      c     9.950249E-03     2.932593E-01     1.990050E+00           
          7  traj.accel.states:velocity_0                 c     0.000000E+00     5.094015E-01     1.200000E+00           
          8  traj.accel.states:velocity_1                 c     0.000000E+00     8.991905E-01     1.200000E+00           
          9  traj.accel.states:velocity_2                 c     0.000000E+00     1.018243E+00     1.200000E+00           
         10  traj.accel.states:mass_0                     c     0.000000E+00     1.157638E+00     6.666667E+15           
         11  traj.accel.states:mass_1                     c     0.000000E+00     1.157381E+00     6.666667E+15           
         12  traj.accel.states:mass_2                     c     0.000000E+00     1.157022E+00     6.666667E+15           
         13  traj.accel.states:mass_3                     c     0.000000E+00     1.156907E+00     6.666667E+15           
         14  traj.accel.states:distance_0                 c     0.000000E+00     3.426977E-01     3.000000E+01           
         15  traj.accel.states:distance_1                 c     0.000000E+00     5.211404E-01     3.000000E+01           
         16  traj.accel.states:distance_2                 c     0.000000E+00     5.834464E-01     3.000000E+01           
         17  traj.ascent.states:flight_path_angle_0       c    -9.999996E+00    -4.629714E-03     1.999999E+01           
         18  traj.ascent.states:flight_path_angle_1       c    -9.999996E+00    -1.092972E-02     1.999999E+01           
         19  traj.ascent.states:flight_path_angle_2       c    -9.999996E+00    -1.682805E-02     1.999999E+01           
         20  traj.ascent.states:flight_path_angle_3       c    -9.999996E+00    -3.733104E-03     1.999999E+01           
         21  traj.ascent.states:flight_path_angle_4       c    -9.999996E+00     3.556699E-02     1.999999E+01           
         22  traj.ascent.states:flight_path_angle_5       c    -9.999996E+00     4.650806E-02     1.999999E+01           
         23  traj.ascent.states:flight_path_angle_6       c    -9.999996E+00     6.076386E-02     1.999999E+01           
         24  traj.ascent.states:flight_path_angle_7       c    -9.999996E+00     7.400593E-02     1.999999E+01           
         25  traj.ascent.states:flight_path_angle_8       c    -9.999996E+00     8.149275E-02     1.999999E+01           
         26  traj.ascent.states:flight_path_angle_9       c    -9.999996E+00     1.076052E-01     1.999999E+01           
         27  traj.ascent.states:flight_path_angle_10      c    -9.999996E+00     1.243980E-01     1.999999E+01           
         28  traj.ascent.states:flight_path_angle_11      c    -9.999996E+00     1.148223E-01     1.999999E+01           
         29  traj.ascent.states:altitude_0                c     0.000000E+00     5.811335E-01     7.000000E-01           
         30  traj.ascent.states:altitude_1                c     0.000000E+00     5.770878E-01     7.000000E-01           
         31  traj.ascent.states:altitude_2                c     0.000000E+00     5.577523E-01     7.000000E-01           
         32  traj.ascent.states:altitude_3                c     0.000000E+00     5.481762E-01     7.000000E-01           
         33  traj.ascent.states:altitude_4                c     0.000000E+00     5.220809E-01     7.000000E-01           
         34  traj.ascent.states:altitude_5                c     0.000000E+00     4.968251E-01     7.000000E-01           
         35  traj.ascent.states:altitude_6                c     0.000000E+00     4.943967E-01     7.000000E-01           
         36  traj.ascent.states:altitude_7                c     0.000000E+00     4.902336E-01     7.000000E-01           
         37  traj.ascent.states:altitude_8                c     0.000000E+00     4.854764E-01     7.000000E-01           
         38  traj.ascent.states:altitude_9                c     0.000000E+00     4.832352E-01     7.000000E-01           
         39  traj.ascent.states:altitude_10               c     0.000000E+00     4.859963E-01     7.000000E-01           
         40  traj.ascent.states:altitude_11               c     0.000000E+00     4.974812E-01     7.000000E-01           
         41  traj.ascent.states:altitude_12               c     0.000000E+00     5.000000E-01     7.000000E-01           
         42  traj.ascent.states:velocity_0                c     0.000000E+00     7.739338E-01     3.500000E+00           
         43  traj.ascent.states:velocity_1                c     0.000000E+00     7.907358E-01     3.500000E+00           
         44  traj.ascent.states:velocity_2                c     0.000000E+00     7.960919E-01     3.500000E+00           
         45  traj.ascent.states:velocity_3                c     0.000000E+00     8.112733E-01     3.500000E+00           
         46  traj.ascent.states:velocity_4                c     0.000000E+00     8.275859E-01     3.500000E+00           
         47  traj.ascent.states:velocity_5                c     0.000000E+00     8.308342E-01     3.500000E+00           
         48  traj.ascent.states:velocity_6                c     0.000000E+00     8.381708E-01     3.500000E+00           
         49  traj.ascent.states:velocity_7                c     0.000000E+00     8.479422E-01     3.500000E+00           
         50  traj.ascent.states:velocity_8                c     0.000000E+00     8.510685E-01     3.500000E+00           
         51  traj.ascent.states:velocity_9                c     0.000000E+00     8.551861E-01     3.500000E+00           
         52  traj.ascent.states:velocity_10               c     0.000000E+00     8.587320E-01     3.500000E+00           
         53  traj.ascent.states:velocity_11               c     0.000000E+00     8.606683E-01     3.500000E+00           
         54  traj.ascent.states:mass_0                    c     0.000000E+00     1.157766E+00     6.666667E+15           
         55  traj.ascent.states:mass_1                    c     0.000000E+00     1.157754E+00     6.666667E+15           
         56  traj.ascent.states:mass_2                    c     0.000000E+00     1.157739E+00     6.666667E+15           
         57  traj.ascent.states:mass_3                    c     0.000000E+00     1.157734E+00     6.666667E+15           
         58  traj.ascent.states:mass_4                    c     0.000000E+00     1.157723E+00     6.666667E+15           
         59  traj.ascent.states:mass_5                    c     0.000000E+00     1.157707E+00     6.666667E+15           
         60  traj.ascent.states:mass_6                    c     0.000000E+00     1.157702E+00     6.666667E+15           
         61  traj.ascent.states:mass_7                    c     0.000000E+00     1.157691E+00     6.666667E+15           
         62  traj.ascent.states:mass_8                    c     0.000000E+00     1.157675E+00     6.666667E+15           
         63  traj.ascent.states:mass_9                    c     0.000000E+00     1.157670E+00     6.666667E+15           
         64  traj.ascent.states:mass_10                   c     0.000000E+00     1.157659E+00     6.666667E+15           
         65  traj.ascent.states:mass_11                   c     0.000000E+00     1.157643E+00     6.666667E+15           
         66  traj.ascent.states:mass_12                   c     0.000000E+00     1.157638E+00     6.666667E+15           
         67  traj.ascent.states:distance_0                c     0.000000E+00     5.742292E-01     1.500000E+00           
         68  traj.ascent.states:distance_1                c     0.000000E+00     5.919617E-01     1.500000E+00           
         69  traj.ascent.states:distance_2                c     0.000000E+00     5.977528E-01     1.500000E+00           
         70  traj.ascent.states:distance_3                c     0.000000E+00     6.109927E-01     1.500000E+00           
         71  traj.ascent.states:distance_4                c     0.000000E+00     6.294097E-01     1.500000E+00           
         72  traj.ascent.states:distance_5                c     0.000000E+00     6.351760E-01     1.500000E+00           
         73  traj.ascent.states:distance_6                c     0.000000E+00     6.480857E-01     1.500000E+00           
         74  traj.ascent.states:distance_7                c     0.000000E+00     6.653070E-01     1.500000E+00           
         75  traj.ascent.states:distance_8                c     0.000000E+00     6.706166E-01     1.500000E+00           
         76  traj.ascent.states:distance_9                c     0.000000E+00     6.823280E-01     1.500000E+00           
         77  traj.ascent.states:distance_10               c     0.000000E+00     6.973190E-01     1.500000E+00           
         78  traj.ascent.states:distance_11               c     0.000000E+00     7.017758E-01     1.500000E+00           
         79  traj.ascent.controls:alpha_0                 c    -6.000000E+00     2.019022E+00     6.000000E+00           
         80  traj.ascent.controls:alpha_1                 c    -6.000000E+00     1.207215E+00     6.000000E+00           
         81  traj.ascent.controls:alpha_2                 c    -6.000000E+00     8.923152E-01     6.000000E+00           
         82  traj.ascent.controls:alpha_3                 c    -6.000000E+00     9.871789E-01     6.000000E+00           
         83  traj.ascent.controls:alpha_4                 c    -6.000000E+00     1.216348E+00     6.000000E+00           
         84  traj.ascent.controls:alpha_5                 c    -6.000000E+00     1.220122E+00     6.000000E+00           
         85  traj.ascent.controls:alpha_6                 c    -6.000000E+00     1.145838E+00     6.000000E+00           
         86  traj.ascent.controls:alpha_7                 c    -6.000000E+00     1.005389E+00     6.000000E+00           
         87  traj.ascent.controls:alpha_8                 c    -6.000000E+00     1.045409E+00     6.000000E+00           
         88  traj.ascent.controls:alpha_9                 c    -6.000000E+00     1.114561E+00     6.000000E+00           
         89  traj.ascent.controls:alpha_10                c    -6.000000E+00     1.165205E+00     6.000000E+00           
         90  traj.ascent.controls:alpha_11                c    -6.000000E+00     7.793180E-01     6.000000E+00           
         91  traj.climb1.t_duration_0                     c     1.818182E-01     1.021397E+00     1.818182E+00           
         92  traj.climb1.states:altitude_0                c     4.000000E-02     7.167701E-02     1.100000E+00           
         93  traj.climb1.states:altitude_1                c     4.000000E-02     4.140071E-01     1.100000E+00           
         94  traj.climb1.states:altitude_2                c     4.000000E-02     8.659357E-01     1.100000E+00           
         95  traj.climb1.states:altitude_3                c     4.000000E-02     1.000000E+00     1.100000E+00           
         96  traj.climb1.states:mass_0                    c     0.000000E+00     1.156907E+00     6.666667E+15           
         97  traj.climb1.states:mass_1                    c     0.000000E+00     1.155449E+00     6.666667E+15           
         98  traj.climb1.states:mass_2                    c     0.000000E+00     1.153566E+00     6.666667E+15           
         99  traj.climb1.states:mass_3                    c     0.000000E+00     1.153000E+00     6.666667E+15           
        100  traj.climb1.states:distance_0                c     0.000000E+00     7.183512E-01     5.000000E+01           
        101  traj.climb1.states:distance_1                c     0.000000E+00     1.344616E+00     5.000000E+01           
        102  traj.climb1.states:distance_2                c     0.000000E+00     1.552378E+00     5.000000E+01           
        103  traj.climb2.t_duration_0                     c     2.325581E-02     1.299466E-01     1.976744E+00           
        104  traj.climb2.states:altitude_0                c     3.000000E-01     5.441597E-01     1.333333E+00           
        105  traj.climb2.states:altitude_1                c     3.000000E-01     7.757551E-01     1.333333E+00           
        106  traj.climb2.states:altitude_2                c     3.000000E-01     8.374224E-01     1.333333E+00           
        107  traj.climb2.states:altitude_3                c     3.000000E-01     9.594748E-01     1.333333E+00           
        108  traj.climb2.states:altitude_4                c     3.000000E-01     1.081075E+00     1.333333E+00           
        109  traj.climb2.states:altitude_5                c     3.000000E-01     1.113365E+00     1.333333E+00           
        110  traj.climb2.states:altitude_6                c     3.000000E-01     1.180970E+00     1.333333E+00           
        111  traj.climb2.states:altitude_7                c     3.000000E-01     1.239994E+00     1.333333E+00           
        112  traj.climb2.states:altitude_8                c     3.000000E-01     1.250000E+00     1.333333E+00           
        113  traj.climb2.states:mass_0                    c     0.000000E+00     1.153000E+00     6.666667E+15           
        114  traj.climb2.states:mass_1                    c     0.000000E+00     1.150352E+00     6.666667E+15           
        115  traj.climb2.states:mass_2                    c     0.000000E+00     1.147244E+00     6.666667E+15           
        116  traj.climb2.states:mass_3                    c     0.000000E+00     1.146355E+00     6.666667E+15           
        117  traj.climb2.states:mass_4                    c     0.000000E+00     1.144377E+00     6.666667E+15           
        118  traj.climb2.states:mass_5                    c     0.000000E+00     1.141936E+00     6.666667E+15           
        119  traj.climb2.states:mass_6                    c     0.000000E+00     1.141199E+00     6.666667E+15           
        120  traj.climb2.states:mass_7                    c     0.000000E+00     1.139536E+00     6.666667E+15           
        121  traj.climb2.states:mass_8                    c     0.000000E+00     1.137417E+00     6.666667E+15           
        122  traj.climb2.states:mass_9                    c     0.000000E+00     1.136796E+00     6.666667E+15           
        123  traj.climb2.states:distance_0                c     2.000000E-02     5.537061E-02     2.000000E+00           
        124  traj.climb2.states:distance_1                c     2.000000E-02     9.306663E-02     2.000000E+00           
        125  traj.climb2.states:distance_2                c     2.000000E-02     1.059079E-01     2.000000E+00           
        126  traj.climb2.states:distance_3                c     2.000000E-02     1.368422E-01     2.000000E+00           
        127  traj.climb2.states:distance_4                c     2.000000E-02     1.825136E-01     2.000000E+00           
        128  traj.climb2.states:distance_5                c     2.000000E-02     1.975584E-01     2.000000E+00           
        129  traj.climb2.states:distance_6                c     2.000000E-02     2.332554E-01     2.000000E+00           
        130  traj.climb2.states:distance_7                c     2.000000E-02     2.819304E-01     2.000000E+00           
        131  traj.climb2.states:distance_8                c     2.000000E-02     2.961779E-01     2.000000E+00           
        132  traj.cruise.t_initial_0                      c     0.000000E+00     1.705194E+00     1.000000E+02           
        133  traj.cruise.t_duration_0                     c    -2.000000E+02    -6.751920E-01    -2.000000E-05           
        134  traj.desc1.t_duration_0                      c     5.000000E-01     9.562037E-01     1.500000E+00           
        135  traj.desc1.states:altitude_0                 c     3.333333E-02     1.253135E+00     1.333333E+00           
        136  traj.desc1.states:altitude_1                 c     3.333333E-02     1.162537E+00     1.333333E+00           
        137  traj.desc1.states:altitude_2                 c     3.333333E-02     1.031382E+00     1.333333E+00           
        138  traj.desc1.states:altitude_3                 c     3.333333E-02     9.874374E-01     1.333333E+00           
        139  traj.desc1.states:altitude_4                 c     3.333333E-02     8.811771E-01     1.333333E+00           
        140  traj.desc1.states:altitude_5                 c     3.333333E-02     7.157101E-01     1.333333E+00           
        141  traj.desc1.states:altitude_6                 c     3.333333E-02     6.572832E-01     1.333333E+00           
        142  traj.desc1.states:altitude_7                 c     3.333333E-02     5.354824E-01     1.333333E+00           
        143  traj.desc1.states:altitude_8                 c     3.333333E-02     3.798934E-01     1.333333E+00           
        144  traj.desc1.states:altitude_9                 c     3.333333E-02     3.333333E-01     1.333333E+00           
        145  traj.desc1.states:mass_0                     c     0.000000E+00     9.768559E-01     7.142857E+15           
        146  traj.desc1.states:mass_1                     c     0.000000E+00     9.767944E-01     7.142857E+15           
        147  traj.desc1.states:mass_2                     c     0.000000E+00     9.766960E-01     7.142857E+15           
        148  traj.desc1.states:mass_3                     c     0.000000E+00     9.766614E-01     7.142857E+15           
        149  traj.desc1.states:mass_4                     c     0.000000E+00     9.765741E-01     7.142857E+15           
        150  traj.desc1.states:mass_5                     c     0.000000E+00     9.764323E-01     7.142857E+15           
        151  traj.desc1.states:mass_6                     c     0.000000E+00     9.763814E-01     7.142857E+15           
        152  traj.desc1.states:mass_7                     c     0.000000E+00     9.762578E-01     7.142857E+15           
        153  traj.desc1.states:mass_8                     c     0.000000E+00     9.760685E-01     7.142857E+15           
        154  traj.desc1.states:mass_9                     c     0.000000E+00     9.760038E-01     7.142857E+15           
        155  traj.desc1.states:distance_0                 c     8.163265E-01     9.735273E-01     1.360544E+00           
        156  traj.desc1.states:distance_1                 c     8.163265E-01     9.768182E-01     1.360544E+00           
        157  traj.desc1.states:distance_2                 c     8.163265E-01     9.778738E-01     1.360544E+00           
        158  traj.desc1.states:distance_3                 c     8.163265E-01     9.803113E-01     1.360544E+00           
        159  traj.desc1.states:distance_4                 c     8.163265E-01     9.837353E-01     1.360544E+00           
        160  traj.desc1.states:distance_5                 c     8.163265E-01     9.848364E-01     1.360544E+00           
        161  traj.desc1.states:distance_6                 c     8.163265E-01     9.872033E-01     1.360544E+00           
        162  traj.desc1.states:distance_7                 c     8.163265E-01     9.902551E-01     1.360544E+00           
        163  traj.desc1.states:distance_8                 c     8.163265E-01     9.911761E-01     1.360544E+00           
        164  traj.desc2.t_duration_0                      c     3.921569E-02     1.691248E-01     1.960784E+00           
        165  traj.desc2.states:altitude_0                 c    -5.555556E-02     9.201194E-01     1.111111E+00           
        166  traj.desc2.states:altitude_1                 c    -5.555556E-02     7.517188E-01     1.111111E+00           
        167  traj.desc2.states:altitude_2                 c    -5.555556E-02     5.349885E-01     1.111111E+00           
        168  traj.desc2.states:altitude_3                 c    -5.555556E-02     3.160796E-01     1.111111E+00           
        169  traj.desc2.states:altitude_4                 c    -5.555556E-02     1.365429E-01     1.111111E+00           
        170  traj.desc2.states:altitude_5                 c    -5.555556E-02     2.674544E-02     1.111111E+00           
        171  traj.desc2.states:altitude_6                 c    -5.555556E-02    -4.801256E-13     1.111111E+00           
        172  traj.desc2.states:mass_0                     c     0.000000E+00     9.109369E-01     6.666667E+15           
        173  traj.desc2.states:mass_1                     c     0.000000E+00     9.108778E-01     6.666667E+15           
        174  traj.desc2.states:mass_2                     c     0.000000E+00     9.107457E-01     6.666667E+15           
        175  traj.desc2.states:mass_3                     c     0.000000E+00     9.105599E-01     6.666667E+15           
        176  traj.desc2.states:mass_4                     c     0.000000E+00     9.103523E-01     6.666667E+15           
        177  traj.desc2.states:mass_5                     c     0.000000E+00     9.101656E-01     6.666667E+15           
        178  traj.desc2.states:mass_6                     c     0.000000E+00     9.100436E-01     6.666667E+15           
        179  traj.desc2.states:mass_7                     c     0.000000E+00     9.100129E-01     6.666667E+15           
        180  traj.desc2.states:distance_0                 c     0.000000E+00     1.041457E+00     1.428571E+00           
        181  traj.desc2.states:distance_1                 c     0.000000E+00     1.042990E+00     1.428571E+00           
        182  traj.desc2.states:distance_2                 c     0.000000E+00     1.044981E+00     1.428571E+00           
        183  traj.desc2.states:distance_3                 c     0.000000E+00     1.047016E+00     1.428571E+00           
        184  traj.desc2.states:distance_4                 c     0.000000E+00     1.048704E+00     1.428571E+00           
        185  traj.desc2.states:distance_5                 c     0.000000E+00     1.049745E+00     1.428571E+00           
        186  traj.desc2.states:distance_6                 c     0.000000E+00     1.050000E+00     1.428571E+00           
        187  traj.groundroll.t_duration_0                 c     1.980198E-02     7.178305E-01     1.980198E+00           
        188  traj.groundroll.states:velocity_0            c     0.000000E+00     3.913909E-01     6.666667E+00           
        189  traj.groundroll.states:velocity_1            c     0.000000E+00     8.499864E-01     6.666667E+00           
        190  traj.groundroll.states:velocity_2            c     0.000000E+00     9.751807E-01     6.666667E+00           
        191  traj.groundroll.states:mass_0                c     0.000000E+00     1.158717E+00     6.666667E+15           
        192  traj.groundroll.states:mass_1                c     0.000000E+00     1.158408E+00     6.666667E+15           
        193  traj.groundroll.states:mass_2                c     0.000000E+00     1.157978E+00     6.666667E+15           
        194  traj.groundroll.states:mass_3                c     0.000000E+00     1.157840E+00     6.666667E+15           
        195  traj.groundroll.states:distance_0            c     0.000000E+00     2.206578E-01     3.333333E+00           
        196  traj.groundroll.states:distance_1            c     0.000000E+00     1.171255E+00     3.333333E+00           
        197  traj.groundroll.states:distance_2            c     0.000000E+00     1.606237E+00     3.333333E+00           
        198  traj.rotation.t_duration_0                   c     1.980198E-02     6.003200E-02     1.980198E+00           
        199  traj.rotation.states:alpha_0                 c     0.000000E+00     1.251155E-02     1.000000E+00           
        200  traj.rotation.states:alpha_1                 c     0.000000E+00     2.977483E-02     1.000000E+00           
        201  traj.rotation.states:alpha_2                 c     0.000000E+00     3.523865E-02     1.000000E+00           
        202  traj.rotation.states:velocity_0              c     0.000000E+00     1.497119E+00     1.000000E+01           
        203  traj.rotation.states:velocity_1              c     0.000000E+00     1.543392E+00     1.000000E+01           
        204  traj.rotation.states:velocity_2              c     0.000000E+00     1.557679E+00     1.000000E+01           
        205  traj.rotation.states:mass_0                  c     0.000000E+00     1.157840E+00     6.666667E+15           
        206  traj.rotation.states:mass_1                  c     0.000000E+00     1.157814E+00     6.666667E+15           
        207  traj.rotation.states:mass_2                  c     0.000000E+00     1.157777E+00     6.666667E+15           
        208  traj.rotation.states:mass_3                  c     0.000000E+00     1.157766E+00     6.666667E+15           
        209  traj.rotation.states:distance_0              c     0.000000E+00     1.018393E+00     2.000000E+00           
        210  traj.rotation.states:distance_1              c     0.000000E+00     1.096194E+00     2.000000E+00           
        211  traj.rotation.states:distance_2              c     0.000000E+00     1.121373E+00     2.000000E+00           

   Constraints (i - inequality, e - equality)
      Index  Name                                                          Type          Lower           Value           Upper    Status  Lagrange Multiplier (N/A)
          0  mission:constraints:range_residual                               e   0.000000E+00    1.339004E-09    0.000000E+00              9.00000E+100
          1  h_fit.h_init_gear                                                e   1.000000E+00   -3.086317E+02    1.000000E+00         E    9.00000E+100
          2  h_fit.h_init_flaps                                               e   1.000000E+00   -1.980033E+01    1.000000E+00         E    9.00000E+100
          3  groundroll_boundary.velocity                                     e   0.000000E+00    3.244496E-06    0.000000E+00         E    9.00000E+100
          4  gtow_constraint.GTOW                                             e   0.000000E+00   -4.457469E-08    0.000000E+00              9.00000E+100
          5  mission:constraints:mass_residual                                e   0.000000E+00    8.804876E-07    0.000000E+00              9.00000E+100
          6  taxi_groundroll_mass_constraint.mass                             e   0.000000E+00   -4.656002E-10    0.000000E+00              9.00000E+100
          7  traj.linkages.groundroll:mass_final|rotation:mass_initial        e   0.000000E+00   -2.144952E-08    0.000000E+00              9.00000E+100
          8  traj.linkages.rotation:mass_final|ascent:mass_initial            e   0.000000E+00   -1.214386E-06    0.000000E+00         E    9.00000E+100
          9  traj.linkages.rotation:alpha_final|ascent:alpha_initial          e   0.000000E+00    3.162786E-07    0.000000E+00              9.00000E+100
         10  traj.linkages.ascent:mass_final|accel:mass_initial               e   0.000000E+00   -1.930515E-06    0.000000E+00         E    9.00000E+100
         11  traj.linkages.accel:mass_final|climb1:mass_initial               e   0.000000E+00   -1.932698E-06    0.000000E+00         E    9.00000E+100
         12  traj.linkages.climb1:mass_final|climb2:mass_initial              e   0.000000E+00   -1.948385E-06    0.000000E+00         E    9.00000E+100
         13  traj.linkages.climb2:mass_final|cruise:mass_initial              e   0.000000E+00   -9.277519E-08    0.000000E+00              9.00000E+100
         14  traj.linkages.cruise:mass_final|desc1:mass_initial               e   0.000000E+00   -1.152684E-07    0.000000E+00              9.00000E+100
         15  traj.linkages.desc1:mass_final|desc2:mass_initial                e   0.000000E+00   -2.591754E-06    0.000000E+00         E    9.00000E+100
         16  traj.phases.accel->final_boundary_constraint->EAS                e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         17  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00    1.464598E-03    0.000000E+00         E    9.00000E+100
         18  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00   -1.330296E-04    0.000000E+00         E    9.00000E+100
         19  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00    2.192564E-05    0.000000E+00         E    9.00000E+100
         20  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00   -8.102454E-07    0.000000E+00              9.00000E+100
         21  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00   -3.731624E-07    0.000000E+00              9.00000E+100
         22  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00    9.680293E-08    0.000000E+00              9.00000E+100
         23  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    5.887266E-03    0.000000E+00         E    9.00000E+100
         24  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    3.266219E-03    0.000000E+00         E    9.00000E+100
         25  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    1.203153E-06    0.000000E+00         E    9.00000E+100
         26  traj.phases.ascent->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         27  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    1.380484E-02    0.000000E+00         E    9.00000E+100
         28  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    4.441568E-02    0.000000E+00         E    9.00000E+100
         29  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    3.853038E-02    0.000000E+00         E    9.00000E+100
         30  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    4.923615E-02    0.000000E+00         E    9.00000E+100
         31  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.504203E-02    0.000000E+00         E    9.00000E+100
         32  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.822982E-02    0.000000E+00         E    9.00000E+100
         33  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.088487E-02    0.000000E+00         E    9.00000E+100
         34  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    5.663389E-02    0.000000E+00         E    9.00000E+100
         35  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    6.085314E-02    0.000000E+00         E    9.00000E+100
         36  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    6.841690E-02    0.000000E+00         E    9.00000E+100
         37  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.407530E-02    0.000000E+00         E    9.00000E+100
         38  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    3.129158E-02    0.000000E+00         E    9.00000E+100
         39  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00    7.941375E-09    0.000000E+00              9.00000E+100
         40  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.066456E-02    0.000000E+00         E    9.00000E+100
         41  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -2.618508E-02    0.000000E+00         E    9.00000E+100
         42  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -3.386167E-02    0.000000E+00         E    9.00000E+100
         43  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -3.403067E-02    0.000000E+00         E    9.00000E+100
         44  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.996331E-02    0.000000E+00         E    9.00000E+100
         45  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.652913E-02    0.000000E+00         E    9.00000E+100
         46  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.616438E-02    0.000000E+00         E    9.00000E+100
         47  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -2.054141E-02    0.000000E+00         E    9.00000E+100
         48  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.960978E-02    0.000000E+00         E    9.00000E+100
         49  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.111819E-02    0.000000E+00         E    9.00000E+100
         50  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.410874E-02    0.000000E+00         E    9.00000E+100
         51  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00   -3.529673E-02    0.000000E+00         E    9.00000E+100
         52  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00   -3.524609E-03    0.000000E+00         E    9.00000E+100
         53  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    7.048000E-03    0.000000E+00         E    9.00000E+100
         54  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    9.630944E-03    0.000000E+00         E    9.00000E+100
         55  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    8.878926E-03    0.000000E+00         E    9.00000E+100
         56  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.750280E-03    0.000000E+00         E    9.00000E+100
         57  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    2.755662E-03    0.000000E+00         E    9.00000E+100
         58  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    2.923903E-03    0.000000E+00         E    9.00000E+100
         59  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.858739E-03    0.000000E+00         E    9.00000E+100
         60  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.184044E-03    0.000000E+00         E    9.00000E+100
         61  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    7.213077E-07    0.000000E+00              9.00000E+100
         62  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    1.857785E-03    0.000000E+00         E    9.00000E+100
         63  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.604698E-07    0.000000E+00              9.00000E+100
         64  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    6.959095E-07    0.000000E+00              9.00000E+100
         65  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    6.342592E-07    0.000000E+00              9.00000E+100
         66  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.631446E-07    0.000000E+00              9.00000E+100
         67  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.752087E-07    0.000000E+00              9.00000E+100
         68  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    4.688257E-07    0.000000E+00              9.00000E+100
         69  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.942549E-07    0.000000E+00              9.00000E+100
         70  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    4.524239E-07    0.000000E+00              9.00000E+100
         71  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.771071E-07    0.000000E+00              9.00000E+100
         72  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    2.963147E-07    0.000000E+00              9.00000E+100
         73  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.842102E-07    0.000000E+00              9.00000E+100
         74  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.512698E-07    0.000000E+00              9.00000E+100
         75  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    2.386550E-03    0.000000E+00         E    9.00000E+100
         76  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    7.109714E-04    0.000000E+00         E    9.00000E+100
         77  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    4.040022E-04    0.000000E+00         E    9.00000E+100
         78  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    3.363939E-04    0.000000E+00         E    9.00000E+100
         79  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    3.499745E-04    0.000000E+00         E    9.00000E+100
         80  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -1.230211E-04    0.000000E+00         E    9.00000E+100
         81  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -4.322351E-04    0.000000E+00         E    9.00000E+100
         82  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -1.096509E-03    0.000000E+00         E    9.00000E+100
         83  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -1.994956E-03    0.000000E+00         E    9.00000E+100
         84  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -2.299998E-03    0.000000E+00         E    9.00000E+100
         85  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -3.334338E-03    0.000000E+00         E    9.00000E+100
         86  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -4.775712E-03    0.000000E+00         E    9.00000E+100
         87  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00   -3.828503E-09    0.000000E+00              9.00000E+100
         88  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00    3.551076E-09    0.000000E+00              9.00000E+100
         89  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00    2.607846E-09    0.000000E+00              9.00000E+100
         90  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00   -7.028954E-09    0.000000E+00              9.00000E+100
         91  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00    2.381165E-08    0.000000E+00              9.00000E+100
         92  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00    6.787690E-09    0.000000E+00              9.00000E+100
         93  traj.phases.climb1->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         94  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00    1.198847E-04    0.000000E+00         E    9.00000E+100
         95  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00    3.475093E-04    0.000000E+00         E    9.00000E+100
         96  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00   -2.463892E-04    0.000000E+00         E    9.00000E+100
         97  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -2.263555E-06    0.000000E+00         E    9.00000E+100
         98  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -2.437732E-06    0.000000E+00         E    9.00000E+100
         99  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -8.823583E-07    0.000000E+00              9.00000E+100
        100  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -5.521311E-04    0.000000E+00         E    9.00000E+100
        101  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -7.363968E-04    0.000000E+00         E    9.00000E+100
        102  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -3.938234E-04    0.000000E+00         E    9.00000E+100
        103  traj.phases.climb2->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        104  traj.phases.climb2->final_boundary_constraint->mach              e   8.000000E-01    8.000000E-01    8.000000E-01              9.00000E+100
        105  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -4.871115E-04    0.000000E+00         E    9.00000E+100
        106  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -5.387985E-03    0.000000E+00         E    9.00000E+100
        107  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -6.567229E-03    0.000000E+00         E    9.00000E+100
        108  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    3.936598E-03    0.000000E+00         E    9.00000E+100
        109  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -1.230941E-02    0.000000E+00         E    9.00000E+100
        110  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -1.322212E-02    0.000000E+00         E    9.00000E+100
        111  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    4.136420E-03    0.000000E+00         E    9.00000E+100
        112  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    7.590919E-03    0.000000E+00         E    9.00000E+100
        113  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -2.659460E-03    0.000000E+00         E    9.00000E+100
        114  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    2.887250E-08    0.000000E+00              9.00000E+100
        115  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    5.231343E-05    0.000000E+00         E    9.00000E+100
        116  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.342698E-04    0.000000E+00         E    9.00000E+100
        117  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00   -1.637094E-05    0.000000E+00         E    9.00000E+100
        118  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    8.493104E-05    0.000000E+00         E    9.00000E+100
        119  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.046990E-04    0.000000E+00         E    9.00000E+100
        120  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.982160E-05    0.000000E+00         E    9.00000E+100
        121  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00   -2.366292E-05    0.000000E+00         E    9.00000E+100
        122  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    2.083190E-05    0.000000E+00         E    9.00000E+100
        123  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    4.946557E-07    0.000000E+00              9.00000E+100
        124  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    2.211217E-04    0.000000E+00         E    9.00000E+100
        125  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.229913E-04    0.000000E+00         E    9.00000E+100
        126  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.269905E-04    0.000000E+00         E    9.00000E+100
        127  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.529940E-04    0.000000E+00         E    9.00000E+100
        128  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.775890E-04    0.000000E+00         E    9.00000E+100
        129  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    9.319001E-04    0.000000E+00         E    9.00000E+100
        130  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.276918E-03    0.000000E+00         E    9.00000E+100
        131  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00   -2.488483E-04    0.000000E+00         E    9.00000E+100
        132  traj.phases.desc1->final_boundary_constraint->altitude           e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        133  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00    2.808183E-05    0.000000E+00         E    9.00000E+100
        134  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -7.600279E-05    0.000000E+00         E    9.00000E+100
        135  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -3.269526E-04    0.000000E+00         E    9.00000E+100
        136  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -3.709200E-04    0.000000E+00         E    9.00000E+100
        137  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -5.637960E-04    0.000000E+00         E    9.00000E+100
        138  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -2.315669E-03    0.000000E+00         E    9.00000E+100
        139  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00    1.035239E-04    0.000000E+00         E    9.00000E+100
        140  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -6.066107E-05    0.000000E+00         E    9.00000E+100
        141  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -1.161643E-04    0.000000E+00         E    9.00000E+100
        142  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -1.135168E-08    0.000000E+00              9.00000E+100
        143  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -1.591237E-07    0.000000E+00              9.00000E+100
        144  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -1.956528E-07    0.000000E+00              9.00000E+100
        145  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.728901E-07    0.000000E+00              9.00000E+100
        146  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.701770E-07    0.000000E+00              9.00000E+100
        147  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.386989E-07    0.000000E+00              9.00000E+100
        148  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -2.657785E-07    0.000000E+00              9.00000E+100
        149  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00    7.638579E-08    0.000000E+00              9.00000E+100
        150  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00    1.386148E-07    0.000000E+00              9.00000E+100
        151  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    1.101909E-03    0.000000E+00         E    9.00000E+100
        152  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    2.771372E-04    0.000000E+00         E    9.00000E+100
        153  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    2.753606E-05    0.000000E+00         E    9.00000E+100
        154  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    9.004347E-05    0.000000E+00         E    9.00000E+100
        155  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    9.312668E-05    0.000000E+00         E    9.00000E+100
        156  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    4.864784E-04    0.000000E+00         E    9.00000E+100
        157  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00   -1.577391E-04    0.000000E+00         E    9.00000E+100
        158  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00   -3.902908E-05    0.000000E+00         E    9.00000E+100
        159  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    5.141498E-06    0.000000E+00         E    9.00000E+100
        160  traj.phases.desc2->final_boundary_constraint->altitude           e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        161  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    6.764003E-05    0.000000E+00         E    9.00000E+100
        162  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    1.652023E-04    0.000000E+00         E    9.00000E+100
        163  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.141564E-04    0.000000E+00         E    9.00000E+100
        164  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.382631E-04    0.000000E+00         E    9.00000E+100
        165  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.230187E-04    0.000000E+00         E    9.00000E+100
        166  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    1.843114E-04    0.000000E+00         E    9.00000E+100
        167  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    5.853690E-05    0.000000E+00         E    9.00000E+100
        168  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00    7.052591E-09    0.000000E+00              9.00000E+100
        169  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -7.706893E-08    0.000000E+00              9.00000E+100
        170  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -2.656773E-07    0.000000E+00              9.00000E+100
        171  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -1.363871E-07    0.000000E+00              9.00000E+100
        172  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -3.510428E-07    0.000000E+00              9.00000E+100
        173  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -2.024110E-07    0.000000E+00              9.00000E+100
        174  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -3.130180E-08    0.000000E+00              9.00000E+100
        175  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00    5.384401E-08    0.000000E+00              9.00000E+100
        176  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -2.639568E-05    0.000000E+00         E    9.00000E+100
        177  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -6.606589E-05    0.000000E+00         E    9.00000E+100
        178  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -8.598270E-05    0.000000E+00         E    9.00000E+100
        179  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -7.245823E-05    0.000000E+00         E    9.00000E+100
        180  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -3.830249E-05    0.000000E+00         E    9.00000E+100
        181  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -8.277211E-06    0.000000E+00         E    9.00000E+100
        182  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -6.128665E-05    0.000000E+00         E    9.00000E+100
        183  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -9.822673E-04    0.000000E+00         E    9.00000E+100
        184  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -8.990368E-04    0.000000E+00         E    9.00000E+100
        185  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -1.621972E-09    0.000000E+00              9.00000E+100
        186  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -4.273467E-08    0.000000E+00              9.00000E+100
        187  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -7.767439E-08    0.000000E+00              9.00000E+100
        188  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00   -1.325550E-09    0.000000E+00              9.00000E+100
        189  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00    8.288562E-03    0.000000E+00         E    9.00000E+100
        190  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00    6.935584E-03    0.000000E+00         E    9.00000E+100
        191  traj.phases.rotation->final_boundary_constraint->normal_force    e   0.000000E+00    0.000000E+00    0.000000E+00              9.00000E+100
        192  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -7.133575E-08    0.000000E+00              9.00000E+100
        193  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -3.269037E-07    0.000000E+00              9.00000E+100
        194  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -2.369138E-07    0.000000E+00              9.00000E+100
        195  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00   -3.855136E-04    0.000000E+00         E    9.00000E+100
        196  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00   -5.231107E-04    0.000000E+00         E    9.00000E+100
        197  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00    2.307653E-04    0.000000E+00         E    9.00000E+100
        198  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.272794E-08    0.000000E+00              9.00000E+100
        199  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.277555E-08    0.000000E+00              9.00000E+100
        200  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.694396E-08    0.000000E+00              9.00000E+100
        201  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.055516E-03    0.000000E+00         E    9.00000E+100
        202  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.404154E-03    0.000000E+00         E    9.00000E+100
        203  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.802190E-03    0.000000E+00         E    9.00000E+100
        204  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    6.817723E-01    1.100000E+00              9.00000E+100
        205  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.279978E-01    1.100000E+00              9.00000E+100
        206  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.472447E-01    1.100000E+00              9.00000E+100
        207  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.823346E-01    1.100000E+00              9.00000E+100
        208  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.823346E-01    1.100000E+00              9.00000E+100
        209  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.730250E-01    1.100000E+00              9.00000E+100
        210  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.938300E-01    1.100000E+00              9.00000E+100
        211  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.720765E-01    1.100000E+00              9.00000E+100
        212  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.720766E-01    1.100000E+00              9.00000E+100
        213  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.310326E-01    1.100000E+00              9.00000E+100
        214  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.559812E-01    1.100000E+00              9.00000E+100
        215  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.848612E-01    1.100000E+00              9.00000E+100
        216  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.848612E-01    1.100000E+00              9.00000E+100
        217  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    5.093960E-01    1.100000E+00              9.00000E+100
        218  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.704541E-01    1.100000E+00              9.00000E+100
        219  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    2.847477E-01    1.100000E+00              9.00000E+100
        220  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.595111E+00    1.500000E+01              9.00000E+100
        221  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    4.270811E+00    1.500000E+01              9.00000E+100
        222  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.335349E+00    1.500000E+01              9.00000E+100
        223  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.471718E+00    1.500000E+01              9.00000E+100
        224  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.471718E+00    1.500000E+01              9.00000E+100
        225  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    4.367851E+00    1.500000E+01              9.00000E+100
        226  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.638448E+00    1.500000E+01              9.00000E+100
        227  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.893906E+00    1.500000E+01              9.00000E+100
        228  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.893907E+00    1.500000E+01              9.00000E+100
        229  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.008459E+00    1.500000E+01              9.00000E+100
        230  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.967276E+00    1.500000E+01              9.00000E+100
        231  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.741995E+00    1.500000E+01              9.00000E+100
        232  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.741995E+00    1.500000E+01              9.00000E+100
        233  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    1.049135E+01    1.500000E+01              9.00000E+100
        234  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    9.524072E+00    1.500000E+01              9.00000E+100
        235  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.814242E+00    1.500000E+01              9.00000E+100
        236  traj.phases.climb2->final_boundary_constraint->altitude_rate     i   1.000000E-01    3.496109E+02    1.000000E+30              9.00000E+100

--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.

The argument models/test_aircraft/aircraft_for_bench_GwGm.csv shows where the aircraft csv file lives. This argument must be one of the following three options:

  • an absolute path,

  • a relative path relative to the current working directory,

  • a relative path relative to the Aviary package.

Aviary searches for the given dataset in this order. If a dataset file exists in multiple directories, the first one Aviary finds will be used.

Note

When you run Aviary all of the outputs will be saved in the location where you run the command or run script. Outputs files such as optimization history and Dymos output files are saved at that location. A folder named ‘report’ will also be created there, and contains all of the associated report files for your runs. These reports include optimization and trajectory reports generated by OpenMDAO and Dymos.

Warning messages#

During your first run you may notice some warning messages in the Aviary output. Frequently seen warnings are:

  • PromotionWarning: Issued when there is ambiguity due to variable promotion (an OpenMDAO warning).

  • RuntimeWarning: Issued for warnings about dubious runtime features (a Python warning).

  • UserWarning: Issued for warnings about potential OpenMDAO, Dymos, and/or Aviary problems.

  • DerivativesWarning: Issued when the approximated partials or coloring cannot be evaluated as expected (an OpenMDAO warning).

Some of these warnings are expected and can be ignored. For example, the following warnings are expected when running the aircraft_for_bench_GwGm.csv model:

~/workspace/OpenMDAO/openmdao/core/group.py:326: PromotionWarning:<class Group>: Setting input defaults for input 'mission:takeoff:ascent_duration' which override previously set defaults for ['auto', 'prom', 'src_shape', 'val'].
~/workspace/dymos/dymos/phase/phase.py:2007: RuntimeWarning: Invalid options for non-optimal control 'throttle' in phase 'groundroll': lower, upper
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
~/workspace/dymos/dymos/phase/phase.py:1967: UserWarning: Phase time options have no effect because fix_initial=True or input_initial=True for phase 'traj.phases.rotation': initial_ref
  warnings.warn(f'Phase time options have no effect because fix_initial=True '

For now, we can ignore the warning messages and continue.

Level 1 run options#

In addition to a model input file, this option has additional options that are important, as seen here:

!aviary run_mission
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
usage: aviary run_mission [-h] [-o OUTDIR]
                          [--optimizer {SNOPT,IPOPT,SLSQP,None}]
                          [--phase_info PHASE_INFO] [--max_iter MAX_ITER]
                          [--shooting] [--verbosity {0,1,2,3}]
                          indeck
aviary run_mission: error: the following arguments are required: indeck

Let us discuss these in more detail:

  • -o ––outdir: Use specified directory to write output. The default is the current directory.

  • --optimizer: Name of optimizer. Choices are: SNOPT, IPOPT, SLSQP, and None. The default is SNOPT. If optimizer is None, it will be set to IPOPT or SNOPT depending on the analysis scheme. The optimization objective is fuel burn for level 1 runs. The objective is

    • mission:objectives:fuel if mission_method is GASP

    • fuel_burned if mission_method is FLOPS.

  • --phase_info: Path to phase info file. If not provided, it is default_phase_info/gasp.py if Mission origin is 2DOF and default_phase_info/flops.py for simple.

  • --max_iter: Maximum number of iterations. Default is 50.

  • shooting: Use shooting instead of collocation.

For the aircraft_for_bench_GwGm example so far, we have used the IPOPT optimizer because it is publicly available through pyOptSparse.

Aviary run_mission command#

The Level 1 CLI (i.e. run_mission option) is designed to offer the lowest barrier to entry for new users. Analysts who have experience with legacy tools, such as FLOPS and GASP, should find the switch from FORTRAN namelists to csv-based input decks to be straightforward. Aviary input decks allow the user the ability to set aircraft characteristics and basic mission parameters, such as cruise Mach number and altitude, in a simple text-based format that does not require any familiarity with Python or OpenMDAO.

Aviary can then be executed by calling aviary run_mission with the csv input deck. Although the order of the variables in the input deck is not important, you might find it helpful to separate the variables based on if they’re used as initial guesses or in the metadata.

First aircraft model#

We have a few sample aircraft csv files in aviary/models. They are all .csv files. For example, aircraft_for_bench_GwGm.csv (in test_models) is a large single aisle aircraft mission model.

Open aircraft_for_bench_GwGm.csv using your favorite editor (e.g., MS Excel or VS Code). If you choose to use MS Excel, you need to take extra steps to make sure that it is in the original format. We recommend running dos2unix first (e.g. dos2unix aircraft_for_bench_GwGm.csv).

The file aircraft_for_bench_GwGm.csv is a data input file modeling large single aisle aircraft which was converted from a GASP input deck. Here is a snippet of the file:

aircraft:air_conditioning:mass_coefficient,1.65,unitless
aircraft:anti_icing:mass,551,lbm
aircraft:apu:mass,928,lbm
......
mission:design:cruise_altitude,37500,ft
mission:design:gross_mass,175400,lbm
......
reserves,0
rotation_mass,0.99
time_to_climb,0
......
INGASP.ALR,1.11,
INGASP.BENGOB,0.05,
INGASP.CINP,0.11,
......

The input .csv file defines the aircraft and mission. In an Aviary input file you see a list of Aviary variables, their values, and units. Aviary variables are colon delimited strings (e.g. aircraft:wing:aspect_ratio). An Aviary variable name usually has three words. The first word is aircraft or mission, the second word is a subsystem name and the third is the variable name. Each variable has a value and units. Aviary requires units for every variable because it reduces errors and is a good engineering practice.

Note

If you have used the fortran_to_aviary tool to create a .csv file you must check the outputted variable values and units to ensure the file is valid and what you expect.

Be aware that some variables do not follow the standard naming format (e.g. reserves). These are used for the initial guessing of the trajectory. They are intentionally separate to prevent conflicts with some similarly named variables. They are only used once to set up the problem and then are discarded.

Finally we see a list of GASP variables that are not converted. Not all variables are converted to Aviary right now. They may represent some features in FLOPS and GASP that we haven’t implemented in Aviary yet.

To find information about a variable (e.g. description, data type, etc.), users should read the Variable Metadata Doc.

There are special rules for the mapping from the input file variable names to the metadata. For example, variable aircraft:wing:aspect_ratio in aircraft_for_bench_GwGm.csv is mapped to Aircraft.Wing.ASPECT_RATIO in aviary/variable_info/variable_meta_data.py. So, the first part (e.g., aircraft or mission) is mapped to the same word but with the first letter capitalized (e.g., Aircraft or Mission). The third word is all caps (e.g., ASPECT_RATIO). The middle part (e.g., wing) is a little more complicated. In most cases, this part capitalizes the first letter (e.g., Wing). The following words have special mappings:

  • air_conditioning -> AirConditioning

  • anti_icing -> AntiIcing

  • blended_wing_body -> BWB

  • crew_and_payload -> CrewPayload

  • horizontal_tail -> HorizontalTail

  • landing_gear -> LandingGear

  • tail_boom -> TailBoom

  • vertical_tail -> VerticalTail

This can be summarized as to capitalize the leading letter and to capitalize the first letter after special character “_”.

File aviary/variable_info/variables.py is a variable hierarchy that is for a single mission. Each mission gets a copy of this hierarchy. Below is a snippet of this file:

'''
This is a variable hierarchy that is for a single mission. Each mission
gets a copy of this hierarchy.
'''


class Aircraft:
    """Aircraft data hierarchy"""

    class AirConditioning:
        MASS = 'aircraft:air_conditioning:mass'
        MASS_COEFFICIENT = 'aircraft:air_conditioning:mass_coefficient'
        MASS_SCALER = 'aircraft:air_conditioning:mass_scaler'

    class AntiIcing:
        MASS = 'aircraft:anti_icing:mass'
        MASS_SCALER = 'aircraft:anti_icing:mass_scaler'

Aviary variables are always set to default values before the input file is read in. Then Aviary will update the values based on the user-provided .csv input file. If you want to set different values, you can set them in the .csv input file.

Other available input files#

Other input files reside in aviary/models. They include:

aviary/models/validation_data/N3CC/N3CC_generic_low_speed_polars_FLOPSinp.csv
aviary/models/large_single_aisle_1/large_single_aisle_1_GwGm.csv
aviary/models/small_single_aisle/small_single_aisle_GwGm.csv
aviary/models/test_aircraft/converter_configuration_test_data_GwGm.csv

For example, to run large_single_aisle_1 model, you execute the following command:

!aviary run_mission models/large_single_aisle_1/large_single_aisle_1_GwGm.csv --max_iter 1 --optimizer IPOPT
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
Loaded default phase_info for 2dof equations of motion
User has specified Design.NUM_PASSENGERS but CrewPayload.NUM_PASSENGERS has been left blank or set to zero.
Assuming they are equal to maintain backwards compatibility with GASP and FLOPS output files.
If you intended to have no passengers on this flight, please set Aircraft.CrewPayload.TOTAL_PAYLOAD_MASS to zero in aviary_values.
/home/runner/work/Aviary/Aviary/aviary/utils/preprocessors.py:375: UserWarning: Mount location for engines of type <turbofan_23k_1> not specified. Wing-mounted engines are assumed.
  warnings.warn(
The following variables have been overridden:
  'aircraft:fuselage:wetted_area  4000  ft**2
  'aircraft:horizontal_tail:form_factor  1.25  unitless
  'aircraft:nacelle:form_factor  [1.5]  unitless
  'aircraft:strut:fuselage_interference_factor  0  unitless
  'aircraft:vertical_tail:form_factor  1.25  unitless
  'aircraft:wing:form_factor  1.25  unitless
  'aircraft:wing:slat_span_ratio  0.9  unitless
  'aircraft:wing:slat_span_ratio  0.9  unitless
  'aircraft:wing:slat_span_ratio  0.9  unitless
--- Constraint Report [traj] ---
    --- groundroll ---
        None
    --- rotation ---
        [final]   0.0000e+00 == normal_force [lbf]
    --- ascent ---
        [final]   5.0000e+02 == altitude [ft]
        [path]    0.0000e+00 <= load_factor <= 1.1000e+00  [unitless]
        [path]    0.0000e+00 <= fuselage_pitch <= 1.5000e+01  [deg]
    --- accel ---
        [final]   2.5000e+02 == EAS [kn]
    --- climb1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- climb2 ---
        [final]   3.7500e+04 == altitude [ft]
        [final]   1.0000e-01 <= altitude_rate  [ft/min]
        [final]   8.0000e-01 == mach [unitless]
    --- cruise ---
        None
    --- desc1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- desc2 ---
        [final]   1.0000e+03 == altitude [ft]
'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.ascent' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.


'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.climb1' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.


'rhs_checking' is disabled for 'DirectSolver in 'traj.phases.desc1' <class Phase>' but that solver has redundant adjoint solves. If it is expensive to compute derivatives for this solver, turning on 'rhs_checking' may improve performance.
Full total jacobian for problem 'large_single_aisle_1_GwGm' was computed 3 times, taking 1.8603167610000355 seconds.
Total jacobian shape: (238, 212) 
Jacobian shape: (238, 212)  (8.73% nonzero)
FWD solves: 51   REV solves: 0
Total colors vs. total size: 51 vs 212  (75.94% improvement)

Sparsity computed using tolerance: 1e-25
Time to compute sparsity:   1.8603 sec
Time to compute coloring:   0.3301 sec
Memory to compute coloring:   1.7500 MB
Coloring created on: 2024-12-17 21:05:47
This is Ipopt version 3.14.17, running with linear solver MUMPS 5.7.3.

Number of nonzeros in equality constraint Jacobian...:     7490
Number of nonzeros in inequality constraint Jacobian.:     1592
Number of nonzeros in Lagrangian Hessian.............:        0
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/total_jac.py:1728: DerivativesWarning:Constraints or objectives [('traj.phases.climb2->final_boundary_constraint->mach', inds=[(11, 0)])] cannot be impacted by the design variables of the problem.
Total number of variables............................:      212
                     variables with only lower bounds:        2
                variables with lower and upper bounds:      210
                     variables with only upper bounds:        0
Total number of equality constraints.................:      204
Total number of inequality constraints...............:       33
        inequality constraints with only lower bounds:        1
   inequality constraints with lower and upper bounds:       32
        inequality constraints with only upper bounds:        0

iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls
   0  5.2731333e+00 1.74e+03 1.65e+01  -5.0 0.00e+00    -  0.00e+00 0.00e+00   0
Number of Iterations....: 1

                                   (scaled)                 (unscaled)
Objective...............:   4.4320662856342086e+00    4.4320662856342086e+00
Dual infeasibility......:   1.6269237200068577e+01    1.6269237200068577e+01
Constraint violation....:   5.3952849193409811e+01    3.1262254997489958e+02
Variable bound violation:   0.0000000000000000e+00    0.0000000000000000e+00
Complementarity.........:   1.0000016687732632e+05    1.0000016687732632e+05
Overall NLP error.......:   1.0000016687732632e+05    1.0000016687732632e+05


Number of objective function evaluations             = 2
Number of objective gradient evaluations             = 2
Number of equality constraint evaluations            = 2
Number of inequality constraint evaluations          = 2
Number of equality constraint Jacobian evaluations   = 2
Number of inequality constraint Jacobian evaluations = 2
Number of Lagrangian Hessian evaluations             = 0
Total seconds in IPOPT                               = 2.396

EXIT: Maximum Number of Iterations Exceeded.
minimal_print is not available for this solution


Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
    Objective Function: _objfunc

    Solution: 
--------------------------------------------------------------------------------
    Total Time:                    2.3970
       User Objective Time :       0.6921
       User Sensitivity Time :     1.6692
       Interface Time :            0.0224
       Opt Solver Time:            0.0134
    Calls to Objective Function :       3
    Calls to Sens Function :            3


   Objectives
      Index  Name                               Value
          0  mission:objectives:fuel     4.432066E+00

   Variables (c - continuous, i - integer, d - discrete)
      Index  Name                                      Type      Lower Bound            Value      Upper Bound     Status
          0  mission:design:gross_mass_0                  c     5.714286E-05     9.949129E-01     1.000000E+30           
          1  mission:summary:gross_mass_0                 c     5.714286E-05     9.949129E-01     1.000000E+30           
          2  mission:takeoff:ascent_t_initial_0           c     0.000000E+00     5.640273E-01     3.333333E+00           
          3  mission:takeoff:ascent_duration_0            c     1.000000E-01     5.382008E-01     1.000000E+02           
          4  tau_gear_0                                   c     1.000000E-02     2.004342E-01     1.000000E+00           
          5  tau_flaps_0                                  c     1.000000E-02     9.684392E-01     1.000000E+00           
          6  traj.accel.t_duration_0                      c     9.950249E-03     2.936943E-01     1.990050E+00           
          7  traj.accel.states:velocity_0                 c     0.000000E+00     5.087915E-01     1.200000E+00           
          8  traj.accel.states:velocity_1                 c     0.000000E+00     8.990488E-01     1.200000E+00           
          9  traj.accel.states:velocity_2                 c     0.000000E+00     1.018243E+00     1.200000E+00           
         10  traj.accel.states:mass_0                     c     0.000000E+00     1.158263E+00     6.666667E+15           
         11  traj.accel.states:mass_1                     c     0.000000E+00     1.158006E+00     6.666667E+15           
         12  traj.accel.states:mass_2                     c     0.000000E+00     1.157646E+00     6.666667E+15           
         13  traj.accel.states:mass_3                     c     0.000000E+00     1.157531E+00     6.666667E+15           
         14  traj.accel.states:distance_0                 c     0.000000E+00     3.426106E-01     3.000000E+01           
         15  traj.accel.states:distance_1                 c     0.000000E+00     5.213065E-01     3.000000E+01           
         16  traj.accel.states:distance_2                 c     0.000000E+00     5.837040E-01     3.000000E+01           
         17  traj.ascent.states:flight_path_angle_0       c    -9.999996E+00    -4.662254E-03     1.999999E+01           
         18  traj.ascent.states:flight_path_angle_1       c    -9.999996E+00    -1.101022E-02     1.999999E+01           
         19  traj.ascent.states:flight_path_angle_2       c    -9.999996E+00    -1.692525E-02     1.999999E+01           
         20  traj.ascent.states:flight_path_angle_3       c    -9.999996E+00    -3.826068E-03     1.999999E+01           
         21  traj.ascent.states:flight_path_angle_4       c    -9.999996E+00     3.551656E-02     1.999999E+01           
         22  traj.ascent.states:flight_path_angle_5       c    -9.999996E+00     4.646857E-02     1.999999E+01           
         23  traj.ascent.states:flight_path_angle_6       c    -9.999996E+00     6.073000E-02     1.999999E+01           
         24  traj.ascent.states:flight_path_angle_7       c    -9.999996E+00     7.396081E-02     1.999999E+01           
         25  traj.ascent.states:flight_path_angle_8       c    -9.999996E+00     8.144627E-02     1.999999E+01           
         26  traj.ascent.states:flight_path_angle_9       c    -9.999996E+00     1.075661E-01     1.999999E+01           
         27  traj.ascent.states:flight_path_angle_10      c    -9.999996E+00     1.243560E-01     1.999999E+01           
         28  traj.ascent.states:flight_path_angle_11      c    -9.999996E+00     1.147667E-01     1.999999E+01           
         29  traj.ascent.states:altitude_0                c     0.000000E+00     5.818720E-01     7.000000E-01           
         30  traj.ascent.states:altitude_1                c     0.000000E+00     5.778146E-01     7.000000E-01           
         31  traj.ascent.states:altitude_2                c     0.000000E+00     5.584218E-01     7.000000E-01           
         32  traj.ascent.states:altitude_3                c     0.000000E+00     5.488168E-01     7.000000E-01           
         33  traj.ascent.states:altitude_4                c     0.000000E+00     5.226454E-01     7.000000E-01           
         34  traj.ascent.states:altitude_5                c     0.000000E+00     4.972917E-01     7.000000E-01           
         35  traj.ascent.states:altitude_6                c     0.000000E+00     4.948390E-01     7.000000E-01           
         36  traj.ascent.states:altitude_7                c     0.000000E+00     4.906178E-01     7.000000E-01           
         37  traj.ascent.states:altitude_8                c     0.000000E+00     4.857664E-01     7.000000E-01           
         38  traj.ascent.states:altitude_9                c     0.000000E+00     4.834891E-01     7.000000E-01           
         39  traj.ascent.states:altitude_10               c     0.000000E+00     4.861666E-01     7.000000E-01           
         40  traj.ascent.states:altitude_11               c     0.000000E+00     4.975256E-01     7.000000E-01           
         41  traj.ascent.states:altitude_12               c     0.000000E+00     5.000000E-01     7.000000E-01           
         42  traj.ascent.states:velocity_0                c     0.000000E+00     7.735769E-01     3.500000E+00           
         43  traj.ascent.states:velocity_1                c     0.000000E+00     7.903587E-01     3.500000E+00           
         44  traj.ascent.states:velocity_2                c     0.000000E+00     7.957103E-01     3.500000E+00           
         45  traj.ascent.states:velocity_3                c     0.000000E+00     8.108858E-01     3.500000E+00           
         46  traj.ascent.states:velocity_4                c     0.000000E+00     8.271845E-01     3.500000E+00           
         47  traj.ascent.states:velocity_5                c     0.000000E+00     8.304254E-01     3.500000E+00           
         48  traj.ascent.states:velocity_6                c     0.000000E+00     8.377463E-01     3.500000E+00           
         49  traj.ascent.states:velocity_7                c     0.000000E+00     8.474992E-01     3.500000E+00           
         50  traj.ascent.states:velocity_8                c     0.000000E+00     8.506211E-01     3.500000E+00           
         51  traj.ascent.states:velocity_9                c     0.000000E+00     8.547275E-01     3.500000E+00           
         52  traj.ascent.states:velocity_10               c     0.000000E+00     8.582598E-01     3.500000E+00           
         53  traj.ascent.states:velocity_11               c     0.000000E+00     8.601932E-01     3.500000E+00           
         54  traj.ascent.states:mass_0                    c     0.000000E+00     1.158391E+00     6.666667E+15           
         55  traj.ascent.states:mass_1                    c     0.000000E+00     1.158380E+00     6.666667E+15           
         56  traj.ascent.states:mass_2                    c     0.000000E+00     1.158364E+00     6.666667E+15           
         57  traj.ascent.states:mass_3                    c     0.000000E+00     1.158359E+00     6.666667E+15           
         58  traj.ascent.states:mass_4                    c     0.000000E+00     1.158348E+00     6.666667E+15           
         59  traj.ascent.states:mass_5                    c     0.000000E+00     1.158333E+00     6.666667E+15           
         60  traj.ascent.states:mass_6                    c     0.000000E+00     1.158328E+00     6.666667E+15           
         61  traj.ascent.states:mass_7                    c     0.000000E+00     1.158316E+00     6.666667E+15           
         62  traj.ascent.states:mass_8                    c     0.000000E+00     1.158301E+00     6.666667E+15           
         63  traj.ascent.states:mass_9                    c     0.000000E+00     1.158296E+00     6.666667E+15           
         64  traj.ascent.states:mass_10                   c     0.000000E+00     1.158284E+00     6.666667E+15           
         65  traj.ascent.states:mass_11                   c     0.000000E+00     1.158268E+00     6.666667E+15           
         66  traj.ascent.states:mass_12                   c     0.000000E+00     1.158263E+00     6.666667E+15           
         67  traj.ascent.states:distance_0                c     0.000000E+00     5.739882E-01     1.500000E+00           
         68  traj.ascent.states:distance_1                c     0.000000E+00     5.916628E-01     1.500000E+00           
         69  traj.ascent.states:distance_2                c     0.000000E+00     5.974352E-01     1.500000E+00           
         70  traj.ascent.states:distance_3                c     0.000000E+00     6.106318E-01     1.500000E+00           
         71  traj.ascent.states:distance_4                c     0.000000E+00     6.289877E-01     1.500000E+00           
         72  traj.ascent.states:distance_5                c     0.000000E+00     6.347341E-01     1.500000E+00           
         73  traj.ascent.states:distance_6                c     0.000000E+00     6.475972E-01     1.500000E+00           
         74  traj.ascent.states:distance_7                c     0.000000E+00     6.647522E-01     1.500000E+00           
         75  traj.ascent.states:distance_8                c     0.000000E+00     6.700405E-01     1.500000E+00           
         76  traj.ascent.states:distance_9                c     0.000000E+00     6.817021E-01     1.500000E+00           
         77  traj.ascent.states:distance_10               c     0.000000E+00     6.966229E-01     1.500000E+00           
         78  traj.ascent.states:distance_11               c     0.000000E+00     7.010570E-01     1.500000E+00           
         79  traj.ascent.controls:alpha_0                 c    -6.000000E+00     2.021233E+00     6.000000E+00           
         80  traj.ascent.controls:alpha_1                 c    -6.000000E+00     1.208358E+00     6.000000E+00           
         81  traj.ascent.controls:alpha_2                 c    -6.000000E+00     8.931579E-01     6.000000E+00           
         82  traj.ascent.controls:alpha_3                 c    -6.000000E+00     9.882097E-01     6.000000E+00           
         83  traj.ascent.controls:alpha_4                 c    -6.000000E+00     1.217799E+00     6.000000E+00           
         84  traj.ascent.controls:alpha_5                 c    -6.000000E+00     1.221602E+00     6.000000E+00           
         85  traj.ascent.controls:alpha_6                 c    -6.000000E+00     1.147195E+00     6.000000E+00           
         86  traj.ascent.controls:alpha_7                 c    -6.000000E+00     1.006488E+00     6.000000E+00           
         87  traj.ascent.controls:alpha_8                 c    -6.000000E+00     1.046473E+00     6.000000E+00           
         88  traj.ascent.controls:alpha_9                 c    -6.000000E+00     1.115690E+00     6.000000E+00           
         89  traj.ascent.controls:alpha_10                c    -6.000000E+00     1.166394E+00     6.000000E+00           
         90  traj.ascent.controls:alpha_11                c    -6.000000E+00     7.801499E-01     6.000000E+00           
         91  traj.climb1.t_duration_0                     c     1.818182E-01     1.021879E+00     1.818182E+00           
         92  traj.climb1.states:altitude_0                c     4.000000E-02     7.169730E-02     1.100000E+00           
         93  traj.climb1.states:altitude_1                c     4.000000E-02     4.140227E-01     1.100000E+00           
         94  traj.climb1.states:altitude_2                c     4.000000E-02     8.659405E-01     1.100000E+00           
         95  traj.climb1.states:altitude_3                c     4.000000E-02     1.000000E+00     1.100000E+00           
         96  traj.climb1.states:mass_0                    c     0.000000E+00     1.157531E+00     6.666667E+15           
         97  traj.climb1.states:mass_1                    c     0.000000E+00     1.156072E+00     6.666667E+15           
         98  traj.climb1.states:mass_2                    c     0.000000E+00     1.154189E+00     6.666667E+15           
         99  traj.climb1.states:mass_3                    c     0.000000E+00     1.153623E+00     6.666667E+15           
        100  traj.climb1.states:distance_0                c     0.000000E+00     7.186850E-01     5.000000E+01           
        101  traj.climb1.states:distance_1                c     0.000000E+00     1.345250E+00     5.000000E+01           
        102  traj.climb1.states:distance_2                c     0.000000E+00     1.553110E+00     5.000000E+01           
        103  traj.climb2.t_duration_0                     c     2.325581E-02     1.301056E-01     1.976744E+00           
        104  traj.climb2.states:altitude_0                c     3.000000E-01     5.442344E-01     1.333333E+00           
        105  traj.climb2.states:altitude_1                c     3.000000E-01     7.758903E-01     1.333333E+00           
        106  traj.climb2.states:altitude_2                c     3.000000E-01     8.375682E-01     1.333333E+00           
        107  traj.climb2.states:altitude_3                c     3.000000E-01     9.596309E-01     1.333333E+00           
        108  traj.climb2.states:altitude_4                c     3.000000E-01     1.081211E+00     1.333333E+00           
        109  traj.climb2.states:altitude_5                c     3.000000E-01     1.113490E+00     1.333333E+00           
        110  traj.climb2.states:altitude_6                c     3.000000E-01     1.181064E+00     1.333333E+00           
        111  traj.climb2.states:altitude_7                c     3.000000E-01     1.240021E+00     1.333333E+00           
        112  traj.climb2.states:altitude_8                c     3.000000E-01     1.250000E+00     1.333333E+00           
        113  traj.climb2.states:mass_0                    c     0.000000E+00     1.153623E+00     6.666667E+15           
        114  traj.climb2.states:mass_1                    c     0.000000E+00     1.150971E+00     6.666667E+15           
        115  traj.climb2.states:mass_2                    c     0.000000E+00     1.147859E+00     6.666667E+15           
        116  traj.climb2.states:mass_3                    c     0.000000E+00     1.146969E+00     6.666667E+15           
        117  traj.climb2.states:mass_4                    c     0.000000E+00     1.144988E+00     6.666667E+15           
        118  traj.climb2.states:mass_5                    c     0.000000E+00     1.142544E+00     6.666667E+15           
        119  traj.climb2.states:mass_6                    c     0.000000E+00     1.141806E+00     6.666667E+15           
        120  traj.climb2.states:mass_7                    c     0.000000E+00     1.140141E+00     6.666667E+15           
        121  traj.climb2.states:mass_8                    c     0.000000E+00     1.138019E+00     6.666667E+15           
        122  traj.climb2.states:mass_9                    c     0.000000E+00     1.137398E+00     6.666667E+15           
        123  traj.climb2.states:distance_0                c     2.000000E-02     5.541482E-02     2.000000E+00           
        124  traj.climb2.states:distance_1                c     2.000000E-02     9.315566E-02     2.000000E+00           
        125  traj.climb2.states:distance_2                c     2.000000E-02     1.060120E-01     2.000000E+00           
        126  traj.climb2.states:distance_3                c     2.000000E-02     1.369826E-01     2.000000E+00           
        127  traj.climb2.states:distance_4                c     2.000000E-02     1.827075E-01     2.000000E+00           
        128  traj.climb2.states:distance_5                c     2.000000E-02     1.977700E-01     2.000000E+00           
        129  traj.climb2.states:distance_6                c     2.000000E-02     2.335091E-01     2.000000E+00           
        130  traj.climb2.states:distance_7                c     2.000000E-02     2.822427E-01     2.000000E+00           
        131  traj.climb2.states:distance_8                c     2.000000E-02     2.965078E-01     2.000000E+00           
        132  traj.cruise.t_initial_0                      c     0.000000E+00     1.706097E+00     1.000000E+02           
        133  traj.cruise.t_duration_0                     c    -2.000000E+02    -6.754130E-01    -2.000000E-05           
        134  traj.desc1.t_duration_0                      c     5.000000E-01     9.561960E-01     1.500000E+00           
        135  traj.desc1.states:altitude_0                 c     3.333333E-02     1.252952E+00     1.333333E+00           
        136  traj.desc1.states:altitude_1                 c     3.333333E-02     1.162367E+00     1.333333E+00           
        137  traj.desc1.states:altitude_2                 c     3.333333E-02     1.031227E+00     1.333333E+00           
        138  traj.desc1.states:altitude_3                 c     3.333333E-02     9.872877E-01     1.333333E+00           
        139  traj.desc1.states:altitude_4                 c     3.333333E-02     8.810401E-01     1.333333E+00           
        140  traj.desc1.states:altitude_5                 c     3.333333E-02     7.156067E-01     1.333333E+00           
        141  traj.desc1.states:altitude_6                 c     3.333333E-02     6.571990E-01     1.333333E+00           
        142  traj.desc1.states:altitude_7                 c     3.333333E-02     5.354315E-01     1.333333E+00           
        143  traj.desc1.states:altitude_8                 c     3.333333E-02     3.798822E-01     1.333333E+00           
        144  traj.desc1.states:altitude_9                 c     3.333333E-02     3.333333E-01     1.333333E+00           
        145  traj.desc1.states:mass_0                     c     0.000000E+00     9.774216E-01     7.142857E+15           
        146  traj.desc1.states:mass_1                     c     0.000000E+00     9.773601E-01     7.142857E+15           
        147  traj.desc1.states:mass_2                     c     0.000000E+00     9.772617E-01     7.142857E+15           
        148  traj.desc1.states:mass_3                     c     0.000000E+00     9.772271E-01     7.142857E+15           
        149  traj.desc1.states:mass_4                     c     0.000000E+00     9.771397E-01     7.142857E+15           
        150  traj.desc1.states:mass_5                     c     0.000000E+00     9.769980E-01     7.142857E+15           
        151  traj.desc1.states:mass_6                     c     0.000000E+00     9.769470E-01     7.142857E+15           
        152  traj.desc1.states:mass_7                     c     0.000000E+00     9.768235E-01     7.142857E+15           
        153  traj.desc1.states:mass_8                     c     0.000000E+00     9.766342E-01     7.142857E+15           
        154  traj.desc1.states:mass_9                     c     0.000000E+00     9.765694E-01     7.142857E+15           
        155  traj.desc1.states:distance_0                 c     8.163265E-01     9.735258E-01     1.360544E+00           
        156  traj.desc1.states:distance_1                 c     8.163265E-01     9.768168E-01     1.360544E+00           
        157  traj.desc1.states:distance_2                 c     8.163265E-01     9.778724E-01     1.360544E+00           
        158  traj.desc1.states:distance_3                 c     8.163265E-01     9.803099E-01     1.360544E+00           
        159  traj.desc1.states:distance_4                 c     8.163265E-01     9.837339E-01     1.360544E+00           
        160  traj.desc1.states:distance_5                 c     8.163265E-01     9.848349E-01     1.360544E+00           
        161  traj.desc1.states:distance_6                 c     8.163265E-01     9.872017E-01     1.360544E+00           
        162  traj.desc1.states:distance_7                 c     8.163265E-01     9.902534E-01     1.360544E+00           
        163  traj.desc1.states:distance_8                 c     8.163265E-01     9.911744E-01     1.360544E+00           
        164  traj.desc2.t_duration_0                      c     3.921569E-02     1.691568E-01     1.960784E+00           
        165  traj.desc2.states:altitude_0                 c    -5.555556E-02     9.201203E-01     1.111111E+00           
        166  traj.desc2.states:altitude_1                 c    -5.555556E-02     7.517211E-01     1.111111E+00           
        167  traj.desc2.states:altitude_2                 c    -5.555556E-02     5.349917E-01     1.111111E+00           
        168  traj.desc2.states:altitude_3                 c    -5.555556E-02     3.160824E-01     1.111111E+00           
        169  traj.desc2.states:altitude_4                 c    -5.555556E-02     1.365445E-01     1.111111E+00           
        170  traj.desc2.states:altitude_5                 c    -5.555556E-02     2.674580E-02     1.111111E+00           
        171  traj.desc2.states:altitude_6                 c    -5.555556E-02    -4.733064E-13     1.111111E+00           
        172  traj.desc2.states:mass_0                     c     0.000000E+00     9.114648E-01     6.666667E+15           
        173  traj.desc2.states:mass_1                     c     0.000000E+00     9.114057E-01     6.666667E+15           
        174  traj.desc2.states:mass_2                     c     0.000000E+00     9.112736E-01     6.666667E+15           
        175  traj.desc2.states:mass_3                     c     0.000000E+00     9.110877E-01     6.666667E+15           
        176  traj.desc2.states:mass_4                     c     0.000000E+00     9.108801E-01     6.666667E+15           
        177  traj.desc2.states:mass_5                     c     0.000000E+00     9.106933E-01     6.666667E+15           
        178  traj.desc2.states:mass_6                     c     0.000000E+00     9.105714E-01     6.666667E+15           
        179  traj.desc2.states:mass_7                     c     0.000000E+00     9.105407E-01     6.666667E+15           
        180  traj.desc2.states:distance_0                 c     0.000000E+00     1.041456E+00     1.428571E+00           
        181  traj.desc2.states:distance_1                 c     0.000000E+00     1.042988E+00     1.428571E+00           
        182  traj.desc2.states:distance_2                 c     0.000000E+00     1.044980E+00     1.428571E+00           
        183  traj.desc2.states:distance_3                 c     0.000000E+00     1.047015E+00     1.428571E+00           
        184  traj.desc2.states:distance_4                 c     0.000000E+00     1.048704E+00     1.428571E+00           
        185  traj.desc2.states:distance_5                 c     0.000000E+00     1.049745E+00     1.428571E+00           
        186  traj.desc2.states:distance_6                 c     0.000000E+00     1.050000E+00     1.428571E+00           
        187  traj.groundroll.t_duration_0                 c     1.980198E-02     7.178866E-01     1.980198E+00           
        188  traj.groundroll.states:velocity_0            c     0.000000E+00     3.911834E-01     6.666667E+00           
        189  traj.groundroll.states:velocity_1            c     0.000000E+00     8.495790E-01     6.666667E+00           
        190  traj.groundroll.states:velocity_2            c     0.000000E+00     9.747302E-01     6.666667E+00           
        191  traj.groundroll.states:mass_0                c     0.000000E+00     1.159342E+00     6.666667E+15           
        192  traj.groundroll.states:mass_1                c     0.000000E+00     1.159034E+00     6.666667E+15           
        193  traj.groundroll.states:mass_2                c     0.000000E+00     1.158603E+00     6.666667E+15           
        194  traj.groundroll.states:mass_3                c     0.000000E+00     1.158466E+00     6.666667E+15           
        195  traj.groundroll.states:distance_0            c     0.000000E+00     2.205426E-01     3.333333E+00           
        196  traj.groundroll.states:distance_1            c     0.000000E+00     1.170685E+00     3.333333E+00           
        197  traj.groundroll.states:distance_2            c     0.000000E+00     1.605474E+00     3.333333E+00           
        198  traj.rotation.t_duration_0                   c     1.980198E-02     6.009772E-02     1.980198E+00           
        199  traj.rotation.states:alpha_0                 c     0.000000E+00     1.252525E-02     1.000000E+00           
        200  traj.rotation.states:alpha_1                 c     0.000000E+00     2.980743E-02     1.000000E+00           
        201  traj.rotation.states:alpha_2                 c     0.000000E+00     3.527723E-02     1.000000E+00           
        202  traj.rotation.states:velocity_0              c     0.000000E+00     1.496460E+00     1.000000E+01           
        203  traj.rotation.states:velocity_1              c     0.000000E+00     1.542757E+00     1.000000E+01           
        204  traj.rotation.states:velocity_2              c     0.000000E+00     1.557051E+00     1.000000E+01           
        205  traj.rotation.states:mass_0                  c     0.000000E+00     1.158466E+00     6.666667E+15           
        206  traj.rotation.states:mass_1                  c     0.000000E+00     1.158439E+00     6.666667E+15           
        207  traj.rotation.states:mass_2                  c     0.000000E+00     1.158402E+00     6.666667E+15           
        208  traj.rotation.states:mass_3                  c     0.000000E+00     1.158391E+00     6.666667E+15           
        209  traj.rotation.states:distance_0              c     0.000000E+00     1.017953E+00     2.000000E+00           
        210  traj.rotation.states:distance_1              c     0.000000E+00     1.095781E+00     2.000000E+00           
        211  traj.rotation.states:distance_2              c     0.000000E+00     1.120970E+00     2.000000E+00           

   Constraints (i - inequality, e - equality)
      Index  Name                                                          Type          Lower           Value           Upper    Status  Lagrange Multiplier (N/A)
          0  mission:constraints:range_residual                               e   0.000000E+00    1.334183E-09    0.000000E+00              9.00000E+100
          1  h_fit.h_init_gear                                                e   1.000000E+00   -3.116225E+02    1.000000E+00         E    9.00000E+100
          2  h_fit.h_init_flaps                                               e   1.000000E+00   -2.005189E+01    1.000000E+00         E    9.00000E+100
          3  groundroll_boundary.velocity                                     e   0.000000E+00    3.969998E-06    0.000000E+00         E    9.00000E+100
          4  gtow_constraint.GTOW                                             e   0.000000E+00   -4.370230E-08    0.000000E+00              9.00000E+100
          5  mission:constraints:mass_residual                                e   0.000000E+00    7.861158E-07    0.000000E+00              9.00000E+100
          6  taxi_groundroll_mass_constraint.mass                             e   0.000000E+00   -5.205045E-10    0.000000E+00              9.00000E+100
          7  traj.linkages.groundroll:mass_final|rotation:mass_initial        e   0.000000E+00   -8.876668E-09    0.000000E+00              9.00000E+100
          8  traj.linkages.rotation:mass_final|ascent:mass_initial            e   0.000000E+00   -1.203705E-06    0.000000E+00         E    9.00000E+100
          9  traj.linkages.rotation:alpha_final|ascent:alpha_initial          e   0.000000E+00    3.167567E-07    0.000000E+00              9.00000E+100
         10  traj.linkages.ascent:mass_final|accel:mass_initial               e   0.000000E+00   -1.920998E-06    0.000000E+00         E    9.00000E+100
         11  traj.linkages.accel:mass_final|climb1:mass_initial               e   0.000000E+00   -1.923239E-06    0.000000E+00         E    9.00000E+100
         12  traj.linkages.climb1:mass_final|climb2:mass_initial              e   0.000000E+00   -1.938926E-06    0.000000E+00         E    9.00000E+100
         13  traj.linkages.climb2:mass_final|cruise:mass_initial              e   0.000000E+00   -9.241749E-08    0.000000E+00              9.00000E+100
         14  traj.linkages.cruise:mass_final|desc1:mass_initial               e   0.000000E+00   -1.148250E-07    0.000000E+00              9.00000E+100
         15  traj.linkages.desc1:mass_final|desc2:mass_initial                e   0.000000E+00   -2.581917E-06    0.000000E+00         E    9.00000E+100
         16  traj.phases.accel->final_boundary_constraint->EAS                e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         17  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00    1.550687E-03    0.000000E+00         E    9.00000E+100
         18  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00   -7.076577E-05    0.000000E+00         E    9.00000E+100
         19  traj.accel.collocation_constraint.defects:velocity               e   0.000000E+00    7.856566E-05    0.000000E+00         E    9.00000E+100
         20  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00   -8.188246E-07    0.000000E+00              9.00000E+100
         21  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00   -3.794817E-07    0.000000E+00              9.00000E+100
         22  traj.accel.collocation_constraint.defects:mass                   e   0.000000E+00    9.533276E-08    0.000000E+00              9.00000E+100
         23  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    5.946486E-03    0.000000E+00         E    9.00000E+100
         24  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    3.302874E-03    0.000000E+00         E    9.00000E+100
         25  traj.accel.collocation_constraint.defects:distance               e   0.000000E+00    7.711292E-06    0.000000E+00         E    9.00000E+100
         26  traj.phases.ascent->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         27  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    1.373030E-02    0.000000E+00         E    9.00000E+100
         28  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    4.434143E-02    0.000000E+00         E    9.00000E+100
         29  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    3.844836E-02    0.000000E+00         E    9.00000E+100
         30  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    4.917956E-02    0.000000E+00         E    9.00000E+100
         31  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.504408E-02    0.000000E+00         E    9.00000E+100
         32  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.824375E-02    0.000000E+00         E    9.00000E+100
         33  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.088463E-02    0.000000E+00         E    9.00000E+100
         34  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    5.660427E-02    0.000000E+00         E    9.00000E+100
         35  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    6.082644E-02    0.000000E+00         E    9.00000E+100
         36  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    6.840249E-02    0.000000E+00         E    9.00000E+100
         37  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    7.407479E-02    0.000000E+00         E    9.00000E+100
         38  traj.ascent.collocation_constraint.defects:flight_path_angle     e   0.000000E+00    3.124106E-02    0.000000E+00         E    9.00000E+100
         39  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00    7.947686E-09    0.000000E+00              9.00000E+100
         40  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.069356E-02    0.000000E+00         E    9.00000E+100
         41  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -2.625838E-02    0.000000E+00         E    9.00000E+100
         42  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -3.395165E-02    0.000000E+00         E    9.00000E+100
         43  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -3.412323E-02    0.000000E+00         E    9.00000E+100
         44  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -2.002821E-02    0.000000E+00         E    9.00000E+100
         45  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.658741E-02    0.000000E+00         E    9.00000E+100
         46  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.622422E-02    0.000000E+00         E    9.00000E+100
         47  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -2.061930E-02    0.000000E+00         E    9.00000E+100
         48  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.969082E-02    0.000000E+00         E    9.00000E+100
         49  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.119609E-02    0.000000E+00         E    9.00000E+100
         50  traj.ascent.collocation_constraint.defects:altitude              e   0.000000E+00   -1.419851E-02    0.000000E+00         E    9.00000E+100
         51  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00   -3.536633E-02    0.000000E+00         E    9.00000E+100
         52  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00   -3.545162E-03    0.000000E+00         E    9.00000E+100
         53  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    7.048725E-03    0.000000E+00         E    9.00000E+100
         54  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    9.636142E-03    0.000000E+00         E    9.00000E+100
         55  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    8.881796E-03    0.000000E+00         E    9.00000E+100
         56  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.741589E-03    0.000000E+00         E    9.00000E+100
         57  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    2.744853E-03    0.000000E+00         E    9.00000E+100
         58  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    2.913446E-03    0.000000E+00         E    9.00000E+100
         59  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.851907E-03    0.000000E+00         E    9.00000E+100
         60  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    3.177131E-03    0.000000E+00         E    9.00000E+100
         61  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00   -8.791620E-06    0.000000E+00         E    9.00000E+100
         62  traj.ascent.collocation_constraint.defects:velocity              e   0.000000E+00    1.852124E-03    0.000000E+00         E    9.00000E+100
         63  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.626273E-07    0.000000E+00              9.00000E+100
         64  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    6.976613E-07    0.000000E+00              9.00000E+100
         65  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    6.361549E-07    0.000000E+00              9.00000E+100
         66  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.653630E-07    0.000000E+00              9.00000E+100
         67  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    5.767307E-07    0.000000E+00              9.00000E+100
         68  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    4.706536E-07    0.000000E+00              9.00000E+100
         69  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.964275E-07    0.000000E+00              9.00000E+100
         70  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    4.540411E-07    0.000000E+00              9.00000E+100
         71  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.788343E-07    0.000000E+00              9.00000E+100
         72  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    2.983472E-07    0.000000E+00              9.00000E+100
         73  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.856516E-07    0.000000E+00              9.00000E+100
         74  traj.ascent.collocation_constraint.defects:mass                  e   0.000000E+00    3.527655E-07    0.000000E+00              9.00000E+100
         75  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    2.361309E-03    0.000000E+00         E    9.00000E+100
         76  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    6.816899E-04    0.000000E+00         E    9.00000E+100
         77  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    3.731243E-04    0.000000E+00         E    9.00000E+100
         78  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    3.052153E-04    0.000000E+00         E    9.00000E+100
         79  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00    3.184056E-04    0.000000E+00         E    9.00000E+100
         80  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -1.559671E-04    0.000000E+00         E    9.00000E+100
         81  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -4.659108E-04    0.000000E+00         E    9.00000E+100
         82  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -1.131660E-03    0.000000E+00         E    9.00000E+100
         83  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -2.031716E-03    0.000000E+00         E    9.00000E+100
         84  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -2.337166E-03    0.000000E+00         E    9.00000E+100
         85  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -3.372726E-03    0.000000E+00         E    9.00000E+100
         86  traj.ascent.collocation_constraint.defects:distance              e   0.000000E+00   -4.815500E-03    0.000000E+00         E    9.00000E+100
         87  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00   -3.835667E-09    0.000000E+00              9.00000E+100
         88  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00    3.558047E-09    0.000000E+00              9.00000E+100
         89  traj.ascent.continuity_comp.defect_control_rates:alpha_rate      e   0.000000E+00    2.610795E-09    0.000000E+00              9.00000E+100
         90  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00   -7.039482E-09    0.000000E+00              9.00000E+100
         91  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00    2.385292E-08    0.000000E+00              9.00000E+100
         92  traj.ascent.continuity_comp.defect_controls:alpha                e   0.000000E+00    6.792369E-09    0.000000E+00              9.00000E+100
         93  traj.phases.climb1->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
         94  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00    2.146880E-04    0.000000E+00         E    9.00000E+100
         95  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00    4.426926E-04    0.000000E+00         E    9.00000E+100
         96  traj.climb1.collocation_constraint.defects:altitude              e   0.000000E+00   -1.554989E-04    0.000000E+00         E    9.00000E+100
         97  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -2.269606E-06    0.000000E+00         E    9.00000E+100
         98  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -2.444011E-06    0.000000E+00         E    9.00000E+100
         99  traj.climb1.collocation_constraint.defects:mass                  e   0.000000E+00   -8.848115E-07    0.000000E+00              9.00000E+100
        100  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -5.556437E-04    0.000000E+00         E    9.00000E+100
        101  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -7.399681E-04    0.000000E+00         E    9.00000E+100
        102  traj.climb1.collocation_constraint.defects:distance              e   0.000000E+00   -3.962516E-04    0.000000E+00         E    9.00000E+100
        103  traj.phases.climb2->final_boundary_constraint->altitude          e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        104  traj.phases.climb2->final_boundary_constraint->mach              e   8.000000E-01    8.000000E-01    8.000000E-01              9.00000E+100
        105  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -5.253283E-04    0.000000E+00         E    9.00000E+100
        106  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -5.406960E-03    0.000000E+00         E    9.00000E+100
        107  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -6.563124E-03    0.000000E+00         E    9.00000E+100
        108  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    3.956191E-03    0.000000E+00         E    9.00000E+100
        109  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -1.229840E-02    0.000000E+00         E    9.00000E+100
        110  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -1.319060E-02    0.000000E+00         E    9.00000E+100
        111  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    4.177165E-03    0.000000E+00         E    9.00000E+100
        112  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00    7.622489E-03    0.000000E+00         E    9.00000E+100
        113  traj.climb2.collocation_constraint.defects:altitude              e   0.000000E+00   -2.664431E-03    0.000000E+00         E    9.00000E+100
        114  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    2.873291E-08    0.000000E+00              9.00000E+100
        115  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    5.206070E-05    0.000000E+00         E    9.00000E+100
        116  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.339751E-04    0.000000E+00         E    9.00000E+100
        117  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00   -1.686830E-05    0.000000E+00         E    9.00000E+100
        118  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    8.453192E-05    0.000000E+00         E    9.00000E+100
        119  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.043711E-04    0.000000E+00         E    9.00000E+100
        120  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    1.935743E-05    0.000000E+00         E    9.00000E+100
        121  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00   -2.407175E-05    0.000000E+00         E    9.00000E+100
        122  traj.climb2.collocation_constraint.defects:mass                  e   0.000000E+00    2.069119E-05    0.000000E+00         E    9.00000E+100
        123  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    5.348640E-07    0.000000E+00              9.00000E+100
        124  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    2.191252E-04    0.000000E+00         E    9.00000E+100
        125  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.186274E-04    0.000000E+00         E    9.00000E+100
        126  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.220398E-04    0.000000E+00         E    9.00000E+100
        127  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.471591E-04    0.000000E+00         E    9.00000E+100
        128  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.720010E-04    0.000000E+00         E    9.00000E+100
        129  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    9.315702E-04    0.000000E+00         E    9.00000E+100
        130  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00    3.279655E-03    0.000000E+00         E    9.00000E+100
        131  traj.climb2.collocation_constraint.defects:distance              e   0.000000E+00   -2.490166E-04    0.000000E+00         E    9.00000E+100
        132  traj.phases.desc1->final_boundary_constraint->altitude           e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        133  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00    2.487153E-05    0.000000E+00         E    9.00000E+100
        134  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -7.843395E-05    0.000000E+00         E    9.00000E+100
        135  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -3.287796E-04    0.000000E+00         E    9.00000E+100
        136  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -3.728108E-04    0.000000E+00         E    9.00000E+100
        137  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -5.654851E-04    0.000000E+00         E    9.00000E+100
        138  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -2.294016E-03    0.000000E+00         E    9.00000E+100
        139  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00    9.664605E-05    0.000000E+00         E    9.00000E+100
        140  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -6.626720E-05    0.000000E+00         E    9.00000E+100
        141  traj.desc1.collocation_constraint.defects:altitude               e   0.000000E+00   -1.206803E-04    0.000000E+00         E    9.00000E+100
        142  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -8.733842E-09    0.000000E+00              9.00000E+100
        143  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -1.571074E-07    0.000000E+00              9.00000E+100
        144  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -1.952676E-07    0.000000E+00              9.00000E+100
        145  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.701606E-07    0.000000E+00              9.00000E+100
        146  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.679159E-07    0.000000E+00              9.00000E+100
        147  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -7.321711E-07    0.000000E+00              9.00000E+100
        148  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00   -2.644429E-07    0.000000E+00              9.00000E+100
        149  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00    7.666227E-08    0.000000E+00              9.00000E+100
        150  traj.desc1.collocation_constraint.defects:mass                   e   0.000000E+00    1.383872E-07    0.000000E+00              9.00000E+100
        151  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00   -5.888077E-05    0.000000E+00         E    9.00000E+100
        152  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    3.225136E-05    0.000000E+00         E    9.00000E+100
        153  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    8.597211E-05    0.000000E+00         E    9.00000E+100
        154  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    8.971158E-05    0.000000E+00         E    9.00000E+100
        155  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    9.282590E-05    0.000000E+00         E    9.00000E+100
        156  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    4.790127E-04    0.000000E+00         E    9.00000E+100
        157  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00   -1.569916E-04    0.000000E+00         E    9.00000E+100
        158  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00   -3.865078E-05    0.000000E+00         E    9.00000E+100
        159  traj.desc1.collocation_constraint.defects:distance               e   0.000000E+00    5.202309E-06    0.000000E+00         E    9.00000E+100
        160  traj.phases.desc2->final_boundary_constraint->altitude           e   1.000000E+00    1.000000E+00    1.000000E+00              9.00000E+100
        161  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    8.182355E-05    0.000000E+00         E    9.00000E+100
        162  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    1.791134E-04    0.000000E+00         E    9.00000E+100
        163  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.274549E-04    0.000000E+00         E    9.00000E+100
        164  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.507400E-04    0.000000E+00         E    9.00000E+100
        165  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    3.345589E-04    0.000000E+00         E    9.00000E+100
        166  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    1.950581E-04    0.000000E+00         E    9.00000E+100
        167  traj.desc2.collocation_constraint.defects:altitude               e   0.000000E+00    6.876125E-05    0.000000E+00         E    9.00000E+100
        168  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00    7.025804E-09    0.000000E+00              9.00000E+100
        169  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -7.706998E-08    0.000000E+00              9.00000E+100
        170  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -2.654911E-07    0.000000E+00              9.00000E+100
        171  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -1.360880E-07    0.000000E+00              9.00000E+100
        172  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -3.507615E-07    0.000000E+00              9.00000E+100
        173  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -2.022791E-07    0.000000E+00              9.00000E+100
        174  traj.desc2.collocation_constraint.defects:mass                   e   0.000000E+00   -3.132357E-08    0.000000E+00              9.00000E+100
        175  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00    6.481460E-08    0.000000E+00              9.00000E+100
        176  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -2.635014E-05    0.000000E+00         E    9.00000E+100
        177  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -6.596916E-05    0.000000E+00         E    9.00000E+100
        178  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -8.586096E-05    0.000000E+00         E    9.00000E+100
        179  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -7.235470E-05    0.000000E+00         E    9.00000E+100
        180  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -3.824407E-05    0.000000E+00         E    9.00000E+100
        181  traj.desc2.collocation_constraint.defects:distance               e   0.000000E+00   -8.258734E-06    0.000000E+00         E    9.00000E+100
        182  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -9.730877E-05    0.000000E+00         E    9.00000E+100
        183  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -1.016834E-03    0.000000E+00         E    9.00000E+100
        184  traj.groundroll.collocation_constraint.defects:velocity          e   0.000000E+00   -9.231414E-04    0.000000E+00         E    9.00000E+100
        185  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -1.804814E-09    0.000000E+00              9.00000E+100
        186  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -4.340941E-08    0.000000E+00              9.00000E+100
        187  traj.groundroll.collocation_constraint.defects:mass              e   0.000000E+00   -7.753064E-08    0.000000E+00              9.00000E+100
        188  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00   -1.323597E-09    0.000000E+00              9.00000E+100
        189  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00    8.249485E-03    0.000000E+00         E    9.00000E+100
        190  traj.groundroll.collocation_constraint.defects:distance          e   0.000000E+00    6.865935E-03    0.000000E+00         E    9.00000E+100
        191  traj.phases.rotation->final_boundary_constraint->normal_force    e   0.000000E+00    0.000000E+00    0.000000E+00              9.00000E+100
        192  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -7.145365E-08    0.000000E+00              9.00000E+100
        193  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -3.274435E-07    0.000000E+00              9.00000E+100
        194  traj.rotation.collocation_constraint.defects:alpha               e   0.000000E+00   -2.373049E-07    0.000000E+00              9.00000E+100
        195  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00   -3.968917E-04    0.000000E+00         E    9.00000E+100
        196  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00   -5.343371E-04    0.000000E+00         E    9.00000E+100
        197  traj.rotation.collocation_constraint.defects:velocity            e   0.000000E+00    2.226118E-04    0.000000E+00         E    9.00000E+100
        198  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.238554E-08    0.000000E+00              9.00000E+100
        199  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.309410E-08    0.000000E+00              9.00000E+100
        200  traj.rotation.collocation_constraint.defects:mass                e   0.000000E+00   -2.702920E-08    0.000000E+00              9.00000E+100
        201  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.031322E-03    0.000000E+00         E    9.00000E+100
        202  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.379952E-03    0.000000E+00         E    9.00000E+100
        203  traj.rotation.collocation_constraint.defects:distance            e   0.000000E+00    1.778089E-03    0.000000E+00         E    9.00000E+100
        204  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    6.818859E-01    1.100000E+00              9.00000E+100
        205  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.279411E-01    1.100000E+00              9.00000E+100
        206  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.471723E-01    1.100000E+00              9.00000E+100
        207  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.822907E-01    1.100000E+00              9.00000E+100
        208  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.822907E-01    1.100000E+00              9.00000E+100
        209  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.730462E-01    1.100000E+00              9.00000E+100
        210  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.938526E-01    1.100000E+00              9.00000E+100
        211  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.720754E-01    1.100000E+00              9.00000E+100
        212  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.720755E-01    1.100000E+00              9.00000E+100
        213  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.309773E-01    1.100000E+00              9.00000E+100
        214  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.558824E-01    1.100000E+00              9.00000E+100
        215  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.847562E-01    1.100000E+00              9.00000E+100
        216  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    4.847562E-01    1.100000E+00              9.00000E+100
        217  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    5.092822E-01    1.100000E+00              9.00000E+100
        218  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    3.703512E-01    1.100000E+00              9.00000E+100
        219  traj.phases.ascent->path_constraint->load_factor                 i   0.000000E+00    2.846567E-01    1.100000E+00              9.00000E+100
        220  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.606163E+00    1.500000E+01              9.00000E+100
        221  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    4.274662E+00    1.500000E+01              9.00000E+100
        222  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.334950E+00    1.500000E+01              9.00000E+100
        223  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.471303E+00    1.500000E+01              9.00000E+100
        224  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    2.471302E+00    1.500000E+01              9.00000E+100
        225  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    4.369775E+00    1.500000E+01              9.00000E+100
        226  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.642958E+00    1.500000E+01              9.00000E+100
        227  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.898426E+00    1.500000E+01              9.00000E+100
        228  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    6.898427E+00    1.500000E+01              9.00000E+100
        229  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.012013E+00    1.500000E+01              9.00000E+100
        230  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.970007E+00    1.500000E+01              9.00000E+100
        231  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.744976E+00    1.500000E+01              9.00000E+100
        232  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    8.744977E+00    1.500000E+01              9.00000E+100
        233  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    1.049505E+01    1.500000E+01              9.00000E+100
        234  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    9.525826E+00    1.500000E+01              9.00000E+100
        235  traj.phases.ascent->path_constraint->theta                       i   0.000000E+00    7.814118E+00    1.500000E+01              9.00000E+100
        236  traj.phases.climb2->final_boundary_constraint->altitude_rate     i   1.000000E-01    3.485365E+02    1.000000E+30              9.00000E+100

--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.

For FLOPS-derived models, let us run:

!aviary run_mission models/test_aircraft/aircraft_for_bench_FwFm.csv --max_iter 1 --optimizer IPOPT
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
Loaded default phase_info for height_energy equations of motion
User has specified Design.NUM_* passenger values but CrewPyaload.NUM_* has been left blank or set to zero.
Assuming they are equal to maintain backwards compatibility with GASP and FLOPS output files.
If you intended to have no passengers on this flight, please set Aircraft.CrewPayload.TOTAL_PAYLOAD_MASS to zero in aviary_values.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:897: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'climb': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'climb': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'cruise': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'cruise': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'descent': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2323: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'descent': lower, upper, ref
  warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
The following variables have been overridden:
  'aircraft:design:touchdown_mass  152800  lbm
  'aircraft:engine:mass  [7400.]  lbm
  'aircraft:fins:mass  0  lbm
  'aircraft:fuel:auxiliary_fuel_capacity  0  lbm
  'aircraft:fuel:fuselage_fuel_capacity  0  lbm
  'aircraft:fuel:total_capacity  45694  lbm
  'aircraft:fuselage:planform_area  1578.24  ft**2
  'aircraft:fuselage:wetted_area  4158.62  ft**2
  'aircraft:horizontal_tail:wetted_area  592.65  ft**2
  'aircraft:landing_gear:main_gear_oleo_length  102  inch
  'aircraft:landing_gear:nose_gear_oleo_length  67  inch
  'aircraft:vertical_tail:wetted_area  581.13  ft**2
  'aircraft:wing:aspect_ratio  11.22091  unitless
  'aircraft:wing:control_surface_area  137  ft**2
  'aircraft:wing:wetted_area  2396.56  ft**2
--- Constraint Report [traj] ---
    --- climb ---
        [path]    0.0000e+00 <= throttle <= 1.0000e+00  [unitless]
    --- cruise ---
        [initial] 0.0000e+00 <= throttle <= 1.0000e+00  [unitless]
        [final]   0.0000e+00 <= throttle <= 1.0000e+00  [unitless]
    --- descent ---
        [path]    0.0000e+00 <= throttle <= 1.0000e+00  [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:177: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:177: SolverWarning:DirectSolver in 'traj.phases.descent.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
Full total jacobian for problem 'aircraft_for_bench_FwFm' was computed 3 times, taking 0.5960236790000408 seconds.
Total jacobian shape: (137, 96) 
Jacobian shape: (137, 96)  (4.12% nonzero)
FWD solves: 5   REV solves: 0
Total colors vs. total size: 5 vs 96  (94.79% improvement)

Sparsity computed using tolerance: 1e-25
Time to compute sparsity:   0.5960 sec
Time to compute coloring:   0.1031 sec
Memory to compute coloring:   0.1250 MB
Coloring created on: 2024-12-17 21:06:02
This is Ipopt version 3.14.17, running with linear solver MUMPS 5.7.3.

Number of nonzeros in equality constraint Jacobian...:     2947
Number of nonzeros in inequality constraint Jacobian.:     1344
Number of nonzeros in Lagrangian Hessian.............:        0
Total number of variables............................:       96
                     variables with only lower bounds:        2
                variables with lower and upper bounds:       94
                     variables with only upper bounds:        0
Total number of equality constraints.................:       94
Total number of inequality constraints...............:       42
        inequality constraints with only lower bounds:        0
   inequality constraints with lower and upper bounds:       42
        inequality constraints with only upper bounds:        0

iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls
   0  3.0000000e-01 1.90e+02 1.91e-02  -5.0 0.00e+00    -  0.00e+00 0.00e+00   0
Number of Iterations....: 1

                                   (scaled)                 (unscaled)
Objective...............:   3.3334852721650099e-01    3.3334852721650099e-01
Dual infeasibility......:   1.0071386876135500e+00    1.0071386876135500e+00
Constraint violation....:   1.8729457092734259e+02    1.8729457092734259e+02
Variable bound violation:   0.0000000000000000e+00    0.0000000000000000e+00
Complementarity.........:   1.0000000827403710e+05    1.0000000827403710e+05
Overall NLP error.......:   1.0000000827403710e+05    1.0000000827403710e+05


Number of objective function evaluations             = 2
Number of objective gradient evaluations             = 2
Number of equality constraint evaluations            = 2
Number of inequality constraint evaluations          = 2
Number of equality constraint Jacobian evaluations   = 2
Number of inequality constraint Jacobian evaluations = 2
Number of Lagrangian Hessian evaluations             = 0
Total seconds in IPOPT                               = 1.058

EXIT: Maximum Number of Iterations Exceeded.
minimal_print is not available for this solution


Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
    Objective Function: _objfunc

    Solution: 
--------------------------------------------------------------------------------
    Total Time:                    1.0594
       User Objective Time :       0.7549
       User Sensitivity Time :     0.2878
       Interface Time :            0.0086
       Opt Solver Time:            0.0082
    Calls to Objective Function :       3
    Calls to Sens Function :            3


   Objectives
      Index  Name                               Value
          0  mission:objectives:fuel     3.333485E-01

   Variables (c - continuous, i - integer, d - discrete)
      Index  Name                              Type      Lower Bound            Value      Upper Bound     Status
          0  mission:design:gross_mass_0          c     5.714286E-05     1.001024E+00     1.000000E+30           
          1  mission:summary:gross_mass_0         c     5.714286E-05     1.001024E+00     1.000000E+30           
          2  traj.climb.t_duration_0              c     5.000000E-01     5.000001E-01     1.500000E+00          l
          3  traj.climb.states:mass_0             c     0.000000E+00     7.945996E+00     1.000000E+17           
          4  traj.climb.states:mass_1             c     0.000000E+00     7.946041E+00     1.000000E+17           
          5  traj.climb.states:mass_2             c     0.000000E+00     7.946087E+00     1.000000E+17           
          6  traj.climb.states:mass_3             c     0.000000E+00     7.946100E+00     1.000000E+17           
          7  traj.climb.states:mass_4             c     0.000000E+00     7.946147E+00     1.000000E+17           
          8  traj.climb.states:mass_5             c     0.000000E+00     7.946194E+00     1.000000E+17           
          9  traj.climb.states:mass_6             c     0.000000E+00     7.946208E+00     1.000000E+17           
         10  traj.climb.states:mass_7             c     0.000000E+00     7.946243E+00     1.000000E+17           
         11  traj.climb.states:mass_8             c     0.000000E+00     7.946293E+00     1.000000E+17           
         12  traj.climb.states:mass_9             c     0.000000E+00     7.946311E+00     1.000000E+17           
         13  traj.climb.states:mass_10            c     0.000000E+00     7.946346E+00     1.000000E+17           
         14  traj.climb.states:mass_11            c     0.000000E+00     7.946400E+00     1.000000E+17           
         15  traj.climb.states:mass_12            c     0.000000E+00     7.946417E+00     1.000000E+17           
         16  traj.climb.states:mass_13            c     0.000000E+00     7.946437E+00     1.000000E+17           
         17  traj.climb.states:mass_14            c     0.000000E+00     7.946465E+00     1.000000E+17           
         18  traj.climb.states:mass_15            c     0.000000E+00     7.946475E+00     1.000000E+17           
         19  traj.climb.states:distance_0         c     0.000000E+00     9.795121E-03     1.000000E+15           
         20  traj.climb.states:distance_1         c     0.000000E+00     9.703745E-03     1.000000E+15           
         21  traj.climb.states:distance_2         c     0.000000E+00     9.672563E-03     1.000000E+15           
         22  traj.climb.states:distance_3         c     0.000000E+00     9.512481E-03     1.000000E+15           
         23  traj.climb.states:distance_4         c     0.000000E+00     9.253729E-03     1.000000E+15           
         24  traj.climb.states:distance_5         c     0.000000E+00     9.162869E-03     1.000000E+15           
         25  traj.climb.states:distance_6         c     0.000000E+00     8.894159E-03     1.000000E+15           
         26  traj.climb.states:distance_7         c     0.000000E+00     8.473438E-03     1.000000E+15           
         27  traj.climb.states:distance_8         c     0.000000E+00     8.328541E-03     1.000000E+15           
         28  traj.climb.states:distance_9         c     0.000000E+00     8.033991E-03     1.000000E+15           
         29  traj.climb.states:distance_10        c     0.000000E+00     7.595005E-03     1.000000E+15           
         30  traj.climb.states:distance_11        c     0.000000E+00     7.448407E-03     1.000000E+15           
         31  traj.climb.states:distance_12        c     0.000000E+00     7.279755E-03     1.000000E+15           
         32  traj.climb.states:distance_13        c     0.000000E+00     7.039741E-03     1.000000E+15           
         33  traj.climb.states:distance_14        c     0.000000E+00     6.962037E-03     1.000000E+15           
         34  traj.cruise.t_duration_0             c     5.000000E-01     5.202490E-01     1.500000E+00           
         35  traj.cruise.states:mass_0            c     0.000000E+00     7.946170E+00     1.000000E+17           
         36  traj.cruise.states:mass_1            c     0.000000E+00     7.945750E+00     1.000000E+17           
         37  traj.cruise.states:mass_2            c     0.000000E+00     7.945617E+00     1.000000E+17           
         38  traj.cruise.states:mass_3            c     0.000000E+00     7.944996E+00     1.000000E+17           
         39  traj.cruise.states:mass_4            c     0.000000E+00     7.944145E+00     1.000000E+17           
         40  traj.cruise.states:mass_5            c     0.000000E+00     7.943876E+00     1.000000E+17           
         41  traj.cruise.states:mass_6            c     0.000000E+00     7.943147E+00     1.000000E+17           
         42  traj.cruise.states:mass_7            c     0.000000E+00     7.942148E+00     1.000000E+17           
         43  traj.cruise.states:mass_8            c     0.000000E+00     7.941833E+00     1.000000E+17           
         44  traj.cruise.states:mass_9            c     0.000000E+00     7.941228E+00     1.000000E+17           
         45  traj.cruise.states:mass_10           c     0.000000E+00     7.940397E+00     1.000000E+17           
         46  traj.cruise.states:mass_11           c     0.000000E+00     7.940135E+00     1.000000E+17           
         47  traj.cruise.states:mass_12           c     0.000000E+00     7.939842E+00     1.000000E+17           
         48  traj.cruise.states:mass_13           c     0.000000E+00     7.939438E+00     1.000000E+17           
         49  traj.cruise.states:mass_14           c     0.000000E+00     7.939310E+00     1.000000E+17           
         50  traj.cruise.states:distance_0        c     0.000000E+00     8.043272E-03     1.000000E+15           
         51  traj.cruise.states:distance_1        c     0.000000E+00     9.534506E-03     1.000000E+15           
         52  traj.cruise.states:distance_2        c     0.000000E+00     1.000632E-02     1.000000E+15           
         53  traj.cruise.states:distance_3        c     0.000000E+00     1.221181E-02     1.000000E+15           
         54  traj.cruise.states:distance_4        c     0.000000E+00     1.525221E-02     1.000000E+15           
         55  traj.cruise.states:distance_5        c     0.000000E+00     1.621383E-02     1.000000E+15           
         56  traj.cruise.states:distance_6        c     0.000000E+00     1.883020E-02     1.000000E+15           
         57  traj.cruise.states:distance_7        c     0.000000E+00     2.243641E-02     1.000000E+15           
         58  traj.cruise.states:distance_8        c     0.000000E+00     2.357684E-02     1.000000E+15           
         59  traj.cruise.states:distance_9        c     0.000000E+00     2.577210E-02     1.000000E+15           
         60  traj.cruise.states:distance_10       c     0.000000E+00     2.879838E-02     1.000000E+15           
         61  traj.cruise.states:distance_11       c     0.000000E+00     2.975553E-02     1.000000E+15           
         62  traj.cruise.states:distance_12       c     0.000000E+00     3.082835E-02     1.000000E+15           
         63  traj.cruise.states:distance_13       c     0.000000E+00     3.230796E-02     1.000000E+15           
         64  traj.cruise.states:distance_14       c     0.000000E+00     3.277610E-02     1.000000E+15           
         65  traj.descent.t_duration_0            c     5.000000E-01     5.503103E-01     1.500000E+00           
         66  traj.descent.states:mass_0           c     0.000000E+00     7.938939E+00     1.000000E+17           
         67  traj.descent.states:mass_1           c     0.000000E+00     7.938426E+00     1.000000E+17           
         68  traj.descent.states:mass_2           c     0.000000E+00     7.938263E+00     1.000000E+17           
         69  traj.descent.states:mass_3           c     0.000000E+00     7.937510E+00     1.000000E+17           
         70  traj.descent.states:mass_4           c     0.000000E+00     7.936481E+00     1.000000E+17           
         71  traj.descent.states:mass_5           c     0.000000E+00     7.936153E+00     1.000000E+17           
         72  traj.descent.states:mass_6           c     0.000000E+00     7.935266E+00     1.000000E+17           
         73  traj.descent.states:mass_7           c     0.000000E+00     7.934014E+00     1.000000E+17           
         74  traj.descent.states:mass_8           c     0.000000E+00     7.933601E+00     1.000000E+17           
         75  traj.descent.states:mass_9           c     0.000000E+00     7.932865E+00     1.000000E+17           
         76  traj.descent.states:mass_10          c     0.000000E+00     7.931945E+00     1.000000E+17           
         77  traj.descent.states:mass_11          c     0.000000E+00     7.931675E+00     1.000000E+17           
         78  traj.descent.states:mass_12          c     0.000000E+00     7.931383E+00     1.000000E+17           
         79  traj.descent.states:mass_13          c     0.000000E+00     7.930992E+00     1.000000E+17           
         80  traj.descent.states:mass_14          c     0.000000E+00     7.930869E+00     1.000000E+17           
         81  traj.descent.states:distance_0       c     0.000000E+00     3.425172E-02     1.000000E+15           
         82  traj.descent.states:distance_1       c     0.000000E+00     3.625133E-02     1.000000E+15           
         83  traj.descent.states:distance_2       c     0.000000E+00     3.687523E-02     1.000000E+15           
         84  traj.descent.states:distance_3       c     0.000000E+00     3.973432E-02     1.000000E+15           
         85  traj.descent.states:distance_4       c     0.000000E+00     4.351553E-02     1.000000E+15           
         86  traj.descent.states:distance_5       c     0.000000E+00     4.467136E-02     1.000000E+15           
         87  traj.descent.states:distance_6       c     0.000000E+00     4.771510E-02     1.000000E+15           
         88  traj.descent.states:distance_7       c     0.000000E+00     5.165885E-02     1.000000E+15           
         89  traj.descent.states:distance_8       c     0.000000E+00     5.284307E-02     1.000000E+15           
         90  traj.descent.states:distance_9       c     0.000000E+00     5.503469E-02     1.000000E+15           
         91  traj.descent.states:distance_10      c     0.000000E+00     5.786060E-02     1.000000E+15           
         92  traj.descent.states:distance_11      c     0.000000E+00     5.870595E-02     1.000000E+15           
         93  traj.descent.states:distance_12      c     0.000000E+00     5.962520E-02     1.000000E+15           
         94  traj.descent.states:distance_13      c     0.000000E+00     6.084330E-02     1.000000E+15           
         95  traj.descent.states:distance_14      c     0.000000E+00     6.121655E-02     1.000000E+15           

   Constraints (i - inequality, e - equality)
      Index  Name                                                      Type          Lower           Value           Upper    Status  Lagrange Multiplier (N/A)
          0  mission:constraints:range_residual                           e   0.000000E+00    1.872946E+02    0.000000E+00         E    9.00000E+100
          1  gtow_constraint.GTOW                                         e   0.000000E+00    0.000000E+00    0.000000E+00              9.00000E+100
          2  link_climb_mass.mass                                         e   0.000000E+00   -2.910383E-16    0.000000E+00              9.00000E+100
          3  mission:constraints:mass_residual                            e   0.000000E+00   -3.626541E-01    0.000000E+00         E    9.00000E+100
          4  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    2.982414E-04    0.000000E+00         E    9.00000E+100
          5  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    2.624315E-04    0.000000E+00         E    9.00000E+100
          6  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    2.297180E-04    0.000000E+00         E    9.00000E+100
          7  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    4.535355E-04    0.000000E+00         E    9.00000E+100
          8  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    4.024720E-04    0.000000E+00         E    9.00000E+100
          9  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    3.661309E-04    0.000000E+00         E    9.00000E+100
         10  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    4.281977E-04    0.000000E+00         E    9.00000E+100
         11  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    4.174994E-04    0.000000E+00         E    9.00000E+100
         12  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    4.169244E-04    0.000000E+00         E    9.00000E+100
         13  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    3.515382E-04    0.000000E+00         E    9.00000E+100
         14  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    3.547700E-04    0.000000E+00         E    9.00000E+100
         15  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    3.610187E-04    0.000000E+00         E    9.00000E+100
         16  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    1.778012E-04    0.000000E+00         E    9.00000E+100
         17  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    1.788765E-04    0.000000E+00         E    9.00000E+100
         18  traj.climb.collocation_constraint.defects:mass               e   0.000000E+00    1.811353E-04    0.000000E+00         E    9.00000E+100
         19  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00    9.214773E-05    0.000000E+00         E    9.00000E+100
         20  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -1.182565E-04    0.000000E+00         E    9.00000E+100
         21  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -2.043722E-04    0.000000E+00         E    9.00000E+100
         22  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -4.058403E-04    0.000000E+00         E    9.00000E+100
         23  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -4.701549E-04    0.000000E+00         E    9.00000E+100
         24  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -5.566144E-04    0.000000E+00         E    9.00000E+100
         25  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -6.936294E-04    0.000000E+00         E    9.00000E+100
         26  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -7.787535E-04    0.000000E+00         E    9.00000E+100
         27  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -8.922004E-04    0.000000E+00         E    9.00000E+100
         28  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -7.798135E-04    0.000000E+00         E    9.00000E+100
         29  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -8.352428E-04    0.000000E+00         E    9.00000E+100
         30  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -9.092146E-04    0.000000E+00         E    9.00000E+100
         31  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -4.563565E-04    0.000000E+00         E    9.00000E+100
         32  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -4.686959E-04    0.000000E+00         E    9.00000E+100
         33  traj.climb.collocation_constraint.defects:distance           e   0.000000E+00   -4.854155E-04    0.000000E+00         E    9.00000E+100
         34  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.270193E-04    0.000000E+00         E    9.00000E+100
         35  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.269269E-04    0.000000E+00         E    9.00000E+100
         36  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.267778E-04    0.000000E+00         E    9.00000E+100
         37  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.588151E-04    0.000000E+00         E    9.00000E+100
         38  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.583520E-04    0.000000E+00         E    9.00000E+100
         39  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.577400E-04    0.000000E+00         E    9.00000E+100
         40  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    3.062049E-04    0.000000E+00         E    9.00000E+100
         41  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    3.056177E-04    0.000000E+00         E    9.00000E+100
         42  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    3.048525E-04    0.000000E+00         E    9.00000E+100
         43  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.562208E-04    0.000000E+00         E    9.00000E+100
         44  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.558582E-04    0.000000E+00         E    9.00000E+100
         45  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    2.553844E-04    0.000000E+00         E    9.00000E+100
         46  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.249781E-04    0.000000E+00         E    9.00000E+100
         47  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.248948E-04    0.000000E+00         E    9.00000E+100
         48  traj.cruise.collocation_constraint.defects:mass              e   0.000000E+00    1.247799E-04    0.000000E+00         E    9.00000E+100
         49  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.331650E-04    0.000000E+00         E    9.00000E+100
         50  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.330063E-04    0.000000E+00         E    9.00000E+100
         51  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.327872E-04    0.000000E+00         E    9.00000E+100
         52  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.837327E-04    0.000000E+00         E    9.00000E+100
         53  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.830699E-04    0.000000E+00         E    9.00000E+100
         54  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.821545E-04    0.000000E+00         E    9.00000E+100
         55  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -1.048450E-03    0.000000E+00         E    9.00000E+100
         56  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -1.047511E-03    0.000000E+00         E    9.00000E+100
         57  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -1.046214E-03    0.000000E+00         E    9.00000E+100
         58  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.796386E-04    0.000000E+00         E    9.00000E+100
         59  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.789728E-04    0.000000E+00         E    9.00000E+100
         60  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -8.780532E-04    0.000000E+00         E    9.00000E+100
         61  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.297942E-04    0.000000E+00         E    9.00000E+100
         62  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.296343E-04    0.000000E+00         E    9.00000E+100
         63  traj.cruise.collocation_constraint.defects:distance          e   0.000000E+00   -4.294135E-04    0.000000E+00         E    9.00000E+100
         64  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    2.834206E-05    0.000000E+00         E    9.00000E+100
         65  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    2.808045E-05    0.000000E+00         E    9.00000E+100
         66  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    2.767147E-05    0.000000E+00         E    9.00000E+100
         67  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    5.620727E-05    0.000000E+00         E    9.00000E+100
         68  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    5.534016E-05    0.000000E+00         E    9.00000E+100
         69  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    5.330763E-05    0.000000E+00         E    9.00000E+100
         70  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    6.249366E-05    0.000000E+00         E    9.00000E+100
         71  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    5.980619E-05    0.000000E+00         E    9.00000E+100
         72  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    5.467559E-05    0.000000E+00         E    9.00000E+100
         73  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    4.448944E-05    0.000000E+00         E    9.00000E+100
         74  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    4.124631E-05    0.000000E+00         E    9.00000E+100
         75  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    3.964390E-05    0.000000E+00         E    9.00000E+100
         76  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    1.935859E-05    0.000000E+00         E    9.00000E+100
         77  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    1.950722E-05    0.000000E+00         E    9.00000E+100
         78  traj.descent.collocation_constraint.defects:mass             e   0.000000E+00    1.991093E-05    0.000000E+00         E    9.00000E+100
         79  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -2.202914E-04    0.000000E+00         E    9.00000E+100
         80  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -2.170380E-04    0.000000E+00         E    9.00000E+100
         81  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -2.124420E-04    0.000000E+00         E    9.00000E+100
         82  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -4.308439E-04    0.000000E+00         E    9.00000E+100
         83  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -4.163778E-04    0.000000E+00         E    9.00000E+100
         84  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -3.955410E-04    0.000000E+00         E    9.00000E+100
         85  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -4.621715E-04    0.000000E+00         E    9.00000E+100
         86  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -4.395629E-04    0.000000E+00         E    9.00000E+100
         87  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -4.069757E-04    0.000000E+00         E    9.00000E+100
         88  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -3.333599E-04    0.000000E+00         E    9.00000E+100
         89  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -3.157475E-04    0.000000E+00         E    9.00000E+100
         90  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -2.906597E-04    0.000000E+00         E    9.00000E+100
         91  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -1.383413E-04    0.000000E+00         E    9.00000E+100
         92  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -1.338275E-04    0.000000E+00         E    9.00000E+100
         93  traj.descent.collocation_constraint.defects:distance         e   0.000000E+00   -1.275099E-04    0.000000E+00         E    9.00000E+100
         94  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    5.643333E-01    1.000000E+00              9.00000E+100
         95  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    5.241742E-01    1.000000E+00              9.00000E+100
         96  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.898411E-01    1.000000E+00              9.00000E+100
         97  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.825034E-01    1.000000E+00              9.00000E+100
         98  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.825034E-01    1.000000E+00              9.00000E+100
         99  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.658966E-01    1.000000E+00              9.00000E+100
        100  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.783067E-01    1.000000E+00              9.00000E+100
        101  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.871277E-01    1.000000E+00              9.00000E+100
        102  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    4.871277E-01    1.000000E+00              9.00000E+100
        103  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    5.273592E-01    1.000000E+00              9.00000E+100
        104  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    6.055021E-01    1.000000E+00              9.00000E+100
        105  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    6.360401E-01    1.000000E+00              9.00000E+100
        106  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    6.360401E-01    1.000000E+00              9.00000E+100
        107  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    6.994589E-01    1.000000E+00              9.00000E+100
        108  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    8.069827E-01    1.000000E+00              9.00000E+100
        109  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    8.450170E-01    1.000000E+00              9.00000E+100
        110  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    8.450170E-01    1.000000E+00              9.00000E+100
        111  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    8.920666E-01    1.000000E+00              9.00000E+100
        112  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    9.560546E-01    1.000000E+00              9.00000E+100
        113  traj.phases.climb->path_constraint->throttle                 i   0.000000E+00    9.771954E-01    1.000000E+00              9.00000E+100
        114  traj.phases.cruise->initial_boundary_constraint->throttle    i   0.000000E+00    7.736965E-01    1.000000E+00              9.00000E+100
        115  traj.phases.cruise->final_boundary_constraint->throttle      i   0.000000E+00    8.431198E-01    1.000000E+00              9.00000E+100
        116  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.780837E-01    1.000000E+00              9.00000E+100
        117  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.502843E-01    1.000000E+00              9.00000E+100
        118  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.172868E-01    1.000000E+00              9.00000E+100
        119  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.079391E-01    1.000000E+00              9.00000E+100
        120  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.079391E-01    1.000000E+00              9.00000E+100
        121  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    2.647135E-01    1.000000E+00              9.00000E+100
        122  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    2.112276E-01    1.000000E+00              9.00000E+100
        123  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.957721E-01    1.000000E+00              9.00000E+100
        124  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.957721E-01    1.000000E+00              9.00000E+100
        125  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.592641E-01    1.000000E+00              9.00000E+100
        126  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.194467E-01    1.000000E+00              9.00000E+100
        127  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.082926E-01    1.000000E+00              9.00000E+100
        128  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    1.082926E-01    1.000000E+00              9.00000E+100
        129  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    8.760524E-02    1.000000E+00              9.00000E+100
        130  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    4.650921E-02    1.000000E+00              9.00000E+100
        131  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.450878E-02    1.000000E+00              9.00000E+100
        132  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    3.450878E-02    1.000000E+00              9.00000E+100
        133  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    2.169453E-02    1.000000E+00              9.00000E+100
        134  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    5.537003E-03    1.000000E+00              9.00000E+100
        135  traj.phases.descent->path_constraint->throttle               i   0.000000E+00    6.889462E-04    1.000000E+00              9.00000E+100

--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.

Aviary fortran_to_aviary command#

Note

You only need to use the fortran_to_aviary command if you are converting a FLOPS or GASP deck to Aviary’s csv format.

This fortran_to_aviary tool converts FORTRAN namelists into Aviary’s csv based format using the mappings found in the historical_name section of the variable_meta_data. The resulting csv is automatically sorted into three sections:

  1. Input Values: Any FORTRAN variables that were mapped to input variables in Aviary components

  2. Initial Guesses: Some variables are only used as initial guesses for the trajectory. These are displayed separately from the Input Values because they will not be passed directly to components

  3. Unconverted Values: If the fortran_to_aviary converter can’t find an Aviary variable that matches the FORTRAN variable, it is added to the end of the csv file. We recommend that you check this section after converting a namelist to ensure that there aren’t any variables you expected to be converted here. Many of these unconverted variables represent features or options that are not used in Aviary and can be safely ignored. Variables related to mission definition are important, but Aviary defines mission profiles in a significantly different way. Currently, the user must build a new mission definition file that recreates the mission. Aviary will ignore unconverted variables when loading the csv, so you can safely leave them.

# Testing Cell
from aviary.interface.cmd_entry_points import _command_map
_command_map['fortran_to_aviary'];

Phase info#

Aviary runs depend not only on input *.csv files, but also on phase_info dictionaries. A phase_info dictionary defines the mission and any analysis settings used within Aviary. It is used by Aviary to build a mission trajectory (e.g., from take-off to landing).

A Python dictionary is a set of key-value pairs. Most keys in phase_info are self-explained. Users should read Dymos’ Phase Options. For more details about phase_info keys, especially their example usages, please read onboarding phase information.

Default GASP phases#

The default_phase_info/two_dof.py file is shown below:

phase_info = {
    'groundroll': {
        'user_options': {
            'num_segments': 1,
            'order': 3,
            'connect_initial_mass': False,
            'fix_initial': True,
            'fix_initial_mass': False,
            'duration_ref': (50., 's'),
            'duration_bounds': ((1., 100.), 's'),
            'velocity_lower': (0, 'kn'),
            'velocity_upper': (1000, 'kn'),
            'velocity_ref': (150, 'kn'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'ft'),
            'distance_upper': (10.e3, 'ft'),
            'distance_ref': (3000, 'ft'),
            'distance_defect_ref': (3000, 'ft'),
        },
        'initial_guesses': {
            'time': ([0.0, 40.0], 's'),
            'velocity': ([0.066, 143.1], 'kn'),
            'distance': ([0.0, 1000.], 'ft'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'rotation': {
        'user_options': {
            'num_segments': 1,
            'order': 3,
            'fix_initial': False,
            'duration_bounds': ((1, 100), 's'),
            'duration_ref': (50.0, 's'),
            'velocity_lower': (0, 'kn'),
            'velocity_upper': (1000, 'kn'),
            'velocity_ref': (100, 'kn'),
            'velocity_ref0': (0, 'kn'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'ft'),
            'distance_upper': (10_000, 'ft'),
            'distance_ref': (5000, 'ft'),
            'distance_defect_ref': (5000, 'ft'),
            'angle_lower': (0., 'rad'),
            'angle_upper': (5., 'rad'),
            'angle_ref': (5., 'rad'),
            'angle_defect_ref': (5., 'rad'),
            'normal_ref': (10000, 'lbf'),
        },
        'initial_guesses': {
            'time': ([40.0, 5.0], 's'),
            'alpha': ([0.0, 2.5], 'deg'),
            'velocity': ([143, 150.], 'kn'),
            'distance': ([3680.37217765, 4000], 'ft'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'ascent': {
        'user_options': {
            'num_segments': 4,
            'order': 3,
            'fix_initial': False,
            'velocity_lower': (0, 'kn'),
            'velocity_upper': (700, 'kn'),
            'velocity_ref': (200, 'kn'),
            'velocity_ref0': (0, 'kn'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'ft'),
            'distance_upper': (15_000, 'ft'),
            'distance_ref': (1e4, 'ft'),
            'distance_defect_ref': (1e4, 'ft'),
            'alt_lower': (0, 'ft'),
            'alt_upper': (700, 'ft'),
            'alt_ref': (1000, 'ft'),
            'alt_defect_ref': (1000, 'ft'),
            'final_altitude': (500, 'ft'),
            'alt_constraint_ref': (500, 'ft'),
            'angle_lower': (-10, 'rad'),
            'angle_upper': (20, 'rad'),
            'angle_ref': (57.2958, 'deg'),
            'angle_defect_ref': (57.2958, 'deg'),
            'pitch_constraint_lower': (0., 'deg'),
            'pitch_constraint_upper': (15., 'deg'),
            'pitch_constraint_ref': (1., 'deg'),
        },
        'initial_guesses': {
            'time': ([45., 25.], 's'),
            'flight_path_angle': ([0.0, 8.], 'deg'),
            'alpha': ([2.5, 1.5], 'deg'),
            'velocity': ([150., 185.], 'kn'),
            'distance': ([4.e3, 10.e3], 'ft'),
            'altitude': ([0.0, 500.], 'ft'),
            'tau_gear': (0.2, 'unitless'),
            'tau_flaps': (0.9, 'unitless'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'accel': {
        'user_options': {
            'num_segments': 1,
            'order': 3,
            'fix_initial': False,
            'alt': (500, 'ft'),
            'EAS_constraint_eq': (250, 'kn'),
            'duration_bounds': ((1, 200), 's'),
            'duration_ref': (1000, 's'),
            'velocity_lower': (150, 'kn'),
            'velocity_upper': (270, 'kn'),
            'velocity_ref': (250, 'kn'),
            'velocity_ref0': (150, 'kn'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'NM'),
            'distance_upper': (150, 'NM'),
            'distance_ref': (5, 'NM'),
            'distance_defect_ref': (5, 'NM'),
        },
        'initial_guesses': {
            'time': ([70., 13.], 's'),
            'velocity': ([185., 250.], 'kn'),
            'distance': ([10.e3, 20.e3], 'ft'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'climb1': {
        'user_options': {
            'num_segments': 1,
            'order': 3,
            'fix_initial': False,
            'EAS_target': (250, 'kn'),
            'mach_cruise': 0.8,
            'target_mach': False,
            'final_altitude': (10.e3, 'ft'),
            'duration_bounds': ((30, 300), 's'),
            'duration_ref': (1000, 's'),
            'alt_lower': (400, 'ft'),
            'alt_upper': (11_000, 'ft'),
            'alt_ref': (10.e3, 'ft'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'NM'),
            'distance_upper': (500, 'NM'),
            'distance_ref': (10, 'NM'),
            'distance_ref0': (0, 'NM'),
        },
        'initial_guesses': {
            'time': ([1., 2.], 'min'),
            'distance': ([20.e3, 100.e3], 'ft'),
            'altitude': ([500., 10.e3], 'ft'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'climb2': {
        'user_options': {
            'num_segments': 3,
            'order': 3,
            'fix_initial': False,
            'EAS_target': (270, 'kn'),
            'mach_cruise': 0.8,
            'target_mach': True,
            'final_altitude': (37.5e3, 'ft'),
            'required_available_climb_rate': (0.1, 'ft/min'),
            'duration_bounds': ((200, 17_000), 's'),
            'duration_ref': (5000, 's'),
            'alt_lower': (9000, 'ft'),
            'alt_upper': (40000, 'ft'),
            'alt_ref': (30000, 'ft'),
            'alt_ref0': (0, 'ft'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (10, 'NM'),
            'distance_upper': (1000, 'NM'),
            'distance_ref': (500, 'NM'),
            'distance_ref0': (0, 'NM'),
            'distance_defect_ref': (500, 'NM'),
        },
        'initial_guesses': {
            'time': ([216., 1300.], 's'),
            'distance': ([100.e3, 200.e3], 'ft'),
            'altitude': ([10.e3, 37.5e3], 'ft'),
            'throttle': ([0.956, 0.956], 'unitless'),
        }
    },
    'cruise': {
        'user_options': {
            'alt_cruise': (37.5e3, 'ft'),
            'mach_cruise': 0.8,
        },
        'initial_guesses': {
            # [Initial mass, delta mass] for special cruise phase.
            'mass': ([171481., -35000], 'lbm'),
            'initial_distance': (200.e3, 'ft'),
            'initial_time': (1516., 's'),
            'altitude': (37.5e3, 'ft'),
            'mach': (0.8, 'unitless'),
        }
    },
    'desc1': {
        'user_options': {
            'num_segments': 3,
            'order': 3,
            'fix_initial': False,
            'input_initial': False,
            'EAS_limit': (350, 'kn'),
            'mach_cruise': 0.8,
            'input_speed_type': SpeedType.MACH,
            'final_altitude': (10.e3, 'ft'),
            'duration_bounds': ((300., 900.), 's'),
            'duration_ref': (1000, 's'),
            'alt_lower': (1000, 'ft'),
            'alt_upper': (40_000, 'ft'),
            'alt_ref': (30_000, 'ft'),
            'alt_ref0': (0, 'ft'),
            'alt_constraint_ref': (10000, 'ft'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (140_000, 'lbm'),
            'mass_ref0': (0, 'lbm'),
            'mass_defect_ref': (140_000, 'lbm'),
            'distance_lower': (3000., 'NM'),
            'distance_upper': (5000., 'NM'),
            'distance_ref': (mission_distance, 'NM'),
            'distance_ref0': (0, 'NM'),
            'distance_defect_ref': (100, 'NM'),
        },
        'initial_guesses': {
            'mass': (136000., 'lbm'),
            'altitude': ([37.5e3, 10.e3], 'ft'),
            'throttle': ([0.0, 0.0], 'unitless'),
            'distance': ([.92*mission_distance, .96*mission_distance], 'NM'),
            'time': ([28000., 500.], 's'),
        }
    },
    'desc2': {
        'user_options': {
            'num_segments': 1,
            'order': 7,
            'fix_initial': False,
            'input_initial': False,
            'EAS_limit': (250, 'kn'),
            'mach_cruise': 0.80,
            'input_speed_type': SpeedType.EAS,
            'final_altitude': (1000, 'ft'),
            'duration_bounds': ((100., 5000), 's'),
            'duration_ref': (500, 's'),
            'alt_lower': (500, 'ft'),
            'alt_upper': (11_000, 'ft'),
            'alt_ref': (10.e3, 'ft'),
            'alt_ref0': (1000, 'ft'),
            'alt_constraint_ref': (1000, 'ft'),
            'mass_lower': (0, 'lbm'),
            'mass_upper': (None, 'lbm'),
            'mass_ref': (150_000, 'lbm'),
            'mass_defect_ref': (150_000, 'lbm'),
            'distance_lower': (0, 'NM'),
            'distance_upper': (5000, 'NM'),
            'distance_ref': (3500, 'NM'),
            'distance_defect_ref': (100, 'NM'),
        },
        'initial_guesses': {
            'mass': (136000., 'lbm'),
            'altitude': ([10.e3, 1.e3], 'ft'),
            'throttle': ([0., 0.], 'unitless'),
            'distance': ([.96*mission_distance, mission_distance], 'NM'),
            'time': ([28500., 500.], 's'),
        }
    },
}


def phase_info_parameterization(phase_info, post_mission_info, aviary_inputs):
    """
    Modify the values in the phase_info dictionary to accomodate different values
    for the following mission design inputs: cruise altitude, cruise mach number,
    cruise range, design gross mass.

    Parameters
    ----------
    phase_info : dict
        Dictionary of phase settings for a mission profile
    post_mission_info : dict
        Dictionary of phase settings for a post mission profile
    aviary_inputs : <AviaryValues>
        Object containing values and units for all aviary inputs and options

    Returns
    -------
    dict
        Modified phase_info that has been changed to match the new mission
        parameters
    """

    range_cruise = aviary_inputs.get_val(Mission.Design.RANGE, units='NM')
    alt_cruise = aviary_inputs.get_val(Mission.Design.CRUISE_ALTITUDE, units='ft')
    gross_mass = aviary_inputs.get_val(Mission.Design.GROSS_MASS, units='lbm')
    mach_cruise = aviary_inputs.get_val(Mission.Design.MACH)

    # Range
    old_range_cruise = phase_info['desc2']['initial_guesses']['distance'][0][1]
    range_scale = 1.0
    if range_cruise != old_range_cruise:

        phase_info['desc1']['initial_guesses']['distance'] = \
            ([.92*range_cruise, .96*range_cruise], 'NM')
        phase_info['desc2']['initial_guesses']['distance'] = \
            ([.96*range_cruise, range_cruise], 'NM')
        range_scale = range_cruise / old_range_cruise

    # Altitude
    old_alt_cruise = phase_info['climb2']['user_options']['final_altitude'][0]
    if alt_cruise != old_alt_cruise:

        phase_info['climb2']['user_options']['final_altitude'] = (alt_cruise, 'ft')
        phase_info['climb2']['initial_guesses']['altitude'] = ([10.e3, alt_cruise], 'ft')
        phase_info['cruise']['initial_guesses']['altitude'] = (alt_cruise, 'ft')
        phase_info['desc1']['initial_guesses']['altitude'] = ([alt_cruise, 10.e3], 'ft')

        # TODO - Could adjust time guesses/bounds in climb2 and desc2.

    # Mass
    old_gross_mass = 175400.0
    if gross_mass != old_gross_mass:

        # Note, this requires that the guess for gross mass is pretty close to the
        # compute mass.

        fuel_used = 35000 * range_scale
        phase_info['groundroll']['initial_guesses']['mass'] = \
            ([gross_mass, gross_mass], 'lbm')
        phase_info['rotation']['initial_guesses']['mass'] = \
            ([gross_mass, gross_mass], 'lbm')
        phase_info['accel']['initial_guesses']['mass'] = \
            ([gross_mass, gross_mass], 'lbm')
        phase_info['ascent']['initial_guesses']['mass'] = \
            ([gross_mass, gross_mass], 'lbm')

        phase_info['cruise']['initial_guesses']['mass'] = \
            ([gross_mass, -fuel_used], 'lbm')

        end_mass = gross_mass - fuel_used
        phase_info['desc1']['initial_guesses']['mass'] = (end_mass, 'lbm')
        phase_info['desc2']['initial_guesses']['mass'] = (end_mass, 'lbm')

    # Mach
    old_mach_cruise = phase_info['cruise']['initial_guesses']['mach'][0]
    if mach_cruise != old_mach_cruise:

        phase_info['cruise']['initial_guesses']['mach'] = (mach_cruise, 'unitless')

    return phase_info, post_mission_info

The file default_phase_info/two_dof.py contains the following phases:

groundroll, rotation, ascent, accel, climb1, climb2, cruise, desc1, desc2.

All of the above phases belong to mission. No pre-mission phase is provided. If pre_mission is missing, a default pre_mission is provided:

    'pre_mission': {
        'include_takeoff': True,
        'external_subsystems': [],
    }

Similarly, if no post-mission phase is provided, then a default post_mission is provided:

    'post_mission': {
        'include_landing': True,
        'external_subsystems': [],
    },

For `GASP’ missions, taxi is considered part of pre-mission and landing is considered part of post-mission.

Groundroll phase#

Let us discuss the groundroll phase in detail as an example.

In groundroll phase, we are given duration_bounds of ((1., 100.), 's') and a reference value duration_ref of (50., 's'). Here, duration_bounds is a tuple of (lower, upper) bounds for the duration of the integration variable across the phase and duration_ref is the unit-reference value for the duration of the integration variable across the phase (see Dymos Variables). This implies a time range of (0.2, 2) in ground roll phase, or:

0.02 ≤ traj.groundroll.t_duration ≤ 2.

We see velocity_lower: (0, 'kn'), velocity_upper: (1000, 'kn'), and velocity_ref: (100, 'kn'). They result in TAS as a design variable with a range: (0, 10), or

0 ≤ traj.groundroll.states:TAS ≤ 10.

Similarly, we get:

0 ≤ traj.groundroll.states:distance ≤ 3.3333

The mass parameter is a little bit different because we have mass_upper: (None, 'lbm') (and mass_lower: (0, 'lbm')). In this situation, OpenMDAO provides replaces the None with an actual upper bound that is very large. We get:

0 ≤ traj.groundroll.states:mass ≤ 5.7e+15.

Comparing to TAS and distance, mass is not scaled to the same range. This is okay because aircraft mass is actually a constant in this phase. A huge upper bound will not have an impact.

In ground roll phase, throttle setting is set to maximum (1.0). Aviary sets a phase parameter:

Dynamic.Vehicle.Propulsion.THROTTLE = 1.0

For the COLLOCATION setting, there is one segment ('num_segments': 1) and polynomial interpolation degree is 3 ('order': 3). Increasing the number of segments and/or increasing the degree of polynomial will improve accuracy but will also increase the complexity of computation. For groundroll, it is unnecessary.

Since groundroll is the first phase, connect_initial_mass is set to False. fix_initial is set to True, but fix_initial_mass to False because aircraft total mass is a design parameter.

Parameter mass_defect_ref: (10.e3, 'lbm') (or mass_defect_ref: (1.e4, 'lbm')) lbm. Defect is a residual that measures how accurately the proposed state and control history obeyed the ODE governing the system dynamics. For distance, the defect reference is 100 ft. For TAS, defect_ref is not provided. The default behavior is to make make the defect_ref equal to the ref value. Please see here for more info.

Other phases#

We will not discuss the other phases in detail.

Phase climb1 is for climb up to 10,000 ft and phase climb2 is for climb to cruise phase. Phase desc1 is for descent down to 10,000 ft and phase desc2 is for descent from 10,000 ft down to 1,000 ft.

Default FLOPS phases#

The file default_phase_info/flops.py contains the following phases:

pre_mission, climb, cruise, descent, post_mission

The differences between GASP and FLOPS phases are due to how GASP and FLOPS implement trajectory analysis.

Note

File default_phase_info/flops.py has a pre_mission phase and a post_mission phase. In pre_mission, takeoff is the simplified takeoff and in post_mission, landing is the simplified landing. For FLOPS missions, there are detailed takeoff and landing available. But they are not used in default_phase_info/flops.py. The other phases are mission phases.

Outputted files and how to read them#

Note

Please also see this doc page for more info about the Aviary dashboard and the information is contains. We usually expect that users access these reports through the dashboard, though you can also access the html files directly.

Each standard run generates several output files. Which output files are generated depend on the run options. In this section, we assume that we’ve set max_iter = 50 to ensure convergence. The following screenshots used in this article are all from a run using aircraft_for_bench_GwGm.csv.

First, there is always a sub-folder reports/aviary that contains a few HTML files.

  • driver_scaling_report.html

  • inputs.html

  • n2.html

  • opt_report.html

  • total_coloring.html

  • traj_linkage_report.html

  • traj_results_report.html

File driver_scaling_report.html#

This file is a summary of driver scaling information. After all design variables, objectives, and constraints are declared and the problem has been set up, this report presents all the design variables and constraints in all phases as well as the objectives. The file is divided to three blocks: Design Variables, Constraints, and Objectives. It contains the following columns: name, size, indices, driver value and units, model value and units, ref, ref0, scaler, adder, etc. It also shows Jacobian Info - responses with respect to design variables (DV). A screen shot of design variables is shown here:

Sample driver_scaling_report.html (top lines)

This file is needed when you are debugging. New users can skip it.

File inputs.html#

File inputs.html is a sortable and filterable input report of input variables in different phases. It contains all the input variables (with possibly duplicate names for different phases) but only shows those in Source is IVC (abbreviation for IndepVarComp) on opening up. Users can choose to show other inputs by selecting and deselecting this checkbox. Users can filter their inputs by input name, source name, units, shape, tags, values (Val), or design variables (Source is DV). Here is a screen shot (top part) when it is opened:

Sample inputs.html (top lines)

New users can choose to use input_list.txt instead (see below).

File n2.html#

N2, sometimes referred to as an eXtended Design Structure Matrix (XDSM), is a powerful tool for understanding your model in OpenMDAO. It is an N-squared diagram in the shape of a matrix representing functional or physical interfaces between system elements. It can be used to systematically identify, define, tabulate, design, and analyze functional and physical interfaces.

Note

We strongly recommend that you understand N2 diagrams well, especially if you are debugging a model or adding external subsystems. Here is a doc page featuring an informative video on how to use N2 diagrams effectively. MDO lab has another resource for learning about N2 and XDSM.

Here is a screenshot of an N2 diagram after we run aviary run_mission models/test_aircraft/aircraft_for_bench_GwGm.csv.

Sample n2.html (on start)

At level 1, we are interested in the phases, so let us Hide solvers. So from model, zoom into traj, then into phases. We see the phases from groundroll to desc2 as expected. To see how those phases are connected, we click on Show all connections in view. Here is what we see:

Sample n2.html (on phases)

The solid arrow connections show how outputs of one phase feed as inputs to the next. As you can see, the input of each phase is linked from the previous phase and its output is linked to the next phase. This is pretty much what we expect. The dashed arrows are links to and/or from other places not in the view. If you are curious where those links go, you must zoom out. This can be done by clicking the vertical bar on the left. Let us click it twice. Here is what we get:

Sample n2.html (on phases)

Note that desc2 phase is linked back up to static_analysis. This means that at least one variable must be matched in the two places. There must be a nonlinear solver for this step.

For GASP missions, there is a takeoff subsystem (including taxi phase) within pre_mission. For FLOPS missions, there is a takeoff subsystem within pre_mission in which a takeoff phase is added. It is added after static_analysis in add_pre_mission_systems method.

Similarly, there is a landing phase within post_mission.

File opt_report.html#

This file is OpenMDAO Optimization Report. All values are in unscaled, physical units. On the top is a summary of the optimization, followed by the objective, design variables, constraints, and optimizer settings. Here is a screenshot:

Sample opt_report.html (top lines)

This file is important when dissecting optimal results produced by Aviary.

Coloring files#

There is a sub-folder coloring_files. Those are used internally by OpenMDAO for derivative computation. Users should generally skip those files.

File traj_linkage_report.html#

This is a dymos linkage report in a customized N2 diagram. It provides a report detailing how phases are linked together via constraint or connection. It can be used to identify errant linkages between fixed quantities.

Sample traj_linkage_report.html (top of timeseries)

We clearly see how those mission phases are linked.

File traj_results_report.html#

Note

This is one of the most important files produced by Aviary. It will help you visualize and understand the optimal trajectory produced by Aviary.

This file contains timeseries and phase parameters in different tabs. For our aircraft_for_bench_GwGm run, they are: groundroll, rotation, ascent, accel, climb1, climb2, cruise, desc1, and desc2 parameters. On the timeseries tab, users can select which phases to view. The following are the top of timeseries tab and ascent parameters tab:

Sample traj_results_report.html (top part of timeseries)

Let’s find the altitude chart. Move the cursor to the top of climb2. We see that the aircraft climbs to 37910 feet at 2814 second. Then it enters to cruise phase. At 28290 second, it starts descent from 37660 feet.

Sample traj_results_report.html (altitude)

Let’s switch to ascent tab. We see the following:

Sample traj_results_report.html (top lines) of ascent

This file is quite important. Users should play with it and try to grasp all possible features. For example, you can hover the mouse over the solution points to see solution value; you can save the interesting images; you can zoom into a particular region for details, etc.

Optimizer output#

If IPOPT is the optimizer, IPOPT.out is generated. If SLSQP is the optimizer and pyOptSparseDriver is the driver, SLSQP.out is generated. Generally speaking, IPOPT and SNOPT converge Aviary optimization problems better than SLSQP, but SLSQP is bundled with Scipy by default, making it more widely available.

SQLite database file#

There is a .db file after run. By default, it is aviary_history.db. This is an SQLite database file. In level 2 and level 3, we will be able to choose a different name. Our run is recorded into this file. You generally shouldn’t need to parse through this file on your own, but it is available if you’re seeking additional problem information.

Plain text inputs and outputs#

If Settings.VERBOSITY is set to 2 (Verbosity.VERBOSE) or higher, input_list.txt and output_list.txt are generated. We recommend this parameter be set to 1 (Verbosity.BRIEF) for beginners. See Coding Standards for more details.

Let us look at input_list.txt:

Sample input_list.txt (top lines)

In this screenshot, we see a tree structure. There are three columns. The left column is a list of variable names. The middle column is the value and the right column is the promoted variable name. pre_mission is a phase, static_analysis is a subgroup which contains other subcomponents (e.g. propulsion_group). engine is a component under which is a list of input variables (e.g. aircraft:engine:scaled_sls_thrust).

An input variable can appear under different phases and within different components. Note that its values can be different because its value has been updated during the computation. On the top-left corner is the total number of inputs. That number counts the duplicates because one variable can appear in different phases. Aviary variable structure are discussed in Understanding the Variable Hierarchy and Understanding the Variable Metadata.

If you zoom the N2 diagram into propulsion_group and gasp_based_geom, you see the exact same tree structure:

Sample N2 Diagram (top blocks)

File output_list.txt follows the same pattern. But there is another tree for implicit outputs. That helps when debugging is needed.

Additional messages#

When Aviary is run, some messages are printed on the command line and they are important. For example, the following constraint report tells us whether our desired constraints are met:

--- Constraint Report [traj] ---
    --- groundroll ---
        None
    --- rotation ---
        [final]   0.0000e+00 == normal_force [lbf]
    --- ascent ---
        [final]   5.0000e+02 == altitude [ft]
        [path]    0.0000e+00 <= load_factor <= 1.1000e+00  [unitless]
        [path]    0.0000e+00 <= fuselage_pitch <= 1.5000e+01  [deg]
    --- accel ---
        [final]   2.5000e+02 == EAS [kn]
    --- climb1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- climb2 ---
        [final]   3.7500e+04 == altitude [ft]
        [final]   1.0000e-01 <= altitude_rate  [ft/min]
        [final]   8.0000e-01 == mach [unitless]
    --- cruise ---
        None
    --- desc1 ---
        [final]   1.0000e+04 == altitude [ft]
    --- desc2 ---
        [final]   1.0000e+03 == altitude [ft]

Beyond level 1#

Aviary Level 1 is quite capable yet limited in how much is exposed to the user. If you need additional flexibility or want to add your own subsystems, you will need to move to Levels 2 or 3. Let us briefly discuss some of the possible reasons for you to move to Levels 2 or 3.

If you examine the Aviary metadata, you find that Aviary has two choices of objectives: Mission.Objectives.FUEL and Mission.Objectives.RANGE. But in level 1, you are limited to minimizing fuel. Users will be able to make other choices in level 2 and level 3.

Level 1 uses the default phase_info dictionaries, though there is an option to use a customized phase_info. For more details, please read drawing and running simple missions.

There may be more or fewer key-value pairs in each of the phases. For example, if you have an external subsystem, say, a battery with BatteryBuilder(), you can add a key external_subsystems with value BatteryBuilder(). We will go into details later.

No matter what mission method users choose, they can create their own phase sequences and run Aviary using either level 2 or 3 interface.

In the aviary/interface/default_phase_info subdirectory, you see another phase_info file solved.py. It defines the phase_info object for the “solved” mission method which is not explained in level 1. The “solved” mission method is a more simplified and robust method than the other methods. In just a few words, it’s a way of implementing a trajectory where the aircraft only flies in straight lines in the altitude-range space and usually at constant velocity. The idea is to simplify the design space by reducing the flexibility the optimizer has to choose a trajectory profile. It still uses a collocation method to solve the mission dynamics, but the complexity of the design problem is lessened.

We are ready to move on to Level 2.