Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized linear model”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Chlorophyll-a Algorithms for Oligotrophic Oceans: A Novel Approach Based on Three-Band Reflectance Difference

A new empirical algorithm is proposed to estimate surface chlorophyll-a concentrations (Chl) in the global ocean for Chl less than or equal to 0.25 milligrams per cubic meters (approximately 77% of the global ocean area). The algorithm is based on a color index (CI), defined as the difference between remote sensing reflectance (R(sub rs), sr(sup -1) in the green and a reference formed linearly between R(sub rs) in the blue and red. For low Chl waters, in situ data showed a tighter (and therefore better) relationship between CI and Chl than between traditional band-ratios and Chl, which was further validated using global data collected concurrently by ship-borne and SeaWiFS satellite instruments. Model simulations showed that for low Chl waters, compared with the band-ratio algorithm, the CI-based algorithm (CIA) was more tolerant to changes in chlorophyll-specific backscattering coefficient, and performed similarly for different relative contributions of non-phytoplankton absorption. Simulations using existing atmospheric correction approaches further demonstrated that the CIA was much less sensitive than band-ratio algorithms to various errors induced by instrument noise and imperfect atmospheric correction (including sun glint and whitecap corrections). Image and time-series analyses of SeaWiFS and MODIS/Aqua data also showed improved performance in terms of reduced image noise, more coherent spatial and temporal patterns, and consistency between the two sensors. The reduction in noise and other errors is particularly useful to improve the detection of various ocean features such as eddies. Preliminary tests over MERIS and CZCS data indicate that the new approach should be generally applicable to all existing and future ocean color instruments.

Hu, Chuanmin↗
NASA NTRS2017

Evaluation of Sentinel-1A Data For Above Ground Biomass Estimation in Different Forests in India

Use of remote sensing data for mapping and monitoring of forest biomass across large spatial scales can aid in addressing uncertainties in carbon cycle. Earlier, several researchers reported on the use of Synthetic Aperture Radar (SAR) data for characterizing forest structural parameters and the above ground biomass estimation. However, these studies cannot be generalized and the algorithms cannot be applied to all types of forests without additional information on the forest physiognomy, stand structure and biomass characteristics. The radar backscatter signal also saturates as forest parameters such as biomass and the tree height increase. It is also not clear how different polarizations (VV versus VH) impact the backscatter retrievals in different forested regions. Thus, it is important to evaluate the potential of SAR data in different landscapes for characterizing forest structural parameters. In this study, the SAR data from Sentinel-1A has been used to characterize forest structural parameters including the above ground biomass from tropical forests of India. Ground based data on tree density, basal area and above ground biomass data from thirty-eight different forested sites has been collected to relate to SAR data. After the pre-processing of Sentinel 1-A data for radiometric calibration, geo-correction, terrain correction and speckle filtering, the variability in the backscatter signal in relation tree density, basal area and above biomass density has been investigated. Results from the curve fitting approach suggested exponential model between the Sentinel-1A backscatter versus tree density and above ground biomass whereas the relationship was almost linear with the basal area in the VV polarization mode. Of the different parameters, tree density could explain most of the variations in backscatter. Both VV and VH backscatter signals could explain only thirty and thirty three percent of variation in above biomass in different forest sites of India. Results also suggested saturation of the Sentinel-1A backscatter signal around hundred tonnes per hectare for VV polarization and one hundred and forty five tonnes per hectare for VH polarization. The presentation will highlight the above results in addition to potentials and limitations of Sentinel-1A data for retrieving forest structural parameters. Also, background information on different forest types of India, biomass variations and forest type mapping efforts in the region will be presented.

Data↗
NASA NTRS1993

Status of LDEF ionizing radiation measurements and analysis

At this symposium significant new data and analyses were reported in cosmic ray research, radiation dosimetry, induced radioactivity, and radiation environment modeling. Measurements of induced radioactivity and absorbed dose are nearly complete, but much analysis and modeling remains. Measurements and analyses of passive nuclear track detectors (PNTD), used to derive the cosmic ray composition and spectra, and linear energy transfer (LET) spectra, are only a few percent complete, but important results have already emerged. As one might expect at this stage of the research, some of the new information has produced questions rather than answers. Low-energy heavy nuclei detected by two experiments are not compatible with known solar or cosmic components. Various data sets on absorbed dose are not consistent, and a new trapped proton environment model does not match the absorbed dose data. A search for cosmogenic nuclei other than Be-7 on Long Duration Exposure Facility (LDEF) surfaces has produced an unexpected result, and some activation data relating to neutrons is not yet understood. Most of these issues will be resolved by the analysis of further experiment data, calibrations, or the application of the large LDEF data set that offers alternate data or analysis techniques bearing on the same problem. The scope of the papers at this symposium defy a compact technical summary. I have attempted to group the new information that I noted into the following groups: induced radioactivity; absorbed dose measurements; LET spectra and heavy ion dosimetry; environment modeling and three dimensional shielding effects; cosmogenic nuclei; and cosmic rays and other heavy ions. The papers generally are expository and have excellent illustrations, and I refer to their figures rather than reproduce them here. The general program and objectives of ionizing radiation measurements and analyses on LDEF has been described previously.

Parnell, Thomas A.↗
NASA NTRS2016

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Baurle, R. A.↗
NASA NTRS2006

M(sub W) = 7.2-7.4 Estimated for A.D. 900 Seattle Fault Earthquake by Modeling the Uplift of a Lidar-Mapped Marine Terrace

Inverse modeling of slip on the Seattle fault system, constrained by elevations of uplifted marine terraces, provides a well-constrained estimate of the magnitude of the largest known upper-crust earthquake in the Puget Sound region within the past 2500 years. The terrace elevations that constrain the slip inversion are extracted from elevation and slope images generated from LIDAR surveys of the Puget Sound collected in 1996-2002. The images reveal a single uplifted terrace, dated to 1000 cal yr B.P. near Restoration Point, which is morphologically continuous along the southern shoreline of Bainbridge Island and is visible at comparable elevations within a 25 km by 12 km region encompassing coastlines of West Seattle, Bremerton, East Bremerton, Port Orchard, and Waterman Point. Considering sea level changes since A.D. 900, the maximum uplift magnitudes of shoreline inner edges approach 9 m and are located at the southernmost coastline of Bainbridge Island and the northern tip of Waterman Point, while tilt magnitudes are modest - approaching 0.1 degrees. For each of several different Seattle fault geometry interpretations, we use a linear inversion code to solve for distributed slip on the fault surfaces. Moment magnitudes of 7.2 to 7.4 are calculated directly from the different slip solutions. In general, the greatest slip of the A.D. 900 event was confined to the frontal thrust of the Seattle fault system and was centered beneath Puget Sound between Restoration Point and Alki Point.

Jordan R. Muller↗
NASA NTRS1998

The Generation of Smooth High Speed Solar Wind from Plume-Interplume Mixing

Plumes and rays are magnetic field aligned density striations in coronal holes with different values of plasma beta. The overall plasma beta is very small in the low corona but exceeds unity beyond 15-20 solar radius. High speed solar wind reported beyond 0.3 AU is relatively smooth and uniform and known to originate from the much filamented coronal hole. Thus the obvious question is how to generate a smooth solar wind from seemingly filamentary structure. Hence one has to find a mechanism to substantiate this apparent observed (Ulysses) phenomenon. To do this we model plumes as jets (or wakes) of plasma emitted from the solar surface. The shear between a jet and its ambient is known to become unstable to the MHD Kelvin-Helmholtz ("KH") instability if the Alfven Mach number of the jet is greater than one and the uniform external magnetic field is small. Starting with a simple configuration we consider a jet of half thickness R, having uniform density and uniform internal magnetic field. The external medium has also a uniform density and uniform magnetic field. The jet is perturbed at the boundary with a linear amplitude and fixed frequency. We simulate the coronal jet using the 3D ZEUS code. The first results indicate the slab jet is unstable to the MHD KH instability at 5-10 solar radius for some angle of wave propagation. The propagating instability may smooth the filamented flow. It may also produce the entrained Alfvenic fluctuations observed by Ulysses in the high speed wind. We are at present determining the parameters which induce large growth rate. This may clarify the mystery behind the emergence of fast smooth solar wind from very filamentary structures in coronal holes. Also, using the dispersion relation already available for such a flow we obtain some general description of the instability criteria for the KH instability at a jet interface.

Parhi, Shyamsundar↗
NASA NTRS2003

Mesospheric Non-Migrating Tides Generated With Planetary Waves: Characteristics - 1

We discuss results from a modeling study with our Numerical Spectral Model (NSM) that specifically deals with the non-migrating tides generated in the mesosphere. The NSM extends from the ground to the thermosphere, incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GWs), and it describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. Accounting solely for the excitation sources of the solar migrating tides, the NSM generates through dynamical interactions also non-migrating tides in the mesosphere that are comparable in magnitude to those observed. Large non-migrating tides are produced in the diurnal and semi-diurnal oscillations for the zonal mean (m = 0) and in the semidiurnal oscillation for m = 1. In general, significant eastward and westward propagating tides are generated for all the zonal wave numbers m = 1 to 4. To identify the cause, the NSM is run without the solar heating for the zonal mean (m = 0), and the amplitudes of the resulting non-migrating tides are then negligibly small. In this case, the planetary waves are artificially suppressed, which are generated in the NSM through instabilities. This leads to the conclusion that the non-migrating tides are generated through non-linear interactions between planetary waves and migrating tides, as Forbes et al. and Talaat and Liberman had proposed. In an accompanying paper, we present results from numerical experiments, which indicate that gravity wave filtering contributes significantly to produce the non-linear coupling that is involved.

Mayr, H. G.↗
NASA NTRS2016

A Quasi-Global Approach to Improve Day-Time Satellite Surface Soil Moisture Anomalies through the Land Surface Temperature Input

Passive microwave observations from various spaceborne sensors have been linked to the soil moisture of the Earth's surface layer. A new generation of passive microwave sensors are dedicated to retrieving this variable and make observations in the single theoretically optimal L-band frequency (1-2 GHz). Previous generations of passive microwave sensors made observations in a range of higher frequencies, allowing for simultaneous estimation of additional variables required for solving the radiative transfer equation. One of these additional variables is land surface temperature, which plays a unique role in the radiative transfer equation and has an influence on the final quality of retrieved soil moisture anomalies. This study presents an optimization procedure for soil moisture retrievals through a quasi-global precipitation-based verification technique, the so-called Rvalue metric. Various land surface temperature scenarios were evaluated in which biases were added to an existing linear regression, specifically focusing on improving the skills to capture the temporal variability of soil moisture. We focus on the relative quality of the day-time (01:30 pm) observations from the Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E), as these are theoretically most challenging due to the thermal equilibrium theory, and existing studies indicate that larger improvements are possible for these observations compared to their night-time (01:30 am) equivalent. Soil moisture data used in this study were retrieved through the Land Parameter Retrieval Model (LPRM), and in line with theory, both satellite paths show a unique and distinct degradation as a function of vegetation density. Both the ascending (01:30 pm) and descending (01:30 am) paths of the publicly available and widely used AMSR-E LPRM soil moisture products were used for benchmarking purposes. Several scenarios were employed in which the land surface temperature input for the radiative transfer was varied by imposing a bias on an existing regression. These scenarios were evaluated through the Rvalue technique, resulting in optimal bias values on top of this regression. In a next step, these optimal bias values were incorporated in order to re-calibrate the existing linear regression, resulting in a quasi-global uniform LST relation for day-time observations. In a final step, day-time soil moisture retrievals using the re-calibrated land surface temperature relation were again validated through the Rvalue technique. Results indicate an average increasing Rvalue of 16.5%, which indicates a better performance obtained through the re-calibration. This number was confirmed through an independent Triple Collocation verification over the same domain, demonstrating an average root mean square error reduction of 15.3%. Furthermore, a comparison against an extensive in situ database (679 stations) also indicates a generally higher quality for the re-calibrated dataset. Besides the improved day-time dataset, this study furthermore provides insights on the relative quality of soil moisture retrieved from AMSR-E's day- and night-time observations.

soil moisture↗
NASA NTRS1992

Loss of 'complexity' and aging. Potential applications of fractals and chaos theory to senescence

The concept of "complexity," derived from the field of nonlinear dynamics, can be adapted to measure the output of physiologic processes that generate highly variable fluctuations resembling "chaos." We review data suggesting that physiologic aging is associated with a generalized loss of such complexity in the dynamics of healthy organ system function and hypothesize that such loss of complexity leads to an impaired ability to adapt to physiologic stress. This hypothesis is supported by observations showing an age-related loss of complex variability in multiple physiologic processes including cardiovascular control, pulsatile hormone release, and electroencephalographic potentials. If further research supports this hypothesis, measures of complexity based on chaos theory and the related geometric concept of fractals may provide new ways to monitor senescence and test the efficacy of specific interventions to modify the age-related decline in adaptive capacity.

Non-NASA Center↗
NASA NTRS2020

Assessments of S-NPP and N20 VIIRS DNB and M bands calibration stability and consistency using a homogeneous ground target

The S-NPP and N20 VIIRS Day-Night band (DNB) and M bands top-of-atmosphere (TOA) radiance and reflectance are calculated and a kernel-driven Bidirectional Reflectance Distribution Function (BRDF) correction model is used to correct the surface and atmosphere combined BRDF influence. Due to degradation in the DNB modulated Relative Spectral Response (RSR), the S-NPP VIIRS observed TOA DNB reflectance indicates a decrease of 1.89% and the SCIAMACHY spectra derived TOA DNB reflectance has a decrease of 1.63% for the past 8.5 years. The N20 VIIRSTOA DNB reflectance decreased 0.50% in the past 2.5 years. The DNB radiance is also compared with the integral of the M bands radiance from M4, M5, and M7. The fitting trends of DNB to the integral of M bands ratios indicate a0.48% decrease for S-NPP VIIRS and 0.14% increase for N20 VIIRS. The N20 VIIRS DNB to integral M bands ratios are closer to 1 than the S-NPP VIIRS data. The BRDF corrected reflectance comparisons show that the N20 VIIRS data is slightly lower than those of S-NPP VIIRS. The averages of the linear fit values of the N20 to S-NPP VIIRS from 2018to 2020 June are -1.97%, -4.99%, -3.36%, and -0.85% for M4, M5, M7, and DNB, respectively. Our results indicate that the S-NPP and N20 VIIRS DNB and M bands calibration have been stable. The cross-sensor differences in DNB and M bands are generally consistent with other independent studies using similar and different approaches.

Libya 4↗
NASA NTRS1996

Dust in the small Magellanic cloud. 1: Interstellar polarization and extinction data

The typical extinction curve for the Small Magellanic Cloud (SMC), in contrast to that for the Galaxy, has no bump at 2175 A and has a steeper rise into the far ultraviolet. For the Galaxy the interpretation of the extinction and, therefore, the dust content of the interstellar medium has been greatly assisted by measurements of the wavelength dependence of the polarization. For the SMC no such measurements existed. Therefore, to further elucidate the dust properties in the SMC we have for the first time measured linear polarization with five colors in the optical region of the spectrum for a sample of reddened stars. For two of these stars, for which there were no existing UV spectrophotometric measurements, but for which we measured a relatively large polarization, we have also obtained data from the International Ultraviolet Explorer (IUE) in order to study the extinction. We also attempt to correlate the SMC extinction and polarization data. The main results are: the wavelength of maximum polarization, lambda(sub max), in the SMC is typically smaller than that in the Galaxy; however, AZC 456, which shows the UV extinction bump, has a lambda(sub max) typical of that in the Galaxy, but its polarization curve is narrower and its bump is shifted to shorter wavelengths as compared to the Galaxy; and from an analysis of both the extinction and polarization data it appears that the SMC has typically smaller grains than those in the Galaxy. The absence of the extinction bump in the SMC has generally been thought to imply a lower carbon abundance in the SMC compared to the Galaxy. We interpret our results to mean that te size distribution of the interstellar grains, and not only the carbon abundance, is different in the SMC as compared to the Galaxy. In Paper 2 we present dust model fits to these observations.

Magalhaes, A. M.↗
NASA NTRS2020

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations Over the Indian Ocean

Observational analyses and ocean general circulation model experiments were performed to understand the influence of oceanic processes on warm sea surface temperature anomalies (SSTAs) prior to convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden‐Julian Oscillation (MJO), in the tropical Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool (WP) region during the November–April season of the 2001–2012 period. Seventeen of them (44%) initiated over the Seychelles‐Chagos Thermocline Ridge (SCTR) before propagating eastward to the WP and Maritime Continent. Including only global‐scale MJOs, 71% (24%) of the remaining MJO events initiated from the SCTR (WP). Four (seven) SCTR (WP) ISO events were preceded by SSTAs that were strongly influenced by wind stress‐driven oceanic processes; all four (2/7) SCTR (WP) events were MJOs. Composite analyses of the four oceanic process influenced SCTR MJO events showed that wind stress‐driven oceanic processes contributed about 30% (between 0.1 and 0.2 °C) of the preconvection warming in the SCTR, which coincides with the Intertropical Convergence Zone with mean SST > 28 °C. Reduced upwelling and entrainment played a more important role than horizontal advection in elevating the +SSTAs. Case studies revealed that the two ocean process‐influenced primary MJO SCTR events were also associated with oceanic equatorial Rossby waves. Furthermore, results of a linear atmospheric mixed layer model indicated that boundary layer processes related to SSTAs contributed approximately half of the total convergence during MJO initiation, relative to the effect of the free troposphere, for the Rossby wave‐influenced events.

MJO Initiation↗
NASA NTRSDate not supplied

TROPOMI Geometry-dependent Lambertian-Equivalent surface Reflectivity (GLER) Product for Improved Trace-Gas Retrieval

Accurate information about the reflectivity of the Earth's surface is required for most satellite retrievals of atmospheric composition, and this information is generally taken from monthly surface reflectivity climatology that neglects angular dependence. Previously we introduced Geometry-dependent Lambertian-equivalent surface reflectivity (GLER) which captures solar and satellite viewing angle dependence as well as daily and seasonal changes. GLER is calculated from simulations of Rayleigh-only top-of-atmosphere (TOA) radiances over non-Lambertian surfaces. We use NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) product over land and the wind-dependent Cox–Munk wave-facet slope distribution including water-leaving radiance over water to accounts for surface BRDF. We have developed global GLER product, previously for the Ozone Monitoring Instrument (OMI) and recently for Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) with several new improvements and updates. We have implemented the near real time daily V006 MODIS MCD43C1 BRDF data and gap-filled with a daily BRDF coefficient climatology created from 2002-2017 V006 MCD43GF data. The NASA’s Global Modelling Initiative hourly 0.25 x 0.25 deg Replay simulations are used for more accurate determination of pixel specific terrain pressure. To improve detection of seasonal snow/ice scenes, we use the 4-km snow cover product from the Interactive Multi-sensor Snow and Ice Mapping System (IMS). Finally, we use an improved version of the vector linearized discrete ordinate radiative transfer (VLIDORT) for update of the top-of-atmosphere (TOA) radiance look-up-tables (LUTs). We demonstrate how the use of GLER is beneficial to TROPOMI’s high spatial resolution (up to 3.5 km x 3.5 km) measurements to monitor atmospheric trace gas pollutants down to the sub-city scale.

TROPOMI↗
NASA NTRS1990

Granularity in the magnetic field structure of M83

Researchers recently reported Very Large Array (VLA) 20 cm continuum polarization observations of the bright, nearly face-on southern spiral galaxy M83 (NGC 5236) at a spatial resolution of 2 kpc (Sukumar and Allen 1989). The strongest linearly-polarized emission is found in two giant arcs, with typical lengths of about 30 kpc, which are situated roughly opposite each other in the dark outer regions of the galaxy at a radius of 12 kpc from the center. These regions of high polarized intensity (and hence highly-uniform magnetic field) do not coincide with any prominent spiral-arm tracers, in contrast to the expectations of simple models for the large-scale compression of magnetic field in density-wave shock fronts. From a comparison of the data with previous results at 6 cm, the authors concluded that the low polarization in the central regions of the galaxy is a result of disorder in the interstellar magnetic field. The most likely cause of this disorder is the greater star formation activity observed in the inner parts of the galaxy. The intrinsic direction of the magnetic field in the outer parts of the galaxy has also recently been determined on a length scale of 6.5 kpc from a comparison of the VLA 20 cm results with 6.3 cm observations obtained earlier with the Effelsberg telescope (Sukumar et al. 1989). There is very little Faraday rotation in the regions of the highly-polarized arcs of emission. The magnetic field in these polarized arcs is parallel to the general spiral arm structure seen in the usual optical tracers (dust, HII regions) in the bright inner parts of the galaxy disk. The maximum observed polarization at 2 kpc resolution is about 50 percent.

Allen, R. J.↗
NASA NTRS2009

A Mixed Integer Linear Program for Airport Departure Scheduling

Aircraft departing from an airport are subject to numerous constraints while scheduling departure times. These constraints include wake-separation constraints for successive departures, miles-in-trail separation for aircraft bound for the same departure fixes, and time-window or prioritization constraints for individual flights. Besides these, emissions as well as increased fuel consumption due to inefficient scheduling need to be included. Addressing all the above constraints in a single framework while allowing for resequencing of the aircraft using runway queues is critical to the implementation of the Next Generation Air Transport System (NextGen) concepts. Prior work on airport departure scheduling has addressed some of the above. However, existing methods use pre-determined runway queues, and schedule aircraft from these departure queues. The source of such pre-determined queues is not explicit, and could potentially be a subjective controller input. Determining runway queues and scheduling within the same framework would potentially result in better scheduling. This paper presents a mixed integer linear program (MILP) for the departure-scheduling problem. The program takes as input the incoming sequence of aircraft for departure from a runway, along with their earliest departure times and an optional prioritization scheme based on time-window of departure for each aircraft. The program then assigns these aircraft to the available departure queues and schedules departure times, explicitly considering wake separation and departure fix restrictions to minimize total delay for all aircraft. The approach is generalized and can be used in a variety of situations, and allows for aircraft prioritization based on operational as well as environmental considerations. We present the MILP in the paper, along with benefits over the first-come-first-serve (FCFS) scheme for numerous randomized problems based on real-world settings. The MILP results in substantially reduced delays as compared to FCFS, and the magnitude of the savings depends on the queue and departure fix structure. The MILP assumes deterministic aircraft arrival times at the runway queues. However, due to taxi time uncertainty, aircraft might arrive either earlier or later than these deterministic times. Thus, to incorporate this uncertainty, we present a method for using the MILP with "overlap discounted rolling planning horizon". The approach is based on valuing near-term decision results more than future ones. We develop a model of taxitime uncertainty based on real-world data, and then compare the baseline FCFS delays with delays using the above MILP in a simple rolling-horizon method and in the overlap discounted scheme.

Gupta, Gautam↗
NASA NTRS2003

Task II: Three-dimensional Rotating Stall Inception and Effects of Rotating Tip Clearance Asymmetry in Axial Compressors

The effects of two types of flow non-uniformity on stall inception behavior were assessed with linearized stability analyses of two compressor flow models. Response to rotating tip clearance asymmetries induced by a whirling rotor shaft or rotor height variations were investigated with a two-dimensional flow model. A 3-D compressor model was also developed to study the stability of both full-span and part-span rotating stall modes in annular geometries with radial flow variations. The studies focussed on (1) understanding what compressor designs were sensitive to these types of circumferential and spanwise flow non-uniformities, and (2) situations where 2-D stability theories were inadequate because of 3-D flow effects. Rotating tip clearance non-uniformity caused the greatest performance loss for shafts whirling at the rotating stall frequency. A whirling shaft displacement of 1 percent chord caused the stalling mass flow to rise by as much as 10 percent and the peak pressure rise to decrease by 6 percent. These changes were an order of magnitude larger than for equivalent-sized stationary or rotor-locked clearance asymmetries. Spanwise flow non-uniformity always destabilized the compressor, so that 2-D models over-predicted that stall margin compared to 3-D theory. The difference increased for compressors with larger spanwise variations of characteristic slope and reduced characteristic curvature near the peak. Differences between 2-D and 3-D stall point predictions were generally unacceptable (2 - 4 percent of flow coefficient) for single-stage configurations, but were less than 1 percent for multistage compressors. 2-D analyses predicted the wrong stall mode for specific cases of radial inlet flow distortion, mismatching and annulus area contraction, where higher-order radial modes led to stall. The stability behavior of flows with circumferential or radial non-uniformity was unified through a single stability criterion. The stall point for both cases was set by the integral around the annulus of the pressure rise characteristic slope, weighted by the amplitude of the mode shape. For the case of steady circumferential variations, this criterion reduced to the integrated mean slope (IMS) condition associated with steady inlet distortions. The rotating tip clearance asymmetry model was also used to demonstrate the feasibility of actively controlling the shaft position to suppress rotating stall. In axisymmetric mean flow, this method only stabilized the first harmonic mode, increasing the operating range until surge or higher harmonic modes became unstable.

Suder, Kenneth↗
NASA NTRS2014

High-Lift OVERFLOW Analysis of the DLR-F11 Wind Tunnel Model

In response to the 2nd AIAA CFD High Lift Prediction Workshop, the DLR-F11 wind tunnel model is analyzed using the Reynolds-averaged Navier-Stokes flow solver OVERFLOW. A series of overset grids for a bracket-off landing configuration is constructed and analyzed as part of a general grid refinement study. This high Reynolds number (15.1 million) analysis is done at multiple angles-of-attack to evaluate grid resolution effects at operational lift levels as well as near stall. A quadratic constitutive relation recently added to OVERFLOW for improved solution accuracy is utilized for side-of-body separation issues at low angles-of-attack and outboard wing separation at stall angles. The outboard wing separation occurs when the slat brackets are added to the landing configuration and is a source of discrepancy between the predictions and experimental data. A detailed flow field analysis is performed at low Reynolds number (1.35 million) after pressure tube bundles are added to the bracket-on medium grid system with the intent of better understanding bracket/bundle wake interaction with the wing's boundary layer. Localized grid refinement behind each slat bracket and pressure tube bundle coupled with a time accurate analysis are exercised in an attempt to improve stall prediction capability. The results are inconclusive and suggest the simulation is missing a key element such as boundary layer transition. The computed lift curve is under-predicted through the linear range and over-predicted near stall, and the solution from the most complete configuration analyzed shows outboard wing separation occurring behind slat bracket 6 where the experiment shows it behind bracket 5. These results are consistent with most other participants of this workshop.

High-lift↗
NASA NTRS2016

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point-source and extended-source studies rely on the modeling of this diffuse emission for accurate characterization. Here, we describe the development of the Galactic Interstellar Emission Model (GIEM),which is the standard adopted by the LAT Collaboration and is publicly available. This model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse-Compton emission produced in the Galaxy. In the GIEM, we also include large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20deg and we observe an enhanced emission toward their base extending in the north and south Galactic directions and located within approximately 4deg of the Galactic Center.

gamma rays: diffuse background – gamma rays: gen↗