Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random variables”

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 379 records · Page 21

Uncertainty Analysis of the CE-22 Advanced Nozzle Test Facility

This paper presents methods and results of a detailed measurement uncertainty analysis that was performed for the Advanced Nozzle Test Facility, CE-22, located at the NASA Glenn Research Center. Results for the uncertainty in thrust and flow coefficients, in addition to other variables of interest, are provided. Results are presented separately as random uncertainty (characterizing errors inherent to instrument or measurement environment, impacting repeatability), systematic uncertainty (capturing inaccuracies due to measurement process, calibration, installation effects or other sources which may introduce bias), and total combined uncertainty. The statistical methods and engineering judgments used to estimate elemental uncertainties are described. The analysis tool MANTUS (Measurement Analysis Tool for Uncertainty in Systems) was used to quantify instrumentation system uncertainty, and statistical analysis and engineering judgment were used to quantify other random and systematic uncertainty sources. The Monte Carlo method was used to propagate systematic and random elemental uncertainties to determine the uncertainties of various calculated variables of interest.

data quality↗

Combining remotely sensed and other measurements for hydrologic areal averages

A method is described for combining measurements of hydrologic variables of various sampling geometries and measurement accuracies to produce an estimated mean areal value over a watershed and a measure of the accuracy of the mean areal value. The method provides a means to integrate measurements from conventional hydrological networks and remote sensing. The resulting areal averages can be used to enhance a wide variety of hydrological applications including basin modeling. The correlation area method assigns weights to each available measurement (point, line, or areal) based on the area of the basin most accurately represented by the measurement. The statistical characteristics of the accuracy of the various measurement technologies and of the random fields of the hydrologic variables used in the study (water equivalent of the snow cover and soil moisture) required to implement the method are discussed.

Johnson, E. R.↗

Probabilistic and Possibilistic Analyses of the Strength of a Bonded Joint

The effects of uncertainties on the predicted strength of a single lap shear joint are examined. Probabilistic and possibilistic methods are used to account for uncertainties. A total of ten variables are assumed to be random, with normal distributions. Both Monte Carlo Simulation and the First Order Reliability Method are used to determine the probability of failure. Triangular membership functions with upper and lower bounds located at plus or minus three standard deviations are used to model uncertainty in the possibilistic analysis. The alpha cut (or vertex) method is used to evaluate the possibility of failure. Linear and geometrically nonlinear finite element analyses are used calculate the response of the joint; fracture in the adhesive and material strength failure in the strap are used to evaluate its strength. Although probabilistic and possibilistic analyses provide significantly more information than do conventional deterministic analyses, they are computationally expensive. A novel scaling approach is developed and used to substantially reduce the computational cost of the probabilistic and possibilistic analyses. The possibilistic approach for treating uncertainties appears to be viable during the conceptual and preliminary design stages when limited data are available and high accuracies are not needed. However, this viability is mixed with several cautions that are discussed herein.

Stroud, W. Jefferson↗

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.↗

Vegetation Monitoring Optimization with Normalized Difference Vegetation Index and Evapotranspiration Using Remote Sensing Measurements and Land Surface Models Over East Africa

The majority of people in East Africa rely on the agro-pastoral system for their livelihood, which is highly vulnerable to droughts and flooding. Agro-pastoral droughts are endemic to the region and are considered the main natural hazard that contributes to food insecurity. Drought begins with rainfall deficit, gradually leading to soil moisture deficit, higher land surface temperature, and finally impacts to vegetation growth. Therefore, monitoring vegetation conditions is essential in understanding the progression of drought, potential effects on food security, and providing early warning information needed for drought mitigation decisions. Because vegetation processes couple the land and atmosphere, monitoring of vegetation conditions requires consideration of both water provision and demand. While there is consensus in using either the Normalized Difference Vegetation Index (NDVI) or evapotranspiration(ET) for vegetation monitoring, a comprehensive assessment optimizing the use of both has not yet been done. Moreover, the evaluation methods for understanding the relationships between NDVI and ET for vegetation monitoring are also limited. Taking these gaps into account we have developed a framework to optimize vegetation monitoring using both NDVI and ET by identifying where they perform the best by using triple collocation and cross-correlation methods. We estimated the random error structure in Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI; ET from the Operational Simplified Surface Energy Balance (SSEBop) model; and ET from land surface models (LSMs). LSM ET and SSEBop ET have been found to be better indicators for vegetation monitoring during extreme drought events, while NDVI could provide better information on vegetation condition during wetter than normal conditions. The random error structures of these variables suggest that LSM ET is most likely to provide important information for vegetation monitoring over low and high ends of the vegetation fraction areas. Over moderate vegetative areas, any of these variables could provide important vegetation information for drought characterization and food security assessments. While this study provides a framework for optimizing vegetation monitoring for drought and food security assessments over East Africa, the framework can be adopted to optimize vegetation monitoring over any other drought and food insecure region of the world.

triple collocation↗

A probabilistic fatigue analysis of multiple site damage

The variability in initial crack size and fatigue crack growth is incorporated in a probabilistic model that is used to predict the fatigue lives for unstiffened aluminum alloy panels containing multiple site damage (MSD). The uncertainty of the damage in the MSD panel is represented by a distribution of fatigue crack lengths that are analytically derived from equivalent initial flaw sizes. The variability in fatigue crack growth rate is characterized by stochastic descriptions of crack growth parameters for a modified Paris crack growth law. A Monte-Carlo simulation explicitly describes the MSD panel by randomly selecting values from the stochastic variables and then grows the MSD cracks with a deterministic fatigue model until the panel fails. Different simulations investigate the influences of the fatigue variability on the distributions of remaining fatigue lives. Six cases that consider fixed and variable conditions of initial crack size and fatigue crack growth rate are examined. The crack size distribution exhibited a dominant effect on the remaining fatigue life distribution, and the variable crack growth rate exhibited a lesser effect on the distribution. In addition, the probabilistic model predicted that only a small percentage of the life remains after a lead crack develops in the MSD panel.

Rohrbaugh, S. M.↗

Gravitational Core-Mantle Coupling and the Acceleration of the Earth

Gravitational core-mantle coupling may be the cause of the observed variable acceleration of the Earth's rotation on the 1000 year timescale. The idea is that density inhomogeneities which randomly come and go in the liquid outer core gravitationally attract density inhomogeneities in the mantle and crust, torquing the mantle and changing its rotation state. The corresponding torque by the mantle on the core may also explain the westward drift of the magnetic field of 0.2 deg per year. Gravitational core-mantle coupling would stochastically affect the rate of change of the Earth's obliquity by just a few per cent. Its contribution to polar wander would only be about 0.5% the presently observed rate. Tidal friction is slowing down the rotation of the Earth, overwhelming a smaller positive acceleration from postglacial rebound. Coupling between the liquid outer core of the Earth and the mantle has long been a suspected reason for changes in the length-of-day. The present investigation focuses on the gravitational coupling between the density anomalies in the convecting liquid outer core and those in the mantle and crust as a possible cause for the observed nonsecular acceleration on the millenial timescale. The basic idea is as follows. There are density inhomogeneities caused by blobs circulating in the outer core like the blobs in a lava lamp; thus the outer core's gravitational field is not featureless. Moreover, these blobs will form and dissipate somewhat randomly. Thus there will be a time variability to the fields. These density inhomogeneities will gravitationally attract the density anomalies in the mantle.

Rubincam, David Parry↗

Colorado Front Range Disasters: Understanding the Impact of Forest Management on the Cameron Peak and CalWood Fire

Along the Colorado Front Range, forest management has gained significant attention due to uncharacteristically large fires that burned late in 2020. The Cameron Peak Fire (largest in Colorado recorded history) and the CalWood Fire collectively burned an estimated 219,019 acres from August through December of 2020. Project partners at the Coalition for the Poudre River Watershed, Colorado State Forest Service, Ben Delatour Scout Ranch, The Nature Conservancy, Colorado Forest Restoration Institute, and Colorado State University were interested in understanding the effectiveness of previous forest treatments in reducing burn severity within the Cameron Peak Fire and the CalWood Fire. We first collated a forest treatment dataset from pre-existing datasets by reclassifying over 29,000 treatments, which occurred across the Northern Colorado Front Range between 1970-2020. Secondly, we mapped three burn severity indices using Landsat 8 OLI and Sentinel-2 MSI Earth observations and compared them to soil burn severity field data. Thirdly, a total of 35 topographic, disturbance, forest structure, and treatment predictor variables were generated across the fires. Finally, we assessed relationships between these predictor variables and burn severity using the random forest algorithm. Model results indicate that the primary drivers of burn severity were elevation and distance to treatment edge for the Cameron Peak Fire and fire area and forest canopy cover for the CalWood Fire. Further analysis of these variables paired with field data is necessary to understand the relationship between burn severity and treatments to guide future restoration efforts, improve forest resiliency, and mitigate fire risks.

Neal Swayze↗

Forest Potential Productivity Mapping by Linking Remote-Sensing-Derived Metrics to Site Variables

A fine-resolution region-wide map of forest site productivity is an essential need for effective large-scale forestry planning and management. In this study, we incorporated Sentinel-2 satellite data into an increment-based measure of forest productivity (biomass growth index (BGI)) derived from climate, lithology, soils, and topographic metrics to map improved BGI (iBGI) in parts of North American Acadian regions. Initially, several Sentinel-2 variables including nine single spectral bands and 12 spectral vegetation indices (SVIs) were used in combination with forest management variables to predict tree volume/ha and height using Random Forest. The results showed a 10–12 % increase in out of bag (OOB) r2when Sentinel-2 variables were included in the prediction of both volume and height together with BGI. Later, selected Sentinel-2 variables were used for biomass growth prediction in Maine, USA and New Brunswick, Canada using data from 7738 provincial permanent sample plots. The Sentinel-2 red-edge position (S2REP) index was identified as the most important variable over others to have known influence on site productivity. While a slight improvement in the iBGI accuracy occurred compared to the base BGI model (~2%), substantial changes to coefficients of other variables were evident and some site variables became less important when S2REP was included.

site productivity↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit (CPU) on the order of ~30s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Wave theory for microwave remote sensing of a half-space random medium with three-dimensional variations

The two-variable expansion technique is used to solve for the mean Green's functions from the Dyson equation under the nonlinear approximation for a half-space random medium with three-dimensional correlation functions. The Bethe-Salpeter equations are then solved under the ladder approximation. The radiative transfer equations, which have been applied extensively in the study of microwave remote sensing problems, are derived under these approximations. The limiting cases of large and small horizontal correlation lengths are discussed. It is found that there is only one propagation constant except for the case of large horizontal correlation lengths, in which there are two propagation constants. It is shown that boundary layer appears in first-order solutions and does not appear in zeroth-order solutions.

Tsang, L.↗

Propagation in a shearing plasma. II - Turbulence and the frequency range of pulsar microstructure

Numerical calculations are presented which explore the propagation of radio waves and pulses through a turbulent shearing plasma. A stochastic model is used to describe the turbulent velocity field, in which random fluctuations are superposed on a mean profile. Both Gaussian wave packets and shot-noise pulses polarized in the plane of shearing acquire quasi-periodic intensity modulations above the rest-frame plasma frequency. The frequency range over which these modulations appear is shown to depend on the spatial correlation length of the turbulent fluctuations and on their standard deviation from the mean velocity. The period of the modulations is variable from one realization of the random process to the next, and the average period is frequency dependent. The results of these calculations lend further support to a model for pulsar microstructure in which periodic micropulses are pure temporal modulations of the emitted radiation due to propagation effects in magnetospheric shearing regions. It is shown that turbulence on a scale of approximately 10 cm can produce these quasi-periodic modulations over a frequency range comparable to that of observed microstructure.

Harding, A. K.↗

A global reference atmospheric model for surface to orbital altitudes

An empirical atmospheric model has been developed which generates values for pressure, density, temperature and winds from surface levels to orbital altitudes. The output parameters consist of components for: (1) latitude, longitude, and altitude dependent monthly means; (2) quasi-biennial oscillations; and (3) random perturbations to partially simulate the variability due to synoptic, diurnal, planetary wave and gravity wave variations. The monthly mean models consist of: NASA's four dimensional worldwide model for height, latitude, and longitude dependent monthly means from the surface to 25 km; and a newly developed latitude-longitude dependent model which is an extension of the Groves latitude dependent model for the region between 25 and 90 km. The Jacchia 1970 model is used above 90 km and is faired with the modified Groves values between 90 and 115 km. Quasi-biennial and random variation perturbations are computed from parameters determined from various empirical studies, and are added to the monthly mean values.

Justus, C. G.↗

Computational simulation of probabilistic lifetime strength for aerospace materials subjected to high temperature, mechanical fatigue, creep, and thermal fatigue

The results of a fourth year effort of a research program conducted for NASA-LeRC by The University of Texas at San Antonio (UTSA) are presented. The research included on-going development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic material strength degradation model, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subjected to a number of effects or primitive variables. These primitive variables may include high temperature, fatigue, or creep. In most cases, strength is reduced as a result of the action of a variable. This multifactor interaction strength degradation equation was randomized and is included in the computer program, PROMISC. Also included in the research is the development of methodology to calibrate the above-described constitutive equation using actual experimental materials data together with regression analysis of that data, thereby predicting values for the empirical material constants for each effect or primitive variable. This regression methodology is included in the computer program, PROMISC. Actual experimental materials data were obtained from industry and the open literature for materials typically for applications in aerospace propulsion system components. Material data for Inconel 718 was analyzed using the developed methodology.

Boyce, Lola↗

Quantum Adiabatic Optimization and Combinatorial Landscapes

In this paper we analyze the performance of the Quantum Adiabatic Evolution (QAE) algorithm on a variant of Satisfiability problem for an ensemble of random graphs parametrized by the ratio of clauses to variables, gamma = M / N. We introduce a set of macroscopic parameters (landscapes) and put forward an ansatz of universality for random bit flips. We then formulate the problem of finding the smallest eigenvalue and the excitation gap as a statistical mechanics problem. We use the so-called annealing approximation with a refinement that a finite set of macroscopic variables (verses only energy) is used, and are able to show the existence of a dynamic threshold gamma = gammad, beyond which QAE should take an exponentially long time to find a solution. We compare the results for extended and simplified sets of landscapes and provide numerical evidence in support of our universality ansatz.

Smelyanskiy, V. N.↗