Benefit is claimed under 35 U.S.C. 119(a)-(d) to Foreign application Serial No. 1478/CHE/2011 filed in INDIA entitled “SYSTEM AND METHOD FOR AIRCRAFT PERFORMANCE PREDICTIONS FOR DESCENT AND APPROACH PHASES” by AIRBUS ENGINEERING CENTRE INDIA, filed on Apr. 29, 2011, which is herein incorporated in its entirety by reference for all purposes.
Embodiments of the present subject matter relate to flight management systems. More particularly, embodiments of the present subject matter relate to aircraft performance predictions for descent and approach phases.
Modern aircrafts generally include an onboard computing system called flight management system (FMS). The FMS provides the flight crew with information on the future states of the aircraft along a selected flight path. This information includes predictions based on the performance of the aircraft during various phases of the flight. Typically, the flight plan includes a descent phase, which normally starts from the end of cruise (also referred to as top-of-descent distance) and is considered to be completed once the aircraft reaches a runway threshold through an approach phase, which is mainly performed at idle thrust. Predictions of time, fuel consumption, flying altitude, aircraft speed and other parameters at en route points until the runway threshold provide useful information about the future performance of the aircraft. Typically, such predictions are displayed to the flight crew and used in the guidance of the aircraft along the vertical trajectory.
Currently, the criteria used for obtaining the information are preset or the parameter values may be assumed in obtaining the information. For example, climb angle, fuel flow and true air speed (TAS) may be considered to vary linearly with the altitude. Further, constant values may be considered for parameters, such as air speed, lapse rate, temperature and the like to obtain the information. This can result in providing inaccurate information to the flight crew during the descent phase.
A system and method for aircraft performance predictions for descent and approach phases are disclosed. According to one aspect of the present subject matter, a method of aircraft performance predictions for descent and approach phases in a flight management system (FMS) includes stopping cruise computation substantially around a default distance from a destination. Further, current predicted aircraft state is determined using a total energy of the aircraft starting from the default distance. Furthermore, descent and approach phase profiles are computed using the determined current predicted aircraft state. In addition, the aircraft performance predictions are obtained using the computed descent and approach phase profiles.
According to another aspect of the present subject matter, a non-transitory computer-readable storage medium for the aircraft performance predictions for the descent and approach phases, having instructions that, when executed by a computing device causes the computing device to perform the method described above.
According to yet another aspect of the present subject matter, the aircraft includes the FMS. Further, the FMS includes a processor and memory coupled to the processor. Furthermore, the memory includes an aircraft performance predictions module. In one embodiment, the aircraft performance predictions module includes instructions to stop cruise computation substantially around the default distance from the destination. Moreover, the aircraft performance predictions module determines the current predicted aircraft state using the total energy of the aircraft starting from the default distance. In addition, the aircraft performance predictions module computes the descent and approach phase profiles using the determined current predicted aircraft state. Also, the aircraft performance predictions module obtains the aircraft performance predictions using the computed descent and approach phase profiles.
The methods and systems disclosed herein may be implemented in any means for achieving various aspects. Other features will be apparent from the accompanying drawings and from the detailed description that follow.
Various embodiments are described herein with reference to the drawings, wherein:
The drawings described herein are for illustration purposes only and are not intended to limit the scope of the present disclosure in any way.
A system and method for aircraft performance predictions for descent and approach phases are disclosed. In the following detailed description of the embodiments of the present subject matter, reference is made to the accompanying drawings that form a part hereof, and in which are shown by way of illustration specific embodiments in which the present subject matter may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the present subject matter, and it is to be understood that other embodiments may be utilized and that changes may be made without departing from the scope of the present subject matter. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present subject matter is defined by the appended claims.
At block 106, descent and approach phase profiles are computed using the determined current predicted aircraft state. This is explained in more detail with reference to
At block 108, the aircraft performance predictions are obtained using the computed descent and approach phase profiles. In one embodiment, a first top-of-descent distance is obtained from the obtained descent and approach phase profiles. Further, it is determined whether a difference between the starting top-of-descent distance and the first top-of-descent distance is within a design value. If the difference between the starting top-of-descent distance and the first top-of-descent distance is within the design value, the aircraft performance predictions are determined, using the obtained descent and approach phase profiles, in a forward computation. Furthermore, the aircraft performance predictions are determined for each way point using the descent and approach phase profiles in the forward computation.
Also in this embodiment, if the difference between the starting top-of-descent distance and the first top-of-descent distance is not within the design value, the obtained first top-of-descent distance is considered as the starting top-of-descent distance. Further, the steps of obtaining the first top-of-descent distance from the obtained descent and approach phase profiles and determining whether the difference between the starting top-of-descent distance and the obtained first top-of-descent distance is within the design value are repeated until the difference is within the design value. This is explained in more detail with reference to
Referring now to
Typically, the approach phase starts from the point 212, shown in
As shown in
In this embodiment, the descent and approach phase profiles, shown in
In addition in this embodiment, the computation of the descent phase starts from the point 212, shown in
Also in this embodiment, the descent and approach phase profiles are computed using the current predicted aircraft state using the total energy of the aircraft. This is explained in more detail with reference to
Referring now to
In an exemplary first iteration, the starting top-of-descent distance, point 308, as shown in
Further in this embodiment, if the difference between the starting top-of-descent distance 308 and the first top-of-descent distance 310 is within the design value, the aircraft performance predictions are determined using the obtained descent and approach phase profiles associated with the first top-of-descent distance 310. Furthermore, if the difference between the starting top-of-descent distance 308 and the first top-of-descent distance 310 is not within the design value, the iteration is repeated.
In an exemplary second iteration, the first top-of-descent distance 310, obtained from the first iteration, is considered as the starting top-of-descent distance. Further, the backward computation is repeated to obtain another first top-of-descent distance, say point 312, as shown in
Referring now to
At block 406, the starting top-of-descent distance is selected while the aircraft is in the cruise altitude. At block 408, the descent and approach phase profiles are computed from the destination by performing the iterative backward computation. This is explained in more detail with reference to
At block 414, a check is made to determine whether the computed difference is within the design value. If the difference is not within the design value, at block 416, the obtained first top-of-descent distance is considered as the starting top-of-descent distance. Furthermore, the process flow 400 performs the steps in block 408. If the difference is within the design value, at block 418, the obtained first top-of-descent distance is considered as the top-of-descent distance for the flight plan
Referring now to
Generally, the total energy gained by the aircraft 500 when excess force is acting on the aircraft 500 is expressed using equation:
ΔETOTAL=ΔEPE+ΔEKE (1)
wherein,
ΔETOTAL is a total energy gained by the aircraft 500;
ΔEPE is a potential energy gained by the aircraft 500; and
ΔEKE is a kinetic energy gained by the aircraft 500.
The excess force available in the aircraft 500 is expressed using equation:
f
excess
x=mgh+½mv2 (2)
wherein,
fexcess is an excess force available in the aircraft 500;
x is a distance covered by the aircraft 500;
m is a mass of the aircraft 500;
g is acceleration due to gravity;
h is altitude of the aircraft 500; and
v is a true air speed (TAS) 510, shown in
Further, differentiating the equation (2) results in the following equation:
f
excess
dx/dt=mg(dh/dt)+½m(dv2/dt) (3)
Applying, dv2=2vdv and dividing both sides of equation (3) by the term ‘m’, we obtain the following equation:
f
excess
v/m=g(dh/dt)+v(dv/dt) (4)
wherein,
dh/dt is a rate of climb/descent of the aircraft 500; and
dv/dt is an acceleration of the aircraft 500.
Typically, force is a combination of a thrust 506, shown in
excess specific energy=(thrust−drag)*TAS/m (5)
In the equation (5),
(thrust−drag)=fexcess (6)
As shown in
drag=½ρS(TAS)2CD (7)
wherein,
ρ is a density of air;
S is a surface area of the wing of the aircraft 500; and
CD is a drag coefficient.
In this embodiment, the equation (4) defined above is used for computing the descent and approach phase profiles as explained below. The steps involved in computing the approach phase includes computation in the constant FPA CAS segment 222, the constant FPA segment 218 and the second decel segment 214 shown in
dh/dt=TAS*sin γ≈v*γ (8)
wherein, γ, as shown in
Further, a change in TAS (dv) is obtained using equation:
dv=TAS
n
−TAS
n-1 (9)
wherein,
n is a time step (e.g., n=0, 1, 2 . . . ).
Furthermore, the fuel flow is obtained by computing the thrust 506 for the constant FPA CAS segment 222. Using the equation (4), thrust 506 can be derived using equation:
drag−thrust=mg(dh/dt)/v+m(dv/dt)
The thrust 506 is computed in the constant FPA CAS segment 222 and in all other segments, such as the constant FPA segment 218, the decel segment 214, the constant CAS descent segment 210 and the constant MACH descent segment 206 the thrust 506 remains idle.
Since dh/dt=v*γ, we can obtain thrust 506 using equation:
thrust=drag−(mgγ+m(dv/dt))
In addition, drag 502 is computed using the equation (7).
The computation in the constant FPA segment 218, shown in
dv/dt=1/v*(fexcessv/m−g(dh/dt))
Further, dh/dt can be obtained using the equation (8).
The computation in the second decel segment 214, shown in
(dv/dt)=fexcess/m
The rate of climb for the second decel segment 214 is equal to zero, which is expressed as:
dh/dt=0
Further in this embodiment, the computation in the descent phase includes the computations in the constant CAS descent segment 210 and the constant MACH descent segment 206. Generally, in the descent phase, dv/dt is a small value and the CAS/MACH of the aircraft 500 is held constant. The steps involved in computing the descent phase are explained below.
In an exemplary first step, the speed change (dv/dt) is computed by subtracting the previous time step TAS 510 from the current time step TAS 510. Using the equation (9), dv can be derived as shown below:
dv=TASn−TASn−1
wherein,
n is a time step (e.g., n=0, 1, 2 . . . ).
Further in an exemplary second step, the kinetic energy (KE) change (vdv/dt) required is computed. In the constant CAS descent segment 210 and the constant MACH descent segment 206 the MACH of the aircraft 500 is held constant. However, for a given CAS or MACH the TAS 510, shown in
TAS=√(ρ0/ρ)*Ka*CAS
wherein,
ρ0 is an air density at mean sea level;
ρ is a air density at the given altitude; and
Ka is a compressibility correction.
Since the TAS 510 of the aircraft 500 increases with the altitude of the aircraft 500, some portion of the total energy (ΔETOTAL) is used to maintain the required speed. Therefore, a portion of the ΔETOTAL is converted to KE. The KE is represented using equation:
KE=vdv/dt
Furthermore in an exemplary third step, the remaining energy which is used for potential energy (PE) change is computed by subtracting the value of KE required, obtained in the second step, from fexcess v/m obtained using equation (4). Therefore, the energy available for PE change is expressed using equation:
PEchange=fexcessv/m−v(dv/dt)
In addition in an exemplary fourth step, the rate of climb (dh/dt) is computed using equation:
dh/dt=(fexcessv/m−v(dv/dt))/g
Using the above described computations the descent and approach phase profiles for the aircraft 500 are obtained. Further, the top-of-descent distance is obtained using the obtained descent and approach phase profiles. This is explained in more detail with reference to
In this embodiment, the aircraft performance predictions for the next time step are computed as shown below. Exemplary aircraft performance predictions includes flight parameters, such as flying altitude, aircraft speed, distance covered by the aircraft 500, time, fuel consumed and the aircraft gross weight. In this embodiment, the altitude of the aircraft 500 for the next time step is computed using equation:
ALT
n
=ALT
n-1+∫(dh/dt)dt (10)
wherein,
n is a time step (e.g., n=0, 1, 2, . . . );
ALTn is an altitude at the nth time step; and
ALTn-1 is an altitude at the (n−1)th time step.
Further in this embodiment, the TAS 510 of the aircraft 500 for the next time step is computed as shown below. In the constant CAS descent segment 210, shown in
V
TASn
=V
TASn-1+∫(dv/dt)dt
wherein,
n is a time step (e.g., n=0, 1, 2, . . . );
VTASn is a TAS at the nth time step; and
VTASn-1 is a TAS at the (n−1)th time step.
Also, TAS 510 is computed in the constant CAS descent segment 210 and the constant MACH descent segment 206 for a given altitude, CAS and delta international standard atmosphere (DISA). The altitude is obtained using the equation (10) show above. The CAS of the aircraft 500 is held constant. The DISA is the temperature difference with respect to international standard atmosphere which is obtained by computing a difference between the actual temperature and an ISA temperature. The ISA temperature is computed using equation:
ISA
temperature
=T
0−1.98*[ALT(feet)/1000]
wherein,
T0 is a temperature at mean sea level.
Further, the MACH number is computed using standard ISA equations and the TAS 510 is computed using the equation:
V
TAS
=M*a
wherein,
M is a Mach number; and
a is a speed of sound at the flying altitude.
Furthermore in this embodiment, the gross weight for the next time step is obtained using equation:
GW
n
=GW
n-1−∫fuelflowrate*dt
wherein,
n is a time step (e.g., n=0, 1, 2 . . . );
GWn is a gross weight at nth time step; and
GWn−1 is a gross weight at (n−1)th time step.
In addition, fuel flow rate is computed using engine data stored in the performance database.
Moreover in this embodiment, the distance for the next time step is obtained using equation:
DISTn=DISTn−1+∫VGS*dt
wherein,
n is a time step (e.g., n=0, 1, 2 . . . );
DISTn is a distance at nth time step; and
DISTn−1 is a distance at (n−1)th time step.
In addition, VGS is obtained using the equation:
V
GS
=V
TAS+Windspeed
Also, time for the next time step is obtained using equation:
TIMEn=TIMEn-1+dt
wherein,
n is a time step (e.g., n=0, 1, 2 . . . );
TIMEn is a time at nth time step; and
TIMEn−1 is a time at (n−1)th time step.
The above described computations provide the aircraft performance predictions for each way point in the forward computation. The computations are same for both the backward computation and the forward computation. The aircraft performance predictions thus obtained are presented to the flight crew for information and also to the flight guidance computer for guidance through the descent and approach phases.
Now referring to
The FMS 602 includes a processor 604, memory 606, a removable storage 618, and a non-removable storage 620. The FMS 602 additionally includes a bus 614 and a network interface 616. As shown in
Exemplary user input devices 622 include a digitizer screen, a stylus, trackball, keyboard, keypad, mouse and the like. Exemplary output devices 624 include a display unit of the personal computer, the mobile device, the FMS, and the like. Exemplary communication connections 626 include a local area network, a wide area network, and/or other networks.
The memory 606 further includes volatile memory 608 and non-volatile memory 610. A variety of computer-readable storage media are stored in and accessed from the memory elements of the FMS 602, such as the volatile memory 608 and the non-volatile memory 610, the removable storage 618 and the non-removable storage 620. The memory elements include any suitable memory device(s) for storing data and machine-readable instructions, such as read only memory, random access memory, erasable programmable read only memory, electrically erasable programmable read only memory, hard drive, removable media drive for handling compact disks, digital video disks, diskettes, magnetic tape cartridges, memory cards, Memory Sticks™, and the like.
The processor 604, as used herein, means any type of computational circuit, such as, but not limited to, a microprocessor, a microcontroller, a complex instruction set computing microprocessor, a reduced instruction set computing microprocessor, a very long instruction word microprocessor, an explicitly parallel instruction computing microprocessor, a graphics processor, a digital signal processor, or any other type of processing circuit. The processor 604 also includes embedded controllers, such as generic or programmable logic devices or arrays, application specific integrated circuits, single-chip computers, smart cards, and the like.
Embodiments of the present subject matter may be implemented in conjunction with program modules, including functions, procedures, data structures, and application programs, for performing tasks, or defining abstract data types or low-level hardware contexts. Machine-readable instructions stored on any of the above-mentioned storage media may be executable by the processor 604 of the FMS 602. For example, a computer program 612 includes machine-readable instructions capable of providing aircraft performance predictions for descent and approach phase profiles in the FMS 602, according to the teachings and herein described embodiments of the present subject matter. In one embodiment, the computer program 612 is included on a compact disk-read only memory (CD-ROM) and loaded from the CD-ROM to a hard drive in the non-volatile memory 610. The machine-readable instructions cause the FMS 602 to encode according to the various embodiments of the present subject matter.
As shown, the computer program 612 includes an aircraft performance predictions module 628. For example, the aircraft performance predictions module 628 can be in the form of instructions stored on a non-transitory computer-readable storage medium. The non-transitory computer-readable storage medium having the instructions that, when executed by the FMS 602, causes the FMS 602 to perform the one or more methods described in
The systems and methods described herein enables determining the top-of-descent distance by performing an iterative backward computation from the destination to the cruise altitude. The iterative backward computation is performed based on the current predicted aircraft state without making any assumptions with respect to the current predicted aircraft state parameters. Further, the above mentioned embodiments enable determining aircraft performance predictions along the computed flight path without making any simplifying assumptions to the flight path. In addition, the gross weight and the time predictions are computed using a forward computation after determining the top-of-descent distance.
Although the present embodiments have been described with reference to specific example embodiments, it will be evident that various modifications and changes may be made to these embodiments without departing from the broader spirit and scope of the various embodiments. Furthermore, the various devices, modules, analyzers, generators, and the like described herein may be enabled and operated using hardware circuitry, for example, complementary metal oxide semiconductor based logic circuitry, firmware, software and/or any combination of hardware, firmware, and/or software embodied in a machine readable medium. For example, the various electrical structure and methods may be embodied using transistors, logic gates, and electrical circuits, such as application specific integrated circuit.
Number | Date | Country | Kind |
---|---|---|---|
1478/CHE/2011 | Apr 2011 | IN | national |