Development of a Coupled Fluid/Structure Aeroelastic Solver with Applications to Vortex Breakdown-Induced Twin Tail Buffeting

Simulation of tail buffet is studied for several delta wing-vertical tail configurations. Flow conditions are chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on the vertical tail. The dimensions and material properties of the vertical tails are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. This multi-disciplinary problem is solved sequentially for the fluid flow, the elastic deformations and the grid displacements. The flow is simulated by time accurately solving the laminar, unsteady, compressible, Navier-Stokes equations using an implicit, upwind, flux-difference splitting, finite volume scheme. The elastic vibrations of the tail are modeled by coupled bending and torsion beam equations. These equations are solved accurately in time using the Galerkin method and a five-stage, Runge-Kutta-Verner scheme. The grid for the fluid dynamics calculations is continuously deformed using interpolation functions to smoothly disperse the displacements throughout the computational domain. Tail buffet problems are solved for single tail cases, twin F/A-18 tail cases and twin highly swept generic tail cases. The use of an apex flap for buffet control is also computationally studied. The results demonstrate the effects of inertial structural coupling, Reynolds number, aft fuselage geometry and spanwise tail location on the tail buffet loads and response. Favorable comparisons with experimental data indicate that the present aeroelastic method is well suited to providing qualitative insight into the tail buffet problem, as well as quantitative data for refined long duration simulations

Computational Aeroelastic Analysis of Ares Crew Launch Vehicle Bi-Modal Loading

A Reynolds averaged Navier-Stokes analysis, with and without dynamic aeroelastic effects, is presented for the Ares I-X launch vehicle at transonic Mach numbers and flight Reynolds numbers for two grid resolutions and two angles of attack. The purpose of the study is to quantify the force and moment increment imparted by the sudden transition from fully separated flow around the crew module - service module junction to that of the bi-modal flow state in which only part of the flow reattaches. The bi-modal flow phenomenon is of interest to the guidance, navigation and control community because it causes a discontinuous jump in forces and moments. Computations with a rigid structure at zero zero angle of attack indicate significant increases in normal force and pitching moment. Dynamic aeroelastic computations indicate the bi-modal flow state is insensitive to vehicle flexibility due to the resulting deflections imparting only very small changes in local angle of attack. At an angle of attack of 2.5deg, the magnitude of the pitching moment increment resulting from the bi-modal state nearly triples, while occurring at a slightly lower Mach number. Significant grid induced variations between the solutions indicate that further grid refinement is warranted

PAB3D Simulations for the CAWAPI F-16XL

Numerical simulations of the flow around F-16XL are performed as a contribution to the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) using the PAB3D CFD code. Two turbulence models are used in the calculations: a standard k-! model, and the Shih-Zhu-Lumley (SZL) algebraic stress model. Seven flight conditions are simulated for the flow around the F-16XL where the free stream Mach number varies from 0.242 to 0.97. The range of angles of attack varies from 0deg to 20deg. Computational results, surface static pressure, boundary layer velocity profiles, and skin friction are presented and compared with flight data. Numerical results are generally in good agreement with flight data, considering that only one grid resolution is utilized for the different flight conditions simulated in this study. The ASM results are closer to the flight data than the k-! model results. The ASM predicted a stronger primary vortex, however, the origin of the vortex and footprint is approximately the same as in the k-! predictions

Coupled CFD/CSD Analysis of an Active-Twist Rotor in a Wind Tunnel with Experimental Validation

An unsteady Reynolds averaged Navier-Stokes analysis loosely coupled with a comprehensive rotorcraft code is presented for a second-generation active-twist rotor. High fidelity Navier-Stokes results for three configurations: an isolated rotor, a rotor with fuselage, and a rotor with fuselage mounted in a wind tunnel, are compared to lifting-line theory based comprehensive rotorcraft code calculations and wind tunnel data. Results indicate that CFD/CSD predictions of flapwise bending moments are in good agreement with wind tunnel measurements for configurations with a fuselage, and that modeling the wind tunnel environment does not significantly enhance computed results. Actuated rotor results for the rotor with fuselage configuration are also validated for predictions of vibratory blade loads and fixed-system vibratory loads. Varying levels of agreement with wind tunnel measurements are observed for blade vibratory loads, depending on the load component (flap, lag, or torsion) and the harmonic being examined. Predicted trends in fixed-system vibratory loads are in good agreement with wind tunnel measurements

Coupled CFD/CSD Computation of Airloads of an Active-Twist Rotor

An unsteady Reynolds averaged Navier-Stokes analysis loosely coupled with a comprehensive rotorcraft code for blade trim and aeroelastic effects is presented for a second-generation Active-Twist Rotor. Mesh and temporal sensitives of computed airloads are evaluated. In the final paper, computed airloads will be compared with wind tunnel data for the Active-Twist Rotor test that is currently underway

Aerostructural Vortical Flow Interactions with Applications to F/A-18 and F-117 Tail Buffet

The buffet response of the flexible twin-tail configuration of the F/A-18 and a generic F-111 aircraft are computationally simulated and experimentally validated. The problem is a multidisciplinary one which requires the sequential solution of three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational models consist of a 76 deg. swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 or F-117 twin-tail. The configuration is pitched at 30 deg. angle-of-attack. The problem is solved for the initial flow conditions with the twin tails kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the bending and torsional tails responses due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores of the delta wing. Several spanwise locations of the twin tails are investigated. The computational results are validated using several existing experimental data

Temperature-Corrected Model of Turbulence in Hot Jet Flows

An improved correction has been developed to increase the accuracy with which certain formulations of computational fluid dynamics predict mixing in shear layers of hot jet flows. The CFD formulations in question are those derived from the Reynolds-averaged Navier-Stokes equations closed by means of a two-equation model of turbulence, known as the k-epsilon model, wherein effects of turbulence are summarized by means of an eddy viscosity. The need for a correction arises because it is well known among specialists in CFD that two-equation turbulence models, which were developed and calibrated for room-temperature, low Mach-number, plane-mixing-layer flows, underpredict mixing in shear layers of hot jet flows. The present correction represents an attempt to account for increased mixing that takes place in jet flows characterized by high gradients of total temperature. This correction also incorporates a commonly accepted, previously developed correction for the effect of compressibility on mixing

Computation and Validation of Fluid/Structure Twin Tail Buffet Response

The buffet response of the flexible twin-tail/delta wing configuration, a multidisciplinary problem, is solved using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations which are used for obtaining the flow-field vector and the aerodynamic loads on the twin tails. The second set is the coupled aeroelastic equations which are used for obtaining the bending and torsional deflections of the twin tails. The third set is the grid-displacement equations which are used for updating the grid coordinates due to the tail deflections. The computational model is similar to the one used by Washburn et. al. which consists of a delta wing of aspect ratio one and twin tails with taper ratios of 0.23. The vortex of the twin tails are located at the wing trailing edge. The configuration is pitched at 30 deg angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tails fixed as rigid surfaces, the problem is solved for the initial flow conditions. Next, the problem is solved for the twin tail response for uncoupled bending and torsional vibrations due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores. The configuration is investigated for three spanwise positions of the twin tails; inboard, midspan and outboard locations. The computational results are validated and are in very good agreement with the experimental data of Washburn, et. al

Focused and Steady-State Characteristics of Shaped Sonic Boom Signatures: Prediction and Analysis

The objective of this study is to examine the effect of flight, at off-design conditions, on the propagated sonic boom pressure signatures of a small "low-boom" supersonic aircraft. The amplification, or focusing, of the low magnitude "shaped" signatures produced by maneuvers such as the accelerations from transonic to supersonic speeds, climbs, turns, pull-up and pushovers is the concern. To analyze these effects, new and/or improved theoretical tools have been developed, in addition to the use of existing methodology. Several shaped signatures are considered in the application of these tools to the study of selected maneuvers and off-design conditions. The results of these applications are reported in this paper as well as the details of the new analytical tools. Finally, the magnitude of the focused boom problem for "low boom" supersonic aircraft designs has been more accurately quantified and potential "mitigations" suggested. In general, "shaped boom" signatures, designed for cruise flight, such as asymmetric and symmetric flat-top and initial-shock ramp waveforms retain their basic shape during transition flight. Complex and asymmetric and symmetric initial shock ramp waveforms provide lower magnitude focus boom levels than N-waves or asymmetric and symmetric flat-top signatures

Active Flutter Suppression Controllers Derived from Linear and Nonlinear Aerodynamics: Application to a Transport Aircraft Model

Active flutter suppression has been demonstrated in simulation by many researchers, generally using methods based on linear aerodynamics and often with simplistic geometries. In this paper, active flutter suppression is demonstrated in a simulation using a Navier-Stokes aerodynamics code, FUN3D (Fully Unstructured Navier-Stokes Three-Dimensional), and a realistic transport aircraft configuration. This is accomplished using simple observer-feedback controllers derived from linear aeroelastic models, including reduced order models built via FUN3D data. The development of these reduced order models is described here. It is shown that controllers derived from reduced order models of the nonlinear aerodynamics outperform controllers based on linear aerodynamics