Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “standard normal variate”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large streamwise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.↗
DOE OSTI2020

Generating realistic building electrical load profiles through the Generative Adversarial Network (GAN)

Building electrical load profiles can improve understanding of building energy efficiency, demand flexibility, and building-grid interactions. Current approaches to generating load profiles are time-consuming and not capable of reflecting the dynamic and stochastic behaviors of real buildings; some approaches also trigger data privacy concerns. In this study, we proposed a novel approach for generating realistic electrical load profiles of buildings through the Generative Adversarial Network (GAN), a machine learning technique that is capable of revealing an unknown probability distribution purely from data. The proposed approach has three main steps: (1) normalizing the daily 24-hour load profiles, (2) clustering the daily load profiles with the k-means algorithm, and (3) using GAN to generate daily load profiles for each cluster. The approach was tested with an open-source database – the Building Data Genome Project. We validated the proposed method by comparing the mean, standard deviation, and distribution of key parameters of the generated load profiles with those of the real ones. The KL divergence of the generated and real load profiles are within 0.3 for majority of parameters and clusters. Additionally, results showed the load profiles generated by GAN can capture not only the general trend but also the random variations of the actual electrical loads in buildings. We report the proposed GAN approach can be used to generate building electrical load profiles, verify other load profile generation models, detect changes to load profiles, and more importantly, anonymize smart meter data for sharing, to support research and applications of grid-interactive efficient buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2019

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large stream wise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.↗
NASA NTRS1973

On the normality and accuracy of simulated random processes.

Efficient methods of simulating stationary and nonstationary random processes and envelopes, by using a series of sine or cosine functions or by using the fast Fourier transform, have been proposed previously. Without applying the central limit theorem, it is shown in this paper that the simulated random processes are asymptotically Gaussian processes as the number of terms, N, of sine or cosine functions approaches infinity. The accuracy of the first-order probability densities of the simulated random processes is investigated by using the fast Fourier transform. Numerical results are computed with respect to the variation of the number of terms, N, of sine or cosine functions used for simulation. It is shown that within the practical range of N, such as 500, the accuracy is remarkably satisfactory even outside the region of 3 standard deviations.

Yang, J.-N.↗
NASA NTRS2019

Quantifying the Storm Time Thermospheric Neutral Density Variations Using Model and Observations

Accurate determination of thermospheric neutral density holds crucial importance for satellited rag calculations. The problem is twofold and involves the correct estimation of the quiet time climatology and storm time variations. In this work, neutral density estimations from two empirical and threephysics based models of the ionosphere thermosphere are compared with the neutral densities along the Challenging MicroSatellite Payload satellite track for six geomagnetic storms. Storm time variations are extracted from neutral density by (1) subtracting the mean difference between model and observation (bias),(2) setting climatological variations to zero, and (3) multiplying model data with the quiet time ratio between the model and observation. Several metrics are employed to evaluate the model performances. We find that the removal of bias or climatology reveals actual performance of the model in simulating the storm time variations. When bias is removed, depending on event and model, storm time errors in neutral density can decrease by an amount of 113% or can increase by an amount of 12% with respect to error in models with quiet time bias. It is shown that using only average and maximum values of neutral density to determine the model performances can be misleading since a model can estimate the averages fairly well but may not capture the maximum value or vice versa. Since each of the metrics used for determining model performances provides different aspects of the error, among these, we suggest employing mean absolute error, prediction efficiency, and normalized root mean square error together as a standard set ofmetrics for the neutral density.

Eyiguler, E. Ceren Kalafatoglu↗
DOE OSTI2024

Forest carbon sequestration on the west coast, USA: Role of species, productivity, and stockability

Forest ecosystems store large amounts of carbon and can be important sources, or sinks, of the atmospheric carbon dioxide that is contributing to global warming. Understanding the carbon storage potential of different forests and their response to management and disturbance events are fundamental to developing policies and scenarios to partially offset greenhouse gas emissions. Projections of live tree carbon accumulation are handled differently in different models, with inconsistent results. We developed growth-and-yield style models to predict stand-level live tree carbon density as a function of stand age in all vegetation types of the coastal Pacific region, US (California, Oregon, and Washington), from 7,523 national forest inventory plots. We incorporated site productivity and stockability within the Chapman-Richards equation and tested whether intensively managed private forests behaved differently from less managed public forests. We found that the best models incorporated stockability in the equation term controlling stand carrying capacity, and site productivity in the equation terms controlling the growth rate and shape of the curve. RMSEs ranged from 10 to 137 Mg C/ha for different vegetation types. There was not a significant effect of ownership over the standard industrial rotation length (~50 yrs) for the productive Douglas-fir/western hemlock zone, indicating that differences in stockability and productivity captured much of the variation attributed to management intensity. Our models suggest that doubling the rotation length on these intensively managed lands from 35 to 70 years would result in 2.35 times more live tree carbon stored on the landscape. These findings are at odds with some studies that have projected higher carbon densities with stand age for the same vegetation types, and have not found an increase in yields (on an annual basis) with longer rotations. We suspect that differences are primarily due to the application of yield curves developed from fully-stocked, undisturbed, single-species, “normal” stands without accounting for the substantial proportion of forests that don’t meet those assumptions. The carbon accumulation curves developed here can be applied directly in growth-and-yield style projection models, and used to validate the predictions of ecophysiological, cohort, or single-tree style models being used to project carbon futures for forests in the region. Our approach may prove useful for developing robust models in other forest types.

Chisholm, Paul J. (ORCID:0000000238784707)↗
NASA NTRS2010

Global Patterns of Lightning Properties Derived by LIS

The Lightning Imaging Sensor LIS aboard the TRMM satellite provides unmatched empirical data of the global lightning distribution (up to approx.35deg S/N) since end of 1997. Climatological flash rate densities derived from LIS are standard references, e.g. for flash rate parameterizations used in GCMs. It is known that flash characteristics are quite variable, and that various quantities (like the flash energy or the NOx production per flash) vary considerably, statistically as well as systematically on regional and seasonal scales. LIS provides information beyond flash counts, in particular radiance and flash footprint. Here we present an analysis of global patterns of various lightning properties derived from LIS, in relation to the number of flashes. These normalized flash characteristics show consistent spatial patterns of regions with "strong" versus regions with "weak" lightning. Most striking is a clear land-ocean contrast, with oceanic flashes being "stronger" than continental flashes. But also over continents, flash strength shows systematic variations. Highest continental values are found over the US, while values over South America and India are quite low. These regional variations cannot be simply parameterized as function of latitude. Information on spatial patterns of mean flash "strength", though rather qualitative up to now, is potentially a valuable input for improving empirical parameterizations based on flash counts (like precipitation or lightning NOx). Further investigation is in progress to come to a more physical and quantitative understanding of the spatial patterns of the different LIS properties. In particular, it has to be checked how far they could be related to established lightning properties (like energy or the fraction of intra-cloud to cloud-to-ground flashes) or to meteorological quantities (like CAPE).

Beirle, S.↗
NASA NTRS2010

Development of an Uncertainty Model for the National Transonic Facility

This paper introduces an uncertainty model being developed for the National Transonic Facility (NTF). The model uses a Monte Carlo technique to propagate standard uncertainties of measured values through the NTF data reduction equations to calculate the combined uncertainties of the key aerodynamic force and moment coefficients and freestream properties. The uncertainty propagation approach to assessing data variability is compared with ongoing data quality assessment activities at the NTF, notably check standard testing using statistical process control (SPC) techniques. It is shown that the two approaches are complementary and both are necessary tools for data quality assessment and improvement activities. The SPC approach is the final arbiter of variability in a facility. Its result encompasses variation due to people, processes, test equipment, and test article. The uncertainty propagation approach is limited mainly to the data reduction process. However, it is useful because it helps to assess the causes of variability seen in the data and consequently provides a basis for improvement. For example, it is shown that Mach number random uncertainty is dominated by static pressure variation over most of the dynamic pressure range tested. However, the random uncertainty in the drag coefficient is generally dominated by axial and normal force uncertainty with much less contribution from freestream conditions.

Walter, Joel A.↗
NASA NTRS1997

Experimental Investigation of Unsteady Shock Wave Turbulent Boundary Layer Interactions About a Blunt Fin

A series of experiments were performed to investigate the effects of Mach number variation on the characteristics of the unsteady shock wave/turbulent boundary layer interaction generated by a blunt fin. A single blunt fin hemicylindrical leading edge diameter size was used in all of the experiments which covered the Mach number range from 2.0 to 5.0. The measurements in this investigation included surface flow visualization, static and dynamic pressure measurements, both on centerline and off-centerline of the blunt fin axis. Surface flow visualization and static pressure measurements showed that the spatial extent of the shock wave/turbulent boundary layer interaction increased with increasing Mach number. The maximum static pressure, normalized by the incoming static pressure, measured at the peak location in the separated flow region ahead of the blunt fin was found to increase with increasing Mach number. The mean and standard deviations of the fluctuating pressure signals from the dynamic pressure transducers were found to collapse to self-similar distributions as a function of the distance perpendicular to the separation line. The standard deviation of the pressure signals showed initial peaked distribution, with the maximum standard deviation point corresponding to the location of the separation line at Mach number 3.0 to 5.0. At Mach 2.0 the maximum standard deviation point was found to occur significantly upstream of the separation line. The intermittency distributions of the separation shock wave motion were found to be self-similar profiles for all Mach numbers. The intermittent region length was found to increase with Mach number and decrease with interaction sweepback angle. For Mach numbers 3.0 to 5.0 the separation line was found to correspond to high intermittencies or equivalently to the downstream locus of the separation shock wave motion. The Mach 2.0 tests, however, showed that the intermittent region occurs significantly upstream of the separation line. Power spectral densities measured in the intermittent regions were found to have self-similar frequency distributions when compared as functions of a Strouhal number for all Mach numbers and interaction sweepback angles. The maximum zero-crossing frequencies were found to correspond with the peak frequencies in the power spectra measured in the intermittent region.

Barnhart, Paul J.↗
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear-and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute time-scale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

irradiance↗
DOE OSTI2025

Thermal oxidation of nuclear graphite and pyrolytic carbon coatings

The oxidation of pyrolytic carbon (PyC) deposited via fluidized bed chemical vapor deposition was characterized and compared with that of standard nuclear-grade graphite. The materials were heated at 700 to 1000 °C in a thermogravimetric analysis system under 20% v/v O 2 flow, allowing for direct comparison of dynamic oxidative mass change in each material. Further, three different PyC samples fabricated under different conditions exhibited variation in total mass loss and mass loss rate, varying by as much as 709 mg/cm 2 in total mass loss and 14.2 (mg/cm 2 )/min in mass loss rate at a single temperature. These variations highlight the correlation between PyC microstructure/defect density and oxidation susceptibility. Additionally, changes in the microstructure and composition between PyC and graphite were characterized via scanning electron microscopy and correlated to the mass loss results. The results of this work have implications toward the safety of tristructural isotropic (TRISO) and other coated particle fuels, especially under off-normal conditions, given the limited information that exists about the oxidation behavior of PyC.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Variability of prediction of maximal oxygen concumption on the cycle ergometer using standard equations

Several investigations within the Exercise Countermeasures Project at the NASA Johnson Space Center focused on the assessment of maximum oxygen consumption (VO2(sub max)) within the Astronaut Corps pre- and postspace flight. Investigations during the Apollo era suggested that there was a significant decrease in postflight VO2(sub max) when compared to preflight values, and current studies have documented that this trend continues in the Space Shuttle era. It is generally accepted and was confirmed in our laboratory that VO2(sub max) can be predicted from submaximal measures taken during graded exercise tests on the cycle ergometer with respect to populations. However, previous work had not examined the effect of day-to-day variations in the physiologic responses that might alter these predictions for individuals. Stability of individual submaximal data over serial tests is important so that predicted changes in VO2(sub max) are reflective of actual VO2(sub max) changes. Therefore, the purpose of this investigation was to determine which of the accepted equations to predict VO2(sub max) would be less affected by normal daily physiologic changes.

Greenisen, Michael C.↗
DOE OSTI2025

Analyzing the Quantum Approximate Optimization Algorithm: Ansätze, Symmetries, and Lie Algebras

The quantum approximate optimization algorithm (QAOA) has been proposed as a method to obtain approximate solutions for combinatorial optimization tasks. In this work, we study the underlying algebraic properties of three QAOA ansätze for the maximum-cut problem on connected graphs, while focusing on the generated Lie algebras as well as their invariant subspaces. Specifically, we analyze the standard QAOA ansatz as well as the orbit and multiangle ansätze. We are able to fully characterize the Lie algebras of the multiangle ansatz across arbitrary connected graphs, finding that they only fall into one of just six families. Aside from the cycle and path graphs, the Lie dimensions for every graph are exponentially large in the system size, meaning that multiangle ansätze are extremely prone to exhibiting barren plateaus. Then, a similar quasi-graph-independent Lie-algebraic characterization beyond the multiangle ansatz is impeded as the circuit exhibits additional “hidden” symmetries besides those naturally arising from a certain parity-superselection operator and all automorphisms of the considered graph. Disregarding the “hidden” symmetries, we can upper bound the dimensions of the orbit and the standard Lie algebras, and the dimensions of the associated invariant subspaces are determined via explicit character formulas. To finish, we conjecture that (for most graphs) the standard Lie algebras have only components that are either exponential or that grow, at most, polynomially with the system size. This would imply that the QAOA is either prone to barren plateaus or classically simulable. More generally, our work provides a symmetry framework and tools to analyze any desired variational quantum algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Heat and Momentum Transfer Studies in High Reynolds Number Wavy Films at Normal and Reduced Gravity Conditions

We examined the effect of the gas flow on the liquid film when the gas flows in the countercurrent direction in a vertical pipe at normal gravity conditions. The most dramatic effect of the simultaneous flow of gas and liquid in pipes is the greatly increased transport rates of heat, mass, and momentum. In practical situations this enhancement can be a benefit or it can result in serious operational problems. For example, gas-liquid flow always results in substantially higher pressure drop and this is usually undesirable. However, much higher heat transfer coefficients can be expected and this can obviously be of benefit for purposes of design. Unfortunately, designers know so little of the behavior of such two phase systems and as a result these advantages are not utilized. Due to the complexity of the second order boundary model as well as the fact that the pressure variation across the film is small compared to the imposed gas phase pressure, the countercurrent gas flow affect was studied for the standard boundary layer model. A different stream function that can compensate the shear stress affect was developed and this stream function also can predict periodic solutions. The discretized model equations were transformed to a traveling wave coordinate system. A stability analysis of these sets of equations showed the presence of a Hopf bifurcation for certain values of the traveling wave velocity and the shear stress. The Hopf celerity was increased due to the countercurrent shear. For low flow rate the increases of celerity are more than for the high flow rate, which was also observed in experiments. Numerical integration of a traveling wave simplification of the model also predicts the existence of chaotic large amplitude, nonperiodic waves as observed in the experiments. The film thickness was increased by the shear.

Balakotaiah, V.↗
NASA NTRS2005

Crop Surveillance Demonstration Using a Near-Daily MODIS Derived Vegetation Index Time Series

Effective response to crop disease outbreaks requires rapid identification and diagnosis of an event. A near-daily vegetation index product, such as a Normalized Difference Vegetation Index (NDVI), at moderate spatial resolution may serve as a good method for monitoring quick-acting diseases. NASA s Moderate Resolution Imaging Spectroradiometer (MODIS) instrument flown on the Terra and Aqua satellites has the temporal, spatial, and spectral properties to make it an excellent coarse-resolution data source for rapid, comprehensive surveillance of agricultural areas. A proof-of-concept wide area crop surveillance system using daily MODIS imagery was developed and tested on a set of San Joaquin cotton fields over a growing season. This area was chosen in part because excellent ground truth data were readily available. Preliminary results indicate that, at least in the southwestern part of the United States, near-daily NDVI products can be generated that show the natural variations in the crops as well as specific crop practices. Various filtering methods were evaluated and compared with standard MOD13 NDVI MODIS products. We observed that specific chemical applications that produce defoliation, which would have been missed using the standard 16-day product, were easily detectable with the filtered daily NDVI products.

McKellip, Rodney↗
DOE OSTI2020

Hierarchical Bayesian Thermonuclear Rate for the 7 Be(n,p) 7 Li Big Bang Nucleosynthesis Reaction

Big Bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than 1000 seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the universe. This work derives new 7 Be(n,p) 7 Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of 7 Be and 7 Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel, R -matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the 7 Be(n,p) 7 Li thermonuclear rates. From the resulting posteriors, we extract R -matrix parameters ($E_r, γ^2_n, γ^2_p$) and derive excitation energies and partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our 7 Be(n,p) 7 Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of ≤1 GK. Finally, we compare our rates to previous results and find that the 7 Be(n,p) 7 Li rates most commonly used in big bang simulations have uncertainties that are too optimistic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Raman intensity as a probe of concentration near a crystal growing in solution

The feasibility of using Raman spectral scattering signals for measurements of concentration profiles near a crystal interface during growth or dissolution is discussed. With KH2PO4 (KDP) as a test material, optical multichannel analyzer (OMA) detection of a solute Raman vibrational band provided direct quantification of solute concentration with band intensity. The intersection of incident laser and Raman collection optics provided 3-D selective point measurements of the solution concentration field. Unlike many other techniques, the Raman band intensity is not sensitive to the typical temperature variations. Precision calibration of Raman intensity versus KDP concentration with less than 1 pct standard deviation error levels was demonstrated. A fiber optic, which sampled incident laser intensity and coupled it to the OMA, provided a fully synchronized monitor of fluctuations in laser power to correlate with observed Raman signals. With 1 W of laser power at the sample, good data statistics required eight repeated data collections at approximately 2.5 min collection. The accumulated time represents the concentration measurement time at one spatial location. Photomicroscopy documented a 30 micrometer diameter by 200 micrometer of laser Raman scattering region in the solution near the crystal surface. The laser beam was able to approach up to 25 micrometer from the crystal surface. However, a crystal surface reflected intensity contribution was weakly detectable. Nucleated microcrystals were seen in the crystal-growing solution. These microcrystals convect right up to the crystal surface and indicate no quiet diffusion region under normal gravity conditions. Translation of the solution cell with respect to the optics caused systematic intensity errors.

Wilkinson, R. Allen↗