Development of Aircraft Mission Performance Analysis Program
 Author: Lee Hyunseok, Lee Hyungjoon, Kwak Einkeun, Lee Seungsoo, Bae Seungho
 Organization: Lee Hyunseok; Lee Hyungjoon; Kwak Einkeun; Lee Seungsoo; Bae Seungho
 Publish: International Journal Aeronautical and Space Sciences Volume 14, Issue2, p162~171, 30 June 2013

ABSTRACT
A general purpose aircraft mission performance analysis program has been developed. The program can be used in design mode or in analysis mode. Fuel weight for a given mission profile can be estimated when the design mode is chosen, while mission time or mission range for a given fuel can be estimated when the analysis mode is chosen. The mission analysis program is written with JAVA and includes GUI(Graphic User Interface) for users’ conveniences. With a proper combination of databases for propulsion, aerodynamics and weight, the program can be configured to compute the performance of any type of aircraft. The program is validated by comparing its results with the results of a well known performance analysis program by ADD(Agency for Defense Development)

KEYWORD
Aircraft , Mission Profile , Mission Performance

Nomenclature
CD Drag coefficient
CL Lift coefficient
CLmax Maximum lift coefficient
M Mach number
p Pressure
Pset Power setting
T Net thrust
V Airspeed
W Aircraft weight
a Speed of sound
g Gravitational acceleration
x Range
△CD Delta drag coefficient
α Angle of attack
ρ Density
μ Viscosity
θ Pitch angle
1. Introduction
As all other machines, inflight performance of an airplane is one of the key factors to determine whether the airplane is worth buying over other airplanes. Many aspects of inflight performance can be addressed. While point performance deals with the airplanes ability to fly at a specific point in flight, mission performance deals with the overall performance throughout the entire flight. With the mission performance analysis, engineers can determine the fuel required for a given flight mission with a specific payload. The mission performance analysis, therefore, is necessary to compute not only the design fuel weight for aircraft design, but also the fuel weight for a commercial flight. The mission performance analysis, therefore, must be done by aircraft manufacturers as well as aviation companies.
Mission performance analysis programs have been developed by research institutes, and aircraft developers. Program NSEG is a rapid mission analysis code based on the use of approximate flight path equations of motion[1]. It was originally developed by U.S. Air Force Aeronautical System Division of WrightPatterson AirForce Base. Equations vary with the segment types; accelerations, climbs, cruises, descents, and deceleration. DARcorperation developed APP(Aircraft Performance Program) for aircraft performance analysis which is used as a performance module for the company’s aircraft design program, AAA(Advanced Aircraft Analysis)[2]. RDSprofessional and RDSstudent were developed and maintained by Conceptual Research Corporation[3]. They also have a common mission performance analysis module inside. The inhouse aircraft synthesis tool developed by Lockheed Martin is called the Adaptable Design Synthesis Tool(ADST)[4]. In this code, the aircraft design is perturbed by changing the wing geometric parameters from the baseline, and the aircraft is then fuel sized to meet the design mission. ADST equips with a performance module called MAPS(Mission Analysis and Performance System).
Most of the performance analysis programs and the design codes, however, were written with high level programming languages, which means they are highly machinedependent. A mission performance analysis program called IMAP(Interactive Mission Performance Analysis Program) is developed with JAVA programming language. Java is a generalpurpose, concurrent, classbased, objectoriented computer programming language that is specifically designed such that WORA(write once, run anywhere) is possible. WORA means that a code which operates on one platform does not need to be recompiled to run on another. JAVA programs are complied to bytecode that can run on any JAVA virtual machine regardless of computer architecture[5]. GUI(Graphic User Interface) is developed along with the program, so that users can use IMAP easily. Furthermore, IMAP is designed to be generalpurpose. With a proper combination of aerodynamics, propulsion, and weight databases, it can be used from lowspeed propeller driven aircraft to highspeed jet propelled aircraft.
The objective of the paper is to introduce the outline of IMAP, and to show its validity and its usefulness. To this end, we will first discuss the structure of IMAP, the databases needed to run the program, namely atmospheric database, aerodynamic database, weight database, and propulsion database. We will describe equations of motion for each mission segment and the numerical methods used in the program. Next, we will present three examples to demonstrate the validity and the usefulness of IMAP. Finally, we will summarize our work, and describe future research plans.
2. IMAP
As stated in the introduction, IMAP is written in JAVA so that it can run on any operating system. User interface with JAVA graphic library is also included in IMAP for user conveniences. One of features IMAP has is that users have freedom to generate a mission profile by sequencing mission segments in order. The mission segments can be specified as outbound(increasing range), or inbound(decreasing range). Also, the users can designate some of the segments as reserve. During the reserve segments, the aircraft weight decreases while the mission range or the mission time does not increase.
IMAP can be used either design mode or analysis mode. With the design mode, the mission fuel weight can be computed for a given mission profile. With the analysis mode, on the other hand, the users are allowed to specify one or two mission parameters as unknowns. The unknown parameters are determined in order to use the fuel completely at the end of the mission. For example, the cruise range or the loiter time can be set as an unknown for a given fuel. If both the cruise range and the loiter time are to be determined, the user should specify the time or the position of the aircraft at the end of the mission.
Figure 1 depicts the structure of IMAP. The users use GUI to input the mission profile, and to construct databases for an aircraft. The input data is then converted to SI unit via unit converter, and stored for later use. The solver consists of trim routine, optimization routine, segment analysis routine, and root finding routine. The databases that IMAP uses are atmospheric database, aerodynamic database, weight database, and propulsion database. The database can be generated by IMAP, text editors, or spreadsheet programs. Once the databases are set up, IMAP loads the database from files whose names are specified by the users.
2.1. Equations of Motion
Figure 2 depicts the free body diagram of an aircraft. It is
assumed that the pitching moment acting on the aircraft is trimmed with horizontal trim devices. It is equivalent to assume that the aircraft is modeled as a pointmass for the performance evaluation. As a consequence, the user should use the trimmed aerodynamic data. The equations of motion for the aircraft in the stability axis can be found,
where
γ ,R , and？_{T} are the flight path angle, the radius of the flight path arc, and the engine installation angle with respect to the longitudinal body axis, respectively. Equation (1) is supplemented with the kinetic relationswhere
x andh are the range and the altitude of the aircraft.2.2. Trim Method
Forces acting on the aircraft must be trimmed for all steady flights. The trim condition can be determined from the following equations, which are derived from equations (1) and (2)
where
n andP_{set} are load factor and power setting, respectively. Depending on the type of the segment, two parameters are computed with NewtonRaphson method. For a constant Mach number and constant altitude cruise, for example, theangle of attack and the power setting are the parameters to be determined from Eq. (4) and (5) with
n =1, andγ =0. The unknown vector,X can be computed iteratively by using the NewtonRaphson method,where
X ,F , andJ are given byThe derivatives are computed using the Newton quotient,
where
ε is a small parameter. In this paper,ε =10^{5} is used for all parameters. The value ofε is determined numerically such that it does not change the converged trimmed condition.2.3 Time Integration Method
From Eq. (1), (2), the equations of motion for the acceleration segment can be found,
From Eq. (10), the angle of attack can be found with a full throttle. The Euler explicit time integration method, then, is used to integrate Eq. (9) to obtain the time histories.
2.4. Optimization Methods
BFGS(BroydenFletcherGoldfarbShanno) method[6] is adopted for segment optimization. The BFGS method belongs to a class of quasiNewton methods. In the quasiNewton methods, the Hessian matrix of second derivatives does not need to be evaluated directly. Instead, the Hessian matrix is approximated using rankone updates specified by gradient evaluations. The quasiNewton methods are a generalization of the secant method to find the root of the first derivative for multidimensional problems. The BFGS method is one of the most popular methods for nonconstrained optimization problems. Fibonacci search[7] is used as a onedimensional search algorithm in conjunction with the BFGS method.
2.5. Mission Profile and Mission Segments
Mission profile is defined as an idealized diagram that describes the mission of the aircraft from takeoff to landing. The purpose of the mission performance analysis is to determine the fuel required to complete the mission. The mission profile, therefore, is divided into number of mission segments, or mission legs. The mission segment can be modeled in terms of the states and control variables as function of time. Table 1 summarizes the types of mission segments with which the users can model the mission profile. Fuel consumption during landing is ignored because it is small compared to other flight phases.
2.5.1. WUTO(Warmup and Takeoff)
The fuel burned during the warmup is computed from the power setting and the warmup time that the users specify. The takeoff is modeled as aircraft accelerating to a takeoff speed at a fixed angle of attack. The aero database for the WUTO segment should include the ground effects as well as landing gears and flaps.
2.5.2. Climb
The minimum time climb schedule gives a path defined by altitude and Mach number for transitioning from one level of energy to a higher level of energy in the shortest time. The time required to climb can be written as
where the specific excess power and the energy height are defined by
In order to climb in the shortest time, therefore, we can choose the climb condition where the specific excess power is the maximum at a constant energy height.
The weight change during the climb can be written as
where
？_{f} is the fuel flow rate. Thus, the minimum fuel climb can be achieved by following the route to maximizeP_{s} /？_{f} at the constant energy height.If a climb segment is followed by a cruise segment or a loiter segment, the user can couple the climb segment with the next cruise segment or loiter segment. The final altitude of the cruise segment is determined iteratively such that the next cruise segment or loiter segment is optimal in terms of the range or the loiter time.
2.5.3. Cruise and Loiter
The distance during the cruise can be written as a function of weight as
Here, it is assumed that the wind speed is zero. To maximize the range, therefore, we need to follow the flight path along which
V /？_{f} is maximal.The flight condition for the optimum loiter time can be computed by maximizing 1/
？_{f} , which can be seen from the following equation2.6. Databases
To use IMAP, the users should supply 4 databases, namely the atmospheric database, the aerodynamic database, the weight database, and the propulsion database. Figure 3 describes the details of the databases. The atmospheric properties are computed using ISA(International Standard Atmosphere) model. According to ISA model, density, speed of sound, viscosity and pressure are functions of the altitude and the deviation temperature from ISA condition, that is
where △
T is the deviation from ISA condition.As stated in the previous section, IMAP uses the trimmed aerodynamic data, so that the aero database contains the lift and the drag coefficients as functions of static margin, and configurations of the aircraft. The aero database also includes the delta drag due to external stores, the maximum lift coefficient. The basic aero database is given in the form of
where
SM ,I_{GE} ,δ_{F} , andI_{G} are the static margin, the ground effect flag, the flap angle, and the landing gear flag, respectively. Thedelta drag due to stores is a function of the store index,
m and Mach number. The maximum lift coefficient is a function of Mach number, Reynolds number, the ground effect, and the flap angle,The net thrust and the fuel flow rate for a jet engine are provided as following,
The shaft horse power, and the fuel flow rate for a propeller driven aircraft are given by
The weight database consists of the static margin, the store weight,
W_{s} and the c.g. movement, △_{cg} due to store,3. Numerical Examples
To demonstrate the validity and the usefulness of IMAP, three exemplary mission analyses involving two types of aircraft are presented. The results of IMAP are compared with the results of a wellknown mission performance analysis program, AMAP(ADD Mission Performance Analysis Program)
[8] for comparison and validation. Examples are analyzed with both IMAP and AMAP using the analysis mode. In addition, IMAP is run with the initial and final conditions of all segments taken from the results of AMAP. Therefore, the computational methods used for mission segment analysis can be validated. The result with the analysis mode can validate IMAP for integrating segments.
3.1. Jet Powered Aircraft
The first example is a HiLoLoHi strike mission of a jet powered aircraft. The mission consists of a high altitude cruise, low altitude dashin, weight drop, low altitude dashout, and high altitude cruise. Here, the total mission radius is
an unknown, and is to be determined from the total fuel. The details of the segments are presented in Fig. 4. The detailed descriptions of the aircraft including the weight of fuel tanks are listed in Table 2. The drag decrease and the weight decrease due to the tank release are included in the computation. Figures 5 and 6 present the two analysis results of IMAP. Figure 5 presents the results in terms of range and altitude, while Fig. 6 the results in terms of range and weight (wedge plot). The result of AMAP is also presented for comparison. As stated earlier, the initial and final conditions for the validation case are taken from s of AMAP. Three results match well each other, indicating IMAP is accurate.
The mission profile of the second example is depicted in Fig. 7.
The mission comprises an outbound cruise, weight drop, and inbound cruise besides takeoff and landing (HiHiHi Attack). The same aircraft as the first example is used for the second example. The mission radius is going to be determined from the total fuel. Figure 8 and 9 summarize the mission analysis results. Figures show the validation result matches well with the results of AMAP, which indicates that the numerical methods used to solve the segments are correct. The mission range obtained with IMAP, however, is wider than that with AMAP. In addition, the wedge plot obtained with IMAP differs from that with AMAP. The reason for this discrepancy can be understood from Fig. 10. Figure 10 presents the optimization results of four segments for the mission profile with two programs.
P_{s} /？_{f} must be maximal along the climb schedule. For the outbound and inbound cruise segments, on the other hand,
V /？_{f} should be maximal. For the constant altitude loiter, the flight Mach number should be chosen such that 1/？_{f} should be maximal. Figure 10 shows that during 4 segments, IMAP does better job in finding the optimal point for all segments compared to AMAP.3.2. Propeller Driven Aircraft
The mission profile of the final example is presented in Fig. 11. The aircraft used in this example is a propeller driven aircraft.
The mission starts with WUTO, is followed by a min time climb, fixed range outbound cruise, loiter, fixed range inbound cruise, and landing. In this example, the loiter time is to be determined from the given fuel. Furthermore, the final altitude of the min time climb segment is chosen such that the next cruise altitude is optimal. The specifications of the aircraft are summarized in Table 3. As before, the validation result is also presented for comparison. Unlike IMAP and AMAP cases, the initial climb segment for the validation is not linked to the next cruise segment. Moreover, the endurance for the validation is set with the result of AMAP.
Figure 12 and 13 present the analysis results with IMAP and AMAP. The results with IMAP show the discrepancy with those of AMAP as well as the validation. Also, the difference occurred during the climb is maintained for the rest of the mission. IMAP finds the optimal cruise altitude of the outbound cruise segment lower than AMAP does. According the results of the validation case, the aircraft lowers its altitude after the min time climb in order to meet the optimal condition of the next cruise segment. The difference in the final altitude of the initial climb persists throughout the rest of the mission. In Fig. 14, the optimization results of four segments in regard to the mission are presented. The specific excess power with IMAP is higher than that with AMAP. As a result of higher climb rate, the distance for the climb segment is shorter. The difference in this segment persists over
the rest of the mission. Table 4 summarizes the elapsed times calculated from two programs during the climb segment. The table also confirms that IMAP finds the shortest climb schedule. The value of 1/
？_{f} calculated with IMAP is higher than that with AMAP during the loiter segment, which indicates IMAP finds the optimal loiter condition. Table 5 presents the loiter times obtained from two programs. IMAP, however, fails to find the optimal cruise condition from the inbound cruise. However, the difference is small compared to the rest of the mission. As inthe second examples, the differences in the results of IMAP and AMAP for the third example come from the optimization results. Therefore, it is safe to say that the differences in the results of the analysis come from the differences in optimization methods for the two programs.
4. Summary and Concluding Remarks
In this paper, an aircraft mission performance analysis program called IMAP was developed. IMAP was written in JAVA so that it is machineindependent. Also, a graphic user interface was developed with IMAP for users’ conveniences. IMAP can estimate fuel required for a given mission profile, or compute mission time or mission radius for a given fuel. Three exemplary mission analyses were chosen to demonstrate IMAP’s accuracy, efficiency and usefulness. IMAP is now being upgraded to include the point performance analysis capability. It is intended to further expand IMAP’s capability to include the
aircraft synthesis in near future.

2.

3.

5.

[Fig. 1.] Structure of IMAP

[Fig. 2.] Freebody Diagram for Aircraft Performance Evaluation

[Table 1.] Mission segment types

[Fig. 3.] Databases for IMAP

[Table 2.] Specifications of Jet Powered Aircraft

[Fig. 4.] Mission Profile of Sample 1

[Fig. 5.] Analysis Results of Sample 1

[Fig. 6.] Weight Variation of Sample 1

[Fig. 7.] Mission Profile of Sample 2

[Fig. 8.] Analysis Results of Sample 2

[Fig. 9.] Weight Variation of Sample 2

[Fig. 10.] Optimization Results

[Fig. 11.] Mission Profile of Sample 3

[Table 3.] Specifications of Propeller Driven Aircraft

[Fig. 12.] Analysis Results of Sample 3

[Fig. 14.] Optimization Results

[Fig. 13.] Weight Variation of Sample 3

[Table 4.] Times Elapsed during Min Time Climb

[Table 5.] Optimum Loiter Times