SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This paper presents view-graphs and notes on sonic boom aging. Topic covered include sonic boom propagation, George's minimized F-function, final minimum shock boom, amplitude and age parameters, off-track aging, scaling flight test experiments, the potential for thin shocks, and results of a Boomfile flight test that showed significant waveform distortion.
The topics discussed include the following: introduction to marine mammals; marine mammal hearing; noise effects; sonic booms; sonic boom effects on marine mammals; predictions; and recommendations.
From the early-1960s until the mid-1980s, there was very little theoretical development for sonic fatigue prediction. Design nomographs based on simple theoretical models and results of specially designed tests were developed for most common aircraft structures. The use of advanced composites in the 1980s, however, generated an increased interest in development of more sophisticated theoretical models because of the possibilities for a much wider range of structural designs. The purpose of this report is to review sonic fatigue technology and, in particular, to assess recent developments. It also suggests a plan for a coordinated program of theoretical and experimental work to meet the anticipated needs of future aerospace vehicles.
Our objectives are to assess the effect of turbulence and molecular absorption (which is now known to be a key factor in sonic boom shock structure) on shaped sonic booms. Today I will discuss the combination of physical mechanisms for idealized turbulence. In parallel, we are reviewing models for mixed layer turbulence, and these physical effects will eventually be generalized.
To examine the building and human response to sonic boom in the range 3 Hz to 30 Hz, Georgia Institute of Technology is building a special acoustic driver system to simulate sonic boom. To support the NASA LaRC program on building and human response, this simulator's capability has been extended to an upper frequency of 4 KHz. A residential test house was made available by Georgia Tech for these tests. At the time of preparation of this document, most of the acoustic drivers and the associated electronics have been built and assembled. The system has, however, not been fully tested. The following pages provide an overview of the progress to date. The acoustic driver systems, and the principle of their operation together with the test house are described. Future plans are also summarized.
Three experiments were conducted to determine subjective equivalence of aircraft subsonic flyover noise and sonic booms. Two of the experiments were conducted in a loudspeaker-driven sonic boom simulator, and the third in a large room containing conventional loudspeakers. The sound generation system of the boom simulator had a frequency response extending to very low frequencies (about 1 Hz) whereas the large room loudspeakers were limited to about 20 Hz. Subjective equivalence between booms and flyovers was quantified in terms of the difference between the noise level of a boom and that of a flyover when the two were judged equally annoying. Noise levels were quantified in terms of the following noise descriptors: Perceived Level (PL), Perceived Noise Level (PNL), C-weighted sound exposure level (SELC), and A-weighted sound exposure level (SELA). Results from the present study were compared, where possible, to similar results obtained in other studies. Results showed that noise level differences depended upon the descriptor used, specific boom and aircraft noise events being compared and, except for the PNL descriptor, varied between the simulator and large room. Comparison of noise level differences obtained in the present study with those of other studies indicated good agreement across studies only for the PNL and SELA descriptors. Comparison of the present results with assessments of community response to high-energy impulsive sounds made by Working Group 84 of the National Research Council's Committee on Hearing, Bioacoustics, and Biomechanics (CHABA) showed good agreement when boom/flyover noise level differences were based on SELA. However, noise level differences obtained by CHABA using SELA for aircraft flyovers and SELC for booms were not in agreement with results obtained in the present study.
Raspet, Bass, and others have carried out a number of studies of sonic boom propagation using algorithms that can be traced back to a 1973 doctoral thesis written by Pestorius at the University of Texas. The algorithm was originally set forth without a formal derivation from a governing set of partial differential equations. Perusal of Pestorius's report and of subsequent reports from the Texas nonlinear accoustics group suggests that such a derivation has not yet appeared in the archival literature. The author and his colleagues, on the other hand, have been working with an explicit set of approximate partial differential equations analogous to Burgers' equation. The present paper traces through the pertinent statements in the Pestorius reports and shows that the derived algorithm, at least in the limit of sufficiently short steps, is equivalent to the numerical solution of a similar (and nearly equivalent) explicit set of partial differential equations. The tying-together of the two approaches is important because it is possible that some of the basic ideas embodied in the Pestorius algorithm may be incorporated in other contexts where more nearly general models of sonic boom propagation through realistic atmospheres are to be solved numerically.
The potential impact of supersonic operations includes structural damage from the sonic boom overpressure. This paper describes a study of how structural-acoustic modeling and testing techniques may be used to assess the potential for such damage in the absence of actual flyovers. Procedures are described whereby transfer functions relating structural response to sonic boom signature may be obtained with a stationary acoustic source and appropriate data processing. Further, by invoking structural-acoustic reciprocity, these transfer functions may also be acquired by measuring the radiated sound from the structure under a mechanical drive. The approach is based on the fundamental assumption of linearity, both with regard to the (acoustic) propagation of the boom in the vicinity of the structure and to the structure's response. Practical issues revolve around acoustic far field and source directivity requirements. The technique was implemented on a specially fabricated test structure at Edwards AFB, CA with the support of Wyle Laboratories, Inc. Blank shots from a cannon served as our acoustic source and taps from an instrumented hammer generated the mechanical drive. Simulated response functions were constructed. Results of comparisons with corresponding measurements recorded during dedicated supersonic flyovers with F-15 aircraft are presented for a number of sensor placements.
The sonic boom propagation codes reviewed in this study, SHOCKN and ZEPHYRUS, implement current theory on air absorption using different computational concepts. Review of the codes with a realistic atmosphere model confirm the agreement of propagation results reported by others for idealized propagation conditions. ZEPHYRUS offers greater flexibility in propagation conditions and is thus preferred for practical aircraft analysis. The ZEPHYRUS code was used to propagate sonic boom waveforms measured approximately 1000 feet away from an SR-71 aircraft flying at Mach 1.25 to 5000 feet away. These extrapolated signatures were compared to measurements at 5000 feet. Pressure values of the significant shocks (bow, canopy, inlet and tail) in the waveforms are consistent between extrapolation and measurement. Of particular interest is that four (independent) measurements taken under the aircraft centerline converge to the same extrapolated result despite differences in measurement conditions. Agreement between extrapolated and measured signature duration is prevented by measured duration of the 5000 foot signatures either much longer or shorter than would be expected. The duration anomalies may be due to signature probing not sufficiently parallel to the aircraft flight direction.
A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.
A diamond wing and body combination was designed to have an area distribution which would result in near optimum zero-lift wave-drag coefficients at a Mach number of 1.00, and decreasing wave-drag coefficient with increasing Mach number up to near sonic leading-edge conditions for the wing. The airfoil section were computed by varying their shape along with the body radii (blending process) to match the selected area distribution and the given plan form. The exposed wing section had an average maximum thickness of about 3 percent of the local chords, and the maximum thickness of the center-line chord was 5.49 percent. The wing had an aspect ratio of 2 and a leading-edge sweep of 45 deg. Test data were obtained throughout the Mach number range from 0.20 to 3.50 at Reynolds numbers based on the mean aerodynamic chord of roughly 6,000,000 to 9,000,000. The zero-lift wave-drag coefficients of the diamond model satisfied the design objectives and were equal to the low values for the Mach number 1.00 equivalent body up to the limit of the transonic tests. From the peak drag coefficient near M = 1.00 there was a gradual decrease in wave-drag coefficient up to M = 1.20. Above sonic leading-edge conditions of the wing there was a rise in the wave-drag coefficient which was attributed in part to the body contouring as well as to the wing geometry. The diamond model had good lift characteristics, in spite of the prediction from low-aspect-ratio theory that the rear half of the diamond wing would carry little lift. The experimental lift-curve slope obtained at supersonic speeds were equal to or greater than the values predicted by linear theory. Similarly the other basic aerodynamic parameters, aerodynamic center position, and maximum lift-drag ratios were satisfactorily predicted at supersonic speeds.
This paper presents a discussion of the sensitivity of supersonic-transport design and operation to sonic-boom considerations and shows the necessity for a study of these problems early in the development program. Methods of predicting pressure signatures are outlined and examples of the correlation of these estimates with wind-tunnel and flight measurements are shown. Estimates of sonic-boom characteristics for a representative supersonic transport show that in the critical transonic acceleration portion of the flight, overpressures somewhat lower than estimated by the use of far-field assumptions may be expected. Promising design possibilities for the achievement of further overpressure reductions are explored.
Flow field solutions over the Mars Pathfinder Probe spanning the trajectory through the Martian atmosphere at angles of attack from 0 to 11 degrees are obtained. Aerodynamic coefficients derived from these solutions reveal two regions where the derivative of pitching moment with respect to angle of attack is positive at small angles of attack. The behavior is associated with the transition of the sonic line location between the blunted nose and the windside shoulder of the 70 degree half-angle cone in a gas with a low effective ratio of specific heats. The transition first occurs as the shock layer gas chemistry evolves from highly nonequilibrium to near equilibrium, above approximately 6.5 km/s and 40 km altitude, causing the effective specific heat ratio to decrease. The transition next occurs in an equilibrium flow regime as velocities decrease through 3.5 km/s and the specific heat ratio increases again with decreasing enthalpy. The effects of the expansion over the shoulder into the wake are more strongly felt on the fustrum when the sonic line sits on the shoulder. The transition also produces a counter-intuitive trend in which windside heating levels decrease with increasing angle of attack resulting from an increase in the effective radius of curvature. Six-degree-of-freedom trajectory analyses utilizing the computed aerodynamic coefficients predict a moderate, 3 to 4 degree increase in total angle of attack as the probe, spinning at approximately 2 revolutions per minute, passes through these regions.
A method for estimating the sonic-boom overpressures from a conceptual aircraft where the lift is carried by both a canard and a wing during supersonic cruise is presented and discussed. Computer codes used for the prediction of the aerodynamic performance of the wing, the canard-wing interference, the nacelle-wing interference, and the sonic-boom overpressures are identified and discussed as the procedures in the method are discussed. A canard-wing supersonic-cruise concept was used as an example to demonstrate the application of the method.
A computational study of transverse sonic injection of air and helium into a Mach 1.98 cross-flow is presented. A hybrid large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) turbulence model is used, with the two-equation Menter baseline (Menter-BSL) closure for the RANS part of the flow and a Smagorinsky-type model for the LES part of the flow. A time-dependent blending function, dependent on modeled turbulence variables, is used to shift the closure from RANS to LES. Turbulent structures are initiated and sustained through the use of a recycling / rescaling technique. Two higher-order discretizations, the Piecewise Parabolic Method (PPM) of Colella and Woodward, and the SONIC-A ENO scheme of Suresh and Huyhn are used in the study. The results using the hybrid model show reasonably good agreement with time-averaged Mie scattering data and with experimental surface pressure distributions, even though the penetration of the jet into the cross-flow is slightly over-predicted. The LES/RANS results are used to examine the validity of commonly-used assumptions of constant Schmidt and Prandtl numbers in the intense mixing zone downstream of the injection location.
43 sonic booms generated (a few were evanescent waves) a) Overpressures of 0.08 to 2.20 lbf/sq ft; b) Rise-times of about 0.7 to 50 ms. Objectives: a) Structural response of a house of modern construction; b) Sonic boom propagation code validation. Approach: a) Measure shockwave directionality; b) Determine effect of height above ground on acoustic level; c) Generate atmospheric turbulence filter functions.
A cut-cell approach to Computational Fluid Dynamics (CFD) that utilizes the median dual of a tetrahedral background grid is described. The discrete adjoint is also calculated, which permits adaptation based on improving the calculation of a specified output (off-body pressure signature) in supersonic inviscid flow. These predicted signatures are compared to wind tunnel measurements on and off the configuration centerline 10 body lengths below the model to validate the method for sonic boom prediction. Accurate mid-field sonic boom pressure signatures are calculated with the Euler equations without the use of hybrid grid or signature propagation methods. Highly-refined, shock-aligned anisotropic grids were produced by this method from coarse isotropic grids created without prior knowledge of shock locations. A heuristic reconstruction limiter provided stable flow and adjoint solution schemes while producing similar signatures to Barth-Jespersen and Venkatakrishnan limiters. The use of cut-cells with an output-based adaptive scheme completely automated this accurate prediction capability after a triangular mesh is generated for the cut surface. This automation drastically reduces the manual intervention required by existing methods.
A detailed geometric description, in wave drag format, has been developed for the Convair B-58 and North American XB-70-1 delta wing airplanes. These descriptions have been placed on electronic files, the contents of which are described in this paper They are intended for use in wave drag and sonic boom calculations. Included in the electronic file and in the present paper are photographs and 3-view drawings of the two airplanes, tabulated geometric descriptions of each vehicle and its components, and comparisons of the electronic file outputs with existing data. The comparisons include a pictorial of the two airplanes based on the present geometric descriptions, and cross-sectional area distributions for both the normal Mach cuts and oblique Mach cuts above and below the vehicles. Good correlation exists between the area distributions generated in the late 1950s and 1960s and the present files. The availability of these electronic files facilitates further validation of sonic boom prediction codes through the use of two existing data bases on these airplanes, which were acquired in the 1960s and have not been fully exploited.