ALBERT

Description: This report describes the transitional activities of the JPL Analysis Center.

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: 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: Early observations with ERS-1 SAR image data revealed a large ice stream in North East Greenland (Fahnestock 1993). The ice stream has a number of the characteristics of the more closely studied ice streams in Antarctica, including its large size and gross geometry. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. These features can be seen in the balance velocities for the ice stream (Joughin 1997) and its outlets. The ice stream is identifiable for more than 700 km, making it much longer than any other flow feature in Greenland. Our research goals are to gain a greater understanding of the ice flow in the northeast Greenland ice stream and its outlet glaciers in order to assess their impact on the past, present, and future mass balance of the ice sheet. We will accomplish these goals using a combination of remotely sensed data and ice sheet models. We are using satellite radar interferometry data to produce a complete maps of velocity and topography over the entire ice stream. We are in the process of developing methods to use these data in conjunction with existing ice sheet models similar to those that have been used to improve understanding of the mechanics of flow in Antarctic ice streams.

Description: This research is focusing on two related areas that are fundamental to the NASA PARCA (Program for Arctic Regional Climate Assessment) program. The primary research area is the determination of the amount, rate, and timing of accumulation at distributed sites in the dry snow zone of Greenland and evaluation of these results in light of accumulation modeling results. The secondary research area is the calibration of the isotope "thermometer" at these ice sheet sites as well as the determination of long-term temperature trends in Greenland.

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: This report summarizes the current and future plans of the Crustal Dynamics Data Information System (CDDIS) with respect to the International VLBI Service for Geodesy and Astrometry (IVS). Included are background information about the CDDIS, the computer architecture, staffing supporting the system, archive contents, and future plans for the CDDIS within the IVS.

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: 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: This year repeat elevation surveys in the southern half of Greenland were made using the Airborne Topographic Mapper (ATM). The intent of these surveys is to compare present elevations to those measured in 1993 and determine the magnitude and spatial distributions of thickening and thinning rates. In order to effectively interpret any observed changes, it is important to understand the processes that affect these changes. Moreover, because the surveys are made over a brief period (2-4 weeks) during the spring or summer, it is also important to understand the effects of seasonal and interannual elevation variability, in relation to the timing of these surveys. Toward that end we are examining data from weather stations along the coast of Greenland along with data from GC-Net automatic weather stations (AWS's) on the ice sheet. The objectives are to assess: a) the importance of the timing of the flights in relation to natural processes that affect surface heights, namely accumulation and melt, and b) the temperature characteristics of the region in the five years that separated the two sets of surveys (1993-1998), in relation to the past 19 years.

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 report gives an overall view of the CORE program at Goddard Space Flight Center (GSFC). It summarizes the different CORE sessions and gives information about the technical staff. The outlook summarizes the evolution of the different CORE programs.

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: Airborne laser-altimeter flight lines from 1993 over southern parts of the ice sheet were resurveyed with almost complete repeat coverage. In 1993 and 1994, NASA surveyed the entire Greenland ice sheet by airborne laser altimeter, obtaining surface-elevation profiles with root mean square (rms) accuracies of 10 cm or better (Krabill 1995) along flight lines that crossed all the major catchment basins. In 1998, the ten flight lines flown in 1993 in the south of Greenland were resurveyed with about 99% repeat coverage; flight lines in the north will be resurveyed in 1999. Additional flights in 1998 were over glaciers, identified by E. Rignot, where existing SAR data give information on ice motion.

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: Geosynchronous Synthetic Aperture Radar (GeoSAR) is a consortium project consisting of The Jet Propulsion Laboratory (JPL), Calgis (a small GIS company based in Fresno, CA) and the California Department of Conservation with funding provided by Defense Research Projects Agency (DARPA) started in November 1996. The two main objectives of the GeoSAR Program are: 1) To develop a state of the art dual frequency interferometric radar mapping instrument capable of mapping the true ground surface height beneath the vegetation canopy; and 2) To transition this mapping technology to a commercial company, Calgis. JPL, the technical lead, has the following program deliverables at program completion in November 1999 include radar design and radar hardware for X-band (3 cm) and P-band (83 cm) radars, processor software, hardware and documentation, and calibrated X-band radar.

Description: Data in the wavelength range 0.545 - 1.652 microns from the Moderate Resolution Imaging Spectroradiometer (MODIS), to be launched aboard the Earth Observing System (EOS) Terra in the fall of 1999, will be used to map daily global snow cover at 500m resolution. However, during darkness, or when the satellite's view of the surface is obscured by cloud, snow cover cannot be mapped using MODIS data. We show that during these conditions, it is possible to supplement the MODIS product by mapping the snow cover using passive microwave data from the Special Sensor Microwave Imager (SSM/I), albeit with much poorer resolution. For a 7-day time period in March 1999, a prototype MODIS snow-cover product was compared with a prototype MODIS-SSM/I product for the same area in the mid-western United States. The combined MODIS-SSM/I product mapped 9% more snow cover than the MODIS-only product.

Description: The Laser Vegetation Imaging Sensor (LVIS) is an airborne, scanning laser altimeter designed and developed at NASA's Goddard Space Flight Center. LVIS operates at altitudes up to 10 km above ground, and is capable of producing a data swath up to 1000 m wide nominally with 25 m wide footprints. The entire time history of the outgoing and return pulses is digitized, allowing unambiguous determination of range and return pulse structure. Combined with aircraft position and attitude knowledge, this instrument produces topographic maps with decimeter accuracy and vertical height and structure measurements of vegetation. The laser transmitter is a diode-pumped Nd:YAG oscillator producing 1064 nm, 10 nsec, 5 mJ pulses at repetition rates up to 500 Hz. LVIS has recently demonstrated its ability to determine topography (including sub-canopy) and vegetation height and structure on flight missions to various forested regions in the U.S. and Central America. The LVIS system is the airborne simulator for the Vegetation Canopy Lidar (VCL) mission (a NASA Earth remote sensing satellite due for launch in 2000), providing simulated data sets and a platform for instrument proof-of-concept studies. The topography maps and return waveforms produced by LVIS provide Earth scientists with a unique data set allowing studies of topography, hydrology, and vegetation with unmatched accuracy and coverage.

Description: The authors propose a grouped threshold method for scene identification in Advanced Very High Resolution Radiometer imagery that may contain clouds, fire, smoke, or snow. The philosophy of the approach is to build modules that contain groups of spectral threshold tests that are applied concurrently, not sequentially, to each pixel in an image. The purpose of each group of tests is to identify uniquely a specific class in the image, such as smoke. A strength of this approach is that insight into the limits used in the threshold tests may be gained through the use of radiative transfer theory. Methodology and examples are provided for two different scenes, one containing clouds, forest fires, and smoke; and the other containing clouds over snow in the central United States. For both scenes, a limited amount of supporting information is provided by surface observers.

Description: The upcoming generation of laser altimeters record the interaction of emitted laser radiation with terrestrial surfaces in the form of a digitized waveform. We model these laser altimeter return waveforms as the sum of the reflections from individual surfaces within laser footprints, accounting for instrument-specific properties. We compare over 1000 modeled and recorded waveform pairs using the Pearson correlation. We show that we reliably synthesize the vertical structure information for vegetation canopies contained in a medium-large diameter laser footprint from a high-resolution elevation data set.

Description: Radar data from vegetated land surfaces depend on many structural and compositional parameters describing the terrain. Because early, noninterferometric radar systems usually constituted an insufficient observation set from which to estimate parameters of the terrain, statistical regression techniques were used which incorporated some level of apriori knowledge or field measurements. With the advent of radar interferometry and polarimetric interferometry, potentially at multiple baselines, the observation set is now approaching that required to quantitatively estimate the parameters describing a vegetated land surface. Quantitative estimation entails formulating a physical scattering model relating the radar observations to the vegetation and surface parameters on which they depend. This paper describes the physics of candidate scattering models, and shows how the models determine the estimable parameter set. It also indicates the measurement accuracy of parameters such as vegetation height, height-to-base-of-live-crown, and surface topography with multibaseline polarimetric interferometry.

Description: In this paper, a multiscale approach is introduced to classify the Japanese Research Satellite-1 (JERS-1) mosaic image over the Central African rainforest. A series of texture maps are generated from the 100 m mosaic image at various scales. Using a quadtree model and relating classes at each scale by a Markovian relationship, the multiscale images are classified from course to finer scale. The results are verified at various scales and the evolution of classification is monitored by calculating the error at each stage.

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: This study extends a previous investigation on estimating surface soil moisture using the Special Sensor Microwave/Imager (SSM/I) over a grassland region. Although SSM/I is not optimal for soil moisture retrieval, it can under some conditions provide information. Rigorous analyses over land have been difficult due to the lack of good validation data sets. A scientific objective of the Southern Great Plains 1997 (SGP97) Hydrology Experiment was to investigate whether the retrieval algorithms for surface soil moisture developed at higher spatial resolution using truck-and aircraft-based passive microwave sensors can be extended to the coarser resolutions expected from satellite platform. With the data collected for the SGP97, the objective of this study is to compare the surface soil moisture estimated from the SSM/I data with those retrieved from the L-band Electronically Scanned Thinned Array Radiometer (ESTAR) data, the core sensor for the experiment, using the same retrieval algorithm. The results indicated that an error of estimate of 7.81% could be achieved with SSM/I data as contrasted to 2.82% with ESTAR data over three intensive sampling areas of different vegetation regimes. It confirms the results of previous study that SSM/I data can be used to retrieve surface soil moisture information at a regional scale under certain conditions.

Description: Two ascending European Space Agency (ESA) Earth Resources Satellites (ERS)-1/-2 tandem-mode, synthetic aperture radar (SAR) pairs are used to calculate the surface elevation of Hofsjokull, an ice cap in central Iceland. The motion component of the interferometric phase is calculated using the 30 arc-second resolution USGS GTOPO30 global digital elevation product and one of the ERS tandem pairs. The topography is then derived by subtracting the motion component from the other tandem pair. In order to assess the accuracy of the resultant digital elevation model (DEM), a geodetic airborne laser-altimetry swath is compared with the elevations derived from the interferometry. The DEM is also compared with elevations derived from a digitized topographic map of the ice cap from the University of Iceland Science Institute. Results show that low temporal correlation is a significant problem for the application of interferometry to small, low-elevation ice caps, even over a one-day repeat interval, and especially at the higher elevations. Results also show that an uncompensated error in the phase, ramping from northwest to southeast, present after tying the DEM to ground-control points, has resulted in a systematic error across the DEM.

Description: Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Description: The age of secondary forests in the Amazon will become more critical with respect to the estimation of biomass and carbon budgets as tropical forest conversion continues. Multitemporal Thematic Mapper data were used to develop land cover histories for a 33,000 Square kM area near Ariquemes, Rondonia over a 7 year period from 1989-1995. The age of the secondary forest, a surrogate for the amount of biomass (or carbon) stored above-ground, was found to be unimportant in terms of biomass budget error rates in a forested TM scene which had undergone a 20% conversion to nonforest/agricultural cover types. In such a situation, the 80% of the scene still covered by primary forest accounted for over 98% of the scene biomass. The difference between secondary forest biomass estimates developed with and without age information were inconsequential relative to the estimate of biomass for the entire scene. However, in futuristic scenarios where all of the primary forest has been converted to agriculture and secondary forest (55% and 42% respectively), the ability to age secondary forest becomes critical. Depending on biomass accumulation rate assumptions, scene biomass budget errors on the order of -10% to +30% are likely if the age of the secondary forests are not taken into account. Single-date TM imagery cannot be used to accurately age secondary forests into single-year classes. A neural network utilizing TM band 2 and three TM spectral-texture measures (bands 3 and 5) predicted secondary forest age over a range of 0-7 years with an RMSE of 1.59 years and an R(Squared) (sub actual vs predicted) = 0.37. A proposal is made, based on a literature review, to use satellite imagery to identify general secondary forest age groups which, within group, exhibit relatively constant biomass accumulation rates.

Description: We obtain valuable information on the angular and seasonal variability of surface reflectance using a hand-held spectrometer from a light aircraft. The data is used to test a procedure that allows us to estimate visible surface reflectance from the longer wavelength 2.1 micrometer channel (mid-IR). Estimating or avoiding surface reflectance in the visible is a vital first step in most algorithms that retrieve aerosol optical thickness over land targets. The data indicate that specular reflection found when viewing targets from the forward direction can severely corrupt the relationships between the visible and 2.1 micrometer reflectance that were derived from nadir data. There is a month by month variation in the ratios between the visible and the mid-IR, weakly correlated to the Normalized Difference Vegetation Index (NDVI). If specular reflection is not avoided, the errors resulting from estimating surface reflectance from the mid-IR exceed the acceptable limit of DELTA-rho approximately 0.01 in roughly 40% of the cases, using the current algorithm. This is reduced to 25% of the cases if specular reflection is avoided. An alternative method that uses path radiance rather than explicitly estimating visible surface reflectance results in similar errors. The two methods have different strengths and weaknesses that require further study.

Description: Snow-cover maps generated from surface data are based on direct measurements, however they are prone to interpolation errors where climate stations are sparsely distributed. Snow cover is clearly discernable using satellite-attained optical data because of the high albedo of snow, yet the surface is often obscured by cloud cover. Passive microwave (PM) data is unaffected by clouds, however, the snow-cover signature is significantly affected by melting snow and the microwaves may be transparent to thin snow (less than 3cm). Both optical and microwave sensors have problems discerning snow beneath forest canopies. This paper describes a method that combines ground and satellite data to produce a Multiple-Dataset Snow-Cover Product (MDSCP). Comparisons with current snow-cover products show that the MDSCP draws together the advantages of each of its component products while minimizing their potential errors. Improved estimates of the snow-covered area are derived through the addition of two snow-cover classes ("thin or patchy" and "high elevation" snow cover) and from the analysis of the climate station data within each class. The compatibility of this method for use with Moderate Resolution Imaging Spectroradiometer (MODIS) data, which will be available in 2000, is also discussed. With the assimilation of these data, the resolution of the MDSCP would be improved both spatially and temporally and the analysis would become completely automated.

Description: This paper investigates the sensitivity of potential evapotranspiration to input meteorological variables, viz- surface air temperature and surface vapor pressure. The sensitivity studies have been carried out for a wide range of land surface variables such as wind speed, leaf area index and surface temperatures. Errors in the surface air temperature and surface vapor pressure result in errors of different signs in the computed potential evapotranspiration. This result has implications for use of estimated values from satellite data or analysis of surface air temperature and surface vapor pressure in large scale hydrological modeling. The comparison of cumulative potential evapotranspiration estimates using ground observations and satellite observations over Manhattan, Kansas for a period of several months shows very little difference between the two. The cumulative differences between the ground based and satellite based estimates of potential evapotranspiration amounted to less that 20mm over a 18 month period and a percentage difference of 15%. The use of satellite estimates of surface skin temperature in hydrological modeling to update the soil moisture using a physical adjustment concept is studied in detail including the extent of changes in soil moisture resulting from the assimilation of surface skin temperature. The soil moisture of the surface layer is adjusted by 0.9mm over a 10 day period as a result of a 3K difference between the predicted and the observed surface temperature. This is a considerable amount given the fact that the top layer can hold only 5mm of water.

Description: The effective use of satellite observations of the land surface is limited by the lack of high spatial resolution ground data sets for validation of satellite products. Recent large scale field experiments include FIFE, HAPEX-Sahel and BOREAS which provide us with data sets that have large spatial coverage and long time coverage. It is the objective of this paper to characterize the difference between the satellite estimates and the ground observations. This study and others along similar lines will help us in utilization of satellite retrieved data in large scale modeling studies.

Description: Forest ecosystem dynamics modeling, remote sensing data analysis, and a geographical information system (GIS) were used together to determine the possible growth and development of a northern forest in Maine, USA. Field measurements and airborne synthetic aperture radar (SAR) data were used to produce maps of forest cover type and above ground biomass. These forest attribute maps, along with a conventional soils map, were used to identify the initial conditions for forest ecosystem model simulations. Using this information along with ecosystem model results enabled the development of predictive maps of forest development. The results obtained were consistent with observed forest conditions and expected successional trajectories. The study demonstrated that ecosystem models might be used in a spatial context when parameterized and used with georeferenced data sets.

Description: In anticipation of the launch of the Earth Observing System (EOS) Terra, and the PM-1 spacecraft in 1999 and 2000, respectively, efforts are ongoing to determine errors of satellite-derived snow-cover maps. EOS Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer-E (AMSR-E) snow-cover products will be produced. For this study we compare snow maps covering the same study area acquired from different sensors using different snow- mapping algorithms. Four locations are studied: 1) southern Saskatchewan; 2) a part of New England (New Hampshire, Vermont and Massachusetts) and eastern New York; 3) central Idaho and western Montana; and 4) parts of North and South Dakota. Snow maps were produced using a prototype MODIS snow-mapping algorithm used on Landsat Thematic Mapper (TM) scenes of each study area at 30-m and when the TM data were degraded to 1 -km resolution. National Operational Hydrologic Remote Sensing Center (NOHRSC) 1 -km resolution snow maps were also used, as were snow maps derived from 1/2 deg. x 1/2 deg. resolution Special Sensor Microwave Imager (SSM/1) data. A land-cover map derived from the International Geosphere-Biosphere Program (IGBP) land-cover map of North America was also registered to the scenes. The TM, NOHRSC and SSM/I snow maps, and land-cover maps were compared digitally. In most cases, TM-derived maps show less snow cover than the NOHRSC and SSM/I maps because areas of incomplete snow cover in forests (e.g., tree canopies, branches and trunks) are seen in the TM data, but not in the coarser-resolution maps. The snow maps generally agree with respect to the spatial variability of the snow cover. The 30-m resolution TM data provide the most accurate snow maps, and are thus used as the baseline for comparison with the other maps. Comparisons show that the percent change in amount of snow cover relative to the 3 0-m resolution TM maps is lowest using the TM I -km resolution maps, ranging from 0 to 40%. The highest percent change (less than 100%) is found in the New England study area, probably due to the presence of patchy snow cover. A scene with patchy snow cover is more difficult to map accurately than is a scene with a well-defined snowline such as is found on the North and South Dakota scene where the percent change ranged from 0 to 40%. There are also some important differences in the amount of snow mapped using the two different SSM/I algorithms because they utilize different channels.

Description: Landsat-7 is scheduled for launch on April 15 from the Western Test Range at Vandenberg Air Force Base, Calif., on a Delta-H expendable launch vehicle. The Landsat 7 satellite consists of a spacecraft bus being provided by Lockheed Martin Missiles and Space (Valley Forge, Pa.) and the Enhanced Thematic Mapper Plus instrument built by Raytheon (formerly Hughes) Santa Barbara Remote Sensing (Santa Barbara, Calif.). The instrument on board Landsat 7 is the Enhanced Thematic Mapper Plus (ETM+). ETM+ improves upon the previous Thematic Mapper (TM) instruments on Landsat's 4 and 5 (Fig. la and lb). It includes the previous 7 spectral bands measuring reflected solar radiation and emitted thermal emissions but, in addition, includes a new 15 in panchromatic (visible-near infrared) band. The spatial resolution of the thermal infrared band has also been improved to 60 m. Both the radiometric precision and accuracy of the sensor are also improved from the previous TM sensors. After being launched into a sun-synchronous polar orbit, the satellite will use on-board propulsion to adjust its orbit to a circular altitude of 438 miles (705 kilometers) crossing the equator at approximately 10 a.m. on its southward track. This orbit will place Landsat 7 along the same ground track as previous Landsat satellites. The orbit will be maintained with periodic adjustments for the life of the mission. A three-axis attitude control subsystem will stabilize the satellite and keep the instrument pointed toward the Earth to within 0.05 degrees. Later this year, plans call for the NASA Earth Observation System (EOS) Terra (AM-1) observatory and the experimental EO-1 mission to closely follow Landsat-7's orbit to support synergistic research and applications from this new suite of terrestrial sensor systems. Landsat is the United States' oldest land-surface observation satellite system, with satellites continuously operating since 1972. Although the program has scored numerous successes in scientific and resource-management applications, Landsat has had a tumultuous history of management and funding changes over its nearly 27-year history. Landsat-7 marks a new direction in the program to reduce the cost of data and increase systematic global coverage for use in global change research as well as commercial and regional applications. With the passage of the Land Remote Sensing Policy Act in 1992, oversight of the Landsat program began to shift from the commercial sector to the federal government. NASA integrated Landsat-7 into its EOS science program in 1994. Landsat-7 is managed and operated jointly by NASA and U.S. Geological Survey (USGS). As a result, the costs of acquiring observations from

Description: We describe a model to calculate brightness temperature and surface energy balance for a forest canopy system. The model is an extension of an earlier vegetation only model by inclusion of a simple soil layer. The root mean square error in brightness temperature for a dense forest canopy was 2.5 C. Surface energy balance predictions were also in good agreement. The corresponding root mean square errors for net radiation, latent, and sensible heat were 38.9, 30.7, and 41.4 W/sq m respectively.

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.

Description: Full-scale flight performance predictions can be made using CFD or a combination of CFD and analytical skin-friction predictions. However, no matter what method is used to obtain full-scale flight performance predictions knowledge of the boundary layer state is critical. The implementation of CFD codes solving the Navier-Stokes equations to obtain these predictions is still a time consuming, expensive process. In addition, to ultimately obtain accurate performance predictions the transition location must be fixed in the CFD model. An example, using the M2.4-7A geometry, of the change in Navier-Stokes solution with changes in transition and in turbulence model will be shown. Oil flow visualization using the M2.4-7A 4.0% scale model in the 14'x22' wind tunnel shows that fixing transition at 10% x/c in the CFD model best captures the flow physics of the wing flow field. A less costly method of obtaining full-scale performance predictions is the use of non-linear Euler codes or linear CFD codes, such as panel methods, combined with analytical skin-friction predictions. Again, knowledge of the boundary layer state is critical to the accurate determination of full-scale flight performance. Boundary layer transition detection has been performed at 0.3 and 0.9 Mach numbers over an extensive Reynolds number range using the 2.2% scale Reference H model in the NTF. A temperature sensitive paint system was used to determine the boundary layer state for these conditions. Data was obtained for three configurations: the baseline, undeflected flaps configuration; the transonic cruise configuration; and, the high-lift configuration. It was determined that at low Reynolds number conditions, in the 8 to 10 million Reynolds number range, the baseline configuration has extensive regions of laminar flow, in fact significantly more than analytical skin-friction methods predict. This configuration is fully turbulent at about 30 million Reynolds number for both 0.3 and 0.9, Mach numbers. Both the transonic cruise and the high-lift configurations were fully turbulent aft of the leading-edge flap hingeline at all Reynolds numbers.

Description: The Transonic Nozzle Boattail Drag Study was initiated in 1995 to develop an understanding of how external nozzle transonic aerodynamics effect airplane performance and how strongly those effects are dependent on nozzle configuration (2D vs. axisymmetric). MDC analyzed the axisymmetric nozzle. Boeing subcontracted Northrop-Grumman to analyze the 2D nozzle. AU participants analyzed the AGARD nozzle as a check-out and validation case. Once the codes were checked out and the gridding resolution necessary for modeling the separated flow in this region determined, the analysis moved to the installed wing/body/nacelle/diverter cases. The boat tail drag validation case was the AGARD B.4 rectangular nozzle. This test case offered both test data and previous CFD analyses for comparison. Results were obtained for test cases B.4.1 (M=0.6) and B.4.2 (M=0.938) and compared very well with the experimental data. Once the validation was complete a CFD grid was constructed for the full Ref. H configuration (wing/body/nacelle/diverter) using a combination of patched and overlapped (Chimera) grids. This was done to ensure that the grid topologies and density would be adequate for the full model. The use of overlapped grids allowed the same grids from the full configuration model to be used for the wing/body alone cases, thus eliminating the risk of grid differences affecting the determination of the installation effects. Once the full configuration model was run and deemed to be suitable the nacelle/diverter grids were removed and the wing/body analysis performed. Reference H wing/body results were completed for M=0.9 (a=0.0, 2.0, 4.0, 6.0 and 8.0), M=1.1 (a=4.0 and 6.0) and M=2.4 (a=0.0, 2.0, 4.4, 6.0 and 8.0). Comparisons of the M=0.9 and M=2.4 cases were made with available wind tunnel data and overall comparisons were good. The axi-inlet/2D nozzle nacelle was analyzed isolated. The isolated nacelle data coupled with the wing/body result enabled the interference effects of the installed nacelles to be determined. Isolated nacelle mm were made at M=0.9 and M=1.1 for both the supersonic and transonic nozzle settings. AU of the isolated nacelle cases were run at alpha=0. Full configuration runs were to be made at Mach numbers of 0.9, 1.1, and 2.4 (the same as the wing/body and isolated nacelles). Both the isolated nacelles and installed nacelles were run with inlet conditions designed to give zero spillage. This was to be done in order to isolate the boattail effects as much as possible. Full configuration runs with the supersonic nozzles were completed for M=0.9 and 1.1 at a=4.0 and 6.0 (4 runs total) and with the transonic nozzles at M=0.9 and 1.1 at a=2.0, 4.0 and 6.0 (6 runs total). Drag breakdowns were completed for the M=0.9 and M= 1.1 showing favorable interference drag for both cases.

Description: This paper describes an Overview of Technology Integration. The topics include: 1) Technology Concept Airplane Description; 2) LCAP Overview; and 3) ACE (Aeroelastic Concept Engineering) Overview.

Description: The SLFC (Supersonic Laminar Flow Control) mission is to develop and validate technologies and perform the SLFC aerodynamic design for the HSCT with an assessment of the net benefits and risks.

Description: This paper describes HSR (High Speed Research) Aerodynamic Performance Status and Challenges. The topics include: 1) Aero impact on HSR; 2) Goals and Targets; 3) Progress and Status; and 4) Remaining Challenges. This paper is presented in viewgraph form.