Yu, L.; Song, W. P.
A hybrid prediction method for aerodynamic noise is performed using high order accuracy method in this paper. The method combines a two-dimensional Unsteady Reynolds-Averaged Navier-Stokes(URANS) solver with the acoustic analogy method using Ffowcs Williams-Hawkings equation with penetrable data surface (FW-Hpds). Tandem cylinders are chosen to validate the prediction method. The computations are conducted at a Reynolds number of 1.66 × 105 based on the cylinder diameter. Both the aerodynamic and acoustic results show good agreement with the experimental data, showing a successful application of the hybrid prediction method using two-dimensional URANS simulation.
Strickland, J.H.; Meyer, J.
The purpose of this paper is to provide a brief review of vortex methods with application to parachute aerodynamics. A somewhat generalized discussion of analysis techniques which are applicable to development of both two- and three-dimensional numerical solutions will be presented. A brief review of results from several bluff body simulations will be presented along with very recent results from work being conducted by Sandia National Laboratories in this area. 32 refs.
Benepe, D. B., Sr.
A set of empirical methods was developed to predict low-speed lift, drag and pitching-moment variations with angle of attack for a class of low aspect ratio irregular planform wings suitable for application to advanced aerospace vehicles. The data base, an extensive series of wind-tunnel tests accomplished by the Langley Research Center of the National Aeronautics and Space Administration, is summarized. The approaches used to analyze the wind tunnel data, the evaluation of previously existing methods, data correlation efforts, and the development of the selected methods are presented and discussed. A summary of the methods is also presented to document the equations, computational charts and design guides which have been programmed for digital computer solution. Comparisons of predictions and test data are presented which show that the new methods provide a significant improvement in capability for evaluating the landing characteristics of advanced aerospace vehicles during the preliminary design phase of the configuration development cycle.
Giesing, J. P.; Kalman, T. P.; Rodden, W. P.
A method for correcting discrete element lifting surface theory to reflect given experimental data is presented. Theoretical pressures are modified such that imposed constraints are satisfied while minimizing the changes to the pressures. Several types of correction procedures are presented and correlated; (1) scaling of pressures; (2) scaling of downwash values; and (3) addition of an increment to the downwash that is proportioned to pressure. Some special features are included in these methods and they include: (1) consideration of experimental data from multiple deflection modes, (2) limitation of the amplitudes of the correction factors, and (3) the use of correction factor mode shapes. These methods are correlated for cases involving all three Mach Number ranges using a FORTRAN IV computer program. Subsonically, a wing with an oscillating partial span control surface and a wing with a leading edge droop are presented. Transonically a two-dimensional airfoil with an oscillating flap is considered. Supersonically an arrow wing with and without camber is analyzed. In addition to correction factor methods an investigation is presented dealing with a new simplified transonic modification of the two-dimensional subsonic lifting surface theory. Correlations are presented for an airfoil with an oscillating flap.
Lee-Rausch, Elizabeth M.; Batina, John T.
Modifications to an existing 3D implicit upwind Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the governing flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 deg swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.
Lee-Rausch, Elizabeth M.; Batina, John T.
Modifications to an existing three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1) for the aeroelastic analysis of wings are described. These modifications, which were previously added to CFL3D Version 1.0, include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the government flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 degree swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.
Landrum, E. J.; Miller, D. S.
Trends toward the automation of the design process for airplanes and missiles accentuate the need for analytic techniques for the prediction of aerodynamic characteristics. A number of computer codes have been developed or are under development which show promise of significantly improving the estimation of aerodynamic characteristics for arbitrarily-shaped bodies at supersonic speeds. The programs considered range in complexity from a simple linearized solution employing slender body theory to an exact finite difference solution of the Euler equations. The results from five computer codes are compared with experimental data to determine the accuracy, range of applicability, ease of use, and computer time and cost of the programs. The results provide a useful guide for selecting the appropriate method for treating bodies at the various levels of an automated design process.
Bijleveld, H. A.; Veldman, A. E. P.
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
Durston, Donald A.
LinAir is a vortex lattice aerodynamic prediction method similar to Weissinger's extended lifting-line theory, except that the circulation around a wing is represented by discrete horseshoe vortices, not a continuous distribution of vorticity. The program calculates subsonic longitudinal and lateral/directional aerodynamic forces and moments for arbitrary aircraft geometries. It was originally written by Dr. Ilan Kroo of Stanford University, and subsequently modified by the author to simplify modeling of complex configurations. The Polhamus leading-edge suction analogy was added by the author to extend the range of applicability of LinAir to low aspect ratio (i.e., fighter-type) configurations. A brief discussion of the theory of LinAir is presented, and details on how to run the program are given along with some comparisons with experimental data to validate the code. Example input and output files are given in the appendices to aid in understanding the program and its use. This version of LinAir runs in the VAX/VMS, Cray UNICOS, and Silicon Graphics Iris workstation environments at the time of this writing.
Albright, Alan E.; Dixon, Charles J.; Hegedus, Martin C.
A conceptual/preliminary design level subsonic aerodynamic prediction code HASC (High Angle of Attack Stability and Control) has been improved in several areas, validated, and documented. The improved code includes improved methodologies for increased accuracy and robustness, and simplified input/output files. An engineering method called VTXCHN (Vortex Chine) for prediciting nose vortex shedding from circular and non-circular forebodies with sharp chine edges has been improved and integrated into the HASC code. This report contains a summary of modifications, description of the code, user's guide, and validation of HASC. Appendices include discussion of a new HASC utility code, listings of sample input and output files, and a discussion of the application of HASC to buffet analysis.
Jones, K. M.
A nonlinear aerodynamic prediction technique which solves the conservative full potential equation has been applied to the analysis of three waverider configurations. This technique was selected based on its capability to analyze the off-design characteristics of the waveriders. Very good correlations were achieved with surface pressure data for both the Mach 4 elliptic cone waverider and the Mach 6 caret-wing derivative. Off-design Mach number and angle-of-attack pressure correlations were very good for the elliptic cone waverider. The range of correlation with data exceeded that expected based on the theory limitations. A surface pressure integration routine was demonstrated and agreement between predicted aerodynamic forces and experimental force data for the Mach 4 waverider was excellent. Analysis of a nonconical waverider configuration was initiated where a discrete input option is used to achieve the computational gridding. Preliminary analysis of this configuration indicates the correct shock location will be predicted.
Gilmer, B. R.; Miner, G. A.; Bristow, D. R.
A numerical procedure was developed for the aerodynamic force and moment analysis of V/STOL aircraft operating in the transition regime between hover and conventional forward flight. The trajectories, cross sectional area variations, and mass entrainment rates of the jets are calculated by the Adler-Baron Jet-in-Crossflow Program. The inviscid effects of the interaction between the jets and airframe on the aerodynamic properties are determined by use of the MCAIR 3-D Subsonic properties are determined by use of the MCAIR 3-D Subsonic Potential Flow Program, a surface panel method. In addition, the MCAIR 3-D Geometry influence Coefficient Program is used to calculate a matrix of partial derivatives that represent the rate of change of the inviscid aerodynamic properties with respect to arbitrary changes in the effective wing shape.
Pamadi, Bandu N.; Taylor, Lawrence W., Jr.
A semi-empirical method is presented for the estimation of aerodynamic forces and moments acting on a steadily spinning (rotating) light airplane. The airplane is divided into wing, body, and tail surfaces. The effect of power is ignored. The strip theory is employed for each component of the spinning airplane to determine its contribution to the total aerodynamic coefficients. Then, increments to some of the coefficients which account for centrifugal effect are estimated. The results are compared to spin tunnel rotary balance test data.
Goldman, L. J.; Augler, R. E.
A generalized analysis to predict the two-dimensional aerodynamic losses of film-cooled vanes by using integral boundary-layer parameters is presented. Heat-transfer and trailing-edge injection effects are included in the method. An approximate solution of the generalized equations is also included to show more clearly the effect of the different boundary-layer and cooling parameters on the losses. The analytical predictions agree well with the experimental results, indicating that available boundary-layer calculations for cooled vanes are of sufficient accuracy to use in the prediction method.
Schuster, David M.
Over the past three years, the National Aeronautics and Space Administration (NASA) has initiated design, development, and testing of a new human-rated space exploration system under the Constellation Program. Initial designs within the Constellation Program are scheduled to replace the present Space Shuttle, which is slated for retirement within the next three years. The development of vehicles for the Constellation system has encountered several unsteady aerodynamics challenges that have bearing on more traditional unsteady aerodynamic and aeroelastic analysis. This paper focuses on the synergy between the present NASA challenges and the ongoing challenges that have historically been the subject of research and method development. There are specific similarities in the flows required to be analyzed for the space exploration problems and those required for some of the more nonlinear unsteady aerodynamic and aeroelastic problems encountered on aircraft. The aggressive schedule, significant technical challenge, and high-priority status of the exploration system development is forcing engineers to implement existing tools and techniques in a design and application environment that is significantly stretching the capability of their methods. While these methods afford the users with the ability to rapidly turn around designs and analyses, their aggressive implementation comes at a price. The relative immaturity of the techniques for specific flow problems and the inexperience with their broad application to them, particularly on manned spacecraft flight system, has resulted in the implementation of an extensive wind tunnel and flight test program to reduce uncertainty and improve the experience base in the application of these methods. This provides a unique opportunity for unsteady aerodynamics and aeroelastic method developers to test and evaluate new analysis techniques on problems with high potential for acquisition of test and even flight data against which they
A general low-order surface-singularity panel method is used to predict the aerodynamic characteristics of a problem where a wing-tip vortex from one wing closely interacts with an aft mounted wing in a low Reynolds Number flow; i.e., 125,000. Nonlinear effects due to wake roll-up and the influence of the wings on the vortex path are included in the calculation by using a coupled iterative wake relaxation scheme. The interaction also affects the wing pressures and boundary layer characteristics: these effects are also considered using coupled integral boundary layer codes and preliminary calculations using free vortex sheet separation modelling are included. Calculated results are compared with water tunnel experimental data with generally remarkably good agreement.
Cruz, Josue; Anders, Scott G.
The application of a Computational Fluid Dynamics tool to a jet flap control effector on an elliptical airfoil-section wing was investigated. The study utilized the Tetrahedral Unstructured Software System developed at NASA Langley Research Center. The Reynolds-averaged Navier-Stokes flow solver code used was USM3D. The CFD-based jet flap simulations were compared to experimental results from a wind tunnel test conducted at the NASA Langley Transonic Dynamics Tunnel. The wind tunnel model consisted of a six percent thick elliptical airfoil with a modified trailing edge. The jet flap was located at 95% chord and exited at 90 degrees to the lower surface. The experimental model was designed to promote two-dimensional flow across the wing. It was found that the CFD simulation had to model the three-dimensional geometry of the experiment in order to obtain good agreement. Tests were performed at two Mach numbers at several different jet momentum coefficients. In order to be consistent with the experimental method, the CFD lift and pitching moment values were determined by integrating the pressures over the wing.
Zorumski, W. E. (Editor); Weir, D. S. (Editor)
The prediction sequence used in the aircraft noise prediction program (ANOPP) is described. The elements of the sequence are called program modules. The first group of modules analyzes the propeller geometry, the aerodynamics, including both potential and boundary-layer flow, the propeller performance, and the surface loading distribution. This group of modules is based entirely on aerodynamic strip theory. The next group of modules deals with the first group. Predictions of periodic thickness and loading noise are determined with time-domain methods. Broadband noise is predicted by a semiempirical method. Near-field predictions of fuselage surface pressrues include the effects of boundary layer refraction and scattering. Far-field predictions include atmospheric and ground effects.
Lamar, J. E.
The development of a nonplanar lifting surface method having a continuous distribution of singularities and satisfying the tangent flow boundary condition on the mean camber surface is given. The method predicts some incompressible longitudinal aerodynamic coefficients of rectangular wings which have circular-arc camber. The solution method is of the integral-equation type and the resulting surface integrals are evaluated by either using numerical or analytical techniques, as are appropriate. Applications are made and the results compared with those from an exact two-dimensional circular-arc camber solution, a three-dimensional flat-wing solution which represents the camber by a projected slope onto the flat surface, and a flat-wing experiment. From these comparisons, the present method is found to predict well the flat-wing experiment and limiting values, in addition to the center of pressure variation at an angle of attack of zero for any camber. For wings having camber ratios larger than about 1.25% and moderate to high aspect ratios, the results deterioriate due to the inadequacy of lifting pressure modes employed.
Bidwell, Colin; Potapczuk, Mark; Addy, Gene; Wright, William
LEWICE is a computer program that predicts the accumulation of ice on two-dimensional aerodynamic surfaces under conditions representative of the flight of an aircraft through an icing cloud. The software first calculates the airflow surrounding the body of interest, then uses the airflow to compute the trajectories of water droplets that impinge on the surface of the body. The droplet trajectories are also used to compute impingement limits and local collection efficiencies, which are used in subsequent ice-growth calculations and are also useful for designing systems to protect against icing. Next, the software predicts the shape of accumulating ice by modeling transfers of mass and energy in small control volumes. The foregoing computations are repeated over several computational time steps until the total icing exposure time is reached. Results of computations by LEWICE have been compared with an extensive database of measured ice shapes obtained from experiments, and have been shown to closely approximate those shapes under most conditions of interest to the aviation community.
Griffin, Lisa W.; Mcconnaughey, H. V.
A numerical study of the aerodynamic and thermal environment associated with axial turbine stages is presented. Computations were performed using a modification of the unsteady NASA Ames viscous code, ROTOR1, and an improved version of the NASA Lewis steady inviscid cascade system MERIDL-TSONIC coupled with boundary layer codes BLAYER and STAN5. Two different turbine stages were analyzed: the first stage of the United Technologies Research Center Large Scale Rotating Rig (LSRR) and the first stage of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump turbine. The time-averaged airfoil midspan pressure and heat transfer profiles were predicted for numerous thermal boundary conditions including adiabatic wall, prescribed surface temperature, and prescribed heat flux. Computed solutions are compared with each other and with experimental data in the case of the LSRR calculations. Modified ROTOR1 predictions of unsteady pressure envelopes and instantaneous contour plots are also presented for the SSME geometry. Relative merits of the two computational approaches are discussed.
Hrach, Frank J. (Inventor)
In the present method, boundary layer thickening is combined with laminar flow control to reduce drag. An aerodynamic body is accelerated enabling a ram turbine on the body to receive air at velocity V sub 0. The discharge air is directed over an aft portion of the aerodynamic body producing boundary layer thickening. The ram turbine also drives a compressor by applying torque to a shaft connected between the ram turbine and the compressor. The compressor sucks in lower boundary layer air through inlets in the shell of the aircraft producing laminar flow control and reducing drag. The discharge from the compressor is expanded in a nozzle to produce thrust.
Mcmillan, O. J.; Mendenhall, M. R.; Perkins, S. C., Jr.
Work is described dealing with two areas which are dominated by the nonlinear effects of vortex flows. The first area concerns the stall/spin characteristics of a general aviation wing with a modified leading edge. The second area concerns the high-angle-of-attack characteristics of high performance military aircraft. For each area, the governing phenomena are described as identified with the aid of existing experimental data. Existing analytical methods are reviewed, and the most promising method for each area used to perform some preliminary calculations. Based on these results, the strengths and weaknesses of the methods are defined, and research programs recommended to improve the methods as a result of better understanding of the flow mechanisms involved.
Kroll, R. I.; Miller, R. D.
The results of a study are presented in which the applicability at subsonic speeds of several aerodynamic methods for predicting dynamic gust loads on aircraft, including active control systems, was examined and compared. These aerodynamic methods varied from steady state to an advanced unsteady aerodynamic formulation. Brief descriptions of the structural and aerodynamic representations and of the motion and load equations are presented. Comparisons of numerical results achieved using the various aerodynamic methods are shown in detail. From these results, aerodynamic representations for dynamic gust analyses are identified. It was concluded that several aerodynamic methods are satisfactory for dynamic gust analyses of configurations having either controls fixed or active control systems that primarily affect the low frequency rigid body aircraft response.
Lamar, J. E.; Gloss, B. B.
Because the potential flow suction along the leading and side edges of a planform can be used to determine both leading- and side-edge vortex lift, the present investigation was undertaken to apply the vortex-lattice method to computing side-edge suction force for isolated or interacting planforms. Although there is a small effect of bound vortex sweep on the computation of the side-edge suction force, the results obtained for a number of different isolated planforms produced acceptable agreement with results obtained from a method employing continuous induced-velocity distributions. By using the method outlined, better agreement between theory and experiment was noted for a wing in the presence of a canard than was previously obtained.
Debraux, Pierre; Grappe, Frederic; Manolova, Aneliya V; Bertucci, William
When cycling on level ground at a speed greater than 14 m/s, aerodynamic drag is the most important resistive force. About 90% of the total mechanical power output is necessary to overcome it. Aerodynamic drag is mainly affected by the effective frontal area which is the product of the projected frontal area and the coefficient of drag. The effective frontal area represents the position of the cyclist on the bicycle and the aerodynamics of the cyclist-bicycle system in this position. In order to optimise performance, estimation of these parameters is necessary. The aim of this study is to describe and comment on the methods used during the last 30 years for the evaluation of the effective frontal area and the projected frontal area in cycling, in both laboratory and actual conditions. Most of the field methods are not expensive and can be realised with few materials, providing valid results in comparison with the reference method in aerodynamics, the wind tunnel. Finally, knowledge of these parameters can be useful in practice or to create theoretical models of cycling performance. PMID:21936289
Heeg, Jennifer; Chawlowski, Pawel
The AIAA Aeroelastic Prediction Workshop (AePW) was held in April 2012, bringing together communities of aeroelasticians, computational fluid dynamicists and experimentalists. The extended objective was to assess the state of the art in computational aeroelastic methods as practical tools for the prediction of static and dynamic aeroelastic phenomena. As a step in this process, workshop participants analyzed unsteady aerodynamic and weakly-coupled aeroelastic cases. Forced oscillation and unforced system experiments and computations have been compared for three configurations. This paper emphasizes interpretation of the experimental data, computational results and their comparisons from the perspective of validation of unsteady system predictions. The issues examined in detail are variability introduced by input choices for the computations, post-processing, and static aeroelastic modeling. The final issue addressed is interpreting unsteady information that is present in experimental data that is assumed to be steady, and the resulting consequences on the comparison data sets.
Gennaretti, M.; Testa, C.; Bernardini, G.
An aerodynamic/aeroacoustic solution methodology for predction of tonal noise emitted by helicopter rotors and propellers is presented. It is particularly suited for configurations dominated by localized, high-frequency inflow velocity fields as those generated by blade-vortex interactions. The unsteady pressure distributions are determined by the sectional, frequency-domain Küssner-Schwarz formulation, with downwash including the wake inflow velocity predicted by a three-dimensional, unsteady, panel-method formulation suited for the analysis of rotors operating in complex aerodynamic environments. The radiated noise is predicted through solution of the Ffowcs Williams-Hawkings equation. The proposed approach yields a computationally efficient solution procedure that may be particularly useful in preliminary design/multidisciplinary optimization applications. It is validated through comparisons with solutions that apply the airloads directly evaluated by the time-marching, panel-method formulation. The results are provided in terms of blade loads, noise signatures and sound pressure level contours. An estimation of the computational efficiency of the proposed solution process is also presented.
Rajkumar, T.; Bardina, Jorge; Clancy, Daniel (Technical Monitor)
Basic aerodynamic coefficients are modeled as functions of angles of attack and sideslip with vehicle lateral symmetry and compressibility effects. Most of the aerodynamic parameters can be well-fitted using polynomial functions. In this paper a fast, reliable way of predicting aerodynamic coefficients is produced using a neural network. The training data for the neural network is derived from wind tunnel test and numerical simulations. The coefficients of lift, drag, pitching moment are expressed as a function of alpha (angle of attack) and Mach number. The results produced from preliminary neural network analysis are very good.
Pinier, Jeremy T.
A novel method for implementing aerodynamic data dispersion analysis is herein introduced. A general mathematical approach combined with physical modeling tailored to the aerodynamic quantity of interest enables the generation of more realistically relevant dispersed data and, in turn, more reasonable flight simulation results. The method simultaneously allows for the aerodynamic quantities and their derivatives to be dispersed given a set of non-arbitrary constraints, which stresses the controls model in more ways than with the traditional bias up or down of the nominal data within the uncertainty bounds. The adoption and implementation of this new method within the NASA Ares I Crew Launch Vehicle Project has resulted in significant increases in predicted roll control authority, and lowered the induced risks for flight test operations. One direct impact on launch vehicles is a reduced size for auxiliary control systems, and the possibility of an increased payload. This technique has the potential of being applied to problems in multiple areas where nominal data together with uncertainties are used to produce simulations using Monte Carlo type random sampling methods. It is recommended that a tailored physics-based dispersion model be delivered with any aerodynamic product that includes nominal data and uncertainties, in order to make flight simulations more realistic and allow for leaner spacecraft designs.
Buning, Pieter G.; Wong, Tin-Chee; Dilley, Arthur D.; Pao, Jenn L.
The NASA X-43 "Hyper-X" hypersonic research vehicle will be boosted to a Mach 7 flight test condition mounted on the nose of an Orbital Sciences Pegasus launch vehicle. The separation of the research vehicle from the Pegasus presents some unique aerodynamic problems, for which computational fluid dynamics has played a role in the analysis. This paper describes the use of several CFD methods for investigating the aerodynamics of the research and launch vehicles in close proximity. Specifically addressed are unsteady effects, aerodynamic database extrapolation, and differences between wind tunnel and flight environments.
Carlson, John R.; Pao, S. Paul; Abdol-Hamid, Khaled S.; Jones, William T.
The multiblock three-dimensional Navier-Stokes method PAB3D was utilized by the Component Integration Branch (formerly Propulsion Aerodynamics Branch) at the NASA-Langley Research Center in an international study sponsored by AGARD Working Group #17 for the assessment of the state-of-the-art of propulsion-airframe integration testing techniques and CFD prediction technologies. Three test geometries from ONERA involving fundamental flow physics and four geometries from NASA-LaRC involving realistic flow interactions of wing, body, tail, and jet plumes were chosen by the Working Group. An overview of results on four (1 ONERA and 3 LaRC) of the seven test cases is presented. External static pressures, integrated pressure drag and total drag were calculated for the Langley test cases and jet plume velocity profiles and turbulent viscous stresses were calculated for the ONERA test case. Only selected data from these calculations are presented in this paper. The complete data sets calculated by the participants will be presented in an AGARD summary report. Predicted surface static pressures compared favorably with experimental data for the Langley geometries. Predicted afterbody drag compared well with experiment. Predicted nozzle drag was typically low due to over-compression of the flow near the trailing edge. Total drag was typically high. Predicted jet plume quantities on the ONERA case compared generally well with data.
Sullivan, Anders James
The design of missile body shapes often requires a compromise between aero-dynamic and electromagnetic performance goals. In general, the missile shape producing the lowest radar signature will be different from the preferred aero-dynamic shape. Interdisciplinary shape optimization is utilized to combine multiple disciplines to determine the best possible shape for a hybrid aerodynamic-electromagnetic problem. A composite missile body consisting of an axisymmetric body of revolution (BOR) and two thin flat plate attachments is considered. The goal is to minimize the drag and backscatter associated with this composite shape. The body is assumed to be perfectly conducting, and flying at zero degrees angle of attack. The variable nose shape serves as the optimization design parameter. To characterize the system performance, a cost function is defined which is comprised of weighted values of the drag and backscatter. To solve the electromagnetic problem, methods to treat electrically large complex bodies are investigated. Hybrid methods which combine the method of moments (MoM) with physical optics (PO) are developed to calculate the scattering from simple two-dimensional bodies. A surface-wave hybrid approach is shown to effectively approximate the traveling wave currents on the smooth interior portions of a BOR. Asymptotic methods are used to solve the resulting integral equations more efficiently. The hybrid methods are shown to produce MoM-quality results, while requiring less computational resources. To solve the composite body problem, an iterative technique is developed that preserves the simplicity of the original BOR scheme. In this formulation, the current over each part of the composite body is solved independently. The results from one part of the body are used to update the fields incident on the other part of the body. This procedure is repeated until the solution converges. To solve the aerodynamic problem, slender body theory is used to calculate the
Verdon, Joseph M.; Barnett, M.; Ayer, T. C.
Applications for unsteady aerodynamics analysis in this report are: (1) aeroelastic: blade flutter and forced vibration; (2) aeroacoustic: noise generation; (3) vibration and noise control; and (4) effects of unsteadiness on performance. This requires that the numerical simulations and analytical modeling be accurate and efficient and contain realistic operating conditions and arbitrary modes of unsteady excitation. The assumptions of this application contend that: (1) turbulence and transition can be modeled with the Reynolds averaged and using Navier-Stokes equations; (2) 'attached' flow with high Reynolds number will require thin-layer Navier-Stokes equations, or inviscid/viscid interaction analyses; (3) small-amplitude unsteady excitations will need nonlinear steady and linearized unsteady analyses; and (4) Re to infinity will concern inviscid flow. Several computer programs (LINFLO, CLT, UNSVIS, AND SFLOW-IVI) are utilized for these analyses. Results and computerized grid examples are shown. This report was given during NASA LeRC Workshop on Forced Response in Turbomachinery in August of 1993.
A simple method for an unsteady aerodynamic model tuning is proposed in this study. This method is based on the direct modification of the aerodynamic influence coefficient matrices. The aerostructures test wing 2 flight-test data is used to demonstrate the proposed model tuning method. The flutter speed margin computed using only the test validated structural dynamic model can be improved using the additional unsteady aerodynamic model tuning, and then the flutter speed margin requirement of 15 percent in military specifications can apply towards the test validated aeroelastic model. In this study, unsteady aerodynamic model tunings are performed at two time invariant flight conditions, at Mach numbers of 0.390 and 0.456. When the Mach number for the unsteady aerodynamic model tuning approaches to the measured fluttering Mach number, 0.502, at the flight altitude of 9,837 ft, the estimated flutter speed is approached to the measured flutter speed at this altitude. The minimum flutter speed difference between the estimated and measured flutter speed is -0.14 percent.
Abid, Ridha; Morrison, Joseph H.; Gatski, Thomas B.; Speziale, Charles G.
An explicit algebraic stress equation, developed by Gatski and Speziale, is used in the framework of K-epsilon formulation to predict complex aerodynamic turbulent flows. The nonequilibrium effects are modeled through coefficients that depend nonlinearly on both rotational and irrotational strains. The proposed model was implemented in the ISAAC Navier-Stokes code. Comparisons with the experimental data are presented which clearly demonstrate that explicit algebraic stress models can predict the correct response to nonequilibrium flow.
Beatty, T. D.; Worthey, M. K.
The V/STOL Aircraft Propulsive Effects (VAPE) computerized prediction method is evaluated. The program analyzes viscous effects, various jet, inlet, and Short TakeOff and Landing (STOL) models, and examines the aerodynamic configurations of V/STOL aircraft.
Wolfe, W.P.; Ochs, S.S.
An investigation was conducted into the capabilities and accuracy of a representative computational fluid dynamics code to predict the flow field and aerodynamic characteristics of typical wind-turbine airfoils. Comparisons of the computed pressure and aerodynamic coefficients were made with wind tunnel data. This work highlights two areas in CFD that require further investigation and development in order to enable accurate numerical simulations of flow about current generation wind-turbine airfoils: transition prediction and turbulence modeling. The results show that the laminar-to turbulent transition point must be modeled correctly to get accurate simulations for attached flow. Calculations also show that the standard turbulence model used in most commercial CFD codes, the k-e model, is not appropriate at angles of attack with flow separation. 14 refs., 28 figs., 4 tabs.
Mcconnaughey, H. V.; Griffin, L. W.
A preliminary assessment is made of two NASA-developed unsteady turbine stage computer codes. The methodology and previous partial validation of the codes are briefly outlined. Application of these codes to a Space Shuttle main engine turbine for two sets of operating conditions is then described. Steady and unsteady, two and three-dimensional results are presented, compared, and discussed. These results include time-mean and instantaneous airfoil pressure distributions and pressure fluctuations, streamlines on the airfoil surfaces and endwalls, and relative total pressure contours at different axial locations in the rotor passage. Although not available at the time of this writing, experimental data for one of the operating conditions simulated is forthcoming and will be used to assess the accuracy of the unsteady, as well as, the steady predictions presented. Issues related to code usage and resource requirements of the two codes are also discussed.
A simple method for an unsteady aerodynamic model tuning is proposed in this study. This method is based on the direct modification of the aerodynamic influence coefficient matrices. The aerostructures test wing 2 flight-test data is used to demonstrate the proposed model tuning method. The flutter speed margin computed using only the test validated structural dynamic model can be improved using the additional unsteady aerodynamic model tuning, and then the flutter speed margin requirement of 15 % in military specifications can apply towards the test validated aeroelastic model. In this study, unsteady aerodynamic model tunings are performed at two time invariant flight conditions, at Mach numbers of 0.390 and 0.456. When the Mach number for the unsteady model tuning approaches to the measured fluttering Mach number, 0.502, at the flight altitude of 9,837 ft, the estimated flutter speed is approached to the measured flutter speed at this altitude. The minimum flutter speed difference between the estimated and measured flutter speed is -.14 %.
Rajkumar, T.; Bardina, Jorge; Clancy, Daniel (Technical Monitor)
Wind tunnels use scale models to characterize aerodynamic coefficients, Wind tunnel testing can be slow and costly due to high personnel overhead and intensive power utilization. Although manual curve fitting can be done, it is highly efficient to use a neural network to define the complex relationship between variables. Numerical simulation of complex vehicles on the wide range of conditions required for flight simulation requires static and dynamic data. Static data at low Mach numbers and angles of attack may be obtained with simpler Euler codes. Static data of stalled vehicles where zones of flow separation are usually present at higher angles of attack require Navier-Stokes simulations which are costly due to the large processing time required to attain convergence. Preliminary dynamic data may be obtained with simpler methods based on correlations and vortex methods; however, accurate prediction of the dynamic coefficients requires complex and costly numerical simulations. A reliable and fast method of predicting complex aerodynamic coefficients for flight simulation I'S presented using a neural network. The training data for the neural network are derived from numerical simulations and wind-tunnel experiments. The aerodynamic coefficients are modeled as functions of the flow characteristics and the control surfaces of the vehicle. The basic coefficients of lift, drag and pitching moment are expressed as functions of angles of attack and Mach number. The modeled and training aerodynamic coefficients show good agreement. This method shows excellent potential for rapid development of aerodynamic models for flight simulation. Genetic Algorithms (GA) are used to optimize a previously built Artificial Neural Network (ANN) that reliably predicts aerodynamic coefficients. Results indicate that the GA provided an efficient method of optimizing the ANN model to predict aerodynamic coefficients. The reliability of the ANN using the GA includes prediction of aerodynamic
Romere, P. O.; Young, J. C.
An overview of the aerodynamic performance, longitudinal trim, and control surface effectiveness as developed from the second Shuttle flight is presented, along with comparisons with preflight predictions. The Orbiter enters the atmosphere at Mach 29 and descends through all the Mach speeds until touchdown. Control mechanisms for the Orbiter are outlined, noting the lessening of angle of attack as the terminal configuration is reached. Methods of developing the wind tunnel predictions are reviewed, including performance in different altitudes and speed regimes, and the impetus to decouple the flight control system from aerodynamics with respect to uncertainties. Preflight predictions agreed well with the lift/drag ratio above Mach 1, but higher L/D was exhibited than predicted below Mach 1. A deficiency in the wind tunnel environment was concluded for the inaccuracies in the predictions for hypersonic trim characteristics.
Maksimov, F. A.; Churakov, D. A.; Shevelev, Yu. D.
Complex-geometry design and grid generation are addressed. The gasdynamic equations are solved, and the numerical results are compared with experimental data. For aerodynamic problems, a suite of mathematical and information technology tools is proposed for the support and management of geometric models of actual objects. Based on the mathematical modeling methods developed, numerical experiments can be performed for a wide class of geometric forms and the aerodynamic properties of aircraft can be predicted with allowance for the viscosity effects.
Mendenhall, M.R. )
The present volume discusses tactical missile aerodynamic drag, drag-prediction methods for axisymmetric missile bodies, an aerodynamic heating analysis for supersonic missiles, a component buildup method for engineering analysis of missiles at low-to-high angles of attack, experimental and analytical methods for missiles with noncircular fuselages, and a vortex-cloud model for body vortex shedding and tracking. Also discussed are panel methods with vorticity effects and corrections for nonlinear compressibility, supersonic full-potential methods for missile body analysis, space-marching Euler solvers, the time-asymptotic Euler/Navier-Stokes methods for subsonic and transonic flows, 3D boundary layers on missiles, Navier-Stokes analyses of flows over slender airframes, and the interaction of exhaust plumes with missile airframes.
Jeng, D. R.; Keith, T. G.; Aliakbarkhanafjeh, A.
A new method for calculating the aerodynamic performance of horizontal axis wind turbines is described. The method, entitled the helical vortex method, directly calculates the local induced velocity due to helical vortices that originate at the rotor blade. Furthermore, the method does not require a specified circulation distribution. Results of the method are compared to similar results obtained from Wilson PROP code methods as well as to existing experimental data taken from a Mod-O wind turbine. It is shown that results of the proposed method agree well with experimental values of the power output both near cut-in and at rated wind speeds. Further, it is found that the method does not experience some of the numerical difficulties encountered by the PROP code when run at low wind velocities.
Gnoffo, Peter A.; Braun, Robert D.; Weilmuenster, K. James; Mitcheltree, Robert A.; Engelund, Walter C.; Powell, Richard W.
Postflight analysis of the Mars Pathfinder hypersonic, continuum aerodynamic data base is presented. Measured data include accelerations along the body axis and axis normal directions. Comparisons of preflight simulation and measurements show good agreement. The prediction of two static instabilities associated with movement of the sonic line from the shoulder to the nose and back was confirmed by measured normal accelerations. Reconstruction of atmospheric density during entry has an uncertainty directly proportional to the uncertainty in the predicted axial coefficient. The sensitivity of the moment coefficient to freestream density, kinetic models and center-of-gravity location are examined to provide additional consistency checks of the simulation with flight data. The atmospheric density as derived from axial coefficient and measured axial accelerations falls within the range required for sonic line shift and static stability transition as independently determined from normal accelerations.
Dillenius, Marnix F. E.; Johnson, David L.; Lesieutre, Daniel J.
The computer program DEMON3 was developed for the aerodynamic analysis of nonconventional supersonic configurations comprising a body with noncircular cross section and up to two wing or fin sections. Within a wing or fin section, the lifting surfaces may be cruciform, triform, planar, or low profile layouts; the planforms of the lifting surfaces allow for breaks in sweep. The body and fin sections are modeled by triplet and constant u-velocity panels, respectively, accounting for mutual body-fin interference. Fin thickness effects are included for the use of supersonic planar source panels. One of the unique features of DEMON3 is the modeling of high angle of attack vortical effects associated with the lifting surfaces and the body. In addition, shock expansion and Newtonian pressure calculation methods can be optionally engaged. These two dimensional nonlinear methods are augmented by aerodynamic interference determined from the linear panel methods. Depending on the geometric details of the body, the DEMON3 program can be used to analyze nonconventional configurations at angles of attack up to 25 degrees for Mach numbers from 1.1 to 6. Calculative results and comparisons with experimental data demonstrate the capabilities of DEMON3. Limitations and deficiencies are listed.
Caughey, David A.
Work has been directed at the development of efficient multigrid methods for the solution of aerodynamic problems involving complex geometries, including the development of computational methods for the solution of both inviscid and viscous transonic flow problems. The emphasis is on problems of complex, three-dimensional geometry. The methods developed are based upon finite-volume approximations to both the Euler and the Reynolds-Averaged Navier-Stokes equations. The methods are developed for use on multi-block grids using diagonalized implicit multigrid methods to achieve computational efficiency. The work is focused upon aerodynamic problems involving complex geometries, including advanced engine inlets.
Brentner, Kenneth S.
A great deal of progress has been made in the modeling of aerodynamically generated sound for rotors over the past decade. The Ffowcs Williams-Hawkings (FW-H ) equation has been the foundation for much of the development. Both subsonic and supersonic quadrupole noise formulations have been developed for the prediction of high-speed impulsive noise. In an effort to eliminate the need to compute the quadrupole contribution, the FW-H has also been utilized on permeable surfaces surrounding all physical noise sources. Comparison of the Kirchhoff formulation for moving surfaces with the FW-H equation have shown that the Kirchhoff formulation for moving surfaces can give erroneous results for aeroacoustic problems.
Tang, Deman; Kholodar, Denis; Juang, Jer-Nan; Dowell, Earl H.
The representation of unsteady aerodynamic flow fields in terms of global aerodynamic modes has proven to be a useful method for reducing the size of the aerodynamic model over those representations that use local variables at discrete grid points in the flow field. Eigenmodes and Proper Orthogonal Decomposition (POD) modes have been used for this purpose with good effect. This suggests that system identification models may also be used to represent the aerodynamic flow field. Implicit in the use of a systems identification technique is the notion that a relative small state space model can be useful in describing a dynamical system. The POD model is first used to show that indeed a reduced order model can be obtained from a much larger numerical aerodynamical model (the vortex lattice method is used for illustrative purposes) and the results from the POD and the system identification methods are then compared. For the example considered, the two methods are shown to give comparable results in terms of accuracy and reduced model size. The advantages and limitations of each approach are briefly discussed. Both appear promising and complementary in their characteristics.
Sirenko V.; Rohatgi U.
A small scale model (length 1710 mm) of General Motor SUV was built and tested in the wind tunnel for expected wind conditions and road clearance. Two passive devices, rear screen which is plate behind the car and rear fairing where the end of the car is aerodynamically extended, were incorporated in the model and tested in the wind tunnel for different wind conditions. The conclusion is that rear screen could reduce drag up to 6.5% and rear fairing can reduce the drag by 26%. There were additional tests for front edging and rear vortex generators. The results for drag reduction were mixed. It should be noted that there are aesthetic and practical considerations that may allow only partial implementation of these or any drag reduction options.
Peeters, M. F.
Five methods to increase the computational efficiency of aerodynamic design using numerical optimization, by reducing the computer time required to perform gradient calculations, are examined. The most promising method consists of drastically reducing the size of the computational domain on which aerodynamic calculations are made during gradient calculations. Since a gradient calculation requires the solution of the flow about an airfoil whose geometry was slightly perturbed from a base airfoil, the flow about the base airfoil is used to determine boundary conditions on the reduced computational domain. This method worked well in subcritical flow.
Wissink, Andrew M.; Lyrintzis, Anastasios S.; Strawn, Roger C.; Oliker, Leonid; Biswas, Rupak
This paper presents parallel implementations of two codes used in a combined CFD/Kirchhoff methodology to predict the aerodynamics and aeroacoustics properties of helicopters. The rotorcraft Navier-Stokes code, TURNS, computes the aerodynamic flowfield near the helicopter blades and the Kirchhoff acoustics code computes the noise in the far field, using the TURNS solution as input. The overall parallel strategy adds MPI message passing calls to the existing serial codes to allow for communication between processors. As a result, the total code modifications required for parallel execution are relatively small. The biggest bottleneck in running the TURNS code in parallel comes from the LU-SGS algorithm that solves the implicit system of equations. We use a new hybrid domain decomposition implementation of LU-SGS to obtain good parallel performance on the SP-2. TURNS demonstrates excellent parallel speedups for quasi-steady and unsteady three-dimensional calculations of a helicopter blade in forward flight. The execution rate attained by the code on 114 processors is six times faster than the same cases run on one processor of the Cray C-90. The parallel Kirchhoff code also shows excellent parallel speedups and fast execution rates. As a performance demonstration, unsteady acoustic pressures are computed at 1886 far-field observer locations for a sample acoustics problem. The calculation requires over two hundred hours of CPU time on one C-90 processor but takes only a few hours on 80 processors of the SP2. The resultant far-field acoustic field is analyzed with state of-the-art audio and video rendering of the propagating acoustic signals.
Griswold, M.; Roskam, J.
An analytical method is presented for predicting lateral-directional aerodynamic characteristics of light twin engine propeller-driven airplanes. This method is applied to the Advanced Technology Light Twin Engine airplane. The calculated characteristics are correlated against full-scale wind tunnel data. The method predicts the sideslip derivatives fairly well, although angle of attack variations are not well predicted. Spoiler performance was predicted somewhat high but was still reasonable. The rudder derivatives were not well predicted, in particular the effect of angle of attack. The predicted dynamic derivatives could not be correlated due to lack of experimental data.
Verdon, Joseph M.; Hall, Kenneth C.
A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.
Homicz, G. F.; Moselle, J. R.
A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.
Gebert, G. A.
Current aerodynamic interest has turned to the study of supermaneuverable fighters and weapon performance when launched in extreme flight conditions. The evaluation of design missile performance requires multiple runs of six degree-of-freedom (6-DOF) simulations, analyzing the missile behavior for a variety of launch and flight conditions. Before wind-tunnel tests, it is necessary to produce the aerodynamic loading of candidate missiles for 6-DOF analyses. Since semi-empirical formulas fail in regions of nonlinear aerodynamics, and solutions to the full Navier-Stokes equations are too costly and time consuming, an alternative method of discrete vortex analysis is re-examined. The present theory examines the three-dimensional nature of the shed vorticity and generalizes previous discrete vortex analyses. Consequently, the results demonstrate relative user independence in determining all slender-body loading at angles of attack from 0 to 70 deg. The rapid calculations of the discrete vortex method makes it a prime candidate for the determinations of high angle-of-attack aerodynamic databases.
Luo, JiaQi; Xiong, JunTao; Liu, Feng
This paper presents the fundamentals of a continuous adjoint method and the applications of this method to the aerodynamic design optimization of both external and internal flows. General formulation of the continuous adjoint equations and the corresponding boundary conditions are derived. With the adjoint method, the complete gradient information needed in the design optimization can be obtained by solving the governing flow equations and the corresponding adjoint equations only once for each cost function, regardless of the number of design parameters. An inverse design of airfoil is firstly performed to study the accuracy of the adjoint gradient and the effectiveness of the adjoint method as an inverse design method. Then the method is used to perform a series of single and multiple point design optimization problems involving the drag reduction of airfoil, wing, and wing-body configuration, and the aerodynamic performance improvement of turbine and compressor blade rows. The results demonstrate that the continuous adjoint method can efficiently and significantly improve the aerodynamic performance of the design in a shape optimization problem.
Ibrahim, A. H.; Hou, G. J.-W.; Tiwari, S. N. (Principal Investigator)
Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite
Keys, C. N.; Stephniewski, W. Z. (Editor)
Application of theories, as well as, special methods of procedures applicable to performance prediction are illustrated first, on an example of the conventional helicopter and then, winged and tandem configurations. Performance prediction of conventional helicopters in hover and vertical ascent are investigated. Various approaches to performance prediction in forward translation are presented. Performance problems are discussed only this time, a wing is added to the baseline configuration, and both aircraft are compared with respect to their performance. This comparison is extended to a tandem. Appendices on methods for estimating performance guarantees and growth of aircraft concludes this volume.
Mercer, J. E.; Weber, J. A.; Lesferd, E. P.
A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.
Mercer, J. E.; Weber, J. A.; Lesferd, E. P.
A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.
Korsmeyer, David (Technical Monitor); Rajkumar, T.; Bardina, Jorge
Basic aerodynamic coefficients are modeled as functions of angle of attack, speed brake deflection angle, Mach number, and side slip angle. Most of the aerodynamic parameters can be well-fitted using polynomial functions. We previously demonstrated that a neural network is a fast, reliable way of predicting aerodynamic coefficients. We encountered few under fitted and/or over fitted results during prediction. The training data for the neural network are derived from wind tunnel test measurements and numerical simulations. The basic questions that arise are: how many training data points are required to produce an efficient neural network prediction, and which type of transfer functions should be used between the input-hidden layer and hidden-output layer. In this paper, a comparative study of the efficiency of neural network prediction based on different transfer functions and training dataset sizes is presented. The results of the neural network prediction reflect the sensitivity of the architecture, transfer functions, and training dataset size.
Mendenhall, M. R.; Spangler, S. B.
A theoretical method has been developed to predict the longitudinal aerodynamic characteristics of engine-wing-flap combinations with externally blown flaps (EBF) and upper surface blowing (USB) high lift devices. Potential flow models of the lifting surfaces and the jet wake are combined to calculate the induced interference of the engine wakes on the lifting surfaces. The engine wakes may be circular, elliptic, or rectangular cross-sectional jets, and the lifting surfaces are comprised of a wing with multiple-slotted trailing-edge flaps or a deflected trailing-edge Coanda surface. Results are presented showing comparisons of measured and predicted forces, pitching moments, span-load distributions, and flow fields.
Kamel, Ouari; Mohand, Ouhrouche; Toufik, Rekioua; Taib, Nabil
In order to improvement of the performances for wind energy conversions systems (WECS), an advanced control techniques must be used. In this paper, as an alternative to conventional PI-type control methods, a nonlinear predictive control (NPC) approach is developed for DFIG-based wind turbine. To enhance the robustness of the controller, a disturbance observer is designed to estimate the aerodynamic torque which is considered as an unknown perturbation. An explicitly analytical form of the optimal predictive controller is given consequently on-line optimization is not necessary The DFIG is fed through the rotor windings by a back-to-back converter controlled by Pulse Width Modulation (PWM), where the stator winding is directly connected to the grid. The presented simulation results show a good performance in trajectory tracking of the proposed strategy and rejection of disturbances is successfully achieved.
Chi, x.; Li, Y.; Chen, H.; Addy, H. E.; Choo, Y. K.; Shih, T. I-P.
WIND, Fluent, and PowerFLOW were used to predict the lift, drag, and moment coefficients of a business-jet airfoil with a rime ice (rough and jagged, but no protruding horns) and with a glaze ice (rough and jagged end has two or more protruding horns) for angles of attack from zero to and after stall. The performance of the following turbulence models were examined by comparing predictions with available experimental data. Spalart-Allmaras (S-A), RNG k-epsilon, shear-stress transport, v(sup 2)-f, and a differential Reynolds stress model with and without non-equilibrium wall functions. For steady RANS simulations, WIND and FLUENT were found to give nearly identical results if the grid about the iced airfoil, the turbulence model, and the order of accuracy of the numerical schemes used are the same. The use of wall functions was found to be acceptable for the rime ice configuration and the flow conditions examined. For rime ice, the S-A model was found to predict accurately until near the stall angle. For glaze ice, the CFD predictions were much less satisfactory for all turbulence models and codes investigated because of the large separated region produced by the horns. For unsteady RANS, WIND and FLUENT did not provide better results. PowerFLOW, based on the Lattice Boltzmann method, gave excellent results for the lift coefficient at and near stall for the rime ice, where the flow is inherently unsteady.
Egolf, T. Alan; Anderson, Olof L.; Edwards, David E.; Landgrebe, Anton J.
A user's manual for the computer program developed for the prediction of propeller-nacelle aerodynamic performance reported in, An Analysis for High Speed Propeller-Nacelle Aerodynamic Performance Prediction: Volume 1 -- Theory and Application, is presented. The manual describes the computer program mode of operation requirements, input structure, input data requirements and the program output. In addition, it provides the user with documentation of the internal program structure and the software used in the computer program as it relates to the theory presented in Volume 1. Sample input data setups are provided along with selected printout of the program output for one of the sample setups.
Hoyniak, D.; Fleeter, S.
A mathematical model is developed to predict the enhanced coupled bending-torsion unstalled supersonic flutter stability due to alternate circumferential spacing aerodynamic detuning of a turbomachine rotor. The translational and torsional unsteady aerodynamic coefficients are developed in terms of influence coefficients, with the coupled bending-torsion stability analysis developed by considering the coupled equations of motion together with the unsteady aerodynamic loading. The effect of this aerodynamic detuning on coupled bending-torsion unstalled supersonic flutter as well as the verification of the modeling are then demonstrated by considering an unstable 12 bladed rotor, with Verdon's uniformly spaced Cascade B flow geometry as a baseline. However, with the elastic axis and center of gravity at 60 percent of the chord, this type of aerodynamic detuning has a minimal effect on stability. For both uniform and nonuniform circumferentially space rotors, a single degree of freedom torsion mode analysis was shown to be appropriate for values of the bending-torsion natural frequency ratio lower than 0.6 and higher 1.2. When the elastic axis and center of gravity are not coincident, the effect of detuning on cascade stability was found to be very sensitive to the location of the center of gravity with respect to the elastic axis. In addition, it was determined that when the center of gravity was forward of an elastic axis located at midchord, a single degree of freedom torsion model did not accurately predict cascade stability.
Dowell, E. H.
A simple, direct procedure was developed for converting frequency domain aerodynamics into indicial aerodynamics. The data required for aerodynamic forces in the frequency domain may be obtained from any available (linear) theory. The method retains flexibility for the analyst and is based upon the particular character of the frequency domain results. An evaluation of the method was made for incompressible, subsonic, and transonic two dimensional flows.
Faller, W E; Schreck, S J
The capability to control unsteady separated flow fields could dramatically enhance aircraft agility. To enable control, however, real-time prediction of these flow fields over a broad parameter range must be realized. The present work describes real-time predictions of three-dimensional unsteady separated flow fields and aerodynamic coefficients using neural networks. Unsteady surface-pressure readings were obtained from an airfoil pitched at a constant rate through the static stall angle. All data sets were comprised of 15 simultaneously acquired pressure records and one pitch angle record. Five such records and the associated pitch angle histories were used to train the neural network using a time-series algorithm. Post-training, the input to the network was the pitch angle (alpha), the angular velocity (dalpha/dt), and the initial 15 recorded surface pressures at time (t (0)). Subsequently, the time (t+Deltat) network predictions, for each of the surface pressures, were fed back as the input to the network throughout the pitch history. The results indicated that the neural network accurately predicted the unsteady separated flow fields as well as the aerodynamic coefficients to within 5% of the experimental data. Consistent results were obtained both for the training set as well as for generalization to both other constant pitch rates and to sinusoidal pitch motions. The results clearly indicated that the neural-network model could predict the unsteady surface-pressure distributions and aerodynamic coefficients based solely on angle of attack information. The capability for real-time prediction of both unsteady separated flow fields and aerodynamic coefficients across a wide range of parameters in turn provides a critical step towards the development of control systems targeted at exploiting unsteady aerodynamics for aircraft manoeuvrability enhancement. PMID:18263439
Potsdam, Mark; Yeo, Hyeonsoo; Johnson, Wayne
This work couples a computational fluid dynamics (CFD) code and rotorcraft computational structural dynamics (CSD) code to calculate helicopter rotor airloads across a range of flight conditions. An iterative loose (weak) coupling methodology is used to couple the CFD and CSD codes on a per revolution, periodic basis. The CFD uses a high fidelity, Navier-Stokes, overset grid methodology with first principles-based wake capturing. Modifications are made to the CFD code for aeroelastic analysis. For a UH-60A Blackhawk helicopter, four challenging level flight conditions are computed: 1) low speed (u = 0.15) with blade-vortex interaction, 2) high speed (u = 0.37) with advancing blade negative lift, 3) high thrust with dynamic stall (u = 0.24), and 4) hover. Results are compared with UH-60A Airloads Program fight test data. Most importantly, for all cases the loose coupling methodology is shown to be stable, convergent, and robust with full coupling of normal force, pitching moment, and chord force. In comparison with flight test data, normal force and pitching moment magnitudes are in good agreement. For the high speed and dynamic stall cases a phase lag in comparison with the data is seen, nonetheless, the shapes of the curves are very good. Overall, the results are noteworthy improvement over lifting line aerodynamics used in rotorcraft comprehensive codes.
Murman, Scott M.; Aftosmis, Michael J.; Rogers, Stuart E.
An automated Computational Fluid Dynamics process for determining the aerodynamic Characteristics of debris shedding from the Space Shuttle Launch Vehicle during ascent is presented. This process uses Cartesian fully-coupled, six-degree-of-freedom simulations of isolated debris pieces in a Monte Carlo fashion to produce models for the drag and crossrange behavior over a range of debris shapes and shedding scenarios. A validation of the Cartesian methods against ballistic range data for insulating foam debris shapes at flight conditions, as well as validation of the resulting models, are both contained. These models are integrated with the existing shuttle debris transport analysis software to provide an accurate and efficient engineering tool for analyzing debris sources and their potential for damage.
Rajkumar, T.; Aragon, Cecilia; Bardina, Jorge; Britten, Roy
A fast, reliable way of predicting aerodynamic coefficients is produced using a neural network optimized by a genetic algorithm. Basic aerodynamic coefficients (e.g. lift, drag, pitching moment) are modelled as functions of angle of attack and Mach number. The neural network is first trained on a relatively rich set of data from wind tunnel tests of numerical simulations to learn an overall model. Most of the aerodynamic parameters can be well-fitted using polynomial functions. A new set of data, which can be relatively sparse, is then supplied to the network to produce a new model consistent with the previous model and the new data. Because the new model interpolates realistically between the sparse test data points, it is suitable for use in piloted simulations. The genetic algorithm is used to choose a neural network architecture to give best results, avoiding over-and under-fitting of the test data.
Rizzi, Arthur; Lucking, James M.
The CAWAPI-2 coordinated project has been underway to improve CFD predictions of slender airframe aerodynamics. The work is focused on two flow conditions and leverages a unique flight data set obtained with the F-16XL aircraft for comparison and verification. These conditions, a low-speed high angle-of-attack case and a transonic low angle-of-attack case, were selected from a prior prediction campaign wherein the CFD failed to provide acceptable results. In re-visiting these two cases, approaches for improved results include better, denser grids using more grid adaptation to local flow features as well as unsteady higher-fidelity physical modeling like hybrid RANS/URANS-LES methods. The work embodies predictions from multiple numerical formulations that are contributed from multiple organizations where some authors investigate other possible factors that could explain the discrepancies in agreement, e.g. effects due to deflected control surfaces during the flight tests, as well as static aeroelastic deflection of the outer wing. This paper presents the synthesis of all the results and findings and draws some conclusions that lead to an improved understanding of the underlying flow physics, and finally making the connections between the physics and aircraft features.
Maskew, B.; Strash, D.; Nathman, J.; Dvorak, F. A.
A computer program, VSAERO, has been applied to a number of V/STOL configurations with a view to advancing prediction techniques for the low-speed aerodynamic characteristics. The program couples a low-order panel method with surface streamline calculation and integral boundary layer procedures. The panel method--which uses piecewise constant source and doublet panels-includes an iterative procedure for wake shape and models boundary layer displacement effect using the source transpiration technique. Certain improvements to a basic vortex tube jet model were installed in the code prior to evaluation. Very promising results were obtained for surface pressures near a jet issuing at 90 deg from a flat plate. A solid core model was used in the initial part of the jet with a simple entrainment model. Preliminary representation of the downstream separation zone significantly improve the correlation. The program accurately predicted the pressure distribution inside the inlet on the Grumman 698-411 design at a range of flight conditions. Furthermore, coupled viscous/potential flow calculations gave very close correlation with experimentally determined operational boundaries dictated by the onset of separation inside the inlet. Experimentally observed degradation of these operational boundaries between nacelle-alone tests and tests on the full configuration were also indicated by the calculation. Application of the program to the General Dynamics STOL fighter design were equally encouraging. Very close agreement was observed between experiment and calculation for the effects of power on pressure distribution, lift and lift curve slope.
Sim, B. W.; Lim, J. W.
Predictions of blade-vortex interaction (BVI) noise, using blade airloads obtained from a coupled aerodynamic and structural methodology, are presented. This methodology uses an iterative, loosely-coupled trim strategy to cycle information between the OVERFLOW-2 (CFD) and CAMRAD-II (CSD) codes. Results are compared to the HART-II baseline, minimum noise and minimum vibration conditions. It is shown that this CFD/CSD state-of-the-art approach is able to capture blade airload and noise radiation characteristics associated with BVI. With the exception of the HART-II minimum noise condition, predicted advancing and retreating side BVI for the baseline and minimum vibration conditions agrees favorably with measured data. Although the BVI airloads and noise amplitudes are generally under-predicted, this CFD/CSD methodology provides an overall noteworthy improvement over the lifting line aerodynamics and free-wake models typically used in CSD comprehensive analysis codes.
Young, J. C.; Perez, L. F.; Romere, P. O.; Kanipe, D. B.
Results of comparisons of predictions of aerodynamic performance, longitudinal trim, and reaction control jet interaction with data from the initial Shuttle flight are presented. The Shuttle's control surfaces are described, and it is noted that the flight plan contained no provisions for maneuvering capability tests. Wind tunnel testing totaling 35,000 hr were used to replace graduated flight testing, and calculations were made to allow for nonsimulated structural deformation, flowfield parameters, and profile drag. The goal was to desensitize the flight control system with respect to the aerodynamics by adding variations to the predictions. Lift/drag agreed well above Mach 1, while lower drag was encountered below Mach 1. Trim characteristics were predicted accurately between Mach 2-10, and less than satisfactorily outside that range. Discrepancies were also observed for jet interaction effects for the aft yaw jets at Mach numbers greater than 10.
Riley, Donald C.
This paper contains a collection of some results of four individual studies presenting calculated numerical values for airfoil aerodynamic stability derivatives in unseparated inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory motion were considered. Variables include the number of vortices representing the airfoil, the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and circular or parabolic arc camber. Comparisons with some experimental and other theoretical information are included. The calculated aerodynamic numerical results obtained using a limited number of vortices provided in each study compared favorably with thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated herein (such as induced drag) that are not readily available elsewhere.
Hollis, Brian R.; Thompson, Richard A.; Berry, Scott A.; Horvath, Thomas J.; Murphy, Kelly J.; Nowak, Robert J.; Alter, Stephen J.
This report details a computational fluid dynamics study conducted in support of the phase II development of the X-33 vehicle. Aerodynamic and aeroheating predictions were generated for the X-33 vehicle at both flight and wind-tunnel test conditions using two finite-volume, Navier-Stokes solvers. Aerodynamic computations were performed at Mach 6 and Mach 10 wind-tunnel conditions for angles of attack from 10 to 50 with body-flap deflections of 0 to 20. Additional aerodynamic computations were performed over a parametric range of free-stream conditions at Mach numbers of 4 to 10 and angles of attack from 10 to 50. Laminar and turbulent wind-tunnel aeroheating computations were performed at Mach 6 for angles of attack of 20 to 40 with body-flap deflections of 0 to 20. Aeroheating computations were performed at four flight conditions with Mach numbers of 6.6 to 8.9 and angles of attack of 10 to 40. Surface heating and pressure distributions, surface streamlines, flow field information, and aerodynamic coefficients from these computations are presented, and comparisons are made with wind-tunnel data.
Maskew, B.; Rao, B. M.; Dvorak, F. A.
The development of a simple yet effective technique for modelling the effects of trailing edge separation is discussed. The model encloses the low energy region with free vortex sheets coupled with a potential flow panel method. The technique includes an interation cycle between viscous and potential flow routines and its development from the two dimensional case to the three dimensional case is discussed. A description of the potential flow panel method, which is based on an internal Dirichlet boundary condition, is included.
Rubbert, P. E.
A small-disturbance theory is developed for predicting the aerodynamics of an airplane in sideslip. Second-order terms involving the interaction between sideslip angle and angle of attack, sideslip angle and wing camber, etc., are retained. It is found that the second-order terms can produce the dominant sideslip effects when the dihedral of the lifting surfaces is small. Numerical implementation of the theory requires a solution procedure capable of producing accurate velocity gradients in the first-order solution.
Larson, R. S.; Nelson, D. P.; Stevens, B. S.
Five co-annular nozzle models, covering a systematic variation of nozzle geometry, were tested statically over a range of exhaust conditions including inverted velocity profile (IVP) (fan to primary stream velocity ratio 1) and non IVP profiles. Fan nozzle pressure ratio (FNPR) was varied from 1.3 to 4.1 at primary nozzle pressure ratios (PNPR) of 1.53 and 2.0. Fan stream temperatures of 700 K (1260 deg R) and 1089 K(1960 deg R) were tested with primary stream temperatures of 700 K (1260 deg R), 811 K (1460 deg R), and 1089 K (1960 deg R). At fan and primary stream velocities of 610 and 427 m/sec (2000 and 1400 ft/sec), respectively, increasing fan radius ratio from 0.69 to 0.83 reduced peak perceived noise level (PNL) 3 dB, and an increase in primary radius ratio from 0 to 0.81 (fan radius ratio constant at 0.83) reduced peak PNL an additional 1.0 dB. There were no noise reductions at a fan stream velocity of 853 m/sec (2800 ft/sec). Increasing fan radius ratio from 0.69 to 0.83 reduced nozzle thrust coefficient 1.2 to 1.5% at a PNPR of 1.53, and 1.7 to 2.0% at a PNPR of 2.0. The developed acoustic prediction procedure collapsed the existing data with standard deviation varying from + or - 8 dB to + or - 7 dB. The aerodynamic performance prediction procedure collapsed thrust coefficient measurements to within + or - .004 at a FNPR of 4.0 and a PNPR of 2.0.
Borggaard, Jeff; Burns, John
The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.
Sawyer, W. C.; Allen, J. M.; Hernandez, G.; Dillenius, M. F. E.; Hemsch, M. J.
This paper presents a survey of engineering computational methods and experimental programs used for estimating the aerodynamic characteristics of missile configurations. Emphasis is placed on those methods which are suitable for preliminary design of conventional and advanced concepts. An analysis of the technical approaches of the various methods is made in order to assess their suitability to estimate longitudinal and/or lateral-directional characteristics for different classes of missile configurations. Some comparisons between the predicted characteristics and experimental data are presented. These comparisons are made for a large variation in flow conditions and model attitude parameters. The paper also presents known experimental research programs developed for the specific purpose of validating analytical methods and extending the capability of data-base programs.
Romere, P. O.; Whitnah, A. M.
The Space Shuttle orbiter flight test program has required the aerodynamicist to take a new approach in determining flight characteristics. A conventional flight test program of slowly and cautiously approaching more severe flight conditions was not possible with the orbiter. On the first four orbital flights, The orbiter entered the atmosphere at approximately Mach 29 and decelerated through the Mach range (the subsonic portion of flight had also been flown by another orbiter vehicle during the Approach and Landing Test Program). Certification for these flights was achieved by an extensive wind tunnel test and analysis program. The initial series of flights of the orbiter were heavily instrumented for the purpose of obtaining accurate aerodynamic data. The flight data derived from the entry Mach range provided comparisons between flight and wind tunnel derived predicted data in the areas of both aerodynamic performance and longitudinal trim.
Carlson, Leland A.
Research conducted during the period from July 1991 through December 1992 is covered. A method based upon the quasi-analytical approach was developed for computing the aerodynamic sensitivity coefficients of three dimensional wings in transonic and subsonic flow. In addition, the method computes for comparison purposes the aerodynamic sensitivity coefficients using the finite difference approach. The accuracy and validity of the methods are currently under investigation.
Adenan, Nur Hamiza; Md Noorani, Mohd Salmi
The analysis and prediction of time series data have been addressed by researchers. Many techniques have been developed to be applied in various areas, such as weather forecasting, financial markets and hydrological phenomena involving data that are contaminated by noise. Therefore, various techniques to improve the method have been introduced to analyze and predict time series data. In respect of the importance of analysis and the accuracy of the prediction result, a study was undertaken to test the effectiveness of the improved nonlinear prediction method for data that contain noise. The improved nonlinear prediction method involves the formation of composite serial data based on the successive differences of the time series. Then, the phase space reconstruction was performed on the composite data (one-dimensional) to reconstruct a number of space dimensions. Finally the local linear approximation method was employed to make a prediction based on the phase space. This improved method was tested with data series Logistics that contain 0%, 5%, 10%, 20% and 30% of noise. The results show that by using the improved method, the predictions were found to be in close agreement with the observed ones. The correlation coefficient was close to one when the improved method was applied on data with up to 10% noise. Thus, an improvement to analyze data with noise without involving any noise reduction method was introduced to predict the time series data.
Margason, Richard; Kuhn, Richard
Available prediction methods applied to problems of aero/propulsion interactions for short takeoff and vertical landing (STOVL) aircraft are critically reviewed and an assessment of their strengths and weaknesses provided. The first two problems deal with aerodynamic performance effects during hover: (1) out-of-ground effect, and (2) in-ground effect. The first can be evaluated for some multijet cases; however, the second problem is very difficult to evaluate for multijets. The ground-environment effects due to wall jets and fountain flows directly affect hover performance. In a related problem: (3) hot-gas ingestion affects the engine operation. Both of these problems as well as jet noise affect the ability of people to work near the aircraft and the ability of the aircraft to operate near the ground. Additional problems are: (4) the power-augmented lift due to jet-flap effects (both in- and out-of-ground effects), and (5) the direct jet-lift effects during short takeoff and landing (STOL) operations. The final problem: (6) is the aerodynamic/propulsion interactions in transition between hover and wing-borne flight. Areas where modern CFD methods can provide improvements to current computational capabilities are identified.
Petot, D.; Loiseau, H.
Unsteady aerodynamic methods adopted for the study of aeroelasticity in helicopters are considered with focus on the development of a semiempirical model of unsteady aerodynamic forces acting on an oscillating profile at high incidence. The successive smoothing algorithm described leads to the model's coefficients in a very satisfactory manner.
Hennes, Christopher C.
Many computational tools are used when developing a modern helicopter. As the design space is narrowed, more accurate and time-intensive tools are brought to bear. These tools are used to determine the effect of a design decision on the performance, handling, stability and efficiency of the aircraft. One notable parameter left out of this process is acoustics. This is due in part to the difficulty in making useful acoustics calculations that reveal the differences between various design configurations. This thesis presents a new approach designed to bridge the gap in prediction capability between fast but low-fidelity Lagrangian particle methods, and slow but high-fidelity Eulerian computational fluid dynamics simulations. A multi-pronged approach is presented. First, a simple flow solver using well-understood and tested flow solution methodologies is developed specifically to handle bodies in arbitrary motion. To this basic flow solver two new technologies are added. The first is an Immersed Boundary technique designed to be tolerant of geometric degeneracies and low-resolution grids. This new technique allows easy inclusion of complex fuselage geometries at minimal computational cost, improving the ability of a solver to capture the complex interactional aerodynamic effects expected in modern rotorcraft design. The second new technique is an extension of a concept from flow visualization where the motion of tip vortices are tracked through the solution using massless particles convecting with the local flow. In this extension of that concept, the particles maintain knowledge of the expected and actual vortex strength. As a post-processing step, when the acoustic calculations are made, these particles are used to augment the loading noise calculation and reproduce the highly-impulsive character of blade-vortex interaction noise. In combination these new techniques yield a significant improvement to the state of the art in rotorcraft blade-vortex interaction noise
Zhao, Lin; Ge, Yaojun
This study proposes a new identification algorithm about the admittance function, which can estimate the full set of six aerodynamic admittance functions considering cross power spectral density functions about the forces and the turbulence components. The method was first numerically validated through Monte Carlo simulations, and then adopted to estimate the aerodynamic admittance of a streamlined bridge deck. The identification method was further validated through a comparison between the numerical calculation and wind tunnel tests on a moving bridge section.
Berry, John D.; Letnikov, Victor; Bavykina, Irena; Chaffin, Mark S.
Recent advances in computing subsonic flow have been applied to helicopter configurations with various degrees of success. This paper is a comparison of two specific methods applied to a particularly challenging regime of helicopter flight, very low speeds, where the interaction of the rotor wake and the fuselage are most significant. Comparisons are made between different methods of predicting the interactional aerodynamics associated with a simple generic helicopter configuration. These comparisons are made using fuselage pressure data from a Mach-scaled powered model helicopter with a rotor diameter of approximately 3 meters. The data shown are for an advance ratio of 0.05 with a thrust coefficient of 0.0066. The results of this comparison show that in this type of complex flow both analytical techniques have regions where they are more accurate in matching the experimental data.
Goldman, L. J.
The prediction of the cooled aerodynamic performance, for both stators and turbines, at actual primary to coolant inlet total temperature ratios from the results obtained at a reduced total temperature ratio is described. Theoretical and available experimental results were compared for convection film and transpiration cooled stator vanes and for a film cooled, single stage core turbine. For these tests the total temperature ratio varied from near 1.0 to about 2.7. The agreement between the theoretical and the experimental results was, in general, reasonable.
Pamadi, Bandu N.; Taylor, Larry W.; Leary, Terrance O.
A method is proposed for the reduction of the aerodynamic drag of bluff bodies, particularly for application to road transport vehicles. This technique consists of installation of panels on the forward surface of the vehicle facing the airstream. With the help of road tests, it was demonstrated that the attachment of proposed panels can reduce aerodynamic drag of road vehicles and result in significant fuel cost savings and conservation of energy resources.
Punta, Marco; Forrest, Lucy R.; Bigelow, Henry; Kernytsky, Andrew; Liu, Jinfeng; Rost, Burkhard
We survey computational approaches that tackle membrane protein structure and function prediction. While describing the main ideas that have led to the development of the most relevant and novel methods, we also discuss pitfalls, provide practical hints and highlight the challenges that remain. The methods covered include: sequence alignment, motif search, functional residue identification, transmembrane segment and protein topology predictions, homology and ab initio modeling. Overall, predictions of functional and structural features of membrane proteins are improving, although progress is hampered by the limited amount of high-resolution experimental information available. While predictions of transmembrane segments and protein topology rank among the most accurate methods in computational biology, more attention and effort will be required in the future to ameliorate database search, homology and ab initio modeling. PMID:17367718
Luckring, James M.; Rizzi, Arthur; Davis, M. Bruce
A coordinated project has been underway to improve CFD predictions of slender airframe aerodynamics. The work is focused on two flow conditions and leverages a unique flight data set obtained with an F-16XL aircraft. These conditions, a low-speed high angleof- attack case and a transonic low angle-of-attack case, were selected from a prior prediction campaign wherein the CFD failed to provide acceptable results. In this paper the background, objectives and approach to the current project are presented. The work embodies predictions from multiple numerical formulations that are contributed from multiple organizations, and the context of this campaign to other multi-code, multiorganizational efforts is included. The relevance of this body of work toward future supersonic commercial transport concepts is also briefly addressed.
Mohrfeld-Halterman, J. A.; Uddin, M.
We described in this paper the development of a high fidelity vehicle aerodynamic model to fit wind tunnel test data over a wide range of vehicle orientations. We also present a comparison between the effects of this proposed model and a conventional quasi steady-state aerodynamic model on race vehicle simulation results. This is done by implementing both of these models independently in multi-body quasi steady-state simulations to determine the effects of the high fidelity aerodynamic model on race vehicle performance metrics. The quasi steady state vehicle simulation is developed with a multi-body NASCAR Truck vehicle model, and simulations are conducted for three different types of NASCAR race tracks, a short track, a one and a half mile intermediate track, and a higher speed, two mile intermediate race track. For each track simulation, the effects of the aerodynamic model on handling, maximum corner speed, and drive force metrics are analysed. The accuracy of the high-fidelity model is shown to reduce the aerodynamic model error relative to the conventional aerodynamic model, and the increased accuracy of the high fidelity aerodynamic model is found to have realisable effects on the performance metric predictions on the intermediate tracks resulting from the quasi steady-state simulation.
Egolf, T. Alan; Anderson, Olof L.; Edwards, David E.; Landgrebe, Anton J.
A computer program, the Propeller Nacelle Aerodynamic Performance Prediction Analysis (PANPER), was developed for the prediction and analysis of the performance and airflow of propeller-nacelle configurations operating over a forward speed range inclusive of high speed flight typical of recent propfan designs. A propeller lifting line, wake program was combined with a compressible, viscous center body interaction program, originally developed for diffusers, to compute the propeller-nacelle flow field, blade loading distribution, propeller performance, and the nacelle forebody pressure and viscous drag distributions. The computer analysis is applicable to single and coaxial counterrotating propellers. The blade geometries can include spanwise variations in sweep, droop, taper, thickness, and airfoil section type. In the coaxial mode of operation the analysis can treat both equal and unequal blade number and rotational speeds on the propeller disks. The nacelle portion of the analysis can treat both free air and tunnel wall configurations including wall bleed. The analysis was applied to many different sets of flight conditions using selected aerodynamic modeling options. The influence of different propeller nacelle-tunnel wall configurations was studied. Comparisons with available test data for both single and coaxial propeller configurations are presented along with a discussion of the results.
VanZante, Dale E.; Koch, L. Danielle; Wernet, Mark P.; Podboy, Gary G.
Driven by the need for low production costs, electronics cooling fans have evolved differently than the bladed components of gas turbine engines which incorporate multiple technologies to enhance performance and durability while reducing noise emissions. Drawing upon NASA Glenn's experience in the measurement and prediction of gas turbine engine aeroacoustic performance, tests have been conducted to determine if these tools and techniques can be extended for application to the aerodynamics and acoustics of electronics cooling fans. An automated fan plenum installed in NASA Glenn's Acoustical Testing Laboratory was used to map the overall aerodynamic and acoustic performance of a spaceflight qualified 80 mm diameter axial cooling fan. In order to more accurately identify noise sources, diagnose performance limiting aerodynamic deficiencies, and validate noise prediction codes, additional aerodynamic measurements were recorded for two operating points: free delivery and a mild stall condition. Non-uniformities in the fan s inlet and exhaust regions captured by Particle Image Velocimetry measurements, and rotor blade wakes characterized by hot wire anemometry measurements provide some assessment of the fan aerodynamic performance. The data can be used to identify fan installation/design changes which could enlarge the stable operating region for the fan and improve its aerodynamic performance and reduce noise emissions.
Riley, Christopher J.; DeJarnette, Fred R.
A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.
Buffum, Daniel H.; Fleeter, Sanford
The steady and unsteady aerodynamics of a linear oscillating cascade are investigated using experimental and computational methods. Experiments are performed to quantify the torsion mode oscillating cascade aerodynamics of the NASA Lewis Transonic Oscillating Cascade for subsonic inlet flowfields using two methods: simultaneous oscillation of all the cascaded airfoils at various values of interblade phase angle, and the unsteady aerodynamic influence coefficient technique. Analysis of these data and correlation with classical linearized unsteady aerodynamic analysis predictions indicate that the wind tunnel walls enclosing the cascade have, in some cases, a detrimental effect on the cascade unsteady aerodynamics. An Euler code for oscillating cascade aerodynamics is modified to incorporate improved upstream and downstream boundary conditions and also the unsteady aerodynamic influence coefficient technique. The new boundary conditions are shown to improve the unsteady aerodynamic influence coefficient technique. The new boundary conditions are shown to improve the unsteady aerodynamic predictions of the code, and the computational unsteady aerodynamic influence coefficient technique is shown to be a viable alternative for calculation of oscillating cascade aerodynamics.
Griffin, L. W.; Mcconnaughey, H. V.
A numerical study of the aerodynamic and thermal environment associated with axial turbine stages is presented. Computations were performed using a modification of the unsteady viscous code, ROTORI, and an improved version of the steady inviscid cascade system, MERIDL-TSONIC, coupled with boundary layer codes, BLAYER and STAN5. Two different turbine stages were analyzed: the first stage of the United Technologies Research Center Large Scale Rotating Rig (LSRR) and the first stage of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump turbine. The time-averaged airfoil midspan pressure and heat transfer profiles were predicted for numerous thermal boundary conditions including adiabatic wall, prescribed surface temperature, and prescribed heat flux. Computed solutions are compared with each other and with experimental data in the case of the LSRR calculations. Modified ROTORI predictions of unsteady pressure envelopes and instantaneous contour plots are also presented. Relative merits of the two computational approaches are discussed.
Griffin, L. W.; McConnaughey, H. V.
A numerical study of the aerodynamic and thermal environment associated with axial turbine stages is presented. Computations were performed using a modification of the unsteady viscous code, ROTORI, and an improved version of the steady inviscid cascade system, MERIDL-TSONIC, coupled with boundary layer codes, BLAYER and STAN5. Two different turbine stages were analyzed: the first stage of the United Technologies Research Center Large Scale Rotating Rig (LSRR) and the first stage of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump turbine. The time-averaged airfoil midspan pressure and heat transfer profiles were predicted for numerous thermal boundary conditions including adiabatic wall, prescribed surface temperature, and prescribed heat flux. Computed solutions are compared with each other and with experimental data in the case of the LSRR calculations. Modified ROTORI predictions of unsteady pressure envelopes and instantaneous contour plots are also presented. Relative merits of the two computational approaches are discussed.
Labbe, Steven G.; Gilbert, Michael G.; Kehoe, Michael W.
The initial X-43A flight test, June 2, 2001, resulted in a mishap and loss of the vehicle. A mishap investigation board (MIB) report and findings, including the established root cause, were publicly released on July, 23, 2003. The X-43A Flight 1 Hyper-X Launch Vehicle (HXLV) failed because the vehicle control system design was deficient for the trajectory flown due to inaccurate analytical models (Pegasus heritage and HXLV specific), which overestimated the (control) system margin ? X-43A Mishap Investigation Report, Vol. I. ? included as Reference 1. Several specific errors were noted, 1) HXLV aerodynamics ? failure to model changes to wing, fin and rudder airfoil shapes due to addition of thermal protection system (TPS); 2) Fin actuation system (FAS) modeling ? under prediction of the control surface hinge moments and FAS compliance; and 3) Parametric uncertainties ? insufficient variation in the aerodynamic, FAS and control system models. In response to the MIB findings, the X-43A program has been working RTF through an approved Corrective Action Plan (CAP) over the last two years.
Draxl, C.; Churchfield, M.; Mirocha, J.; Lee, S.; Lundquist, J.; Michalakes, J.; Moriarty, P.; Purkayastha, A.; Sprague, M.; Vanderwende, B.
Wind plant aerodynamics are influenced by a combination of microscale and mesoscale phenomena. Incorporating mesoscale atmospheric forcing (e.g., diurnal cycles and frontal passages) into wind plant simulations can lead to a more accurate representation of microscale flows, aerodynamics, and wind turbine/plant performance. Our goal is to couple a numerical weather prediction model that can represent mesoscale flow [specifically the Weather Research and Forecasting model] with a microscale LES model (OpenFOAM) that can predict microscale turbulence and wake losses.
Lee, J. H.; Timmermans, J.; Su, Z.; Mancini, M.
Aerodynamic roughness height (Zom) is a key parameter required in land surface hydrological model, since errors in heat flux estimations are largely dependent on accurate optimization of this parameter. Despite its significance, it remains an uncertain parameter that is not easily determined. This is mostly because of non-linear relationship in Monin-Obukhov Similarity (MOS) and unknown vertical characteristic of vegetation. Previous studies determined aerodynamic roughness using traditional wind profile method, remotely sensed vegetation index, minimization of cost function over MOS relationship or linear regression. However, these are complicated procedures that presume high accuracy for several other related parameters embedded in MOS equations. In order to simplify a procedure and reduce the number of parameters in need, this study suggests a new approach to extract aerodynamic roughness parameter via Ensemble Kalman Filter (EnKF) that affords non-linearity and that requires only single or two heat flux measurement. So far, to our knowledge, no previous study has applied EnKF to aerodynamic roughness estimation, while a majority of data assimilation study has paid attention to land surface state variables such as soil moisture or land surface temperature. This approach was applied to grassland in semi-arid Tibetan area and maize on moderately wet condition in Italy. It was demonstrated that aerodynamic roughness parameter can inversely be tracked from data assimilated heat flux analysis. The aerodynamic roughness height estimated in this approach was consistent with eddy covariance result and literature value. Consequently, this newly estimated input adjusted the sensible heat overestimated and latent heat flux underestimated by the original Surface Energy Balance System (SEBS) model, suggesting better heat flux estimation especially during the summer Monsoon period. The advantage of this approach over other methodologies is that aerodynamic roughness height
Batina, John T.; Lee, Elizabeth M.; Kleb, William L.; Rausch, Russ D.
The current status of unstructured grid methods development in the Unsteady Aerodynamics Branch at NASA-Langley is described. These methods are being developed for unsteady aerodynamic and aeroelastic analyses. The flow solvers are highlighted which were developed for the solution of the unsteady Euler equations and selected results are given which show various features of the capability. The results demonstrate 2-D and 3-D applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.
Batina, John T.; Lee, Elizabeth M.; Kleb, William L.; Rausch, Russ D.
The current status of unstructured grid methods developed in the Unsteady Aerodynamics Branch at NASA Langley Research Center is described. These methods are being developed for unsteady aerodynamic and aeroelastic analyses. Flow solvers that have been developed for the solution of unsteady Euler equations are highlighted. The results demonstrate two and three dimensional applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.
In this article two integral computational fluid dynamics methods for steady-state and transient vehicle aerodynamic simulations are described using a Chevrolet Corvette ZR-1 surface panel model. In the last decade, road-vehicle aerodynamics have become an important design consideration. Originally, the design of low-drag shapes was given high priority due to worldwide fuel shortages that occurred in the mid-seventies. More recently, there has been increased interest in the role aerodynamics play in vehicle stability and passenger safety. Consequently, transient aerodynamics and the aerodynamics of vehicle in yaw have become important issues at the design stage. While there has been tremendous progress in Navier-Stokes methodology in the last few years, the physics of bluff-body aerodynamics are still very difficult to model correctly. Moreover, the computational effort to perform Navier-Stokes simulations from the geometric stage to complete flow solutions requires much computer time and impacts the design cycle time. In the short run, therefore, simpler methods must be used for such complicated problems. Here, two methods are described for the simulation of steady-state and transient vehicle aerodynamics.
Hartill, W. R.
A hypersonic wind tunnel test method for obtaining credible aerodynamic data on a complete hypersonic vehicle (generic X-24c) with scramjet exhaust flow simulation is described. The general problems of simulating the scramjet exhaust as well as accounting for scramjet inlet flow and vehicle forces are analyzed, and candidate test methods are described and compared. The method selected as most useful makes use of a thrust-minus-drag flow-through balance with a completely metric model. Inlet flow is diverted by a fairing. The incremental effect of the fairing is determined in the testing of two reference models. The net thrust of the scramjet module is an input to be determined in large-scale module tests with scramjet combustion. Force accounting is described, and examples of force component levels are predicted. Compatibility of the test method with candidate wind tunnel facilities is described, and a preliminary model mechanical arrangement drawing is presented. The balance design and performance requirements are described in a detailed specification. Calibration procedures, model instrumentation, and a test plan for the model are outlined.
Adenan, N. H.; Noorani, M. S. M.
The estimation of river flow is significantly related to the impact of urban hydrology, as this could provide information to solve important problems, such as flooding downstream. The nonlinear prediction method has been employed for analysis of four years of daily river flow data for the Langat River at Kajang, Malaysia, which is located in a downstream area. The nonlinear prediction method involves two steps; namely, the reconstruction of phase space and prediction. The reconstruction of phase space involves reconstruction from a single variable to the m-dimensional phase space in which the dimension m is based on optimal values from two methods: the correlation dimension method (Model I) and false nearest neighbour(s) (Model II). The selection of an appropriate method for selecting a combination of preliminary parameters, such as m, is important to provide an accurate prediction. From our investigation, we gather that via manipulation of the appropriate parameters for the reconstruction of the phase space, Model II provides better prediction results. In particular, we have used Model II together with the local linear prediction method to achieve the prediction results for the downstream area with a high correlation coefficient. In summary, the results show that Langat River in Kajang is chaotic, and, therefore, predictable using the nonlinear prediction method. Thus, the analysis and prediction of river flow in this area can provide river flow information to the proper authorities for the construction of flood control, particularly for the downstream area.
Tang, D.L.; Chang, M.F.; Sultan, M.C.
This patent describes a method of determining spark time in a spark timing system of an internal combustion engine having a plurality of cylinders and a spark period for each cylinder in which a spark occurs. It comprises: generating at least one crankshaft position reference pulse for each spark firing event, the reference pulse nearest the next spark being set to occur within a same cylinder event as the next spark; measuring at least two reference periods between recent reference pulses; calculating the spark timing synchronously with crankshaft position by performing the calculation upon receipt of the reference pulse nearest the next spark; predicting the engine speed for the next spark period from at least two reference periods including the most recent reference period; and based on the predicted speed, calculating a spark time measured from the the reference pulse nearest the next spark.
Olson, Erik D.
This paper lays out a comprehensive methodology for computing a low-speed, high-lift polar, without requiring additional details about the aircraft design beyond what is typically available at the conceptual design stage. Introducing low-order, physics-based aerodynamic analyses allows the methodology to be more applicable to unconventional aircraft concepts than traditional, fully-empirical methods. The methodology uses empirical relationships for flap lift effectiveness, chord extension, drag-coefficient increment and maximum lift coefficient of various types of flap systems as a function of flap deflection, and combines these increments with the characteristics of the unflapped airfoils. Once the aerodynamic characteristics of the flapped sections are known, a vortex-lattice analysis calculates the three-dimensional lift, drag and moment coefficients of the whole aircraft configuration. This paper details the results of two validation cases: a supercritical airfoil model with several types of flaps; and a 12-foot, full-span aircraft model with slats and double-slotted flaps.
Nguyen, Nhan T. (Inventor)
An aerodynamic control apparatus for an air vehicle improves various aerodynamic performance metrics by employing multiple spanwise flap segments that jointly form a continuous or a piecewise continuous trailing edge to minimize drag induced by lift or vortices. At least one of the multiple spanwise flap segments includes a variable camber flap subsystem having multiple chordwise flap segments that may be independently actuated. Some embodiments also employ a continuous leading edge slat system that includes multiple spanwise slat segments, each of which has one or more chordwise slat segment. A method and an apparatus for implementing active control of a wing shape are also described and include the determination of desired lift distribution to determine the improved aerodynamic deflection of the wings. Flap deflections are determined and control signals are generated to actively control the wing shape to approximate the desired deflection.
Rodný, Marek; Nolz, Reinhard; Novák, Viliam; Hlaváčiková, Hana; Loiskandl, Willibald; Himmelbauer, Margarita
The aim of this study was to present and validate an alternative evapotranspiration calculation procedure that includes specific expression for the aerodynamic resistance. Calculated daily potential evapotranspiration totals were compared to the results of FAO56 procedure application and to the results of measurements taken with a precision weighing lysimeter permanently grown with irrigated, short grass. For the examination period from March 17 through October 31, 2011, it was found that daily potential evapotranspiration estimates obtained by both calculation procedures fitted well to the lysimeter measurements. Potential evapotranspiration daily totals calculated with the use of the proposed aerodynamic resistance calculation procedure gave better results for days with higher evapotranspiration, compared to the FAO56 method. The most important is that the approach based on the proposed alternative aerodynamic resistance could be effectively used even for a wide variety of crops, because it is not limited to any particular crop.
Yates, E. Carson, Jr.
Progress in the development of computational methods for steady and unsteady aerodynamics has perennially paced advancements in aeroelastic analysis and design capabilities. Since these capabilities are of growing importance in the analysis and design of high-performance aircraft, considerable effort has been directed toward the development of appropriate aerodynamic methodology. The contributions to those efforts from the integral-equations research program at the NASA Langley Research Center is reviewed. Specifically, the current scope, progress, and plans for research and development for inviscid and viscous flows are discussed, and example applications are shown in order to highlight the generality, versatility, and attractive features of this methodology.
Dhruv, Akash V.
Flight has been one of the greatest realizations of human imagination, revolutionizing communication and transportation over the years. This has greatly influenced the growth of technology itself, enabling researchers to communicate and share their ideas more effectively, extending the human potential to create more sophisticated systems. While the end product of a sophisticated technology makes our lives easier, its development process presents an array of challenges in itself. In last decade, scientists and engineers have turned towards bio-inspiration to design more efficient and robust aerodynamic systems to enhance the ability of Unmanned Aerial Vehicles (UAVs) to be operated in cluttered environments, where tight maneuverability and controllability are necessary. Effective use of UAVs in domestic airspace will mark the beginning of a new age in communication and transportation. The design of such complex systems necessitates the need for faster and more effective tools to perform preliminary investigations in design, thereby streamlining the design process. This thesis explores the implementation of numerical panel methods for aerodynamic analysis of bio-inspired morphing wings. Numerical panel methods have been one of the earliest forms of computational methods for aerodynamic analysis to be developed. Although the early editions of this method performed only inviscid analysis, the algorithm has matured over the years as a result of contributions made by prominent aerodynamicists. The method discussed in this thesis is influenced by recent advancements in panel methods and incorporates both viscous and inviscid analysis of multi-flap wings. The surface calculation of aerodynamic coefficients makes this method less computationally expensive than traditional Computational Fluid Dynamics (CFD) solvers available, and thus is effective when both speed and accuracy are desired. The morphing wing design, which consists of sequential feather-like flaps installed
Schwing, Alan M.
Exploration Flight Test One (EFT-1) was an incredible milestone in the development NASA's Orion spacecraft. It incorporated hundreds of articles of flight test instrumentation and returned with a wealth of data. Aerodynamic surface pressures were collected during launch vehicle ascent and capsule reentry and descent. These discrete surface pressure measurements enable comparisons to computational results and ground test data. This paper details the comparisons between pre-test predictions and flight test data for the Orion MPCV Crew Module (CM) and Launch Abort Tower (LAT) during all phases of flight. Regions with strong comparisons, poor predictions, and lessons learned are discussed. 38 pressure measurements were made on the LAT during ascent. Nine of the gauges were Honeywell PPTs and the remainder were Kulite pressure transducers. In order to address bias in the Kulites, a two-point linear calibration was used and the details are discussed. Results from the flight are compared to existing database products. 44 pressure measurements were made on the CM during reentry and descent. Nine of the gauges were Honeywell PPTs and the remainder were Kulite pressure transducers. In order to address bias in the Kulites, a tare was made against the vacuum measurements as described below. Once the bias was removed from the gauges, comparisons between predicted loading and the measured results are compared.
The present volume discusses the original development of the panel method, the mapping solutions and singularity distributions of linear potential schemes, the capabilities of full-potential, Euler, and Navier-Stokes schemes, the use of the grid-generation methodology in applied aerodynamics, subsonic airfoil design, inverse airfoil design for transonic applications, the divergent trailing-edge airfoil innovation in CFD, Euler and potential computational results for selected aerodynamic configurations, and the application of CFD to wing high-lift systems. Also discussed are high-lift wing modifications for an advanced-capability EA-6B aircraft, Navier-Stokes methods for internal and integrated propulsion system flow predictions, the use of zonal techniques for analysis of rotor-stator interaction, CFD applications to complex configurations, CFD applications in component aerodynamic design of the V-22, Navier-Stokes computations of a complete F-16, CFD at supersonic/hypersonic speeds, and future CFD developments.
Erickson, Gary E.; Deloach, Richard
A collection of statistical and mathematical techniques referred to as response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration using data obtained on small-scale models at supersonic speeds in the NASA Langley Research Center Unitary Plan Wind Tunnel. The simulated Mach 3 staging was dominated by multiple shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. This motivated a partitioning of the overall inference space into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using cuboidal and spherical central composite designs capable of fitting full second-order response functions. The primary goal was to approximate the underlying overall aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle using relatively simple, lower-order polynomial functions that were piecewise-continuous across the full independent variable ranges of interest. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. The potential benefits of augmenting the central composite designs to full third order using computer-generated D-optimality criteria were also evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting low-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.
Erickson, Gary E.
Response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration at supersonic speeds in the NASA LaRC Unitary Plan Wind Tunnel. The Mach 3 staging was dominated by shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. The inference space was partitioned into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using central composite designs capable of fitting full second-order response functions. The underlying aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle were estimated using piecewise-continuous lower-order polynomial functions. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. Augmenting the central composite designs to full third-order using computer-generated D-optimality criteria was evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting lower-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.
Capon, Christopher; Boyce, Russell; Brown, Melrose
Interactions between orbiting bodies and the charged space environment are complex. The large variation in passive body parameters e.g. size, geometry and materials, makes the plasma-body interaction in Low Earth Orbit (LEO) a region rich in fundamental physical phenomena. The aerodynamic interaction of LEO orbiting bodies with the neutral environment constitutes the largest non-conservative force on the body. However in general, study of the LEO plasma-body interaction has not been concerned with external flow physics, but rather with the effects on surface charging. The impact of ionospheric flow physics on the forces on space debris (and active objects) is not well understood. The work presented here investigates the contribution that plasma-body interactions have on the flow structure and hence on the total atmospheric force vector experienced by a polar orbiting LEO body. This work applies a hybrid Particle-in-Cell (PIC) - Direct Simulation Monte Carlo (DSMC) code, pdFoam, to self-consistently model the electrostatic flowfield about a cylinder with a uniform, fixed surface potential. Flow conditions are representative of the mean conditions experienced by the Earth Observing Satellite (EOS) based on the International Reference Ionosphere model (IRI-86). The electron distribution function is represented by a non-linear Boltzmann electron fluid and ion gas-surface interactions are assumed to be that of a neutralising, conducting, thermally accommodating solid wall with diffuse reflections. The variation in flowfield and aerodynamic properties with surface potential at a fixed flow condition is investigated, and insight into the relative contributions of charged and neutral species to the flow physics experienced by a LEO orbiting body is provided. This in turn is intended to help improve the fidelity of physics-based orbit predictions for space debris and other near-Earth space objects.
Arnold, J.R.; Kelly, J.F.; Delzingaro, M.J.
Motor Operated Valve (MOV) squirrel cage AC motor rotors are susceptible to degradation under certain conditions. Premature failure can result due to high humidity/temperature environments, high running load conditions, extended periods at locked rotor conditions (i.e. > 15 seconds) or exceeding the motor`s duty cycle by frequent starts or multiple valve stroking. Exposure to high heat and moisture due to packing leaks, pressure seal ring leakage or other causes can significantly accelerate the degradation. ComEd and Liberty Technologies have worked together to provide and validate a non-intrusive method using motor power diagnostics to evaluate MOV rotor condition and predict failure. These techniques have provided a quick, low radiation dose method to evaluate inaccessible motors, identify degradation and allow scheduled replacement of motors prior to catastrophic failures.
Rhode, M. N.; Engelund, Walter C.; Mendenhall, Michael R.
Experimental longitudinal and lateral-directional aerodynamic characteristics were obtained for the Pegasus and Pegasus XL configurations over a Mach number range from 1.6 to 6 and angles of attack from -4 to +24 degrees. Angle of sideslip was varied from -6 to +6 degrees, and control surfaces were deflected to obtain elevon, aileron, and rudder effectiveness. Experimental data for the Pegasus configuration are compared with engineering code predictions performed by Nielsen Engineering & Research, Inc. (NEAR) in the aerodynamic design of the Pegasus vehicle, and with results from the Aerodynamic Preliminary Analysis System (APAS) code. Comparisons of experimental results are also made with longitudinal flight data from Flight #2 of the Pegasus vehicle. Results show that the longitudinal aerodynamic characteristics of the Pegasus and Pegasus XL configurations are similar, having the same lift-curve slope and drag levels across the Mach number range. Both configurations are longitudinally stable, with stability decreasing towards neutral levels as Mach number increases. Directional stability is negative at moderate to high angles of attack due to separated flow over the vertical tail. Dihedral effect is positive for both configurations, but is reduced 30-50 percent for the Pegasus XL configuration because of the horizontal tail anhedral. Predicted longitudinal characteristics and both longitudinal and lateral-directional control effectiveness are generally in good agreement with experiment. Due to the complex leeside flowfield, lateral-directional characteristics are not as well predicted by the engineering codes. Experiment and flight data are in good agreement across the Mach number range.
Aerodynamic investigation of an air-cooled axial-flow turbine. Part 2: Rotor blade tip-clearance effects on overall turbine performance and internal gas flow conditions: Experimental results and prediction methods
Yamamoto, A.; Takahara, K.; Nouse, H.; Mimura, F.; Inoue, S.; Usui, H.
Total turbine blade performance was investigated while changing the blade tip clearance in three ways. The internal flow at the moving blade outlet point was measured. Experimental results were compared with various theoretical methods. Increased blade clearance leads to decreased turbine efficiency.
Meakin, Robert L.
A prototype method for time-accurate simulation of multiple aerodynamic bodies in relative motion is presented. The method is general and features unsteady chimera domain decomposition techniques and an implicit approximately factored finite-difference procedure to solve the time-dependent thin-layer Navier-Stokes equations. The method is applied to a set of two- and three- dimensional test problems to establish spatial and temporal accuracy, quantify computational efficiency, and begin to test overall code robustness.
Takizawa, Kenji; Tezduyar, Tayfun E.; Boswell, Cody; Tsutsui, Yuki; Montel, Kenneth
The space-time fluid-structure interaction (STFSI) methods for 3D parachute modeling are now at a level where they can bring reliable, practical analysis to some of the most complex parachute systems, such as spacecraft parachutes. The methods include the Deforming-Spatial-Domain/Stabilized ST method as the core computational technology, and a good number of special FSI methods targeting parachutes. Evaluating the stability characteristics of a parachute based on how the aerodynamic moment varies as a function of the angle of attack is one of the practical analyses that reliable parachute FSI modeling can deliver. We describe the special FSI methods we developed for this specific purpose and present the aerodynamic-moment data obtained from FSI modeling of NASA Orion spacecraft parachutes and Japan Aerospace Exploration Agency (JAXA) subscale parachutes.
Lacasse, James M.
A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations. Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies. The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from the CFD discretized flow equations, recast in adjoint form, and solved using direct matrix inversion techniques. This investigation is performed to demonstrate an aerodynamic shape optimization technique on a multiblock domain and its applicability to complex geometries. The objectives are accomplished by shape optimizing two aerodynamic configurations. First, the shape optimization of a transonic airfoil is performed to investigate the behavior of the method in highly nonlinear flows and the effect of different grid blocking strategies on the procedure. Secondly, shape optimization of a two-element configuration in subsonic flow is completed. Cases are presented for this configuration to demonstrate the effect of simultaneously reshaping interfering elements. The aerodynamic shape optimization is shown to produce supercritical type airfoils in the transonic flow from an initially symmetric airfoil. Multiblocking effects the path of optimization while providing similar results at the conclusion. Simultaneous reshaping of elements is shown to be more effective than individual element reshaping due to the inclusion of mutual interference effects.
Pao, S. Paul
The role of acoustics and unsteady aerodynamics research in understanding the fundamental physics of time-dependent fluid phenomena is reviewed. The key issues are illustrated by considering the sound radiation of turbulent jets and the aeroacoustics of rotating bodies such as helicopter rotors. The importance of computational methods as a link between aerodynamics and acoustics is also discussed. It is noted that where acoustic analogy techniques are sufficiently accurate, unsteady aerodynamics can be used for acoustic prediction. In supersonic problems where acoustics and aerodynamics are coupled, an integrated nonlinear analysis can provide an accurate problem solution.
Milholen, William E., II
Computational Fluid Dynamics based design methods are maturing to the point that they are beginning to be used in the aircraft design process. Many design methods however have demonstrated deficiencies in the leading edge region of airfoil sections. The objective of the present research is to develop an efficient inverse design method which is valid in the leading edge region. The new design method is a streamline curvature method, and a new technique is presented for modeling the variation of the streamline curvature normal to the surface. The new design method allows the surface coordinates to move normal to the surface, and has been incorporated into the Constrained Direct Iterative Surface Curvature (CDISC) design method. The accuracy and efficiency of the design method is demonstrated using both two-dimensional and three-dimensional design cases.
Yates, E. Carson, Jr.; Whitlow, Woodrow, Jr.
The current scope, recent progress, and plans for research and development of computational methods for unsteady aerodynamics at the NASA Langley Research Center are reviewed. Both integral equations and finite difference methods for inviscid and viscous flows are discussed. Although the great bulk of the effort has focused on finite difference solution of the transonic small perturbation equation, the integral equation program is given primary emphasis here because it is less well known.
Yates, E. Carson, Jr.; Whitlow, Woodrow, Jr.
The current scope, recent progress, and plans for research and development of computational methods for unsteady aerodynamics at the NASA Langley Research Center are reviewed. Both integral-equations and finite-difference method for inviscid and viscous flows are discussed. Although the great bulk of the effort has focused on finite-difference solution of the transonic small-perturbation equation, the integral-equation program is given primary emphasis here because it is less well known.
Bennett, R. M.
The method of integral relations is applied in a one-strip approximation to the perturbation equations governing small motions of an inclined, sharp-edged, flat surface about the mean supersonic steady flow. Algebraic expressions for low reduced-frequency aerodynamics are obtained and a set of ordinary differential equations are obtained for general oscillatory motion. Results are presented for low reduced-frequency aerodynamics and for the variation of the unsteady forces with frequency. The method gives accurate results for the aerodynamic forces at low reduced frequency which are in good agreement with available experimental data. However, for cases in which the aerodynamic forces vary rapidly with frequency, the results are qualitatively correct, but of limited accuracy. Calculations indicate that for a range of inclination angles near shock detachment such that the flow in the shock layer is low supersonic, the aerodynamic forces vary rapidly both with inclination angle and with reduced frequency.
Edwards, John W.; Malone, John B.
The current status of computational methods for unsteady aerodynamics and aeroelasticity is reviewed. The key features of challenging aeroelastic applications are discussed in terms of the flowfield state: low-angle high speed flows and high-angle vortex-dominated flows. The critical role played by viscous effects in determining aeroelastic stability for conditions of incipient flow separation is stressed. The need for a variety of flow modeling tools, from linear formulations to implementations of the Navier-Stokes equations, is emphasized. Estimates of computer run times for flutter calculations using several computational methods are given. Applications of these methods for unsteady aerodynamic and transonic flutter calculations for airfoils, wings, and configurations are summarized. Finally, recommendations are made concerning future research directions.
Lyle, Karen H.
Acceptance of new spacecraft structural architectures and concepts requires validated design methods to minimize the expense involved with technology validation via flighttesting. This paper explores the implementation of probabilistic methods in the sensitivity analysis of the structural response of a Hypersonic Inflatable Aerodynamic Decelerator (HIAD). HIAD architectures are attractive for spacecraft deceleration because they are lightweight, store compactly, and utilize the atmosphere to decelerate a spacecraft during re-entry. However, designers are hesitant to include these inflatable approaches for large payloads or spacecraft because of the lack of flight validation. In the example presented here, the structural parameters of an existing HIAD model have been varied to illustrate the design approach utilizing uncertainty-based methods. Surrogate models have been used to reduce computational expense several orders of magnitude. The suitability of the design is based on assessing variation in the resulting cone angle. The acceptable cone angle variation would rely on the aerodynamic requirements.
Edwards, John W.; Malone, John B.
The status of computational methods for unsteady aerodynamics and aeroelasticity is reviewed. The key features of challenging aeroelastic applications is discussed in terms of the flowfield state - low angle high speed flows and high angle vortex dominated flows. The critical role played by viscous effects in determining aeroelastic stability for conditions of incipient flow separation is stressed. The need for a variety of flow modeling tools, from linear formulations to implementations of the Navier-Stokes equations, is emphasized. Estimates of computer run times for flutter calculations using several computational methods are given. Applications of these methods for unsteady aerodynamic and transonic flutter calculations for airfoils, wings, and configurations are summarized. Finally, recommendations are made concerning future research directions.
Edwards, J. W.; Malone, J. B.
The current status of computational methods for unsteady aerodynamics and aeroelasticity is reviewed. The key features of challenging aeroelastic applications are discussed in terms of the flowfield state: low-angle high speed flows and high-angle vortex-dominated flows. The critical role played by viscous effects in determining aeroelastic stability for conditions of incipient flow separation is stressed. The need for a variety of flow modeling tools, from linear formulations to implementations of the Navier-Stokes equations, is emphasized. Estimates of computer run times for flutter calculations using several computational methods are given. Application of these methods for unsteady aerodynamic and transonic flutter calculations for airfoils, wings, and configurations are summarized. Finally, recommendations are made concerning future research directions.
Lai, Y. G. J.; Przekwas, A. J.; So, R. M. C.
In the past, most aerodynamic flow calculations were carried out with density-based numerical methods and zero-equation turbulence models. However, pressure-based methods and more advanced turbulence models have been routinely used in industry for many internal flow simulations and for incompressible flows. Unfortunately, their usefulness in calculating aerodynamic flows is still not well demonstrated and accepted. In this study, an advanced pressure-based numerical method and a recently proposed near-wall compressible two-equation turbulence model are used to calculate external aerodynamic flows. Several TVD-type schemes are extended to pressure-based method to better capture discontinuities such as shocks. Some improvements are proposed to accelerate the convergence of the numerical method. A compressible near-wall two-equation turbulence model is then implemented to calculate transonic turbulent flows over NACA 0012 and RAE 2822 airfoils with and without shocks. The calculated results are compared with wind tunnel data as well as with results obtained from the Baldwin-Lomax model. The performance of the two-equation turbulence model is evaluated and its merits or lack thereof are discussed.
Anderson, L. R.
The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.
Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian method in which the domain is represented by particles. Each particle is assigned properties such as mass, pressure, density, temperature, and velocity. These properties are then evaluated at the particle positions using a smoothing kernel that integrates over the values of the surrounding particles. In the present study the SPH method is first used to obtain numerical solutions for fluid flows over a cylinder and then we are going to apply the same principle over an airfoil obstacle.
Chen, Zuoyi; Guo, Jingrong
Progress in the use of the Hodograph method of aerodynamic design is discussed. It was found that there are some restricted conditions in the application of Hodograph design to transonic turbine and compressor cascades. The Hodograph method is suitable not only to the transonic turbine cascade but also to the transonic compressor cascade. The three dimensional Hodograph method will be developed after obtaining the basic equation for the three dimensional Hodograph method. As an example of the Hodograph method, the use of the method to design a transonic turbine and compressor cascade is discussed.
Abhari, R.S.; Giles, M.
An unsteady, compressible, two-dimensional, thin shear layer Navier-Stokes solver is modified to predict the motion-dependent unsteady flow around oscillating airfoils in a cascade. A quasi-three-dimensional formulation is used to account for the stream-wise variation of streamtube height. The code uses Ni`s Lax-Wendroff algorithm in the outer region an implicit ADI method in the inner region, conservative coupling at the interface, and the Baldwin-Lomax turbulence model. The computational mesh consists of an O-grid around each blade plus an unstructured outer grid of quadrilateral or triangular cells. The unstructured computational grid was adapted to the flow to better resolve shocks and wakes. Motion of each airfoil was simulated at each time step by stretching and compressing the mesh within the O-grid. This imposed motion consists of harmonic solid body translation in two directions and rotation, combined with the correct interblade phase angles. The validity of the code is illustrated by comparing its predictions to a number of test cases, including an axially oscillating flat plate in laminar flow, the Aeroelasticity of Turbomachines Symposium Fourth Standard Configuration (a transonic turbine cascade), and the Seventh Standard Configuration (a transonic compressor cascade). The overall comparison between the predictions and the test data is reasonably good. A numerical study on a generic transonic compressor rotor was performed in which the impact of varying the amplitude of the airfoil oscillation on the normalized predicted magnitude and phase of the unsteady pressure around the airfoil was studied. It was observed that for this transonic compressor, the nondimensional aerodynamic damping was influenced by the amplitude of the oscillation.
Carlson, Leland A.
The primary accomplishments of the project are as follows: (1) Using the transonic small perturbation equation as a flowfield model, the project demonstrated that the quasi-analytical method could be used to obtain aerodynamic sensitivity coefficients for airfoils at subsonic, transonic, and supersonic conditions for design variables such as Mach number, airfoil thickness, maximum camber, angle of attack, and location of maximum camber. It was established that the quasi-analytical approach was an accurate method for obtaining aerodynamic sensitivity derivatives for airfoils at transonic conditions and usually more efficient than the finite difference approach. (2) The usage of symbolic manipulation software to determine the appropriate expressions and computer coding associated with the quasi-analytical method for sensitivity derivatives was investigated. Using the three dimensional fully conservative full potential flowfield model, it was determined that symbolic manipulation along with a chain rule approach was extremely useful in developing a combined flowfield and quasi-analytical sensitivity derivative code capable of considering a large number of realistic design variables. (3) Using the three dimensional fully conservative full potential flowfield model, the quasi-analytical method was applied to swept wings (i.e. three dimensional) at transonic flow conditions. (4) The incremental iterative technique has been applied to the three dimensional transonic nonlinear small perturbation flowfield formulation, an equivalent plate deflection model, and the associated aerodynamic and structural discipline sensitivity equations; and coupled aeroelastic results for an aspect ratio three wing in transonic flow have been obtained.
Batina, John T.
The current status of the development of unstructured grid methods in the Unsteady Aerodynamics Branch at NASA-Langley is described. These methods are being developed for steady and unsteady aerodynamic applications. The flow solvers that were developed for the solution of the unsteady Euler and Navier-Stokes equations are highlighted and selected results are given which demonstrate various features of the capability. The results demonstrate 2-D and 3-D applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.
Batina, John T.
The current status of the development of unstructured grid methods in the Unsteady Aerodynamic Branch at NASA-Langley is described. These methods are being developed for steady and unsteady aerodynamic applications. The flow solvers that were developed for the solution of the unsteady Euler and Navier-Stokes equations are highlighted and selected results are given which demonstrate various features of the capability. The results demonstrate 2-D and 3-D applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.
Raj, P.; Long, L. N.
The current capabilities and the future plans for a three dimensional Euler Aerodynamic Method are described. The basic solution algorithm is based on the finite volume, Runge-Kutta pseudo-time-stepping scheme of FLO-57. Several modifications to improve accuracy and computational efficiency were incorporated and others are being investigated. The computer code is used to analyze a cropped delta wing at 0.6 Mach number and an arrow wing at 0.85 Mach number. Computed aerodynamic parameters are compared with experimental data. In all cases, the configuration is impulsively started and no Kutta condition is applied at sharp edges. The results indicate that with additional development and validation, the present method will be a useful tool for engineering analysis of high speed aircraft.
Jorgensen, L. H.
An engineering-type method is presented for computing normal-force and pitching-moment coefficients for slender bodies of circular and noncircular cross section alone and with lifting surfaces. In this method, a semi-empirical term representing viscous-separation crossflow is added to a term representing potential-theory crossflow. For many bodies of revolution, computed aerodynamic characteristics are shown to agree with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. The angles of attack extend from 0 deg to 180 deg for M = 2.9 from 0 deg to 60 deg for M = 0.6 to 2.0. For several bodies of elliptic cross section, measured results are also predicted reasonably well over the investigated Mach number range from 0.6 to 2.0 and at angles of attack from 0 deg to 60 deg. As for the bodies of revolution, the predictions are best for supersonic Mach numbers. For body-wing and body-wing-tail configurations with wings of aspect ratios 3 and 4, measured normal-force coefficients and centers are predicted reasonably well at the upper test Mach number of 2.0. Vapor-screen and oil-flow pictures are shown for many body, body-wing and body-wing-tail configurations. When spearation and vortex patterns are asymmetric, undesirable side forces are measured for the models even at zero sideslip angle. Generally, the side-force coefficients decrease or vanish with the following: increase in Mach number, decrease in nose fineness ratio, change from sharp to blunt nose, and flattening of body cross section (particularly the body nose).
Jorgensen, L. H.
An engineering-type method is presented for computing normal-force and pitching-moment coefficients for slender bodies of circular and noncircular cross section alone and with lifting surfaces. In this method, a semi-empirical term representing viscous-separation crossflow is added to a term representing potential-theory crossflow. For many bodies of revolution, computed aerodynamic characteristics are shown to agree with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. For several bodies of elliptic cross section, measured results are also predicted reasonably well over the investigated Mach number range from 0.6 to 2.0 and at angles of attack from 0 to 60 deg. As for the bodies of revolution, the predictions are best for supersonic Mach numbers. For body-wing and body-wing-tail configurations with wings of aspect ratios 3 and 4, measured normal-force coefficients and centers are predicted reasonably well at the upper test Mach number of 2.0. However, with a decrease in Mach number to 0.6, the agreement for C sub N rapidly deteriorates, although the normal-force centers remain in close agreement. Vapor-screen and oil-flow pictures are shown for many body, body-wing, and body-wing-tail configurations. When separation and vortex patterns are asymmetric, undesirable side forces are measured for the models even at zero sideslip angle. Generally, the side-force coefficients decrease or vanish with the following: increase in Mach number, decrease in nose fineness ratio, change from sharp to blunt nose, and flattening of body cross section (particularly the body nose).
Loomis, Mark P.; Venkatapathy, Ethiraj; Papadopoulos, Periklis; Davies, Carol B.; Berry, Scott; Horvath, Tom; Campbell, Charles
This paper presents the results of a computational fluid dynamics validation exercise using the General Aerodynamic Simulation Program (GASP) code for vehicles under consideration for the X-38 program. The scope of the exercise involved simulation and comparison with wind tunnel aeroheating and aerodynamic data in the hypersonic regime. In addition, simulations were performed and results are presented for the full-scale flight vehicle along an expected re- entry trajectory.
Streett, C. L.
Fundamental concepts underlying spectral collocation methods, especially pertaining to their use in the solution of partial differential equations, are outlined. Theoretical accuracy results are reviewed and compared with results from test problems. A number of practical aspects of the construction and use of spectral methods are detailed, along with several solution schemes which have found utility in applications of spectral methods to practical problems. Results from a few of the successful applications of spectral methods to problems of aerodynamic and fluid mechanic interest are then outlined, followed by a discussion of the problem areas in spectral methods and the current research under way to overcome these difficulties.
Harrison, B. A.; Richard, M.
The information necessary for execution of the digital computer program L216 on the CDC 6600 is described. L216 characteristics are based on the doublet lattice method. Arbitrary aerodynamic configurations may be represented with combinations of nonplanar lifting surfaces composed of finite constant pressure panel elements, and axially summetric slender bodies composed of constant pressure line elements. Program input consists of configuration geometry, aerodynamic parameters, and modal data; output includes element geometry, pressure difference distributions, integrated aerodynamic coefficients, stability derivatives, generalized aerodynamic forces, and aerodynamic influence coefficient matrices. Optionally, modal data may be input on magnetic field (tape or disk), and certain geometric and aerodynamic output may be saved for subsequent use.
Axelson, J. A.
The AEROX program estimates lift, induced-drag and pitching moments to high angles (typ. 60 deg) for wings and for wingbody combinations with or without an aft horizontal tail. Minimum drag coefficients are not estimated, but may be input for inclusion in the total aerodynamic parameters which are output in listed and plotted formats. The theory, users' guide, test cases, and program listing are presented.
Various control analysis, design, and simulation techniques of aeroservoelastic systems require the equations of motion to be cast in a linear, time-invariant state-space form. In order to account for unsteady aerodynamics, rational function approximations must be obtained to represent them in the first order equations of the state-space formulation. A computer program, MIST, has been developed which determines minimum-state approximations of the coefficient matrices of the unsteady aerodynamic forces. The Minimum-State Method facilitates the design of lower-order control systems, analysis of control system performance, and near real-time simulation of aeroservoelastic phenomena such as the outboard-wing acceleration response to gust velocity. Engineers using this program will be able to calculate minimum-state rational approximations of the generalized unsteady aerodynamic forces. Using the Minimum-State formulation of the state-space equations, they will be able to obtain state-space models with good open-loop characteristics while reducing the number of aerodynamic equations by an order of magnitude more than traditional approaches. These low-order state-space mathematical models are good for design and simulation of aeroservoelastic systems. The computer program, MIST, accepts tabular values of the generalized aerodynamic forces over a set of reduced frequencies. It then determines approximations to these tabular data in the LaPlace domain using rational functions. MIST provides the capability to select the denominator coefficients in the rational approximations, to selectably constrain the approximations without increasing the problem size, and to determine and emphasize critical frequency ranges in determining the approximations. MIST has been written to allow two types data weighting options. The first weighting is a traditional normalization of the aerodynamic data to the maximum unit value of each aerodynamic coefficient. The second allows weighting the
Cheung, Samson H.; Edwards, Thomas A.
This paper presents a new methodology for the optimization of supersonic airplane designs to meet the dual design objectives of low sonic boom and high aerodynamic performance. Two sets of design parameters are used on an existing High Speed Civil Transport (HSCT) configuration to maximize the aerodynamic performance and minimize the sonic boom under the flight track. One set of the parameters perturbs the camber line of the wing sections to maximize the lift-over-drag ratio (L/D). A preliminary optimization run yielded a 3.75 percent improvement in L/D over a baseline low-boom configuration. The other set of parameters modifies the fuselage area to achieve a target F-function. Starting from an initial configuration with strong bow, wing, and tail shocks, a modified design with a flat-top signature is obtained. The methods presented can easily incorporate other design variables and objective functions. Extensions to the present capability in progress are described.
Polanský, Jiří; Kalmár, László; Gášpár, Roman
The main aim of this paper is determine the centrifugal fan with forward curved blades aerodynamic characteristics based on numerical modeling. Three variants of geometry were investigated. The first, basic "A" variant contains 12 blades. The geometry of second "B" variant contains 12 blades and 12 semi-blades with optimal length . The third, control variant "C" contains 24 blades without semi-blades. Numerical calculations were performed by CFD Ansys. Another aim of this paper is to compare results of the numerical simulation with results of approximate numerical procedure. Applied approximate numerical procedure  is designated to determine characteristics of the turbulent flow in the bladed space of a centrifugal-flow fan impeller. This numerical method is an extension of the hydro-dynamical cascade theory for incompressible and inviscid fluid flow. Paper also partially compares results from the numerical simulation and results from the experimental investigation. Acoustic phenomena observed during experiment, during numerical simulation manifested as deterioration of the calculation stability, residuals oscillation and thus also as a flow field oscillation. Pressure pulsations are evaluated by using frequency analysis for each variant and working condition.
Reuther, James; Alonso, Juan Jose; Vassberg, John C.; Jameson, Antony; Martinelli, Luigi
The work presented in this paper describes the application of a multiblock gridding strategy to the solution of aerodynamic design optimization problems involving complex configurations. The design process is parallelized using the MPI (Message Passing Interface) Standard such that it can be efficiently run on a variety of distributed memory systems ranging from traditional parallel computers to networks of workstations. Substantial improvements to the parallel performance of the baseline method are presented, with particular attention to their impact on the scalability of the program as a function of the mesh size. Drag minimization calculations at a fixed coefficient of lift are presented for a business jet configuration that includes the wing, body, pylon, aft-mounted nacelle, and vertical and horizontal tails. An aerodynamic design optimization is performed with both the Euler and Reynolds Averaged Navier-Stokes (RANS) equations governing the flow solution and the results are compared. These sample calculations establish the feasibility of efficient aerodynamic optimization of complete aircraft configurations using the RANS equations as the flow model. There still exists, however, the need for detailed studies of the importance of a true viscous adjoint method which holds the promise of tackling the minimization of not only the wave and induced components of drag, but also the viscous drag.
Soohoo, P.; Morino, L.; Noll, R. B.; Ham, N. D.
The Green's function method was used to study tilting proprotor aircraft aerodynamics with particular application to the problem of the mutual interference of the wing-fuselage-tail-rotor wake configuration. While the formulation is valid for fully unsteady rotor aerodynamics, attention was directed to steady state aerodynamics, which was achieved by replacing the rotor with the actuator disk approximation. The use of an actuator disk analysis introduced a mathematical singularity into the formulation; this problem was studied and resolved. The pressure distribution, lift, and pitching moment were obtained for an XV-15 wing-fuselage-tail rotor configuration at various flight conditions. For the flight configurations explored, the effects of the rotor wake interference on the XV-15 tilt rotor aircraft yielded a reduction in the total lift and an increase in the nose-down pitching moment. This method provides an analytical capability that is simple to apply and can be used to investigate fuselage-tail rotor wake interference as well as to explore other rotor design problem areas.
This conference was held December 4--8, 1994 in Asilomar, California. The purpose of this meeting was to provide a forum for exchange of state-of-the-art information concerning the prediction of protein structure. Attention if focused on the following: comparative modeling; sequence to fold assignment; and ab initio folding.
Campbell, John P; Mckinney, Marion O
A summary of methods for making dynamic lateral stability and response calculations and for estimating the aerodynamic stability derivatives required for use in these calculations is presented. The processes of performing calculations of the time histories of lateral motions, of the period and damping of these motions, and of the lateral stability boundaries are presented as a series of simple straightforward steps. Existing methods for estimating the stability derivatives are summarized and, in some cases, simple new empirical formulas are presented. Detailed estimation methods are presented for low-subsonic-speed conditions but only a brief discussion and a list of references are given for transonic and supersonic speed conditions.
Pilon, Anthony R.; Lyrintzis, Anastasios S.
The accurate prediction of aerodynamically generated noise has become an important goal over the past decade. Aeroacoustics must now be an integral part of the aircraft design process. The direct calculation of aerodynamically generated noise with CFD-like algorithms is plausible. However, large computer time and memory requirements often make these predictions impractical. It is therefore necessary to separate the aeroacoustics problem into two parts, one in which aerodynamic sound sources are determined, and another in which the propagating sound is calculated. This idea is applied in acoustic analogy methods. However, in the acoustic analogy, the determination of far-field sound requires the solution of a volume integral. This volume integration again leads to impractical computer requirements. An alternative to the volume integrations can be found in the Kirchhoff method. In this method, Green's theorem for the linear wave equation is used to determine sound propagation based on quantities on a surface surrounding the source region. The change from volume to surface integrals represents a tremendous savings in the computer resources required for an accurate prediction. This work is concerned with the development of enhancements of the Kirchhoff method for use in a wide variety of aeroacoustics problems. This enhanced method, the modified Kirchhoff method, is shown to be a Green's function solution of Lighthill's equation. It is also shown rigorously to be identical to the methods of Ffowcs Williams and Hawkings. This allows for development of versatile computer codes which can easily alternate between the different Kirchhoff and Ffowcs Williams-Hawkings formulations, using the most appropriate method for the problem at hand. The modified Kirchhoff method is developed primarily for use in jet aeroacoustics predictions. Applications of the method are shown for two dimensional and three dimensional jet flows. Additionally, the enhancements are generalized so that
Holmes, Bruce J.; Schairer, Edward; Hicks, Gary; Wander, Stephen; Blankson, Isiaiah; Rose, Raymond; Olson, Lawrence; Unger, George
Presented here is a comprehensive review of the following aerodynamics elements: computational methods and applications, computational fluid dynamics (CFD) validation, transition and turbulence physics, numerical aerodynamic simulation, drag reduction, test techniques and instrumentation, configuration aerodynamics, aeroacoustics, aerothermodynamics, hypersonics, subsonic transport/commuter aviation, fighter/attack aircraft and rotorcraft.
Williams, Louis J.; Hessenius, Kristin A.; Corsiglia, Victor R.; Hicks, Gary; Richardson, Pamela F.; Unger, George; Neumann, Benjamin; Moss, Jim
The annual accomplishments is reviewed for the Aerodynamics Division during FY 1991. The program includes both fundamental and applied research directed at the full spectrum of aerospace vehicles, from rotorcraft to planetary entry probes. A comprehensive review is presented of the following aerodynamics elements: computational methods and applications; CFD validation; transition and turbulence physics; numerical aerodynamic simulation; test techniques and instrumentation; configuration aerodynamics; aeroacoustics; aerothermodynamics; hypersonics; subsonics; fighter/attack aircraft and rotorcraft.
Agarwal, Ramesh K.; Desse, Jerry E.
The current state of the art in computational aerodynamics is described. Recent advances in the discretization of surface geometry, grid generation, and flow simulation algorithms have led to flowfield predictions for increasingly complex and realistic configurations. As a result, computational aerodynamics is emerging as a crucial enabling technology for the development and design of flight vehicles. Examples illustrating the current capability for the prediction of aircraft, launch vehicle and helicopter flowfields are presented. Unfortunately, accurate modeling of turbulence remains a major difficulty in the analysis of viscosity-dominated flows. In the future inverse design methods, multidisciplinary design optimization methods, artificial intelligence technology and massively parallel computer technology will be incorporated into computational aerodynamics, opening up greater opportunities for improved product design at substantially reduced costs.
Srinivasan, G. R.; Baeder, J. D.
This paper outlines some recent advances in the application of the Euler and Navier-Stokes computational fluid dynamics methods to analyze nonlinear problems of helicopter aerodynamics and acoustics. A complete flowfield simulation of helicopters is currently not feasible with these methods. However, the use of the state-of-the-art numerical algorithms in conjunction with powerful supercomputers, like the Cray-2, have enabled notable progress to be made in modeling several individual components of this complex flow in hover and forward flight.
Wood, R. M.; Miller, D. S.; Brentner, K. S.
A theoretical and experimental investigation has been conducted to evaluate the fundamental supersonic aerodynamic characteristics of a generic twin-body model at a Mach number of 2.70. Results show that existing aerodynamic prediction methods are adequate for making preliminary aerodynamic estimates.
Carney, Paul R.; Myers, Stephen; Geyer, James D.
Epilepsy, one of the most common neurological diseases, affects over 50 million people worldwide. Epilepsy can have a broad spectrum of debilitating medical and social consequences. Although antiepileptic drugs have helped treat millions of patients, roughly a third of all patients have seizures that are refractory to pharmacological intervention. The evolution of our understanding of this dynamic disease leads to new treatment possibilities. There is great interest in the development of devices that incorporate algorithms capable of detecting early onset of seizures or even predicting them hours before they occur. The lead time provided by these new technologies will allow for new types of interventional treatment. In the near future, seizures may be detected and aborted before physical manifestations begin. In this chapter we discuss the algorithms that make these devices possible and how they have been implemented to date. We also compare and contrast these measures, and review their individual strengths and weaknesses. Finally, we illustrate how these techniques can be combined in a closed-loop seizure prevention system. PMID:22078526
Race car performance depends on elements such as the engine, tires, suspension, road, aerodynamics, and of course the driver. In recent years, however, vehicle aerodynamics gained increased attention, mainly due to the utilization of the negative lift (downforce) principle, yielding several important performance improvements. This review briefly explains the significance of the aerodynamic downforce and how it improves race car performance. After this short introduction various methods to generate downforce such as inverted wings, diffusers, and vortex generators are discussed. Due to the complex geometry of these vehicles, the aerodynamic interaction between the various body components is significant, resulting in vortex flows and lifting surface shapes unlike traditional airplane wings. Typical design tools such as wind tunnel testing, computational fluid dynamics, and track testing, and their relevance to race car development, are discussed as well. In spite of the tremendous progress of these design tools (due to better instrumentation, communication, and computational power), the fluid dynamic phenomenon is still highly nonlinear, and predicting the effect of a particular modification is not always trouble free. Several examples covering a wide range of vehicle shapes (e.g., from stock cars to open-wheel race cars) are presented to demonstrate this nonlinear nature of the flow field.
Bristow, D. R.
Information on sensitivity analysis in computational aerodynamics is given in outline, graphical, and chart form. The prediction accuracy if the MCAERO program, a perturbation analysis method, is discussed. A procedure for calculating perturbation matrix, baseline wing paneling for perturbation analysis test cases and applications of an inviscid sensitivity matrix are among the topics covered.
The aerodynamic energy method is used in this paper to synthesize control laws for NASA's Drone for Aerodynamic and Structural Testing-Aerodynamic Research Wing 1 (DAST-ARW1) mathematical model. The performance of these control laws in terms of closed-loop flutter dynamic pressure, control surface activity, and robustness is compared against other control laws that appear in the literature and relate to the same model. A control law synthesis technique that makes use of the return difference singular values is developed in this paper. it is based on the aerodynamic energy approach and is shown to yield results superior to those given in the literature and based on optimal control theory. Nyquist plots are presented together with a short discussion regarding the relative merits of the minimum singular value as a measure of robustness, compared with the more traditional measure of robustness involving phase and gain margins.
The aerodynamic energy method is used to synthesize control laws for NASA's drone for aerodynamic and structural testing-aerodynamic research wing 1 (DAST-ARW1) mathematical model. The performance of these control laws in terms of closed-loop flutter dynamic pressure, control surface activity, and robustness is compared with other control laws that relate to the same model. A control law synthesis technique that makes use of the return difference singular values is developed. It is based on the aerodynamic energy approach and is shown to yield results that are superior to those results given in the literature and are based on optimal control theory. Nyquist plots are presented, together with a short discussion regarding the relative merits of the minimum singular value as a measure of robustness as compared with the more traditional measure involving phase and gain margins.
Young, William R.
Small resonant targets used in conjunction with a microwave reflectometer to determine the deformation of the Hypersonic Inflatable Aerodynamic Decelerator (HIAD) during reentry are investigated. The reflectometer measures the distance to the targets and from this the HIAD deformation is determined. The HIAD is used by the Inflatable Reentry Vehicle Experiment (IRVE) which investigates the use of inflatable heat shields for atmospheric reentry. After several different microwave reflectometer systems were analyzed and compared it was determined that the most desirable for this application is the Frequency Doubling Target method.
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
A modular framework for aerodynamic optimization of complex geometries is developed. By working directly with a parametric CAD system, complex-geometry models are modified nnd tessellated in an automatic fashion. The use of a component-based Cartesian method significantly reduces the demands on the CAD system, and also provides for robust and efficient flowfield analysis. The optimization is controlled using either a genetic or quasi-Newton algorithm. Parallel efficiency of the framework is maintained even when subject to limited CAD resources by dynamically re-allocating the processors of the flow solver. Overall, the resulting framework can explore designs incorporating large shape modifications and changes in topology.
Tiffany, Sherwood H.; Adams, William M., Jr.
The approximation of unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft are discussed. Two methods of formulating these approximations are extended to include the same flexibility in constraining the approximations and the same methodology in optimizing nonlinear parameters as another currently used extended least-squares method. Optimal selection of nonlinear parameters is made in each of the three methods by use of the same nonlinear, nongradient optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is lower order than that required when no optimization of the nonlinear terms is performed. The free linear parameters are determined using the least-squares matrix techniques of a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from different approaches are described and results are presented that show comparative evaluations from application of each of the extended methods to a numerical example.
Glassman, A. J.; Moffitt, T. P.
Cursory review of some recent work that has been done in turbine aerodynamic research. Topics discussed include the aerodynamic effect of turbine coolant, high work-factor (ratio of stage work to square of blade speed) turbines, and computer methods for turbine design and performance prediction. Experimental cooled-turbine aerodynamics programs using two-dimensional cascades, full annular cascades, and cold rotating turbine stage tests are discussed with some typical results presented. Analytically predicted results for cooled blade performance are compared to experimental results. The problems and some of the current programs associated with the use of very high work factors for fan-drive turbines of high-bypass-ratio engines are discussed. Computer programs have been developed for turbine design-point performance, off-design performance, supersonic blade profile design, and the calculation of channel velocities for subsonic and transonic flowfields. The use of these programs for the design and analysis of axial and radial turbines is discussed.
Lamar, John E.; Cronin, Catherine K.; Scott, Laura E.
A review of the steps taken to establish an international virtual laboratory (VL) at the NASA Langley Research Center for aerodynamic prediction and comparison of flight data in the post-09/11/2001 cyber-terrorist environment is detailed here. The key features of the VL include an intuitive, web-based user interface for ease of access, a secure high-speed Internet connection between browser and server, a relational database architecture for data and information search, and a secure file-storage system. The detailed planning and handling of such issues as security, computer firewall access and legal protection of data are provided.
Hess, J. R.; Bear, R. L.
A viable, single engine, supersonic V/STOL fighter/attack aircraft concept was defined. This vectored thrust, canard wing configuration utilizes an advanced technology separated flow engine with fan stream burning. The aerodynamic characteristics of this configuration were estimated and performance evaluated. Significant aerodynamic and aerodynamic propulsion interaction uncertainties requiring additional investigation were identified. A wind tunnel model concept and test program to resolve these uncertainties and validate the aerodynamic prediction methods were defined.
Adamov, N. P.; Kharitonov, A. M.; Chasovnikov, E. A.
The paper deals with an experimental study of stationary and nonstationary aerodynamic characteristics of a circular cone in the range of Mach numbers 1.75-7. The experimental equipment and the method of determining the aerodynamic characteristics are briefly described. The integral aerodynamic characteristics of the model in tests with force measurements and the aerodynamic derivatives of the pitching moment in dynamic tests on a setup with free oscillations are obtained. The experimental data are compared with numerical predictions.
Using the two-dimensional ARC2D Navier-Stokes flow solver analyses were conducted to predict the sectional aerodynamic characteristics of the flapped NACA-0015 airfoil section. To facilitate the analyses and the generation of the computational grids, the airfoil with the deflected trailing edge flap was treated as a single element airfoil with no allowance for a gap between the flap's leading edge and the base of the forward portion of the airfoil. Generation of the O-type computational grids was accomplished using the HYGRID hyperbolic grid generation program. Results were obtained for a wide range of Mach numbers, angles of attack and flap deflections. The predicted sectional lift, drag and pitching moment values for the airfoil were then cast in tabular format (C81) to be used in lifting-line helicopter rotor aerodynamic performance calculations. Similar were also generated for the flap. Mathematical expressions providing the variation of the sectional lift and pitching moment coefficients for the airfoil and for the flap as a function of flap chord length and flap deflection angle were derived within the context of thin airfoil theory. The airfoil's sectional drag coefficient were derived using the ARC2D drag predictions for equivalent two dimensional flow conditions.
Grauer, Jared A.; Morelli, Eugene A.
A new formulation of the filter-error method for estimating aerodynamic parameters in nonlinear aircraft dynamic models during turbulence was developed and demonstrated. The approach uses an estimate of the measurement noise covariance to identify the model parameters, their uncertainties, and the process noise covariance, in a relaxation method analogous to the output-error method. Prior information on the model parameters and uncertainties can be supplied, and a post-estimation correction to the uncertainty was included to account for colored residuals not considered in the theory. No tuning parameters, needing adjustment by the analyst, are used in the estimation. The method was demonstrated in simulation using the NASA Generic Transport Model, then applied to the subscale T-2 jet-engine transport aircraft flight. Modeling results in different levels of turbulence were compared with results from time-domain output error and frequency- domain equation error methods to demonstrate the effectiveness of the approach.
Seginer, A.; Rusak, Z.; Wasserstrom, E.
Nonlinear panel methods have no proof for the existence and uniqueness of their solutions. The convergence characteristics of an iterative, nonlinear vortex-lattice method are, therefore, carefully investigated. The effects of several parameters, including (1) the surface-paneling method, (2) an integration method of the trajectories of the wake vortices, (3) vortex-grid refinement, and (4) the initial conditions for the first iteration on the computed aerodynamic coefficients and on the flow-field details are presented. The convergence of the iterative-solution procedure is usually rapid. The solution converges with grid refinement to a constant value, but the final value is not unique and varies with the wing surface-paneling and wake-discretization methods within some range in the vicinity of the experimental result.
Carlson, H. W.
Sonic boom overpressures and signature duration may be predicted for the entire affected ground area for a wide variety of supersonic airplane configurations and spacecraft operating at altitudes up to 76 km in level flight or in moderate climbing or descending flight paths. The outlined procedure relies to a great extent on the use of charts to provide generation and propagation factors for use in relatively simple expressions for signature calculation. Computational requirements can be met by hand-held scientific calculators, or even by slide rules. A variety of correlations of predicted and measured sonic-boom data for airplanes and spacecraft serve to demonstrate the applicability of the simplified method.
Wells, William L.
Two scaled models of the Aeroassist Flight Experiment (AFE) vehicle were tested in two air wind tunnels and one CF4 tunnel. The tests were to determine the static longitudinal aerodynamic characteristics, and shock shapes for the configuration in hypersonic continuum flow. The tests were conducted with a range of angle of attack to evaluate the effects of Mach number, Reynolds numbers, and normal shock density ratio.
Mei, Ren-Wei; Shyy, Wei; Yu, Da-Zhi; Luo, Li-Shi; Rudy, David (Technical Monitor)
The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is
Rowe, W. S.; Sebastian, J. D.; Petrarca, J. R.
Results of theoretical and numerical investigations conducted to develop economical computing procedures were applied to an existing computer program that predicts unsteady aerodynamic loadings caused by leading and trailing edge control surface motions in subsonic compressible flow. Large reductions in computing costs were achieved by removing the spanwise singularity of the downwash integrand and evaluating its effect separately in closed form. Additional reductions were obtained by modifying the incremental pressure term that account for downwash singularities at control surface edges. Accuracy of theoretical predictions of unsteady loading at high reduced frequencies was increased by applying new pressure expressions that exactly satisified the high frequency boundary conditions of an oscillating control surface. Comparative computer result indicated that the revised procedures provide more accurate predictions of unsteady loadings as well as providing reduction of 50 to 80 percent in computer usage costs.
Isakova, N. P.; Kraiko, A. A.; P'yankov, K. S.
A direct optimization method for a broad class of three-dimensional aerodynamic shapes based on the approximation of the desired geometry by Bernstein-Bézier surfaces is developed. The high efficiency of the method is demonstrated by applying it to the design of an optimal supersonic section of an axisymmetric maximum-thrust de Laval nozzle. The method is also tested as applied to the design of a three-dimensional supersonic nozzle section in a dense multi-nozzle setup. In addition to three-dimensional supersonic nozzle sections with a circular throat, nozzles with a varying throat shape are considered. The results suggest that the method can be applied to various problems of 3D shape optimization.
Cao, Chunqin; Xiang, Jingbo; Zhang, Xiaoyang; Wang, Weiqiang
The infrared (IR) imaging missile's dome will be heated when fly at high speed in the atmosphere because of the friction of the air flow blocking. The detector's performance will be decline if the dome surface is heated to a certain temperature. In this paper, we find a right way to evaluate the aerothermal effects in the imaging and information processing algorithm. Which have three steps including the aerothermal radiation calculation, quantization and image reconstruction. Firstly, the aerothermal radiation is calculated by using a combination of both methods of theoretical analysis and experiment data. Secondly, the relationship between aerothermal radiation and IR images background mean gray and noise can be calculated through the analysis of the experiment data. At last, we can rebuild an aerodynamic heating effect of infrared images fusion with target and decoy, which can be used for virtual prototyping platform missile trajectory simulation. It can be found that the above constructed images have good agreements with the actual image according to comparison between the simulation data and experiment data. It is an economic method that can solve the lab aerodynamic heating simulation and modeling problems.
Gard, Eric; Riot, Vincent; Coffee, Keith; Woods, Bruce; Tobias, Herbert; Birch, Jim; Weisgraber, Todd
A computerized method and system for designing an aerodynamic focusing lens stack, using input from a designer related to, for example, particle size range to be considered, characteristics of the gas to be flowed through the system, the upstream temperature and pressure at the top of a first focusing lens, the flow rate through the aerodynamic focusing lens stack equivalent at atmosphere pressure; and a Stokes number range. Based on the design parameters, the method and system determines the total number of focusing lenses and their respective orifice diameters required to focus the particle size range to be considered, by first calculating for the orifice diameter of the first focusing lens in the Stokes formula, and then using that value to determine, in iterative fashion, intermediate flow values which are themselves used to determine the orifice diameters of each succeeding focusing lens in the stack design, with the results being output to a designer. In addition, the Reynolds numbers associated with each focusing lens as well as exit nozzle size may also be determined to enhance the stack design.
Mohammad, A. F.; Zaki, S. A.; Hagishima, A.; Ali, M. S. M.
Estimates of aerodynamic parameters, in particular roughness length z 0 and displacement height d, are important for the analysis of the roughness of an urban surface, which affects processes that occur within the urban boundary layer such as pollutant dispersion and urban ventilation. Findings regarding the aerodynamic effects of various configurations of urban arrays were compiled from various studies. Several experimental, numerical and semi-empirical studies to estimate z 0 and d were reviewed and compared with each other. The results can be summarized as follows: (1) the influence of the frontal area index ( λ f ) on z 0 is significant and their relationship has been confirmed by both experimental and numerical data; (2) compared to one-parameter and two-parameter fitting methods, the three-parameter fitting method is the least accurate; (3) the physical meaning of d remains vague because its definition as the height where surface drag acts may not be accurate for sharp-edged roughness blocks and (4) the peak values of z 0 for uniform and heterogeneous block heights indicate presence of skimming or wake-interference flow effects, which may influence surface roughness. Finally, the semi-empirical models were found to be limited to cases derived from available experimental data, which normally involve uniform arrays of cubes.
Mendenhall, Michael R.; Lesieutre, Daniel J.; Whittaker, C. H.; Curry, Robert E.; Moulton, Bryan
Pegasus, a three-stage, air-launched, winged space booster was developed to provide fast and efficient commercial launch services for small satellites. The aerodynamic design and analysis of Pegasus was conducted without benefit of wind tunnel tests using only computational aerodynamic and fluid dynamic methods. Flight test data from the first two operational flights of Pegasus are now available, and they provide an opportunity to validate the accuracy of the predicted pre-flight aerodynamic characteristics. Comparisons of measured and predicted flight characteristics are presented and discussed. Results show that the computational methods provide reasonable aerodynamic design information with acceptable margins. Post-flight analyses illustrate certain areas in which improvements are desired.
Mehta, Daryush D.; Hillman, Robert E.
Purpose of review This paper describes recent advances in perceptual, acoustic, aerodynamic, and endoscopic imaging methods for assessing voice production. Recent findings Perceptual assessment Speech-language pathologists are being encouraged to use the new CAPE-V inventory for auditory perceptual assessment of voice quality, and recent studies have provided new insights into listener reliability issues that have plagued subjective perceptual judgments of voice quality. Acoustic assessment Progress is being made on the development of algorithms that are more robust for analyzing disordered voices, including the capability to extract voice quality-related measures from running speech segments. Aerodynamic assessment New devices for measuring phonation threshold air pressures and air flows have the potential to serve as sensitive indices of glottal phonatory conditions, and recent developments in aeroacoustic theory may provide new insights into laryngeal sound production mechanisms. Endoscopic imaging The increased light sensitivity of new ultra high-speed color digital video processors is enabling high-quality endoscopic imaging of vocal fold tissue motion at unprecedented image capture rates, which promises to provide new insights into mechanisms of normal and disordered voice production. Summary Some of the recent research advances in voice quality assessment could be more readily adopted into clinical practice, while others will require further development. PMID:18475073
Lores, M. E.
Experience with advanced computational methods being used at the Lockheed-Georgia Company to aid in the evaluation and design of new and modified aircraft indicates that large and specialized computers will be needed to make advanced three-dimensional viscous aerodynamic computations practical. The Numerical Aerodynamic Simulation Facility should be used to provide a tool for designing better aerospace vehicles while at the same time reducing development costs by performing computations using Navier-Stokes equations solution algorithms and permitting less sophisticated but nevertheless complex calculations to be made efficiently. Configuration definition procedures and data output formats can probably best be defined in cooperation with industry, therefore, the computer should handle many remote terminals efficiently. The capability of transferring data to and from other computers needs to be provided. Because of the significant amount of input and output associated with 3-D viscous flow calculations and because of the exceedingly fast computation speed envisioned for the computer, special attention should be paid to providing rapid, diversified, and efficient input and output.
Pandya, Mohagna J.; Baysal, Oktay
A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.
Nelson, Andrea; Aftosmis, Michael J.; Nemec, Marian; Pulliam, Thomas H.
Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the
Farassat, F.; Myers, M. K.
Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.
Farassat, F.; Myers, M. K.
Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.
Wing, L. D.
A method is presented for calculating the aerodynamic heating and shear stresses at the wall for tangent ogive noses that are slender enough to maintain an attached nose shock through that portion of flight during which heat transfer from the boundary layer to the wall is significant. The lower entropy of the attached nose shock combined with the inclusion of the streamwise pressure gradient yields a reasonable estimate of the actual flow conditions. Both laminar and turbulent boundary layers are examined and an approximation of the effects of (up to) moderate angles-of-attack is included in the analysis. The analytical method has been programmed in FORTRAN 4 for an IBM 360/91 computer.
Luckring, J. M.; Hemsch, M. J.; Morrison, J. H.
An approach is presented to treat computational aerodynamics as a process, subject to the fundamental quality assurance principles of process control and process improvement. We consider several aspects affecting uncertainty for the computational aerodynamic process and present a set of stages to determine the level of management required to meet risk assumptions desired by the customer of the predictions.
Brahimi, M.T.; Paraschivoiu, I.
The earlier aerodynamic models for studying vertical axis wind turbines (VAWT`s) are based on constant incident wind conditions and are thus capable of predicting only periodic variations in the loads. The purpose of the present study is to develop a model capable of predicting the aerodynamic loads on the Darrieus rotor in a turbulent wind. This model is based on the double-multiple streamtube method (DMS) and incorporates a stochastic wind model. The method used to simulate turbulent velocity fluctuations is based on the power spectral density. The problem consists in generating a region of turbulent flow with a relevant spectrum and spatial correlation. The first aerodynamic code developed is based on a one-dimensional turbulent wind model. However, since this model ignores the structure of the turbulence in the crossflow plane, an extension to three dimensions has been made. The computer code developed, CARDAAS, has been used to predict aerodynamic loads for the Sandia-17m rotor and compared to CARDAAV results and experimental data. Results have shown that the computed aerodynamic loads have been improved by including stochastic wind into the aerodynamic model.
A parametric study of planform and aeroelastic effects on aerodynamic center, alpha- and q- stability derivatives. Appendix C: Method for computing the aerodynamic influence coefficient matrix of nonplanar wing-body-tail configurations
Expressions are derived for computing the aerodynamic influence coefficient matrix for nonplanar wing-body-tail configurations. An aerodynamic influence coefficient is defined as the load in lbs. induced on a panel as a result of a unit angle of attack on another panel. Fuselage, wing and tail thickness are assumed to be small with the result that the thickness effect on the flow-field is negligible. The method for determining the aerodynamic influence coefficient matrix is based on the lifting solution to the small perturbation, steady potential flow equation.
Fujii, Kozo; Takanashi, Susumu
An overview of aerodynamic aircraft design methods and their recent applications in Japan is presented. A design code which was developed at the National Aerospace Laboratory (NAL) and is in use now is discussed, hence, most of the examples are the result of the collaborative work between heavy industry and the National Aerospace Laboratory. A wide variety of applications in transonic to supersonic flow regimes are presented. Although design of aircraft elements for external flows are the main focus, some of the internal flow applications are also presented. Recent applications of the design code, using the Navier Stokes and Euler equations in the analysis mode, include the design of HOPE (a space vehicle) and Upper Surface Blowing (USB) aircraft configurations.
Ishimatsu, Takuto; Morishita, Etsuo; Okunuki, Takeo; Koyama, Hisao
Flows over two circular cylinders in tandem, side-by-side, and staggered arrangements were analyzed using the overset grid method, which is capable of handling a variety of sizes and arrangements. The Reynolds number was 100 based on the cylinder diameter. The present computation code was validated by comparison with benchmark solutions for flow around a single cylinder. Wind-tunnel experiments were conducted for the side-by-side cylinder flow for comparison with numerical simulations. Calculation showed two critical spacings in the tandem arrangement where the aerodynamic forces and Strouhal number change discontinuously. Three critical spacings and four distinct flow patterns were found numerically in the side-by-side arrangement. Similar critical spacings were found in the staggered arrangement calculation and formed critical lines. Furthermore, a pocket region was found for a staggered arrangement surrounded by the critical line.
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
Automated, high-fidelity tools for aerodynamic design face critical issues in attempting to optimize real-life geometry arid in permitting radical design changes. Success in these areas promises not only significantly shorter design- cycle times, but also superior and unconventional designs. To address these issues, we investigate the use of a parmetric-CAD system in conjunction with an embedded-boundary Cartesian method. Our goal is to combine the modeling capabilities of feature-based CAD with the robustness and flexibility of component-based Cartesian volume-mesh generation for complex geometry problems. We present the development of an automated optimization frame-work with a focus on the deployment of such a CAD-based design approach in a heterogeneous parallel computing environment.
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
The objective for this paper is to present the development of an optimization capability for the Cartesian inviscid-flow analysis package of Aftosmis et al. We evaluate and characterize the following modules within the new optimization framework: (1) A component-based geometry parameterization approach using a CAD solid representation and the CAPRI interface. (2) The use of Cartesian methods in the development Optimization techniques using a genetic algorithm. The discussion and investigations focus on several real world problems of the optimization process. We examine the architectural issues associated with the deployment of a CAD-based design approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute nodes. In addition, we study the influence of noise on the performance of optimization techniques, and the overall efficiency of the optimization process for aerodynamic design of complex three-dimensional configurations. of automated optimization tools. rithm and a gradient-based algorithm.
Yates, E. Carson, Jr.; Desmarais, Robert N.
The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. The present process appears to be readily adaptable to more general surface-panel methods.
Wang, K.D.; Amano, R.S.
In designing a turbine airfoil, a continuous second derivative is desirable to prevent local flow separation and recirculation and a minimum curvature on the airfoil surface helps in decreasing turbulence level that causes potential energy loss in turbine efficiency. To meet this requirement, the design program uses the higher-order Bernstein polynomials to define for both the pressure-side and the suction-side profile curves. Two sets of data are required by the design program. The first is relevant to the geometric parameters that include radius of nose and tail, inlet and exit metal angles, length and orientation of chord, throat and pitch, bend angle, and wedge angle. The second set of data consists of five control variables to determine the suction-side profile and three control variables to determine the pressure-side profile. The program was successfully applied to modify an existing industrial airfoil. Geometric parameters were obtained by direct measurement. Trial and error procedures were employed to determine the values of those eight control variables so that the airfoil could be curve-fitted by running the design program. The shape of the airfoil reproduced by the present model was continuous and smooth. This was utilized in a subsequent computational fluid dynamics performance analysis and was revealed to be satisfactory as it eliminates certain discontinuities in the measured data. A more efficient shape of an airfoil could then be conformed by adjusting geometrical parameters and the control variables. Results obtained by using this program indicate that the present technique can not only effectively optimize an airfoil shape to achieve better aerodynamic performance but also can preserve the original design intents. A viscous flow program was used to aerodynamically analyze the original and the new designs. The produced airfoil was used to perform frequency response computations. The model developed in this study was shown to be useful for flutter
Tweedt, Daniel L.
Computational Aerodynamic simulations of an 840 ft/sec tip speed, Advanced Ducted Propulsor fan system were performed at five different operating points on the fan operating line, in order to provide detailed internal flow field information for use with fan acoustic prediction methods presently being developed, assessed and validated. The fan system is a sub-scale, lownoise research fan/nacelle model that has undergone extensive experimental testing in the 9- by 15- foot Low Speed Wind Tunnel at the NASA Glenn Research Center, resulting in quality, detailed aerodynamic and acoustic measurement data. Details of the fan geometry, the computational fluid dynamics methods, the computational grids, and various computational parameters relevant to the numerical simulations are discussed. Flow field results for three of the five operating conditions simulated are presented in order to provide a representative look at the computed solutions. Each of the five fan aerodynamic simulations involved the entire fan system, excluding a long core duct section downstream of the core inlet guide vane. As a result, only fan rotational speed and system bypass ratio, set by specifying static pressure downstream of the core inlet guide vane row, were adjusted in order to set the fan operating point, leading to operating points that lie on a fan operating line and making mass flow rate a fully dependent parameter. The resulting mass flow rates are in good agreement with measurement values. The computed blade row flow fields for all five fan operating points are, in general, aerodynamically healthy. Rotor blade and fan exit guide vane flow characteristics are good, including incidence and deviation angles, chordwise static pressure distributions, blade surface boundary layers, secondary flow structures, and blade wakes. Examination of the computed flow fields reveals no excessive boundary layer separations or related secondary-flow problems. A few spanwise comparisons between
Birckelbaw, Lourdes G.; Mcneil, Walter E.; Wardwell, Douglas A.
This report describes the aerodynamics model used in a simulation model of an advanced short takeoff and vertical landing (ASTOVL) lift-fan fighter aircraft. The simulation model was developed for use in piloted evaluations of transition and hover flight regimes, so that only low speed (M approximately 0.2) aerodynamics are included in the mathematical model. The aerodynamic model includes the power-off aerodynamic forces and moments and the propulsion system induced aerodynamic effects, including ground effects. The power-off aerodynamics data were generated using the U.S. Air Force Stability and Control Digital DATCOM program and a NASA Ames in-house graphics program called VORVIEW which allows the user to easily analyze arbitrary conceptual aircraft configurations using the VORLAX program. The jet-induced data were generated using the prediction methods of R. E. Kuhn et al., as referenced in this report.
Rumsey, Christopher L.; Wahls, Richard A.
Several recent workshops and studies are used to make an assessment of the current status of CFD for subsonic fixed wing aerodynamics. Uncertainty quantification plays a significant role in the assessment, so terms associated with verification and validation are given and some methodology and research areas are highlighted. For high-subsonic-speed cruise through buffet onset, the series of drag prediction workshops and NASA/Boeing buffet onset studies are described. For low-speed flow control for high lift, a circulation control workshop and a synthetic jet flow control workshop are described. Along with a few specific recommendations, gaps and needs identified through the workshops and studies are used to develop a list of broad recommendations to improve CFD capabilities and processes for this discipline in the future.
Abdelrahim, Mohamed E; Chrystyn, Henry
Characterization of the aerosolized dose emitted from a nebulized system can be determined using CEN (prEN13544-1) methodology and more recently with a Next Generation Impactor (NGI), but evaporative effects can influence the results. We have investigated these characteristics using different flows and cooling with the NGI and compared the results to the standard CEN method using two different nebulizer systems. The NGI was operated using flows of 15 and 30 L min(-1) at room (ROOM) temperature and immediately after cooling at 5 degrees C for 90 min (COLD). Two nebulizer systems, the Sidestream jet nebulizer (SIDE) and the Aeroneb Pro (AERO), were used to nebulize terbutaline sulphate respiratory solution. The CEN method was also used to provide the aerodynamic characteristics of the aerosolized dose from these two nebulizer systems. The mean (SD) mass median aerodynamic diameter (MMAD) using 15COLD, 15ROOM, 30COLD, 30ROOM, and CEN for AERO was 5.0(0.1), 4.1(0.3), 4.4(0.2), 2.0(0.3), and 3.0(1.1) microm, respectively, and 4.2(0.4), 2.6(0.4), 3.5(0.1), 1.7(0.1), and 3.2(0.3) microm for SIDE. The fine particle fraction (FPF), using the NGI, followed the expected trend associated with the corresponding MMAD values, ranging from 48.1 to 70.5% from AERO and 57.3 to 87.8% for SIDE. The mean FPF for AERO and SIDE using the CEN methodology was 72.5 and 63.6%. Overall there was a highly significant difference (p < 0.001) between the different operating conditions for the FPF and MMAD of both nebulizer systems. All methods revealed a significant difference between AERO and SIDE except CEN. Both nebulizer systems were prone to evaporation effects during in vitro testing. Cooling and using a slow flow minimizes evaporation effects with the NGI and should be adopted as the recommended compendial method. The CEN method provides different values to those of the NGI operating conditions and could not differentiate between the two nebulizers. PMID:19392586
Dillenius, Marnix F. E.
Program LRCDM2 was developed for supersonic missiles with axisymmetric bodies and up to two finned sections. Predicted are pressure distributions and loads acting on a complete configuration including effects of body separated flow vorticity and fin-edge vortices. The computer program is based on supersonic panelling and line singularity methods coupled with vortex tracking theory. Effects of afterbody shed vorticity on the afterbody and tail-fin pressure distributions can be optionally treated by companion program BDYSHD. Preliminary versions of combined shock expansion/linear theory and Newtonian/linear theory have been implemented as optional pressure calculation methods to extend the Mach number and angle-of-attack ranges of applicability into the nonlinear supersonic flow regime. Comparisons between program results and experimental data are given for a triform tail-finned configuration and for a canard controlled configuration with a long afterbody for Mach numbers up to 2.5. Initial tests of the nonlinear/linear theory approaches show good agreement for pressures acting on a rectangular wing and a delta wing with attached shocks for Mach numbers up to 4.6 and angles of attack up to 20 degrees.
reduction of bluff-body noise. Xiaoyu Wang and Xiaofeng Sun discuss the interaction of fan stator and acoustic treatments using the transfer element method. S Saito and his colleagues in JAXA report the development of active devices for reducing helicopter noise. The paper by A Tamura and M Tsutahara proposes a brand new methodology for aerodynamic sound by applying the lattice Boltzmann finite difference method. As the method solves the fluctuation of air density directly, it has the advantage of not requiring modeling of the sound generation. M A Langthjem and M Nakano solve the hole-tone feedback cycle in jet flow by a numerical method. Y Ogami and S Akishita propose the application of a line-vortex method to the three-dimensional separated flow from a bluff body. I hope that a second issue on aerodynamic sound will be published in FDR in the not too distant future.
Munro, C. D.; Krus, P.; Jouannet, C.
Scale effect is of significance in all experimental aerodynamics testing. At high angles of attack such issues as Reynolds and Mach number scaling are further complicated by the presence of complex, unsteady separated flow structures. The understanding and quantification of such flow effects remains a substantial challenge to the combat aircraft designer. The challenge, especially during conceptual and preliminary design stages, is to obtain sufficiently accurate information to make relevant design decisions, addressing potential weaknesses before proceeding to detail design. The focus in the present review is on the application of the water tunnel in such cases. The combination of qualitative and quantitative data obtainable from the water tunnel provides a useful complement to other tools during the early design stages. The effect of Reynolds number in particular is however critical; and appreciation of its effects are critical to the effective use of such a tool. It is shown however that for certain configurations such effects can be small in comparison to other experimental effects.
Iungo, Giacomo Valerio; Viola, Francesco; Camarri, Simone; Porté-Agel, Fernando; Gallaire, Francois
Instability of the hub vortex, which is a vorticity structure present in wind turbine near-wake and mainly oriented along the streamwise direction, is predicted from wake velocity measurements. In this work, stability analysis is performed on wind tunnel velocity measurements acquired in the wake produced from a wind turbine model immersed in a uniform flow. Turbulence effects on wake dynamics are taken into account by modeling the Reynolds stresses through eddy-viscosity models, which are calibrated on the wind tunnel data. This formulation leads to the identification of one dominant mode associated with the hub vortex instability, which is characterized by a counter-winding single-helix mode. Moreover, this analysis also predicts accurately the frequency of the hub vortex instability observed experimentally. The hub vortex instability is also investigated by considering incoming wind fields with different turbulence characteristics, different turbine aerodynamic designs and operational regimes, which affect the morphology of the wake vorticity structures and their dynamics. The ultimate goal of this work consists in providing useful information for predicting wind turbine wake dynamics and their effects on downstream wake recovery, thus to maximize wind power harvesting.
Cumming, Stephen B.; Smith, Mark S.; Ali, Aliyah N.; Bui, Trong T.; Ellsworth, Joel C.; Garcia, Christian A.
The aerodynamic effects of compliant flaps installed onto a modified Gulfstream III airplane were investigated. Analyses were performed prior to flight to predict the aerodynamic effects of the flap installation. Flight tests were conducted to gather both structural and aerodynamic data. The airplane was instrumented to collect vehicle aerodynamic data and wing pressure data. A leading-edge stagnation detection system was also installed. The data from these flights were analyzed and compared with predictions. The predictive tools compared well with flight data for small flap deflections, but differences between predictions and flight estimates were greater at larger deflections. This paper describes the methods used to examine the aerodynamics data from the flight tests and provides a discussion of the flight-test results in the areas of vehicle aerodynamics, wing sectional pressure coefficient profiles, and air data.
Glassman, A. J.; Moffitt, T. P.
A cursory review is presented of some of the recent work that has been done in turbine aerodynamic research at NASA-Lewis Research Center. Topics discussed include the aerodynamic effect of turbine coolant, high work-factor (ratio of stage work to square of blade speed) turbines, and computer methods for turbine design and performance prediction. An extensive bibliography is included. Experimental cooled-turbine aerodynamics programs using two-dimensional cascades, full annular cascades, and cold rotating turbine stage tests are discussed with some typical results presented. Analytically predicted results for cooled blade performance are compared to experimental results. The problems and some of the current programs associated with the use of very high work factors for fan-drive turbines of high-bypass-ratio engines are discussed. Turbines currently being investigated make use of advanced blading concepts designed to maintain high efficiency under conditions of high aerodynamic loading. Computer programs have been developed for turbine design-point performance, off-design performance, supersonic blade profile design, and the calculation of channel velocities for subsonic and transonic flow fields. The use of these programs for the design and analysis of axial and radial turbines is discussed.
Carlson, Harry W.; McElroy, Marcus O.; Lessard, Wendy B.; McCullers, L. Arnold
Prediction of the loss of wing leading-edge thrust and the accompanying increase in drag due to lift, when flow is not completely attached, presents a difficult but commonly encountered problem. A method (called the previous method) for the prediction of attainable leading-edge thrust and the resultant effect on airplane aerodynamic performance has been in use for more than a decade. Recently, the method has been revised to enhance its applicability to current airplane design and evaluation problems. The improved method (called the present method) provides for a greater range of airfoil shapes from very sharp to very blunt leading edges. It is also based on a wider range of Reynolds numbers than was available for the previous method. The present method, when employed in computer codes for aerodynamic analysis, generally results in improved correlation with experimental wing-body axial-force data and provides reasonable estimates of the measured drag.
Nemec, M.; Aftosmis, M. J.; Pulliam, T. H.
The objective for this paper is to present the development of an optimization capability for Curt3D, a Cartesian inviscid-flow analysis package. We present the construction of a new optimization framework and we focus on the following issues: 1) Component-based geometry parameterization approach using parametric-CAD models and CAPRI. A novel geometry server is introduced that addresses the issue of parallel efficiency while only sparingly consuming CAD resources; 2) The use of genetic and gradient-based algorithms for three-dimensional aerodynamic design problems. The influence of noise on the optimization methods is studied. Our goal is to create a responsive and automated framework that efficiently identifies design modifications that result in substantial performance improvements. In addition, we examine the architectural issues associated with the deployment of a CAD-based approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute engines. We demonstrate the effectiveness of the framework for a design problem that features topology changes and complex geometry.
Gregory, T. J.
While faster computers will be needed to make solution of the Navier-Stokes equations practical and useful, most all of the other aerodynamic solution techniques can benefit from faster computers. There is a wide variety of computational and measurement techniques, the prospect of more powerful computers permits extension and an enhancement across all aerodynamic methods, including wind-tunnel measurement. It is expected that, as in the past, a blend of methods will be used to predict aircraft aerodynamics in the future. These will include methods based on solution of the Navier-Stokes equations and the potential flow equations as well as those based on empirical and measured results. The primary flows of interest in aircraft aerodynamics are identified, the predictive methods currently in use and/or under development are reviewed and two of these methods are analyzed in terms of the computational resources needed to improve their usefulness and practicality.
Sherman, L.L.; Taylor, A.C. III; Hou, G.W.; Korivi, V.M.
The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained. The basic equations for second-order sensitivity derivatives are presented, which results in a comparison of four different methods. Each of these four schemes for second-order derivatives requires that large systems are solved first for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first-order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent-flow example. In each of these two sample problems, three dependent variables (coefficients of lift, drag, and pitching-moment) and six independent variables (three geometric-shape and three flow-condition design variables) are considered. Several different procedures are tested, and results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method. Furthermore, it is at least two to four times faster than central finite differences, without an overwhelming penalty in computer memory. 23 refs., 14 tabs.
Unstructured grid Euler computations, performed at supersonic cruise speed, are presented for a High Speed Civil Transport (HSCT) configuration, designated as the Technology Concept Airplane (TCA) within the High Speed Research (HSR) Program. The numerical results are obtained for the complete TCA cruise configuration which includes the wing, fuselage, empennage, diverters, and flow through nacelles at M (sub infinity) = 2.4 for a range of angles-of-attack and sideslip. Although all the present computations are performed for the complete TCA configuration, appropriate assumptions derived from the fundamental supersonic aerodynamic principles have been made to extract aerodynamic predictions to complement the experimental data obtained from a 1.675%-scaled truncated (aft fuselage/empennage components removed) TCA model. The validity of the computational results, derived from the latter assumptions, are thoroughly addressed and discussed in detail. The computed surface and off-surface flow characteristics are analyzed and the pressure coefficient contours on the wing lower surface are shown to correlate reasonably well with the available pressure sensitive paint results, particularly, for the complex flow structures around the nacelles. The predicted longitudinal and lateral/directional performance characteristics for the truncated TCA configuration are shown to correlate very well with the corresponding wind-tunnel data across the examined range of angles-of-attack and sideslip. The complementary computational results for the longitudinal and lateral/directional performance characteristics for the complete TCA configuration are also presented along with the aerodynamic effects due to empennage components. Results are also presented to assess the computational method performance, solution sensitivity to grid refinement, and solution convergence characteristics.
Vandam, C. P. G.; Griswold, M.; Roskam, J.
An analytical method is presented for predicting the lift coefficient, the pitching moment coefficient, and the drag coefficient of light, twin-engine, propeller-driven airplanes. The method was applied to the Advanced Technology Light Twin-Engine airplane. The calculated characteristics were then correlated against full scale wind tunnel data. The analytical method was found to predict the drag and pitching moment fairly well. However, the lift prediction was extremely poor.
Batina, John T.
The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.
Boyd, David D. Jr.
Preliminary aerodynamic and performance predictions for an active twist rotor for a HART-II type of configuration are performed using a computational fluid dynamics (CFD) code, OVERFLOW2, and a computational structural dynamics (CSD) code, CAMRAD -II. These codes are loosely coupled to compute a consistent set of aerodynamics and elastic blade motions. Resultant aerodynamic and blade motion data are then used in the Ffowcs-Williams Hawkins solver, PSU-WOPWOP, to compute noise on an observer plane under the rotor. Active twist of the rotor blade is achieved in CAMRAD-II by application of a periodic torsional moment couple (of equal and opposite sign) at the blade root and tip at a specified frequency and amplitude. To provide confidence in these particular active twist predictions for which no measured data is available, the rotor system geometry and computational set up examined here are identical to that used in a previous successful Higher Harmonic Control (HHC) computational study. For a single frequency equal to three times the blade passage frequency (3P), active twist is applied across a range of control phase angles at two different amplitudes. Predicted results indicate that there are control phase angles where the maximum mid-frequency noise level and the 4P non -rotating hub vibrations can be reduced, potentially, both at the same time. However, these calculated reductions are predicted to come with a performance penalty in the form of a reduction in rotor lift-to-drag ratio due to an increase in rotor profile power.
Ding, Lei; Zhang, JinBi; Zheng, HaiYang; Wang, YingPing; Fang, Li
For the purpose of classification of airborne particles, this paper describes an experimental apparatus for simultaneously measuring shape characteristics and aerodynamic size at single particle level. The shape of a particle is indicated through near forward scattering light collected by 3 PMTs placed at 120-degree offset azimuthal angles and the aerodynamic diameter is obtained by time-of-flight that a particle takes to traverse double laser beams. Laboratory experiments are performed on sampled aerosol particles in spherical, cuboid and elongated shape, and preliminary results indicate that the experimental apparatus has a good capability of discriminating between spherical and irregular particles. A variance factor of scattered light related to shape of ambient airborne particles under different conditions are also presented, which can be modeled using lognormal probability density distribution. Combined with aerodynamic size information, these results suggest potential uses in environmental aerosol monitoring for characterizing constituents of particles.
Rege, Alok Ashok
Insect flight comes with a lot of intricacies that cannot be explained by conventional aerodynamics. Even with their small-size, insects have the ability to generate the required aerodynamic forces using high frequency flapping motion of their wings to perform different maneuvers. The maneuverability obtained by these flyers using flapping motion belies the classical aerodynamics theory and calls for a new approach to study this highly unsteady aerodynamics. Research is on to find new ways to realize the flight capabilities of these insects and engineer a micro-flyer which would have various applications, ranging from autonomous pollination of crop fields and oil & gas exploration to area surveillance and detection & rescue missions. In this research, a parametric study of flapping trajectories is performed using a two-dimensional wing to identify the factors that affect the force production. These factors are then non-dimensionalized and used in a design of experiments set-up to conduct sensitivity analysis. A procedure to determine an aerodynamic model comprising cycle-averaged force coefficients is described. This aerodynamic model is then used in a nonlinear dynamics framework to perform flight dynamics analysis using a micro-flyer with model properties based on Drosophila. Stability analysis is conducted to determine different steady state flight conditions that could achieved by the micro-flyer with the given model properties. The effect of scaling the mass properties is discussed. An LQR design is used for closed-loop control. Open and closed-loop simulations are performed. The results show that nonlinear dynamics framework can be used to determine values for model properties of a micro-flyer that would enable it to perform different flight maneuvers.
Gangwani, S. T.
A detailed analysis of dynamic stall experiments has led to a set of relatively compact analytical expressions, called synthesized unsteady airfoil data, which accurately describe in the time-domain the unsteady aerodynamic characteristics of stalled airfoils. An analytical research program was conducted to expand and improve this synthesized unsteady airfoil data method using additional available sets of unsteady airfoil data. The primary objectives were to reduce these data to synthesized form for use in rotor airload prediction analyses and to generalize the results. Unsteady drag data were synthesized which provided the basis for successful expansion of the formulation to include computation of the unsteady pressure drag of airfoils and rotor blades. Also, an improved prediction model for airfoil flow reattachment was incorporated in the method. Application of this improved unsteady aerodynamics model has resulted in an improved correlation between analytic predictions and measured full scale helicopter blade loads and stress data.
Tweedt, Daniel L.
Computational Aerodynamic simulations of a 1215 ft/sec tip speed transonic fan system were performed at five different operating points on the fan operating line, in order to provide detailed internal flow field information for use with fan acoustic prediction methods presently being developed, assessed and validated. The fan system is a sub-scale, low-noise research fan/nacelle model that has undergone extensive experimental testing in the 9- by 15-foot Low Speed Wind Tunnel at the NASA Glenn Research Center. Details of the fan geometry, the computational fluid dynamics methods, the computational grids, and various computational parameters relevant to the numerical simulations are discussed. Flow field results for three of the five operating points simulated are presented in order to provide a representative look at the computed solutions. Each of the five fan aerodynamic simulations involved the entire fan system, which for this model did not include a split flow path with core and bypass ducts. As a result, it was only necessary to adjust fan rotational speed in order to set the fan operating point, leading to operating points that lie on a fan operating line and making mass flow rate a fully dependent parameter. The resulting mass flow rates are in good agreement with measurement values. Computed blade row flow fields at all fan operating points are, in general, aerodynamically healthy. Rotor blade and fan exit guide vane flow characteristics are good, including incidence and deviation angles, chordwise static pressure distributions, blade surface boundary layers, secondary flow structures, and blade wakes. Examination of the flow fields at all operating conditions reveals no excessive boundary layer separations or related secondary-flow problems.
Nielsen, Jack N.
The fundamental aerodynamics of slender bodies is examined in the reprint edition of an introductory textbook originally published in 1960. Chapters are devoted to the formulas commonly used in missile aerodynamics; slender-body theory at supersonic and subsonic speeds; vortices in viscid and inviscid flow; wing-body interference; downwash, sidewash, and the wake; wing-tail interference; aerodynamic controls; pressure foredrag, base drag, and skin friction; and stability derivatives. Diagrams, graphs, tables of terms and formulas are provided.
Hamilton, J. T.; Wallace, R. O.; Dill, C. C.
The aerodynamic development plan for the Space Shuttle integrated vehicle had three major objectives. The first objective was to support the evolution of the basic configuration by establishing aerodynamic impacts to various candidate configurations. The second objective was to provide continuing evaluation of the basic aerodynamic characteristics in order to bring about a mature data base. The third task was development of the element and component aerodynamic characteristics and distributed air loads data to support structural loads analyses. The complexity of the configurations rendered conventional analytic methods of little use and therefore required extensive wind tunnel testing of detailed complex models. However, the ground testing and analyses did not predict the aerodynamic characteristics that were extracted from the Space Shuttle flight test program. Future programs that involve the use of vehicles similar to the Space Shuttle should be concerned with the complex flow fields characteristics of these types of complex configurations.
Greer, Donald; Harmory, Phil; Krake, Keith; Drela, Mark
A sailplane being developed at NASA Dryden Flight Research Center will support a high-altitude flight experiment. The experiment will measure the performance parameters or an airfoil at high altitudes (70,000 - 100,000 ft), low Reynolds numbers (2 x 10(exp 5) - 7 x 10(exp 5)), and high subsonic Mach numbers (0.5 and 0.65). The airfoil section lift and drag are determined from pilot and static pressure measurements. The locations of the separation bubble, Tollmien-Schlichting boundary-layer instability frequencies, and vortex shedding are measured from a hot-film strip. The details of the planned flight experiment are presented as well as several predictions of the airfoil performance.
Greer, Donald; Hamory, Phil; Krake, Keith; Drela, Mark
A sailplane being developed at NASA Dryden Flight Research Center will support a high-altitude flight experiment. The experiment will measure the performance parameters of an airfoil at high altitudes (70,000 to 100,000 ft), low Reynolds numbers (200,000 to 700,000), and high subsonic Mach numbers (0.5 and 0.65). The airfoil section lift and drag are determined from pitot and static pressure measurements. The locations of the separation bubble, Tollmien-Schlichting boundary layer instability frequencies, and vortex shedding are measured from a hot-film strip. The details of the planned flight experiment are presented. Several predictions of the airfoil performance are also presented. Mark Drela from the Massachusetts Institute of Technology designed the APEX-16 airfoil, using the MSES code. Two-dimensional Navier-Stokes analyses were performed by Mahidhar Tatineni and Xiaolin Zhong from the University of California, Los Angeles, and by the authors at NASA Dryden.
Pearson, H A
An empirical method is given for estimating the aerodynamic effect of ordinary and split flaps on airfoils similar to the Clark Y. The method is based on a series of charts that have been derived from an analysis of existing wind-tunnel data. Factors are included by which such variables as flap location, flap span, wing aspect ratio, and wing taper may be taken into account. A series of comparisons indicate that the method would be suitable for use in making preliminary performance calculations and in structural design.
Jorgensen, L. H.
An engineering-type procedure is presented for computing normal-force, axial-force, and pitching-moment coefficients for bodies at angles of attack from 0 deg to 180 deg. The procedure is ideally suited for estimating the aerodynamic characteristics of space shuttle booster-like bodies because of the wide range of angles of attack, Mach numbers, and Reynolds numbers to be encountered. The analytical formulas, plots, and references given are also applicable for shuttle orbiter, missile, and aircraft-like bodies of both circular and noncircular cross section. The method for computing normal-force and pitching-moment coefficients is based upon the original proposal of Allen that the crossflow or lift distribution over a body can be expressed as the sum of slender-body potential term and an empirical viscous crossflow term. Although experimental data from which to verify the procedure at very high angles of attack are extremely limited, the comparisons made thus far of computed with experimental results are good. In this report the procedure has been shown to be capable of predicting reasonably well the experimental variation of C sub N, C sub A, C sub m, and x sub ac/l with angle of attack for nine bodies of revolution at a free-stream Mach number of 2.86.
Gurbacki, Phillip M.
In this work, a feasibility study of a novel technique for the real-time identification of uncertain nonlinear aircraft aerodynamic coefficients has been conducted. The major objective of this paper is to investigate the feasibility of a system for parameter identification in a real-time flight environment. This system should be able to calculate aerodynamic coefficients and derivative information using typical pilot inputs while ensuring robust, stable, and rapid convergence. The parameter estimator investigated is based upon the nonlinear sliding mode control schema; one of the main advantages of the sliding mode estimator is the ability to guarantee a stable and robust convergence. Stable convergence is ensured by choosing a sliding surface and function that satisfies the Lyapunov stability criteria. After a proper sliding surface has been chosen, the nonlinear equations of motion for an F-16 aircraft are substituted into the sliding surface yielding an estimator capable of identifying a single aircraft parameter. Multiple sliding surfaces are then developed for each of the different flight parameters that will be identified. Sliding surfaces and parameter estimators have been developed and simulated for the pitching moment, lift force, and drag force coefficients of the F-16 aircraft. Comparing the estimated coefficients with the reference coefficients shows rapid and stable convergence for a variety of pilot inputs. Starting with simple doublet and sin wave commands, and followed by more complicated continuous pilot inputs, estimated aerodynamic coefficients have been shown to match the actual coefficients with a high degree of accuracy. This estimator is also shown to be superior to model reference or adaptive estimators, it is able to handle positive and negative estimated parameters and control inputs along with guaranteeing Lyapunov stability during convergence. Accurately estimating these aerodynamic parameters in real-time during a flight is essential
The Hypersonic Arbitrary-Body Aerodynamics (HABA) computer program predicts static and dynamic aerodynamic derivatives at hypersonic speeds for any vehicle geometry. It is intended to be used during conceptual design studies where fast computational speed is required. It uses the same geometry and hypersonic aerodynamic methods as the Mark IV Supersonic/Hypersonic Arbitrary-Body Program (SHABP) developed under sponsorship of the Air Force Flight Dynamics Laboratory; however, the input and output formats have been improved to make it easier to use. This program is available as part of the Department 9140 CAE software.
Henderson, W. P.
Powerplant installations involve complex flows, strongly influenced by viscous effects and often with important aerodynamic interactions between the airframe and propulsion system. The introduction of vehicle propulsion concepts, and points of emphasis in aircraft and missile design requirements, provide an expanding range of aerodynamic problems which cal for both experimental and theoretical study. Aerodynamic problems in powerplant installation are surveyed and work which has improved basic understanding or has enhanced prediction and design methods in this field is reviewed. Powerplant installation effects for both combat and transport aircraft are emphasized.
Meakin, Robert L.; Nygaard, Tor A.
A generic spinning missile with dithering canards is used to demonstrate the utility of an overset structured grid approach for simulating the aerodynamics of rolling airframe missile systems. The approach is used to generate a modest aerodynamic database for the generic missile. The database is populated with solutions to the Euler and Navier-Stokes equations. It is used to evaluate grid resolution requirements for accurate prediction of instantaneous missile loads and the relative aerodynamic significance of angle-of-attack, canard pitching sequence, viscous effects, and roll-rate effects. A novel analytical method for inter- and extrapolation of database results is also given.
Ammara, I.; Masson, C.; Paraschivoiu, I.
In most existing performance-analysis methods, wind turbines are considered isolated so that interference effects caused by other rotors or by the site topography are neglected. The main objective of this paper is to propose a practical 3-D method suitable for the study of these effects, in order to optimize the arrangement and the positioning of Horizontal-Axis Wind Turbines (HAWTs) in a wind farm. In the proposed methodology, the flow field around isolated HAWTs is predicted by solving the 3-D, time-averaged, steady-state, incompressible, Navier-Stokes equations in which the turbines are represented by distributions of momentum sources. The resulting governing equations are solved using a Control-Volume Finite Element Method (CVFEM). The fundamental aspects related to the development of a practical 3-D method are discussed in this paper, with an emphasis on some of the challenges that arose during its implementation. The current implementation is limited to the analysis of isolated HAWTs. Preliminary results have indicated that, the proposed 3-D method reaches the same level of accuracy, in terms of performance predictions, that the previously developed 2-D axisymmetric model and the well-known momentum-strip theory, while still using reasonable computers resources. It can be considered as a useful tool for the design of HAWTs. Its main advantages, however, are its intrinsic capacity to predict the details of the flow in the wake, and its capabilities of modelling arbitrary wind-turbine arrangements and including ground effects.
Duque, Earl P. N.; Johnson, Wayne; vanDam, C. P.; Chao, David D.; Cortes, Regina; Yee, Karen
Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.
Yates, E. C., Jr.; Cunningham, H. J.; Desmarais, R. N.; Silva, W. A.; Drobenko, B.
The SOUSSA (steady, oscillatory, and unsteady subsonic and supersonic aerodynamics) program is the computational implementation of a general potential flow analysis (by the Green's function method) that can generate pressure distributions on complete aircraft having arbitrary shapes, motions and deformations. Some applications of the initial release version of this program to several wings in steady and oscillatory motion, including flutter are presented. The results are validated by comparisons with other calculations and experiments. Experiences in using the program as well as some recent improvements are described.
Maughmer, M.; Straussfogel, D.; Long, L.; Ozoroski, L.
This work examines the ability of the aerodynamic analysis methods contained in an industry standard conceptual design code, the Aerodynamic Preliminary Analysis System (APAS II), to estimate the forces and moments generated through control surface deflections from low subsonic to high hypersonic speeds. Predicted control forces and moments generated by various control effectors are compared with previously published wind-tunnel and flight-test data for three vehicles: the North American X-15, a hypersonic research airplane concept, and the Space Shuttle Orbiter. Qualitative summaries of the results are given for each force and moment coefficient and each control derivative in the various speed ranges. Results show that all predictions of longitudinal stability and control derivatives are acceptable for use at the conceptual design stage.
Gern, Frank H.
This paper provides an overview of ongoing efforts to develop, evaluate, and validate different tools for improved aerodynamic modeling and systems analysis of Hybrid Wing Body (HWB) aircraft configurations. Results are being presented for the evaluation of different aerodynamic tools including panel methods, enhanced panel methods with viscous drag prediction, and computational fluid dynamics. Emphasis is placed on proper prediction of aerodynamic loads for structural sizing as well as viscous drag prediction to develop drag polars for HWB conceptual design optimization. Data from transonic wind tunnel tests at the Arnold Engineering Development Center s 16-Foot Transonic Tunnel was used as a reference data set in order to evaluate the accuracy of the aerodynamic tools. Triangularized surface data and Vehicle Sketch Pad (VSP) models of an X-48B 2% scale wind tunnel model were used to generate input and model files for the different analysis tools. In support of ongoing HWB scaling studies within the NASA Environmentally Responsible Aviation (ERA) program, an improved finite element based structural analysis and weight estimation tool for HWB center bodies is currently under development. Aerodynamic results from these analyses are used to provide additional aerodynamic validation data.
Gao, YingYing; He, Feng; Shen, MengYu
Based on the idea of adjoint method and the dynamic evolution method, a new optimum aerodynamic design technique is presented in this paper. It can be applied to the optimum problems with a large number of design variables and is time saving. The key of the new method lies in that the optimization process is regarded as an unsteady evolution, i.e., the optimization is executed, simultaneously with solving the unsteady flow governing equations and adjoint equations. Numerical examples for both the inverse problem and drag minimization using Euler equations have been presented, and the results show that the method presented in this paper is more efficient than the optimum methods based on the steady flow solution and the steady solution of adjoint equations.
Stone, J. R.
A method is provided for predicting jet noise for a wide range of nozzle geometries and operating conditions of interest for aircraft engines. Jet noise theory, data and existing prediction methods was reviewed, and based on this information a interim method of jet noise prediction is proposed. Problem areas are idenified where further research is needed to improve the prediction method. This method predicts only the noise generated by the exhaust jets mixing with the surrounding air and does not include other noises emanating from the engine exhaust, such as combustion and machinery noise generated inside the engine (i.e., core noise). It does, however, include thrust reverser noise. Prediction relations are provided for conical nozzles, plug nozzles, coaxial nozzles and slot nozzles.
Hemsch, M. J.; Nielsen, J. N.
A method has been developed for estimating the nonlinear aerodynamic characteristics of missile wing and control surfaces. The method is based on the following assumption: if a fin on a body has the same normal-force coefficient as a wing alone composed of two of the same fins joined together at their root chords, then the other force and moment coefficients of the fin and the wing alone are the same including the nonlinearities. The method can be used for deflected fins at arbitrary bank angles and at high angles of attack. In the paper, a full derivation of the method is given, its accuracy demonstrated and its use in extending missile data bases is shown.
Lesieutre, Daniel J.; Mendenhall, Michael R.; Nazario, Susana M.; Hemsch, Michael J.
An aerodynamic prediction method for missile aerodynamic performance and preliminary design has been developed to utilize a newly available systematic fin data base and an improved equivalent angle of attack methodology. The method predicts total aerodynamic loads and individual fin forces and moments for body-tail (wing-body) and canard-body-tail configurations with cruciform fin arrangements. The data base and the prediction method are valid for angles of attack up to 45 deg, arbitrary roll angles, fin deflection angles between -40 deg and 40 deg, Mach numbers between 0.6 and 4.5, and fin aspect ratios between 0.25 and 4.0. The equivalent angle of attack concept is employed to include the effects of vorticity and geometric scaling.
Kirk, Donn B.
A method is presented for obtaining the nonlinear aerodynamic stability characteristics of bodies of revolution from free-flight test.s The necessary conditions for the application of this method are: (1) that the roll rate and damping encountered in a single cycle of oscillation be small, and (2) that the resulting motion be reasonably planar. Four approximations to the nonlinear restoring moment are considered and solutions are obtained in closed form: 1. A single-term polynomial in an arbitrary power of the angle of attack. 2. A two-term polynomial having linear and cubic terms. 3. A three-term polynomial having linear, quadratic, and cubic terms. 4. A three-term polynomial having linear, quadratic, and cubic terms. An iteration procedure is formulated to allow the use of each of these approximations for obtaining the aerodynamic coefficients of bodies of revolution from free-flight test data. It is found that although the equations that are solved pertain strictly to planar motion, the solutions are applicable to motions that deviate to a fairly large degree from planar motion.
In the field of computational structural proteomics, contact predictions have shown new prospects of solving the longstanding problem of ab initio protein structure prediction. In the last few years, application of deep learning algorithms and availability of large protein sequence databases, combined with improvement in methods that derive contacts from multiple sequence alignments, have shown a huge increase in the precision of contact prediction. In addition, these predicted contacts have also been used to build three-dimensional models from scratch. In this chapter, we briefly discuss many elements of protein residue–residue contacts and the methods available for prediction, focusing on a state-of-the-art contact prediction tool, DNcon. Illustrating with a case study, we describe how DNcon can be used to make ab initio contact predictions for a given protein sequence and discuss how the predicted contacts may be analyzed and evaluated. PMID:27115648
Adhikari, Badri; Cheng, Jianlin
In the field of computational structural proteomics, contact predictions have shown new prospects of solving the longstanding problem of ab initio protein structure prediction. In the last few years, application of deep learning algorithms and availability of large protein sequence databases, combined with improvement in methods that derive contacts from multiple sequence alignments, have shown a huge increase in the precision of contact prediction. In addition, these predicted contacts have also been used to build three-dimensional models from scratch.In this chapter, we briefly discuss many elements of protein residue-residue contacts and the methods available for prediction, focusing on a state-of-the-art contact prediction tool, DNcon. Illustrating with a case study, we describe how DNcon can be used to make ab initio contact predictions for a given protein sequence and discuss how the predicted contacts may be analyzed and evaluated. PMID:27115648
Hemsch, Michael J. (Editor); Nielsen, Jack N. (Editor)
The present conference on tactical missile aerodynamics discusses autopilot-related aerodynamic design considerations, flow visualization methods' role in the study of high angle-of-attack aerodynamics, low aspect ratio wing behavior at high angle-of-attack, supersonic airbreathing propulsion system inlet design, missile bodies with noncircular cross section and bank-to-turn maneuvering capabilities, 'waverider' supersonic cruise missile concepts and design methods, asymmetric vortex sheding phenomena from bodies-of-revolution, and swept shock wave/boundary layer interaction phenomena. Also discussed are the assessment of aerodynamic drag in tactical missiles, the analysis of supersonic missile aerodynamic heating, the 'equivalent angle-of-attack' concept for engineering analysis, the vortex cloud model for body vortex shedding and tracking, paneling methods with vorticity effects and corrections for nonlinear compressibility, the application of supersonic full potential method to missile bodies, Euler space marching methods for missiles, three-dimensional missile boundary layers, and an analysis of exhaust plumes and their interaction with missile airframes.