ALBERT

Beschreibung: The test capabilities of the Stability Wind Tunnel of the Virginia Polytechnic Institute and State University are described, and calibrations for curved and rolling flow techniques are given. Oscillatory snaking tests to determine pure yawing derivatives are considered. Representative aerodynamic data obtained for a current fighter configuration using the curved and rolling flow techniques are presented. The application of dynamic derivatives obtained in such tests to the analysis of airplane motions in general, and to high angle of attack flight conditions in particular, is discussed.

Beschreibung: Solar variability is examined in search of implications for global change. The topics covered include the following: solar variation modification of global surface temperature; the significance of solar variability with respect to future climate change; and methods of reducing the uncertainty of the potential amplitude of solar variability on longer time scales.

Beschreibung: The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Beschreibung: The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Beschreibung: A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Beschreibung: We determined the infrared optical constants of nitric acid trihydrate, nitric acid dihydrate, nitric acid monohydrate, and solid amorphous nitric acid solutions which crystallize to form these hydrates. We have also found the infrared optical constants of H2O ice. We measured the transmission of infrared light throught thin films of varying thickness over the frequency range from about 7000 to 500/cm at temperatures below 200 K. We developed a theory for the transmission of light through a substrate that has thin films on both sides. We used an iterative Kramers-Kronig technique to determine the optical constants which gave the best match between measured transmission spectra and those calculated for a variety of films of different thickness. These optical constants should be useful for calculations of the infrared spectrum of polar stratospheric clouds.

Beschreibung: Meteorological conditions, extremely conducive to fire development and spread in the spring of 1987, resulted in forest fires burning over extremely large areas in the boreal forest zone in northeastern China and the southeastern region of Siberia. The great China fire, one of the largest and most destructive forest fires in recent history, occurred during this period in the Heilongjiang Province of China. Satellite imagery is used to examine the development and areal distribution of 1987 forest fires in this region. Overall trace gas emissions to the atmosphere from these fires are determined using a satellite-derived estimate of area burned in combination with fuel consumption figures and carbon emission ratios for boreal forest fires.

Beschreibung: In November to December 1991, a substantial number of remote sensors and in situ instruments were operated together in Coffeyville, Kansas, during the climate experiment First ISCCP Regional Experiment Phase 2 (FIRE 2). Includede in the suite of instruments were (1) the NOAA Environmental Technology Laboratory (ETL) three-channel microwave radiometer, (2) the NASA GSFC Raman lidar, (3) ETL radio acoustic sounding system (RASS), and (4) frequent, research-quality radiosondes. The Raman lidar operated only at night and the focus of this portion of the experiment concentrated on clear conditions. The lidar data, together with frequent radiosondes and measurements of temperature profiles (every 15 min) by RASS allowed profiles of temperature and absolute humidity to be estimated every minute. We compared 20 min measurements of brightness temperature (T(sub b) with calculations of T(sub b) that were based on the Liebe and Layton (1987) and Liebe et al. (1993) microwave propagation models, as well as the Waters (1976) model. The comparisons showed the best agreement at 20.6 GHz with the Waters model, with the Liebe et al. (1993) model being best at 31.65 GHz. The results at 90 GHz gave about equal success with the Liebe and Layton (1987) and Liebe et al. (1993) models. Comparisons of precipitable water vapor derived independently from the two instruments also showed excellent agreement, even for averages as short as 2 min. The rms difference between Raman and radiometric determinations of precipitable water vapor was 0.03 cm which is roughly 2%. The experiments clearly demonstrate the potential of simultaneous operation of radiometers and Raman lidars for fundamental physical studies of water vapor.

Beschreibung: This paper describes the validation of the International Satellite Cloud Climatology project (ISCCP)cloud detections by verifying the accuracy of the infrared clear-sky radiances. Comparison of retrieved surface temperatures to other measurements shows that bias errors are less than 2 K and random errors are about 2 K for sea surface (monthly means at 280-km scales) and that bias errors are less than 2 K and random errors are about 4 K for land surfaces (3 hourly at 280-km scales). Bias errors over a few persistently cloudy locations are sometimes -(2-4) K and over winter sea ice may be about +2 K. Surface reflectances are confirmed to be within 3% of other measurements and models for ocean, except for sun glint geometries, and to be within 3%-5% for land surfaces. Sufficiently accurate validation data are not avilable for visible reflectances of sea ice and snow-covered land, but some tests of specific cases suggest that errors are approximately 10%. These errors in clear-sky radiances suggest uncertainties in the ISCCP cloud detections of about 10% with a small (3%-6%) negative bias over land. Some specific regions exhibit both larger rms uncertainties and somewhat larger biases in cloud amount approaching 10%. ISCCP cloud detections are more in error over the polar regions than anywhere else. Based on comparisons with an anlysis of radiances measured at other wavelengths, the ISCCP analysis appears to miss 15%-25% of the clouds in summer but only 5%-10% of the winter clouds.

Beschreibung: Water vapor concentrations obtained by the Stratospheric Aerosol and Gas Experiment 2 (SAGE 2) and collocated temperatures provided by the National Meteorological Center (NMC) from 1986 to 1990 are used to deduce seasonally and zonally averaged acidity, density, and refractive index of stratospheric aerosols. It is found that the weight percentage of sulfuric acid in the aerosols increases from about 60 just above the tropopause to about 86 at 35 km. The density increases from about 1.55 to 1.85 g/cu cm between the same altitude limits. Some seasonal variations of composition and density are evident at high latitudes. The refractive indices at 1.02, 0.694, and 0.532 micrometers increase, respectively, from about 1.425, 1.430, and 1.435 just above the tropopause to about 1.445, 1.455, and 1.458 at altitudes above 27 km, depending on the season and latitude. The aerosol properties presented can be used in models to study the effectiveness of heterogeneous chemistry, the mass loading of stratospheric aerosols, and the extinction and backscatter of aerosols at different wavelengths. Computed aerosol surface areas, rate coefficients for the heterogeneous reaction ClONO2 + H2O yields HOCl + HNO3 and aerosol mass concentrations before and after the Pinatubo eruption in June 1991 are shown as sample applications.

Beschreibung: The global distribution of ozone during the October 31, 1978, to May 6, 1993, observing lifetime of the Nimbus 7 total ozone mapping spectrometer (TOMS) is described, with emphasis on the low ozone amounts observed during 1992 and 1993. Ozone amount time series are extended beyong May 6, 1993, to the end of July 1993 using preliminary Meteor 3 TOMS data. Time series for zonally averaged ozone amounts show that there has not been a significant shift in the seasonal patterns of ozone maxima and minima caused by the Mount Pinatubo eruption or by the onset of very low ozone values during 1992 and 1993. There has been a relatively slow, nearly linear decrease in the amount of ozone over the entire globe from 1979 to the end of 1991, with rates ranging from no change at the equator to a 4 - 6% decrease per decade at midlatitudes and a 10 - 12% decrease per decade at higher latitudes. After the eruption of Mount Pinatubo during June 1991, the ozone amount decreased in the equatorial latitudes (10 deg S to 10 deg N) for about 6 months (-10 Dobson units (DU) between 0 deg and 10 deg S and -3 DU between 0 deg and 10 deg N). During 1992 and continuing into 1993, the rate of ozone decrease deviated from the previously linear trend with the onset of changes that were large in comparison with the historical range of ozone values from 1979 to 1991. The first of the large decreases in ozone amount occurred earlier, in February 1990 to May 1990, at 50 deg - 70 deg N. At high northern latitudes, the 1993 decreased ozone amounts were about 12.5% below the envelope of historical values; at midlatitudes they were about 7% lower; and at low latitudes they were about 4% lower. Area-weighted averages in the northern and southern hemispheres show that most of the 1992 - 1993 ozone losses have occurred in the northern hemisphere. The 1993 global average (70 deg S to 70 deg N) ozone amount is 3% below the 1979 to 1991 minimum, 5% below the historical envelope in the northern hemisphere, and near the lower boundary of the historical envelope in the southern hemisphere. In the 70 deg - 60 deg S latitude band, the ozone losses between 1979 and 1993 have reduced the annual minimum amount to values below those seen in the equatorial regions.

Beschreibung: A new 8-year global cloud climatology has been produced by the International Satellite Cloud Climatology Project (ISCCP) that provides information every 3 h at 280-km spatial resolution covering the period from July 1983 through June 1991. If cloud detection errors and differences in area sampling are neglected, individual ISCCP cloud amounts agree with individual surface observations to within 15% rms with biases of only a few percent. When measurements of small-scale, broken clouds are isolated in the comparison, the rms differences between satellite and surface cloud amounts are about 25%, similar to the rms difference between ISCCP and Landsat determinations of cloud amount. For broken clouds, the average ISCCP cloud amounts are about 5% smaller than estimated by surface observers (difference between earth cover and sky cover), but about 5% larger than estimated from very high spatial resolution satellite observations (overestimate due to low spatial resolution offset by underestimate due to finite radiance thresholds). Detection errors caused by errors in the ISCCP clear-sky radiances or incorrect radiance threshold magnitudes, are the dominant source of error in monthly average cloud amounts. The ISCCP cloud amounts appear to be too low over land by about 10%, somewhat less in summer and somewhat more in winter, and about right (maybe slightly low) over oceans. In polar regions, ISCCP cloud amounts are probably too low by about 15%-25% in summer and 5%-10% in winter. Comparison of the ISCCP climatology to three other cloud climatologies shows excellent agreement in the geographic distribution and seasonal variation of cloud amounts; there is little agreement about day/night contrasts in cloud amount. Notable results from ISCCP are that the gobal annual mean cloud amount is about 63%, being about 23% higher over oceans than over land, that it varies by less than 1% rms from month to month, and that it has varied by about 4% on a time scale approximately equal to 2-4 years. The magnitude of interannual variations of local (280-km scale) monthly mean cloud amounts is about 7%-9%. Longitudinal contrasts in cloud amount are just as large as latitudinal contrasts. The largest seasonal variation of cloud amount occurs in the tropics, being larger in summer than in winter; the seasonal variation in the middle latitudes has the opposite phase. Polar regions may have little seasonable variability in cloud amount. The ISCCP results show slightly more nightime than daytime cloud amount over oceans and more daytime than nightime cloud amount over land.

Beschreibung: This paper, the first of three, describes the cloud detection part of the International Satellite Cloud Climatology Project (ISCCP) analysis. Key features of the cloud detection alogrithm are (1) use of space and time radiance variation tests over several different space and time domains to account for the global variety of cloudy and clear characteristics, (2) estimation of clear radiance values for every time and place, and, (3) use of radiance thresholds that vary with the type of surface and climate regime. Design of the detection algorithm was supported by global, multiyear surveys of the statistical behavior of satellite-measured infrared and visible radiances to determine those characteristics that differentiate cloudy and clear scenes and how these characteristics vary among climate regimes. A summary of these statistical results is presented to illustrate how the cloud detection method works in a variety of circumstances. The sensitivity of the results to changing test parameter values is determined to provide a first estimate of the uncertainty of ISCCP cloud amounts. These test results (which exclude polar regions) suggest detection uncertainties of about 10% with possible negative biases of 5% (especially at night).

Beschreibung: The Chicxulub impact crater in Mexico is the site of the impact purported to have caused mass extinctions at the Cretaceous/Tertiary (K/T) boundary. 2-D hydrocode modeling of the impact, coupled with studies of the impact site geology, indiate that between 0.4 and 7.0 x 10(exp 17) g of sulfur were vaporized by the impact into anhydrite target rocks. A small portion of the sulfur was released as SO3 or SO4, which converted rapidly into H2SO4 aerosol and fell as acid rain. A radiative transfer model, coupled with a model of coagulation indicates that the aerosol prolonged the initial blackout period caused by impact dust only if the aerosol contained impurities. A larger portion of sulfur was released as SO2, which converted to aerosol slowly, due to the rate-limiting oxidation of SO2. Our radiative transfer calculations, combined with rates of acid production, coagulation, and diffusion indicate that solar transmission was reduced to 10-20% of normal for a period of 8-13 yr. This reduction produced a climate forcing (cooling) of -300 W/sq.m, which far exceeded the +8 W/sq.m greenhouse warming, caused by the CO2 released through the vaporization of carbonates, and therefore produced a decade of freezing and near-freezing temperatures. Several decades of moderate warming followed the decade of severe cooling due to the long residence time of CO2. The prolonged impact winter may have been a major cause of the K/T extinctions.

Beschreibung: Rotor noise prediction codes predict the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predicitions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and cyclic flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the chordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of Central Processing Unit (CPU) time necessary for the various approximations.

Beschreibung: Clear-sky albedos and outgoing longwave radiation (OLR) determined from Earth Radiation Budget Experiment (ERBE) scanners on board the earth radiation budget satellite and NOAA-9 spacecraft were analyzed for three target sites for the months February 1985-January 1987. The targets were oceans, deserts, and a multiscene site covering half the earth's surface. Year-to-year ratios of the monthly albedos and OLR were within the 0.98-1.02 range with a standard error of about 1%. The data indicate that ERBE scanner measurements were stable to within a few tenths of a percent for the two-year periods.

Beschreibung: A set of system simulations has been performed to evaluate candidate scanner designs for an Earth Radiation Budget Instrument (ERBI) for the Earth Observing System (EOS) of the late 1990s. Five different instruments are considered: (1) the Active Cavity Array (ACA), (2) the Clouds and Earth's Radiant Energy System-Instrument (CERES-I), (3) the Conically Scanning Radiometer (CSR), (4) the Earth Radiation Budget Experiment Cross-Track Scanner (ERBE), and (5) the Nimbus-7 Biaxial Scanner (N7). Errors in instantaneous, top-of-the-atmosphere (TOA) satellite flux estimates are assumed to arise from two measurement problems: the sampling of space over a given geographic domain, and sampling in angle about a given spatial location. When angular sampling errors vanish due to the application of correct angular dependence models (ADMs) during inversion, the accuracy of each scanner design is determined by the instrument's ability to map the TOA radiance field in a uniform manner. In this regard, the instruments containing a cross-track scanning component (CERES-I and ERBE) do best. As errors in ADMs are encountered, cross-track instruments incur angular sampling errors more rapidly than biaxial instruments (N7, ACA, and CSR) and eventually overtake the biaxial designs in their total error amounts. A latitude bias (north-south error gradient) in the ADM error of cross-track instruments also exists. This would be objectionable when ADM errors are systematic over large areas of the globe. For instantaneous errors, however, cross-track scanners outperform biaxial or conical scanners for 2.5 deg latitude x 2.5 deg longitude target areas, providing that the ADM error is less than or equal to 30%. A key issue is the amount of systematic ADM error (departures from the mean models) that is present at the 2.5 deg resolution of the ERBE target areas. If this error is less than 30%, then the CERES-I, ERBE, and CSR, in order of increasing error, provide the most accurate instantaneous flux estimates, within 2-3 W/sq m of each other in reflected shortwave flux. The magnitude of this error is near the 10 W/sq m accuracy requirement of the user community. Longwave flux errors have been found to have the same space and time characteristics as errors in shortwave radiation, but only about 25% as large.

Beschreibung: This paper briefly reviews several single-frequency rain profiling methods for an airborne or spaceborne radar. The authors describe the different methods from a unified point of view starting from the basic differential equation. This facilitates the comparisons between the methods and also provides a better understanding of the physical and mathematical basis of the methods. The application of several methods to airborne radar data taken during the Convective and Precipitation/Electrification Experiment is shown. Finally, the authors consider a hybrid method that provides a smooth transition between the Hitschfeld-Bordan method, which performs well at low attenuations, and the surface reference method, for which the relative error decreases with increasing path attenuation.

Beschreibung: An analysis of a spontaneous sudden stratospheric warming that occurred during a 2-year integration of the Langley Research Center (LaRC) Atmospheric Simulation Model is presented. The simulated warming resembles observed 'wave 1' warmings in the Northern Hemisphere stratosphere and provides an opportunity to investigate the radiative and dynamical processes occurring during the warming event. Isentropic analysis of potential vorticity sources and sinks indicates that dynamically induced departures from radiative equilibrium play an important role in the warming event. Enhanced radiative cooling associated with a series of upper stratospheric warm pools leads to radiative dampening within the polar vortex. Within the 'surf zone' large-scale radiative cooling leads to diabatic advection of high potential vorticity air from aloft. Lagrangian area diagnostics of the simulated warming agree well with Limb Infrared Monitor of the Stratosphere (LIMS) analyses. Dynamical mixing is shown to account for the majority of the decrease in the size of the polar vortex during the simulated warming. An investigation of the nonlinear deformation of material lines that are initially coincident with diagnosed potential vorticity isopleths is conducted to clarify the relationship between the Lagrangian area diagnostics and potential vorticity advection during wave breaking events.

Beschreibung: Calculations of the quasi biennial oscillation (QBO) signal in Stratospheric Aerosol and Gas Experiment (SAGE) II O3 and NO2 data between 1984 and 1991 are presented and have been investigated by using a two-dimensional model. The isentropic 2D model is a fully interactive radiative-dynamical-chemical model in which the eddy fluxes of chemical species are calculated in a consistent manner. The QBO in the model has been forced by relaxing the equatorial zonal wind toward the observations at Singapore allowing the comparison of the model with observations from specific years. The model reproduces the observed vertical structure of the equatorial ozone anomaly with the well-known transition from dynamical to photochemical control at around 28km. The model also reproduces the observed vertical structure of the SAGE II observed NO2 anomaly. The model studies have shown that it is the QBO modulation of NO2 which the main cause of QBO signal in O3 above 30 km. The model also reproduces the observed latitudinal structure of the QBO signals in O3 and NO2. Due to the differing horizontal distribution of O3 and NO(y) the ozone signal shows a distinct phase change in the subtropics whereas the NO2 anomaly gives a broader signal.

Beschreibung: A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Beschreibung: A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.

Beschreibung: Two stratospheric warmings during February and March 1993 are described using United Kingdom Meteorological Office (UKMO) analyses, calculated potential vorticity (PV) and diabetic heating, and N2O observed by the Cryogenic Limb Array Etalon Spectrometer (CLAES) instrument on the Upper Atmosphere Research Satellite (UARS). The first warming affected temperatures over a larger region, while the second produced a larger region of reversed zonal winds. Tilted baroclinic zones formed in the temperature field, and the polar vortex tilted westward with height. Narrow tongues of high PV and low N2O were drawn off the polar vortex, and irreversibly mixed. Tongues of material were drawn from low latitudes into the region between the polar vortex and the anticyclone; diabatic descent was also strongest in this region. Increased N2O over a broad region near the edge of the polar vortex indicates the importance of horizontal transport. N2O decreased in the vortex, consistent with enhanced diabatic descent during the warmings.

Beschreibung: The U.S. National Aeronautics and Space Administration (NASA) Balloon Program has been highly successful since recovering from the catastrophic balloon failure problems of the early to mid 1980s. Balloons have continued to perform at unprecedented success rates. The comprehensive research and development (R&D) effort has continued with advances being made across the spectrum of balloon related disciplines. The long duration balloon project will be transitioning from a development effort to an operational capability this year. Recently, emphasis has been placed on the development and implementation of new support systems and facilities. A new permanent launch facility at Fort Sumner, New Mexico has been established. New ground station support equipment is being implemented, and a new heavy load launch vehicle is scheduled to be implemented in 1992. The progress, status and future plans for these and other aspects of the NASA program will be presented.

Beschreibung: In this paper, preliminary results in using orthogonal and continuous wavelet transform (WT) to identify period doubling and time-frequency localization in both synthetic and real data are presented. First, the Haar WT is applied to synthetic time series derived from a simple nonlinear dynamical system- a first-order quadratic difference equation. Second, the complex Morlet WT is used to study the time-frequency localization of tropical convection based on a high-resolution Japanese Geostationary Meteorological Satellite infrared (IR) radiance dataset. The Haar WT of the synthetic time series indicates the presence and distinct separation of multiple frequencies in a period-doubling sequence. The period-doubling process generates a multiplicity of intermediate frequencies, which are manifested in the nonuniformity in time with respect to the phase of oscillations in the lower frequencies. Wavelet transform also enables the detection of extremely weak signals in high-order subharmonics resulting from the period-doubling bifurcations. These signals are either undetected or considered statistically insignificant by traditional Fourier analysis. The Morlet WT of the IR radiance dataset indicates the presence of multiple timescales, which are localized in both frequency and time. There are two regimes in the variation of IR radiance, corresponding to the wet and dry periods. Multiple timescales, ranging from semidiurnal, diurnal, synoptic, to intraseasonal with embedding structures, are active in the wet regime. In particular, synoptic variability is more prominent during the wet phase of an intensive intraseasonal cycle. These are not only consistent with, but also show more details than, previous findings by using other techniques. The phase-locking relationships among the oscillations with different time-scales suggest that both synoptic and intraseasonal variations may be mixed oscillations due to the interaction of self-excited oscillations in the tropical atmosphere and external forcings such as annual and diurnal solar radiation variations. Both examples show that WT is a powerful tool for analysis of phenomena involving multiscale interactions that exhibit localization in both frequency and time. A discussion on the caveats in the use of WT in geophysical data analysis is also presented.

Beschreibung: The catastrophic balloon failure during the first half of the 1980's identified the need for a comprehensive and continuing balloon research and development (R&D) commitment by NASA. Technical understanding was lacking in many of the disciplines and processes associated with scientific ballooning. A comprehensive balloon R&D plan was developed in 1986 and implemented in 1987. The objectives were to develop the understanding of balloon system performance, limitations, and failure mechanisms. The program consisted of five major technical areas: structures, performance and analysis, materials, chemistry and processing, and quality control. Research activitites have been conducted at NASA/Goddard Space Flight Center (GSFC)-Wallops Flight Facility (WFF), other NASA centers and government facilities, universities, and the balloon manufacturers. Several new and increased capabilities and resources have resulted from this activity. The findings, capabilities, and plan of the balloon R&D program are presented.

Beschreibung: Caps have been used to structurally reinforce scientific research balloons since the late 1950's. The scientific research balloons used by the National Aeronautics and Space Administration (NASA) use internal caps. A NASA cap placement specification does not exist since no empirical information exisits concerning cap placement. To develop a cap placement specification, NASA has completed two in-hangar inflation tests comparing the structural contributions of internal caps and external caps. The tests used small scale test balloons designed to develop the highest possible stresses within the constraints of the hangar and balloon materials. An externally capped test balloon and an internally capped test balloon were designed, built, inflated and simulated to determine the structural contributions and benefits of each. The results of the tests and simulations are presented.

Beschreibung: Observations of upper tropospheric relative humidity obtained from Raman lidar and Cross-chain Loran Atmospheric Sounding System (CLASS) sonde instruments obtained during the First ISCCP Regional Experiment (FIRE) Cirrus-II field program are compared with satellite measurements from the GOES 6.7-micron channel. The 6.7-micron channel is sensitive to water vapor integrated over a broad layer in the upper troposphere (roughly 500-200 mbar). Instantaneous measurements of the upper tropospheric relative humidity from GOES are shown to agree to within roughly 6% of the nearest lidar observations and 9% of the nearest CLASS observations. The CLASS data exhibit a slight yet systematic dry bias in upper tropospheric humidity, a result which is consistent with previous radiosonde intercomparisons. Temporal stratification of the CLASS data indicates that the magnitude of the bias is dependent upon the time of day, suggesting a solar heating effect in the radiosonde sensor. Using CLASS profiles, the impact of vertical variability in relative humidity upon the GOES upper tropospheric humidity measurements is also examined. The upper tropospheric humidity inferred from the GOES 6.7-micron channel is demonstrated to agree to within roughly 5% of the relative humidity vertically averaged over the depth of atmosphere to which the 6.7-micron channel is sensitive. The results of this study encourage the use of satellite measurements in the 6.7-micron channel to quantitatively describe the distribution and temporal evolution of the upper tropospheric humidity field.

Beschreibung: Published estimates of cloud liquid water path (LWP) from satellite-measured microwave radiation show little agreement, even about the relative magnitudes of LWP in the tropics and midlatitudes. To understand these differences and to obtain more reliable estimate, optical and microwave LWP retrieval methods are compared using the International Satellite Cloud Climatology Project (ISCCP) and special sensor microwave/imager (SSM/I) data. Errors in microwave LWP retrieval associated with uncertainties in surface, atmosphere, and cloud properties are assessed. Sea surface temperature may not produce great LWP errors, if accurate contemporaneous measurements are used in the retrieval. An uncertainty of estimated near-surface wind speed as high as 2 m/s produces uncertainty in LWP of about 5 mg/sq cm. Cloud liquid water temperature has only a small effect on LWP retrievals (rms errors less than 2 mg/sq cm), if errors in the temperature are less than 5 C; however, such errors can produce spurious variations of LWP with latitude and season. Errors in atmospheric column water vapor (CWV) are strongly coupled with errors in LWP (for some retrieval methods) causing errors as large as 30 mg/sq cm. Because microwave radiation is much less sensitive to clouds with small LWP (less than 7 mg/sq cm) than visible wavelength radiation, the microwave results are very sensitive to the process used to separate clear and cloudy conditions. Different cloud detection sensitivities in different microwave retrieval methods bias estimated LWP values. Comparing ISCCP and SSM/I LWPs, we find that the two estimated values are consistent in global, zonal, and regional means for warm, nonprecipitating clouds, which have average LWP values of about 5 mg/sq cm and occur much more frequently than precipitating clouds. Ice water path (IWP) can be roughly estimated from the differences between ISCCP total water path and SSM/I LWP for cold, nonprecipitating clouds. IWP in the winter hemisphere is about 3 times the LWP but only half the LWP in the summer hemisphere. Precipitating clouds contribute significantly to monthly, zonal mean LWP values determined from microwave, especially in the intertropical convergence zone (ITCZ), because they have almost 10 times the liquid water (cloud plus precipitation) of nonprecipitating clouds on average. There are significant differences among microwave LWP estimates associated with the treatment of precipitating clouds.

Beschreibung: Cloud radiative forcing (CRF) is the radiative impact of clouds on the Earth's radiation budget. This study examines the diurnal variations of CRF using the Earth Radiation Budget Experiment (ERBE) monthly hourly flux data and the flux data derived from the International Satellite Cloud Climatology Project (ISCCP) using the Goddard Institute for Space Studies general circulation model radiation code. The results for the months of April, July, and October 1985 and January 1986 are analyzed. We found that, in general, two data sets agreed. For longwave (LW) CRF the diurnal range over land is generally greater than that observed over oceans. For the 4-month averages the ERBE values are 15.8 W/sq m and 6.8 W/sq m for land and ocean, respectively, compared with the ISCCP calculated values of 18.4 W/sq m and 8.0 W/sq m, respectively. The land/ocean contrast is largely associated with changes in cloud amount and the temperature difference between surface and cloud top. It would be more important to note that the clear-sky flux (i.e., surface temperature) variabilities are shown to be a major contributor to the large variabilities over land. The maximum diurnal range is found to be in the summer hemisphere, and the minimum values in the winter hemisphere. It is also shown that the daytime maximum and the nighttime minimum are seen over large portions of land, whereas they occur at any local hour over most oceans. For shortwave (SW) CRF the daytime maximum values are about twice as large as monthly averages, and their highest frequency occurs at local noon, indicating that solar insolation is a primary factor for the diurnal variation of SW CRF. However, the comparison of the ERBE data with the ISCCP results demonstrated that the largest differences in the diurnal range and monthly mean of LW CRF were associated with tropical convergence zones, where clear-sky fluxes could be easily biased by persistent cloudiness and the inadequate treatment of the atmospheric water vapor.

Beschreibung: Temperature variations in the stratosphere from 1979 to 1992 are investigated using 365-day running mean of the National Meteorological Center gridded analysis temperature data. Significant variations are seen at all levels between 70 and 1 mbar. The middle stratosphere shows temperature peaks during 1982 and 1983. The upper stratosphere has significant temperature declines between 1 and 10 mbar from 1981 to 1984. Temperatures at all levels recover to near their prior values after 1984, with the 5-mbar temperatures requiring the greatest time to fully recover. The temperature declines at 1 mbar occur in both hemispheres, over all longitudes, and in every month of the year. The decreases are largest in the middle latitudes and the polar regions and during the fall and the winter months. Such temperature variations, which appear to be of natural origin, must be taken into consideration when searching for temperature trends caused by the increase of CO2 or other greenhouse gases which affect the radiative balance of the Earth-atmosphere system or stratospheric ozone.

Beschreibung: This paper presents a comparison of the water vapor distribution obtained from two general circulation models, the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM), with satellite observations of total precipitable water (TPW) from Special Sensor Microwave/Imager (SSM/I) and upper tropospheric relative humidity (UTH) from GOES. Overall, both models are successful in capturing the primary features of the observed water vapor distribution and its seasonal variation. For the ECMWF model, however, a systematic moist bias in TPW is noted over well-known stratocumulus regions in the eastern subtropical oceans. Comparison with radiosonde profiles suggests that this problem is attributable to difficulties in modeling the shallowness of the boundary layer and large vertical water vapor gradients which characterize these regions. In comparison, the CCM is more successful in capturing the low values of TPW in the stratocumulus regions, although it tends to exhibit a dry bias over the eastern half of the subtropical oceans and a corresponding moist bias in the western half. The CCM also significantly overestimates the daily variability of the moisture fields in convective regions, suggesting a problem in simulating the temporal nature of moisture transport by deep convection. Comparison of the monthly mean UTH distribution indicates generally larger discrepancies than were noted for TPW owing to the greater influence of large-scale dynamical processes in determining the distribution of UTH. In particular, the ECMWF model exhibits a distinct dry bias along the Intertropical Convergence Zone (ITCZ) and a moist bias over the subtropical descending branches of the Hadley cell, suggesting an underprediction in the strength of the Hadley circulation. The CCM, on the other hand, demonstrates greater discrepancies in UTH than are observed for the ECMWF model, but none that are as clearly correlated with well-known features of the large-scale circulation.

Beschreibung: In this paper we describe fundamental properties of an 'off-line' three-dimensional transport model, that is, a model which uses prescribed rather than predicted winds. The model is currently used primarily for studying problems of the middle atmosphere because we have not (yet) incorporated a formulation for the convective transport of trace species, a prerequisite for many tropospheric problems. The off-line model is simpler and less expensive than a model which predicts the wind and mass evolution (an 'on-line' model), but it is more complex than the two-dimensional (2-D) zonally averaged transport models often used in the study of chemistry and transport in the middle atmosphere. It thus serves as a model of intermediate complexity and can fill a useful niche for the study of transport and chemistry. We compare simulations of four tracers, released in the lower stratosphere, in both the on- and off-line models to document the difference resulting from differences in modeling the same problem with this intermediate model. These differences identify the price to be paid in going to a cheaper and simpler calculation. The off-line model transports a tracer in three dimensions. For this reason, it requires fewer approximations than 2-D transport model, which must parameterize the effects of mixing by transient and zonally asymmetric wind features. We compare simulations of the off-line model with simulations of a 2-D model for two problems. First, we compare 2-D and three-dimensional (3-D) models by simulating the emission of an NO(x)-like tracer by a fleet of high-speed aircraft. The off-line model is then used to simulate the transport of C-14 and to contrast its simulation properties to that of the host of 2-D models which participated in an identical simulation in a recent NASA model intercomparison. The off-line model is shown to be somewhat sensitive to the sampling strategy for off-line winds. Simulations with daily averaged winds are in very good qualitative agreement but are less diffusive than when driven with instantaneous winds sampled at half-hour intervals. Simulations with the off-line and 2-D models are quite similar in the middle and upper stratosphere but behave quite differently in the lower stratosphere, where the 3-D model has a substantially more vigorous circulation. The off-line model is quite realistic in its simulation of C-14. While there are still systematic differences between the 3-D calculation and the observations, the differences seem to be substantially reduced when compared with the body of 2-D simulations documented in the above mentioned NASA intercomparison, particularly at 31 deg N.

Beschreibung: SEASAT synthetic aperture radar (SAR) echoes from the sea have previously been shown to be the result of rain and winds produced by convective stroms; rain damps the surface waves and causes ech-free holes, while the diverging winds associated with downdraft generate waves and associated echoes surrounding the holes. Gust fronts are also evident. Such a snapshot from 8 July 1978 has been examined in conjunction with ground-based radar. This leads to the conclusion that the SAR storm footprints resulted from storm processes that occurred up to an hour or more prior to the snapshot. A sequence of events is discerned from the SAR imagery in which new cell growth is triggered in between the converging outflows of two preexisting cells. In turn, the new cell generates a mini-squall line along its expanding gust front. While such phenomena are well known over land, the spaceborne SAR now allows important inferences to be made about the nature and frequency of convective storms over the oceans. The storm effects on the sea have significant implications for spaceborne wind scatterometry and rainfall measurements. Some of the findings herein remain speculative because of the great distance to the Miami weather radar-the only source of corroborative data.

Beschreibung: Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Beschreibung: Climatologies of convective precipitation were derived from passive microwave observations from the Special Sensor Microwave Imager using a scattering-based algorithm of Adler et al. Data were aggregated over periods of 3-5 months using data from 4 to 5 years. Data were also stratified by satellite overpass times (primarily 06 00 and 18 00 local time). Four regions (Mexico, Amazonia, western Africa, and the western equatorial Pacific Ocean (TOGA COARE area) were chosen for their meteorological interest and relative paucity of conventional observations. The strong diurnal variation over Mexico and the southern United States was the most striking aspect of the climatologies. Pronounced morning maxima occured offshore, often in concativities in the coastline, the result of the increased convergence caused by the coastline shape. The major feature of the evening rain field was a linear-shaped maximum along the western slope of the Sierra Madre Occidental. Topography exerted a strong control on the rainfall in other areas, particularly near the Nicaragua/Honduras border and in Guatemala, where maxima in excess of 700 mm/month were located adjacent to local maxima in terrain. The correlation between the estimates and monthly gage data over the southern United States was low (0.45), due mainly to poor temporal sampling in any month and an inadequate sampling of the diurnal cycle. Over the Amazon Basin the differences in morning versus evening rainfall were complex, with an alternating series of morning/evening maxima aligned southwest to northeast from the Andes to the northeast Brazilian coast. A real extent of rainfall in Amazonia was slightly higher in the evening, but a maximum in morning precipitation was found on the Amazon River just east of Manaus. Precipitation over the water in the intertropical convergence zone (ITCZ) north of Brazil was more pronounced in the morning, and a pronounced land-/sea-breeze circulation was found along the northeast coast of Brazil. Inter-comparison of four years revealed 1992 to be the driest over Amazonia, with about a 23% decrease in mean rate compared to the 4-year mean estimated rain rate.

Beschreibung: The purpose of this Note is to present results from an analytic/experimental study that investigated the potential for passively changing blade twist through the use of extension-twist coupling. A set of composite model rotor blades was manufactured from existing blade molds for a low-twist metal helicopter rotor blade, with a view toward establishing a preliminary proof concept for extension-twist-coupled rotor blades. Data were obtained in hover for both a ballasted and unballasted blade configuration in sea-level atmospheric conditions. Test data were compared with results obtained from a geometrically nonlinear analysis of a detailed finite element model of the rotor blade developed in MSC/NASTRAN.

Beschreibung: Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Beschreibung: The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high-temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet total temperatures from ambient to 1370 K. The results of this comparison show that the first-order supersonic instability wave and the Kelvin-Hemlhlotz first-, second-, and third-order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiatio process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Beschreibung: The objective of the present work is to study the mixing characteristics of a linear array of supersonic rectangular jets under conditions of screech synchronization. The screech synchronization at a fully expanded jet Mach number of 1.61 is achieved by a precise adjustment of the internozzle spacing. To our knowledge, such an experiment on the resonant mixing of screech synchronized multiple rectangular jets has not been reported before. The results are compared with the case where the screech was suppressed in the multijet configuration.

Beschreibung: The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Beschreibung: The meteorological satellite program began in the United States as the result of the actions taken by a very small but dedicated group of people from the late 1940s to 1960. This paper provides firsthand accounts by two of these dedicated individuals. Their remarks provide an insight into the trials and tribulations they and the program encountered during these very early years. Those now active in the program, many of whom do not recall this time, might appreciate the effort of these pioneers and the legacy they left for us.

Beschreibung: The Experimental Cloud Lidar Pilot Study (ECLIPS) was initiated to obtain statistics on cloud-base height, extinction, optical depth, cloud brokenness, and surface fluxes. Two observational phases have taken place, in October-December 1989 and April-July 1991, with intensive 30-day periods being selected within the two time intervals. Data are being archived at NASA Langley Research Center and, once there, are readily available to the international scientific community. This article describes the scale of the study in terms of its international involvement and in the range of data being recorded. Lidar observations of cloud height and backscatter coefficient have been taken from a number of ground-based stations spread around the globe. Solar shortwave and infrared longwave fluxes and infrared beam radiance have been measured at the surface wherever possible. The observations have been tailored to occur around the overpass times of the NOAA weather satellites. This article describes in some detail the various retrieval methods used to obtain results on cloud-base height, extinction coefficient, and infrared emittance, paying particular attention to the uncertainties involved.

Beschreibung: The benefits of using a hypersonic waverider for spacecraft trajectory modification are presented. A waverider is a hypersonic vehicle specifically designed so that the undersurface bow shock is attached to the leading edge, which provides for the highest known lift-to-drag ratios achievable at high Mach number flight. Several viable space missions are suggested which could use such configurations for low-drag aero-assisted maneuvers in planetary atmospheres. It is shown that large changes in the spacecraft velocity vector can be accomplished with acceptably small losses in energy due to drag using a waverider aeroshell. The primary advantage of an aero-assist maneuver is suggested by comparison to a traditional gravity-assist trajectory. Some scaling laws are presented for comparing waveriders designed for different planetary atmospheres, and it is shown that the compositional differences between the terrestrial planets has a minimal impact on waverider design.

Beschreibung: In this paper, we investigate the relative importance of local vs remote control on cloud radiative forcing using a cumulus ensemble model. It is found that cloud and surface radiation forcings are much more sensitive to the mean vertical motion assoicated with large scale tropical circulation than to the local SST (sea surface temperature). When the local SST is increased with the mean vertical motion held constant, increased surface latent and sensible heat flux associated with enhanced moisture recycling is found to be the primary mechanism for cooling the ocean surface. Large changes in surface shortwave fluxes are related to changes in cloudiness induced by changes in the large scale circulation. These results are consistent with a number of earlier empirical studies, which raised concerns regarding the validity of the cirrus-thermostat hypothesis (Ramanathan and Collins, 1991). It is argued that for a better understanding of cloud feedback, both local and remote controls need to be considered and that a cumulus ensemble model is a powerful tool that should be explored for such purpose.

Beschreibung: The probability of polar stratospheric cloud (PSC) occurrence in the Antarctic and Arctic has been estimated using Stratospheric Aerosol Measurement (SAM) II aerosol extinction data from 1978 to 1989. Antarctic PSCs are typically observed by SAM II from mid-May to early November, with a maximum zonal average probability of about 0.6 at 18-20 km in August. The typical Arctic PSC season extends only from late November to early March, with a peak zonal average probability of about 0.1 in early February at 20-22 km. There is considerable year-to-year variability in Arctic PSC sightings because of changes in the dynamics of the northern polar vortex. Year-to-year variability in Antarctic sightings is most prominent in the number of late season clouds. Maximum PSC sighting probabilities in both polar regions occur in the region from 90 deg W through the Greenwich meridian to 90 deg E, where temperatures are coldest on average. Arctic sighting probabilities approach zero outside this region, but clouds have been sighted in the Antarctic at all longitudes during most months. Inferred PSC formation temperatures remain constant throughout the Arctic winter and are similar to those in early Antarctic winter. PSC formation temperatures in the Antarctic drop markedly in the 15 to 20-km region by September, a pattern consistent with the irreversible loss of HNO3 and H2O vapor in sedimenting PSC particles.

Beschreibung: A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.

Beschreibung: Five years of meteorological and hydrological data from a typical New England watershed where winter snow cover is significant were used to drive and validate two off-line land surface schemes suitable for use in the Goddard Institute for Space Studies (GISS) general circulation model (GCM): a baseline scheme that does not model the physics of a snowpack and therefore, neglects the insulating properties of snow cover; and a modified scheme in which a three-layer snowpack is modeled. Comparing baseline model results with validation data reveals several model deficiencies. Surface radiation temperatures could not adequately be modeled and the ground froze to unreasonable depths. Furthermore, because of ground cooling resulting from large surface heat fluxes to the atmosphere from the uninsulated surface, deeper model layers did not unfreeze until midsummer. As such, the normal hydrologic processes of runoff, ground water infiltration, and movement, etc., are compromised for a good part of the year. With the inclusion of a simple three-layer snow model into the baseline model, not only are the ground and surface radiation temperatures adequately modeled but all the features of snowpack ripening that characterize pack growth/ablation are simulated.

Beschreibung: A three-year climatology of satellite-estimated rainfall for the warm season for the southwest United States and Mexico has been derived from data from the Special Sensor Microwave Imager (SSM/I). The microwave data have been stratified by month (June, July, August), year (1988, 1989, 1990), and time of day (morning and evening orbits). A rain algorithm was employed that relates 86-GHz brightness temperatures to rain rate using a coupled cloud-radiative transfer model. Results identify an early evening maximum in rainfall along the western slope of the Sierra Madre Occidental during all three months. A prominent morning rainfall maximum was found off the western Mexican coast near Mazatlan in July and August. Substantial differences between morning and evening estimates were noted. To the extent that three years constitute a climatology, results of interannual variability are presented. Results are compared and contrasted to high-resolution (8 km, hourly) infrared cloud climatologies, which consist of the frequency of occurrence of cloud colder than -38 C and -58 C. This comparison has broad implications for the estimation of rainfall by simple (cloud threshold) techniques. By sampling the infrared data to approximate the time and space resolution of the microwave, we produce ratios (or adjustment factors) by which we can adjust the infrared rain estimation schemes. This produces a combined microwave/infrared rain algorithm for monthly rainfall. Using a limited set of raingage data as ground truth, an improvement (lower bias and root-mean-square error) was demonstrated by this combined technique when compared to either method alone. The diurnal variability of convection during July 1990 was examined using hourly rain estimates from the Geostationary Operational Environmental Satellite (GOES) precipitation index and the convective stratiform technique, revealing a maximum in estimated rainfall from 1800 to 2100 local time. It is in this time period when the SSM/I evening orbit occurs. A high-resolution topographic database was available to aid in interpreting the influence of topography on the rainfall patterns.

Beschreibung: An advanced Microwave Precipitation Radiometer (AMPR) has been developed and flown in the NASA ER-2-high-altitude aircraft for imaging various atmospheric and surface processes, primarily the internal structure of rain clouds. The AMPR is a scanning four-frequency total power microwave radiometer that is externally calibrated with high-emissivity warm and cold loads. Separate antenna systems allow the sampling of the 10.7- and 19.35-GHz channels at the same spatial resolution, while the 37.1- and 85.5-GHz channels utilize the same multifrequency feedhorn as the 19.35-GHz channel. Spatial resolutions from an aircraft altitude of 20-km range from 0.6 km at 85.5 GHz to 2.8 km at 19.35 and 10.7 GHz. All channels are sampled every 0.6 km in both along-track and cross-track directions, leading to a contiguous sampling pattern of the 85.5-GHz 3-dB beamwidth footprints, 2.3X oversampling of the 37.1-GHz data, and 4.4X oversampling of the 19.35- and 10.7-GHz data. Radiometer temperature sensitivities range from 0.2 to 0.5 C. Details of the system are described, including two different calibration systems and their effect on the data collected. Examples of oceanic rain systems are presented from Florida and the tropical west Pacific that illustrate the wide variety of cloud water, rainwater, and precipitation-size ice combinations that are observable from aircraft altitudes.

Beschreibung: A general circulation model (GCM) is used to model global lightning distributions and frequencies. Both total and cloud-to-ground lightning frequencies are modeled using parameterizations that relate the depth of convective clouds to lightning frequencies. The model's simulations of lightning distributions in time and space show good agreement with available observations. The model's annual mean climatology shows a global lightning frequency of 77 flashes per second, with cloud-to-ground lightning making up 25% of the total. The maximum lightning activity in the GCM occurs during the Northern Hemisphere summer, with approximately 91% of all lightning occurring over continental and coastal regions.

Beschreibung: This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Beschreibung: A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Beschreibung: Radar or satellite observations of an area generate sequences of rain-rate maps. From a gridded map a histogram of rain rates can be obtained representing the relative areas occupied by rain rates of various strengths. The histograms vary with time as precipitating systems in the area evolve and decay and amounts of convective and stratiform rain in the area change. A method of decomposing the histograms into linear combinations of a few empirical distributions with time-dependent coefficients is developed, using principal component analysis as a starting point. When applied to a tropical Atlantic dataset (GATE), two distributions emerge naturally from the analysis, resembling stratiform and convective rain-rate distributions in that they peak at low and high rain rates, respectively. The two 'modes' have different timescales and only the high-rain-rate mode has a statistically significant diurnal cycle. The ability of just two modes to describe rain variabiltiy over an area can explain why methods of estimating area-averaged rain rate from the area covered by rain rates above a certain threshold are so successful.

Beschreibung: Carbon monoxide (CO) and methane (CH4) were measured in the 0.15- to 6-km portion of the troposphere over subarctic and boreal landscapes of midcontinent and eastern Canada during July - August 1990. In the mid-continent region, Arctic air entering the region was characterized by relatively uniform CO concentrations (86-108 parts per billion by volume (ppbv)) and CH4 concentrations (1729-1764 ppbv). Local biomass burning and long-range transport of CO into the area from industrial/urban sources and distant fires did frequently produce enhanced and variable concentrations. Emissions of CH4 from the Hudson Bay lowlands was the primary source for enhanced and variable concentrations, especially at altitudes of 0.15-1 km. In eastern Canada, most of the observed variability in CO and CH4 was similar in origin to the phenomena described for the midcontinent region. However, unexpectedly low concentrations of CO (51 ppbv) and CH4 (1688 ppbv) were measured in the midtroposphere on several flights. Combined meteorological and chemical data indicated that the low CO-CH4 events were the result of long-range transport of tropical Pacific marine air to subarctic latitudes.

Beschreibung: We report features of acidic gases in the troposphere from 9 to 5000 m altitude above ground over the Canadian taiga in the summer of 1990. The measurements were conducted at a 30-m meteorological tower and from the NASA Wallops Electra aircraft as part of the joint U.S.-Canadian Arctic Boundary Layer Expedition (ABLE) 3B Northern Wetland Studies (NOWES). We sampled air for acidic gases using the mist chamber collector coupled with subsequent analysis using ion chromatography. At the tower we collected samples at two heights during a 13-day period, including diurnal studies. Using eddy flux and profile data, we estimated the biosphere/troposphere fluxes of nitric, formic, and acetic acids and sulfur dioxide. For the organic acids, emissions from the taiga in the afternoon hours and deposition during the predawn morning hours were observed. The flux intensities alone were however not high enough to explain the observed changes in mixing ratios. The measured deposition fluxes of nitric acid were high enough to have a significant influence on its mixing ratio in the boundary layer. On three days we measured vertical profiles of nitric, formic, and acetic acids through the lower to midtroposphere. We found that the chemical composition of the troposphere was extremely heterogenous. Pronounced layers of polluted air were readily apparent from our measurements. Local photochemical production and episodic long-range transport of trace components, originating from biomass burning and possibly industrial emissions, appear to have a strong influence on the composition of the troposphere and biosphere/troposphere fluxes of acidic gases at this site.

Beschreibung: In this study we perform an error analysis for cloud-top pressure retrieval using the High-Resolution Infrared Radiometric Sounder (HIRS/2) 15-microns CO2 channels for the two-layer case of transmissive cirrus overlying an overcast, opaque stratiform cloud. This analysis includes standard deviation and bias error due to instrument noise and the presence of two cloud layers, the lower of which is opaque. Instantaneous cloud pressure retrieval errors are determined for a range of cloud amounts (0.1-1.0) and cloud-top pressures (850-250 mb). Large cloud-top pressure retrieval errors are found to occur when a lower opaque layer is present underneath an upper transmissive cloud layer in the satellite field of view (FOV). Errors tend to increase with decreasing upper-cloud effective cloud amount and with decreasing cloud height (increasing pressure). Errors in retrieved upper-cloud pressure result in corresponding errors in derived effective cloud amount. For the case in which a HIRS FOV has two distinct cloud layers, the difference between the retrieved and actual cloud-top pressure is positive in all cases, meaning that the retrieved upper-cloud height is lower than the actual upper-cloud height. In addition, errors in retrieved cloud pressure are found to depend upon the lapse rate between the low-level cloud top and the surface. We examined which sounder channel combinations would minimize the total errors in derived cirrus cloud height caused by instrument noise and by the presence of a lower-level cloud. We find that while the sounding channels that peak between 700 and 1000 mb minimize random errors, the sounding channels that peak at 300-500 mb minimize bias errors. For a cloud climatology, the bias errors are most critical.

Beschreibung: We examine the influence of aggregation errors on developing estimates of regional soil-CO2 flux from temperate forests. We find daily soil-CO2 fluxes to be more sensitive to changes in soil temperatures (Q(sub 10) = 3.08) than air temperatures (Q(sub 10) = 1.99). The direct use of mean monthly air temperatures with a daily flux model underestimates regional fluxes by approximately 4%. Temporal aggregation error varies with spatial resolution. Overall, our calibrated modeling approach reduces spatial aggregation error by 9.3% and temporal aggregation error by 15.5%. After minimizing spatial and temporal aggregation errors, mature temperate forest soils are estimated to contribute 12.9 Pg C/yr to the atmosphere as carbon dioxide. Georeferenced model estimates agree well with annual soil-CO2 fluxes measured during chamber studies in mature temperate forest stands around the globe.

Beschreibung: This study examines the role of synoptic-scale eddies during the development of persistent anticyclonic height anomalies over the central North Pacific in a general circulation model under perpetual January conditions. The General Circulation Model (GCM) replicates the basic characteristics of the evolution of the anomaly patterns found in observations. The life cycle is characterized by the rapid establishment of the major anomaly center and considerably longer maintenance and decay phases, which include the development of downstream anomaly centers. The simulation also shows a realistic evolution of synoptic-scale activity beginning with enhanced activity off the east coast of Asia prior to onset, followed by a northward shift of the Pacific storm track, which lasts throughout the maintenance phase. The initial enhancement of synoptic-scale eddy activity is associated with a large-scale cyclonic anomaly that developes over Siberia several days prior to the onset of the main anticyclonic anomaly over the central North Pacific. The observations, however, show considerable interdecadel variability in the details of the composite onset behavior; it is unclear whether this variability is real or whether it reflects differences in the data assimilation systems. The role of the time mean flow and synoptic-scale eddies in the development of the persistent Pacific anomalies is studied within the context of a kinetic energy budget in which the flow is decomposed into the time-mean, low-frequency (timescales longer than 10 days), and synoptic (timescales less than 6 days) components. The budget, which is carried out for the simulation at 500 mb, shows that the initial growth of the persistent anticyclonic anomalies is associated with barotropic conversions of energy, with approximately equal contributions coming from the mean flow and the synoptic-scale eddies. After onset the barotropic conversion from the mean flow dominates, whereas the decay phase is associated with baroclinic processes within the low-frequency flow.

Beschreibung: Winds derived from a stratospheric and tropospheric data assimilation system (STRATAN) are compared with balance winds derived from National Meteorological Center/Climate Analysis Center (NMC/CAC) heights. At middle latitudes in the lower stratosphere, the results show that STRATAN winds are comparable to the balance winds. In addition STRATAN winds provide useful horizontal divergence analyses, and hence, vertical velocity fields. More generally, the STRATAN winds are useful in a more extended domain than the balanced winds. In particular, they are useful in the Tropics and the upper stratosphere where the balanced winds fail. The assimilation also captures the quasi-biennial oscillation, but does not do a good job of representing tropical waves.

Beschreibung: The Institute for Naval Oceanography, in cooperation with Naval Research Laboratories and universities, executed the Data Assimilation and Model Evaluation Experiment (DAMEE) for the Gulf Stream region during fiscal years 1991-1993. Enormous effort has gone into the preparation of several high-quality and consistent datasets for model initialization and verification. This paper describes the preparation process, the temporal and spatial scopes, the contents, the structure, etc., of these datasets. The goal of DAMEE and the need of data for the four phases of experiment are briefly stated. The preparation of DAMEE datasets consisted of a series of processes: (1) collection of observational data; (2) analysis and interpretation; (3) interpolation using the Optimum Thermal Interpolation System package; (4) quality control and re-analysis; and (5) data archiving and software documentation. The data products from these processes included a time series of 3D fields of temperature and salinity, 2D fields of surface dynamic height and mixed-layer depth, analysis of the Gulf Stream and rings system, and bathythermograph profiles. To date, these are the most detailed and high-quality data for mesoscale ocean modeling, data assimilation, and forecasting research. Feedback from ocean modeling groups who tested this data was incorporated into its refinement. Suggestions for DAMEE data usages include (1) ocean modeling and data assimilation studies, (2) diagnosis and theoretical studies, and (3) comparisons with locally detailed observations.

Beschreibung: A new framework for interpreting the origin of the tropical intraseasonal oscillation (TISO), which avoids the speed and scale selection problems in the previous theories, is proposed in this study. In this interpretation TISO is viewed as an oscillation driven by an eastward moving convective region. This convective region consists of one or more super cloud clusters originating in the Indian Ocean and terminating in mid-Pacific, and is then followed by another convective region arising in the Indian Ocean in a period of 40-50 days. Additionally, a formal analogy is pointed out between super cloud clusters and the middle-latitude baroclinic wave packets. This study includes a simulation of TISO in a 2D model to support our interpretation. Experiments were conducted with four different convection schemes. The authors advocate that the successful simulation of TISO depends on the successful simulation of super clouds cluster, which in turn depends on the successful simulation of the life cycle of cloud clusters, which further in turn depends on the choice of cumulus convection scheme. What makes a cumulus convection scheme successful in simulating TISO is discussed.

Beschreibung: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Beschreibung: Panel method computer programs are software tools of moderate cost used for solving a wide range of engineering problems. The panel code PMARC_12 (Panel Method Ames Research Center, version 12) can compute the potential flow field around complex three-dimensional bodies such as complete aircraft models. PMARC_12 is a well-documented, highly structured code with an open architecture that facilitates modifications and the addition of new features. Adjustable arrays are used throughout the code, with dimensioning controlled by a set of parameter statements contained in an include file; thus, the size of the code (i.e. the number of panels that it can handle) can be changed very quickly. This allows the user to tailor PMARC_12 to specific problems and computer hardware constraints. In addition, PMARC_12 can be configured (through one of the parameter statements in the include file) so that the code's iterative matrix solver is run entirely in RAM, rather than reading a large matrix from disk at each iteration. This significantly increases the execution speed of the code, but it requires a large amount of RAM memory. PMARC_12 contains several advanced features, including internal flow modeling, a time-stepping wake model for simulating either steady or unsteady (including oscillatory) motions, a Trefftz plane induced drag computation, off-body and on-body streamline computations, and computation of boundary layer parameters using a two-dimensional integral boundary layer method along surface streamlines. In a panel method, the surface of the body over which the flow field is to be computed is represented by a set of panels. Singularities are distributed on the panels to perturb the flow field around the body surfaces. PMARC_12 uses constant strength source and doublet distributions over each panel, thus making it a low order panel method. Higher order panel methods allow the singularity strength to vary linearly or quadratically across each panel. Experience has shown that low order panel methods can provide nearly the same accuracy as higher order methods over a wide range of cases with significantly reduced computation times; hence, the low order formulation was adopted for PMARC_12. The flow problem is solved by modeling the body as a closed surface dividing space into two regions: the region external to the surface in which an unknown velocity potential exists representing the flow field of interest, and the region internal to the surface in which a known velocity potential (representing a fictitious flow) is prescribed as a boundary condition. Both velocity potentials are required to satisfy Laplace's equation. A surface integral equation for the unknown potential external to the surface can be written by applying Green's Theorem to the external region. Using the internal potential and zero flow through the surface as boundary conditions, the unknown potential external to the surface can be solved for. When the internal flow option, which allows the analysis of closed ducts, wind tunnels, and similar internal flow problems, is selected, the geometry is modeled such that the flow field of interest is inside the geometry and the fictitious flow is outside the geometry. Items such as wings, struts, or aircraft models can be included in the internal flow problem. The time-stepping wake model gives PMARC_12 the ability to model both steady and unsteady flow problems. The wake is convected downstream from the wake-separation line by the local velocity field. With each time step, a new row of wake panels is added to the wake at the wake-separation line. Time stepping can start from time t=0 (no initial wake) or from time t=t0 (an initial wake is specified). A wide range of motions can be prescribed, including constant rates of translation, constant rate of rotation about an arbitrary axis, oscillatory translation, and oscillatory rotation about any of the three coordinate axes. Investigators interested in a visual representation of the phenomenon they are studying with PMARC_12 may want to consider obtaining the program GVS (ARC-13361), the General Visualization System. GVS is a Silicon Graphics IRIS program which was created for the purpose of supporting the scientific visualization needs of PMARC_12. GVS is available separately from COSMIC. PMARC_12 is written in standard FORTRAN 77, with the exception of the NAMELIST extension used for input. This makes the code fairly machine independent. A compiler which supports the NAMELIST extension is required. The amount of free disk space and RAM memory required for PMARC_12 will vary depending on how the code is dimensioned using the parameter statements in the include file. The recommended minimum requirements are 20Mb of free disk space and 4Mb of RAM. PMARC_12 has been successfully implemented on a Macintosh II running System 6.0.7 or 7.0 (using MPW/Language Systems Fortran 3.0), a Sun SLC running SunOS 4.1.1, an HP 720 running HP-UX 8.07, an SGI IRIS running IRIX 4.0 (it will not run under IRIX 3.x.x without modifications), an IBM RS/6000 running AIX, a DECstation 3100 running ULTRIX, and a CRAY-YMP running UNICOS 6.0 or later. Due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The standard distribution medium for PMARC_12 is a set of three 3.5 inch 800K Macintosh format diskettes and one 3.5 inch 1.44Mb Macintosh format diskette which contains an electronic copy of the documentation in MS Word 5.0 format for the Macintosh. Alternate distribution media and formats are available upon request, but these will not include the electronic version of the document. No executables are included on the distribution media. This program is an update to PMARC version 11, which was released in 1989. PMARC_12 was released in 1993. It is available only for use by United States citizens.

Beschreibung: This program determines the supersonic flowfield surrounding three-dimensional wing-body configurations of a delta wing. It was designed to provide the numerical computation of three dimensional inviscid, flowfields of either perfect or real gases about supersonic or hypersonic airplanes. The governing equations in conservation law form are solved by a finite difference method using a second order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically. The flowfield between the body and outermost shock is treated in a shock capturing fashion and therefore allows for the correct formation of secondary internal shocks . The program operates in batch mode, is in CDC update format, has been implemented on the CDC 7600, and requires more than 140K (octal) word locations.

Beschreibung: This theoretical aerodynamics program, TAD, was developed to predict the aerodynamic characteristics of vehicles with sounding rocket configurations. These slender, axisymmetric finned vehicle configurations have a wide range of aeronautical applications from rockets to high speed armament. Over a given range of Mach numbers, TAD will compute the normal force coefficient derivative, the center-of-pressure, the roll forcing moment coefficient derivative, the roll damping moment coefficient derivative, and the pitch damping moment coefficient derivative of a sounding rocket configured vehicle. The vehicle may consist of a sharp pointed nose of cone or tangent ogive shape, up to nine other body divisions of conical shoulder, conical boattail, or circular cylinder shape, and fins of trapezoid planform shape with constant cross section and either three or four fins per fin set. The characteristics computed by TAD have been shown to be accurate to within ten percent of experimental data in the supersonic region. The TAD program calculates the characteristics of separate portions of the vehicle, calculates the interference between separate portions of the vehicle, and then combines the results to form a total vehicle solution. Also, TAD can be used to calculate the characteristics of the body or fins separately as an aid in the design process. Input to the TAD program consists of simple descriptions of the body and fin geometries and the Mach range of interest. Output includes the aerodynamic characteristics of the total vehicle, or user-selected portions, at specified points over the mach range. The TAD program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 123K of 8 bit bytes. The TAD program was originally developed in 1967 and last updated in 1972.

Beschreibung: This program, which is called 'AOFA', determines the complete viscous and inviscid flow around a body of revolution at a given angle of attack and traveling at supersonic speeds. The viscous calculations from this program agree with experimental values for surface and pitot pressures and with surface heating rates. At high speeds, lee-side flows are important because the local heating is difficult to correlate and because the shed vortices can interact with vehicle components such as a canopy or a vertical tail. This program should find application in the design analysis of any high speed vehicle. Lee-side flows are difficult to calculate because thin-boundary-layer theory is not applicable and the concept of matching inviscid and viscous flow is questionable. This program uses the parabolic approximation to the compressible Navier-Stokes equations and solves for the complete inviscid and viscous regions of flow, including the pressure. The parabolic approximation results from the assumption that the stress derivatives in the streamwise direction are small in comparison with derivatives in the normal and circumferential directions. This assumption permits the equation to be solved by an implicit finite difference marching technique which proceeds downstream from the initial data point, provided the inviscid portion of flow is supersonic. The viscous cross-flow separation is also determined as part of the solution. To use this method it is necessary to first determine an initial data point in a region where the inviscid portion of the flow is supersonic. Input to this program consists of two parts. Problem description is conveyed to the program by namelist input. Initial data is acquired by the program as formatted data. Because of the large amount of run time this program can consume the program includes a restart capability. Output is in printed format and magnetic tape for further processing. This program is written in FORTRAN IV and has been implemented on a CDC 7600 with a central memory requirement of approximately 35K (octal) of 60 bit words.

Beschreibung: GEMPAK is a general meteorological software package developed at NASA/Goddard Space Flight Center. It includes programs to analyze and display surface, upper-air, and gridded data, including model output. There are very general programs to list, edit, and plot data on maps, to display profiles and time series, to draw and fill contours, to draw streamlines, to plot symbols for clouds, sky cover, and pressure tendency, and draw cross sections in the case of gridded data and sounding data. In addition, there are Barnes objective analysis programs to grid surface and upper-air data. The programs include the capabilities to derive meteorological parameters from those found in the dataset, to perform vertical interpolations of sounding data to different coordinate systems, and to compute an extensive set of gridded diagnostic quantities by specifying various nested combinations of scalars and vector arithmetic, algebraic, and differential operators. The GEMPAK 5.1 graphics/transformation subsystem, GEMPLT, provides device-independent graphics. GEMPLT also has the capability to display output in a variety of map projections or overlaid on satellite imagery. GEMPAK 5.1 is written in FORTRAN 77 and C-language and has been implemented on VAX computers under VMS and on computers running the UNIX operating system. During installation and normal use, this package occupies approximately 100Mb of hard disk space. The UNIX version of GEMPAK includes drivers for several graphic output systems including MIT's X Window System (X11,R4), Sun GKS, PostScript (color and monochrome), Silicon Graphics, and others. The VMS version of GEMPAK also includes drivers for several graphic output systems including PostScript (color and monochrome). The VMS version is delivered with the object code for the Transportable Applications Environment (TAE) program, version 4.1 which serves as a user interface. A color monitor is recommended for displaying maps on video display devices. Data for rendering regional maps is included with this package. The standard distribution medium for the UNIX version of GEMPAK 5.1 is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the VMS version of GEMPAK 5.1 is a 6250 BPI 9-track magnetic tape in DEC VAX BACKUP format. The VMS version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. This program was developed in 1985. The current version, GEMPAK 5.1, was released in 1992. The package is delivered with source code. An extensive collection of subroutine libraries allows users to format data for use by GEMPAK, to develop new programs, and to enhance existing ones.

Beschreibung: The Comprehensive Analytical Model of Rotorcraft Aerodynamics, CAMRAD, program is designed to calculate rotor performance, loads, and noise; helicopter vibration and gust response; flight dynamics and handling qualities; and system aeroelastic stability. The analysis is a consistent combination of structural, inertial, and aerodynamic models applicable to a wide range of problems and a wide class of vehicles. The CAMRAD analysis can be applied to articulated, hingeless, gimballed, and teetering rotors with an arbitrary number of blades. The rotor degrees of freedom included are blade/flap bending, rigid pitch and elastic torsion, and optionally gimbal or teeter motion. General two-rotor aircrafts can be modeled. Single main-rotor and tandem helicopter and sideby-side tilting proprotor aircraft configurations can be considered. The case of a rotor or helicopter in a wind tunnel can also be modeled. The aircraft degrees of freedom included are the six rigid body motion, elastic airframe motions, and the rotor/engine speed perturbations. CAMRAD calculates the load and motion of helicopters and airframes in two stages. First the trim solution is obtained; then the flutter, flight dynamics, and/or transient behavior can be calculated. The trim operating conditions considered include level flight, steady climb or descent, and steady turns. The analysis of the rotor includes nonlinear inertial and aerodynamic models, applicable to large blade angles and a high inflow ratio, The rotor aerodynamic model is based on two-dimensional steady airfoil characteristics with corrections for three-dimensional and unsteady flow effects, including a dynamic stall model. In the flutter analysis, the matrices are constructed that describe the linear differential equations of motion, and the equations are analyzed. In the flight dynamics analysis, the stability derivatives are calculated and the matrices are constructed that describe the linear differential equations of motion. These equations are analyzed. In the transient analysis, the rigid body equations of motion are numerically integrated, for a prescribed transient gust or control input. The CAMRAD program product is available by license for a period of ten years to domestic U.S. licensees. The licensed program product includes the CAMRAD source code, command procedures, sample applications, and one set of supporting documentation. Copies of the documentation may be purchased separately at the price indicated below. CAMRAD is written in FORTRAN 77 for the DEC VAX under VMS 4.6 with a recommended core memory of 4.04 megabytes. The DISSPLA package is necessary for graphical output. CAMRAD was developed in 1980.

Beschreibung: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Beschreibung: Panel methods are moderate cost tools for solving a wide range of engineering problems. PMARC (Panel Method Ames Research Center) is a potential flow panel code that numerically predicts flow fields around complex three-dimensional geometries. PMARC's predecessor was a panel code named VSAERO which was developed for NASA by Analytical Methods, Inc. PMARC is a new program with many additional subroutines and a well-documented code suitable for powered-lift aerodynamic predictions. The program's open architecture facilitates modifications or additions of new features. Another improvement is the adjustable size code which allows for an optimum match between the computer hardware available to the user and the size of the problem being solved. PMARC can be resized (the maximum number of panels can be changed) in a matter of minutes. Several other state-of-the-art PMARC features include internal flow modeling for ducts and wind tunnel test sections, simple jet plume modeling essential for the analysis and design of powered-lift aircraft, and a time-stepping wake model which allows the study of both steady and unsteady motions. PMARC is a low-order panel method, which means the singularities are distributed with constant strength over each panel. In many cases low-order methods can provide nearly the same accuracy as higher order methods (where the singularities are allowed to vary linearly or quadratically over each panel). Low-order methods have the advantage of a shorter computation time and do not require exact matching between panels. The flow problem is solved by assuming that the body is at rest in a moving flow field. The body is modeled as a closed surface which divides space into two regions -- one region contains the flow field of interest and the other contains a fictitious flow. External flow problems, such as a wing in a uniform stream, have the external region as the flow field of interest and the internal flow as the fictitious flow. This arrangement is reversed for internal flow problems where the internal region contains the flow field of interest and the external flow field is fictitious. In either case it is assumed that the velocity potentials in both regions satisfy Laplace's equation. PMARC has extensive geometry modeling capabilities for handling complex, three-dimensional surfaces. As with all panel methods, the geometry must be modeled by a set of panels. For convenience, the geometry is usually subdivided into several pieces and modeled with sets of panels called patches. A patch may be folded over on itself so that opposing sides of the patch form a common line. For example, wings are normally modeled with a folded patch to form the trailing edge of the wing. PMARC also has the capability to automatically generate a closing tip patch. In the case of a wing, a tip patch could be generated to close off the wing's third side. PMARC has a simple jet model for simulating a jet plume in a crossflow. The jet plume shape, trajectory, and entrainment velocities are computed using the Adler/Baron jet in crossflow code. This information is then passed back to PMARC. The wake model in PMARC is a time-stepping wake model. The wake is convected downstream from the wake separation line by the local velocity flowfield. With each time step, a new row of wake panels is added to the wake at the wake separation line. PMARC also allows an initial wake to be specified if desired, or, as a third option, no wakes need be modeled. The effective presentation of results for aerodynamics problems requires the generation of report-quality graphics. PMAPP (ARC-12751), the Panel Method Aerodynamic Plotting Program, (Sterling Software), was written for scientists at NASA's Ames Research Center to plot the aerodynamic analysis results (flow data) from PMARC. PMAPP is an interactive, color-capable graphics program for the DEC VAX or MicroVAX running VMS. It was designed to work with a variety of terminal types and hardcopy devices. PMAPP is available separately from COSMIC. PMARC was written in standard FORTRAN77 using adjustable size arrays throughout the code. Redimensioning PMARC will change the amount of disk space and memory the code requires to be able to run; however, due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The program was implemented on an Apple Macintosh (using 2.5 MB of memory) and tested on a VAX/VMS computer. The program is available on a 3.5 inch Macintosh format diskette (standard media) or in VAX BACKUP format on TK50 tape cartridge or 9-track magnetic tape. PMARC was developed in 1989.

Beschreibung: The Transonic Airfoil analysis computer code, TAIR, was developed to employ a fast, fully implicit algorithm to solve the conservative full-potential equation for the steady transonic flow field about an arbitrary airfoil immersed in a subsonic free stream. The full-potential formulation is considered exact under the assumptions of irrotational, isentropic, and inviscid flow. These assumptions are valid for a wide range of practical transonic flows typical of modern aircraft cruise conditions. The primary features of TAIR include: a new fully implicit iteration scheme which is typically many times faster than classical successive line overrelaxation algorithms; a new, reliable artifical density spatial differencing scheme treating the conservative form of the full-potential equation; and a numerical mapping procedure capable of generating curvilinear, body-fitted finite-difference grids about arbitrary airfoil geometries. Three aspects emphasized during the development of the TAIR code were reliability, simplicity, and speed. The reliability of TAIR comes from two sources: the new algorithm employed and the implementation of effective convergence monitoring logic. TAIR achieves ease of use by employing a "default mode" that greatly simplifies code operation, especially by inexperienced users, and many useful options including: several airfoil-geometry input options, flexible user controls over program output, and a multiple solution capability. The speed of the TAIR code is attributed to the new algorithm and the manner in which it has been implemented. Input to the TAIR program consists of airfoil coordinates, aerodynamic and flow-field convergence parameters, and geometric and grid convergence parameters. The airfoil coordinates for many airfoil shapes can be generated in TAIR from just a few input parameters. Most of the other input parameters have default values which allow the user to run an analysis in the default mode by specifing only a few input parameters. Output from TAIR may include aerodynamic coefficients, the airfoil surface solution, convergence histories, and printer plots of Mach number and density contour maps. The TAIR program is written in FORTRAN IV for batch execution and has been implemented on a CDC 7600 computer with a central memory requirement of approximately 155K (octal) of 60 bit words. The TAIR program was developed in 1981.

Beschreibung: GEMPAK is a general meteorological software package developed at NASA/Goddard Space Flight Center. It includes programs to analyze and display surface, upper-air, and gridded data, including model output. There are very general programs to list, edit, and plot data on maps, to display profiles and time series, to draw and fill contours, to draw streamlines, to plot symbols for clouds, sky cover, and pressure tendency, and draw cross sections in the case of gridded data and sounding data. In addition, there are Barnes objective analysis programs to grid surface and upper-air data. The programs include the capabilities to derive meteorological parameters from those found in the dataset, to perform vertical interpolations of sounding data to different coordinate systems, and to compute an extensive set of gridded diagnostic quantities by specifying various nested combinations of scalars and vector arithmetic, algebraic, and differential operators. The GEMPAK 5.1 graphics/transformation subsystem, GEMPLT, provides device-independent graphics. GEMPLT also has the capability to display output in a variety of map projections or overlaid on satellite imagery. GEMPAK 5.1 is written in FORTRAN 77 and C-language and has been implemented on VAX computers under VMS and on computers running the UNIX operating system. During installation and normal use, this package occupies approximately 100Mb of hard disk space. The UNIX version of GEMPAK includes drivers for several graphic output systems including MIT's X Window System (X11,R4), Sun GKS, PostScript (color and monochrome), Silicon Graphics, and others. The VMS version of GEMPAK also includes drivers for several graphic output systems including PostScript (color and monochrome). The VMS version is delivered with the object code for the Transportable Applications Environment (TAE) program, version 4.1 which serves as a user interface. A color monitor is recommended for displaying maps on video display devices. Data for rendering regional maps is included with this package. The standard distribution medium for the UNIX version of GEMPAK 5.1 is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the VMS version of GEMPAK 5.1 is a 6250 BPI 9-track magnetic tape in DEC VAX BACKUP format. The VMS version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. This program was developed in 1985. The current version, GEMPAK 5.1, was released in 1992. The package is delivered with source code. An extensive collection of subroutine libraries allows users to format data for use by GEMPAK, to develop new programs, and to enhance existing ones.

Beschreibung: The TAWFIVE program calculates transonic flow over a transport-type wing and fuselage. Although more complex Euler and Navier-Stokes methods are available, TAWFIVE combines a multi-grid acceleration technique in the iterative solution of the potential equation with the use of integral-form boundary-layer equations to provide a computationally efficient and sufficiently accurate design tool. TAWFIVE simplifies the solution process by breaking the problem into a loosely coupled set of modified equations. The inviscid method, using standard inviscid equations (nonlinear full potential), is valid in the "outer" region away from the wing, whereas the boundary-layer equations are valid in the thin region near the solid surface of the wing. The two types of equations are coupled by a technique of modifying surface boundary conditions for the inviscid equations. This interaction process starts with a solution of the outer flow field. Pressures are computed at the wing surface and are used to calculate the boundary layer. The boundary-layer and wake properties are then computed using a three-dimensional integral method, and the computed displacement thickness is added to the surface of the "hard" geometry. This new displaced wing surface is then regridded and the inviscid flowfield is recomputed. New values of the inviscid pressures are then used by the boundary-layer method to predict a new displacement thickness distribution. An under-relaxed update of the previously predicted displacement thickness is then made to obtain a new displacement thickness correction that is added to the "hard" geometry. These global iterations are continued until suitable convergence is obtained. Input to TAWFIVE is limited to geometric definition of the configuration, free-stream flow quantities, and iteration control parameters. The geometric input consists of the definition of a series of airfoil sections to define the wing and a series of fuselage cross sections to model the fuselage. High-aspect-ratio wings are modeled more accurately than low-aspect-ratio wings since no special provisions are made to accurately model the wing-fuselage juncture or the wingtip region. The user can specify the solution either in terms of lift or in terms of angle of attack. TAWFIVE can produce tabular output and input files for PLOT3D (COSMIC program number ARC-12779). TAWFIVE is written in FORTRAN 77 for CRAY series computers running UNICOS. The main memory requirement is 2.7Mb for execution. This program is available on a 9-track 1600 BPI UNIX tar format magnetic tape. TAWFIVE was under development from 1979 to 1989 and first released by COSMIC in 1991. CRAY and UNICOS are registered trademarks of Cray Research, Inc.

Beschreibung: A summary is presented of vortex control applications and current techniques for the control of longitudinal vortices produced by bodies, leading edges, tips and intersections. Vortex control has up till now been performed by many approaches in an empirical fashion, assisted by the essentially inviscid nature of much of longitudinal vortex behavior. Attention is given to Reynolds number sensitivities, vortex breakdown and interactions, vortex control on highly swept wings, and vortex control in juncture flows.

Beschreibung: A meteorological overview of the Arctic Boundary Layer Expedition (ABLE 3A) flight series is presented. Synoptic analyses of mid-tropospheric circulation patterns are combined with isentropic back trajectory calculations to describe the long-range (400-3000 km) atmospheric transport mechanisms and pathways of air masses to the Arctic and sub-Arctic regions of North America during July and August 1988. Siberia and the northern Pacific Ocean were found to be the two most likely source areas for 3-day transport to the study areas in Alaska. Transport to the Barrow region was frequently influenced by polar vortices and associated short-wave troughs over the Arctic Ocean, while the Bethel area was most often affected by lows migrating across the Bering Sea and the Gulf of Alaska, as well as ridges of high pressure which built into interior Alaska. July 1988 was warmer and dryer than normal over much of Alaska. As a result, the 1988 Alaska fire season was one of the most active of the past decade. Airborne lidar measurements verified the presence of biomass burning plumes on many flights, often trapped in thin subsidence layer temperature inversions. Several cases of stratosphere/troposphere exchange were noted, based upon potential vorticity analyses and aircraft lidar data, especially in the Barrow region and during transit flights to and from Alaska.

Beschreibung: The flow-field within an axial flow turbomachine, such as a turbine or compressor, is extremely complex because of three-dimensional features such as hub-corner stall, tip-leakage flows, and airfoil wakes. These flow features interact with each other and with rotor and stator airfoils inducing time-varying forces on the airfoils. These complicated rotor-stator interactions must be understood in order to design turbomachines that are light and compact as well as reliable and efficient. Two codes, STAGE-2 and STAGE-3, have been developed to compute these unsteady rotor-stator interaction flows in multistage turbomachines. An implicit, thin-layer Euler/Navier-Stokes zonal algorithm is used to compute the unsteady flow-field within both turbine and compressor configurations. Results include surface pressures and wake profiles for two-dimensional turbine and compressor configurations and surface pressures for a three-dimensional single-stage turbine configuration. The results compare well with experimental data and other unsteady computations.

Beschreibung: The aerosol optical depth over the Konza Prairie, near Manhattan, Kansas, was recorded at various locations by five separate teams. These measurements were made in support of the First ISLSCP Field Experiment (FIFE) and used to correct imagery from a variety of satellite and aircraft sensors for the effects of atmospheric scattering and absorption. The results from one instrument are reported here for 26 days in 1987 and for 7 in 1989. Daily averages span a range of 0.05 to 0.28 in the midvisible wavelengths. In addition, diurnal variations are noted in which the afternoon optical depths are greater than those of the morning by as much as 0.07. A comparison between instruments and processing techniques used to determine these aerosol optical depths is provided. The first comparisons are made using summer 1987 data. Differences of as much as 0.05 (midvisible) are observed. Although these data allow reasonable surface reflectance retrievals, they do not agree to within the performance limits typically associated with these types of instruments. With an accuracy goal of 0.02 a preseason calibration/comparison experiment was conducted at a mountain site prior to the final field campaign in 1989. Good calibration data were obtained, and good agreement (0.01, midvisible) was observed in the retrieved optical depth acquired over the Konza. By comparing data from the surface instruments at different locations, spatial inhomogeneities are determined. Then, data from the airborne tracking sunphotometer allow one to determine variations as a function of altitude. Finally, a technique is proposed for using the in situ data to establish an instrument calibration.

Beschreibung: The underlying mean and variance properties of surface net radiation, sensible-latent heat fluxes and soil heat flux are studied over the densely instrumented grassland region encompassing FIFE. Flux variability is discussed together with the problem of scaling up to area-averaged fluxes. Results are compared and contrasted for cloudy and clear situations and examined for the influence of surface-induced biophysical controls (burn and grazing treatments) and topographic controls (aspect ratios and slope factors).

Beschreibung: Airborne measurements of flux densities of carbon dioxide CO2, sensible heat, and latent heat (H2O) obtained over the First ISLSCP Field Experiment (FIFE) site during three intensive field campaigns in 1987 and one in 1989 are examined to characterize the spatial and temporal variations of CO2 and energy transfer processes. These data were collected by the National Research Council Twin Otter using low-level flight patterns, all flown at constant pressure altitude during relatively clear days. The spatial variations are larger in 1989 than in 1987 and a higher correlation is observed between the fluxes and the surface features. The temporal patterns are easier to characterize with the relatively homogeneous situation of 1987. Functional relationships obtained between fluxes of CO2 and latent heat, CO2 fluxes and greenness index, latent heat fluxes and greenness index, and between sensible heat fluxes and surface air temperature differences are presented for one day in 1987 and one in 1989 as an example of the kind of information that can be obtained from grid flights at constant pressure altitude.

Beschreibung: Radiometric measurements at 90 GHz and three sideband frequencies near the peak water vapor absorption line of 183.3 GHz were made with Advanced Microwave Moisture Sounder (AMMS) aboard the NASA DC-8 aircraft during the Global Aerosol Backscatter Experiment (GLOBE) mission over the Pacific Ocean in November 1989. Some of the measurements over the high-latitude regions (greater than 50 deg N or 50 deg S) were analyzed for the retrieval of total precipitable water less than 0.5 g/sq cm both over land and ocean surfaces. The results show that total precipitable water from a relatively dry atmosphere could be estimated with high sensitivity from these radiometric measurements. The retrieved values over ocean surface show a decrease toward the polar region as expected. The retrieved total precipitable water over land correlates positively with the aircraft radar altitude. This positive correlation is expected because the aircraft radar altitude provides a measure of atmospheric water vapor burden above the surface. Retrieved high reflectivities over land surface at 90 GHz and 183 GHz are presumably related to snow cover on the ground. This suggests that radiometric measurements at these frequencies could be used to map snow at high-latitude regions.

Beschreibung: Jet noise and jet-induced structural loads have become key issues in the design of commercial and military aircraft. Computational Fluid Dynamics (CFD) can be of use in predicting the underlying jet shear-layer instabilities and, in conjunction with classical acoustic theory, jet noise. The computational issues involved in the resolution of high Reynolds number unsteady jet flows are addressed in this paper. Once these jet flows can be accurately resolved, it should be possible to use acoustic theory to extract, for example, the far-field jet noise. An assessment of future work and computational resources required for directly computing far-field jet noise is also presented.

Beschreibung: This computer program is designed to calculate the flow fields in two-dimensional and three-dimensional axisymmetric supersonic inlets. The method of characteristics is used to compute arrays of points in the flow field. At each point the total pressure, local Mach number, local flow angle, and static pressure are calculated. This program can be used to design and analyze supersonic inlets by determining the surface compression rates and throat flow properties. The program employs the method of characteristics for a perfect gas. The basic equation used in the program is the compatibility equation which relates the change in stream angle to the change in entropy and the change in velocity. In order to facilitate the computation, the flow field behind the bow shock wave is broken into regions bounded by shock waves. In each region successive rays are computed from a surface to a shock wave until the shock wave intersects a surface or falls outside the cowl lip. As soon as the intersection occurs a new region is started and the previous region continued only in the area in which it is needed, thus eliminating unnecessary calculations. The maximum number of regions possible in the program is ten, which allows for the simultaneous calculations of up to nine shock waves. Input to this program consists of surface contours, free-stream Mach number, and various calculation control parameters. Output consists of printed and/or plotted results. For plotted results an SC-4020 or similar plotting device is required. This program is written in FORTRAN IV to be executed in the batch mode and has been implemented on a CDC 7600 with a central memory requirement of approximately 27k (octal) of 60 bit words.

Beschreibung: This program was developed to predict turbine stage performance taking into account the effects of complex passage geometries. The method uses a quasi-3D inviscid-flow analysis iteratively coupled to calculated losses so that changes in losses result in changes in the flow distribution. In this manner the effects of both the geometry on the flow distribution and the flow distribution on losses are accounted for. The flow may be subsonic or shock-free transonic. The blade row may be fixed or rotating, and the blades may be twisted and leaned. This program has been applied to axial and radial turbines, and is helpful in the analysis of mixed flow machines. This program is a combination of the flow analysis programs MERIDL and TSONIC coupled to the boundary layer program BLAYER. The subsonic flow solution is obtained by a finite difference, stream function analysis. Transonic blade-to-blade solutions are obtained using information from the finite difference, stream function solution with a reduced flow factor. Upstream and downstream flow variables may vary from hub to shroud and provision is made to correct for loss of stagnation pressure. Boundary layer analyses are made to determine profile and end-wall friction losses. Empirical loss models are used to account for incidence, secondary flow, disc windage, and clearance losses. The total losses are then used to calculate stator, rotor, and stage efficiency. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370/3033 under TSS with a central memory requirement of approximately 4.5 Megs of 8 bit bytes. This program was developed in 1985.

Beschreibung: Turbomachinery components are often connected by ducts, which are usually annular. The configurations and aerodynamic characteristics of these ducts are crucial to the optimum performance of the turbomachinery blade rows. The ANDUCT computer program was developed to calculate the velocity distribution along an arbitrary line between the inner and outer walls of an annular duct with axisymmetric swirling flow. Although other programs are available for duct analysis, the use of the velocity gradient method makes the ANDUCT program fast and convenient while requiring only modest computer resources. A fast and easy method of analyzing the flow through a duct with axisymmetric flow is the velocity gradient method, also known as the stream filament or streamline curvature method. This method has been used extensively for blade passages but has not been widely used for ducts, except for the radial equilibrium equation. In ANDUCT, a velocity gradient equation derived from the momentum equation is used to determine the velocity variation along an arbitrary straight line between the inner and outer wall of an annular duct. The velocity gradient equation is used with an assumed variation of meridional streamline curvature. Upstream flow conditions may vary between the inner and outer walls, and an assumed total pressure distribution may be specified. ANDUCT works best for well-guided passages and where the curvature of the walls is small as compared to the width of the passage. The ANDUCT program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 60K of 8 bit bytes. The ANDUCT program was developed in 1982.

Beschreibung: The Panel Code for Planar Cascades was developed as an aid for the designer of turbomachinery blade rows. The effective design of turbomachinery blade rows relies on the use of computer codes to model the flow on blade-to-blade surfaces. Most of the currently used codes model the flow as inviscid, irrotational, and compressible with solutions being obtained by finite difference or finite element numerical techniques. While these codes can yield very accurate solutions, they usually require an experienced user to manipulate input data and control parameters. Also, they often limit a designer in the types of blade geometries, cascade configurations, and flow conditions that can be considered. The Panel Code for Planar Cascades accelerates the design process and gives the designer more freedom in developing blade shapes by offering a simple blade-to-blade flow code. Panel, or integral equation, solution techniques have been used for several years by external aerodynamicists who have developed and refined them into a primary design tool of the aircraft industry. The Panel Code for Planar Cascades adapts these same techniques to provide a versatile, stable, and efficient calculation scheme for internal flow. The code calculates the compressible, inviscid, irrotational flow through a planar cascade of arbitrary blade shapes. Since the panel solution technique is for incompressible flow, a compressibility correction is introduced to account for compressible flow effects. The analysis is limited to flow conditions in the subsonic and shock-free transonic range. Input to the code consists of inlet flow conditions, blade geometry data, and simple control parameters. Output includes flow parameters at selected control points. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 590K of 8 bit bytes. This program was developed in 1982.

Beschreibung: An exact, full-potential-equation model for the steady, irrotational, homoentropic, and homoenergetic flow of a compressible, inviscid fluid through a two-dimensional planar cascade together with its appropriate boundary conditions has been derived. The CAS2D computer program numerically solves an artificially time-dependent form of the actual full-potential-equation, providing a nonrotating blade-to-blade, steady, potential transonic cascade flow analysis code. Comparisons of results with test data and theoretical solutions indicate very good agreement. In CAS2D, the governing equation is discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field is discretized by providing a boundary-fitted, nonuniform computational mesh. This mesh is generated by using a sequence of conformal mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the full-potential equation is solved iteratively by using successive line over relaxation. Possible isentropic shocks are captured by the explicit addition of an artificial viscosity in a conservative form. In addition, a four-level, consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two-dimensional cascade flows. The results from CAS2D are not directly applicable to three-dimensional, potential, rotating flows through a cascade of blades because CAS2D does not consider the effects of the Coriolis force that would be present in the three-dimensional case. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 200K of 8 bit bytes. The CAS2D program was developed in 1980.

Beschreibung: A computer program, QSONIC, has been developed for calculating the full potential, transonic quasi-three-dimensional flow through a rotating turbomachinery blade row. The need for lighter, more efficient turbomachinery components has led to the consideration of machines with fewer stages, each with blades capable of higher speeds and higher loading. As speeds increase, the numerical problems inherent in the transonic regime have to be resolved. These problems include the calculation of imbedded shock discontinuities and the dual nature of the governing equations, which are elliptic in the subcritical flow regions but become hyperbolic for supersonic zones. QSONIC provides the flow analyst with a fast and reliable means of obtaining the transonic potential flow distribution on a blade-to-blade stream surface of a stationary or rotating turbomachine blade row. QSONIC combines several promising transonic analysis techniques. The full potential equation in conservative form is discretized at each point on a body-fitted period mesh. A mass balance is calculated through the finite volume surrounding each point. Each local volume is corrected in the third dimension for any change in stream-tube thickness along the stream tube. The nonlinear equations for all volumes are of mixed type (elliptic or hyperbolic) depending on the local Mach number. The final result is a block-tridiagonal matrix formulation involving potential corrections at each grid point as the unknowns. The residual of each system of equations is solved along each grid line. At points where the Mach number exceeds unity, the density at the forward (sweeping) edge of the volume is replaced by an artificial density. This method calculates the flow field about a cascade of arbitrary two-dimensional airfoils. Three-dimensional flow is approximated in a turbomachinery blade row by correcting for stream-tube convergence and radius change in the through flow direction. Several significant assumptions were made in developing the QSONIC program, including: (1) the flow is inviscid and adiabatic, (2) the flow relative to the blade is steady, (3) the fluid is a perfect gas with constant specific heat, (4) the flow is isentropic and any discontinuities (shocks) are weak enough to be approximated as isentropic jumps, (5) there is no velocity component normal to the stream surface, and (6) the flow relative to a fixed frame in space (absolute velocity) is completely irrotational. These assumptions place some limitations on the application of QSONIC. Sharp leading edges at high incidence and high-Mach-number turbine blade trailing edges with substantial deviation will both cause large velocity peaks on the blade. In addition, the program may have difficulty converging if the passage is nearly choked. Input to QSONIC consists of case control parameters, a geometry description, upstream boundary conditions, and a rotor description. Output includes solution scheme parameters and flow field parameters. A data file is also output which contains data on the solution mesh, surface Mach numbers, surface static pressures, isomachs, and the velocity vector field. This data may be used for further processing or for plotting. The QSONIC is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 500K of 8 bit bytes. QSONIC was developed in 1982.

Beschreibung: This computer program, WIND, was developed to numerically solve the exact, full-potential equation for three-dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three-dimensional, boundary-conforming grid and iteratively solves the full-potential equation while fully accounting for both the rotating and Coriolis effects. WIND is capable of numerically analyzing the flow field about a given blade shape of the horizontal-axis type wind turbine. The rotor hub is assumed representable by a doubly infinite circular cylinder. An arbitrary number of blades may be attached to the hub and these blades may have arbitrary spanwise distributions of taper and of the twist, sweep, and dihedral angles. An arbitrary number of different airfoil section shapes may be used along the span as long as the spanwise variation of all the geometeric parameters is reasonably smooth. The numerical techniques employed in WIND involve rotated, type-dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, WIND is cabable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. Along with the three-dimensional results, WIND provides the results of the two-dimensional calculations to aid the user in locating areas of possible improvement in the aerodynamic design of the blade. Output from WIND includes the chordwise distribution of the coefficient of pressure, the Mach number, the density, and the relative velocity components at spanwise stations along the blade. In addition, the results specify local values of the lift coefficient and the tangent and axial aerodynamic force components. These are also given in integrated form expressing the total torque and the total axial force acting on the shaft. WIND can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the relative oncoming flow is subsonic. The WIND program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 253K of 8 bit bytes. WIND was developed in 1980.

Beschreibung: This computer program calculates the flow field in the supersonic portion of a mixed-compression aircraft inlet at non-zero angle of attack. This approach is based on the method of characteristics for steady three-dimensional flow. The results of this program agree with those produced by the two-dimensional method of characteristics when axisymmetric flow fields are calculated. Except in regions of high viscous interaction and boundary layer removal, the results agree well with experimental data obtained for threedimensional flow fields. The flow field in a variety of axisymmetric mixed compression inlets can be calculated using this program. The bow shock wave and the internal shock wave system are calculated using a discrete shock wave fitting procedure. The internal flow field can be calculated either with or without the discrete fitting of the internal shock wave system. The influence of molecular transport can be included in the calculation of the external flow about the forebody and in the calculation of the internal flow when internal shock waves are not discretely fitted. The viscous and thermal diffussion effects are included by treating them as correction terms in the method of characteristics procedure. Dynamic viscosity is represented by Sutherland's law and thermal conductivity is represented as a quadratic function of temperature. The thermodynamic model used is that of a thermally and calorically perfect gas. The program assumes that the cowl lip is contained in a constant plane and that the centerbody contour and cowl contour are smooth and have continuous first partial derivatives. This program cannot calculate subsonic flow, the external flow field if the bow shock wave does not exist entirely around the forebody, or the internal flow field if the bow flow field is injected into the annulus. Input to the program consists of parameters to control execution, to define the geometry, and the vehicle orientation. Output consists of a list of parameters used, solution planes, and a description of the shock waves. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 series machine with a central memory requirement of 110K (octal) of 60 bit words when it is overlayed. This flow analysis program was developed in 1978.

Beschreibung: This computer program was developed for calculating the subsonic or transonic flow on the hub-shroud mid-channel stream surface of a single blade row of a turbomachine. The design and analysis of blades for compressors and turbines ideally requires methods for analyzing unsteady, three-dimensional, turbulent viscous flow through a turbomachine. Since an exact solution is impossible at present, solutions on two-dimensional surfaces are calculated to obtain a quasi-three dimensional solution. When three-dimensional effects are important, significant information can be obtained from a solution on a cross-sectional surface of the passage normal to the flow. With this program, a solution to the equations of flow on the meridional surface can be carried out. This solution is chosen when the turbomachine under consideration has significant variation in flow properties in the hubshroud direction, especially when input is needed for use in blade-to-blade calculations. The program can also perform flow calculations for annular ducts without blades. This program should prove very useful in the design and analysis of any turbomachine. This program calculates a solution for two-dimensional, adiabatic shockfree flow. The flow must be essentially subsonic, but there may be local areas of supersonic flow. To obtain the solution, this program uses both the finite difference and the quasi-orthogonal (velocity gradient) methods combined in a way that takes maximum advantage of both. The finite-difference method solves a finite-difference equation along the meridional stream surface in a very efficient manner but is limited to subsonic velocities. This approach must be used in cases where the blade aspect ratios are above one, cases where the passage is curved, and cases with low hub-tip-ratio blades. The quasi-orthogonal method solves the velocity gradient equation on the meridional surface and is used if it is necessary to extend the range of solutions into the transonic regime. In general the blade row may be fixed or rotating and the blades may be twisted and leaned. The flow may be axial, radial, or mixed. The upstream and downstream flow conditions can vary from hub to shroud with provisions made for an approximate correction for loss of stagnation pressure. Also, viscous forces are neglected along solution mesh lines running from hub to tip. The capabilities of this program include handling of nonaxial flows without restriction, annular ducts without blades, and specified streamwise loss distributions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 700K of 8 bit bytes. This core requirement can be reduced depending on the size of the problem and the desired solution accuracy. This program was developed in 1977.

Beschreibung: A computer program has been developed for the design of supersonic rotor blades where losses are accounted for by correcting the ideal blade geometry for boundary layer displacement thickness. The ideal blade passage is designed by the method of characteristics and is based on establishing vortex flow within the passage. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal passage, and the final blade geometry is obtained by adding the displacement thicknesses to the ideal nozzle coordinates. The boundary-layer parameters are also used to calculate the aftermixing conditions downstream of the rotor blades assuming the flow mixes to a uniform state. The computer program input consists essentially of the rotor inlet and outlet Mach numbers, upper- and lower-surface Mach numbers, inlet flow angle, specific heat ratio, and total flow conditions. The program gas properties are set up for air. Additional gases require changes to be made to the program. The computer output consists of the corrected rotor blade coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 7094. This program was developed in 1971.

Beschreibung: This program obtains a transonic flow solution on a blade-to-blade surface between blades of a turbomachine. The flow must be essentially subsonic, but there may be locally supersonic flow. The solution is two-dimensional, isentropic, and shock free. The blades may be fixed or rotating. The flow may be axial, radial, or mixed, and there may be a change in stream-channel thickness in the through-flow direction. A loss in relative stagnation pressure may be accounted for. The program input consists of blade and stream-channel geometry, stagnation flow conditions, inlet and outlet flow angles, and blade-to-blade stream-channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. The transonic solution is obtained by a combination of a finite-difference, stream-function solution and a velocity-gradient solution. The finite-difference solution at a reduced weight flow provides information needed to obtain a velocity-gradient solution. This program is written in FORTRAN IV for batch execution and has been implemented on the IBM 360 computer with a central memory requirement of approximately 36K of 8 bit bytes. This program was developed in 1969 and last updated in 1979.

Beschreibung: This program is a revision of an existing program for blade-to-blade aerodynamic analysis of turbomachine blades and it is a simpler program while consistent with related programs. The analysis is for two-dimensional, subsonic, compressible (or incompressible), nonviscous flow in a circular or straight infinite cascade of blades, which may be fixed or rotating. The flow may be axial, radial, or mixed, and the stream channel thickness may change in the through-flow direction. The program input consists of blade and stream channel geometry, total flow conditions, inlet and outlet flow angles, and blade-to-blade stream channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. This program was developed on an IBM 7094/7044 DCS.

Beschreibung: This computer program gives the blade-to-blade solution of the two-dimensional, subsonic, compressible (or incompressible), nonviscous flow problem for a circular or straight infinite cascade of tandem or slotted turbomachine blades. The blades may be fixed or rotating. The flow may be axial, radial , or mixed. The method of solution is based on the stream function using an iterative solution of nonlinear finite-difference equations. These equations are solved using two major levels of iteration. The inner iteration consists of the solution of simultaneous linear equations by successive over-relaxation, using an estimated optimum over-relaxation factor. The outer iteration then changes the coefficients of the simultaneous equations to correct for compressibility. The program input consists of the basic blade geometry, the meridional stream channel coordinates, fluid stagnation conditions, weight flow and flow split through the slot, and inlet and outlet flow angles. The output includes blade surface velocities, velocity magnitude and direction throughout the passage, and the streamline coordinates.

Beschreibung: This FORTRAN IV computer program which incorporates the method of characteristics was written to assist in the design of supersonic inlets. There were two objectives: (1) to study a greater variety of supersonic inlet configurations and (2) to reduce the time required for trial-and-error procedures to arrive at optimum inlet design. The computer program was written with the intention of being able to construct a variety of inlet configurations by interchanging specific subroutines. In this manner, greater flexibility of choice was attained, and the time required to program a specific inlet configuration was greatly reduced. The second objective was accomplished by a reformulation of the boundary value problem for hyperbolic equations. By this reformulation of the boundary data, the engineering design quantities, throat Mach number and flow angle, were introduced as direct input quantities to the computer program. As a consequence of introducing the engineering parameters as input, the computer program will calculate the surface contours required to satisfy the specific throat conditions. Inviscid flow is assumed and the method used to calculate the inlet contour results in minimum distortion to the flow in the throat. This program was developed on an IBM 7094.

Beschreibung: This program represents a subsonic aerodynamic method for determining the mean camber surface of trimmed noncoplaner planforms with minimum vortex drag. With this program, multiple surfaces can be designed together to yield a trimmed configuration with minimum induced drag at some specified lift coefficient. The method uses a vortex-lattice and overcomes previous difficulties with chord loading specification. A Trefftz plane analysis is used to determine the optimum span loading for minimum drag. The program then solves for the mean camber surface of the wing associated with this loading. Pitching-moment or root-bending-moment constraints can be employed at the design lift coefficient. Sensitivity studies of vortex-lattice arrangements have been made with this program and comparisons with other theories show generally good agreement. The program is very versatile and has been applied to isolated wings, wing-canard configurations, a tandem wing, and a wing-winglet configuration. The design problem solved with this code is essentially an optimization one. A subsonic vortex-lattice is used to determine the span load distribution(s) on bent lifting line(s) in the Trefftz plane. A Lagrange multiplier technique determines the required loading which is used to calculate the mean camber slopes, which are then integrated to yield the local elevation surface. The problem of determining the necessary circulation matrix is simplified by having the chordwise shape of the bound circulation remain unchanged across each span, though the chordwise shape may vary from one planform to another. The circulation matrix is obtained by calculating the spanwise scaling of the chordwise shapes. A chordwise summation of the lift and pitching-moment is utilized in the Trefftz plane solution on the assumption that the trailing wake does not roll up and that the general configuration has specifiable chord loading shapes. VLMD is written in FORTRAN for IBM PC series and compatible computers running MS-DOS. This program requires 360K of RAM for execution. The Ryan McFarland FORTRAN compiler and PLINK86 are required to recompile the source code; however, a sample executable is provided on the diskette. The standard distribution medium for VLMD is a 5.25 inch 360K MS-DOS format diskette. VLMD was originally developed for use on CDC 6000 series computers in 1976. It was originally ported to the IBM PC in 1986, and, after minor modifications, the IBM PC port was released in 1993.

Beschreibung: This code was developed to aid design engineers in the selection and evaluation of aerodynamically efficient wing-canard and wing-horizontal-tail configurations that may employ simple hinged-flap systems. Rapid estimates of the longitudinal aerodynamic characteristics of conceptual airplane lifting surface arrangements are provided. The method is particularly well suited to configurations which, because of high speed flight requirements, must employ thin wings with highly swept leading edges. The code is applicable to wings with either sharp or rounded leading edges. The code provides theoretical pressure distributions over the wing, the canard or horizontal tail, and the deflected flap surfaces as well as estimates of the wing lift, drag, and pitching moments which account for attainable leading edge thrust and leading edge separation vortex forces. The wing planform information is specified by a series of leading edge and trailing edge breakpoints for a right hand wing panel. Up to 21 pairs of coordinates may be used to describe both the leading edge and the trailing edge. The code has been written to accommodate 2000 right hand panel elements, but can easily be modified to accommodate a larger or smaller number of elements depending on the capacity of the target computer platform. The code provides solutions for wing surfaces composed of all possible combinations of leading edge and trailing edge flap settings provided by the original deflection multipliers and by the flap deflection multipliers. Up to 25 pairs of leading edge and trailing edge flap deflection schedules may thus be treated simultaneously. The code also provides for an improved accounting of hinge-line singularities in determination of wing forces and moments. To determine lifting surface perturbation velocity distributions, the code provides for a maximum of 70 iterations. The program is constructed so that successive runs may be made with a given code entry. To make additional runs, it is necessary only to add an identification record and the namelist data that are to be changed from the previous run. This code was originally developed in 1989 in FORTRAN V on a CDC 6000 computer system, and was later ported to an MS-DOS environment. Both versions are available from COSMIC. There are only a few differences between the PC version (LAR-14458) and CDC version (LAR-14178) of AERO2S distributed by COSMIC. The CDC version has one main source code file while the PC version has two files which are easier to edit and compile on a PC. The PC version does not require a FORTRAN compiler which supports NAMELIST because a special INPUT subroutine has been added. The CDC version includes two MODIFY decks which can be used to improve the code and prevent the possibility of some infrequently occurring errors while PC-version users will have to make these code changes manually. The PC version includes an executable which was generated with the Ryan McFarland/FORTRAN compiler and requires 253K RAM and an 80x87 math co-processor. Using this executable, the sample case requires about four hours to execute on an 8MHz AT-class microcomputer with a co-processor. The source code conforms to the FORTRAN 77 standard except that it uses variables longer than six characters. With two minor modifications, the PC version should be portable to any computer with a FORTRAN compiler and sufficient memory. The CDC version of AERO2S is available in CDC NOS Internal format on a 9-track 1600 BPI magnetic tape. The PC version is available on a set of two 5.25 inch 360K MS-DOS format diskettes. IBM AT is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. CDC is a registered trademark of Control Data Corporation. NOS is a trademark of Control Data Corporation.

Beschreibung: This program provides a wing design algorithm based on modified linear theory which takes into account the effects of attainable leading-edge thrust. A primary objective of the WINGDES2 approach is the generation of a camber surface as mild as possible to produce drag levels comparable to those attainable with full theoretical leading-edge thrust. WINGDES2 provides both an analysis and a design capability and is applicable to both subsonic and supersonic flow. The optimization can be carried out for designated wing portions such as leading and trailing edge areas for the design of mission-adaptive surfaces, or for an entire planform such as a supersonic transport wing. This program replaces an earlier wing design code, LAR-13315, designated WINGDES. WINGDES2 incorporates modifications to improve numerical accuracy and provides additional capabilities. A means of accounting for the presence of interference pressure fields from airplane components other than the wing and a direct process for selection of flap surfaces to approach the performance levels of the optimized wing surfaces are included. An increased storage capacity allows better numerical representation of those configurations that have small chord leading-edge or trailing-edge design areas. WINGDES2 determines an optimum combination of a series of candidate surfaces rather than the more commonly used candidate loadings. The objective of the design is the recovery of unrealized theoretical leading-edge thrust of the input flat surface by shaping of the design surface to create a distributed thrust and thus minimize drag. The input consists of airfoil section thickness data, leading and trailing edge planform geometry, and operational parameters such as Mach number, Reynolds number, and design lift coefficient. Output includes optimized camber surface ordinates, pressure coefficient distributions, and theoretical aerodynamic characteristics. WINGDES2 is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER computer operating under NOS 2.7.1 with a central memory requirement of approximately 344K (octal) of 60 bit words. This program was developed in 1984, and last updated in 1990. CDC and CYBER are trademarks of Control Data Corporation.

Beschreibung: This study is based on the NCAR King Air aircraft and radiosonde observations on 31 October 1986 during the FIRE in Wisconsin over Oshkosh. Aircraft step-up and spiral descent flights are used to obtain kinematic and thermodynamic data. In the step-up maneuver, six different penetrations were made between 1528 and 1616 UTC. Each penetration was about 30 km long separated in the vertical by about 300 m. The time difference between the two spiral soundings was about 43 min. The aircraft descended at a rate of 1.5 m/s during these spiral soundings. Kinematic, cloud physical, and radiometric observations from various instruments are used to estimate the different terms in the moisture- and heat-budget equations. The results show that the advection terms, estimated using the mean longitudinal wind and vertical velocities, and radiative fluxes are important in forming budgets for the cirrus layers. Ice-crystal growth is significant in the upper layers. The maintenance of cirrus can be attributed to relatively warm and moist air advection, radiative cooling at upper levels, and moisture advection in the vertical. Turbulent heat and moisture fluxes are found to be significant in the low levels of cirrus.