Stability-Constrained Aerodynamic Shape Optimization of a Flying Wing Configuration
|
|
- Bertina Amelia Pierce
- 5 years ago
- Views:
Transcription
1 13th AIAA/ISSMO Multidisciplinary Analysis Optimization Conference September 2010, Fort Worth, Texas AIAA th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference September 13 15, 2010, Fort Worth, Texas Stability-Constrained Aerodynamic Shape Optimization of a Flying Wing Configuration Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim R.R.A. Martins Department of Aerospace Engineering, University of Michigan Ann Arbor, Michigan Stability-constrained aerodynamic shape optimization is conducted using three dimensional, time-spectral computational fluid dynamics. To facilitate this, a method for computing stability constraints with the time-spectral method is introduced. The unique characteristics of this method allow an adjoint technique to be used to generate the gradients of the stability constraints that are required for optimization. These gradients are generated using the automatic differentiation adjoint technique, which provides accurate and efficient results. The resulting stability constraints and their gradients are used in the optimization of a flying wing. I. Introduction Aerodynamic optimization is a rapidly maturing field. Using adjoint methods with Euler and Navier Stokes computational fluid dynamics (CFD), a number of researchers have successfully solved drag minimization problems. 1 5 However, when considering unconventional designs such as a flying wing with no distinct horizontal tail there is a strong coupling between the aerodynamic efficiency of the aircraft and its stability characteristics. This coupling leads to a different optimal solution from a stability-constrained optimization than that which is obtained from a simple aerodynamic optimization. In this work we conduct aerodynamic optimization with basic, static longitudinal stability constraints included. Two approaches are considered. In the first approach, we constrain C m to be approximately zero and C Mα to be negative about the center of gravity (CG) of the aircraft. This combination of parameters ensures that the aerodynamic moments caused by perturbations from the cruise condition are restoring moments, giving a statically stable aircraft. In the second approach, we constrain the value C mα + C Lα Static Margin to be zero. This allows the optimization to return a design with a predefined static margin at the optimal solution. Unlike the force coefficients (C L,C D and C m ), which can be computed directly with any flow solver, stability derivatives, such as C mα and C mq, are themselves derivative quantities and require a special flow solver for their computation. In addition, in order to be able to conduct gradient based optimization, these modifications must be made in a way that allows the gradients of these quantities to be computed in an efficient manner. In this work, the stability derivative, C mα, is calculated using a new method based on a time-spectral CFD solution. The algebraic method used in conjunction with the time spectral solution allows the use of an adjoint method for the derivative computation, which achieves the desired efficiency. The details of this method are provided in Section III. Ph.D Candidate, AIAA Student Member AssociateProfessor, AIAA Senior Member Copyright c 2010 by the authors. Published by the American Institute of Aeronautics and Astronautics, Inc. with permission. 1 of 14 Copyright 2010 by the authors. Published by the American American Institute Institute of Aeronautics of Aeronautics and Astronautics, and Astronautics Inc., with permission. Paper
2 Due to the relatively high computational cost of the CFD analysis and the potentially large number of design variables involved in shape optimization, a gradient-based method is selected for optimization. In this case we use an SQP algorithm (SNOPT 6 ) and an automatic differentiation adjoint method (ADjoint 7 ) to compute the gradients required by the optimizer. These tools are used to demonstrate the effect of static longitudinal stability constraints on an aerodynamic optimization problem. Optimizations have been conducted on cases with and without the stability constraints and the resulting designs are compared to highlight the effects of the stability and trim constraints on the optimum design. II. Time-Spectral Method The underlying CFD algorithm used in this work is a time-spectral Euler solver. The basics of the method are repoduced here. For a more detailed explanation the reader is referred to the original works by van der Weide et al. 8 and Gopinath and Jameson. 9 We start by defining the governing equations for unsteady flow. These can be written as: V w t + f i x i = 0, (1) where x i are the coordinates in the i th direction. For the Euler equations, the state and the fluxes for each cell are: ρ ρu i ρu 1 ρu i u 1 + pδ i1 w = ρu 2, f i = ρu i u 2 + pδ i2. (2) ρu 3 ρu i u 3 + pδ i3 ρe ρu i H To simplify notation, we will rewrite this equation as: V w t + R(w) = 0. (3) To formulate the time-spectral equations, we assume that the flow is periodic in time and express the solution of the states, w, as a Fourier series and its inverse transform: ŵ k = 1 N N 1 n=0 w n e ik 2π T n t, (4) w n = N 2 1 ŵ k e ik 2π T n t, (5) k= N 2 where N is the number of time intervals, T is the period of the solution and t = T/N. We can then evaluate the partial derivative of the states with respect to time as: D t w n = 2π T N 2 1 k= N 2 which can be combined with Equation (4) to form a matrix operator, D t w n = N 1 j=0 ikŵ k e ik 2π T n t, (6) j j nw j, (7) which connects all of the time instances in the solution. simultaneous steady state systems based on the equation: Thus, we end up solving the equivalent of N V D t w n + R(w n ) = 0, (8) 2 of 14
3 where R(w n ) are the steady-state residuals for the n th time instance. For the work in this paper we use the SUmb flow solver 10 that has been developed at Stanford University under the sponsorship of the Department of Energy. SUmb is a cell-centered multiblock solver for the Reynolds-averaged Navier Stokes equations (steady, unsteady, and time-spectral) and it has options for a variety of turbulence models with one, two and four equations. In this work, we limit ourselves to the time-spectral Euler equations. III. Stability Derivatives A key development in this work is the addition of stability constraints to aerodynamic optimization. The computation of stability derivatives using CFD is not a new idea: several techniques have been demonstrated previously, including finite differences, 11 automatic differentiation 12, 13 and analytic sensitivity methods. 14 While these methods have proven to be effective in predicting stability derivatives, they each use well established first derivative methods to compute the stability derivatives. Including these quantities as constraints in an optimization problem requires the computation of their gradients. Since the quantities in question are themselves derivatives, the computation of their gradients would require second derivative information from the CFD solver. This is a computationally expensive proposition and for large numbers of design variables would be prohibitively expensive. The time-spectral approach presented here does not have that limitation. To evaluate stability derivatives from a time-spectral solution, a series of linear regressions are performed on the periodic solution produced by the CFD solver. Because this approach is simple and algebraic, an adjoint technique can be used over the entire process to calculate the gradients for optimization. A similar method for computing the stability derivatives with a periodic steady solution has been demonstrated by Murman, 15 who uses a frequency domain method to produce periodic data for relatively large amplitude oscillations (±2.5 deg.). This data can then be analyzed with the same techniques used to produce stability derivatives from forced oscillation wind tunnel data, which allows the method to take advantage of the large body of knowledge that has been developed in that field. In our case, we work with much smaller amplitude oscillations designed to isolate specific motion variables individually. This allows individual derivatives to be computed directly. The virtual nature of the CFD solutions allows the simulation of these motions, which would be very difficult to reproduce physically. To develop this method, we refer to the section on linear air reactions in Etkin. 16 From this theory we know that for a general motion, the force and moment coefficients of an aircraft for example, the lift coefficient can be approximated as: C L = C L0 + C Lα α + C L α α + C L α α C Lβ β +..., (9) where all motion states α, β, V, p, q, r, E and their associated time derivatives are included. However, if we simplify the motion to consist of a single dynamic state for example α the value of the remaining motion states are zero and they, and their associated derivatives, drop out of the equation. In the case of a pure α motion this leads to: C L = C L0 + C Lα α + C L α α + C L α α +... (10) If we assume that the higher order derivatives are small and neglect them, we get: C L = C L0 + C Lα α + C L α α. (11) Diagrams illustrating pure α and q motions are shown in Figures 1 and 2, respectively. In these diagrams, the grid is moving with respect to a stationary air mass and the arrows in the diagrams indicate the relative motion of the incoming flow with respect to the grid. The α motion, shown in Figure 1, is an oscillating vertical motion with the forward grid velocity modified in sync with the oscillations such that the incoming velocity maintains a constant magnitude, but changes direction. Time instances 2 and 3 are at the top and bottom of the oscillations, respectively, illustrating the points in the cycle where the grid experiences pure forward motion. Time instance 1 illustrates a point midway through the oscillation where the grid is experiencing both forward and vertical motion. The diagram in Figure 2 illustrates a pure q motion. In this case, the grid rotates with an oscillating value of q and the magnitude and direction of the velocity are modified such that the incoming α is constant. 3 of 14
4 Figure 1. Pure angle of attack motion Figure 2. Pure pitching motion When one of these pure motions is simulated with the time-spectral method, the result is a solution that consists of force coefficient (C L,C D ) and moment coefficient (C mx,c my,c mz ) values at N time instances in a periodic solution. A set of solutions corresponding to a pure q motion are shown in Figure 3. To compute stability derivatives from this solution, we compute a linear least squares fit of the output coefficient (e.g., C L ) with respect to the primary motion variable from the solution (e.g., q). The slope of the line resulting from this fit is the stability derivative (C Lq ) while the y-intercept of the line is the value of C L at q = 0. To validate this approach, derivatives for the NACA 0012 airfoil are computed using the time-spectral method and compared to the numeric results published by Limache et al. 14 These reference results are based on computations from an inviscid two-dimensional CFD solver. The solver was modified to include a non-inertial reference frame formulation, allowing the dynamic stability derivatives to be computed from a steady solution. In the reference work, those derivatives were computed using an analytic, direct sensitivity method. The results of these comparisons are shown in Table 1. Figure 3 corresponds to the Mach = 0.5 case in Table 1. Table 1. Time-spectral stability derivatives for a NACA 0012 airfoil at α = 0.0 degrees Derivative Time-Spectral (3 Point) Reference 14 % Difference Mach = 0.1 C Lq % C Mq % Mach = 0.5 C Lq % C Mq % Mach = 0.8 C Lq % C Mq % As is demonstrated in the table above, the derivatives computed with the time spectral results agree well with the reference results: for the Mach = 0.5 case the results match with a relative difference that is less 4 of 14
5 1.0e e e-03 3 point Solution 5 point Solution 7 point Solution Derivative estimate Reference Estimate 4.0e e-03 C L 0.0e e e e e e e e e e e e e e-04 q (rad/s) Figure 3. Time spectral q derivative: Mach = 0.5 than one percent. The maximum difference occurs for Mach = 0.8, where C mq differs by almost 11%. A possible reason for the larger discrepancy in the Mach = 0.8 case is the existence of a shock in the solution. A slight variation in the predicted location of the shock between the two solutions would cause a significant variation in the prediction of C m and any other values derived from it. IV. Time-Spectral ADjoint Method The time-spectral derivatives are computed directly from the coefficient values calculated in the flow solution. Therefore, an adjoint approach applied to the solver is sufficient to compute the gradients of the coefficients that are required for the optimization. We start by writing the vector-valued function, I, that we want to differentiate as: I = I(x, w n (x)), (12) where x represents the vector of design variables and w n is the state variable vector for the n th time instance, where n = 1,..., N. For a given vector x, the solution of the governing equations of the system yields a vector w n, thus establishing the dependence of the state of the system on the design variables. In general, we denote these governing equations by: R (x, w n (x)) = 0. (13) However, since the flow solution is time-spectral in nature, from van der Weide et al., 8 we define: R T S = V D t w n + R(w n ) = 0, (14) where R(w n ) is a normal steady-state residual for the n th time instance,where n = 1,..., N, and D t is a spectral operator that spans all N time instances. Now we have a modified residual: R T S (x, w n (x)) = 0, (15) 5 of 14
6 which can be treated in the same fashion as the steady state residual would be in a normal adjoint formulation. We use the chain rule to write the total sensitivity of the vector-valued function I as: di dx = I x + I dw n w n dx. (16) Further, because the governing equations must always be satisfied at a converged solution, the total derivative of the residuals (15) with respect to any design variable must also be zero. Expanding the total derivative of the governing equations with respect to the design variables we can write: dr T S dx = R T S x + R T S w n dw n dx = 0. (17) This expression provides the means of eliminating the total derivative dw n / dx from the total sensitivity computation for I. Moving the first term of this equation to the right hand side, we get: R T S w n dw n dx Substituting the solution of this system into Equation (16) yields: = R T S x. (18) di dx = I x I [ ] 1 RT S R T S w n w n x. (19) The adjoint approach consists in factorizing the R T S / w n matrix with the term to its left, yielding the adjoint system: [ ] T RT S ψ = I w w. (20) Then, this solution is used in Equation (19) to obtain the total sensitivity: di dx DV = I x DV ψ T R T S x DV. (21) This yields a sensitivity method with a computational cost that is essentially independent of the number of design variables. More details on the ADjoint method can be found in Mader et al. 7 Note that because the time-spectral system is N times the size of the steady-state solutions, the adjoint system is also N times larger than the equivalent steady-state system. To improve the accuracy and reduce the implementation time of the adjoint solver, we use automatic differentiation to compute the values of the partial derivatives. More specifically, we use a reverse mode, source transformation tool (Tapenade 17, 18 ) to generate derivative code for each of the partial derivative terms. This reverse mode approach is very efficient for the objective function partial derivatives, as the number of objectives is far less than the number of flow states and design variables. The efficiency of the reverse mode approach is less apparent for the residual partial derivatives. However, as described in Mader et al. 7 a modified, single cell version of the residual calculation is used, which makes the reverse mode computation extremely efficient for this calculation. To facilitate the solution of the large, sparse linear systems represented by these equations, we use the Portable, Extensible Toolkit for Scientific Computation 19, 20 (PETSc). This is a sparse linear and non-linear solution package that provides a variety of parallel solution and preconditioning algorithms. V. Center of Gravity Calculation Both the moment coefficient, C m, and all of its derivatives in this case C mα are strongly dependent on the center of gravity (CG) location. In order to facilitate the optimizations described in Section C, a simple, wing center of gravity calculation was implemented. The method used is derived from the work of Chai et al. 21 They state that the wing center of gravity for a normal transport wing is located between the fore and aft spar along the wing mean aerodynamic chord (MAC). That idea is used as the basis for the method described below. We start by calculating the MAC and the location of its quarter chord. These calculations are based on the methods presented in ESDU item This provides a basis on which to 6 of 14
7 locate the fore and aft spars for the wing. As part of the problem definition, the locations of the spars are defined. In this work we have assume that the forward spar is at 25% of the MAC and the aft spar is at 75% of the MAC. We can then specify the location of the CG as a percentage of the distance between the spars (CG %MAC ). This gives a reference location about which the moment coefficient, C m, and all of its derivatives can be calculated. VI. Results The results are presented in four sections. Section A describes the test case used and presents some basic solutions for the case. Section B presents the results relating to the performance and accuracy of the timespectral adjoint method used in this work, and Section C presents the results relating to the optimizations conducted on flying wings. Comparisons between the various optimization results are discussed in Section D A. Test Case The test case used in this work is the ONERA M6 wing, commonly used in CFD. The test case is run at Mach = 0.8 with a symmetry condition at the root and inviscid wall boundary conditions. The solver used is SUmb, 10 a three dimensional, structured, multi-block solver with a cell-centered, finite-volume scheme. The grid contains 155,648 cells and is shown in Figure 4(a). A steady state solution is shown in Figure 4(b) and a time-spectral solution is shown in Figure 5. (a) ONERA M6: mesh Figure 4. (b) ONERA M6: surface C p distribution: Mach = 0.8 α=3.0 deg. ONERA M6 test case (a) Time instance 1 (b) Time instance 2 (c) Time instance 3 Figure 5. ONERA M6: time-spectral surface C p distribution, Mach=0.8, α=3.0 deg. 7 of 14
8 B. Performance To demonstrate the performance of the methods developed, accuracy and timing results are shown below. Table 2. Sensitivity verification: Onera M6 test case 155,648 cells: convergence Cost Function Design Variable Adjoint Complex Step C L α C D α E E-2 C m α E E-1 C mα α C L Sweep E E-5 C D Sweep E E-4 C m Sweep E E-4 C mα Sweep E E-1 C L Twist E E-2 C D Twist E E-4 C m Twist E E-3 C mα Twist E E-2 C L Taper E E-3 C D Taper E E-4 C m Taper E E-5 C mα Taper E E-2 The derivatives shown in Table 2 clearly demonstrate the accuracy of the derivative method used. In most cases the derivatives are accurate to at least 8 digits, with the worst case being a 6 digit disagreement and the best case 10 digit agreement. This is extremely good agreement, especially when considering that the solution is converged to a relative level of 10 10, which means that 10 digit agreement is the best that can be expected from an numerically exact derivative method. Steady-State Onera M6 Case Time-Spectral Number of processors Number of cells 155, ,648 Number of time instances 1 3 Flow solution ADjoint Breakdown: Computation of residual matrices Computation of reconditioning matrix Computation of volume to surface partial derivative Computation of RHS Solution of the adjoint equations Computation of total sensitivities Table 3. ADjoint computational cost breakdown (times in seconds) The timing results shown in Table 3 demonstrate the efficiency of the approach and implementation of the ADjoint method. Comparing the ADjoint solvers to the flow solvers, in both cases the ADjoint takes less than half as long as the flow solver (0.4 and 0.44 times respectively). Given that this is a reverse derivative that computes the derivative of each output with respect to all of the input values, this is even 8 of 14
9 more attractive for problems with large numbers of design variables. To illustrate, consider a simple timespectral optimization problem with 50 design variables and two output variables, one objective and one constraint. A finite difference method would require 50 additional flow solutions 11.8 hours for a single gradient evaluation, while the adjoint method would require 2 adjoint solutions 13 minutes for the same computation. The cost advantage of using the adjoint approach is evident. Further, the cost of using the time-spectral computation instead of a steady computation is also reasonable. The time-spectral flow solver is 3.5 times as costly as the steady solver and the time-spectral adjoint is 3.8 times as costly as the steady adjoint. While a factor of four increase in cost is significant, it is a reasonable cost to pay to be able to include stability characteristics in the optimization process. C. Optimization To assess the effects of stability constraints on aerodynamic shape optimization results, four optimization problems are examined. To serve as a baseline case, the following lift-constrained drag minimization problem is solved: minimize C D (x, w n (x)) w.r.t. x = [α, sweep, twist, taper] (22) subject to R(x, w n (x)) = 0 C Ltarget C L (x, w n (x)) 0 As an intermediate case, a trim constraint is added to the problem by limiting the moment coefficient,i.e., 0 C m This gives: minimize C D (x, w n (x)) w.r.t. x = [α, sweep, twist, taper, x CGtarget, CG %MAC ] (23) subject to R(x, w n (x)) = 0 C Ltarget C L (x, w n (x)) 0 0 C m (x, w n (x)) x CG x CGtarget 10 5 (24) For the first stability constrained problem, the C mα constraint is added, giving: minimize C D (x, w n (x)) w.r.t. x = [α, sweep, twist, taper, x CGtarget, CG %MAC ] (25) subject to R(x, w n (x)) = 0 C Ltarget C L (x, w n (x)) 0 0 C m (x, w n (x)) 10 4 C mα (x, w n (x)) x CG x CGtarget 10 5 (26) For the second stability constrained problem, the static margin constraint replaces the C mα constraint, giving: minimize C D (x, w n (x)) w.r.t. x = [α, sweep, twist, taper, x CGtarget, CG %MAC ] (27) subject to R(x, w n (x)) = 0 C Ltarget C L (x, w n (x)) 0 0 C m (x, w n (x)) 10 4 C mα (x, w n (x)) + C Lα (x, w n (x)) Static Margin(x) = x CG x CGtarget 10 5 (28) 9 of 14
10 The last three problems are formulated as multidisciplinary systems. The CG calculation forms one discipline and the aerodynamics forms the second discipline. An interdisciplinary feasible (IDF) approach is used to solve the system. As shown in the problem definitions, the discipline coupling variable is the CG location. A single target variable, x CGtarget, is added as an input to the aerodynamic discipline and an additional constraint that x CG x CGtarget = 0 is added to the problem to enforce consistency at the optimized solution. In addition, because the CG location and neutral point location are both primarily driven by the sweep variable, the optimizer has been given additional freedom to specify the CG location as a percentage of the distance between the spars (CG %MAC ). This variable was allowed to vary anywhere in between the spars as well as one MAC length in front of the forward spar. The extra forward freedom is allowed to take into account the potential CG contributions from non-wing portions of the aircraft, such as systems, payload, etc.. Each of these optimizations has been run at Mach= 0.8 with a target lift coefficient of 0.3. The reference temperature, density and pressure are K, kg/m 3 and Pa respectively. For the static margin constrained optimization, the requested static margin is 5% of the MAC. Table 4 outlines the limits specified on the various design variables used. The results of the three optimizations are Design Variable Bounds Design Variable Lower Bound Upper Bound α(deg.) Sweep (deg.) Tip Twist(deg.) Taper Ratio x CGtarget (m) CG %MAC -1 1 Table 4. Design variable limits shown in Table 5. Note that in the trim constrained and stability constrained optimizations, the target center of gravity location, x CGtarget, was used as the reference center for the computation of the moments and their derivatives. Optimization Cases Parameter ONERA M6 Lift Constrained Trimmed C mα Constrained Static Margin Constrained C L C D C m C mα CG Location(m) α(deg.) Sweep (deg.) Tip Twist(deg.) Taper Ratio Table 5. Optimization results comparison A convergence plot for the five cases shown in Table 5 is shown in Figure 6 D. Discussion As shown in Table 5, all of the optimizations achieved the requested target lift coefficient of 0.3 and have reduced the drag with respect to the baseline case. The lowest drag was achieved in both the lift constrained case and the trim constrained optimizations. The optimizer was able to match the performance of the lift constrained case with the trim constrained case by moving the CG variable to trim the wing instead of altering the aerodynamic shape. As expected, the predicted drag increases with both of the stability constrained 10 of 14
11 Merit Function Convergence ONERA M6 Baseline Lift Constraint Trim Constraint C mα Constraint Static Margin Constraint Iterations Figure 6. Convergence histories for the various optimization problems cases. The penalty of the static margin constraint is larger than the C mα constraint because it forces a more negative value of C mα than the default C mα constraint. Looking at the final shape of the design, in all four cases the bulk of the drag minimization comes from the increase in sweep. This change reduces the wave drag caused by the shock on the wing in the baseline case. The trim and stability constraints are satisfied with a combination of changes to the other variables. The most prominent changes appear to be in twist and CG location. Both of the stability constrained cases have more twist than the other two. The static margin case has significantly more twist than any of the others and the CG is significantly farther forward. The resulting flow solutions for each of the cases are plotted in Figures 7(a)-7(e) and a comparison of the resulting lift distributions is shown in Figure 8. The lift distributions also show what one would expect. In order to counter the typical nose down moment of a symmetric airfoil, the optimizer has shifted lift from the tips inboard. This lift reduction at the tips creates a nose up moment that trims and stabilizes the wing. Note that the baseline case and the trimmed case have the same lift distribution. These two cases have the same aerodynamic shape, but the CG has been moved forward in the trimmed optimization to counteract the aerodynamic moment. Also note that the optimal lift constrained solution is not elliptical. This is primarily because the twist is forced to maintain a linear distribution along the span. VII. Conclusion Stability-constrained optimization is conducted for a flying wing using a time-spectral flow solver and adjoint method. A new time-spectral approach to the computation of stability derivatives is shown to be valid and an ADjoint implementation of the time-spectral adjoint is shown to be computationally efficient and to have accuracy greater than any previous implementation of a time-spectral adjoint. The optimizations conducted with this framework produced sensible and meaningful results. As expected, the optimum stability constrained result has a higher drag coefficient than the baseline lift-constrained case. The trimmed case is similar to the baseline case. The differences evident in the stability constrained case reinforce the fact that for some unconventional designs, such as flying wings, including stability considerations in the design optimization is important. 11 of 14
12 (a) ONERA M6 solution (b) Lift constrained solution (c) Trim constrained solution (d) C mα constrained solution (e) Static margin constrained solution Figure 7. Optimal Onera M6 results 12 of 14
13 C l c ONERA M6 Lift Constraint Trim Constraint C mα Constraint Static Margin Constraint Elliptical Span Location Figure 8. Lift distribution comparison Acknowledgements The authors would like to thank Dr. Edwin van der Weide for his help with the SUmb flow solver. Financial support for this research was provided by the University of Toronto and the Natural Sciences and Engineering Research Council(NSERC) of Canada. Computations were performed on the GPC supercomputer at the SciNet HPC Consortium. SciNet is funded by: the Canada Foundation for Innovation under the auspices of Compute Canada; the Government of Ontario; Ontario Research Fund Research Excellence; and the University of Toronto. 13 of 14
14 References 1 Jameson, A., Aerodynamic Design via Control Theory, Journal of Scientific Computing, Vol. 3, No. 3, sep 1988, pp Nemec, M. and ingg, D. W., Newton Krylov Algorithm for Aerodynamic Design Using the Navier Stokes Equations, AIAA Journal, Vol. 40, No. 6, June 2002, pp Hicken, J. and ingg, D. W., Aerodynamic Optimization Algorithm with Integrated Geometry Parameterization and Mesh Movement, AIAA Journal, Vol. 48, No. 2, February 2010, pp Reuther, J. J., Jameson, A., Alonso, J. J.,, Rimlinger, M. J., and Saunders, D., Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers, Part 1, Journal of Aircraft, Vol. 36, No. 1, 1999, pp Reuther, J. J., Jameson, A., Alonso, J. J.,, Rimlinger, M. J., and Saunders, D., Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers, Part 2, Journal of Aircraft, Vol. 36, No. 1, 1999, pp Gill, P. E., Murray, W., and Saunders, M. A., SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization, SIAM Review, Vol. 47, No. 1, 2005, pp Mader, C. A., Martins, J. R. R. A., Alonso, J. J., and van der Weide, E., ADjoint: An Approach for the Rapid Development of Discrete Adjoint Solvers, AIAA Journal, Vol. 46, No. 4, 2008, pp van der Weide, E., Gopinath, A. K., and Jameson, A., Turbomachinery Applications with the Time Spectral Method, AIAA Paper , Gopinath, A. K. and Jameson, A., Time Spectral Method for Periodic Unsteady Computations over Two- Three Dimensional Bodies, AIAA Paper , van der Weide, E., Kalitzin, G., Schluter, J., and Alonso, J. J., Unsteady Turbomachinery Computations Using Massively Parallel Platforms, Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Babcock, D. A. and Jr., D. A. S. A., Estimating Aircraft Stability Derivatives Through Finite Element Analysis, AIAA Paper , Aug Park, M. A., Green, L. L., Montgomery, R. C., and Raney, D. L., Determination of Stability and Control Derivatives Using Computational Fluid Dynamics and Automatic Differentiation, AIAA Paper , Park, M. A. and Green, L. L., Steady-State Computation of Constant Rotational Rate Dynamic Stability Derivatives, AIAA Paper , Limache, A. and Cliff, E., Aerodynamic Sensitivity Theory for Rotary Stability Derivatives, Journal of Aircraft, Vol. 37, 2000, pp Murman, S. M., Reduced-Frequency Approach for Calculating Dynamic Derivatives, AIAA Journal, Vol. 45, No. 6, 2007, pp Etkin, B., Dynamics of Atmospheric Flight, Dover Publications, Mineola, New ork, Pascual, V. and Hascoët, L., Extension of TAPENADE Towards Fortran 95, Automatic Differentiation: Applications, Theory, and Tools, edited by H. M. Bücker, G. Corliss, P. Hovland, U. Naumann, and B. Norris, Lecture Notes in Computational Science and Engineering, Springer, Hascoët, L. and Pascual, V., TAPENADE 2.1 User s Guide, Technical report 300, INRIA, Balay, S., Buschelman, K., Eijkhout, V., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and hang, H., PETSc Users Manual, Tech. Rep. ANL-95/11 - Revision 2.1.5, Argonne National Laboratory, Chai, S., Crisafulli, P., and W.H.Mason, Aircraft Center of Gravity Estimation in Conceptual/Preliminary Design, Proceedings of the 1st AIAA Aircraft Engineering, Technology, and Operations Congress, Los Angeles, CA, 1995, AIAA ESDU, Geometrical Properties of Cranked and Straight Tapered Wing Planforms, Tech. Rep. Item No , ESDU International plc, London, UK, January of 14
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim
More informationChallenges and Complexity of Aerodynamic Wing Design
Chapter 1 Challenges and Complexity of Aerodynamic Wing Design Kasidit Leoviriyakit and Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA kasidit@stanford.edu and
More informationComputing Stability Derivatives and their Gradients for Aerodynamic Shape Optimization
Computing Stability Derivatives and their Gradients for Aerodynamic Shape Optimization Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim R.R.A. Martins
More informationAERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS
AERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS Joaquim R. R. A. Martins and Juan J. Alonso Department of Aeronautics and Astronautics Stanford University, Stanford, CA
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationCOMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD
COMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD Joaquim R. R. A. Martins Juan J. Alonso James J. Reuther Department of Aeronautics and Astronautics Stanford
More informationFLIGHT dynamic characteristics are an essential factor in aircraft
AIAA JOURNAL Vol. 49, No. 12, December 2011 Computation of Aircraft Stability Derivatives Using an Automatic Differentiation Adjoint Approach Charles A. Mader University of Toronto, Toronto, Ontario M3H
More informationToward Practical Aerodynamic Design Through Numerical Optimization
Toward Practical Aerodynamic Design Through Numerical Optimization David W. Zingg, and Laura Billing, Institute for Aerospace Studies, University of Toronto 4925 Dufferin St., Toronto, Ontario M3H 5T6,
More informationFar Field Noise Minimization Using an Adjoint Approach
Far Field Noise Minimization Using an Adjoint Approach Markus P. Rumpfkeil and David W. Zingg University of Toronto Institute for Aerospace Studies 4925 Dufferin Street, Toronto, Ontario, M3H 5T6, Canada
More informationActive Flutter Control using an Adjoint Method
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
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 informationMultidisciplinary Design Optimization of Aerospace Systems
Please cite this document as: Martins, Joaquim R. R. A., Multidisciplinary Design Optimization of Aerospace Systems, Advances and Trends in Optimization with Engineering Applications, SIAM, 2017, pp. 249
More informationJun 22-25, 2009/San Antonio, TX
19th AIAA Computational Fluid Dynamics 22-25 June 2009, San Antonio, Texas AIAA 2009-4273 AIAA 2009 4273 An Assessment of Dual-Time Stepping, Time Spectral and Artificial Compressibility based Numerical
More informationNUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM)
Journal of Fundamental and Applied Sciences ISSN 1112-9867 Available online at http://www.jfas.info NUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM) M. Y. Habib
More informationA HARMONIC BALANCE APPROACH FOR MODELING THREE-DIMENSIONAL NONLINEAR UNSTEADY AERODYNAMICS AND AEROELASTICITY
' - ' Proceedings of ASME International Mechanical Engineering Conference and Exposition November 17-22, 22, New Orleans, Louisiana, USA IMECE-22-3232 A HARMONIC ALANCE APPROACH FOR MODELING THREE-DIMENSIONAL
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037
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 informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More informationA COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION
A COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION Joaquim Rafael Rost Ávila Martins Department of Aeronautics and Astronautics Stanford University Ph.D. Oral Examination, Stanford
More informationScalable parallel approach for high-fidelity steady-state aeroelastic analysis and adjoint derivative computations
This is a preprint of the following article, which is available from http://mdolab.engin.umich.edu Gaetan K. W. Kenway, Graeme J. Kennedy, and Joaquim R. R. A. Martins, Scalable Parallel Approach for High-
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
More informationApril 15, 2011 Sample Quiz and Exam Questions D. A. Caughey Page 1 of 9
April 15, 2011 Sample Quiz Exam Questions D. A. Caughey Page 1 of 9 These pages include virtually all Quiz, Midterm, Final Examination questions I have used in M&AE 5070 over the years. Note that some
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationEfficient Hessian Calculation using Automatic Differentiation
25th AIAA Applied Aerodynamics Conference, 25-28 June 2007, Miami Efficient Hessian Calculation using Automatic Differentiation Devendra P. Ghate and Michael B. Giles Oxford University Computing Laboratory,
More informationPRINCIPLES OF FLIGHT
1 Considering a positive cambered aerofoil, the pitching moment when Cl=0 is: A infinite B positive (nose-up). C negative (nose-down). D equal to zero. 2 The angle between the aeroplane longitudinal axis
More informationStudy of Preliminary Configuration Design of F-35 using simple CFD
Study of Preliminary Configuration Design of F-35 using simple CFD http://www.aerospaceweb.org/aircraft/research/x35/pics.shtml David Hall Sangeon Chun David Andrews Center of Gravity Estimation.5873 Conventional
More informationAERODYNAMIC COEFFICIENTS FOR EXTENDING AND BENDING PROJECTILES. William G. Reinecke*
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 AERODYNAMIC COEFFICIENTS FOR EXTENDING AND BENDING PROJECTILES William G. Reinecke* Institute for Advanced Technology.The University
More informationAIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION
AIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION AIAA-2004-0689 Serhat Hosder, Joseph A. Schetz, Bernard Grossman and William H. Mason Virginia Tech Work sponsored by NASA
More informationAircraft Performance, Stability and control with experiments in Flight. Questions
Aircraft Performance, Stability and control with experiments in Flight Questions Q. If only the elevator size of a given aircraft is decreased; keeping horizontal tail area unchanged; then the aircraft
More informationAERODYNAMIC CHARACTERIZATION OF A CANARD GUIDED ARTILLERY PROJECTILE
45th AIAA Aerospace Sciences Meeting and Exhibit 8-11 January 27, Reno, Nevada AIAA 27-672 AERODYNAMIC CHARACTERIZATION OF A CANARD GUIDED ARTILLERY PROJECTILE Wei-Jen Su 1, Curtis Wilson 2, Tony Farina
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More informationWind Tunnel Study of a Large Aerostat, CFD Validation
AIAA Lighter-Than-Air Systems Technology (LTA) Conference 25-28 March 2013, Daytona Beach, Florida AIAA 2013-1339 Wind Tunnel Study of a Large Aerostat, CFD Validation Stephen C. Chan 1, Kaleb Shervington
More informationInvestigation potential flow about swept back wing using panel method
INTERNATIONAL JOURNAL OF ENERGY AND ENVIRONMENT Volume 7, Issue 4, 2016 pp.317-326 Journal homepage: www.ijee.ieefoundation.org Investigation potential flow about swept back wing using panel method Wakkas
More informationAPPENDIX C DRAG POLAR, STABILITY DERIVATIVES AND CHARACTERISTIC ROOTS OF A JET AIRPLANE (Lectures 37 to 40)
APPENDIX C DRAG POLAR, STABILITY DERIVATIVES AND CHARACTERISTIC ROOTS OF A JET AIRPLANE (Lectures 37 to 40 E.G. TULAPURKARA YASHKUMAR A. VENKATTRAMAN REPORT NO: AE TR 2007-3 APRIL 2007 (REVISED NOVEMBER
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 information41st AIAA Aerospace Sciences Meeting and Exhibit January 6 9, 2003/Reno, NV
AIAA 3 267 Calculation of Unsteady Transonic Flow by an Euler Method with Small Perturbation Boundary Conditions C. Gao, S. Luo, F. Liu Department of Mechanical and Aerospace Engineering University of
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationLimit Cycle Oscillations of a Typical Airfoil in Transonic Flow
Limit Cycle Oscillations of a Typical Airfoil in Transonic Flow Denis B. Kholodar, United States Air Force Academy, Colorado Springs, CO 88 Earl H. Dowell, Jeffrey P. Thomas, and Kenneth C. Hall Duke University,
More informationExtended longitudinal stability theory at low Re - Application to sailplane models
Extended longitudinal stability theory at low Re - Application to sailplane models matthieu.scherrer@free.fr November 26 C L C m C m W X α NP W X V NP W Lift coefficient Pitching moment coefficient Pitching
More informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
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 informationNUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE
NUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE Wessam Mahfouz Elnaggar, Zhihua Chen and Hui Zhang Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing,
More informationAero-Propulsive-Elastic Modeling Using OpenVSP
Aero-Propulsive-Elastic Modeling Using OpenVSP August 8, 213 Kevin W. Reynolds Intelligent Systems Division, Code TI NASA Ames Research Center Our Introduction To OpenVSP Overview! Motivation and Background!
More informationSome effects of large blade deflections on aeroelastic stability
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-839 Some effects of large blade deflections on aeroelastic stability
More informationIntroduction to Aerospace Engineering
4. Basic Fluid (Aero) Dynamics Introduction to Aerospace Engineering Here, we will try and look at a few basic ideas from the complicated field of fluid dynamics. The general area includes studies of incompressible,
More informationGiven the water behaves as shown above, which direction will the cylinder rotate?
water stream fixed but free to rotate Given the water behaves as shown above, which direction will the cylinder rotate? ) Clockwise 2) Counter-clockwise 3) Not enough information F y U 0 U F x V=0 V=0
More informationINFLUENCE OF GUST MODELLING ON FREE-FLIGHT AEROFOILS
International Forum on Aeroelasticity and Structural Dynamics IFASD 2017 25-28 June 2017 Como, Italy INFLUENCE OF GUST MODELLING ON FREE-FLIGHT AEROFOILS Reik Thormann 1, Giampaolo Pagliuca 1, Sebastian
More informationComparative Analysis of Computational Methods for Limit-Cycle Oscillations
47th AIAA Structures, Structural Dynamics and Materials Conference, May 1-4, 26 Newport, Rhode Island Comparative Analysis of Computational Methods for Limit-Cycle Oscillations Arathi K. Gopinath, Stanford
More informationCFD COMPUTATION OF THE GROUND EFFECT ON AIRPLANE WITH HIGH ASPECT RATIO WING
28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES CFD COMPUTATION OF THE GROUND EFFECT ON AIRPLANE WITH HIGH ASPECT RATIO WING Sun Tae Kim*, Youngtae Kim**, Tae Kyu Reu* *Agency for Defense Development,
More informationA New Approach for the Estimation of Longitudinal Damping Derivatives: CFD Validation on NACA 0012
A New Approach for the Estimation of Longitudinal Damping Derivatives: CFD Validation on NACA PIERO GILI Polytechnic of Turin DIMEAS C.so Duca degli Abruzzi 4 piero.gili@polito.it MICHELE VISONE Blue Engineering
More informationAERO-STRUCTURAL MDO OF A SPAR-TANK-TYPE WING FOR AIRCRAFT CONCEPTUAL DESIGN
1 26th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AERO-STRUCTURAL MDO OF A SPAR-TANK-TYPE WING FOR AIRCRAFT CONCEPTUAL DESIGN Paulo R. Caixeta Jr., Álvaro M. Abdalla, Flávio D. Marques, Fernando
More informationApplications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow
Applications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow Antony Jameson and Kui Ou Aeronautics & Astronautics Department,
More informationAeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior
Aeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior N. Manoj Abstract The aeroelastic behavior of engine nacelle strake when subjected to unsteady aerodynamic flows is
More informationSENSITIVITY ANALYSIS OF THE FACTORS AFFECTING FORCE GENERATION BY WING FLAPPING MOTION
Proceedings of the ASME 2013 International Mechanical Engineering Congress and Exposition IMECE2013 November 15-21, 2013, San Diego, California, USA IMECE2013-65472 SENSITIVITY ANALYSIS OF THE FACTORS
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationComputational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2
Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 1 M. Tech Scholar, 2 Associate Professor Department of Mechanical Engineering, Bipin Tripathi
More informationComplete Configuration Aero-Structural Optimization Using a Coupled Sensitivity Analysis Method
AIAA 2002 5402 Complete Configuration Aero-Structural Optimization Using a Coupled Sensitivity Analysis Method Joaquim R. R. A. Martins, Juan J. Alonso Stanford University, Stanford, CA 94305 James J.
More informationA STUDY ON AERODYNAMIC DESIGN OF HIGH- EFFICIENCY PROPELLERS AT HIGH-ALTITUDE
A STUDY ON AERODYNAMIC DESIGN OF HIGH- EFFICIENCY PROPELLERS AT HIGH-ALTITUDE Jianhua Xu, Wenping Song & Wenqing Yang National Key Laboratory of Aerodynamic Design and Research, Northwestern Polytechnical
More informationAeroelastic Gust Response
Aeroelastic Gust Response Civil Transport Aircraft - xxx Presented By: Fausto Gill Di Vincenzo 04-06-2012 What is Aeroelasticity? Aeroelasticity studies the effect of aerodynamic loads on flexible structures,
More informationTransonic Flutter Prediction of Supersonic Jet Trainer with Various External Store Configurations
Transonic Flutter Prediction of Supersonic Jet Trainer with Various External Store Configurations In Lee * Korea Advanced Institute of Science and Technology, Daejeon, 305-701, Korea Hyuk-Jun Kwon Agency
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
45 th Aerospace Sciences Meeting and Exhibit, January 8 11, 2007, Reno, Nevada An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics
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 informationMechanics of Flight. Warren F. Phillips. John Wiley & Sons, Inc. Professor Mechanical and Aerospace Engineering Utah State University WILEY
Mechanics of Flight Warren F. Phillips Professor Mechanical and Aerospace Engineering Utah State University WILEY John Wiley & Sons, Inc. CONTENTS Preface Acknowledgments xi xiii 1. Overview of Aerodynamics
More informationConsider a wing of finite span with an elliptic circulation distribution:
Question 1 (a) onsider a wing of finite span with an elliptic circulation distribution: Γ( y) Γo y + b = 1, - s y s where s=b/ denotes the wing semi-span. Use this equation, in conjunction with the Kutta-Joukowsky
More informationSupplementary Section D: Additional Material Relating to Helicopter Flight Mechanics Models for the Case Study of Chapter 10.
Supplementary Section D: Additional Material Relating to Helicopter Flight Mechanics Models for the Case Study of Chapter 1. D1 Nonlinear Flight-Mechanics Models and their Linearisation D1.1 Introduction
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationCHAPTER 3 ANALYSIS OF NACA 4 SERIES AIRFOILS
54 CHAPTER 3 ANALYSIS OF NACA 4 SERIES AIRFOILS The baseline characteristics and analysis of NACA 4 series airfoils are presented in this chapter in detail. The correlations for coefficient of lift and
More informationA hybrid adjoint approach applied to RANS equations
Center for urbulence Research Annual Research Briefs 2013 303 A hybrid adjoint approach applied to RANS equations By. W. R. aylor, F. Palacios, K. Duraisamy AND J. J. Alonso 1. Motivation and objectives
More informationTurbomachinery Applications with the Time Spectral Method
35th AIAA Fluid Dynamics Conference and Exhibit, June 6-9,005, Toronto,Ontario Turbomachinery Applications with the Time Spectral Method Edwin van der Weide, Arathi K Gopinath and Antony Jameson Stanford
More informationEVALUATION OF TRANSONIC BUFFETING ONSET BOUNDARY ESTIMATED BY TRAILING EDGE PRESSURE DIVERGENCE AND RMS DATA OF WING VIBRATION
EVALUATION OF TRANSONIC BUFFETING ONSET BOUNDARY ESTIMATED BY TRAILING EDGE PRESSURE DIVERGENCE AND RMS DATA OF WING VIBRATION Rodrigo Sorbilli C. de Sousa, Roberto da Motta Girardi Roberto Gil Annes da
More informationHigh-Fidelity Multidisciplinary Design Using an Integrated Design Environment
High-Fidelity Multidisciplinary Design Using an Integrated Design Environment Antony Jameson and Juan J. Alonso Department of Aeronautics and Astronautics Stanford University, Stanford CA AFOSR Joint Contractors
More informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
More informationEvaluation of the Drag Reduction Potential and Static Stability Changes of C-130 Aft Body Strakes
U.S. Air Force T&E Days 2009 10-12 February 2009, Albuquerque, New Mexico AIAA 2009-1721 Evaluation of the Drag Reduction Potential and Static Stability Changes of C-130 Aft Body Strakes Heather G. Pinsky
More informationOptimum aerodynamic design using the Navier-Stokes equations
Copyright 1997, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1997 A9715189, AF-AFOSR-91-0391, AIAA Paper 97-0101 Optimum aerodynamic design using the Navier-Stokes
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics Lecture 7 Equations of motion Mark Voskuijl Semester 1-2012 Delft University of Technology Challenge the
More informationPrediction of Aircraft s Longitudinal Motion Based on Aerodynamic Coefficients and Derivatives by Surrogate Model Approach
Journal of Mechanics Engineering and Automation 4 (2014) 584-594 D DAVID PUBLISHING Prediction of Aircraft s Longitudinal Motion Based on Aerodynamic Coefficients and Derivatives by Surrogate Model Approach
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationComparative Analysis of Computational Methods for Limit-Cycle Oscillations
47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, May 1-4, 26 Newport, Rhode Island Comparative Analysis of Computational Methods for Limit-Cycle Oscillations Arathi
More informationAE Stability and Control of Aerospace Vehicles
AE 430 - Stability and ontrol of Aerospace Vehicles Static/Dynamic Stability Longitudinal Static Stability Static Stability We begin ith the concept of Equilibrium (Trim). Equilibrium is a state of an
More informationSolution Algorithms for Viscous Flow
Solution Algorithms for Viscous Flow Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA Indian Institute of Science Bangalore, India September 20, 2004 Flo3xx Computational
More informationRoyal Aeronautical Society 2016 Applied Aerodynamics Conference Tuesday 19 th Thursday 21 st July Science Centre, Bristol, UK
Assessment and validation of aerodynamic performance results for a Natural Laminar Flow transonic wing tested in cryogenic conditions via simulation of turbulent wedges in CFD. Royal Aeronautical Society
More informationDrag (2) Induced Drag Friction Drag Form Drag Wave Drag
Drag () Induced Drag Friction Drag Form Drag Wave Drag Outline Nomenclature and Concepts Farfield Drag Analysis Induced Drag Multiple Lifting Surfaces Zero Lift Drag :Friction and Form Drag Supersonic
More informationExperimental Study on Flow Control Characteristics of Synthetic Jets over a Blended Wing Body Configuration
Experimental Study on Flow Control Characteristics of Synthetic Jets over a Blended Wing Body Configuration Byunghyun Lee 1), Minhee Kim 1), Chongam Kim 1), Taewhan Cho 2), Seol Lim 3), and Kyoung Jin
More informationAA 242B/ ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B/ ME 242B: Mechanical Vibrations (Spring 2016) Homework #2 Due April 17, 2016 This homework focuses on developing a simplified analytical model of the longitudinal dynamics of an aircraft during
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationEXPERIMENTAL INVESTIGATION OF THE DYNAMIC STABILITY DERIVATIVES FOR A FIGHTER MODEL
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES EXPERIMENTAL INVESTIGATION OF THE DYNAMIC STABILITY DERIVATIVES FOR A FIGHTER MODEL MR Soltani*, Ali R Davari** *Associate Professor, **PhD Student
More informationFlight Vehicle Terminology
Flight Vehicle Terminology 1.0 Axes Systems There are 3 axes systems which can be used in Aeronautics, Aerodynamics & Flight Mechanics: Ground Axes G(x 0, y 0, z 0 ) Body Axes G(x, y, z) Aerodynamic Axes
More informationAirfoils and Wings. Eugene M. Cliff
Airfoils and Wings Eugene M. Cliff 1 Introduction The primary purpose of these notes is to supplement the text material related to aerodynamic forces. We are mainly interested in the forces on wings and
More informationChapter 5 Wing design - selection of wing parameters 2 Lecture 20 Topics
Chapter 5 Wing design - selection of wing parameters Lecture 0 Topics 5..4 Effects of geometric parameters, Reynolds number and roughness on aerodynamic characteristics of airfoils 5..5 Choice of airfoil
More informationCFD ANALYSIS OF AERODYNAMIC HEATING FOR HYFLEX HIGH ENTHALPY FLOW TESTS AND FLIGHT CONDITIONS
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES CFD ANALYSIS OF AERODYNAMIC HEATING FOR HYFLEX HIGH ENTHALPY FLOW TESTS AND FLIGHT CONDITIONS Keiichi Murakami*, Yukimitsu Yamamoto*, Olivier Rouzand**
More informationUsing Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver
Using Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708 0300 A novel nonlinear
More informationLONGITUDINAL STABILITY AND TRIM OF AN ARIANE 5 FLY-BACK BOOSTER
12th AIAA International Space Planes and Hypersonic Systems and Technologies 1-19 December 23, Norfolk, Virginia AIAA 23-7 LONGITUDINAL STABILITY AND TRIM OF AN ARIANE FLY-BACK BOOSTER Th. Eggers DLR,
More informationHigh-Fidelity Aero-Structural Design Optimization of a Supersonic Business Jet
AIAA 2002 1483 High-Fidelity Aero-Structural Design Optimization of a Supersonic Business Jet Joaquim R. R. A. Martins, Juan J. Alonso Stanford University, Stanford, CA 94305 James J. Reuther NASA Ames
More informationNumerical study of battle damaged two-dimensional wings
Advances in Fluid Mechanics IX 141 Numerical study of battle damaged two-dimensional wings S. Djellal, T. Azzam, M. Djellab & K. Lakkaichi Fluid Mechanics Laboratory Polytechnical School Bordj El Bahri,
More informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationCHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH
82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient
More informationFurther Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow
29th AIAA Applied Aerodynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3509 Further Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow Antony Jameson, John C. Vassberg,
More information