Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable generation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Avionics Configuration Assessment for Flightdeck Interval Management: A Comparison of Avionics and Notification Methods

Flightdeck Interval Management is one of the NextGen operational concepts that FAA is sponsoring to realize requisite National Airspace System (NAS) efficiencies. Interval Management will reduce variability in temporal deviations at a position, and thereby reduce buffers typically applied by controllers - resulting in higher arrival rates, and more efficient operations. Ground software generates a strategic schedule of aircraft pairs. Air Traffic Control (ATC) provides an IM clearance with the IM spacing objective (i.e., the TTF, and at which point to achieve the appropriate spacing from this aircraft) to the IM aircraft. Pilots must dial FIM speeds into the speed window on the Mode Control Panel in a timely manner, and attend to deviations between actual speed and the instantaneous FIM profile speed. Here, the crew is assumed to be operating the aircraft with autothrottles on, with autopilot engaged, and the autoflight system in Vertical Navigation (VNAV) and Lateral Navigation (LNAV); and is responsible for safely flying the aircraft while maintaining situation awareness of their ability to follow FIM speed commands and to achieve the FIM spacing goal. The objective of this study is to examine whether three Notification Methods and four Avionics Conditions affect pilots' performance, ratings on constructs associated with performance (workload, situation awareness), or opinions on acceptability. Three Notification Methods (alternate visual and aural alerts that notified pilots to the onset of a speed target, conformance deviation from the required speed profile, and reminded them if they failed to enter the speed within 10 seconds) were examined. These Notification Methods were: VVV (visuals for all three events), VAV (visuals for all three events, plus an aural for speed conformance deviations), and AAA (visual indications and the same aural to indicate all three of these events). Avionics Conditions were defined by the instrumentation (and location) used to present IM information to crews: (1) Integrated (IM information is embedded in extant PFD (Primary Flight Display), ND (Navigation Display), EICAS (Engine Indicating and Crew Alerting System) displays); (2) EFB_Aft (IM information is only supplied in an EFB and mounted in location similar to that for MITRE's UPS work); (3) EFB_Fore (IM information is only supplied in an EFB which is mounted more forward, under the side window), and (4) EFB_Aft plus use of an AGD (the same IM information is supplied in an EFB and on an AGD, both mounted in locations similar to that in MITRE's UPS work ). Twelve commercial pilot crews flew descent scenarios (VNAV Speed with the mode control panel (MCP) speed window open until flaps extended, then VNAV Path) in a commercial transport flight simulator with realistic visual scene and communications. The results of this study serve three practical aims: (1) contribute to the down-select of avionics configuration for future assessment of the ASTAR spacing algorithm at NASA; (2) provide information useful to the FAA Human Factors Division (ANG-C1)'s mission to identify issues pertinent to flight certification of, and flight standards; (3) identify methodological considerations in support of future FIM human-in-the-loop (HITL) investigations.

Latorella, Kara A.↗

A Method for Snow Reanalysis: The Sierra Nevada (USA) Example

This work presents a state-of-the art methodology for constructing snow water equivalent (SWE) reanalysis. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data from the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlation greater than 0.95 compared with in situ SWE observations. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (31 years), and accuracy provide a unique dataset for investigating snow processes. This presentation illustrates how the reanalysis dataset was used to provide a basic accounting of the stored snowpack water in the Sierra Nevada over the last 31 years and ultimately improve real-time streamflow predictions.

Bayesian↗

Estimation of Seasonal Snow Water Equivalent Using Landsat Observations

This work presents a methodology for estimating seasonal snow water equivalent (SWE) from the use of remotely sensed Visible and Near Infrared observations from the Landsat mission. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data over the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlations with in situ SWE observations of greater than 0.95. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (32 years), and accuracy provide a unique dataset for investigating snow processes. In particular, this presentation illustrates how the reanalysis dataset was used to provide climatology of the seasonal snowfall accumulation rates, distributions, and variability over the last three decades.

SWE↗

Probability Index of Low Stratus and Fog at Dawn using Dual Geostationary Satellite Observations from COMS and FY-2D near the Korean Peninsula

We developed a new remote sensing method for detecting low stratus and fog (LSF) at dawn in terms of probability index (PI) of LSF from simultaneous stereo observations of two geostationary-orbit satellites; the Korean Communication, Ocean, and Meteorological Satellite (COMS; 128.2°E); and the Chinese FengYun satellite (FY-2D; 86.5°E). The algorithm was validated near the Korean Peninsula between the months of April and August from April 2012 to June 2015, by using surface observations at 45 meteorological stations in South Korea. The optical features of LSF were estimated by using satellite retrievals and simulated data from the radiative transfer model. The PI was calculated using the combination of three satellite-observed variables: 1) the reflectance at 0.67 μm (R(0.67)) from COMS, and 2) the FY-2D R(0.67) minus the COMS R(0.67) (△R(0.67)) and (3) the FY-2D-COMS difference in the brightness temperature difference between 3.7 and 11.0 μm (ΔBTD(3.7-11)). The three variables, adopted from the top three probability of detection (POD) scores for their fog detection thresholds: △R(0.67) (0.82) > ΔBTD(3.7-11) (0.73) > R(0.67) (0.70) > BTD(3.7-11) (0.51). The LSF PI for this algorithm was significantly better in the two case studies compared to that using COMS only (i.e., R(0.67) or BTD(3.7-11)), so that this improvement was due to △R(0.67) and ΔBTD(3.7-11). Overall, PI in the LSF spatial distribution has the merits of a high detection rate, a specific probability display, and a low rate of seasonality and variability in detection accuracy. Therefore, PI would be useful for monitoring LSF in near-real-time, and to further its forecast ability, using next-generation satellites.

fog↗

Evaluation of Global Surface Water Temperature Data Sets for use in Passive Remote Sensing of Soil Moisture

Inland open water bodies often pose a systematic error source in the passive remote sensing retrievals of soil moisture. Water temperature is a necessary variable used to compute water emissions that is required to be subtracted from satellite observation to yield actual emissions from the land portion, which in turn generates accurate soil moisture retrievals. Therefore, overestimation of soil moisture can often be corrected using concurrent water temperature data in the overall mitigation procedure. In recent years, several data sets of lake water temperature have become available, but their specifications and accuracy have rarely been investigated in the context of passive soil moisture remote sensing on a global scale. For this reason, three lake temperature products were evaluated against in-situ measurements from 2007 to 2011. The data sets include the lake surface water temperature (LSWT) from Global Observatory of Lake Responses to Environmental Change (GloboLakes), the Copernicus Global Land Operations Cryosphere and Water (C-GLOPS), as well as the lake mix-layer temperature (LMLT) from the European Centers for Medium-Range Weather Forecast (ECMWF) ERA5 Land Reanalysis. GloboLakes, C-GLOPS, and ERA5 Land have overall comparable performance with Pearson correlations (R) of 0.87, 0.92 and 0.88 in comparison with in-situ measurements. LSWT products exhibit negative median biases of −0.27 K (GloboLakes) and −0.31 K (C-GLOPS), whereas the median bias of LMLT is 1.56 K. When mapped from their respective native resolutions to a common 9 km Equal-Area Scalable Earth (EASE) Grid 2.0 projection, similar relative performance was observed. LMLT and LSWT data are closer in performance over the 9 km grid cells that exhibit a small range of lake cover fractions (0.05–0.5). Despite comparable relative performance, ERA5 Land shows great advantages in spatial coverage and temporal resolution. In summary, an integrated evaluation on data accuracy, long-term availability, global coverage, temporal resolution, and regular forward processing with modest data latency led us to conclude that LMLT from the ERA5 Land Reanalysis product represents the most optimal path for use in the development of a long-term soil moisture product.

Runze Zhang↗

Digital computer program for generating dynamic turbofan engine models (DIGTEM)

This report describes DIGTEM, a digital computer program that simulates two spool, two-stream turbofan engines. The turbofan engine model in DIGTEM contains steady-state performance maps for all of the components and has control volumes where continuity and energy balances are maintained. Rotor dynamics and duct momentum dynamics are also included. Altogether there are 16 state variables and state equations. DIGTEM features a backward-differnce integration scheme for integrating stiff systems. It trims the model equations to match a prescribed design point by calculating correction coefficients that balance out the dynamic equations. It uses the same coefficients at off-design points and iterates to a balanced engine condition. Transients can also be run. They are generated by defining controls as a function of time (open-loop control) in a user-written subroutine (TMRSP). DIGTEM has run on the IBM 370/3033 computer using implicit integration with time steps ranging from 1.0 msec to 1.0 sec. DIGTEM is generalized in the aerothermodynamic treatment of components.

Daniele, C. J.↗

A survey of nested grid techniques and their potential for use within the MASS weather prediction model

A survey of various one- and two-way interactive nested grid techniques used in hydrostatic numerical weather prediction models is presented and the advantages and disadvantages of each method are discussed. The techniques for specifying the lateral boundary conditions for each nested grid scheme are described in detail. Averaging and interpolation techniques used when applying the coarse mesh grid (CMG) and fine mesh grid (FMG) interface conditions during two-way nesting are discussed separately. The survey shows that errors are commonly generated at the boundary between the CMG and FMG due to boundary formulation or specification discrepancies. Methods used to control this noise include application of smoothers, enhanced diffusion, or damping-type time integration schemes to model variables. The results from this survey provide the information needed to decide which one-way and two-way nested grid schemes merit future testing with the Mesoscale Atmospheric Simulation System (MASS) model. An analytically specified baroclinic wave will be used to conduct systematic tests of the chosen schemes since this will allow for objective determination of the interfacial noise in the kind of meteorological setting for which MASS is designed. Sample diagnostic plots from initial tests using the analytic wave are presented to illustrate how the model-generated noise is ascertained. These plots will be used to compare the accuracy of the various nesting schemes when incorporated into the MASS model.

Koch, Steven E.↗

Next Generation Aura-OMI SO2 Retrieval Algorithm: Introduction and Implementation Status

We introduce our next generation algorithm to retrieve SO2 using radiance measurements from the Aura Ozone Monitoring Instrument (OMI). We employ a principal component analysis technique to analyze OMI radiance spectral in 310.5-340 nm acquired over regions with no significant SO2. The resulting principal components (PCs) capture radiance variability caused by both physical processes (e.g., Rayleigh and Raman scattering, and ozone absorption) and measurement artifacts, enabling us to account for these various interferences in SO2 retrievals. By fitting these PCs along with SO2 Jacobians calculated with a radiative transfer model to OMI-measured radiance spectra, we directly estimate SO2 vertical column density in one step. As compared with the previous generation operational OMSO2 PBL (Planetary Boundary Layer) SO2 product, our new algorithm greatly reduces unphysical biases and decreases the noise by a factor of two, providing greater sensitivity to anthropogenic emissions. The new algorithm is fast, eliminates the need for instrument-specific radiance correction schemes, and can be easily adapted to other sensors. These attributes make it a promising technique for producing long-term, consistent SO2 records for air quality and climate research. We have operationally implemented this new algorithm on OMI SIPS for producing the new generation standard OMI SO2 products.

Sulfur dioxide↗

The effect of normal atmospheric variability on Space Shuttle Orbiter entry performance

The Space Shuttle Orbiter is controlled during entry by a combination of aerodynamic surfaces composed of elevons, a body flap, and a rudder that flares as a speedbrake, and a reaction control system that uses 16 small rocket thrusters. To determine the effects of realistic atmospheric density variations on reentry, a series of Monte Carlo trajectories were generated for each season. The performance parameters analyzed include maximum hinge moments experienced by the aerodynamic control surfaces, the maximum deflections experienced by the aerodynamic control surfaces, a range of angle-of-sideslip values, maximum dynamic pressure, maximum stagnation point heat rate, and state conditions at touchdown. The results of the investigation indicate that the Orbiter guidance and control systems could tolerate atmospheric variability without failures in the control system hardware.

Powell, R. W.↗

Landsat thematic mapper (TM) soil variability analysis over Webster County, Iowa

Thematic mapper simulator (TMS) data acquired June 7, June 23, and July 31, 1982, and Landsat thematic mapper (TM) data acquired August 2, September 3, and October 21, 1982, over Webster County, Iowa, were examined for within-field soil effects on corn and soybean spectral signatures. It was found that patterns displayed on various computer-generated map products were in close agreement with the detailed soil survey of the area. The difference in spectral values appears to be due to a combination of subtle soil properties and crop growth patterns resulting from the different soil properties. Bands 4 (0.76-.90 micron), 5 (1.55-1.75 micron), and 7 (2.08-2.35 micron) were found to be responding to the within-field soil variability even with increasing ground cover. While these results are preliminary, they do indicate that the soil influence on the vegetation is being detected by TM and should provide improved information relating to crop and soil properties.

Thompson, D. R.↗

Gas Composition Sensing Using Carbon Nanotube Arrays

This innovation is a lightweight, small sensor for inert gases that consumes a relatively small amount of power and provides measurements that are as accurate as conventional approaches. The sensing approach is based on generating an electrical discharge and measuring the specific gas breakdown voltage associated with each gas present in a sample. An array of carbon nanotubes (CNTs) in a substrate is connected to a variable-pulse voltage source. The CNT tips are spaced appropriately from the second electrode maintained at a constant voltage. A sequence of voltage pulses is applied and a pulse discharge breakdown threshold voltage is estimated for one or more gas components, from an analysis of the current-voltage characteristics. Each estimated pulse discharge breakdown threshold voltage is compared with known threshold voltages for candidate gas components to estimate whether at least one candidate gas component is present in the gas. The procedure can be repeated at higher pulse voltages to estimate a pulse discharge breakdown threshold voltage for a second component present in the gas. The CNTs in the gas sensor have a sharp (low radius of curvature) tip; they are preferably multi-wall carbon nanotubes (MWCNTs) or carbon nanofibers (CNFs), to generate high-strength electrical fields adjacent to the tips for breakdown of the gas components with lower voltage application and generation of high current. The sensor system can provide a high-sensitivity, low-power-consumption tool that is very specific for identification of one or more gas components. The sensor can be multiplexed to measure current from multiple CNT arrays for simultaneous detection of several gas components.

Li, Jing↗

Variable-mesh method of solving differential equations

Multistep predictor-corrector method for numerical solution of ordinary differential equations retains high local accuracy and convergence properties. In addition, the method was developed in a form conducive to the generation of effective criteria for the selection of subsequent step sizes in step-by-step solution of differential equations.

Van Wyk, R.↗

Microstructural Characterization of Felsite Fragments From the Apollo Next Generation Sample Analysis (ANGSA) Double Drive Tube 73001/73002

While the Moon’s surface is dominated by mafic igneous products formed through the crystallization of the lunar magma ocean, and the subsequent eruption of mare basalts – felsic magmatic fragments (variably referred to as felsites, granites, rhyolites or granophyres) have been identified from a number of Apollo samples (e.g. 12013[1], 14321[2], 15405[3] and 73215[4]). Magmatic edifices like the Gruithuisen Domes – silicic constructs sometimes found proximal to the basaltic mare provinces filling nearside basins, may represent a petrogenetic origin. However, this is, is complicated by the fact that they apparently predate mare volcanism (e.g. [5]). Additionally, crater-counting indicates that these silicic surface features also postdate radiogenic ages of Apollo felsites [6]. Various modes of felsic magmatism have been invoked to explain the presence of felsic lithologies on the Moon, including: 1) fractional crystallization and silicate liquid immiscibility; 2) partial melting of lunar crust through basaltic underplating; and 3) fractional crystallization of the mare parent magmas. Although these models are plausible, based on known Apollo samples and lunar meteorites, they are complicated by the lack of sample lithologies with intermediate composition between basaltic magmas and the felsic components, and the fact that partial melting of many crustal rock types on the Moon (e.g., anorthosite, troctolites) are unlikely to form granites. Along with other unique magmatic lithologies, new felsites have been identified within the < 1 mm size fractions of the ANGSA double drive core tube (73001/73002) from Apollo 17 station 3, sampling the light mantle landslide deposit from the South Massif [7]. These additional examples of lunar felsite will help us to better constrain the processes that lead to the formation of these evolved lithologies. To do this, we have employed a gambit of high-resolution scanning electron microscopy (SEM) and scanning transmission electron microscopy (TEM) analytical techniques, including electron backscatter diffraction (EBSD), cathodoluminescence and high-angle annular dark field (HAADF) imaging. These data provide unique insights into the felsite mineralogy and microstructures including the coexistence of quartz and tridymite in many of the fragments. In addition, these analyses provide petrological context and assist targeted in situ secondary ion mass spectrometry (SIMS) measurements of the U-Pb systematics of accessory minerals, and the volatile abundances and D/H ratios of apatite grains identified within the clasts.

lunar↗

Split-Step Simulations to Assess the Effects of Atmospheric Boundary Layer Turbulence on the Dose Variability of N-Waves and Shaped Booms

The effects of atmospheric boundary layer turbulence on the loudness variability of a sonic boom N-wave and shaped boom are examined with split-step simulations. The shaped boom is representative of a design iteration of the NASA X-59 aircraft. Inhomogeneous atmospheric boundary layer turbulence is generated in the computational domain by a Fourier synthesis method. The N-wave and shaped boom are propagated through turbulent fields representing eight different convection levels measured at the NASA Kennedy Space Center and the NASA Armstrong Flight Research Center. Probability density functions of the formation of caustic regions along the propagation direction are computed from the N-wave results, and a parameter to collapse the caustic PDFs that accounts for both fluctuation intensities and length scales is proposed. Statistical results concerning loudness metric variability are presented, and the standard deviations of several metrics are shown to collapse across different convection levels of turbulence for small nondimensional propagation distances. The loudness metric distributions are observed to be well approximated by a normal distribution for a given range of propagation distances, and become increasingly skewed as distance increases. A model function for the dose variability is proposed, and the function parameters are found to be related to the convection level of the turbulence. The model for the dose variability distribution is compared to simulation data that were not used to find the regression parameters of the model. At several nondimensional propagation distances, agreement is observed between the model and the simulation data. These results indicate that the model may be suitable for providing quick estimates of noise dose variability in the primary carpet region across a wide range of atmospheric boundary layer conditions.

Sonic boom↗

Drop Size Distribution Measurements Supporting the NASA Global Precipitation Measurement Mission: Infrastructure and Preliminary Results

Global Precipitation Measurement Mission (GPM) retrieval algorithm validation requires datasets that characterize the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (5 -- 50 km). Key to this process is the combined use of disdrometer and polarimetric radar platforms. Here the disdrometer measurements serve as a reference for up-scaling dual-polarimetric radar observations of the PSD to the much larger volumetric sampling domain of the radar. The PSD observations thus derived provide a much larger data set for assessing DSD variability, and satellite-based precipitation retrieval algorithm assumptions, in all three spatial dimensions for a range of storm types and seasons. As one component of this effort, the GPM Ground Validation program recently acquired five 3rd generation 2D Video disdrometers as part of its Disdrometer and Radar Observations of Precipitation Facility (DROP), currently hosted in northern Alabama by the NASA Marshall Space Flight Center and the University of Alabama in Huntsville. These next-generation 2DVDs were operated and evaluated in different phases of data collection under the scanning domain of the UAH ARMOR C-band dual-polarimetric radar. During this period approximately 7500 minutes of PSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSDs were then used as input to T-matrix code for computation of polarimetric radar moments at C-band. The combined dataset was then analyzed with two basic objectives in mind: 1) the investigation of seasonal variability in the rain PSD parameters as observed by the 2DVDs; 2) the use of combined polarimetric moments and observed gamma distribution parameters in a functional form to retrieve PSD parameters in 4-D using the ARMOR radar for precipitation occurring in different seasons and for different rain system types. Preliminary results suggest that seasonal variations in the DSD parameters do occur, but are most pronounced when comparing tropical PSDs to either winter or summer convective precipitation. For example the previously documented shift to relatively smaller drop diameters in higher number concentrations for equivalent rain rate bins was observed in tropical storm rainbands occurring over Huntsville. On a more inter seasonal basis empirical fits between parameters such as D0 and ZDR do not appear to exhibit robust seasonal biases- i.e., one fit seems to work for all seasons within acceptable standard error (O[10%]) for estimates of D0. In polarimetric retrievals of the vertical variability in PSD (rain layer) for a tropical rainband we find that the Do varies with height when partitioned by specified precipitation categories (e.g., convective or stratiform, heavy and light stratiform etc.) but this variation is of order 10-20% and is smaller than the difference in D0 observed between the basic delineation of convective and stratiform precipitation types. Currently we are expanding our analysis of the vertical structure of the PSD to include several seasonally and/or dynamically-different storm system types (e.g., winter convection and stratiform events; summer mid-latitude convective etc.) sampled by ARMOR. The study will present the results of our combined analyses.

Petersen, Walter A.↗

Test results of modified electrical charged particle generator for application to fog dispersal

Modifications to a charged particle generator for use in fog dispersal applications were made and additional testing carried out. The modified nozzle, however, did not work as planned, and reported results are the unmodified nozzle. The addition of a positive displacement pump to supply the liquid water was highly successful. Measurements of the generator output current were made with a cylindrical collector system as well as with the needle probe used in previous studies. Measurements with the cylindrical collector and the needle probe showed identical agreement within the variability of the experiment. A high-voltage prove was purchased, and measurements of the corona voltage as well as the voltage variation in the charged particle jet were made. Electric fields in the vertical direction on the order of 1,000,000 v/m were measured. The voltage distribution along the centerline of the jet was compared with the numerical solutions of the Poisson equation and showed very good agreement. Velocity measurements using a pitot tube were made. The resulting measurements were compared with theoretical and other reported experimental results. The measured data showed the appropriate trends and agreed well with reported results. Based on the measured current-to-mass ratio from the charged particle generator, a calculation of the average droplet size was made. Droplet sizes were estimated to range between 0.8 and 0.4 microns. Using measured data, an analysis of the height to which the droplet can be dispersed by the charged particle generator was made. Although the mathematical model is highly simplified, the results indicated that particles would achieve heights on the order of 80 m.

Frost, W.↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

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

Siclari, Michael J.↗

A three-dimensional dynamical-chemical model of atmospheric ozone

A three-year integration of a global three-dimensional model including dynamics and simple photochemistry is used to predict ozone. Good agreement with observation is obtained for stratospheric motion patterns, meridional circulations, ozone density as a function of height and latitude, eddy transports of ozone, surface destruction of ozone, and correlations of ozone with other variables. The annual cycle of columnar ozone in high latitudes is present, but at a smaller amplitude than observed. Vertical transport of ozone downward from the main generation level at 30 km is accomplished primarily by small-scale eddy diffusion between 20 and 30 km and again near the ground; large-scale vertical transport dominates in between. The model predicts a secondary maximum in ozone mixing ratio at 45 km somewhat equatorward of the winter polar-night zone. This feature, recently observed from satellite measurements, is thought to be caused by the temperature dependence of reaction rates in the Chapman scheme. The principal deficiency of the model is an underprediction of the spring ozone concentration in high latitudes in the lower stratosphere.

Cunnold, D.↗