Nonlinear Adaptive Velocity Controller Design for an Air-breathing Supersonic Engine
Nonlinear Adaptive Velocity Controller Design for an Air-breathing Supersonic Engine
International Journal of Aeronautical and Space Sciences. 2012. Sep, 13(3): 361-368
Copyright ©2012, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( 3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : April 18, 2012
  • Accepted : August 24, 2012
  • Published : September 30, 2012
Export by style
Cited by
About the Authors
Jung-Woo, Park
Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
Ik-Soo, Park
Agency for Defense Development, Daejeon, Republic of Korea
Min-Jea, Tahk
Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea

This paper presents an approach on the design of a nonlinear controller to track a reference velocity for an air-breathing supersonic vehicle. The nonlinear control scheme involves an adaptation of propulsive and aerodynamic characteristics in the equations of motion. In this paper, the coefficients of given thrust and drag functions are estimated and they are used to approximate the equations of motion under varying flight conditions. The form of the function of propulsive thrust is extracted from a thrust database which is given by preliminary engine input/output performance analysis. The aerodynamic drag is approximated as a function of angle of attack and fin deflection. The nonlinear controller, designed by using the approximated nonlinear control model equations, provides engine fuel supply command to follow the desired velocity varying with time. On the other hand, the stabilization of altitude, separated from the velocity control scheme, is done by a classical altitude hold autopilot design. Finally, several simulations are performed in order to demonstrate the relevance of the controller design regarding the vehicle.
1. Introduction
The ramjet-powered engine is regarded as one of the high efficient solutions for high speed and long flight distance from its first appearance during the 1900s [1] . However, due to its strong coupling effects between its flight condition and combustion mechanism, the study on the design of flight guidance/control of the engine has been limited by the complexity of propulsive characteristics affected by varying flight conditions.
The ramjet engine is an air-breathing system which uses incoming air-flow as a combustion oxidant, where the airflow is compressed through the shock generation mechanism with its intake shape (called ram-compression). Related with the nature of the ram-compression, it is known that the performance of the engine depends highly on not only its combustor conditions by fuel and exhaust nozzle actuations but the vehicle’s maneuvering during flight [2 - 4] .
In spite of the difficulties in addressing the complex propulsion mechanism in vehicle model dynamics, there have been pioneer researches to design a flight controller of a hypersonic vehicle [5 , 6] . They assumed that the thrust and aerodynamic effects follow designated parametric equations roughly defined with uncertainties. With the approximated control model equations, they used the adaptive sliding mode controller design method to achieve their objectives to control thrust (or velocity), altitude, etc. Recently, the controller design approach by the nonlinear adaptive method has been popularly used to design a controller of a hypersonic vehicle [7 , 8] with more accurately developed nonlinear dynamic models for an air-breathing hypersonic vehicle as in Ref. [9 , 10] . Although the control methods present the design processes with the dynamic models, which capture the characteristics of hypersonic vehicles (i.e. scramjet), the generalized approach for aerodynamic and propulsive properties approximated in equation forms is applicable in the work for a supersonic vehicle (i.e. ramjet).
In our previous study in Ref. [11] , it is shown that the thrust of a considered ramjet engine is given from a tabular database in a function of Mach, altitude, angle of attack, equivalence ratio, and exhaust nozzle throat actuation, which are obtained by the engine model in Ref. [12] . With the thrust database, an approximated function for the thrust variation can be extracted where effecting variables and coefficients are varying with the flight condition for Mach and altitude. On the other hand, the aerodynamic properties of a supersonic vehicle can also be assumed to follow the functions regarding angle of attack and fin deflection angle for the control-design model (CDM) described in Ref. [8] since the aerodynamic characteristics of hypersonic vehicles are not much different from those of supersonic vehicles.
From the preceding discussions, a control model of a ramjet-powered vehicle can be obtained and it shows explicit dynamic relations between state variables and control inputs. In this paper, with explicit dynamic equations of a ramjet-powered vehicle, a nonlinear adaptive controller design is performed for a vehicle’s velocity to track a reference command, as the vehicle’s altitude is stabilized at a desired level by an altitude hold autopilot system. For the verification of the controller, several simulations are carried out with various types of velocity command functions.
Consequently, this paper is organized as follows. First, the generalized control model of a ramjet-powered vehicle based on previous studies [7 - 9] is defined where its thrust profile is approximated in the form of an equation through an engine performance analysis. Second, an overall procedure for the controller design of the vehicle is introduced, mainly disserting the nonlinear adaptive controller design for the velocity subsystem mainly. Finally, some discussion is followed after several simulations in order to demonstrate the controller.
2. Ramjet-powered Vehicle Dynamic Model
In this section, we derive a control-oriented version of the model for the longitudinal dynamics of a ramjet-powered vehicle dynamic model. The simplification is necessary to overcome the analytical intractability of the original model, which will be henceforth referred to as the true model. For the sake of simplicity, the scheme for a control-oriented model in Ref. [9] is employed in this study. The controloriented model, after the procedure of the simplicity, finally represents the explicit state space equation.
- 2.1 Equations of Motion
In general, the nonlinear longitudinal equations of motion are given as
Lager Image
This model is composed of five rigid-body state variables x=[V, h, γ, Q]T . In this paper, the propulsive and aerodynamic terms for T, D, L , and Q in Eq. (1) are originally available from the true model of the engine. The thrust T is evaluated by a higher-fidelity engine model (integral performance model) introduced in Ref. [12] . In addition, the aerodynamic forces are from the aerodynamic database based on the considered ramjet geometry in the present work. The propulsive and aerodynamic forces by the given engine and aerodynamic models represent the true ramjet-powered vehicle motions governed by Eq. (1).
On the other hand, in order to control the vehicle’s velocity, it is necessary to approximate the terms for T and D in regards to the internal loop control (or state) variables. Firstly, when we adopt the scheme for D in Ref. [9] , the following expression is obtained.
Lager Image
The approximation in Eq. (2) is linear in its fit parameters and the parameters are assumed to vary with flight conditions for Mach number M and air density ρ , noting that the simplified drag equation in Eq. (2) mimic the true aerodynamic drag profile. Also, in this paper, it is additionally assumed that the angle of attack α and fin deflection angle δe are measurable (or estimative) variables so that they can be known when constructing a control law in a later section. The approximation of the term for T is done in the following section with its tendency analysis.
- 2.2 Thrust Variation Tendency and its Approximation
As mentioned, the thrust of the ramjet engine is a function of five parameters (Mach number, altitude, angle of attack, fuel flow rate, and nozzle throat area). Here, the fuel flow rate is replaced by the equivalence ratio which is the actual air fuel ratio/the air fuel ratio for the stoichiometric mixture.
Lager Image
Maximum Thrust according to Change of Angle of Attack (a), Equivalence Ratio (b), and Exhaust Nozzle Throat (c)
Lager Image
Prior to initiating the investigation for the influence of the parameters on thrust, in this paper, the flight conditions for Mach number and altitude are excluded as focusing on the internal loop variables for α, Ɛq , and Ath . Figure 1 presents the results of the tendencies of the maximum thrust according to the considered variables and the fitted functions for each of them of which the details are given in Table 1 .
Approximated Thrust Functions
Lager Image
Approximated Thrust Functions
From the results in Fig. 1 and Table 1 , it is acceptable that the thrust can be defined as a function of α , Ɛq , and Ath as in Eq. (4).
Lager Image
In Eq. (4), there are two main control inputs for thrust, Ɛq and Ath . For the sake of simplicity in regards to the control structure, it is assumed that the exhaust nozzle throat actuation is done by a pre-designed operation schedule. The assumption is generally acceptable because the accurate actuation of the exhaust nozzle of the ramjet engine may be achieved with a more difficult process and the effectiveness of the nozzle actuation is less than those of its fuel actuation [2] . From the assumption for simplicity, Eq. (4) is rewritten as follows.
Lager Image
Lager Image
3. Controller Design
In this section, a procedure on the design of a velocity controller is mainly shown by using an adaptive nonlinear approach. Besides, we briefly mention the process of an altitude hold autopilot design.
The controller design starts from the decomposition of the equations of motion into functional subsystems, namely, the horizontal translation dynamics (the velocity subsystem) and the vertical translation dynamics (the altitude subsystem). Each subsystem is controlled separately using the available inputs at that level and intermediate virtual control command, as shown in Fig. 2 . In particular, as an extension from the preceding discussions, a control law with adaptive compensation is derived for the velocity subsystem by controlling thrust from the input of the fuel equivalence ratio, Ɛq . For the velocity subsystem, the goal pursued in this study is to design a dynamic controller of the form
Lager Image
using feedback from the state x only to steer the output of the system in Eq. (1) from a given set of initial values of velocity and altitude to the desired condition V * along reference trajectories Vref(t) . Here, the vector
Lager Image
indicates the vector of uncertain parameters. On the other hand, the altitude dynamics are controlled by means of the missile fin deflection angle actuation, δe , from a reference altitude, where the dynamic sequences are in charge of the altitude response.
Lager Image
Block Diagram of the Control Architecture
- 3.1 Adaptive Nonlinear Controller Design for the Velocity Subsystem
The control law for the velocity subsystem is derived by using robust adaptive dynamic inversion. By substituting the expressions of D and T from Eqs. (2) and (5) into the first equation of Eq.(1), the velocity error dynamics become
Lager Image
Lager Image
Define the vector of uncertain parameters as
Lager Image
The components in Eq. (10) indicate the whole parameters in Eqs. (2) and (6), and they are to be estimated with sequences in the control law generation algorithm. Then, the error dynamics in Eq. (8) can be rewritten with additional vectors for the simplicity of the equation as follows.
Lager Image
Lager Image
In the adaptive nonlinear controller scheme, we want to regulate the variable,
Lager Image
and the estimation error for the uncertain parameters,
Lager Image
To achieve this goal, we consider the Lyapunov candidate function below.
Lager Image
where k 1 is a positive control gain and Г is a diagonal positive definite matrix. The derivative of Eq. (13) induces the following relation.
Lager Image
Here, the term
Lager Image
on the right hand side of Eq.(14) is from a perception as
Lager Image
In order to satisfy the condition in Eq. (14), the updated law of the given uncertain parameter vector θ is finally given as follows.
Lager Image
The control law for the velocity stabilization is derived from Eq. (15) and the one feasible solution from the quadratic equation can be used as the equivalence ratio command. Here, the fuel supply actuation by this command is assumed to be the 1 st order dynamics, with a cutoff frequency of approximately 0.8 Hz [13] , which is equivalent to a time constant τeq of 0.2 sec, and it is realized by a transfer function.
Lager Image
In addition, in this paper, the diagonal positive definite gain matrix Г for the velocity subsystem is constructed after considering the order of each value regarding uncertain parameters. From the true model of the engine under the given flight conditions, the uncertain parameters can be given roughly, which are used for their initial guess. If the motion of the ramjet engine is trimmed at the flight condition of Mach 2.9 and an altitude of 14km, the initial estimation for the uncertain parameters are shown as
Lager Image
The values in Eq. (18) are approximately obtained by fitting the tendencies of the variations in thrust and drag for the independent variables (angle of attack, equivalence ratio, and fin deflection angle). For example, the sample variations in drag with the change of fin deflection at the trimmed value are fitted as the second order polynomial such that the values of
Lager Image
Lager Image
can be obtained.
Considering the order of the uncertainties, the components of the gain matrix Г are
Lager Image
- 3.2 Altitude Hold Autopilot Design
When a reference altitude command, href is given to be stabilized, it is achieved by using an altitude hold autopilot system for the missile’s fin deflection angle actuation, shown in Fig. 3 .
Lager Image
Control Sequence of Altitude Hold Autopilot
In the control sequence, the pitch rate feedback and altitude feedback are considered to generate fin deflection angle command δe_cmd and it is realized in the missile plant by an actuator, which is assumed to have a 1 st order dynamics bandwidth of about 15Hz. The gains used in the altitude hold autopilot system ( Fig. 3 ) are selected from the stability analysis of longitudinal dynamics, and the details are omitted in this paper.
4. Simulations
In this section, the controllers constructed in the preceding section are examined in several simulations. The simulations begin from an initial trim condition. For the trim condition, it is assumed that ramjet engine maintains normal operation during level flight at Mach 2.9 and an altitude of 14 km. In addition, it is noted that the initial guess and gains defined in Eqs. (18) and (19) are consistently available in the simulation.
To demonstrate the control algorithm, we consider two scenarios, in total, with different velocity commands. The first scenario is made for linear velocity and altitude commands. The results are shown in Fig. 4 . The velocity command is linearly increasing, equivalently, from Mach 2.9 to 3.0. After reaching Mach 3.0 the velocity command is maintained to that level. The altitude reference is also given linearly from 14 km with an increment rate of 5 m/sec. From the result, it is shown that the engine needs to consume more fuel supply with increasing velocity and altitude from trimmed states. After reaching reference commands, the maneuver of the vehicle (by equivalence ratio and missile fin deflection) is both stabilized with little variation.
In order to stabilize the error of the vehicle’s velocity for a reference command, the estimation of nonlinear controller results in fine tracking of net force (or net acceleration) along to the velocity vector. As described, the nonlinear velocity controller estimates the vehicle’s thrust and drag with the approximated functions in Eqs. (2) and (4). However, the nonlinear controller takes only one feedback, the velocity error, such that the number of estimation is not matched with that of the stabilization objectives. Therefore, it may not be guaranteed that the estimations both in thrust and drag are converged to their real values, but only the converged estimation of net force (or net acceleration) is secured (see Fig. 4 .(f )).
The second case is for the sinusoidal velocity reference wherein the other conditions are the same as in the first scenario. This simulation might help to verify the control performance with an arbitrarily varying control command. The reference command is assumed to be a sinusoid Mach function with an amplitude of 0.1, a frequency of 0.1 (rad/s), and a mean value of 2.9. The results are shown in Fig. 5 . From the simulation, it is shown that the velocity and altitude command are well tracked by both controllers for corresponding subsystems. The velocity change is delayed by roughly π/2(rad) to the fuel supply as expected in the motion
Lager Image
Simulation Results for Linear Velocity and Altitude Reference Command
dynamics, and the fin deflection actuation is consistent with the velocity, which needs to stabilize the altitude by the velocity change. Furthermore, Fig. 5 (f) demonstrates that the estimation of the net acceleration is continuously converged with time.
5. Conclusion
In this paper, an approach on the controller design of a ramjet-powered vehicle is introduced. The ramjet-powered vehicle has a coupling effect of the vehicle’s flight condition and propulsive performance such that the thrust is one of the major factors under varying instantaneous flight and engine conditions.
In the scheme for the controller design, the longitudinal motion of the vehicle is assumed and is divided into two subsystems (velocity and altitude) where they are separately stabilized by the controllers for the subsystems. This work concentrates on the nonlinear adaptive velocity controller design considering the ramjet engine propulsive characteristics as the vehicle’s altitude is stabilized with a classical altitude hold autopilot control structure. The control law for the velocity subsystem is derived by using approximation functions for the vehicle’s aerodynamic and propulsive properties, and they are inferred from the previous result of the study in Ref. [9] and the thrust tendency analysis from its database for the ramjet engine, respectively.
From the simulations, it is shown that the velocity tracking performance of the ramjet-powered vehicle is demonstrated even in the case regarding altitude change. Therefore, this work can be extended to control for reference trajectory, which is pre-designed from the results of a trajectory optimization problem.
By considering a nonlinear adaptive control structure, the complexity and difficulty in designing a gain-scheduled
Lager Image
Simulation Results for Sinosoidal Velocity and Linear Altitude Reference Commands
PID controller can be reduced as the nonlinear scheme does not require repetition in verifying the control performance according to certain design points for varying flight conditions and combustion characteristics. Therefore, as a preliminary research, it is expected that this work can be extended to applications for specific analyses on system performance under constraints such as the coupling effects between engine characteristics and flight dynamics.
Fry R. S. 2004 “A Century of Ramjet Propulsion Technology Evolution” Journal of Propulsion and Power 20 (1) 27 - 58
Boksenbom A. S. , Novik D. 1948 “Control Requirements and Control Parameters for A Ram Jet with Variable-Area Exhaust Nozzle” NACA RM-E8H24
Cotright E. M. , Connors J. F. 1952 "Survey of Some Preliminary Investigations of Supersonic Diffusers at High Mach Numbers" NACA RM-E52E20
Himmel S. C. 1952 “Some Control Considerations for Ram-jet Engines” NACA RM-E52F10
Tournes C. , Landrum D. B. , Shtessel Y. , Hawk C. W. 1998 “Ramjet-Powered Reusable Launch Vehicle Control by Sliding Modes” Journal of Guidance, Control, and Dynamics 21 (3) 409 - 415
Xu H. , Mirmirani M. D. , Ioannou P. A. 2004 “Adaptive Sliding Mode Control Design for a Hypersonic Flight Vehicle” Journal of Guidance, Control, and Dynamics 27 (5) 829 - 838
Gibson T. E. , Annaswamy A. M. 2008 “Adaptive Control of Hypersonic Vehicles in the Presence of Thrust and Actuator Uncertainties” AIAA Guidance, Navigation and control Conference and Exhibit St. Louis, HI
Fiorentini L. , Serrani A. , Bolender M. A. , Doman D. B. 2009 "Nonlinear Robust Adaptive Control of Flexible Air-Breathing Hypersonic Vehicles" Journal of Guidance, Control, and Dynamics 32 (2) 401 - 416    DOI : 10.2514/1.39210
Parker J. T. , Serrani A. , Yurkovich S. , Bolender M. A. , Doman D. B. 2007 "Control-Oriented Modeling of an Air- Breathing Hypersonic Vehicle" Journal of Guidance, Control, and Dynamics 30 (3) 856 - 869    DOI : 10.2514/1.27830
Bolender M. A. , Doman D. B. 2007 "Nonlinear Longitudinal Dynamical Model of an Air-Breathing Hypersonic Vehicle" Journal of Spacecraft and Rockets 44 (2) 374 - 387    DOI : 10.2514/1.23370
Park J. , Tahk M. , Sung H. 2009 “Trajectory Optimization for A Supersonic Air-Breathing Missile System Using Pseudo-Spectral Method” International Journal of Aeronautical and Space Sciences 10 (1) 112 - 121
Kim S. , Yeom H. , Jeon J. , Sung H. 2008 “An Integrated Dynamic Model to Control Ramjet Propulsion System” Proceedings of the KSAS-JSASS Joint International Symposium on Aerospace Engineering Jeju, Korea
Chandra K. P. B. , Gupta N. K. , Ananthkrishnan N. , Park I. , Yoon H. 2010 "Modeling, Simulation, and Controller Design for an Air-Breathing Combustion System" Journal of Propulsion and Power 26 (3) 562 - 574    DOI : 10.2514/1.42368