Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

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,225 records · Page 68

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗

Isolating the Surface Type Influence on Arctic Cloud Properties

In response to anthropogenic climate change, substantial declines in sea ice have been observed. An urgent question is whether and how clouds respond to the changing Arctic surface. Due to the important radiative influence of clouds, the response of clouds to sea ice cover change constitutes a potentially important climate feedback within the Arctic. The sign of this feedback is strongly influenced by the seasonality of the cloud response to sea ice loss. Given the importance of this phenomenon, previous research focused on using observations of inter-annual variability to quantify the cloud-sea ice relationship. However, the reliance on inter-annual variability to assess the covariance between clouds and sea ice make it challenging to control for the influence of large-scale meteorology on the clouds. We have devised an approach that applies an in situ, airborne observational strategy to quantify the influence of surface type on Arctic clouds to satellite data. Our event-based method composites cloud property differences for ocean, marginal ice zone, and sea ice surface types for >5000 MIZ crossing events. We find that cloud properties in non-summer months below ~1.5 km are influenced by surface type such that greater cloud fraction and water content are found over ocean relative to sea ice regions. Results are statistically significantly at the 95% confidence level. During summer, the surface type does not influence the cloud properties. The primary cause of the surface type dependent cloud property differences is the thermodynamic profiles differences that occur between the two surface types: namely, ocean footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their sea ice counterparts. We find that the differences in thermodynamic profiles by surface type that correspond to the cloud property differences also explains a significant portion of the variability across events. Our analysis provides evidence that the influence of surface type on cloud properties fundamentally occurs due to the inherent surface temperature differences between the surface types. Significant surface type influences on cloud properties are only found in the presence of surface temperature differences between the surface types and only non-summer months meet this criterion. A conceptual model of this process is as follows, fundamental properties of the ocean and sea ice surface types (e.g., surface albedo, heat capacity/thermal inertia, surface roughness, and thermal conductance) lead to systematic differences in the surface temperature by surface type. In response to this surface temperature difference, atmosphere boundary layer processes (including turbulent mixing and radiation) communicate the surface temperature differences to the lower atmosphere driving the greater stability over sea ice than ocean. These differences in the background thermodynamic state profiles between the surface types yield atmospheric conditions that are more conducive to cloud development over ocean than sea ice. The details of these boundary layer processes are sensitive to the surface type but may not be fundamental to the resulting thermodynamic structure and cloud differences. Thus, our results indicate that to accurately model that cloud response to sea ice loss in climate models the surface type properties that yield the systematic differences in ocean and sea ice temperature (e.g., surface albedo, surface turbulent fluxes parameters, and thermal conductance) must be accurately represented.

Arctic radiation budget↗

Optimal Estimation Framework for Ocean Color Atmospheric Correction and Pixel-level Uncertainty Quantification

Ocean color remote sensing requires compensation for atmospheric scattering and absorption (aerosol, Rayleigh, and trace gases), referred to as atmospheric correction (AC). AC allows inference of parameters such as spectrally resolved remote sensing reflectance ( R rs )(λ) ; sr 1 ) at the ocean surface from the top-of-atmosphere reflectance. Often, the uncertainty of this process is not fully explored. Bayesian inference techniques provide a simultaneous AC and uncertainty assessment via a full posterior distribution of the relevant variables, given the prior distribution of those variables and the radiative transfer (RT) likelihood function. Given uncertainties in the algorithm inputs, the Bayesian framework enables better constraints on the AC process by using the complete spectral information compared to traditional approaches that use only a subset of bands for AC. This paper investigates a Bayesian inference research method (Optimal Estimation, OE) for ocean color AC by simultaneously retrieving atmospheric and ocean properties using all visible and near-infrared spectral bands. The OE algorithm analytically approximates the posterior distribution of parameters based on normality assumptions and provides a potentially viable operational algorithm with a reduced computational expense. We developed a Neural Network (NN) RT forward model look-up-table-based emulator to increase algorithm efficiency further and thus speed up the likelihood computations. We then applied the OE algorithm to synthetic data and observations from the MODerate resolution Imaging Spectroradiometer (MODIS) on NASA’s Aqua spacecraft. We compared the R rs )(λ) retrieval and its uncertainty estimates from the OE method with in-situ validation data from the SeaWiFS Bio-optical Archive and Storage System (SeaBASS) and Aerosol Robotic Network Ocean Color (AERONET-OC) datasets. The OE algorithm improved R rs )(λ) estimates relative to the NASA standard operational algorithm by improving all statistical metrics at 443, 555, and 667 nm. Unphysical negative R rs )(λ) , which often appear in complex water conditions, was reduced by a factor of 3. The OE-derived pixel-level R rs )(λ) uncertainty estimates were also assessed relative to in-situ data and were shown to have skill.

Atmospheric correction↗

Kinetic model of turbulence in an incompressible fluid

A statistical description of turbulence in an incompressible fluid obeying the Navier-Stokes equations is proposed, where pressure is regarded as a potential for the interaction between fluid elements. A scaling procedure divides a fluctuation into three ranks representing the three transport processes of macroscopic evolution, transport property, and relaxation. Closure is obtained by relaxation, and a kinetic equation is obtained for the fluctuation of the macroscopic rank of the distribution function. The solution gives the transfer function and eddy viscosity. When applied to the inertia subrange of the energy spectrum the analysis recovers the Kolmogorov law and its numerical coefficient.

Tchen, C. M.↗

Psychophysiological investigations of the biomedical problems of manned spaceflight

The Final Report on psychophysiological investigations of the biomedical problems of manned spaceflight is presented. In the first project, statistical analyses of human autonomic data were performed. The objectives were the following: to establish a relational data base containing human psychophysiological data obtained from Shuttle flight experiments and ground-based research over a 20 year period; to enable multi-user access and retrieval of these data for subsequent analyses and for possible inclusion in the proposed Life Sciences Data Archive; and to enable/conduct statistical analyses across several experiments on large subject populations which can thereby provide definitive answers to questions on human autonomic and behavioral responses and adaptation to environmental stressors on Earth and in space. The second project studied motion sickness. The objectives were: to test/develop hardware and procedures to be incorporated into preflight training of crewmembers of Spacelab-J; and to examine spin-off applications of AFT. The third project studied orthostatic intolerance. The objective was to test the feasibility of applying autogenic-feedback training as a potential treatment for postflight orthostatic intolerance.

Source record↗

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗

Photoionization Modeling and the K Lines of Iron

We calculate the efficiency of iron K line emission and iron K absorption in photoionized models using a new set of atomic data. These data are more comprehensive than those previously applied to the modeling of iron K lines from photoionized gases, and allow us to systematically examine the behavior of the properties of line emission and absorption as a function of the ionization parameter, density and column density of model constant density clouds. We show that, for example, the net fluorescence yield for the highly charged ions is sensitive to the level population distribution produced by photoionization, and these yields are generally smaller than those predicted assuming the population is according to statistical weight. We demonstrate that the effects of the many strongly damped resonances below the K ionization thresholds conspire to smear the edge, thereby potentially affecting the astrophysical interpretation of absorption features in the 7-9 keV energy band. We show that the centroid of the ensemble of K(alpha) lines, the K(beta) energy, and the ratio of the K(alpha(sub 1)) to K(alpha(sub 2)) components are all diagnostics of the ionization parameter of our model slabs.

Kallman, T. R.↗

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Morin, Lee↗

Space Transportation System Liftoff Debris Mitigation Process Overview

Liftoff debris is a top risk to the Space Shuttle Vehicle. To manage the Liftoff debris risk, the Space Shuttle Program created a team with in the Propulsion Systems Engineering & Integration Office. The Shutt le Liftoff Debris Team harnesses the Systems Engineering process to i dentify, assess, mitigate, and communicate the Liftoff debris risk. T he Liftoff Debris Team leverages off the technical knowledge and expe rtise of engineering groups across multiple NASA centers to integrate total system solutions. These solutions connect the hardware and ana lyses to identify and characterize debris sources and zones contribut ing to the Liftoff debris risk. The solutions incorporate analyses sp anning: the definition and modeling of natural and induced environmen ts; material characterizations; statistical trending analyses, imager y based trajectory analyses; debris transport analyses, and risk asse ssments. The verification and validation of these analyses are bound by conservative assumptions and anchored by testing and flight data. The Liftoff debris risk mitigation is managed through vigilant collab orative work between the Liftoff Debris Team and Launch Pad Operation s personnel and through the management of requirements, interfaces, r isk documentation, configurations, and technical data. Furthermore, o n day of launch, decision analysis is used to apply the wealth of ana lyses to case specific identified risks. This presentation describes how the Liftoff Debris Team applies Systems Engineering in their proce sses to mitigate risk and improve the safety of the Space Shuttle Veh icle.

Mitchell, Michael↗

Low-level liquid cloud properties during ORACLES retrieved using airborne polarimetric measurements and a neural network algorithm

In this study we developed a neural network (NN) that can be used to retrieve cloud microphysical properties from multiangular and multispectral polarimetric remote sensing observations. This effort builds upon our previous work, which explored the sensitivity of neural network input, architecture, and other design requirements for this type of remote sensing problem. In particular this work introduces a framework for appropriately weighting total and polarized reflectances, which have vastly different magnitudes and measurement uncertainties. The NN is trained using an artificial training set and applied to research scanning polarimeter (RSP) data obtained during the ORACLES field campaign (ObseRvations of Aerosols above CLouds and their intEractionS). The polarimetric RSP observations are unique in that they observe the same cloud from a very large number of angles within a variety of spectral bands, resulting in a large dataset that can be explored rapidly with a NN approach. The usefulness of applying a NN to a dataset such as this one stems from the possibility of rapidly obtaining a retrieval that could be subsequently applied as a first guess for slower but more rigorous physical-based retrieval algorithms. This approach could be particularly advantageous for more complicated atmospheric retrievals – such as when an aerosol layer lies above clouds like in ORACLES. For RSP observations obtained during ORACLES 2016, comparisons between the NN and standard parametric polarimetric (PP) cloud retrieval give reasonable results for droplet effective radius (r(e): R=0.756, RMSE=1.74 µm) and cloud optical thickness (τ: R=0.950, RMSE=1.82). This level of statistical agreement is shown to be similar to comparisons between the two most well-established cloud retrievals, namely, the polarimetric and the bispectral total reflectance cloud retrievals. The NN retrievals from the ORACLES 2017 dataset result in retrievals of r(e) (R=0.54, RMSE=4.77 µm) and τ (R=0.785, RMSE=5.61) that behave much more poorly. In particular we found that our NN retrieval approach does not perform well for thin (τ<3), inhomogeneous, or broken clouds. We also found that correction for above-cloud atmospheric absorption improved the NN retrievals moderately – but retrievals without this correction still behaved similarly to existing cloud retrievals with a slight systematic offset.

SPECTRAL-RESOLUTION LIDAR↗

Analysis of the effect of numbers of aircraft operations on community annoyance

The general validity of the equivalent-energy concept as applied to community annoyance to aircraft noise has been recently questioned by investigators using a peak-dBA concept. Using data previously gathered around nine U.S. airports, empirical tests of both concepts are presented. Results show that annoyance response follows neither concept, that annoyance increases steadily with energy-mean level for constant daily operations and with numbers of operations up to 100-199 per day (then decreases for higher numbers), and that the behavior of certain response descriptors is dependent upon the statistical distributions of numbers and levels.

Connor, W. K.↗

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗

Spectral characteristics and the extent of paleosols of the Palouse formation

Thematic mapping data was analyzed and verified by comparison to previously gathered transect samples and to aerial photographs. A bare-soil field with exposed paleosols characterized by slight enrichment of iron was investigated. Spectral relationships were first investigated statistically by creating a data set with DN values spatially matched as nearly as possible to field sample points. Chemical data for each point included organic carbon, free iron oxide, and amorphous iron content. The chemical data, DN values, and various band ratios were examined with the program package Statistix in order to find the combinations of reflectance data most likely to show a relationship which would dependably separate the exposed paleosols from the other soils. Cluster analysis and Fastclas classification procedures were applied to the most promising of the band ratio combinations.

Frazier, B. E.↗

Numerical simulation of two-dimensional spatially-developing mixing layers

Two-dimensional, incompressible, spatially developing mixing layer simulations are performed at Re = 10(exp 2) and 10(exp 4) with two classes of perturbations applied at the inlet boundary; combinations of discrete modes from linear stability theory, and a broad spectrum of modes derived from experimentally measured velocity spectra. The effect of the type and strength of inlet perturbations on vortex dynamics and time-averaged properties are explored. Two-point spatial velocity and autocorrelations are used to estimate the size and lifetime of the resulting coherent structures and to explore possible feedback effects. The computed time-averaged properties such as mean velocity profiles, turbulent statistics, and spread rates show good agreement with experimentally measured values. It is shown that by forcing with a broad spectrum of modes derived from an experimental energy spectrum many experimentally observed phenomena can be reproduced by a 2-D simulation. The strength of the forcing merely affected the length required for the dominant coherent structures to become fully-developed. Thus intensities comparable to those of the background turbulence in many wind tunnel experiments produced the same results, given sufficient simulation length.

Wilson, R. V.↗

Soft Decision Analyzer

The Soft Decision Analyzer (SDA) is an instrument that combines hardware, firmware, and software to perform realtime closed-loop end-to-end statistical analysis of single- or dual- channel serial digital RF communications systems operating in very low signal-to-noise conditions. As an innovation, the unique SDA capabilities allow it to perform analysis of situations where the receiving communication system slips bits due to low signal-to-noise conditions or experiences constellation rotations resulting in channel polarity in versions or channel assignment swaps. SDA s closed-loop detection allows it to instrument a live system and correlate observations with frame, codeword, and packet losses, as well as Quality of Service (QoS) and Quality of Experience (QoE) events. The SDA s abilities are not confined to performing analysis in low signal-to-noise conditions. Its analysis provides in-depth insight of a communication system s receiver performance in a variety of operating conditions. The SDA incorporates two techniques for identifying slips. The first is an examination of content of the received data stream s relation to the transmitted data content and the second is a direct examination of the receiver s recovered clock signals relative to a reference. Both techniques provide benefits in different ways and allow the communication engineer evaluating test results increased confidence and understanding of receiver performance. Direct examination of data contents is performed by two different data techniques, power correlation or a modified Massey correlation, and can be applied to soft decision data widths 1 to 12 bits wide over a correlation depth ranging from 16 to 512 samples. The SDA detects receiver bit slips within a 4 bits window and can handle systems with up to four quadrants (QPSK, SQPSK, and BPSK systems). The SDA continuously monitors correlation results to characterize slips and quadrant change and is capable of performing analysis even when the receiver under test is subjected to conditions where its performance degrades to high error rates (30 percent or beyond). The design incorporates a number of features, such as watchdog triggers that permit the SDA system to recover from large receiver upsets automatically and continue accumulating performance analysis unaided by operator intervention. This accommodates tests that can last in the order of days in order to gain statistical confidence in results and is also useful for capturing snapshots of rare events.

Steele, Glen↗

Regional Land Use Mapping: the Phoenix Pilot Project

The Phoenix Pilot Program has been designed to make effective use of past experience in making land use maps and collecting land use information. Conclusions reached from the project are: (1) Land use maps and accompanying statistical information of reasonable accuracy and quality can be compiled at a scale of 1:250,000 from orbital imagery. (2) Orbital imagery used in conjunction with other sources of information when available can significantly enhance the collection and analysis of land use information. (3) Orbital imagery combined with modern computer technology will help resolve the problem of obtaining land use data quickly and on a regular basis, which will greatly enhance the usefulness of such data in regional planning, land management, and other applied programs. (4) Agreement on a framework or scheme of land use classification for use with orbital imagery will be necessary for effective use of land use data.

Anderson, J. R.↗

Venus global radar reflectivity and correlations with elevation

A geological analysis of the global radar reflectivity data for Venus is presented. The statistical and map pattern distribution of the radar reflectivity rho is first analyzed in the context of the porosity structure of the surface. Rho is then correlated with elevation in order to examine the degree of topographic control on material properties which could depend on PT-sensitive chemical weathering phenomena. Comparison of the rms slope alpha and rho correlations with elevation demonstrates that there are distinct topographic zones defined on the basis of the rho and alpha trends. Finally, hierarchical clustering analysis is applied to the global rho, alpha, and 100-km baseline regional slope data sets in order to define radar property-related subregions.

Garvin, J. B.↗