Ascent Phase Trajectory Optimization for a Hypersonic Vehicle Using Nonlinear Programming
|
|
- Berenice Harper
- 5 years ago
- Views:
Transcription
1 Ascent Phase Trajectory Optimization for a Hypersonic Vehicle Using Nonlinear Programming H.M. Prasanna,D.Ghose, M.S. Bhat, C. Bhattacharyya, and J. Umakant 3 Department of Aerospace Engineering Department of Computer Science and Automation 3 Aerodynamics Division, DRDL, Hyderabad, India, Indian Institute of Science, Bangalore 56, India prasannahm@rediffmail.com, {dghose, msbdcl}@aero.iisc.ernet.in, chiru@csa.iisc.ernet.in, umakantj@yahoo.com Abstract. In this paper we present a nonlinear programming solution to one of the most challenging problems in trajectory optimization. Unlike most aerospace trajectory optimization problems the ascent phase of a hypersonic vehicle has to undergo large changes in altitude and associated aerodynamic conditions. As a result, its aerodynamic characteristics, as well as its propulsion parameters, undergo drastic changes. Further, the data available through wind tunnel tests are not always smooth. The challenge in solving such problems lies both in the preprocessing of the data as well as in the judicious use of optimization techniques. In this paper we advocate approximation of the infinite dimensional optimal control problem, derived from practical considerations of a hypersonic vehicle ascending from an altitude of 6 kms to an altitude of 3 kms with specified mach numbers, into a set of finite dimensional nonlinear programming problems. This finite dimensional approximation is shown to produce acceptable optimized results in terms of angle-of-attack control histories and state behaviour. A modification, that exploits the ultimate scheme of linear interpolation of the optimal discrete history, is proposed and is shown to produce accurate results with smaller number of grid points. Introduction Ascent phase trajectory optimization of hypersonic vehicles is one of the most challenging real-life optimization problems []. Numerical algorithms for these problems are difficult to develop due to the inherent nature of the problem []. Our paper is an effort in this direction. Numerical solution of optimal control problems has two important components: optimization and numerical integration of differential equations [3, 4]. A naive approach to solve such problems is a cut-and-paste strategy that uses numerical integration and optimization packages in the form shown in Figure (a). A sophisticated O. Gervasi et al. (Eds.): ICCSA 5, LNCS 3483, pp , 5. c Springer-Verlag Berlin Heidelberg 5
2 Ascent Phase Trajectory Optimization 549 u u Integration package u J Optimization package J x~ i+ x i t t i i+ x~ i+ x i+ x i+ t i+ (a) (b) Fig.. (a) The naive cut-and-paste method (b) Direct transcription method approach exploits the intimate relation that exists between the methods used for solving differential equations and the methods used for optimization [3]. In this paper we will describe such an approach based on nonlinear programming (NLP) and illustrate the method with an example of optimization of the trajectory for an experimental hypersonic vehicle during its ascent phase. The basic approach in solving a trajectory optimization problem using NLP is to represent an infinite dimensional optimal control problem by one or more simpler finite dimensional subproblems. While using this technique, several practical issues arise that delineate these problems from textbook problems. In fact, taking care of these issues requires a sound knowledge of the physics behind the problem. Apart from the well understood issues of infeasible constraints, rank deficient constraints, constraint redundancy, and scaling, there are other issues of missing data in the aerodynamic and propulsion parameters as well as non-smoothness (which may be due to engineering errors or real) which may confuse standard spline interpolation methods. These problems require careful pre-processing of data before they can be used in the optimization process. Methods for Solving Optimal Control Problems. The Optimal Control Problem An optimal control problem may be formulated to find the control vector u(t) to minimize the performance index [3, 5] φ[y(t f ),t f ] () evaluated at the final time t f, which may be fixed or free. The dynamics of the system are given by the state equations and initial and final conditions, ẏ = f[y(t),u(t),t], χ[y(t ),u(t ),t ] χ =,ψ[y(t f ),u(t f ),t f ] ψ f = () where, y is the state vector. Also, we have state and control constraints, ξ L ξ[y(t),u(t),t] ξ U, y L y(t) y U, u L u(t) u U (3)
3 55 H.M. Prasanna et al.. Transcription Method The optimal control problem is an infinite dimensional extension of a finite dimensional NLP problem. This infinite dimensional problem can be transcribed to a finite dimensional approximation. The transcription approach treats the state and control variable at discrete times as optimization variables. The time points are defined as, t I = t <t <...<t M = t F. These points are referred to as node, mesh, or grid points. The values of the state and control variables at these intermediate steps are treated as a set of NLP variables. The differential equation is thus replaced by a finite set of constraints. The transcription method decomposes the dynamic trajectory into several phases or intervals. If there are n number of intervals, the number of grid points are n +. This method uses the value of both state and control at each grid point as NLP (nonlinear programming) variables. We may include initial or final time as one of the NLP variables. The continuity in the states across the phase boundary is achieved by forming an NLP constraint such that, the propagated value of the state from previous phase matches with the value of the state variable at current phase. In the direct transcription method the propagation is done by using algebraic formulas derived from numerical integration schemes. This makes the propagation function approximate (with increasing levels of accuracy depending on the order of the integration scheme). However, as the interval size decreases the accuracy will increase. Let the propagation scheme be given by a function P (.) so that the propagated state at t i+ be defined as x i+. Then P (x i )= x i+ (4) where, x i is the assumed state at time t i. The constraint at t i+ is, x i+ x i+ = x i+ P (x i ) = (5) Figure (b) explains the concept of NLP formulation. The choice of P is important. The control which is obtained from solving NLP problem is interpolated to make it continuous. The accuracy of discrete trajectory, obtained from solving the NLP problem, is assessed by using the interpolated control to integrate the state equations (ODEs). It is expected that as the number of intervals increases, the error between the simulation result and optimization will reduce..3 Discretization Methods Euler Method: This is a first order method and the constraints are given by, =y k+ y k h k f k c(x) (6) Trapezoidal Method: This is a second order method and the constraints are, =y k+ y k (h k /)[f k + f k+ ] c(x) (7) Hermite-Simpson (Compressed) (HSC): This is of order four and the constraints are given by, =y k+ y k (h k /6)[f k+ +4f((y k+ + y k )/+h k (f k f k+ )/8, ū k+,t k + h k /) + f k ] c(x) (8)
4 Ascent Phase Trajectory Optimization 55 The NLP variables are given as x T =(t k,y,u, ū,y,u, ū 3,...,ū M,y M,u M ). Hermite-Simpson (Separated) (HSS): The constraints are given by, =ȳ k+ (/)(y k+ + y k ) (/8)h k (f k f k+ ) (9) =y k+ y k (/6)h k (f k+ +4 f k+ + f k ) () Eqn. (9) defines the hermite interpolation for the state at the interval midpoint, while () enforces the Simpson quadrature over the interval. The NLP variables are, x T =(t k,y,u, ū,y,u, ȳ 3, ū 3,...,ȳ M, ū M,y M,u M ). In general, the duration of a phase may be variable, in which case, the set of NLP variables must be augmented to include the variable time(s) t I and/or t F. Also, we must alter the discretization such that the step size is h k = τ k (t F t I )= τ k t, where, t =(t F t I ) with constants <τ k < chosen so that the grid points are located at fixed fractions of the total phase duration. 3 Basic Models for Ascent Phase Trajectory Optimization of a Hypersonic Vehicle The experimental hypersonic vehicle is air dropped from a carrier aircraft and boosted to Mach number 3.5 at an altitude of 6 km using solid rocket motors. The boosters are separated upon achieving this condition. Then the ascent phase commences and the final condition of Mach number 6. and altitude 3 km is achieved. The problem is to optimize the ascent phase trajectory so that the least fuel is expended (so that maximum remaining fuel can be utilized during the subsequent cruise phase). The constraints and other parameters of the problem are angle of attack α 8, maximum range during ascent phase 5 km, maximum time to reach the cruise conditions 5 s, maximum longitudinal acceleration during ascent phase < 4g, maximum lateral acceleration < 5g, weight of empty vehicle= kg (excluding fuel), fuel weight = 6 kg. The aerodynamic characteristics [6] of the vehicle are usually available in tabular form as values of the normal force coefficient (c N ) and axial force coefficient (c A ) which are functions of Mach number (M) and angle of attack (α). The lift and drag coefficients are computed using the relations c L = c N cos α c A sin α and c D = c N sin α + c A cos α. The resulting variation of lift and drag coefficients with mach number and alpha are shown in Figure. Mass Flow Rate: The mass flow rate of air (ṁ air ) as function of Mach number and α is given in Table (a), for density corresponding to altitude 3.5 km. The mass flow rate of air at any altitude is, ṁ air = ρva, where, A is nozzle area, v is free stream velocity at the nozzle, and ρ is density at given altitude. At 3.5 km, let the density be ρ 3.5, speed of sound be S 3.5, and the speed corresponding to a given mach numbers M is V 3.5 = MS 3.5. We non-dimensionalize Table (a) and determine ṁ air as,
5 55 H.M. Prasanna et al..5 c L.5 α α = = 8 8 α = 6 c D α = 8 α = 6 α = 4 α = 4.5 α =.5 α = α =.4.3 α = Mach number Mach number.5 M = 3.5 M= 4. M = 5. M=6. M = M = 4. M= c L c D M = 6.5 M = 6. M = Fig.. Variation of lift and drag coefficients with α and Mach number Table. (a) Mass flow rate of air (-kg/s) corresponding to altitude 3.5 km (b) Specific impulse (I sp in seconds) M/α (a) Altitude (km) M (b) ṁ nond = ṁ air /(ρ 3.5 V 3.5 A) () ṁ }{{} air = ṁ nond ρva = ṁ }{{} air /(ρ 3.5 V 3.5 ) ρv () }{{} at given altitude at given altitude at given altitude Specific Impulse: Specific impulse is a function of free stream Mach number and altitude (Table (b)). For Mach number between 3. to 4., specific impulse
6 Ascent Phase Trajectory Optimization x Specific Impulse (s) Altitude (m) x 4 Altitude (m) Mach number (a) Temperature (k) (b) Fig. 3. (a) Specific impulse using spline interpolation (b) Temperature variation in Indian standard atmosphere is constant. So linear interpolation is used for this region and spline interpolation for the rest of the data (Figure 3(a)). Thrust: Thrust is calculated as T = I sp ṁ f g. Assuming equivalence ratio φ= and stoichiometric ratio for fuel as 5, we have ṁ f = ṁ air /5, where, ṁ air at given altitude is obtained from the procedure described above. Standard Atmosphere: The Indian standard atmosphere is assumed. The variation of temperature with altitude is shown in Figure 3(b). The atmospheric pressure at sea level is P = 5 (N/m ) and specific gas constant is R= (J/kgK). 4 Problem Formulation and Computational Results During ascent phase the objective is to minimize fuel consumption (or maximize mass at final time). Thus, the optimization problem is formulated as, max m f (3) s.t., ḣ = v sin γ, ẋ = v cos γ, ṁ = T/(I sp g) (4) v =(/m)[t cos α D] µ sin γ/(r e + h) (5) γ =(/(mv))[t sin α + L] + cos γ[v/(r e + h) µ/v(r e + h) ] (6) where, h is the altitude (m), v the velocity (m/sec), α the angle of attack (rad), γ the flight path angle (rad), m the mass (kg), µ the gravitational constant, g the gravitational acceleration, and R e the radius of the earth. The trajectory related constraints t f 5 s, x f 5 m, α 8, a long 4g, a lat 5g
7 554 H.M. Prasanna et al. are also imposed. where, a long and a lat are the longitudinal and lateral accelerations, respectively. defined as a lat =(/m)[t sin α + L] + cos γ[v /(R e + h) µ/(r e + h) ] (7) a long =(/m)[t cos α D] µ sin γ/(r e + h) (8) The lift (L) and drag (D) forces are defined by, D =(/)ρsv c D, L =(/)ρsv c L (9) where, S is the aerodynamic reference area of the vehicle, and ρ is the atmospheric density. The following constants complete the formulation of the problem: h() = 6. m, h(t f ) = 3. m, v() = m/s v(t f ) = 847. m/s, γ() =.698 rad, γ(t f )=rad m() = 36. kg, x() = m, µ = m 3 /s g =9.8 m/s, R e = m, S =m The optimization was started with 5 grid points initially and subsequently higher number of grid points (i.e., 9, 7, 33, 65) were considered. For 5 and 9 grid points, Table. Results of ascent phase trajectory optimization Grid method h f (m) v f (m/s) γ f ( o ) m f (kg) m f(fuel) x f ( 5 m) time (s) points (kg) 5 NLP Step Linear Hermite NLP Step Linear Hermite NLP Step Linear Hermite NLP Step Linear Hermite NLP Step Linear Hermite
8 Ascent Phase Trajectory Optimization 555 the NLP constraints are formed using trapezoidal discretization and for higher number of grid points Hermite-Simpson discretization is used. Each variable is scaled by its upper bound to make the ranges, over which they vary, similar in their order of magnitude. This helps better numerical conditioning of the NLP. The optimization is carried out only at the grid points and the shape of the trajectory between the grid points is not known. To assess the accuracy of the result, we integrated (4) - (6) using control histories obtained from interpolation of the NLP generated discrete control. We used three interpolation schemes (step, linear, and hermite). The initial conditions for n grid points were obtained by linearly interpolating the optimized values from n grid points. The results of optimization and integration (using different interpolation schemes for control) are given in Table. From Table, we observe that the error between the optimization and integration results are significantly large when number of grid points are few and the error reduces as the number of grid points increases. The integration using control generated by step interpolation has large error as compared to linear and hermite interpolations. The results for 33, and 65 grid points show very little improvement in the optimization results (see Figure 4). Figure 5 shows variation of altitude, velocity, flight path angle, mass, and range with time using linearly interpolated control for 65 grid points (a) (b) (c) (d) (e) (f) Fig. 4. (a) (e) Angle of attack vs. time for 5, 9, 7, 33, and 65 grid points in the original NLP (f) Angle of attack vs. time for 65 grid points with modified method
9 556 H.M. Prasanna et al x Altitude (m) Velocity (m/s) x 5 Flight path angle ( γ in deg) 6 4 Overall mass (kg) Horizontal range (m) Fig. 5. State histories (for linearly interpolated control from 65 grid points) 5 Improving Computational Accuracy Since the final step in the transcription method is implemented by linearly interpolating the control history between grid points, we use this knowledge by assuming a linear dynamics for control variable. So, we assume, α = β () Table 3. Results of ascent phase trajectory optimization (with modified method) Grid method h f (m) v f (m/s) γ f ( o ) m f (kg) m f(fuel) x f ( 5 m) time (s) points (kg) 5 NLP Linear NLP Linear NLP Linear NLP Linear NLP Linear e
10 Ascent Phase Trajectory Optimization 557 where, β is now used as the control variable. Then, the equations of motion will include () in addition to (4) - (6). The set of NLP variables will consequently include an additional variable β at each grid point. With this modification, the results of optimization and integration are given in Table 3. Figure 4(f) shows variation of angle of attack against time for 65 grid points. Comparing with the α history in the earlier scheme we see a flat region between 7 and 9 seconds in the new scheme. The modified method also gives better accuracy. In fact it does so for just 33 grid points, whereas the original scheme produces similar accurate results at 65 grid points. However, the additional NLP variable and constraint increases the computational time slightly. 6 Conclusions In this paper, we addressed the challenging application problem of ascent phase trajectory optimization of a hypersonic vehicle. It was shown that the linear interpolation scheme used to create the continuous control history for verification can itself be exploited while forming the dynamical equations to get faster convergence and higher accuracy. Acknowledgements. This work was partially supported by a grant from DRDL. References. P.F. Gath and A.J. Calise, Optimization of Launch Vehicle Ascent Trajectories with Path Constraints and Coast Arcs, Journal of Guidance, Control, and Dynamics, Vol. 4,, pp J.T. Betts, Survey of Numerical Methods for Trajectory Optimization, Journal of Guidance, Control, and Dynamics, Vol., 998, pp J.T. Betts, Practical Methods for Optimal Control Using Nonlinear Programming, The Society for Industrial and Applied Mathematics, Philadelphia,. 4. E. Gill, Murray, and H. Wright, Practical Optimization, Academic Press, A.E. Bryson, and Y.-C. Ho, Applied Optimal Control, John Wiley, New York, J.D. Anderson, Introduction to Flight, McGraw-Hill, Singapore,.
The Direct Transcription Method For Optimal Control. Part 2: Optimal Control
The Direct Transcription Method For Optimal Control Part 2: Optimal Control John T Betts Partner, Applied Mathematical Analysis, LLC 1 Fundamental Principle of Transcription Methods Transcription Method
More informationAEROSPACE ENGINEERING
AEROSPACE ENGINEERING Subject Code: AE Course Structure Sections/Units Topics Section A Engineering Mathematics Topics (Core) 1 Linear Algebra 2 Calculus 3 Differential Equations 1 Fourier Series Topics
More informationMinimum Time Ascent Phase Trajectory Optimization using Steepest Descent Method
IJCTA, 9(39), 2016, pp. 71-76 International Science Press Closed Loop Control of Soft Switched Forward Converter Using Intelligent Controller 71 Minimum Time Ascent Phase Trajectory Optimization using
More informationEXTERNAL-JET (FLUID) PROPULSION ANALOGY FOR PHOTONIC (LASER) PROPULSION By John R. Cipolla, Copyright February 21, 2017
EXTERNAL-JET (FLUID) PROPULSION ANALOGY FOR PHOTONIC (LASER) PROPULSION By John R. Cipolla, Copyright February 21, 2017 ABSTRACT External-jet propulsion uses a narrow jet of high velocity water or conceptually
More informationIntroduction to Flight Dynamics
Chapter 1 Introduction to Flight Dynamics Flight dynamics deals principally with the response of aerospace vehicles to perturbations in their flight environments and to control inputs. In order to understand
More informationCalifornia State Science Fair
California State Science Fair Working Model for Model Rocket Altitude Prediction Edward Ruth drruth@ix.netcom.com This is a complete model of all the forces acting on a model rocket in flight. It calculates
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationStudy of Required Thrust Profile Determination of a Three Stages Small Launch Vehicle
Journal of Physics: Conference Series PAPER OPEN ACCESS Study of Required Thrust Profile Determination of a Three Stages Small Launch Vehicle To cite this article: A Fariz et al 218 J. Phys.: Conf. Ser.
More informationRocket Science 102 : Energy Analysis, Available vs Required
Rocket Science 102 : Energy Analysis, Available vs Required ΔV Not in Taylor 1 Available Ignoring Aerodynamic Drag. The available Delta V for a Given rocket burn/propellant load is ( ) V = g I ln 1+ P
More informationTrajectory Optimization for Ascent and Glide Phases Using Gauss Pseudospectral Method
Trajectory Optimization for Ascent and Glide Phases Using Gauss Pseudospectral Method Abdel Mageed Mahmoud, Chen Wanchun, Zhou Hao, and Liang Yang Abstract The trajectory optimization method for ascent
More informationThe Role of Zero Dynamics in Aerospace Systems
The Role of Zero Dynamics in Aerospace Systems A Case Study in Control of Hypersonic Vehicles Andrea Serrani Department of Electrical and Computer Engineering The Ohio State University Outline q Issues
More informationOptimal Trajectories of Anti-satellite Missiles
Licensed under the Creative Commons Attribution Noncommercial No Derivatives Communications to SIMAI Congress, ISSN 1827-9015, Vol. 3 (2009) 265 (12pp) DOI: 10.1685/CSC09265 Optimal Trajectories of Anti-satellite
More informationInitial Trajectory and Atmospheric Effects
Initial Trajectory and Atmospheric Effects G. Flanagan Alna Space Program July 13, 2011 Introduction A major consideration for an earth-based accelerator is atmospheric drag. Drag loses mean that the gun
More informationReal-time trajectory generation technique for dynamic soaring UAVs
Real-time trajectory generation technique for dynamic soaring UAVs Naseem Akhtar James F Whidborne Alastair K Cooke Department of Aerospace Sciences, Cranfield University, Bedfordshire MK45 AL, UK. email:n.akhtar@cranfield.ac.uk
More informationHow Small Can a Launch Vehicle Be?
UCRL-CONF-213232 LAWRENCE LIVERMORE NATIONAL LABORATORY How Small Can a Launch Vehicle Be? John C. Whitehead July 10, 2005 41 st AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit Tucson, AZ Paper
More informationClassical Numerical Methods to Solve Optimal Control Problems
Lecture 26 Classical Numerical Methods to Solve Optimal Control Problems Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Necessary Conditions
More informationFundamentals of Airplane Flight Mechanics
David G. Hull Fundamentals of Airplane Flight Mechanics With 125 Figures and 25 Tables y Springer Introduction to Airplane Flight Mechanics 1 1.1 Airframe Anatomy 2 1.2 Engine Anatomy 5 1.3 Equations of
More informationChapter 5 Performance analysis I Steady level flight (Lectures 17 to 20) Keywords: Steady level flight equations of motion, minimum power required,
Chapter 5 Performance analysis I Steady level flight (Lectures 17 to 20) Keywords: Steady level flight equations of motion, minimum power required, minimum thrust required, minimum speed, maximum speed;
More informationEXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING. Dilmurat Azimov
EXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING Dilmurat Azimov University of Hawaii at Manoa 254 Dole Street, Holmes 22A Phone: (88)-956-2863, E-mail: azimov@hawaii.edu
More informationBasic Ascent Performance Analyses
Basic Ascent Performance Analyses Ascent Mission Requirements Ideal Burnout Solution Constant & Average Gravity Models Gravity Loss Concept Effect of Drag on Ascent Performance Drag Profile Approximation
More informationTechnology of Rocket
Technology of Rocket Parts of Rocket There are four major parts of rocket Structural system Propulsion system Guidance system Payload system Structural system The structural system of a rocket includes
More informationTurn Performance of an Air-Breathing Hypersonic Vehicle
Turn Performance of an Air-Breathing Hypersonic Vehicle AIAA Aircraft Flight Mechanics Conference Derek J. Dalle, Sean M. Torrez, James F. Driscoll University of Michigan, Ann Arbor, MI 4809 August 8,
More informationGoddard s Problem Steven Finch
Goddard s Problem Steven Finch June 12, 2015 Arocketliftsoff vertically at time =0. Let ( ) be the mass of the rocket (payload and fuel) and ( ) be the altitude. We wish to choose the thrust ( ) and a
More informationModelling Nozzle throat as Rocket exhaust
Vol. 3, Issue. 4, Jul - Aug. 2013 pp-2502-2506 ISSN: 2249-6645 Modelling Nozzle throat as Rocket exhaust Keshava Rao P. 1, Komma Rahul 2, Souda Dinesh 3 1 (Mechanical Engineering, CBIT College, India)
More informationFinite Set Control Transcription for Optimal Control Applications
Finite Set Control Transcription for Optimal Control Stuart A. Stanton 1 Belinda G. Marchand 2 Department of Aerospace Engineering and Engineering Mechanics The University of Texas at Austin 19th AAS/AIAA
More informationMultistage Rocket Performance Project Two
41 Multistage Rocket Performance Project Two Charles R. O Neill School of Mechanical and Aerospace Engineering Oklahoma State University Stillwater, OK 74078 Project Two in MAE 3293 Compressible Flow December
More informationFlight Trajectory Optimization Using Genetic Algorithm Combined with Gradient Method
ITEM, Volume 1, No 1, 2001, Paper 10 e-journal: Information Technology for Economics & Management ISSN 1643-8949 Flight Trajectory Optimization Using Genetic Algorithm Combined with Gradient Method Nobuhiro
More informationCompressible Flow. Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist
Compressible Flow Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist What is Compressible Flow? Compressible Flow is a type of flow in which the density can not be treated as constant.
More informationDynamics and Control Preliminary Examination Topics
Dynamics and Control Preliminary Examination Topics 1. Particle and Rigid Body Dynamics Meirovitch, Leonard; Methods of Analytical Dynamics, McGraw-Hill, Inc New York, NY, 1970 Chapters 1-5 2. Atmospheric
More informationPerformance. 5. More Aerodynamic Considerations
Performance 5. More Aerodynamic Considerations There is an alternative way of looking at aerodynamic flow problems that is useful for understanding certain phenomena. Rather than tracking a particle of
More informationNUMERICAL INVESTIGATIONS ON THE SLENDER AXISYMMETRIC BODIES AERODYNAMICS IN WIDE RANGE OF MACH NUMBERS AND ANGLES OF ATTACK FROM 0 TO 180
NUMERICAL INVESTIGATIONS ON THE SLENDER AXISYMMETRIC BODIES AERODYNAMICS IN WIDE RANGE OF MACH NUMBERS AND ANGLES OF ATTACK FROM 0 TO 180 N.V. Voevodenko*, L.G. Ivanteeva*, V.Ju. Lunin* * TsAGI Central
More informationME 6139: High Speed Aerodynamics
Dr. A.B.M. Toufique Hasan Professor Department of Mechanical Engineering, BUET Lecture-01 04 November 2017 teacher.buet.ac.bd/toufiquehasan/ toufiquehasan@me.buet.ac.bd 1 Aerodynamics is the study of dynamics
More informationLaunch abort trajectory optimisation for reusable launch vehicles
Launch abort trajectory optimisation for reusable launch vehicles Federico Toso and Christie Alisa Maddock Centre for Future Air Space Transportation Technologies University of Strathclyde, Glasgow, G1
More informationStability and Control Some Characteristics of Lifting Surfaces, and Pitch-Moments
Stability and Control Some Characteristics of Lifting Surfaces, and Pitch-Moments The lifting surfaces of a vehicle generally include the wings, the horizontal and vertical tail, and other surfaces such
More information. Frobenius-Perron Operator ACC Workshop on Uncertainty Analysis & Estimation. Raktim Bhattacharya
.. Frobenius-Perron Operator 2014 ACC Workshop on Uncertainty Analysis & Estimation Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. uq.tamu.edu
More informationRobotic Mobility Atmospheric Flight
Gaseous planetary environments (Mars, Venus, Titan) Lighter-than- air (balloons, dirigibles) Heavier-than- air (aircraft, rotorcraft) 1 2018 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
More informationAAE 251 Formulas. Standard Atmosphere. Compiled Fall 2016 by Nicholas D. Turo-Shields, student at Purdue University. Gradient Layer.
AAE 51 Formulas Compiled Fall 016 by Nicholas D. Turo-Shields, student at Purdue University Standard Atmosphere p 0 = 1.0135 10 5 Pascals ρ 0 = 1.5 kg m 3 R = 87 J kg K γ = 1.4 for air p = ρrt ; Equation
More informationI. Introduction. external compression. supersonic flow. II. Design Criteria
Design Optimization of High Speed Inlets Doyle D Knight Dept of Mechanical and Aerospace Engineering Rutgers - The State University of New Jersey New Brunswick, NJ 08903 knight@soemailrutgersedu I Introduction
More informationLecture1: Characteristics of Hypersonic Atmosphere
Module 1: Hypersonic Atmosphere Lecture1: Characteristics of Hypersonic Atmosphere 1.1 Introduction Hypersonic flight has special traits, some of which are seen in every hypersonic flight. Presence of
More informationControl Theory & Applications. Re-entry trajectory optimization using Radau pseudospectral method. HAN Peng, SHAN Jia-yuan
30 8 013 8 DOI: 10.7641/CTA.013.1041 Control Theory & Applications Vol. 30 No. 8 Aug. 013 Radau ( 100081) : Radau. Legendre-Gauss-Radau Lagrange Legendre-Gauss-Radau. (NLP) NLP SNOPT... : ; ; ; Radau :
More informationChapter 1 Lecture 2. Introduction 2. Topics. Chapter-1
Chapter 1 Lecture 2 Introduction 2 Topics 1.4 Equilibrium of airplane 1.5 Number of equations of motion for airplane in flight 1.5.1 Degrees of freedom 1.5.2 Degrees of freedom for a rigid airplane 1.6
More informationTrajectory Optimization for a Mars Ascent Vehicle
AIAA SPACE Forum 13-16 September 216, Long Beach, California AIAA/AAS Astrodynamics Specialist Conference 216. All rights reserved 1.2514/6.216-5441 Trajectory Optimization for a Mars Ascent Vehicle Joel
More informationPreface. 2 Cable space accelerator 39
Contents Abstract Preface xiii xv 1 Space elevator, transport system for space elevator, 1 and tether system 1.1 Brief history 1 1.2 Short description 2 1.3 Transport system for the space elevator 5 1.4
More informationCalculus of Variations Summer Term 2015
Calculus of Variations Summer Term 2015 Lecture 14 Universität des Saarlandes 24. Juni 2015 c Daria Apushkinskaya (UdS) Calculus of variations lecture 14 24. Juni 2015 1 / 20 Purpose of Lesson Purpose
More informationIntroduction to Aerospace Engineering
Introduction to Aerospace Engineering 5. Aircraft Performance 5.1 Equilibrium Flight In order to discuss performance, stability, and control, we must first establish the concept of equilibrium flight.
More informationFlight Test Results for Circular Path Following by Model Predictive Control
Preprints of the 19th World Congress The International Federation of Automatic Control Flight Test Results for Circular Path Following by Model Predictive Control Yoshiro Hamada Taro Tsukamoto Shinji Ishimoto
More informationMissile Interceptor EXTROVERT ADVANCED CONCEPT EXPLORATION ADL P Ryan Donnan, Herman Ryals
EXTROVERT ADVANCED CONCEPT EXPLORATION ADL P- 2011121203 Ryan Donnan, Herman Ryals Georgia Institute of Technology School of Aerospace Engineering Missile Interceptor December 12, 2011 EXTROVERT ADVANCED
More informationA Balance for Measurement of Yaw, Lift and Drag on a Model in a Hypersonic Shock Tunnel
, July 6-8, 2011, London, U.K. A Balance for Measurement of Yaw, Lift and Drag on a Model in a Hypersonic Shock Tunnel S. Trivedi, and V. Menezes Abstract This paper describes the design of an accelerometer
More informationChapter 4 The Equations of Motion
Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations
More informationIntroduction to Flight
l_ Introduction to Flight Fifth Edition John D. Anderson, Jr. Curator for Aerodynamics, National Air and Space Museum Smithsonian Institution Professor Emeritus University of Maryland Me Graw Higher Education
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and orbital mechanics Flight Mechanics practice questions Dr. ir. Mark Voskuijl 20-11-2013 Delft University of Technology Challenge
More informationContinuous Differentiation of Complex Systems Applied to a Hypersonic Vehicle
Continuous of Complex Systems Applied to a Vehicle AIAA Aircraft Flight Mechanics Conference Derek J. Dalle, Sean M. Torrez, James F. Driscoll University of Michigan, Ann Arbor, MI 4819 August 15, 212,
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationRockets, Missiles, and Spacecrafts
36 1 Rockets, Missiles, and Spacecrafts 2 Chinese used rockets in the 12 th century AD against the Mongol attacks. In India Tipu Sultan used rockets against the British army in the 18 th century. The modern
More informationRobotic Mobility Atmospheric Flight
Robotic Mobility Atmospheric Flight Gaseous planetary environments (Mars, Venus, Titan)! Lighter-than- air (balloons, dirigibles)! Heavier-than- air (aircraft, rotorcraft) 1 2014 David L. Akin - All rights
More informationAC : A DESIGN-BY-ANALYSIS PROJECT FOR INTRODUC- TORY STUDENTS IN AEROSPACE ENGINEERING
AC 2012-4116: A DESIGN-BY-ANALYSIS PROJECT FOR INTRODUC- TORY STUDENTS IN AEROSPACE ENGINEERING Dr. Mark Anderson, University of California, San Diego c American Society for Engineering Education, 2012
More informationA trendsetting Micro-Launcher for Europe
A trendsetting Micro-Launcher for Europe Farid Gamgami German Aerospace Center (DLR), Space Launcher Systems Analysis (SART) Bremen, Germany Farid.Gamgami@dlr.de ABSTRACT In this paper we analyse a potential,
More informationUniversity of California at Berkeley Department of Mechanical Engineering ME 163 ENGINEERING AERODYNAMICS FINAL EXAM, 13TH DECEMBER 2005
University of California at Berkeley Department of Mechanical Engineering ME 163 ENGINEERING AERODYNAMICS FINAL EXAM, 13TH DECEMBER 2005 Answer both questions. Question 1 is worth 30 marks and question
More informationPerformance Characterization of Supersonic Retropropulsion for Application to High-Mass Mars Entry, Descent, and Landing
Performance Characterization of Supersonic Retropropulsion for Application to High-Mass Mars Entry, Descent, and Landing Ashley M. Korzun 1 and Robert D. Braun 2 Georgia Institute of Technology, Atlanta,
More informationVerified High-Order Optimal Control in Space Flight Dynamics
Verified High-Order Optimal Control in Space Flight Dynamics R. Armellin, P. Di Lizia, F. Bernelli-Zazzera K. Makino and M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16
More informationFlight Dynamics and Control. Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege
Flight Dynamics and Control Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege Previously on AERO0003-1 We developed linearized equations of motion Longitudinal direction
More informationDesign And Analysis Of Thrust Chamber Of A Cryogenic Rocket Engine S. Senthilkumar 1, Dr. P. Maniiarasan 2,Christy Oomman Jacob 2, T.
Design And Analysis Of Thrust Chamber Of A Cryogenic Rocket Engine S. Senthilkumar 1, Dr. P. Maniiarasan 2,Christy Oomman Jacob 2, T. Vinitha 2 1 Research Scholar, Department of Mechanical Engineering,
More informationTHE METEOROLOGICAL ROCKET SYSTEM FOR ATMOSPHERIC RESEARCH
THE METEOROLOGICAL ROCKET SYSTEM FOR ATMOSPHERIC RESEARCH Komissarenko Alexander I. (1) Kuznetsov Vladimir M. (1) Filippov Valerii V. (1) Ryndina Elena C. (1) (1) State Unitary Enterprise KBP Instrument
More informationFlight Mechanics of Ram-Scram Transition
Flight Mechanics of Ram-Scram Transition Derek J. Dalle and James F. Driscoll University of Michigan, Ann Arbor, MI 48109 Among known air-breathing propulsion systems, the system that can be operated over
More informationOptimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming
Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology
More informationApplied Thermodynamics - II
Gas Turbines Sudheer Siddapureddy sudheer@iitp.ac.in Department of Mechanical Engineering Jet Propulsion - Classification 1. A heated and compressed atmospheric air, mixed with products of combustion,
More informationIAC-13,C4,P,44.p1,x17254 THE DYNAMIC OPERATON OF A HIGH Q EMDRIVE MICROWAVE THRUSTER. Roger Shawyer C.Eng. MIET. FRAeS. SPR Ltd UK
IAC-13,C4,P,44.p1,x1754 THE DYNAMIC OPERATON OF A HIGH Q EMDRIVE MICROWAVE THRUSTER Roger Shawyer C.Eng. MIET. FRAeS SPR Ltd UK sprltd@emdrive.com ABSTRACT The static operation of an EmDrive microwave
More informationResearch on Propeller Characteristics of Tip Induced Loss
4th International Conference on Machinery, Materials and Information Technology Applications (ICMMITA 2016) Research on Propeller Characteristics of Tip Induced Loss Yang Song1, a, Peng Shan2, b 1 School
More information2 Navier-Stokes Equations
1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1
More informationComputational Fluid Dynamics Analysis of Advanced Rocket Nozzle
Research Article International Journal of Current Engineering and Technology E-ISSN 2277 4106, P-ISSN 2347-5161 2014 INPRESSCO, All Rights Reserved Available at http://inpressco.com/category/ijcet Computational
More informationDynamic Inversion Design II
Lecture 32 Dynamic Inversion Design II Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics Summary of Dynamic Inversion Design Advantages Issues
More informationRobotic Mobility Atmospheric Flight
Robotic Mobility Atmospheric Flight Gaseous planetary environments (Mars, Venus, Titan) Lighter-than- air (balloons, dirigibles) Heavier-than- air (aircraft, rotorcraft) 1 2014 David L. Akin - All rights
More informationMinimum-Fuel Powered Descent for Mars Pinpoint Landing
JOURNAL OF SPACECRAFT AND ROCKETS Vol. 44, No. 2, March April 27 Minimum-Fuel Powered Descent for Mars Pinpoint Landing Ufuk Topcu, Jordi Casoliva, and Kenneth D. Mease University of California, Irvine,
More informationAerodynamics Simulation of Hypersonic Waverider Vehicle
Modern Applied Science February, 9 Aerodynamics Simulation of Hypersonic Waverider Vehicle Dingyi Wu & Hong Xiao School of Power and Energy Northwestern Polytechnical University PO box 85, Xian, Shaanxi,
More informationA THREE-DIMENSIONAL PREDICTIVE ENTRY GUIDANCE APPROACH
AIAA 2-3959 A THREE-DIMENSIONAL PREDICTIVE ENTRY GUIDANCE APPROACH K. D. Mease andd.t.chen Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697 S. Tandon, D.H.Young,
More informationEML 6934 Optimal Control Course Project
EML 6934 Optimal Control Course Project Spring 2014 1 Overview of Course Project The project for the course is structured as follows. You are required to choose two problems from the two lists shown below.
More informationAEROSPACE ENGINEERING DEPARTMENT. Second Year - Second Term ( ) Fluid Mechanics & Gas Dynamics
AEROSPACE ENGINEERING DEPARTMENT Second Year - Second Term (2008-2009) Fluid Mechanics & Gas Dynamics Similitude,Dimensional Analysis &Modeling (1) [7.2R*] Some common variables in fluid mechanics include:
More informationMission Scenarios for a Controlled Lunar Impact of a Small Satellite
IAC-4-IAA.4.11.P.5 Mission Scenarios for a Controlled Lunar Impact of a Small Satellite Nikolas Trawny, Michael Graesslin, Rene Laufer and Hans-Peter Roeser Email: n.trawny@gmx.de, {graesslin,laufer,roeser}@irs.uni-stuttgart.de
More informationEnabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty
Engineering, Test & Technology Boeing Research & Technology Enabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty ART 12 - Automation Enrique Casado (BR&T-E) enrique.casado@boeing.com
More informationInvestigation of Combined Airbreathing/Rocket. Air Launch of Micro-Satellites from a Combat Aircraft
6th Responsive Space Conference AIAA-RS6-008-5003 Investigation of Combined Airbreathing/Rocket Propulsion for Air Launch of Micro-Satellites from a Combat Aircraft Avichai Socher and Alon Gany Faculty
More informationEntry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design
Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
More informationNonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu
Nonlinear Filtering With Polynomial Chaos Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu Nonlinear Filtering with PC Problem Setup. Dynamics: ẋ = f(x, ) Sensor Model: ỹ = h(x)
More informationContinuation from a flat to a round Earth model in the coplanar orbit transfer problem
Continuation from a flat to a round Earth model in the coplanar orbit transfer problem M. Cerf 1, T. Haberkorn, Emmanuel Trélat 1 1 EADS Astrium, les Mureaux 2 MAPMO, Université d Orléans First Industrial
More informationSOLUTIONS MANUAL. ROCKET PROPULSION ELEMENTS, 8 th EDITION. By George P. Sutton and Oscar Biblarz. Published by John Wiley & Sons, Inc.
SOLUTIONS MANUAL to ROCKET PROPULSION ELEMENTS, 8 th EDITION By George P. Sutton and Oscar Biblarz Published by John Wiley & Sons, Inc. in 1 This manual is in part an outgrowth of courses taught by Prof.
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationFeedback Control of Aerodynamic Flows
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-843 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Feedback Control of Aerodynamic
More informationContents. Preface... xvii
Contents Preface... xvii CHAPTER 1 Idealized Flow Machines...1 1.1 Conservation Equations... 1 1.1.1 Conservation of mass... 2 1.1.2 Conservation of momentum... 3 1.1.3 Conservation of energy... 3 1.2
More informationEFFECT OF ATMOSPHERIC ALTITUDE ON THE DRAG OF WING AT SUBSONIC AND SUPERSONIC SPEEDS
Journal of Engineering Science and Technology 6 th EURECA 2016 Special Issue May (2017) 71-83 School of Engineering, Taylor s University EFFECT OF ATMOSPHERIC ALTITUDE ON THE DRAG OF WING AT SUBSONIC AND
More informationDynamic trajectory control of gliders
Dynamic trajectory control of gliders Rui DILãO and João FONSECA Institut des Hautes Études Scientifiques 35, route de Chartres 91440 Bures-sur-Yvette (France) Avril 2013 IHES/M/13/09 Dynamic trajectory
More informationL eaving Earth and arriving at another planet or asteroid requires
Designing Interplanetary Transfers L eaving Earth and arriving at another planet or asteroid requires a spacecraft to implement a sequence of manoeuvres. These include changes of velocity needed to escape
More information03. Simple Dynamical Systems
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 03. Simple Dynamical Systems Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationReview of Fundamentals - Fluid Mechanics
Review of Fundamentals - Fluid Mechanics Introduction Properties of Compressible Fluid Flow Basics of One-Dimensional Gas Dynamics Nozzle Operating Characteristics Characteristics of Shock Wave A gas turbine
More informationDynamic trajectory control of gliders
Proceedings of the EuroGNC 2013, 2nd CEAS Specialist Conference on Guidance, Navigation & Control, Delft University of Technology, Delft, The Netherlands, April 10-12, 2013 ThCT1.3 Dynamic trajectory control
More informationAn hp-adaptive pseudospectral method for solving optimal control problems
OPTIMAL CONTROL APPLICATIONS AND METHODS Optim. Control Appl. Meth. 011; 3:476 50 Published online 6 August 010 in Wiley Online Library (wileyonlinelibrary.com)..957 An hp-adaptive pseudospectral method
More informationfor what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory?
1. 5% short answers for what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory? in what country (per Anderson) was the first
More informationExample of Aircraft Climb and Maneuvering Performance. Dr. Antonio A. Trani Professor
Example of Aircraft Climb and Maneuvering Performance CEE 5614 Analysis of Air Transportation Systems Dr. Antonio A. Trani Professor Example - Aircraft Climb Performance Aircraft maneuvering performance
More informationequation 4.1 INTRODUCTION
4 The momentum equation 4.1 INTRODUCTION It is often important to determine the force produced on a solid body by fluid flowing steadily over or through it. For example, there is the force exerted on a
More informationPREDICTION OF SOUND PRESSURE LEVELS ON ROCKET VEHICLES DURING ASCENT Revision E
PREDICTION OF SOUND PRESSURE LEVELS ON ROCKET VEHICLES DURING ASCENT Revision E By Tom Irvine Email: tomirvine@aol.com July 0, 011 Figure 0. Schlieren Photo, Wind Tunnel Test Engineers conducted wind tunnel
More informationDEVELOPMENT OF A COMPRESSED CARBON DIOXIDE PROPULSION UNIT FOR NEAR-TERM MARS SURFACE APPLICATIONS
DEVELOPMENT OF A COMPRESSED CARBON DIOXIDE PROPULSION UNIT FOR NEAR-TERM MARS SURFACE APPLICATIONS Erin Blass Old Dominion University Advisor: Dr. Robert Ash Abstract This work has focused on the development
More informationAugust Progress Report
PATH PREDICTION FOR AN EARTH-BASED DEMONSTRATION BALLOON FLIGHT DANIEL BEYLKIN Mentor: Jerrold Marsden Co-Mentors: Claire Newman and Philip Du Toit August Progress Report. Progress.. Discrete Mechanics
More information