Engineering PapersSearch

Engineering topics

Kevin E. Jacobson

Publications and source records attributed to Kevin E. Jacobson.

Flutter-Constrained Optimization with the Linearized Frequency-Domain Approach

Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.

Aeroelasticity

Fixed- and Free-Mode Flutter Derivatives during Aeroelastic Optimization

Aeroelastic flutter design derivatives may be approximated by ignoring the sensitivity of the structural mode shapes. This “fixed-mode” derivative is less-expensive to compute than the exact “free-mode” derivative (which accommodates the mode shape dependency), but also may provide inaccurate sensitivities for optimization. This work formulates both fixed- and free-mode flutter derivatives, and demonstrates the conditions under which they are equivalent. Two different types of fixed-mode derivatives are also developed in this work. Using a cantilevered-plate demonstration problem with both shape and sizing design parameters, the accuracy of the free-mode derivative and the two fixed-mode derivatives, are all demonstrated via gradient-based optimization.

Bret K. Stanford

Linearized Frequency-Domain Gust Analysis and Adjoint-Based Sensitivities

Gust analysis is added to a linearized frequency-domain method in FUN3D, a NASA computational fluid dynamics solver. The method linearizes about a nonlinear static equilibrium condition and is therefore appropriate for problems with small perturbations such as transonic stochastic gust analysis. In addition to the gust analysis, adjoint-based sensitivities of stochastic gust constraints are implemented for multidisciplinary design optimization. The linearized frequency-domain gust model and adjoint-based sensitivities are described and verified. The method is applied to an optimization for mass minimization of the AGARD 445.6 wing subject to a stochastic gust constraint limiting the displacement of the wing tip.

Aeroelasticity

Flutter and Stochastic Gust Constraint Studies of the Undeflected CRM using MPhys

Aircraft design must ultimately address gust and flutter characteristics to achieve airworthiness. As modern aircraft continue to become more lightweight and flexible, these aeroelastic phenomena are more likely to drive the design process and should thus be considered early. Designing for ride quality due to typical gust spectra is also of interest to ensure a smooth passenger experience. The present work demonstrates the new implementation of stochastic gust constraints in MPhys/OpenMDAO, in addition to flutter and static stress constraints, among others. Design studies are shown for the NASA Common Research Model, including both fuel burn minimization as well as multiobjective optimization of fuel burn and landing gross weight. Two sets of optimization problems are covered, differing mainly by the gust quantity of interest, which include tip displacement and root bending moment. Fuel burn penalties imposed by the flutter and gust constraints are also discussed.

aeroelasticity

Gradient-Based Optimization of the Common Research Model Wing Subject to CFD-Based Gust and Flutter Constraints

The linearized frequency-domain method was recently implemented in the stabilized finite element solver in NASA’s FUN3D code. Previous work by the authors used this method for enforcing flutter constraints during gradient-based optimizations. More recently, the solver was expanded to account for continuous (also known as stochastic) gust responses. This paper expands on recent Common Research Model wing optimization work, which demonstrated gradient-based optimization with flutter and stochastic gust constraints, among others. While that work utilized FUN3D for static aeroelastic solutions but relied on doublet lattice aerodynamics for gust and flutter responses, the present work replaces these unsteady aerodynamic analyses with those of FUN3D’s linearized frequency-domain solver. With analytic derivatives available, gradient-based optimization is performed through the use of the OpenMDAO/MPhys libraries with over 700 shape, structural, and aerodynamic design variables and over 10 nonlinear constraints. Comparisons of analysis results and optimized designs are made between doublet lattice and linearized frequency-domain solutions.

aeroelasticity

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition

Closed-Loop Simulations of Human-Scale Mars Lander Descent Trajectories on Frontier

A computational campaign was performed to run high-fidelity, free-flight simulations of a human-scale Mars lander concept vehicle decelerating under retropropulsion through the Martian atmosphere with closed-loop flight control. A novel approach is used to couple computational fluid dynamics (CFD) software with a mature flight mechanics package, where the two applications communicate in real-time across two geographically-dispersed computational facilities. The CFD is performed on the Frontier exascale system located at Oak Ridge National Laboratory, and the flight mechanics are executed on a system located at NASA Langley Research Center. In the current campaign, CFD is performed using finite-rate chemistry to account for the interactions between the LOXCH 4 engines and the CO 2 Martian atmosphere. A simulation of a closed-loop main engine throttling and RCS actuation is presented, demonstrating that the vehicle and model are able to maintain stability in a long-duration CFD-in-the-loop flight simulation. Comparisons are made to a reduced order model ignoring aero-propulsive interactions.

CFD