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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Solar Signals in CMIP-5 Simulations: The Stratospheric Pathway

The 11 year solar-cycle component of climate variability is assessed in historical simulations of models taken from the Coupled Model Intercomparison Project, phase 5 (CMIP-5). Multiple linear regression is applied to estimate the zonal temperature, wind and annular mode responses to a typical solar cycle, with a focus on both the stratosphere and the stratospheric influence on the surface over the period approximately 1850-2005. The analysis is performed on all CMIP-5 models but focuses on the 13 CMIP-5 models that resolve the stratosphere (high-top models) and compares the simulated solar cycle signature with reanalysis data. The 11 year solar cycle component of climate variability is found to be weaker in terms of magnitude and latitudinal gradient around the stratopause in the models than in the reanalysis. The peak in temperature in the lower equatorial stratosphere (approximately 70 hPa) reported in some studies is found in the models to depend on the length of the analysis period, with the last 30 years yielding the strongest response. A modification of the Polar Jet Oscillation (PJO) in response to the 11 year solar cycle is not robust across all models, but is more apparent in models with high spectral resolution in the short-wave region. The PJO evolution is slower in these models, leading to a stronger response during February, whereas observations indicate it to be weaker. In early winter, the magnitude of the modeled response is more consistent with observations when only data from 1979-2005 are considered. The observed North Pacific high-pressure surface response during the solar maximum is only simulated in some models, for which there are no distinguishing model characteristics. The lagged North Atlantic surface response is reproduced in both high- and low-top models, but is more prevalent in the former. In both cases, the magnitude of the response is generally lower than in observations.

stratosphere↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗
NASA NTRS2003

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗
NASA NTRS2021

Improved Rainfall Data in the Philippines through Concurrent Use of GPM IMERG and Ground-Based Measurements

The availability of accurate and reliable rainfall data that are applicable to various phenomenological, climatological, and modeling studies is important, especially in the Philippines, which is considered to be highly vulnerable to natural hazards and a changing climate. The presented strategy involved constructing a dataset consisting of synoptic data, automatic rain gauge (ARG) measurements, and satellite data that are co-registered, consistent, and formatted in the same manner. Although sparse in number, the synoptic stations provide the most accurate rainfall information and were used as the baseline for creating the dataset. The ARGs that are within a distance of 1 km to the synoptic stations were used to determine the correction factors needed to make the synoptic and ARG data consistent. Subsequently, the corrected ARGs were used to make the satellite IMERG data consistent with both ARG and synoptic data. In case of the latter, only IMERG pixels with at least 10 ARGs within the relatively large footprint of the satellite sensor were used in estimating the required correction parameters derived from a combination of a power transform and linear regression correction techniques. The final results show good agreement of synoptic and corrected ARG data with correlation coefficients of 0.94 and 0.97 for the 10 day and monthly data, respectively, and improvement in the linear regression slope from 0.67 to 0.90 for 10 day data, and 0.70 to 0.94 for monthly data. In addition, the corrected ARG data agree well with the corrected IMERG data, with correlation coefficients of 0.88 and 0.93 for the 10 day and monthly data, respectively, and an improvement in slope from 0.66 to 0.87 for 10 day data, and 0.74 to 0.99 for monthly data. The merit of using a combined dataset is illustrated through comparative analyses of the IMERG data and spatially interpolated synoptic and ARG data. The results show general agreements in spatial patterns of rainfall across the datasets, especially in areas where in situ measurements are recorded. The observed discrepancy when ground data is limited emphasizes the need for satellite IMERG data to obtain the true spatial patterns of rainfall distribution.

rainfall↗
NASA NTRS1998

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Overview - Part 1

Knowledge of the amount of water in the soil is of great importance to many earth science disciplines. Soil moisture is a key variable in controlling the exchange of water and energy between the land surface and the atmosphere. Thus, soil moisture information is valuable in a wide range of applications including weather and climate, runoff potential and flood control, early warning of droughts, irrigation, crop yield forecasting, soil erosion, reservoir management, geotechnical engineering, and water quality. Despite the importance of soil moisture information, widespread and continuous measurements of soil moisture are not possible today. Although many earth surface conditions can be measured from satellites, we still cannot adequately measure soil moisture from space. Research in soil moisture remote sensing began in the mid 1970s shortly after the surge in satellite development. Recent advances in remote sensing have shown that soil moisture can be measured, at least qualitatively, by several methods. Quantitative measurements of moisture in the soil surface layer have been most successful using both passive and active microwave remote sensing, although complications arise from surface roughness and vegetation type and density. Early attempts to measure soil moisture from space-borne microwave instruments were hindered by what is now considered sub-optimal wavelengths (shorter than 5 cm) and the coarse spatial resolution of the measurements. L-band frequencies between 1 and 3 GHz (10-30 cm) have been deemed optimal for detection of soil moisture in the upper few centimeters of soil. The Electronically Steered Thinned Array Radiometer (ESTAR), an aircraft-based instrument operating a 1,4 GHz, has shown great promise for soil moisture determination. Initiatives are underway to develop a similar instrument for space. Existing space-borne synthetic aperture radars (SARS) operating at C- and L-band have also shown some potential to detect surface wetness. The advantage of radar is its much higher resolution than passive microwave systems, but it is currently hampered by surface roughness effects and the lack of a good algorithm based on a single frequency and single polarization. In addition, its repeat frequency is generally low (about 40 days). In the meantime, two new radiometers offer some hope for remote sensing of soil moisture from space. The Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), launched in November 1997, possesses a 10.65 GHz channel and the Advanced Microwave Scanning Radiometer (AMSR) on both the ADEOS-11 and Earth Observing System AM-1 platforms to be launched in 1999 possesses a 6.9 GHz channel. Aside from issues about interference from vegetation, the coarse resolution of these data will provide considerable challenges pertaining to their application. The resolution of TMI is about 45 km and that of AMSR is about 70 km. These resolutions are grossly inconsistent with the scale of soil moisture processes and the spatial variability of factors that control soil moisture. Scale disparities such as these are forcing us to rethink how we assimilate data of various scales in hydrologic models. Of particular interest is how to assimilate soil moisture data by reconciling the scale disparity between what we can expect from present and future remote sensing measurements of soil moisture and modeling soil moisture processes. It is because of this disparity between the resolution of space-based sensors and the scale of data needed for capturing the spatial variability of soil moisture and related properties that remote sensing of soil moisture has not met with more widespread success. Within a single footprint of current sensors at the wavelengths optimal for this application, in most cases there is enormous heterogeneity in soil moisture created by differences in landcover, soils and topography, as well as variability in antecedent precipitation. It is difficult to interpret the meaning of 'mean' soil moisture under such conditions and even more difficult to apply such a value. Because of the non-linear relationships between near-surface soil moisture and other variables of interest, such as surface energy fluxes and runoff, mean soil moisture has little applicability at such large scales. It is for these reasons that the use of remote sensing in conjunction with a hydrologic model appears to be of benefit in capturing the complete spatial and temporal structure of soil moisture. This paper is Part I of a four-part series describing a method for intermittently assimilating remotely-sensed soil moisture information to improve performance of a distributed land surface hydrology model. The method, summarized in section II, involves the following components, each of which is detailed in the indicated section of the paper or subsequent papers in this series: Forward radiative transfer model methods (section II and Part IV); Use of a Kalman filter to assimilate remotely-sensed soil moisture estimates with the model profile (section II and Part IV); Application of a soil hydrology model to capture the continuous evolution of the soil moisture profile within and below the root zone (section III); Statistical aggregation techniques (section IV and Part II); Disaggregation techniques using a neural network approach (section IV and Part III); and Maximum likelihood and Bayesian algorithms for inversely solving for the soil moisture profile in the upper few cm (Part IV).

Crosson, William L.↗
NASA NTRS1980

Changes in plasma vasopressin during motion sickness in cats

Changes in levels of plasma vasopressin (AVP) and cortisol (C) have been shown to be correlated with motion sickness and nausea in man. As part of the research aimed at validation of the cat as an appropriate animal model for motion sickness research, levels of these hormones were investigated in the cat during motion sickness elicited by vertical linear acceleration of approximately 0.6 Hz and 1 +/- 0.6 G. In Study 1, 15 cats previously screened for susceptibility to motion sickness were prepared with indwelling jugular catheters to permit withdrawl of blood with minimal disruption of the stimulus and minimum stress to the animal. AVP and C were measured in blood samples obtained during exposure to vertical linear acceleration and during control sessions in which the animals were placed in the stationary apparatus. 10 min and 1 min prior to duration; 1, 5, 10, and 20 min after start of motion. Total duration of exposure to motion was 20 min. The data indicate that both AVP and C are elevated during exposure to motion if emesis occurs. AVP reaches maximum levels during or about the same time as emesis, while C increases gradually throughout the period of vertical acceleration. In Study 2, four cats were prepared with indwelling catheters and AVP was measured in blood withdrawn during exposure to the vertical linear acceleration. A single pre-motion sample consisting of three samples drawn 5 min prior to motion onset. Two series of samples consisting of three samples drawn at 3-min intervals were obtained during motion. The first series was initiated at emesis, and the second 25 min after emesis. Results show that levels of circulating AVP were elevated (2 to 27 times the control and pre-motion levels) in the samples taken during emesis and decreased, but remained 1 to 6 times above the pre-motion or control levels within 25 min. The results of these two studies indicate that AVP is elevated during motion-produced emesis than is C. These findings are in general agreement with those obtained from humans under motion sickness conditions, and indicate that it is appropriate to continue to use the cat in studies of hormone changes during motion sickness.

R Fox↗
NASA NTRS2000

Laser-Induced Fluorescence Measurements and Modeling of Nitric Oxide in Counterflow Diffusion Flames

The feasibility of making quantitative nonintrusive NO concentration ([NO]) measurements in nonpremixed flames has been assessed by obtaining laser-induced fluorescence (LIF) measurements of [NO] in counterflow diffusion flames at atmospheric and higher pressures. Comparisons at atmospheric pressure between laser-saturated fluorescence (LSF) and linear LIF measurements in four diluted ethane-air counterflow diffusion flames with strain rates from 5 to 48/s yielded excellent agreement from fuel-lean to moderately fuel-rich conditions, thus indicating the utility of a model-based quenching correction technique, which was then extended to higher pressures. Quantitative LIF measurements of [NO] in three diluted methane-air counterflow diffusion flames with strain rates from 5 to 35/s were compared with OPPDIF model predictions using the GRI (version 2.11) chemical kinetic mechanism. The comparisons revealed that the GRI mechanism underpredicts prompt-NO by 30-50% at atmospheric pressure. Based on these measurements, a modified reaction rate coefficient for the prompt-NO initiation reaction was proposed which causes the predictions to match experimental data. Temperature measurements using thin filament pyrometry (TFP) in conjunction with a new calibration method utilizing a near-adiabatic H2-air Hencken burner gave very good comparisons with model predictions in these counterflow diffusion flames. Quantitative LIF measurements of [NO] were also obtained in four methane-air counterflow partially-premixed flames with fuel-side equivalence ratios (phi(sub B)) of 1.45, 1.6, 1.8 and 2.0. The measurements were in excellent agreement with model predictions when accounting for radiative heat loss. Spatial separation between regions dominated by the prompt and thermal NO mechanisms was observed in the phi(sub B) = 1.45 flame. The modified rate coefficient proposed earlier for the prompt-NO initiation reaction improved agreement between code predictions and measurements in the region where prompt-NO dominates. Finally, LIF measurements of NO were obtained in counterflow diffusion flames at 2 to 5 atm. Comparisons between [NO] measurements and predictions show that the GRI mechanism underpredicts prompt-NO by a factor of two to three at all pressures. In general, the results indicate a need for refinement of the CH chemistry, especially the pressure-dependent CH formation and destruction reactions.

Ravikrishna, Rayavarapu V.↗
NASA NTRS1995

Energy Transfer in Rotating Turbulence

The influence or rotation on the spectral energy transfer of homogeneous turbulence is investigated in this paper. Given the fact that linear dynamics, e.g. the inertial waves regime tackled in an RDT (Rapid Distortion Theory) fashion, cannot Affect st homogeneous isotropic turbulent flow, the study of nonlinear dynamics is of prime importance in the case of rotating flows. Previous theoretical (including both weakly nonlinear and EDQNM theories), experimental and DNS (Direct Numerical Simulation) results are gathered here and compared in order to give a self-consistent picture of the nonlinear effects of rotation on tile turbulence. The inhibition of the energy cascade, which is linked to a reduction of the dissipation rate, is shown to be related to a damping due to rotation of the energy transfer. A model for this effect is quantified by a model equation for the derivative-skewness factor, which only involves a micro-Rossby number Ro(sup omega) = omega'/(2(OMEGA))-ratio of rms vorticity and background vorticity as the relevant rotation parameter, in accordance with DNS and EDQNM results fit addition, anisotropy is shown also to develop through nonlinear interactions modified by rotation, in an intermediate range of Rossby numbers (Ro(omega) = (omega)' and Ro(omega)w greater than 1), which is characterized by a marco-Rossby number Ro(sup L) less than 1 and Ro(omega) greater than 1 which is characterized by a macro-Rossby number based on an integral lengthscale L and the micro-Rossby number previously defined. This anisotropy is mainly an angular drain of spectral energy which tends to concentrate energy in tile wave-plane normal to the rotation axis, which is exactly both the slow and the two-dimensional manifold. In Addition, a polarization of the energy distribution in this slow 2D manifold enhances horizontal (normal to the rotation axis) velocity components, and underlies the anisotropic structure of the integral lengthscales. Finally is demonstrated the ability of a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) model to predict the underlying spectral transfer structure and all the subsequent developments of classic anisotropy indicators in physical space, when compared to recent LES results. Even if the applications mainly concern developed strong turbulence, a particular emphasis is given to the strong formal analogy of this EDQNM2 model with recent weakly nonlinear approaches to wave-turbulence.

Cambon, Claude↗
NASA NTRS2011

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗
NASA NTRS1956

On the Contribution of Turbulent Boundary Layers to the Noise Inside a Fuselage

The following report deals in preliminary fashion with the transmission through a fuselage of random noise generated on the fuselage skin by a turbulent boundary layer. The concept of attenuation is abandoned and instead the problem is formulated as a sequence of two linear couplings: the turbulent boundary layer fluctuations excite the fuselage skin in lateral vibrations and the skin vibrations induce sound inside the fuselage. The techniques used are those required to determine the response of linear systems to random forcing functions of several variables. A certain degree of idealization has been resorted to. Thus the boundary layer is assumed locally homogeneous, the fuselage skin is assumed flat, unlined and free from axial loads and the 'cabin' air is bounded only by the vibrating plate so that only outgoing waves are considered. Some of the details of the statistical description have been simplified in order to reveal the basic features of the problem. The results, strictly applicable only to the limiting case of thin boundary layers, show that the sound pressure intensity is proportional to the square of the free stream density, the square of cabin air density and inversely proportional to the first power of the damping constant and to the second power of the plate density. The dependence on free stream velocity and boundary layer thickness cannot be given in general without a detailed knowledge of the characteristics of the pressure fluctuations in the boundary layer (in particular the frequency spectrum). For a flat spectrum the noise intensity depends on the fifth power of the velocity and the first power of the boundary layer thickness. This suggests that boundary layer removal is probably not an economical means for decreasing cabin noise. In general, the analysis presented here only reduces the determination of cabin noise intensity to the measurement of the effect of any one of our variables (free stream velocity, boundary layer thickness, plate thickness or the characteristic velocity of propagation in the plate). The plate generates noise by vibrating in resonance over a wide range of frequencies and increasing the damping constant is consequently an effective method of decreasing noise generation. One of the main features of the results is that the relevant quantities upon which noise intensity depends are non-dimensional numbers in which boundary layer and plate properties enter as ratios. This is taken as an indication that in testing models of structures for boundary layer noise it is not sufficient to duplicate in the model the structural characteristics of the fuselage. One must match properly the characteristics of the exicitng pressure fluctuations to that of the structure.

Corcos, G. M.↗
NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗
NASA NTRS2014

Let's Exploit Available Knowledge on Vegetation Fluorescence

The potential to measure vegetation fluorescence from space (1) and to derive from it direct information on the gross primary productivity (GPP) of terrestrial ecosystems is probably the most thrilling development in remote sensing and global ecology of recent years, as it moves Earth observation techniques from the detection of canopy biophysics (e.g., fraction of absorbed radiation) and biochemistry (chlorophyll and nitrogen content) to the realm of ecosystem function. The existence of a functional relationship between fluorescence and photosynthesis has been elucidated over the last decade by several laboratories, notably as part of the preliminary studies of the European Space Agency Fluorescence Explorer (FLEX) Earth Explorer Mission. The empirical observation presented by Guanter et al. (2) of a linear relationship between fluorescence radiance and GPP, however, provides the first experimental confirmation of the feasibility of the approach— already thoroughly tested at leaf level—at the desired scale, despite the confounding effects associated with the satellite detection of such a faint signal. A word of clarification is needed here. The use of fluorescence as a probe of leaf photochemistry has been a staple of plant ecophysiology for decades, rooted in a sound understanding of photosynthetic energy dissipation. However, most past studies had to rely for the interpretation of results on active (pulse-saturated) techniques, making them unsuitable for remote-sensing applications. Over recent years, however, novel process based models have been developed for the interpretation of steady-state, solar-induced fluorescence at the leaf to canopy scale (3). We are therefore in a position to move beyond the mere empirical observation of an association between GPP and fluorescence radiance. In particular, Guanter et al. (2) base their analysis on the assumption of a constant ratio between photosynthetic and fluorescence light use efficiencies (equation 3 in ref. 2).We know, however, that the ratio is not constant, but changes widely in response to light, CO2, stomatal limitations, and extreme stress (4, 5). What’s more, we can make sense of such changes, thus extracting valuable information from the very scatter that is apparent in their data. However, this process will require the availability of more tailored instruments, such as the one planned for the FLEX mission. As already stressed by Guanter et al. (2), the spatial resolution of the Global Ozone Monitoring Experiment-2 sensor (40 × 80 km) makes it difficult to compare meaningfully the fluorescence signal with ground measurements, when only 60–70% of the footprint consists of the desired land-cover type (table S1 in ref. 2), suggesting that this could be largely responsible for the low signals observed in European grasslands. Moreover, the overpass time of the MetOp-A satellite (9:30 AM) implies that fluorescence is generally measured under light-limiting conditions, when fluorescence is only marginally affected by stomatal closure even under stress conditions. This result could explain the seasonal mismatch with daily GPP observed in natural ecosystems in the absence of irrigation (figure 4 in ref. 2). We hope, therefore, that this welcome contribution to this fast-advancing field will help demonstrate the potential of the new technique, and pave the way for more refined studies under both a technological and scientific point of view.

Space↗
NASA NTRS2018

Background and Lunar Neutron Populations Detected by LEND and Average Concentration of Near-Surface Hydrogen near the Moon's Poles

Neutron flux measurements by the Lunar Exploration Neutron Detector (LEND) on the Lunar Reconnaissance Orbiter (LRO) enable quantifying hydrogen-bearing volatiles in the lunar surface from orbit. Accurately determining hydrogen abundance requires discriminating between the instrument background detection rate and the population of lunar-sourced neutrons that are sensitive to surficial hydrogen. We have investigated the detection rate for lunar and non-lunar (spacecraft-sourced) neutrons in LEND by modeling maps of measured count rate in three LEND detector systems using linear combinations of maps compiled from LEND detectors and from the Lunar Prospector Neutron Spectrometer. We find that 30% of the global-average 24.926 0.020 neutron counts per second (cps) detected by the LEND STN3 thermal-energy neutron sensor are lunar-sourced neutrons in the thermal energy range (E < 0.4 eV), 65% are lunar-sourced neutrons in the epithermal and fast energy range (E > 0.4 eV), and 5% are from spacecraft-sourced background signal. In the SETN epithermal neutron detector, 90% of the 10.622 0.002 neutron detections per second are consistent with a lunar source of epithermal and fast neutrons combined (E > 0.4 eV), with 3% due to lunar-sourced thermal neutron leakage into the detector (E < 0.4 eV), and background signal accounting for 7% of total detections. Background signal due to spacecraft-derived neutrons is substantial in the CSETN collimated detector system, accounting for 57% of the global average detection rate of 5.082 0.001 cps, greater than the 48% estimated from cruise-phase data. Lunar-sourced epithermal and fast neutrons account for 43% of detected neutrons, including neutrons in collimation as well as neutrons that penetrate the collimator wall to reach the detector. We estimate a lower limit of 17% of lunar-sourced neutrons detected by CSETN are epithermal neutrons in collimation (0.37 cps), with an upper limit estimate of 54 11% of lunar-sourced neutrons received in collimation, or 1.2 0.2 cps global average. The pole-to-equator contrast ratio inepithermal and high-energy epithermal neutron flux indicates that the average concentration of hydrogen in the polar regolith above 80 north or south latitude is ~105 ppmw (parts per million by weight), or 0.095 0.01 wt% water-equivalent hydrogen. Above 88 north or south, the concentration increases to ~140 ppmw, or 0.13 0.02 wt% water-equivalent hydrogen. The similar pattern of neutron flux suppression at both poles suggests that hydrogen concentration generally increases nearer the pole and is not closely associated with a specific feature such as Shackleton Crater at the lunar south pole that has no northern counterpart. Epithermal neutron flux decreases with increasing latitude outside the polar regions, consistent with surface hydration that increases with latitude if that hydration extends to ~13-40 cm into the surface.

Livengood, Timothy A.↗
NASA NTRS2017

Background and Lunar Neutron Populations Detected by LEND and Average Concentration of Near-Surface Hydrogen near the Moon's Poles

Neutron flux measurements by the Lunar Exploration Neutron Detector (LEND) on the Lunar Reconnaissance Orbiter (LRO) enable quantifying hydrogen-bearing volatiles in the lunar surface from orbit. Accurately determining hydrogen abundance requires discriminating between the instrument background detection rate and the population of lunar-sourced neutrons that are sensitive to surficial hydrogen. We have investigated the detection rate for lunar and non-lunar (spacecraft-sourced) neutrons in LEND by modeling maps of measured count rate in three LEND detector systems using linear combinations of maps compiled from LEND detectors and from the Lunar Prospector Neutron Spectrometer. We find that 30% of the global-average 24.926 +/- 0.020 neutron counts per second (cps) detected by the LEND STN3 thermal-energy neutron sensor are lunar-sourced neutrons in the thermal energy range (E < 0.4 eV), 65% are lunar-sourced neutrons in the epithermal and fast energy range (E > 0.4 eV), and 5% are from spacecraft-sourced background signal. In the SETN epithermal neutron detector, 90% of the 10.622 +/- 0.002 neutron detections per second are consistent with a lunar source of epithermal and fast neutrons combined (E > 0.4 eV), with 3% due to lunar-sourced thermal neutron leakage into the detector (E < 0.4 eV), and background signal accounting for 7% of total detections. Background signal due to spacecraft-derived neutrons is substantial in the CSETN collimated detector system, accounting for 57% of the global average detection rate of 5.082 +/- 0.001 cps, greater than the 48% estimated from cruise-phase data. Lunar-sourced epithermal and fast neutrons account for 43% of detected neutrons, including neutrons in collimation as well as neutrons that penetrate the collimator wall to reach the detector. We estimate a lower limit of 17% of lunar-sourced neutrons detected by CSETN are epithermal neutrons in collimation (0.37 cps), with an upper limit estimate of 54 +/- 11% of lunar-sourced neutrons received in collimation, or 1.2 +/- 0.2 cps global average. The pole-to-equator contrast ratio in epithermal and high-energy epithermal neutron flux indicates that the average concentration of hydrogen in the polar regolith above 80deg north or south latitude is ∼105 ppmw (parts per million by weight), or 0.095 +/- 0.01 wt% water-equivalent hydrogen. Above 88deg north or south, the concentration increases to ∼140 ppmw, or 0.13 +/- 0.02 wt% water-equivalent hydrogen. The similar pattern of neutron flux suppression at both poles suggests that hydrogen concentration generally increases nearer the pole and is not closely associated with a specific feature such as Shackleton Crater at the lunar south pole that has no northern counterpart. Epithermal neutron flux decreases with increasing latitude outside the polar regions, consistent with surface hydration that increases with latitude if that hydration extends to ∼13-40 cm into the surface.

T A Livengood↗
NASA NTRS1992

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗
NASA NTRS2010

Convective Weather Forecast Accuracy Analysis at Center and Sector Levels

This paper presents a detailed convective forecast accuracy analysis at center and sector levels. The study is aimed to provide more meaningful forecast verification measures to aviation community, as well as to obtain useful information leading to the improvements in the weather translation capacity models. In general, the vast majority of forecast verification efforts over past decades have been on the calculation of traditional standard verification measure scores over forecast and observation data analyses onto grids. These verification measures based on the binary classification have been applied in quality assurance of weather forecast products at the national level for many years. Our research focuses on the forecast at the center and sector levels. We calculate the standard forecast verification measure scores for en-route air traffic centers and sectors first, followed by conducting the forecast validation analysis and related verification measures for weather intensities and locations at centers and sectors levels. An approach to improve the prediction of sector weather coverage by multiple sector forecasts is then developed. The weather severe intensity assessment was carried out by using the correlations between forecast and actual weather observation airspace coverage. The weather forecast accuracy on horizontal location was assessed by examining the forecast errors. The improvement in prediction of weather coverage was determined by the correlation between actual sector weather coverage and prediction. observed and forecasted Convective Weather Avoidance Model (CWAM) data collected from June to September in 2007. CWAM zero-minute forecast data with aircraft avoidance probability of 60% and 80% are used as the actual weather observation. All forecast measurements are based on 30-minute, 60- minute, 90-minute, and 120-minute forecasts with the same avoidance probabilities. The forecast accuracy analysis for times under one-hour showed that the errors in intensity and location for center forecast are relatively low. For example, 1-hour forecast intensity and horizontal location errors for ZDC center were about 0.12 and 0.13. However, the correlation between sector 1-hour forecast and actual weather coverage was weak, for sector ZDC32, about 32% of the total variation of observation weather intensity was unexplained by forecast; the sector horizontal location error was about 0.10. The paper also introduces an approach to estimate the sector three-dimensional actual weather coverage by using multiple sector forecasts, which turned out to produce better predictions. Using Multiple Linear Regression (MLR) model for this approach, the correlations between actual observation and the multiple sector forecast model prediction improved by several percents at 95% confidence level in comparison with single sector forecast.

Wang, Yao↗
NASA NTRS1994

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗
NASA NTRS2013

The Face that Launched a Thousand Slips

Helene, (approximately 17.6 kilometers mean radius) is an L4 Trojan co-orbital of Saturn's moon Dione. Its hemisphere features an unusual morphology consisting of broad depressions and a generally smooth surface patterned with streaks and grooves. The streaks appear to be oriented down-gradient, as are the grooves. This pattern suggests intensive mass-wasting as a dominant process on the leading hemisphere. Kilometer-scale impact craters are very sparse on the leading hemisphere other than the degraded kilometer-scale basins defining the overall satellite shape, and many small craters have a diffuse appearance suggesting ongoing mass wasting. Thus mass wasting must dominate surface-modifying processes at present. In fact, the mass wasting appears to have been sufficient in magnitude to narrow the divides between adjacent basins to narrow septa, similar, but in lower relief, to the honeycomb pattern of Hyperion. The prominent groves occur primarily near topographic divides and appear have cut into a broad, slightly lower albedo surface largely conforming to the present topography but elevated a few meters above the smooth surfaces undergoing mass wasting flow. Low ridges and albedo markings on the surface suggest surface flow of materials traveling up to several kilometers. Diffusive mass wasting produces smooth surfaces - such a pattern characterizes most of the low-lying surfaces. The grooves, however, imply that the transport process is advective at those locations where they occur, that is, erosion tends to concentrate along linear pathways separated by divides. In fact, in many places grooves have a fairly regular spacing of 125-160 meters, defining a characteristic erosional scale. Several questions are prompted by the unusual morphology of Helene: 1) What is the nature of the surface materials? 2) Are the transport processes gradual or catastrophic motion from one or a few events? 3) What mechanisms drive mass wasting and groove development? 4) Have the formative processes been active in the recent past? 5) Finally, is the surface accreting or eroding? The smooth character of the leading edge hemisphere of Helene and the dominance of mass wasting suggest that the surface is composed of fine-grained debris, probably dominated by dust-size to small gravel particles. The Lagrangian points of saturnian satellites are locations where planetesimals might have accreted to form co-orbital satellites such as and they may capture ejecta from their master moon. Published models suggest that Helene is a site of net accretion, but we find no extant explanation for the dominance of fine grain sizes for the surface (and probable subsurface) composition of Helene and the other Lagrangian satellites. Observation of the mass wasting tracks on Helene suggests the presence of well-defined streams of debris with low bordering levees that may be depositional features or remnants of the dissected higher surface. Some flows in grazing illumination appear to have a convex cross-section. This mass-flow morphology might be consistent with the occurrence of large-scale flow events, but which might have occurred through rapid emplacement or slow glacier-like creep. On the other hand, small craters appear to have been "softened" by creep-like processes, indicating ongoing mass wasting.

honeycomb pattern of Hyperion↗