ALBERT

Description: The work to be described was performed at the NASA Langley UPWT (4-ft supersonic), test section #2, during 21-24 May 1996. The configuration being tested was the 1.675% Ref H controls model; test conditions were Ma = 2.40, Re = 3 million/ft. This was an exploration of a new technique, and it was not intended to provide definitive comparison of measured and computed skin friction results. It is, however, hoped that the experience gained will make such a rigorous comparison possible in the future.

Description: To summarize the significant highlights in this report: (1) Data quality, determined by multiple repeat runs performed on the TCA baseline configuration, and long-term repeatability, determined by comparing baseline Reference H data from this test to a previous test, have been shown to be good. (2) The longitudinal stability of the TCA is more non-linear than for the Reference H, and while it is similar at normal lift values, the TCA has considerably more pitch-up at higher lift. (3) Longitudinal control effectiveness of the TCA is similar to the Reference H and the ratio of elevator effectiveness to horizontal tail effectiveness is approximately 0.3. 4) The directional stability of the TCA is improved relative to Reference H at higher angles-of attack. The chine is effective for improving directional stability. (5) The directional control effectiveness 'of the TCA rudder is the same as that of the Reference H rudder at low angles-of-attack, after taking factors, such as number of rudder panels deflected and vertical tail volume into account. However, rudder effectiveness was shown to be reduced at higher angles-of-attack. (6) The lateral stability was shown to be reduced relative to the Reference H, which may be beneficial at low speeds for alleviating lateral control saturation. (7) Lateral control effectiveness for the TCA was shown to be similar to the Reference H for negative trailing-edge flap deflections and was reduced by approximately 25% for positive trailing-edge flap deflections.

Description: This paper gives the results of a grid study, a turbulence model study, and a Reynolds number effect study for transonic flows over a high-speed aircraft using the thin-layer, upwind, Navier-Stokes CFL3D code. The four turbulence models evaluated are the algebraic Baldwin-Lomax model with the Degani-Schiff modifications, the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, and Menter's two-equation Shear-Stress-Transport (SST) model. The flow conditions, which correspond to tests performed in the NASA Langley National Transonic Facility (NTF), are a Mach number of 0.90 and a Reynolds number of 30 million based on chord for a range of angle-of-attacks (1 degree to 10 degrees). For the Reynolds number effect study, Reynolds numbers of 10 and 80 million based on chord were also evaluated. Computed forces and surface pressures compare reasonably well with the experimental data for all four of the turbulence models. The Baldwin-Lomax model with the Degani-Schiff modifications and the one-equation Baldwin-Barth model show the best agreement with experiment overall. The Reynolds number effects are evaluated using the Baldwin-Lomax with the Degani-Schiff modifications and the Baldwin-Barth turbulence models. Five angles-of-attack were evaluated for the Reynolds number effect study at three different Reynolds numbers. More work is needed to determine the ability of CFL3D to accurately predict Reynolds number effects.

Description: The NASA High Speed Research (HSR) Program is intended to establish a technology base enabling industry development of an economically viable and environmentally acceptable second generation high speed civil transport (HSCT). The HSR program consists of work directed towards several broad technology areas, one of which is aerodynamic performance. The objective of the Configuration Aerodynamics task of the Aerodynamic Performance technology area is the development of aerodynamic drag reduction, stability and control, and propulsion airframe integration technologies required to support the HSCT development process. Towards this goal, computational and empirical based aerodynamic design tools are being developed, evaluated, and validated through ground based experimental testing. In addition, methods for ground to flight scaling are being developed and refined. Successful development of validated design and scaling methodologies will result in improved economy of operation for an HSCT and reduce uncertainty in full-scale flight predictions throughout the development process.

Description: The purpose of the chemistry component of the model comparison is to assess to what extent differences in the formulation of chemical processes explain the variance between model results. Observed concentrations of chemical compounds are used to estimate to what degree the various models represent realistic situations. For readability, the materials for the chemistry experiment are reported in three separate sections. This section discussed the data used to evaluate the models in their simulation of the source gases and the Nitrogen compounds (NO(y)) and Chlorine compounds (Cl(y)) species.

Description: It is not unusual when comparing CFD data to experimental data to find discrepancies between the results. Sometimes forces and moments compare well, while surface pressures do not, and vice versa. It is commonplace for the researcher to believe that the flow field has been accurately simulated when these types of measurements compare well. However, being able to routinely predict boundary layer transition and separated flows are not guaranteed. In fact accurate simulation of these types of flow physics has been a challenge to the CFD community. In order to improve Navier-Stokes predictions for complex vortical flow fields, more detailed information about the flow physics is necessary. Unfortunately, the many wind-tunnel tests performed in Langley's NTF and 14x22 facilities as well as in the Ames' 12 ft. Tunnel provided little information about the detailed flow physics, and no priority was given to obtaining any CFD measurements. Using the latest experimental techniques, this information can and should be obtained for present and future use.

Description: Space-based and airborne coherent Doppler lidars designed for measuring global tropospheric wind profiles in cloud-free air rely on backscatter, beta from aerosols acting as passive wind tracers. Aerosol beta distribution in the vertical can vary over as much as 5-6 orders of magnitude. Thus, the design of a wave length-specific, space-borne or airborne lidar must account for the magnitude of 8 in the region or features of interest. The SPAce Readiness Coherent Lidar Experiment under development by the National Aeronautics and Space Administration (NASA) and scheduled for launch on the Space Shuttle in 2001, will demonstrate wind measurements from space using a solid-state 2 micrometer coherent Doppler lidar. Consequently, there is a critical need to understand variability of aerosol beta at 2.1 micrometers, to evaluate signal detection under varying aerosol loading conditions. Although few direct measurements of beta at 2.1 micrometers exist, extensive datasets, including climatologies in widely-separated locations, do exist for other wavelengths based on CO2 and Nd:YAG lidars. Datasets also exist for the associated microphysical and chemical properties. An example of a multi-parametric dataset is that of the NASA GLObal Backscatter Experiment (GLOBE) in 1990 in which aerosol chemistry and size distributions were measured concurrently with multi-wavelength lidar backscatter observations. More recently, continuous-wave (CW) lidar backscatter measurements at mid-infrared wavelengths have been made during the Multicenter Airborne Coherent Atmospheric Wind Sensor (MACAWS) experiment in 1995. Using Lorenz-Mie theory, these datasets have been used to develop a method to convert lidar backscatter to the 2.1 micrometer wavelength. This paper presents comparison of modeled backscatter at wavelengths for which backscatter measurements exist including converted beta (sub 2.1).

Description: The typical fair weather electric field at the ground is between -100 and -300 V/m. At the NASA Kennedy Space Center and US Air Force Cape Canaveral Air Station (KSC) the electric field at the ground sometimes reaches -400 to -1200 V/m within an hour or two after sunrise on days that otherwise seem to be fair weather. We refer to the enhanced negative electric fields as the "sunrise enhancement." To investigate the sunrise enhancement at KSC we measured the electric field (E) in the first few hundred meters above the ground before and during several sunrise enhancements. From these E soundings we can infer the presence of charge layers and determine their thickness and charge density.

Description: This paper focuses on the parallel computation of aerodynamic derivatives via automatic differentiation of the Euler/Navier-Stokes solver CFL3D. The comparison with derivatives obtained by finite differences is presented and the scaling of the time required to obtain the derivatives relative to the number of processors employed for the computation is shown. Finally, the derivative computations are coupled with an optimizer and surface/volume grid deformation tools to perform an optimization to reduce the drag of a three-dimensional wing.

Description: Current parallel computational approaches involve distributed and shared memory paradigms. In the distributed memory paradigm, each processor has its own independent memory. Message passing typically uses a function library such as MPI or PVM. In the shared memory paradigm, such as that used on the SGI Origin 2000 machine, compiler directives are used to instruct the compiler to schedule multiple threads to perform calculations. In this paradigm, it must be assured that processors (threads) do not simultaneously access regions of memory in such away that errors would occur. This paper utilizes the latest version of the SGI MPI function library to combine the two parallelization paradigms to perform aerodynamic shape optimization of a generic wing/body.

Description: The primary objectives of this study were to expand the data base showing the effects of LE radius distribution and corresponding sensitivity to Rn at subsonic and transonic conditions, and to assess the predictive capability of CFD for these effects. Several key elements led to the initiation of this project: 1) the necessity of meeting multipoint design requirements to enable a viable HSCT, 2) the demonstration that blunt supersonic leading-edges can be associated with performance gain at supersonic speeds , and 3) limited data. A test of a modified Reference H model with the TCA planform and 2 LE radius distributions was performed in the NTF, in addition to Navier-Stokes analysis for an additional 3 LE radius distributions. Results indicate that there is a tremendous potential to improve high-lift performance through the use of a blunt LE across the span given an integrated, fully optimized design, and that low Rn data alone is probably not sufficient to demonstrate the benefit.

Description: This paper presents results of a study which attempted to provide some understanding of the relationship between skin friction drag estimates produced by flat plate methods and those produced by Navier-Stokes computations. A brief introduction is followed by analysis, including a flat plate grid study, analysis of the wing flow, an analysis of the fuselage flow. Other results of interest are then presented, including turbulence model sensitivities, and brief analysis of other configurations.

Description: Efforts towards understanding boundary layer transition characteristics on a High Speed Civil Transport (HSCT)-class configuration in the National Transonic Facility (NTF) are ongoing. The majority of the High Speed Research (HSR) data base in the NTF has free transition on the wing, even at low Reynolds numbers (Rn) attainable in conventional facilities. Limited data has been obtained and is described herein showing the effects of a conventional, Braslow method based wing boundary-layer trip on drag. Comparisons are made using force data polars and surface flow visualization at selected angles-of-attack and Mach number. Minimum drag data obtained in this study suggest that boundary layer transition occurred very near the wing leading edge by a chord Rn of 30 million. Sublimating chemicals were used in the air mode of operation only at low Rn and low angles-of-attack with no flap deflections; sublimation results suggest that the forebody and outboard wing panel are the only regions with significant laminar flow. The process and issues related to the sublimating chemical technique as applied in the NTF are discussed. Beyond the existing experience, status of efforts to develop a production transition detection system applicable to both air and cryogenic nitrogen environments is presented.

Description: Model deformation measurement techniques have been investigated and developed at NASA's Langley Research Center. The current technique is based upon a single video camera photogrammetric determination of two dimensional coordinates of wing targets with a fixed (and known) third dimensional coordinate, namely the spanwise location. Variations of this technique have been used to measure wing twist and bending at a few selected spanwise locations near the wing tip on HSR models at the National Transonic Facility, the Transonic Dynamics Tunnel, and the Unitary Plan Wind Tunnel. Automated measurements have been made at both the Transonic Dynamics Tunnel and at Unitary Plan Wind Tunnel during the past year. Automated measurements were made for the first time at the NTF during the recently completed HSR Reference H Test 78 in early 1996. A major problem in automation for the NTF has been the need for high contrast targets which do not exceed the stringent surface finish requirements. The advantages and limitations (including targeting) of the technique as well as the rationale for selection of this particular technique are discussed. Wing twist examples from the HSR Reference H model are presented to illustrate the run-to-run and test-to-test repeatability of the technique in air mode at the NTF. Examples of wing twist in cryogenic nitrogen mode at the NTF are also presented.

Description: To develop full scale flight performance predictions an understanding of Reynolds number effects on HSCT-class configurations is essential. A wind tunnel database utilizing a 2.2% scale Reference H model in NASA Langley Research Centers National Transonic Facility is being developed to assess these Reynolds number effects. In developing this database temperature and aeroelastic corrections to the wind tunnel data have been identified and are being analyzed. Once final corrections have been developed and applied, then pure Reynolds number effects can be determined. In addition, final corrections will yield the data required for CFD validation at q = 0. Presented in this report are the results of seven tests involving the wing/body configuration. This includes summaries of data acquired in these tests, uncorrected Reynolds number effects, and temperature and aeroelastic corrections. The data presented herein illustrates the successes achieved to date as well as the challenges that will be faced in obtaining full scale flight performance predictions.

Description: Experience with afterbody closure effects and accompanying test techniques issues on a High Speed Civil Transport (HSCT)-class configuration is described. An experimental data base has been developed which includes force, moment, and surface pressure data for the High Speed Research (HSR) Reference H configuration with a closed afterbody at subsonic and transonic speeds, and with a cylindrical afterbody at transonic and supersonic speeds. A supporting computational study has been performed using the USM3D unstructured Euler solver for the purposes of computational fluid dynamics (CFD) method assessment and model support system interference assessment with a focus on lower blade mount effects on longitudinal data at transonic speeds. Test technique issues related to a lower blade sting mount strategy are described based on experience in the National Transonic Facility (NTF). The assessment and application of the USM3D code to the afterbody/sting interference problem is discussed. Finally, status and plans to address critical test technique issues and for continuation of the computational study are presented.

Description: The Boeing Reference H configuration was tested in the NASA Ames 9x7 Supersonic Wind Tunnel. A simulated unstarted inlet was evaluated as well as the aerodynamic performance of the configuration with and without nacelle and diverter components. These experimental results were compared with computational results from the unstructured grid Euler flow solver AIRPLANE. The comparisons between computational and experimental results were good, and demonstrated that the Euler code is capable of efficiently and accurately predicting the changes in the aerodynamic coefficients associated with inlet unstart and the effects of the nacelle and diverter components.

Description: This presentation will describe the organization and conduct of the workshops, list the topics discussed, and conclude with a more-detailed examination of a related set of issues dear to the presenters heart. Because the current HSCT configuration is expected to have (mostly) turbulent flow over the wings, and because current CFD predictions assume fully-turbulent flow, the wind tunnel testing to date has attempted to duplicate this condition at the lower Reynolds numbers attainable on the ground. This frequently requires some form of artificial boundary layer trip to induce transition near the wing's leading edge. But this innocent-sounding goal leads to a number of complications, and it is not clear that present-day testing technology is adequate to the task. An description of some of the difficulties, and work underway to address them, forms the "Results" section of this talk. Additional results of the testing workshop will be covered in presentations by other team members.

Description: This paper presents The Propulsion Airframe Integration Advisory report in viewgraph form. The approach of the advisory group is to identify and prioritize technology elements (1.0 Inlet Issues, 2.0 Nozzle Issues, 3.0 Nacelle Design, and 4.0 Airframe Integration).

Description: Preliminary human acceptability studies of sonic booms indicate that supersonic flight is unlikely to be acceptable even at noise levels significantly below 1994 low boom designs (reference 1, p. 288). Further, these low boom designs represent considerable changes to baseline configurations, and changes translate into additional effort and uncertain structural weight penalties that may provide no annoyance benefit, increasing the risk of including low boom technology. Since over land sonic boom designs were so risky (and yet the acceptability studies highlight how annoying sonic booms are), boom softening studies were undertaken to reduce the boom of baseline configurations using minor modifications that would not significantly change the designs. The goal of this work is to reduce boom levels over water. Even though Concorde over water boom has not been found to have any adverse environmental impact, boom levels for baseline HSCT designs are 50% higher in overpressure than the Concorde (due to a doubling in configuration weight with only a 50% increase in length),

Description: Officially, the Tu-144 was the first supersonic-cruise, passenger-carrying aircraft to enter commercial service. Design, construction, and testing were carried out by the Soviet Union, flight certification was by the Soviet Union, and the only regular passenger flights were scheduled and flown across the territory of the Soviet Union. Although it was not introduced to international passenger service, there were many significant engineering accomplishments achieved in the design, production, and flight of this aircraft. Development of the aircraft began with a prototype stage. Systematic testing and redesign led to a production aircraft in discrete stages that measurably improved the performance of the aircraft from the starting concept to final aircraft certification. It flew in competition with the English-French Concorde for a short time, but was withdrawn from national commercial service due to a lack of interest by airlines outside the Soviet Union. NASA became interested in the Tu- 144 aircraft when it was offered for use as a flying "testbed" in the study of operating characteristics of a supersonic-cruise commercial airplane. Since it had been in supersonic-cruise service, the Tu- 144 had operational characteris'tics similar to those anticipated in the conceptual aircraft designs being studied by the United States aircraft companies. In addition to the other operational tests being conducted on the Tu-144 aircraft, it was proposed that two sets of sonic-boom pressure signature measurements be made. The first set would be made on the ground, using techniques and devices similar to those in reference I and many other subsequent studies. A second set would be made in the air with an instrumented aircraft flying close under the Tu-144 in supersonic flight. Such in-flight measurements would require pressure gages that were capable of accurately recording the flow-field overpressures generated by the Tu- 144 at relatively close distances under the vehicle. Therefore, an analysis of the Tu-144 was made to obtain predictions of pressure signature shape and shock strengths at cruise conditions so that the range and characteristics of the required pressure gages could be determined well in advance of the tests. Cancellation of the sonic-boom signature measurement part of the tests removed the need for these pressure gages. Since CFD methods would be used to analyze the aerodynamic performance of the Tu-144 and make similar pressure signature predictions, the relatively quick and simple Whitham-theory pressure signature predictions presented in this paper could be used for comparisons. Pressure signature predictions of sonic-boom disturbances from the Tu- 144 aircraft were obtained from geometry derived from a three-view description of the production aircraft. The geometry was used to calculate aerodynamic performance characteristics at supersonic-cruise conditions. These characteristics and Whitham/Walkden sonic-boom theory were employed to obtain F-functions and flow-field pressure signature predictions at a Mach number of 2.2, at a cruise altitude of 61000 feet, and at a cruise weight of 350000 pounds.

Description: This document contains the details of the thermal analysis of the X-38 aft fin during re-entry. This analysis was performed in order to calculate temperature response of the aft fin components. This would be provided as input to a structural analysis and would also define the operating environment for the electromechanical actuator (EMA). The calculated structural temperature response would verify the performance of the thermal protection system (TPS). The geometric representation of the aft fin was derived from an I-DEAS finite element model that was used for structural analysis. The thermal mass network model was derived from the geometric representation.

Description: To summarize the significant highlights in this report: (1) Data quality, determined by multiple repeat runs performed on the TCA baseline configuration, and long-term repeatability, determined by comparing baseline Reference H data from this test to a previous test, have been shown to be good. (2) The longitudinal stability of the TCA is more non-linear than for the Reference H, and while it is similar at normal lift values, the TCA has considerably more pitch-up at higher lift. (3) Longitudinal control effectiveness of the TCA is similar to the Reference H and the ratio of elevator effectiveness to horizontal tail effectiveness is approximately 0.3. (4) The directional stability of the TCA is improved relative to Reference H at higher angles-of attack. The chine is effective for improving directional stability.

Description: The NASA-industry team has sponsored several studies in the last two years to address the installed nozzle boattail drag issues. Some early studies suggested that nozzle boattail drag could be as much as 25 to 40 percent of the subsonic cruise. As part of this study tests have been conducted at NASA-Langley to determine the uninstalled drag characteristics of a proposed nozzle. The overall objective was to determine the effects of nozzle external flap curvature and sidewall boattail variations. This test would also provide data for validating CFD predictions of nozzle boattail drag.

Description: AIRPLANE (Jameson/Baker) is a steady inviscid unstructured Euler flow solver. It has been validated on many HSR geometries. It is implemented as MESHPLANE, an unstructured mesh generator, and FLOPLANE, an iterative flow solver. The surface description from an Intergraph CAD system goes into MESHPLANE as collections of polygonal curves to generate the 3D mesh. The flow solver uses a multistage time stepping scheme with residual averaging to approach steady state, but R is not time accurate. The flow solver was ported from Cray to IBM SP2 by Wu-Sun Cheng (IBM); it could only be run on 4 CPUs at a time because of memory limitations. Meshes for the four cases had about 655,000 points in the flow field, about 3.9 million tetrahedra, about 77,500 points on the surface. The flow solver took about 23 wall seconds per iteration when using 4 CPUs. It took about eight and a half wall hours to run 1,300 iterations at a time (the queue limit is 10 hours). A revised version of FLOPLANE (Thomas) was used on up to 64 CPUs to finish up some calculations at the end. We had to turn on more communication when using more processors to eliminate noise that was contaminating the flow field; this added about 50% to the elapsed wall time per iteration when using 64 CPUs. This study involved computing lift and drag for a wing/body/nacelle configuration at Mach 0.9 and 4 degrees pitch. Four cases were considered, corresponding to four nacelle mass flow conditions.

Description: The primary objectives of this study were to expand the data base showing the effects of LE radius distribution and corresponding . sensitivity to Rn at subsonic and transonic conditions, and to assess the predictive capability of CFD for these effects. Several key elements led to the initiation of this project: 1) the necessity of meeting multipoint design requirements to enable a viable HSCT, 2) the demonstration that blunt supersonic leading-edges can be associated with performance gain at supersonic speeds , and 3) limited data. A test of a modified Reference H model with the TCA planform and 2 LE radius distributions was performed in the NTF, in addition to Navier-Stokes analysis for an additional 3 LE radius distributions. Results indicate that there is a tremendous potential to improve high-lift performance through the use of a blunt LE across the span given an integrated, fully optimized design, and that low Rn data alone is probably not sufficient to demonstrate the benefit.

Description: In cooperation with personnel from the Boeing ANP Laboratory and NASA Langley, a performance test was conducted using the Reference-H 1.675% model ("NASA Modular Model") without nacelles at the NASA Langley 16-Ft Transonic Tunnel. The main objective of the test was to determine the drag reduction achievable with leading-edge and trailing-edge flaps deflected along the outboard wing span at transonic Mach numbers (M = 0.9 to 1.2) for purpose of preliminary design and for comparison with computational predictions. The obtained drag data with flap deflections for Mach numbers of 1.07 to 1.20 are unique for the Reference H wing. Four leading-edge and two trailing-edge flap deflection angles were tested at a mean-wing chord-Reynolds number of about 5.7 million. An outboard-wing leading-edge flap deflection of 81 provides a 4.5 percent drag reduction at M = 1.2 A = 0.2), and much larger values at lower Mach numbers with larger flap deflections. The present results for the baseline (no flaps deflected) compare reasonably well with previous Boeing and NASA Ref-H tunnel tests, including high-Reynolds number NTF results. Viscous CFD simulations using the OVERFLOW thin-layer N.S. method properly predict the observed trend in drag reduction at M = 1.2 as function of leading-edge flap deflection. Modified linear theory properly predicts the flap effects on drag at subsonic conditions (Aero2S code), and properly predicts the absolute drag for the 40 and 80 leading-edge deflection at M = 1.2 (A389 code).

Description: Unstructured grid Euler computations, performed at supersonic cruise speed, are presented for a proposed high speed civil transport 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 Mach 2.4 for a range of angles-of-attack and sideslip. 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 shock wave structures around the nacelles. The predicted longitudinal and lateral/directional performance characteristics are shown to correlate very well with the measured data across the examined range of angles-of-attack and sideslip. The results from the present effort have been documented into a NASA Controlled-Distribution report which is being presently reviewed for publication.

Description: The objectives of the Cycle 2 Nonlinear Design Optimization Anlaytical Cross Checks are to: 1) Understand the variability in the predicted performance levels of the nonlinear designs arising from the use of different inviscid (full potential/Euler) and viscous (Navier-Stokes) analysis methods; and 2) Provide the information required to allow the performance levels of all three designs to be validated using the data from the NCV (nonlinear Cruise Validation) model test.

Description: During the last cycle of concept design and wind-tunnel testing, the goal of the low-boom- shaped HSCT concepts (the B-935, the LB-16, and the LB- 1 8) was to meet mission requirements and generate shaped, ground-level pressure signatures with nose shock strengths of 1.0 psf or less. The wind-tunnel tests of these concepts produced results that were partially successful and encouraging although not fully up to expectations. In spite of this, however, these conceptual designs were overly optimistic and not acceptable because: the wing planforms had excessive area; the wing structural aspect ratio was too high; one concept had aft-fuselage rather than under-the-wing engines; and the gross takeoff weights were unrealistically low because of engines that were early, high-tech versions of later, revised, more-realistic engines. The need for reducing the ground-level overpressure shock strengths still existed; a need to be met within more restrictive guidelines of mission performance and gross takeoff weight limitations. Therefore, it was decided that the next conceptual design cycle would focus on decreased nose shock strengths, "boom softening," in the signatures of the Boeing and the McDonnell Douglas baseline concepts rather than low-boom concepts with shaped-signature designs. Overly-optimistic results were not the only problem with these low-sonic-boom concepts. Papers given at the 1994 Sonic-Boom Workshop had demonstrated that the problem of successful nacelle integration on HSCT concepts had only been partially solved. Wind-tunnel pressure signature data, from the HSCT-11B (a.k.a. the LB-18) wind-tunnel model, showed that the Langley HSCT design and analysis method had been successful in reducing the nacelle-volume disturbances in the flow field. This was due.to the engine nacelles mounted behind the wing trailing-edge on the aft fuselage so that no nacelle-wing interference-lift flow-field disturbances were generated. While acceptable from a sonic-boom research point of view, this concept was unacceptable from several practical and structural considerations. Preliminary wind-tunnel pressure signature data from the LB-16 wind-tunnel model, which had the engine nacelles mounted under the wings (the usual location), indicated that the application of the Langley nacelle-integration method had been only partially successful in the reduction of the nacelle-volume with nacelle-wing interference-lift pressure disturbances. So, "boom softening" had to also address the task of successful integration of the engine nacelles, with the engines in the required under-the-wing location. Unless this problem was solved, low-sonic-boom and low-drag modifications to the wing planform, the airfoil shape, and the fuselage longitudinal area distribution could be nullified if the nacelle disturbances added increments to the nose-shock strengths that were removed through component tailoring. In this paper, an arrow-wing boom-softened HSC7 concept which incorporated modifications to a baseline McDonnell Douglas concept is discussed. The analysis of the concept's characteristics will include estimates of weight, center of gravity, takeoff field length, mission range, and predictions of its ground-level sonic-boom pressure signature. Additional modifications which enhanced the softened-boom performance of this concept are also described as well as estimates of the performance penalties induced by these modifications.

Description: A 1:300 scale wind-tunnel model of a conceptual High-Speed Civil Transport (HSCT) designed to generate a shaped, low-boom pressure signature on the ground was tested to obtain sonic-boom pressure signatures in the Langley Research Center Unitary Plan Wind Tunnel at a Mach number of 1.8 and a separation distance of about two body lengths or four wing-spans from the model. Two sets of engine nacelles representing two levels of engine technology were used on the model to determine the effects of increased nacelle volume. Pressure signatures were measured for (model lift)/(design lift) ratios of 0.5, 0.63, 0.75, and 1.0 so that the effect of lift on the pressure signature could be determined. The results of these tests were analyzed and used to discuss the agreement between experimental data and design expectations.

Description: This section contains a number of special diagnostics that are designed to examine certain mechanisms. Section 1 reports on the method used to test the photochemical partitioning in the models. Sections 2 and 3 represent efforts to examine the model calculated production and removal rates for ozone and how the values are combined with transport rates in the models to produce the simulated ozone distributions. Sections 4 and 5 concentrate on polar processes including the dynamics aspect of vortex confinement and the chemical aspects of chlorine activation.

Description: MM II defined a series of experiments to better understand and characterize model transport and to assess the realism of this transport by comparison to observations. Measurements from aircraft, balloon, and satellite, not yet available at the time of MM I [Prather and Remsberg, 1993], provide new and stringent constraints on model transport, and address the limits of our transport modeling abilities. Simulations of the idealized tracers the age spectrum, and propagating boundary conditions, and conserved HSCT-like emissions probe the relative roles of different model transport mechanisms, while simulations of SF6 and C02 make the connection to observations. Some of the tracers are related, and transport diagnostics such as the mean age can be derived from more than one of the experiments for comparison to observations. The goals of the transport experiments are: (1) To isolate the effects of transport in models from other processes; (2) To assess model transport for realistic tracers (such as SF6 and C02) for comparison to observations; (3) To use certain idealized tracers to isolate model mechanisms and relationships to atmospheric chemical perturbations; (4) To identify strengths and weaknesses of the treatment of transport processes in the models; (5) To relate evaluated shortcomings to aspects of model formulation. The following section are included:Executive Summary, Introduction, Age Spectrum, Observation, Tropical Transport in Models, Global Mean Age in Models, Source-Transport Covariance, HSCT "ANOY" Tracer Distributions, and Summary and Conclusions.

Description: The NASA High Speed Research (HSR) Program is intended to establish a technology base enabling industry development of an economically viable and environmentally acceptable second generation high speed civil transport (HSCT). The objective of the Configuration Aerodynamics task of the program is the development of aerodynamic drag reduction, stability and control, and propulsion airframe integration technologies required to support the HSCT development process. Aerodynamic design tools are being developed, evaluated, and validated through ground based experimental testing. In addition, methods for ground to flight scaling are being developed and refined.

Description: This paper reports on the model, test, and results from the Langley Supersonic Aftbody Closure wind tunnel test. This project is an experimental evaluation of the 1.5% Technology Concept Aircraft (TCA) aftbody closure model (Model 23) in the Langley Unitary Plan Wind Tunnel. The baseline TCA design is the result of a multidisciplinary, multipoint optimization process and was developed using linear design and analysis methods, supplemented with Euler and Navier-Stokes numerical methods. After a thorough design review, it was decided to use an upswept blade attached to the forebody as the mounting system. Structural concerns dictated that a wingtip support system would not be feasible. Only the aftbody part of the model is metric. The metric break was chosen to be at the fuselage station where prior aft-sting supported models had been truncated. Model 23 is thus a modified version of Model 20. The wing strongback, flap parts, and nacelles from Model 20 were used, whereas new aftbodies, a common forebody, and some new tails were fabricated. In summary, significant differences in longitudinal and direction stability and control characteristics between the ABF and ABB aftbody geometries were measured. Correcting the experimental data obtained for the TCA configuration with the flared aftbody to the representative of the baseline TCA closed aftbody will result in a significant reduction in longitudinal stability, a moderate reduction in stabilizer effectiveness and directional stability, and a moderate to significant reduction in rudder effectiveness. These reductions in the stability and control effectiveness levels of the baseline TCA closed aftbody are attributed to the reduction in carry-over area.

Description: There were two objectives for this test. First, was to assess the reasons why there is approximately 1.5 drag counts (cts) discrepancy between measured and computed drag improvement of the Non-linear Cruise Validation (NCV) over the Technology Concept Airplane (TCA) wing body (WB) configurations. The Navier-Stokes (N-S) pre-test predictions from Boeing Commercial Airplane Group (BCAG) show 4.5 drag cts of improvement for NCV over TCA at a lift coefficient (CL) of 0. I at Mach 2.4. The pre-test predictions from Boeing Phantom Works - Long Beach, BPW-LB, show 3.75 drag cts of improvement. BCAG used OVERFLOW and BPW-LB used CFL3D. The first test entry to validate the improvement was held at the NASA Langley Research Center (LARC) UPV;T, test number 1687. The experimental results showed that the drag improvement was only 2.6 cts, not accounting for laminar run and trip drag. This is approximately 1.5 cts less than predicted computationally. In addition to the low Reynolds Number (RN) test, there was a high RN test in the Boeing Supersonic Wind Tunnel (BSWT) of NCV and TCA. BSV@T test 647 showed that the drag improvement of NCV over TCA was also 2.6 cts, but this did account for laminar run and trip drag. Every effort needed to be done to assess if the improvement measured in LaRC UPWT and BSWT was correct. The second objective, once the first objective was met, was to assess the performance increment of NCV over TCA accounting for the associated laminar run and trip drag corrections in LaRC UPWT. We know that the configurations tested have laminar flow on portions of the wing and have trip drag due to the mechanisms used to force the flow to go from laminar to turbulent aft of the transition location.

Description: An improved laminar run and trip drag correction methodology for supersonic cruise performance testing was derived. This method required more careful analysis of the flow visualization images which revealed delayed transition particularly on the inboard upper surface, even for the largest trip disks. In addition, a new code was developed to estimate the laminar run correction. Once the data were corrected for laminar run, the correct approach to the analysis of the trip drag became evident. Although the data originally appeared confusing, the corrected data are consistent with previous results. Furthermore, the modified approach, which was described in this presentation, extends prior historical work by taking into account the delayed transition caused by the blunt leading edges.

Description: A global lightning model that includes diurnal and annual lightning variation, and total flash density versus latitude for each major land and ocean, has been used as the basis for simulating the global electric circuit charging rate. A particular objective has been to reconcile the difference in amplitude ratios [AR=(max-min)/mean] between global lightning diurnal variation (AR approx. = 0.8) and the diurnal variation of typical atmospheric potential gradient curves (AR approx. = 0.35). A constraint on the simulation is that the annual mean charging current should be about 1000 A. The global lightning model shows that negative ground flashes can contribute, at most, about 10-15% of the required current. For the purpose of the charging rate simulation, it was assumed that each ground flash contributes 5 C to the charging process. It was necessary to assume that all electrified clouds contribute to charging by means other than lightning, that the total flash rate can serve as an indirect indicator of the rate of charge transfer, and that oceanic electrified clouds contribute to charging even though they are relatively inefficient in producing lightning. It was also found necessary to add a diurnally invariant charging current component. By trial and error it was found that charging rate diurnal variation curves in Universal time (UT) could be produced with amplitude ratios and general shapes similar to those of the potential gradient diurnal variation curves measured over ocean and arctic regions during voyages of the Carnegie Institute research vessels.

Description: Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Description: Computations have been performed on the baseline Reference H wing/body configuration, as well as the Wing 704 configuration, an optimized wing and fuselage combination derived from Ref. H through automated optimization. The parabolized Navier-Stokes solver UPS was employed with viscous terms in two directions in an effort to understand the source and level of potential viscous/inviscid interactions. The paper briefly describes the UPS code and the grids used to obtain the solutions before the discussion of results. Results of these computations indicate that viscous/inviscid interaction can contribute increments to both the pressure- and friction-related drag. Computations were performed for wind tunnel conditions-1.675% scale models at a Reynolds number of 4 million per foot. Turbulent flow results were obtained using the Baldwin-Lomax algebraic turbulence model and were compared with laminar flow results. The laminar flow fields were used to obtain upper bounds on potential interaction effects.

Description: The NASA-industry team has sponsored several studies in the last two years to address the installed nozzle boattail drag issues. Some early studies suggested that nozzle boattail drag could be as much as 25 to 40 percent of the subsonic cruise. As part of this study tests have been conducted at NASA-Langley to determine the uninstalled drag characteristics of a proposed nozzle. The overall objective was to determine the effects of nozzle external flap curvature and sidewall boattail variations. This test would also provide data for validating CFD predictions of nozzle boattail drag.

Description: The objective of this milestone is to assess the propulsion/airframe integration characteristics of the Technology Concept Airplane and design variations through computational analysis and experimental subsonic through supersonic wind tunnel testing. The Milestone will generate a comprehensive CFD and wind tunnel data base of the baseline, and design variations. Emphasis will be placed on establishing the propulsion induced effects on the flight performance of the Technology Concept Airplane with all appropriate wind tunnel corrections.

Description: The computational fluid dynamics (CFD) comparisons being presented are compared to each other and to wind tunnel (WT) data on the baseline TCA. Some of the CFD computations were done prior to the tests and others later. Only force data (CL vs CD) from CFD will be presented as part of this report. The WT data presented comes from the testing of the baseline TCA in the Langley Unitary Plan Wind Tunnel (UPWT), Test Section #2. There are 2 sets of wind tunnel data being presented: one from test 1671 of model 2a (flapped wing) and the other from test 1679 of model 2b (solid wing). Most of the plots show only one run from each of the WT tests per configuration. But many repeat runs were taken during the tests. The WT repeat runs showed an uncertainty in the drag of +/- 0.5 count. There were times when the uncertainty in drag was better, +/- 0.25 count. Test 1671 data was of forces and pressures measured from model 2a. The wing had cutouts for installing various leading and trailing edge flaps at lower Mach numbers. The internal duct of the nacelles are not designed and fabricated as defined in the outer mold lines (OML) iges file. The internal duct was fabricated such that a linear transition occurs from the inlet to exhaust. Whereas, the iges definition has a constant area internal duct that quickly transitions from the inlet to exhaust cross sectional shape. The nacelle internal duct was fabricated, the way described, to save time and money. The variation in the cross sectional area is less than 1% from the iges definition. The nacelles were also installed with and without fairings. Fairings are defined as the build up of the nacelles on the upper wing surface so that the nacelles poke through the upper surface as defined in the OML iges file. Test 1679 data was of forces measured from model 2a and 2b. The wing for model 2b was a solid wing. The nacelles were built the same way as for model 2a, except for the nacelle base pressure installation. The nacelles were only tested with the fairings for model 2a and 2b during test 1679.

Description: Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Description: The TetrUSS (Tetrahedral Unstructured Software System), developed at NASA LaRC, enables one to take a vehicle from its surface definition to its analyzed solution. The important parts are the shape definition, accomplished in GRIDTOOL; the initial front and volume grid generation in VGRID; the flow solver USM3D, and the various ways used to post-process the computational results.

Description: The objective of this study is to calibrate a Navier-Stokes code for the TCA (30/10) baseline configuration (partial span leading edge flaps were deflected at 30 degs. and all the trailing edge flaps were deflected at 10 degs). The computational results for several angles of attack are compared with experimental force, moments, and surface pressures. The code used in this study is CFL3D; mesh sequencing and multi-grid were used to full advantage to accelerate convergence. A multi-grid approach was used similar to that used for the Reference H configuration allowing point-to-point matching across all the trailingedge block interfaces. From past experiences with the Reference H (ie, good force, moment, and pressure comparisons were obtained), it was assumed that the mounting system would produce small effects; hence, it was not initially modeled. However, comparisons of lower surface pressures indicated the post mount significantly influenced the lower surface pressures, so the post geometry was inserted into the existing grid using Chimera (overset grids).

Description: The objective of the present study was to address the questions of: 1) how reliably or consistently the Navier-Stokes methods and processes used by the various organizations can predict integrated skin friction drag, and 2) how well the methods can predict trends within a family of optimized configurations. As a first step, all available skin friction drag predictions were accumulated to obtain a mean and standard deviation for the TCA (Technology Concept Airplane) baseline and each of the optimized configurations. It is observed that the optimization process has had little effect on the predicted skin friction drags. The variation in the mean that is observed is dwarfed by the standard deviations. In order to understand the reasons for the relatively large spreads in the computed results, a number of auxiliary computations have been performed using the UPS and OVERFLOW codes in an effort to identify and quantity potential sources of the variations.

Description: The paper presents the recent progress made towards developing an efficient and user-friendly parallel environment for routine analysis of large CFD problems. The coarse-grain parallel version of the CFL3D Euler/Navier-Stokes analysis code, CFL3Dhp, has been ported onto most available parallel platforms. The CFL3Dhp solution accuracy on these parallel platforms has been verified with the CFL3D sequential analyses. User-friendly pre- and post-processing tools that enable a seamless transfer from sequential to parallel processing have been written. Static load balancing tool for CFL3Dhp analysis has also been implemented for achieving good parallel efficiency. For large problems, load balancing efficiency as high as 95% can be achieved even when large number of processors are used. Linear scalability of the CFL3Dhp code with increasing number of processors has also been shown using a large installed transonic nozzle boattail analysis. To highlight the fast turn-around time of parallel processing, the TCA full configuration in sideslip Navier-Stokes drag polar at supersonic cruise has been obtained in a day. CFL3Dhp is currently being used as a production analysis tool.

Description: This paper presents an Unstructured Navier-Stokes Analysis of Full TCA (Technology Concept Airplane) Configuration. The topics include: 1) Motivation; 2) Milestone and approach; 3) Overview of the unstructured-grid system; 4) Results on full TCA W/B/N/D/E configuration; 5) Concluding remarks; and 6) Future directions.

Description: Automatic Grid Generation Wish List Geometry handling, including CAD clean up and mesh generation, remains a major bottleneck in the application of CFD methods. There is a pressing need for greater automation in several aspects of the geometry preparation in order to reduce set up time and eliminate user intervention as much as possible. Starting from the CAD representation of a configuration, there may be holes or overlapping surfaces which require an intensive effort to establish cleanly abutting surface patches, and collections of many patches may need to be combined for more efficient use of the geometrical representation. Obtaining an accurate and suitable body conforming grid with an adequate distribution of points throughout the flow-field, for the flow conditions of interest, is often the most time consuming task for complex CFD applications. There is a need for a clean unambiguous definition of the CAD geometry. Ideally this would be carried out automatically by smart CAD clean up software. One could also define a standard piece-wise smooth surface representation suitable for use by computational methods and then create software to translate between the various CAD descriptions and the standard representation. Surface meshing remains a time consuming, user intensive procedure. There is a need for automated surface meshing, requiring only minimal user intervention to define the overall density of mesh points. The surface mesher should produce well shaped elements (triangles or quadrilaterals) whose size is determined initially according to the surface curvature with a minimum size for flat pieces, and later refined by the user in other regions if necessary. Present techniques for volume meshing all require some degree of user intervention. There is a need for fully automated and reliable volume mesh generation. In addition, it should be possible to create both surface and volume meshes that meet guaranteed measures of mesh quality (e.g. minimum and maximum angle, stretching ratios, etc.).

Description: Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Description: A team was formed to tackle the sonic boom softening issues of the current Boeing HSCT design. The team consisted of personnel from NASA Ames, NASA Langley, and Boeing company. The work described in this paper was done when the first author was at NASA Ames Research Center. This paper presents the sonic boom softening work on two Boeing High Speed Civil Transport (HSCT) baseline configurations, Reference-H and Boeing-1122. This presentation can be divided into two parts: parametric studies and sonic boom minimization by CFD optimization routines.

Description: The objectives of this research are: 1) To determine the effect of geometric variations near the inboard leading-edge flap on high-lift and stability and control performance data; 2) To determine Re effects on TCA (Technology Concept Aircraft) high-lift configuration for optimum high-lift and stability and control performance at takeoff, climbout, approach and landing conditions; and 3) To obtain flow-visualization data on upper surface of wing for CFD validations. This paper is presented in viewgraph form.

Description: This presentation describes the advances being made with the Aerodynamic Shape Optimization (ASO) and high-fidelity Multidisciplinary Optimization (MDO) software used in the High Speed Research Program at NASA Ames Research Center. The description starts with the motivation for continued ASO/MDO development. Objectives of the current work are then presented. A list of ingredients deemed necessary for a flexible design environment is discussed, and the HSR requirement for different geometries at different design points is explained. Multiple design disciplines within a high-fidelity design environment are demonstrated. Finally, progress so far is summarized and planned future work is outlined.

Description: This report considers the effect of canard and horizontal tail vertical position on the aerodynamic characteristics of the PTC configuration without nacelles and diverters. This analysis is followed by three optimization studies using canard and tail incidence as design variables in the first problem followed by an optimization run with canard and tail incidence and wing camber design variables and finally an optimization run with canard incidence and wing camber. The first problem was run at fixed lift while the other two problems were run at fixed angle of attack. The final investigation reported here will show data from a component buildup study using the PTC configuration. This final study will show the aerodynamic interference between the canard, wing and horizontal tail.

Description: This paper presents results of three minor studies into the behavior of the OVERFLOW with respect to the prediction of skin friction drag on wing bodies at cruise Mach number and wind tunnel Reynolds number. The studies include a preliminary assessment of the behavior of the two new 2-equation turbulence models introduced with the latest version of OVERFLOW (v. 1.8f), an investigation into potential improvements in the matrix dissipation scheme currently implemented in OVERFLOW, and an analysis of the observed sensitivity of the code's skin friction predictions to grid stretching at solid surface boundaries.

Description: The computational results of the optimized complete configurations, including nacelles and diverters, are presented in terms of drag count improvement compared with the TCA baseline configuration at Mach 2.4, C(sub L)=0.1. The three candidate designs are designated by the organization from which they were derived. ARC represents the Ames Research Center 1-03 design, BCAG represents the Boeing Commercial Aircraft Group's design from Seattle, and BLB represents the design from Boeing Long Beach. All CFD methods are in unanimous agreement that the Ames 1-03 configuration has the largest performance improvement, followed closely by the BCAG configuration, with a much smaller improvement attained by Boeing Long Beach. The Ames design was obtained using the single-block wing/body code SYN87-SB with its "pseudo" nacelle option-an elaborate technique for incorporating nacelle/diverter effects into the design optimization process. This technique uses AIRPLANE surface pressure coefficient data with and without the nacelles/diverters. Further details of this method are described. It is reasonable to expect that further improvements could be achieved by including the "real" nacelles directly into the optimization process by use of the newly-developed multiblock optimization code, SYN107-MB, which can handle full configurations.

Description: The aim of this work is to demonstrate a simple technique which accounts for aeroelastic deformations experienced by HSR wind-tunnel models within CFD computations. With improved correlations, CFD can become a more effective tool for augmenting the post-test understanding of experimental data. The present technique involves the loose coupling of a low-level structural representation within the ELAPS code, to an unstructured Navier-Stokes flow solver, USM3Dns. The ELAPS model is initially calibrated against bending characteristics of the wind-tunnel model. The strength of this method is that, with a single point calibration of a simple structural representation, the static aeroelastic effects can be accounted for in CFD calculations across a range of test conditions. No prior knowledge of the model deformation during the wind-on test is required. This approach has been successfully applied to the high aspect-ratio planforms of subsonic transports. The current challenge is to adapt the procedure to low aspect-ratio planforms typical of HSR configurations.

Description: This presentation includes three topics: (1) Analysis of isolated boattail drag; (2) Computation of Technology Concept Airplane (TCA)-installed nacelle effects on aerodynamic performance; and (3) Assessment of TCA inlet flow quality.

Description: LaRC conducted a code validation study for the OVERFLOW code to ascertain its accuracy for boattail drag prediction. The OVERFLOW results compared favorably with the LaRC 16-ft. Transonic Wind Tunnel (TWT) data, and prior CFD solutions from PAB3D and CFL3D. The ultimate goal is to investigate the installation drag of the nacelle boattails with powered nozzles at transonic mach numbers. The OVERFLOW solver was chosen because of its ability to accept volume overlapping structured grid for very complex airframe configurations. Structured grid components for representing the transonic nozzle boattail can be added to the BCAG grid for a TCA airframe with 2D bifurcated inlet and flow through nacelle without alteration. The focus of this research was to determine the suitability of the OVERFLOW solver for accomplishing this ultimate goal. This presentation will first introduce the transonic nozzle boattail wind-tunnel model geometry, followed by an examination of aerodynamic features based on the current OVERFLOW solutions and the solutions obtained previously using PAB3D, comparisons of Cp on the flap surface between the OVERFLOW solutions, wind tunnel data, and solutions from other CFD codes, an assessment of boattail drag count prediction, and a work plan for FY99.

Description: A concept of system identification applied to high performance aircraft is introduced followed by a discussion on the identification methodology. Special emphasis is given to model postulation using time invariant and time dependent aerodynamic parameters, model structure determination and parameter estimation using ordinary least squares and mixed estimation methods. At the same time problems of data collinearity detection and its assessment are discussed. These parts of methodology are demonstrated in examples using flight data of the X-29A and X-31A aircraft. In the third example wind tunnel oscillatory data of the F-16XL model are used. A strong dependence of these data on frequency led to the development of models with unsteady aerodynamic terms in the form of indicial functions. The paper is completed by concluding remarks.

Description: This paper reports the predicted M = 2.4 strut-interference effects on a closed aftbody with empennage for the TCA baseline model. The strut mounting technique was needed in order to assess the impact of aft-end shaping, i.e. open for a sting or closed to better represent a flight vehicle. However,this technique can potentially lead to unanticipated effects that are measured on the aft body. Therefore, a set of computations were performed in order to examine the closed aft body with and without strut present, at both zero and non-zero angles of sideslip (AOS). The work was divided into a computational task performed by Javier A. Garriz, using an inviscid (Euler) solver, and a monitoring/reporting task done by John E. Lamar. All this work was performed during FY98 at the NASA Langley Research Center.

Description: Estimates of the effect of pulse stretching on satellite laser altimetry in particular the Geoscience Laser Altimeter System (GLAS), by cloud multiple scattering were made from an analytical method and from Monte Carlo calculations. The path delay of the return pulse was found to be largest for low-level clouds with particle radii (3-20 microns). The magnitude of the path delay was affected by several factors including cloud height, cloud optical depth, cloud particle size, particle shape, and receiver field of view. Polar aerosol and Rayleigh scattering usually added less than 1 cm to the overall path delay. Path delay estimates for all cloud conditions would be less if a simple Gaussian fit of the return pulse peak were used to measure the pulse's centroid. However, care must be taken in determining the centroid as factors such as pulse width, surface slope and the fitting method used will affect the estimate. A planned application for laser altimetry is high precision monitoring of the height change of polar ice sheets. In the absence of a correction algorithm, the required GLAS altimetry accuracies will not be achieved unless atmospheric channel information is used to remove profiles that are likely to be heavily contaminated by multiple scattering.

Description: The zonal mean eddy heat flux is directly proportional to the wave activity that propagates from the troposphere into the stratosphere. This quantity is a simple eddy diagnostic which is easily calculated from conventional meteorological analyses. Because this "wave driving" of the stratosphere has a strong impact on the stratospheric temperature, it is necessary to compare the impact of the flux with respect to stratospheric radiative changes caused by greenhouse gas changes. Hence, we must understand the precision and accuracy of the heat flux derived from our global meteorological analyses. Herein, we quantify the stratospheric heat flux using five different meteorological analyses, and show that there are 30% differences between these analyses during the disturbed conditions of the northern hemisphere winter. Such large differences result from the planetary differences in the stationary temperature and meridional wind fields. In contrast, planetary transient waves show excellent agreement amongst these five analyses, and this transient heat flux appears to have a long term downward trend.

Description: We examine concurrent upper tropospheric measurements of CN (diameter greater than 4 nm). NO, and NO(Y) during the SONEX Experiment over the North Atlantic (Oct.-Nov., 1997). Elevated CN and NO(Y) concentrations observed in the upper troposphere are attributed largely to enhancements in convective outflows. We estimate that less than 7% of observed high-CN plumes (greater than 10000 /cc) may be attributed to aircraft emissions. Dilution of high-CN convective and aircraft plumes appears to be much more rapid than losses of NO(X) and CN by oxidation and coagulation, respectively, and accounts for much of observed CN concentrations. When taking into account of different time scales against dilution for observable aircraft and convective high-CN plumes (estimated to be 1:4), the contribution by aircraft emissions to CN concentrations is significant, about 20% of the convective source. We find no evidence that particle formation in convective plumes is limited by OH oxidation of SO2.

Description: This study examines a unique data set returned by IMP8 and Geotail on January 29, 1995 during a substorm which resulted in the ejection of a plasmoid. The two spacecraft (s/c) were situated in the north lobe of the tail and both observed a traveling compression region (TCR). From single s/c observations only the length of the plasmoid in X and an estimate of its height in Z can be determined. However, we show that dual s/c measurements of TCRs can be used to model all three dimensions of the underlying plasmoid and to estimate of its rate of expansion or contraction. For this event plasmoid dimensions of Delta(X) approximates 18, Delta(Y) approximates 30, and Delta(Z) approximates 10 R(sub e) are inferred from the IMP8 and Geotail lobe magnetic field measurements. The earthward end of the plasmoid was inferred to be near the mean location of the near-earth neutral line, X approximates -26 R(sub e). Its center was underneath IMP 8 at X approximates -34 R(sub e) and its tailward end appeared to be near X approximates -44 R(sub e). Furthermore, a factor of approximately 2 increase in the amplitude of the TCR occurred in the 1.5 min it took to move from IMP 8 to Geotail. Modeled using conservation of the magnetic flux, this increase in lobe compression implies that the underlying plasmoid was expanding at a rate of approximately 140 km/s. Such an expansion is comparable to recently reported V(sub y) speeds in "young" plasmoids in this region of the tail. Finally, the Geotail measurements indicate that a reconfiguration of the lobe magnetic field closely followed the ejection of the plasmoid which moved magnetic flux tubes into the wake behind the plasmoid where they would convect into the near-earth neutral line and reconnect.

Description: Prior to 1991, major warmings (defined by increasing zonal mean temperatures and zonal mean easterly winds from 60degN to the pole at 10 hPa) typically occurred approximately once every two Arctic winters; a major warming in mid-Dec. 1998 was the first since Feb. 1991. The Dec. 1998 warming was also the second earliest on record. The earliest, and the only other major warming on record before the end of Dec. was in early Dec 1987; prior to that, the earliest was in late Dec./early Jan. 1984-85. The 1984-85 and 1987 warmings resulted in the warmest and weakest lower stratospheric polar vortices in the 20 years before 1998-99. Fig. 1 compares temperatures and vortex strength in 1998-99 with those in the previous 20 years, using the US National Center for Environmental Prediction (NCEP) record; 1987-88 and 1984-85 are also highlighted. The Dec. 1998 warming had a more pronounced effect on mid-stratospheric temperatures than the Dec. 1987 warming (Fig. 1a), although smaller than that of warmings later in winter (e.g., 1984-85). 10-hPa temperatures fell well below average again in late Jan. 1999 and remained unusually low until an early final warming began in late Feb. 840 K PV gradients (Fig. 1c) set a record minimum in Jan. 1999, but were near average in Feb before the final warming. The effect of the Dec. 1998 warming on lower stratospheric temperatures was comparable to that of other major warmings; there was a brief period of record-high minimum 46-hPa temperatures in early Jan 1999 (Fig. 1b), and temperatures then fell to near average for a short period in mid-Feb. Lower stratospheric PV gradients were the weakest on record during the 1998-99 winter (Fig. 1d). The evolution of the vortex and minimum temperatures during 1998-99 was remarkably similar to that during 1987-88, the only previous year when a major warming was observed before the end of Dec.

Description: This CFD experiment concludes that the potential difference between the flow between a flight Reynolds number test and a sub-scale wind tunnel test are substantial for this particular nozzle boattail geometry. The early study was performed using a linear k-epsilon turbulence model. The present study was performed using the Girimaji formulation of a algebraic Reynolds stress turbulent simulation.

Description: The Lockheed Martin spillage study was a substantial effort and is worthy of a separate paper. However, since a paper was not submitted a few of the most pertinent results have been pulled out and included in this paper. The reader is urged to obtain a copy of the complete Boeing Configuration Aerodynamics final 1995 contract report for the complete Lockheed documentation of the spillage work. The supersonic cruise studies presented here focus on the bifurcated - axisymmetric inlet drag delta. In the process of analyzing this delta several test/CFD data correlation problems arose that lead to a correction of the measured drag delta from 4.6 counts to 3.1 counts. This study also lead to much better understanding of the OVERFLOW gridding and solution process, and to increased accuracy of the force and moment data. Detailed observations of the CFD results lead to the conclusion that the 3.1 count difference between the two inlet types could be reduced to approximately 2 counts, with an absolute lower bound of 1.2 counts due to friction drag and the bifurcated lip bevel.

Description: The Configuration Aerodynamics (CA) element of the High Speed Research (HSR) program is managed by a joint NASA and Industry team, referred to as the Technology Integration Development (ITD) team. This team is responsible for the development of a broad range of technologies for improved aerodynamic performance and stability and control characteristics at subsonic to supersonic flight conditions. These objectives are pursued through the aggressive use of advanced experimental test techniques and state of the art computational methods. As the HSR program matures and transitions into the next phase the objectives of the Configuration Aerodynamics ITD are being refined to address the drag reduction needs and stability and control requirements of High Speed Civil Transport (HSCT) aircraft. In addition, the experimental and computational tools are being refined and improved to meet these challenges. The presentation will review the work performed within the Configuration Aerodynamics element in 1994 and 1995 and then discuss the plans for the 1996-1998 time period. The final portion of the presentation will review several observations of the HSR program and the design activity within Configuration Aerodynamics.

Description: The NASA High-Speed Research program developed the High-Lift Engine Aeroacoustics Technology (HEAT) program to demonstrate satisfactory interaction between the jet noise suppressor and high-lift system of a High-Speed Civil Transport (HSCT) configuration at takeoff, climb, approach and landing conditions. One scheme for reducing jet exhaust noise generated by an HSCT is the use of a mixer-ejector system which would entrain large quantities of ambient air into the nozzle exhaust flow through secondary inlets in order to cool and slow the jet exhaust before it exits the nozzle. The effectiveness of such a noise suppression device must be evaluated in the presence of an HSCT wing high-lift system before definitive assessments can be made concerning its acoustic performance. In addition, these noise suppressors must provide the required acoustic attenuation while not degrading the thrust efficiency of the propulsion system or the aerodynamic performance of the high-lift devices on the wing. Therefore, the main objective of the HEAT program is to demonstrate these technologies and understand their interactions on a large-scale HSCT model. The HEAT program is a collaborative effort between NASA-Ames, Boeing Commercial Airplane Group, Douglas Aircraft Corp., Lockheed-Georgia, General Electric and NASA - Lewis. The suppressor nozzles used in the tests were Generation 1 2-D mixer-ejector nozzles made by General Electric. The model used was a 13.5%-scale semi-span model of a Boeing Reference H configuration.

Description: The mission of High-Lift Technology is to develop technology allowing the design of practical high lift concepts for the High-Speed Civil Transport (HSCT) in order to: 1) operate safely and efficiently; and 2) reduce terminal control area and community noise. In fulfilling this mission, close and continuous coordination will be maintained with other High-Speed Research (HSR) technology elements in order to support optimization of the overall airplane (rather than just the high lift system).

Description: Nozzle boattail drag is significant for the High Speed Civil Transport (HSCT) and can be as high as 25% of the overall propulsion system thrust at transonic conditions. Thus, nozzle boattail drag has the potential to create a thrust-drag pinch and can reduce HSCT aircraft aerodynamic efficiencies at transonic operating conditions. In order to accurately predict HSCT performance, it is imperative that nozzle boattail drag be accurately predicted. Previous methods to predict HSCT nozzle boattail drag were suspect in the transonic regime. In addition, previous prediction methods were unable to account for complex nozzle geometry and were not flexible enough for engine cycle trade studies. A computational fluid dynamics (CFD) effort was conducted by NASA and McDonnell Douglas to evaluate the magnitude and characteristics of HSCT nozzle boattail drag at transonic conditions. A team of engineers used various CFD codes and provided consistent, accurate boattail drag coefficient predictions for a family of HSCT nozzle configurations. The CFD results were incorporated into a nozzle drag database that encompassed the entire HSCT flight regime and provided the basis for an accurate and flexible prediction methodology.

Description: High-lift system performance will have a large impact on airplane noise and weight. Successful completion of PCD1 activities provided greater understanding of aerodynamic characteristics and configuration features important to high-lift system performance including: 1) Reynolds number effects (Ref. H); 2) Propulsion/airframe integration effects; and 3) Planform effects, canard/3-surface, alternate high-lift concepts, etc. PCD2 plans are aimed at achieving technology development performance goals and increasing technology readiness level for Technology Concept.

Description: Experiments were conducted in the NASA Ames 9-Ft by 7-Ft Supersonic and 11-Ft by 11-Ft Transonic Wind Tunnels of a 2.7% Reference H (Ref. H) Nacelle Airframe Interference (NAI) High Speed Civil Transport (HSCT) model. NASA Ames did the experiment with the cooperation and assistance of Boeing and McDonnell Douglas. The Ref. H geometry was designed by Boeing. The model was built and tested by NASA under a license agreement with Boeing. Detailed forces and pressures of individual components of the configuration were obtained to assess nacelle airframe interference through the transonic and supersonic flight regime. The test apparatus was capable of measuring forces and pressures of the Wing body (WB) and nacelles. Axisymmetric and 2-D inlet nacelles were tested with the WB in both the in-proximity and captive mode. The in-proximity nacelles were mounted to a nacelle support system apparatus and were individually positioned. The right hand nacelles were force instrumented with flow through strain-gauged balances and the left hand nacelles were pressure instrumented. Mass flow ratio was varied to get steady state inlet unstart data. In addition, supersonic spillage data was taken by testing the 2-D inlet nacelles with ramps and the axisymmetric inlet nacelles with an inlet centerbody for the Mach condition of interest. The captive nacelles, both axisymmetric and 2-D, were attached to the WB via diverters. The captive 2-D inlet nacelle was also tested with ramps to get supersonic spillage data. Boeing analyzed the data and showed a drag penalty of four drag counts for the 2-D compared with the axisymmetric inlet nacelle. Two of the four counts were attributable to the external bevel designed into the 2-D inlet contour. Boeing and McDonnell Douglas used these data for evaluating Computational Fluid Dynamic (CFD) codes and for evaluation of nacelle airframe integration problems and solutions.

Description: This paper discusses a method for the identification and application of reduced-order models based on linear and nonlinear aerodynamic impulse responses. The Volterra theory of nonlinear systems and an appropriate kernel identification technique are described. Insight into the nature of kernels is provided by applying the method to the nonlinear Riccati equation in a non-aerodynamic application. The method is then applied to a nonlinear aerodynamic model of an RAE 2822 supercritical airfoil undergoing plunge motions using the CFL3D Navier-Stokes flow solver with the Spalart-Allmaras turbulence model. Results demonstrate the computational efficiency of the technique.

Description: Two different methods for retrieving Upper Tropospheric Humidities (UTH) from the TOVS (TIROS Operational Vertical Sounder) instruments aboard NOAA polar orbiting satellites are presented and compared. The first one, from the Environmental Technology Laboratory, computed by J. Bates and D. Jackson (hereafter BJ method), estimates UTH from a simplified radiative transfer analysis of the upper tropospheric infrared water vapor channel at wavelength measured by HIRS (6.3 micrometer). The second one results from a neural network analysis of the TOVS (HIRS and MSU) data developed at, the Laboratoire de Meteorologie Dynamique (hereafter the 3I (Improved Initialization Inversion) method). Although the two methods give very similar retrievals in temperate regions (30-60 N and S), an absolute bias up to 16% appears in the convective zone of the tropics. The two datasets have also been compared with UTH retrievals from infrared radiance measurements in the 6.3 micrometer channel from the geostationary satellite METEOSAT (hereafter MET method). The METEOSAT retrievals are systematically drier than the TOVS-based results by an absolute bias between 5 and 25%. Despite the biases, the spatial and temporal correlations are very good. The purpose of this study is to explain the deviations observed between the three datasets. The sensitivity of UTH to air temperature and humidity profiles is analysed as are the clouds effects. Overall, the comparison of the three retrievals gives an assessment of the current uncertainties in water vapor amounts in the upper troposphere as determined from NOAA and METEOSAT satellites.

Description: Carbonados are porous polycrystalline (with crystal sizes up to 100 micrometer) diamonds. Carbonado is found only in alluvial deposits in Bahia, Brazil and in the Central African Republic (CAR). Alluvial deposit host is 1.5 Ga while the carbonados are between 2.6 - 3.8 Ga. The process of fusing the carbonado microcrystals together is not fully understood, partly due to fact that the origin of these carbonado, is not known. Several modes of origins are proposed for carbonado. First, a crustal origin. Carbonados have a light carbon and helium isotopic signature. They contain an enrichment of the rare-earth elements (REE). Carbonados have tightly trapped atmospheric noble gases and contain an evidence of high He content despite the carbonado expected depletion of He at mantle temperatures. Carbonados have high porosity incompatible with high pressure mantle conditions. Second, a mantle origin is proposed based on similar REE pattern to kimberlites. The presence of nitrogen platelet (by IR spectra) indicates high temperature origin and syngenetic inclusions of rutile, ilmenite, and magnetite indicates high pressure and high temperature conditions consistent with mantle origin as well. Third, it is proposed that carbonado diamonds are a result of early impacts into crustal rocks. This is indicated by the rare and controversial occurrence of high pressure diamond polymorph, londsdaleite, frequently found in diamonds from meteorite impact sites, and by observation of planar deformation features, possibly indicating shock events. Finally, it is suggested that carbonados have an extraterrestrial origin, as indicated by a long term annealing based on observation of a zero-phonon line, attributed to paired nitrogen atoms in association with a vacancy.

Description: Vertical electron-density (N (sub e)) profiles, deduced from newly-available ISIS-II digital ionospheric topside-sounder data, are used to investigate the "polar-hole" region within the winter, nighttime polar cap ionosphere during solar minimum. The hole region is located around 0200 MLT near the poleward side of the auroral oval. Earlier investigations had revealed very low N (sub e) values in this region (down to 200/cu cm near 300 km). In the present study, such low N, values (approx. 100/cu cm) were only found near the ISIS (International Satellite for Ionospheric Study)-II altitude of 1400 km. The peak ionospheric concentration below the spacecraft remained fairly constant (approx. 10 (exp 5)/cu cm across the hole region but the altitude of the peak dropped dramatically. This peak dropped, surprisingly, to the vicinity of 100 km. These observations suggest that the earlier satellite in situ measurements, interpreted as deep holes in the ionospheric F-region concentration, could have been made during conditions of an extreme decrease in the altitude of the ionospheric N (sub e) peak. The observations, in combination with other data, indicate that the absence of an F-layer peak may be a frequent occurrence at high latitudes.

Description: We have examined the sea level height tide records at seven tide gauge sites in the region of southcentral Alaska that were affected by the 1964 Prince William Sound earthquake to determine the history of crustal uplift subsequent to the earthquake. There is considerable variation in the behavior depending on the location of the site relative to the 1964 rupture. At Seward, on the eastern side of the Kenai Peninsula we find a slow uplift that is consistent with elastic strain accumulation while at Seldovia and Nikiski on the western side of the Kenai we find a persistent rapid uplift of about 1 cm/yr that most likely represents a long term transient response to the earthquake, but which cannot be sustained over the expected recurrence interval for a great earthquake of several hundred years. Further to the southwest, at Kodiak, we find evidence that the rate of uplift, which is still several mm/yr, has slowed significantly over the past three and a half decades. To the east of the Kenai Peninsula we find subsidence at Cordova and an uncertain behavior at Valdez. At both of these sites there is a mathematically significant time-dependence to the uplift behavior, but the data confirming this time dependence are not as convincing as at Kodiak. At Anchorage, to the north there is little evidence of vertical motion since the earthquake. We compare these long term tide gauge records to recent GPS observations. In general there is reasonable consistency except at Anchorage and Cordova where the GPS measurement indicate somewhat more rapid uplift and subsidence, respectively.

Description: The geographic limits of the South Atlantic Anomaly, as defined by radiation damage, are compared to contours of geomagnetic total field intensity, as defined by the 1995 IGRF, for the present and recent past. The most likely secular variation of the geomagnetic field is estimated and used to extrapolate the geomagnetic field to the year 2100. This indicates that radiation damage to spacecraft and humans in space will likely increase and to cover a much larger geographic area.

Description: We have adopted the transport scenarios used in Part 1 to examine the sensitivity of stratospheric aircraft perturbations to transport changes in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric residence time and reduced the magnitude of the negative perturbation response in total ozone. Increasing the stratospheric K(sub yy) increased the residence time and enhanced the global scale negative total ozone response. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and results in a significantly weaker perturbation response, relative to the base case, throughout the stratosphere. We found a relatively minor model perturbation response sensitivity to the magnitude of K(sub yy) in the tropical stratosphere, and only a very small sensitivity to the magnitude of the horizontal mixing across the tropopause and to the strength of the mesospheric gravity wave drag and diffusion. These transport simulations also revealed a generally strong correlation between passive NO(sub y) accumulation and age of air throughout the stratosphere, such that faster transport rates resulted in a younger mean age and a smaller NO(y) mass accumulation. However, specific variations in K(sub yy) and mesospheric gravity wave strength exhibited very little NO(sub y)-age correlation in the lower stratosphere, similar to 3-D model simulations performed in the recent NASA "Models and Measurements" II analysis. The base model transport, which gives the most favorable overall comparison with inert tracer observations, simulated a global/annual mean total ozone response of -0.59%, with only a slightly larger response in the northern compared to the southern hemisphere. For transport scenarios which gave tracer simulations within some agreement with measurements, the annual/globally averaged total ozone response ranged from -0.45% to -0.70%. Our previous 1995 model exhibited overly fast transport rates, resulting in a global/annually averaged perturbation total ozone response of -0.25%, which is significantly weaker compared to the 1999 model. This illustrates how transport deficiencies can bias model simulations of stratospheric aircraft.

Description: According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Description: The subject of neotectonics, covering the structures and structural activity of the last 5 million years (i.e., post-Miocene) is a well-recognized field, including "active tectonics," focussed on the last 500,000 years in a 1986 National Research Council report of that title. However, there is a cartographic gap between tectonic maps, generally showing all features regardless of age, and maps of current seismic or volcanic activity. We have compiled a map intended to bridge this gap, using modern data bases and computer-aided cartographic techniques. The maps presented here are conceptually descended from an earlier map showing tectonic and volcanic activity of the last one million years. Drawn by hand with the National Geographic Society's 1975 "The Physical World" map as a base, the 1981 map in various revisions has been widely reproduced in textbooks and various technical publications. However, two decades of progress call for a completely new map that can take advantage of new knowledge and cartographic techniques. The digital tectonic activity map (DTM), presented in shaded relief (Fig. 1) and schematic (Fig. 2) versions, is the result. The DTM is intended to show tectonism and volcanism of the last one million years, a period long enough to be representative of global activity, but short enough that features such as fault scarps and volcanos are still geomorphically recognizable. Data Sources and Cartographic Methods The DTM is based on a wide range of sources, summarized in Table 1. The most important is the digital elevation model, used to construct a shaded relief map. The bathymetry is largely from satellite altimetry, specifically the marine gravity compilations by Smith and Sandwell (1996). The shaded relief map was designed to match the new National Geographic Society world physical map (1992), although drawn independently, from the digital elevation model. The Robinson Projection is used instead of the earlier Van der Grinten one. Although neither conformal nor equal-area, the Robinson Projection provides a reasonable compromise and retains useful detail at high latitudes.

Description: There are procedures and methods for verification of coding algebra and for validations of models and calculations that are in use in the aerospace computational fluid dynamics (CFD) community. These methods would be efficacious if used by the glacier dynamics modelling community. This paper is a presentation of some of those methods, and how they might be applied to uncertainty management supporting code verification and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modelling are discussed. After establishing sources of uncertainty and methods for code verification, the paper looks at a representative sampling of verification and validation efforts that are underway in the glacier modelling community, and establishes a context for these within overall solution quality assessment. Finally, an information architecture and interactive interface is introduced and advocated. This Integrated Cryospheric Exploration (ICE) Environment is proposed for exploring and managing sources of uncertainty in glacier modelling codes and methods, and for supporting scientific numerical exploration and verification. The details and functionality of this Environment are described based on modifications of a system already developed for CFD modelling and analysis.

Description: The behavior of mesospheric odd nitrogen species during and following relativistic and diffuse auroral precipitation events is simulated, Below 75 km nitric oxide is enhanced in proportion to the ion pair production function associated with the electron precipitation and the length of the event. Nitrogen dioxide and nitric acid are also enhanced. At 65 km the percentage of odd nitrogen for N is 0.1%, HNO3 is 1.6%, NO2 is 15%, and NO is 83.3%. Between 75 and 85 km NO is depleted during particle events due to the faster destruction of NO by N relative to the production of NO by N reacting with O2. Recovery of NO depends on transport from the lower thermosphere, where NO is produced in abundant amounts during particle events.

Description: Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Description: Magnetic remanence of crustal rocks can reside in three common rock-forming magnetic minerals: magnetite, pyrrhotite, and hematite. Thermoremanent magnetization (TRM) of magnetite and pyrrhotite is carried mostly by single domain (SD) grains. The TRM of hematite grains, however, is carried mostly by multidomain (NM) grains. This characteristic is illustrated by TRM acquisition curves for hematite of variable grainsizes. The transition between truly NM behavior and tendency towards SD behavior his been established between hematite grainsizes of 0. 1 and 0.05 mm. Coarse grainsize of lower crustal rocks and the large sensitivity of MD hematite grains to acquire TRM indicates that hematite could be a significant contributor to long-wavelength magnetic anomalies.

Description: We present an extension for the 2D (zonal mean) version of our Numerical Spectral Mode (NSM) that incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations and the semi-annual and quasi-biennial oscillations (SAO and QBO). Our earlier model reproduced the salient features of the mean zonal circulation in the middle atmosphere, including the QBO extension into the upper mesosphere inferred from UARS measurements. In the present model we incorporate also tropospheric heating to reproduce the upwelling at equatorial latitudes associated with the Brewer-Dobson circulation that affects significantly the dynamics of the stratosphere as Dunkerton had pointed out. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. The QBO period in the model is 30 months (mo), which is conducive to synchronize this oscillation with the seasonal cycle of solar forcing. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. Quadratic non-linearities generate interseasonal variations to produce a complicated pattern of variability associated with the QBO. The computed temperature amplitudes for the SAO and QBO are in substantial agreement with observations at equatorial and extratropical latitudes. At high latitudes, however, the observed QBO amplitudes are significantly larger, which may be a signature of propagating planetary waves not included in the present model. The assumption of hydrostatic equilibrium not being imposed, we find that the effects from the vertical Coriolis force associated with the equatorial oscillations are large for the vertical winds and significant for the temperature variations even outside the tropics but are relatively small for the zonal winds.

Description: This paper reviews observations of stratospheric temperatures that have been made over a period of several decades. Those observed temperatures have been used to assess variations and trends in stratospheric temperatures. A wide range of observation datasets have been used, comprising measurements by radiosonde (1940s to the present), satellite (1979 - present), lidar (1979 - present) and rocketsonde (periods varying with location, but most terminating by about the mid-1990s). In addition, trends have also been assessed from meteorological analyses, based on radiosonde and/or satellite data, and products based on assimilating observations into a general circulation model. Radiosonde and satellite data indicate a cooling trend of the annual-mean lower stratosphere since about 1980. Over the period 1979-1994, the trend is 0.6K/decade. For the period prior to 1980, the radiosonde data exhibit a substantially weaker long-term cooling trend. In the northern hemisphere, the cooling trend is about 0.75K/decade in the lower stratosphere, with a reduction in the cooling in mid-stratosphere (near 35 km), and increased cooling in the upper stratosphere (approximately 2 K per decade at 50 km). Model simulations indicate that the depletion of lower stratospheric ozone is the dominant factor in the observed lower stratospheric cooling. In the middle and upper stratosphere both the well-mixed greenhouse gases (such as CO) and ozone changes contribute in an important manner to the cooling.

Description: This paper presents a detailed characterization of seasonal and interannual variability in tropical tropospheric column ozone (TCO). TCO time series are derived from 20 years (1979-1998) of total ozone mapping spectrometer (TOMS) data using the convective cloud differential (CCD) method. Our study identifies three regions in the tropics with distinctly different zonal characteristics related to seasonal and interannual variability. These three regions are the eastern Pacific, Atlantic, and western Pacific. Results show that in both the eastern and western Pacific seasonal-cycle variability of northern hemisphere (NH) TCO exhibits maximum amount during NH spring whereas largest amount in southern hemisphere (SH) TCO occurs during SH spring. In the Atlantic, maximum TCO in both hemispheres occurs in SH spring. These seasonal cycles are shown to be comparable to seasonal cycles present in ground-based ozonesonde measurements. Interannual variability in the Atlantic region indicates a quasi-biennial oscillation (QBO) signal that is out of phase with the QBO present in stratospheric column ozone (SCO). This is consistent with high pollution and high concentrations of mid-to-upper tropospheric O3-producing precursors in this region. The out of phase relation suggests a UV modulation of tropospheric photochemistry caused by the QBO in stratospheric O3. During El Nino events there is anomalously low TCO in the eastern Pacific and high values in the western Pacific, indicating the effects of convectively-driven transport of low-value boundary layer O3 (reducing TCO) and O3 precursors including H2O and OH. A simplified technique is proposed to derive high-resolution maps of TCO in the tropics even in the absence of tropopause-level clouds. This promising approach requires only total ozone gridded measurements and utilizes the small variability observed in TCO near the dateline. This technique has an advantage compared to the CCD method because the latter requires high-resolution footprint measurements of both reflectivity and total ozone in the presence of tropopause-level cloud tops.

Description: Tropospheric column ozone (TCO) and stratospheric column ozone (SCO) gridded data in the tropics for 1979-present are now available from NASA Goddard Space Flight Center via either direct ftp, world-NN,ide-NN,eb, or electronic mail. This note provides a brief overview of the method used to derive the data set including validation and adjustments.

Description: In this study, we examine the sensitivity of long lived tracers to changes in the base transport components in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric mean age, and increased the rate of removal of material from the stratosphere. Increasing the stratospheric K(sub yy) increased the mean age due to the greater recycling of air parcels through the middle atmosphere, via the residual circulation, before returning to the troposphere. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and produces significantly younger ages throughout the stratosphere. Simulations with very small tropical stratospheric K(sub yy) decreased the globally averaged age of air by as much as 25% in the middle and upper stratosphere, and resulted in substantially weaker vertical age gradients above 20 km in the extratropics. We found only very small stratospheric tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and to the strength of the mesospheric gravity wave drag and diffusion used in the model. We also investigated the transport influence on chemically active tracers and found a strong age-tracer correlation, both in concentration and calculated lifetimes. The base model transport gives the most favorable overall comparison with a variety of inert tracer observations, and provides a significant improvement over our previous 1995 model transport. Moderate changes to the base transport were found to provide modest agreement with some of the measurements. Transport scenarios with residence times ranging from moderately shorter to slightly longer relative to the base case simulated N2O lifetimes that were within the observational estimates of Volk et al. [1997]. However, only scenarios with rather fast transport rates were comparable with the Volk et al. estimates of CFCl3 lifetimes. This is inconsistent with model-measurement comparisons of mean age in which the base model or slightly slower transport rates compared the most favorably with balloon SF6 data. For all comparisons shown, large transport changes away from the base case resulted in simulations that were outside the range of measurements, and in many cases, far outside this range.

Description: An intensive field campaign involving measurement of various aerosol physical, chemical, and radiative properties was conducted at Sde Boker in the Negev Desert of Israel, from 18 February to 15 March 1997. Nephelometer measurements gave average background scattering coefficient values of about 25 M/m at 550 nm wavelength, but strong dust events caused the value of this parameter to rise up to about 800 M/m Backscattering fractions did not depend on aerosol loading, and generally fell in the range of 0.1 to 0.25, comparable to values reported for marine and Arctic environments. Chemical analysis of the aerosol revealed that, in the coarse size range (2 - 10 micrometer equivalent aerodynamic diameter (EAD)), calcium (Ca) was by far the most abundant element followed by silicon (Si), both of which are indicators for mineral dust. In the fine size fraction (〈 2 micrometers EAD), sulfur (S) generally was the dominant element, except during high dust episodes when Ca and Si were again the most abundant. Furthermore, fine black carbon (BC) correlates with S, suggesting that they may have originated from the same sources or source regions. An indication of the short-term effect of aerosol loading on radiative forcing was provided by measurements of global and diffuse solar radiation, which showed that during high turbidity periods (strong dust events) almost all of the solar radiation reaching the area is scattered or absorbed.

Description: This presentation provides a brief overview of the scientific activities of the International GPS Service (IGS). This was an approved activity of the International Association of Geodesy (IAG) with official start of service on 1 Jan 1994. The mission of the IGS is "To provide a service to support geodetic and geophysical research activities, through GPS data and data products." The presentation explains the concept of the IGS working group, and pilot projects, and reviews the current working groups and pilot projects.

Description: The objective of this study is to calibrate a Navier-Stokes code for a high-lift Reference H configuration using structured grids. The outline of this presentation will first include a brief description of the grids used and the flow solver. Next the results will be presented in terms of convergence and resources used on the C-90. Predicted force and moment and surface pressure results are compared to experiment and off- and on-surface flow viz. is discussed.

Description: With the continuous progress in hardware and numerical schemes, Computational Unsteady Aerodynamics (CUA), that is, the application of Computational Fluid Dynamics (CFD) to unsteady flowfields, is slowly finding its way as a useful and reliable tool (turbulence and transition modeling permitting) in the aircraft, helicopter, engine and missile design and development process. Before a specific code may be used with confidence it is essential to validate its capability to describe the physics of the flow correctly, or at least to the level of approximation required, for which purpose a comparison with accurate experimental data is needed. Unsteady wind tunnel testing is difficult and expensive; two factors which dramatically limit the number of organizations with the capability and/or resources to perform it. Thus, unsteady experimental data is scarce, often classified and scattered in diverse documents. Additionally, access to the reports does not necessarily assure access to the data itself. The collaborative effort described in this paper was conceived with the aim of collecting into a single easily accessible document as much quality data as possible. The idea is not new. In the early 80's NATO's AGARD (Advisory Group for Aerospace Research & Development) Structures and Material Panel (SMP) produced AGARD Report No. 702 "Compendium of Unsteady Aerodynamic Measurements", which has found and continues to find extensive use within the CUA Community. In 1995 AGARD's Fluid Dynamics Panel (FDP) decided to update and expand the former database with new geometries and physical phenomena, and launched Working Group WG-22 on "Validation Data for Computational Unsteady Aerodynamic Codes". Shortly afterwards AGARD was reorganized as the RTO (Research and Technology Organization) and the WG was renamed as AVT (Applied Vehicle Technolology) WG-003. Contributions were received from AEDC, BAe, DLR, DERA, Glasgow University, IAR, NAL, NASA, NLR, and ONERA. The final publication with the results of the exercise is expected in the second part of 1999. The aim of the present paper is to announce and present the new database to the Aeroelasticity community. It is also intended to identify, together with one of the groups of end users it targets, deficiencies in the compendium that should be addressed by means of new wind tunnel tests or by obtaining access to additionally existing data.

Description: The HSR program moved into phase two with the selection of a new airplane configuration, the Technology Concept Airplane (TCA). The TCA was designed based on the experiences gained while investigating both the Reference H and the Arrow Wing configurations in different wind tunnels and CFD studies. Part of that investigation included performing extensive high Reynolds number testing on the Reference H configuration in the NTF to provide data for predicting full-scale flight performance, as well as developing techniques for testing these types of configurations in the NTF. With the selection of the TCA configuration, a smaller investigation was designed to examine whether or not the scaling characteristics of the TCA configuration are similar to those observed for the Reference H configuration. This presentation will include a description of the 2.2% Modified Reference H model used in this investigation (highlighting the similarities and the differences when compared to the TCA configuration), the testing objectives, and some preliminary findings that are relevant to the current high-lift system.

Description: The codes WINGDES2 and AERO2S were easy to obtain, and technical help was readily available. The codes have a long, well-documented history of successful optimizations of various aircraft configurations. The codes were easy to use, although specification of input data was time-consuming. The Run times were short, allowing the many runs necessary for the Suction Parameter matrix to be accomplished within a day or two. The results of the optimization appear to be reasonable.