Engineering Papers⌕ Search

Engineering topics

Bret K. Stanford

Publications and source records attributed to Bret K. Stanford.

12 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Random Variables with Moment-Matching Staircase Density Functions

This paper proposes a family of random variables for uncertainty modeling. The variables of interest have a bounded support set, and prescribed values for the first four moments. We present the feasibility conditions for the existence of any of such variables, and propose a class of variables that conforms to such constraints. This class is called staircase because the density of its members is a piecewise constant function. Convex optimization is used to calculate their distributions according to several optimality criteria, including maximal entropy and maximal log-likelihood. The flexibility and efficiency of staircases enable modeling phenomena having a possibly skewed and/or multimodal response at a low computational cost. Furthermore, we provide a means to account for the uncertainty in the distribution caused by estimating staircases from data. These ideas are illustrated by generating empirical staircase predictor models. We consider the case in which the predictor matches the sample moments exactly (a setting applicable to large datasets), as well as the case in which the predictor accounts for the sampling error in such moments (a setting applicable to sparse datasets). A predictor model for the dynamics of an aeroelastic airfoil subject to flutter instability is used as an example. The resulting predictor not only describes the system's response accurately, but also enables carrying out a risk analysis for safe flight.

Luis G. Crespo↗
NASA NTRSDate not supplied

Topology Optimization of Low-Speed Aeroelastic Wind Tunnel Models

The goal of this work is to utilize topology optimization methods to obtain low-speed flexible wind tunnel models with simple jig outer mold lines, but highly-tuned aeroelastic flutter behavior. The wing is conceptually constructed with commonly-available 3D printed materials, of which the optimizer may select one of two at every location in the wing; a stiffer material meant to bear the loads, and a softer rubber-like material for space filling, and to form the outer mold line, where needed. Results are shown for three topological parameterizations: 2D topologies with a fixed nondesignable plate through the camber line, 2D parameterizations without the plate, and 3D parameterizations where the topology can vary though the depth of the wing, at a given planform location.

Aeroelasticity↗
NASA NTRSDate not supplied

Aeroservoelastic Control Law Development for the Integrated Adaptive Wing Technology Maturation Wind-Tunnel Test

In this paper, a linear quadratic Gaussian (LQG) regulator was applied to a computational aeroservoelastic (ASE) model of the Integrated Adaptive Wing Technology Maturation (IAWTM) wind-tunnel model to demonstrate its effectiveness for active flutter suppression. The aerodynamics within this computational ASE model were solved using unsteady Reynolds-averaged Navier-Stokes (RANS) equations; the structural dynamics of the model were solved using modal analysis. This RANS-based ASE model was used as the plant for the LQG regulator. The Kalman filter and controller used within the LQG regulator were derived from a reduced-order model (ROM) of the RANS-based ASE model. This ROM-based ASE model was generated using system identification techniques. It was shown that the Kalman filter adequately predicted important states of the plant, and that LQG regulator successfully stabilized an aeroelastically unstable plant without exceeding reasonable control surface limitations. All of this was demonstrated using simulation data, and the methods presented here appear feasible for use in the upcoming IAWTM wind-tunnel tests.

Josiah M. Waite↗
NASA NTRSDate not supplied

Aeroelastic Modeling and CFD Simulation of Wind-Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents an aeroelastic modeling and simulation study of an aspect ratio 13.5 wind-tunnel scale Common Research Model (CRM) with distributed flaps. A vortex-lattice VSPAERO model of the CRM model is developed. A transonic small disturbance/integral boundary layer correction method is implemented in the VSPAERO model to account for the transonic and viscous flow effects. The structural deformation of the CRM model is calculated using a NASTRAN equivalent beam model. The VSPAERO model is coupled to the NASTRAN equivalent beam model to provide a rapid aero-structural analysis. A validation of the VSPAERO aeroelastic model is conducted by comparing the results to FUN3D CFD aeroelastic simulation results. An aerodynamic database is generated using the developed VSPAERO aeroelastic model for the real-time drag optimization and maneuver load alleviation study of the wind-tunnel scale CRM model.

Aeroelastic Modeling↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗