Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Aeroelasticity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Parallel computation of aerodynamic influence coefficients for aeroelastic analysis on a transputer network

Aeroelastic analysis is mult-disciplinary and computationally expensive. Hence, it can greatly benefit from parallel processing. As part of an effort to develop an aeroelastic analysis capability on a distributed-memory transputer network, a parallel algorithm for the computation of aerodynamic influence coefficients is implemented on a network of 32 transputers. The aerodynamic influence coefficients are calculated using a three-dimensional unsteady aerodynamic model and a panel discretization. Efficiencies up to 85 percent are demonstrated using 32 processors. The effects of subtask ordering, problem size and network topology are presented. A comparison to results on a shared-memory computer indicates that higher speedup is achieved on the distributed-memory system.

Janetzke, D. C.↗

Identification of Computational and Experimental Reduced-Order Models

The identification of computational and experimental reduced-order models (ROMs) for the analysis of unsteady aerodynamic responses and for efficient aeroelastic analyses is presented. For the identification of a computational aeroelastic ROM, the CFL3Dv6.0 computational fluid dynamics (CFD) code is used. Flutter results for the AGARD 445.6 Wing and for a Rigid Semispan Model (RSM) computed using CFL3Dv6.0 are presented, including discussion of associated computational costs. Modal impulse responses of the unsteady aerodynamic system are computed using the CFL3Dv6.0 code and transformed into state-space form. The unsteady aerodynamic state-space ROM is then combined with a state-space model of the structure to create an aeroelastic simulation using the MATLAB/SIMULINK environment. The MATLAB/SIMULINK ROM is then used to rapidly compute aeroelastic transients, including flutter. The ROM shows excellent agreement with the aeroelastic analyses computed using the CFL3Dv6.0 code directly. For the identification of experimental unsteady pressure ROMs, results are presented for two configurations: the RSM and a Benchmark Supercritical Wing (BSCW). Both models were used to acquire unsteady pressure data due to pitching oscillations on the Oscillating Turntable (OTT) system at the Transonic Dynamics Tunnel (TDT). A deconvolution scheme involving a step input in pitch and the resultant step response in pressure, for several pressure transducers, is used to identify the unsteady pressure impulse responses. The identified impulse responses are then used to predict the pressure responses due to pitching oscillations at several frequencies. Comparisons with the experimental data are then presented.

Silva, Walter A.↗

Effects of symmetric and asymmetric modes on transonic aeroelastic characteristics of full-span wing-body configurations

To accurately study the transonic aeroelastic characteristics, it is important to model the full aircraft configuration, including asymmetry. Recently, an accurate method of computing unsteady transonic flows on full-span wing-body configurations was developed using the transonic potential flow theory. In this work, the method is further developed to account for the aeroelasticity of full-span wing-body configurations. This is accomplished by simultaneously integrating the unsteady aerodynamic forces and modal structural equations of the wing-body configurations. To validate the method, aeroelastic computations are made for a wing-body configuration with a rectangular wing. The aeroelastic responses of this configuration are correlated with the responses of a similar isolated wing. The comparisons are favorable. Aeroelastic computations associated with symmetric and asymmetric modes are also made to study the influence of modal asymmetry on responses. This new development is further illustrated by computing aeroelastic characteristics of a typical fighter aircraft. The results from this study will be useful in accurately computing the transonic flutter boundaries of aircraft, including those associated with asymmetric modes.

Guruswamy, Guru P.↗

CFD Based Computations of Flexible Helicopter Blades for Stability Analysis

As a collaborative effort among government aerospace research laboratories an advanced version of a widely used computational fluid dynamics code, OVERFLOW, was recently released. This latest version includes additions to model flexible rotating multiple blades. In this paper, the OVERFLOW code is applied to improve the accuracy of airload computations from the linear lifting line theory that uses displacements from beam model. Data transfers required at every revolution are managed through a Unix based script that runs jobs on large super-cluster computers. Results are demonstrated for the 4-bladed UH-60A helicopter. Deviations of computed data from flight data are evaluated. Fourier analysis post-processing that is suitable for aeroelastic stability computations are performed.

Guruswamy, Guru P.↗

Overview of the Aeroelastic Prediction Workshop

The AIAA Aeroelastic Prediction Workshop (AePW) was held in April, 2012, bringing together communities of aeroelasticians and computational fluid dynamicists. The objective in conducting this workshop on aeroelastic prediction was to assess state-of-the-art computational aeroelasticity methods as practical tools for the prediction of static and dynamic aeroelastic phenomena. No comprehensive aeroelastic benchmarking validation standard currently exists, greatly hindering validation and state-of-the-art assessment objectives. The workshop was a step towards assessing the state of the art in computational aeroelasticity. This was an opportunity to discuss and evaluate the effectiveness of existing computer codes and modeling techniques for unsteady flow, and to identify computational and experimental areas needing additional research and development. Three configurations served as the basis for the workshop, providing different levels of geometric and flow field complexity. All cases considered involved supercritical airfoils at transonic conditions. The flow fields contained oscillating shocks and in some cases, regions of separation. The computational tools principally employed Reynolds-Averaged Navier Stokes solutions. The successes and failures of the computations and the experiments are examined in this paper.

Heeg, Jennifer↗

Coupled Aerodynamic and Structural Sensitivity Analysis of a High-Speed Civil Transport

An objective of the High Performance Computing and Communication Program at the NASA Langley Research Center is to demonstrate multidisciplinary shape and sizing optimization of a complete aerospace vehicle configuration by using high-fidelity, finite-element structural analysis and computational fluid dynamics aerodynamic analysis. In a previous study, a multi-disciplinary analysis system for a high-speed civil transport was formulated to integrate a set of existing discipline analysis codes, some of them computationally intensive, This paper is an extension of the previous study, in which the sensitivity analysis for the coupled aerodynamic and structural analysis problem is formulated and implemented. Uncoupled stress sensitivities computed with a constant load vector in a commercial finite element analysis code are compared to coupled aeroelastic sensitivities computed by finite differences. The computational expense of these sensitivity calculation methods is discussed.

Mason, B. H.↗

Comparison between computations and experimental data in unsteady three-dimensional transonic aerodynamics, including aeroelastic applications

Comparisons were made of computed and experimental data in three-dimensional unsteady transonic aerodynamics, including aeroelastic applications. The computer code LTRAN3, which is based on small-disturbance aerodynamic theory, was used to obtain the aerodynamic data. A procedure based on the U-g method was developed to compute flutter boundaries by using the unsteady aerodynamic coefficients obtained from LTRAN3. The experimental data were obtained from available NASA publications. All the studies were conducted for thin, unswept, rectangular wings with circular-arc cross sections. Numerical and experimental steady and unsteady aerodynamic data were compared for a wing with an aspect ratio of 3 and a thickness ratio of 5% at Mach numbers of 0.7 and 0.9. Flutter data were compared for a wing with an aspect ratio of 5. Two thickness ratios, 6% at Mach numbers of 0.715, 0.851, and 0.913, and 4% at Mach number of 0.904, were considered. Based on the unsteady aerodynamic data obtained from LTRAN3, flutter boundaries were computed; they were compared with those obtained from experiments and the code NASTRAN, which uses linear aerodynamics.

Guruswamy, P.↗

Evaluation of Linear, Inviscid, Viscous, and Reduced-Order Modeling Aeroelastic Solutions of the AGARD 445.6 Wing Using Root Locus Analysis

Reduced-order modeling (ROM) methods are applied to the CFD-based aeroelastic analysis of the AGARD 445.6 wing in order to gain insight regarding well-known discrepancies between the aeroelastic analyses and the experimental results. The results presented include aeroelastic solutions using the inviscid CAP-TSD code and the FUN3D code (Euler and Navier-Stokes). Full CFD aeroelastic solutions and ROM aeroelastic solutions, computed at several Mach numbers, are presented in the form of root locus plots in order to better reveal the aeroelastic root migrations with increasing dynamic pressure. Important conclusions are drawn from these results including the ability of the linear CAP-TSD code to accurately predict the entire experimental flutter boundary (repeat of analyses performed in the 1980's), that the Euler solutions at supersonic conditions indicate that the third mode is always unstable, and that the FUN3D Navier-Stokes solutions stabilize the unstable third mode seen in the Euler solutions.

Silva, Walter A.↗

Application of Computational Stability and Control Techniques Including Unsteady Aerodynamics and Aeroelastic Effects

The motivation behind the inclusion of unsteady aerodynamics and aeroelastic effects in the computation of stability and control (S&C) derivatives will be discussed as they pertain to aeroelastic and aeroservoelastic analysis. This topic will be addressed in the context of two applications, the first being the estimation of S&C derivatives for a cable-mounted aeroservoelastic wind tunnel model tested in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The second application will be the prediction of the nonlinear aeroservoelastic phenomenon known as Residual Pitch Oscillation (RPO) on the B-2 Bomber. Techniques and strategies used in these applications to compute S&C derivatives and perform flight simulations will be reviewed, and computational results will be presented.

Schuster, David M.↗

An integrated computer procedure for sizing composite airframe structures

A computerized algorithm to generate cross-sectional dimensions and fiber orientations for composite airframe structures is described, and its application in a wing structural synthesis is established. The algorithm unifies computations of aeroelastic loads, stresses, and deflections, as well as optimal structural sizing and fiber orientations in an open-ended system of integrated computer programs. A finite-element analysis and a mathematical-optimization technique are discussed.

Sobieszczanski-Sobieski, J.↗

Flexible Launch Vehicle Stability Analysis Using Steady and Unsteady Computational Fluid Dynamics

Launch vehicles frequently experience a reduced stability margin through the transonic Mach number range. This reduced stability margin can be caused by the aerodynamic undamping one of the lower-frequency flexible or rigid body modes. Analysis of the behavior of a flexible vehicle is routinely performed with quasi-steady aerodynamic line loads derived from steady rigid aerodynamics. However, a quasi-steady aeroelastic stability analysis can be unconservative at the critical Mach numbers, where experiment or unsteady computational aeroelastic analysis show a reduced or even negative aerodynamic damping.Amethod of enhancing the quasi-steady aeroelastic stability analysis of a launch vehicle with unsteady aerodynamics is developed that uses unsteady computational fluid dynamics to compute the response of selected lower-frequency modes. The response is contained in a time history of the vehicle line loads. A proper orthogonal decomposition of the unsteady aerodynamic line-load response is used to reduce the scale of data volume and system identification is used to derive the aerodynamic stiffness, damping, and mass matrices. The results are compared with the damping and frequency computed from unsteady computational aeroelasticity and from a quasi-steady analysis. The results show that incorporating unsteady aerodynamics in this way brings the enhanced quasi-steady aeroelastic stability analysis into close agreement with the unsteady computational aeroelastic results.

Bartels, Robert E.↗

NACA0012 benchmark model experimental flutter results with unsteady pressure distributions

The Structural Dynamics Division at NASA Langley Research Center has started a wind tunnel activity referred to as the Benchmark Models Program. The primary objective of this program is to acquire measured dynamic instability and corresponding pressure data that will be useful for developing and evaluating aeroelastic type computational fluid dynamics codes currently in use or under development. The program is a multi-year activity that will involve testing of several different models to investigate various aeroelastic phenomena. This paper describes results obtained from a second wind tunnel test of the first model in the Benchmark Models Program. This first model consisted of a rigid semispan wing having a rectangular planform and a NACA 0012 airfoil shape which was mounted on a flexible two degree of freedom mount system. Experimental flutter boundaries and corresponding unsteady pressure distribution data acquired over two model chords located at the 60 and 95 percent span stations are presented.

Rivera, Jose A., Jr.↗

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three-dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first-step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings. The paper describes the modifications that were made to the code and presents results and comparisons which assess the capability.

Robinson, Brian A.↗

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings.

Robinson, Brian A.↗

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗

Functionality Improvements to Overaero

The functionality of the overset, static aeroelasticity, Navier-Stokes flow solver OVERAERO was increased by adding capability to the flow solver and enhancing code performance. Improvements were made to the fluids/structure interface, an MLP version of the parallel OVERAERO code was developed, and the OVERAERO-MPI code was ported to the Cray T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested successfully on the IPG testbed and a means of reducing communication overhead within OVERFLOW-MPI was investigated. To solve an aeroelastic problem computationally, a structures grid surface definition and a fluids grid surface definition are required. Typically, the structures grid surface has a lower fidelity than the fluids grid surface. Thus, the methods developed to transfer data between the two grid systems are vital to the accuracy and efficiency of the aeroelasticity code. The fluids/structures interface developed for the OVERAERO code was improved to more accurately treat fluids surfaces that bridge between two different structural surfaces. For example, the method allowed the forward portion of a flap track fairing to deform with the wing and the aft end of the fairing to deform with the flap. A tightly-coupled version of the code based on OVERFLOW-MLP was developed to improve code performance on the SGI Origin 2000. This required a new parallelization strategy to couple the fluids and structures codes. The OVERAERO-MPI code was ported to the Cray T3E to extend the usability of the code. The port required extensive use of dynamic memory management techniques to fit large problems within the memory limitations of the T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested on the IPG testbed being developed within NASA. For small problems with minimal data transfer between grids, there was little to no performance penalty spreading the computation across two machines. For very large problems, methods were developed to minimize intermachine communication via the grid partitioning scheme. By minimizing the intermachine communication requirements of the problem, it may still be beneficial to run a tightly-coupled flow solver across two machines within the IPG.

Gee, Ken↗

FUN3D Airload Predictions for the Full-Scale UH-60A Airloads Rotor in a Wind Tunnel

An unsteady Reynolds-Averaged Navier-Stokes solver for unstructured grids, FUN3D, is used to compute the rotor performance and airloads of the UH-60A Airloads Rotor in the National Full-Scale Aerodynamic Complex (NFAC) 40- by 80-foot Wind Tunnel. The flow solver is loosely coupled to a rotorcraft comprehensive code, CAMRAD-II, to account for trim and aeroelastic deflections. Computations are made for the 1-g level flight speed-sweep test conditions with the airloads rotor installed on the NFAC Large Rotor Test Apparatus (LRTA) and in the 40- by 80-ft wind tunnel to determine the influence of the test stand and wind-tunnel walls on the rotor performance and airloads. Detailed comparisons are made between the results of the CFD/CSD simulations and the wind tunnel measurements. The computed trends in solidity-weighted propulsive force and power coefficient match the experimental trends over the range of advance ratios and are comparable to previously published results. Rotor performance and sectional airloads show little sensitivity to the modeling of the wind-tunnel walls, which indicates that the rotor shaft-angle correction adequately compensates for the wall influence up to an advance ratio of 0.37. Sensitivity of the rotor performance and sectional airloads to the modeling of the rotor with the LRTA body/hub increases with advance ratio. The inclusion of the LRTA in the simulation slightly improves the comparison of rotor propulsive force between the computation and wind tunnel data but does not resolve the difference in the rotor power predictions at mu = 0.37. Despite a more precise knowledge of the rotor trim loads and flight condition, the level of comparison between the computed and measured sectional airloads/pressures at an advance ratio of 0.37 is comparable to the results previously published for the high-speed flight test condition.

Lee-Rausch, Elizabeth M.↗

Aeroelasticity of wing and wing-body configurations on parallel computers

The objective of this research is to develop computationally efficient methods for solving aeroelasticity problems on parallel computers. Both uncoupled and coupled methods are studied in this research. For the uncoupled approach, the conventional U-g method is used to determine the flutter boundary. The generalized aerodynamic forces required are obtained by the pulse transfer-function analysis method. For the coupled approach, the fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗