Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data Intensive Computing”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.1), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
NASA NTRS1984

Predicted changes in advanced turboprop noise with shaft angle of attack

Advanced turboprop blade designs and new propeller installation schemes motivated an effort to include unsteady loading effects in existing propeller noise prediction computer programs. The present work validates the prediction capability while studing the effects of shaft inclination on the radiated sound field. Classical methods of propeller performance analysis supply the time-dependent blade loading needed to calculate noise. Polar plots of the sound pressure level (SPL) of the first four harmonics and overall SPL are indicative of the change in directivity pattern as a function of propeller angle of attack. Noise predictions are compared with newly available wind tunnel data and the accuracy and applicability of the prediction method are discussed. It is concluded that unsteady blade loading caused by inclining the propeller with respect to the flow changes the directionality and the intensity of the radiated noise. These changes are well modeled by the present quasi-steady prediction method.

Padula, S. L.↗
NASA NTRS2000

X-Ray Astronomy

Dr. S. N. Zhang has lead a seven member group (Dr. Yuxin Feng, Mr. XuejunSun, Mr. Yongzhong Chen, Mr. Jun Lin, Mr. Yangsen Yao, and Ms. Xiaoling Zhang). This group has carried out the following activities: continued data analysis from space astrophysical missions CGRO, RXTE, ASCA and Chandra. Significant scientific results have been produced as results of their work. They discovered the three-layered accretion disk structure around black holes in X-ray binaries; their paper on this discovery is to appear in the prestigious Science magazine. They have also developed a new method for energy spectral analysis of black hole X-ray binaries; four papers on this topics were presented at the most recent Atlanta AAS meeting. They have also carried Monte-Carlo simulations of X-ray detectors, in support to the hardware development efforts at Marshall Space Flight Center (MSFC). These computation-intensive simulations have been carried out entirely on the computers at UAH. They have also carried out extensive simulations for astrophysical applications, taking advantage of the Monte-Carlo simulation codes developed previously at MSFC and further improved at UAH for detector simulations. One refereed paper and one contribution to conference proceedings have been resulted from this effort.

Wu, S. T.↗
NASA NTRS2011

Use of Apollo 17 Epoch Neutron Spectrum as a Benchmark in Testing LEND Collimated Sensor

The Apollo 17 neutron experiment LPNE provided a unique set of data on production of neutrons in the Lunar soil bombarded by Galactic Cosmic Rays (GCR). It serves as valuable "ground-truth" in the age of orbital remote sensing. We used the neutron data attributed to Apollo 17 epoch as a benchmark for testing the LEND's collimated sensor, as introduced by the geometry of collimator and efficiency of He3 counters. The latter is defined by the size of gas counter and pressure inside it. The intensity and energy spectrum of neutrons escaping the lunar surface are dependent on incident flux of Galactic Cosmic Rays (GCR) whose variability is associated with Solar Cycle and its peculiarities. We obtain first the share of neutrons entering through the field of view of collimator as a fraction of the total neutron flux by using the angular distribution of neutron exiting the Moon described by our Monte Carlo code. We computed next the count rate of the 3He sensor by using the neutron energy spectrum from McKinney et al. [JGR, 2006] and by consider geometry and gas pressure of the LEND sensor. Finally the neutron count rate obtained for the Apollo 17 epoch characterized by intermediate solar activity was adjusted to the LRO epoch characterized by low solar activity. It has been done by taking into account solar modulation potential, which affects the GCR flux, and in turn changes the neutron albedo flux.

Chin, Gordon↗
NASA NTRS2018

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗
NASA NTRS2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗
DOE OSTI2023

Power Density Distributions in Electron Beams

Measuring and controlling the power density distribution of electron beams used for welding is critical for producing repeatable welds and for transferring welding parameters between different machines. On any given machine, the power density distribution is controlled by defocusing the beam relative to its sharpest focused condition. However, measuring the power density distribution can be difficult due to the intense nature of welding beams and is further complicated by imperfect electron optics that can distort the beam, making it difficult to quantify. The enhanced modified Faraday cup (EMFC) diagnostic method was used here for beam analysis that utilizes computed tomography to reconstruct the beam’s power density distribution. These results were compared to the International Standards Organization (ISO) method for characterizing laser beams using a second-moment D4σ calculation. For ideal Gaussian-shaped beams, both methods would give the same result. However, for imperfect beams, the calculated D4σ diameter was shown to be about 25% larger relative to the FWe2 diameter measured by the EMFC due to the heavier weighting of data in the tails of the beam by D4σ. Finally, although both methods produce repeatable welds, it is important to understand the differences in the reported beam diameters, divergence angles, and beam parameter products when transferring parameters between machines.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Systems biology and biomarker discovery

Medical practitioners have always relied on surrogate markers of inaccessible biological processes to make their diagnosis, whether it was the pallor of shock, the flush of inflammation, or the jaundice of liver failure. Obviously, the current implementation of biomarkers for disease is far more sophisticated, relying on highly reproducible, quantitative measurements of molecules that are often mechanistically associated with the disease in question, as in glycated hemoglobin for the diagnosis of diabetes [1] or the presence of cardiac troponins in the blood for confirmation of myocardial infarcts [2]. In cancer, where the initial symptoms are often subtle and the consequences of delayed diagnosis often drastic for disease management, the impetus to discover readily accessible, reliable, and accurate biomarkers for early detection is compelling. Yet despite years of intense activity, the stable of clinically validated, cost-effective biomarkers for early detection of cancer is pathetically small and still dominated by a handful of markers (CA-125, CEA, PSA) first discovered decades ago. It is time, one could argue, for a fresh approach to the discovery and validation of disease biomarkers, one that takes full advantage of the revolution in genomic technologies and in the development of computational tools for the analysis of large complex datasets. This issue of Disease Markers is dedicated to one such new approach, loosely termed the 'Systems Biology of Biomarkers'. What sets the Systems Biology approach apart from other, more traditional approaches, is both the types of data used, and the tools used for data analysis - and both reflect the revolution in high throughput analytical methods and high throughput computing that has characterized the start of the twenty first century.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Synchronicity from synchronized chaos

The synchronization of loosely-coupled chaotic oscillators, a phenomenon investigated intensively for the last two decades, may realize the philosophical concept of “synchronicity”—the commonplace notion that related events mysteriously occur at the same time. When extended to continuous media and/or large discrete arrays, and when general (non-identical) correspondences are considered between states, intermittent synchronous relationships indeed become ubiquitous. Meaningful synchronicity follows naturally if meaningful events are identified with coherent structures, defined by internal synchronization between remote degrees of freedom; a condition that has been posited as necessary for synchronizability with an external system. The important case of synchronization between mind and matter is realized if mind is analogized to a computer model, synchronizing with a sporadically observed system, as in meteorological data assimilation. Evidence for the ubiquity of synchronization is reviewed along with recent proposals that: (1) synchronization of different models of the same objective process may be an expeditious route to improved computational modeling and may also describe the functioning of conscious brains; and (2) the nonlocality in quantum phenomena implied by Bell’s theorem may be explained in a variety of deterministic (hidden variable) interpretations if the quantum world resides on a generalized synchronization “manifold”.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗
NASA NTRS2004

Velocity and Temperature Measurement in Supersonic Free Jets Using Spectrally Resolved Rayleigh Scattering

The flow fields of unheated, supersonic free jets from convergent and convergent-divergent nozzles operating at M = 0.99, 1.4, and 1.6 were measured using spectrally resolved Rayleigh scattering technique. The axial component of velocity and temperature data as well as density data obtained from a previous experiment are presented in a systematic way with the goal of producing a database useful for validating computational fluid dynamics codes. The Rayleigh scattering process from air molecules provides a fundamental means of measuring flow properties in a non-intrusive, particle free manner. In the spectrally resolved application, laser light scattered by the air molecules is collected and analyzed using a Fabry-Perot interferometer (FPI). The difference between the incident laser frequency and the peak of the Rayleigh spectrum provides a measure of gas velocity. The temperature is measured from the spectral broadening caused by the random thermal motion and density is measured from the total light intensity. The present point measurement technique uses a CW laser, a scanning FPI and photon counting electronics. The 1 mm long probe volume is moved from point to point to survey the flow fields. Additional arrangements were made to remove particles from the main as well as the entrained flow and to isolate FPI from the high sound and vibration levels produced by the supersonic jets. In general, velocity is measured within +/- 10 m/s accuracy and temperature within +/- 10 K accuracy.

Panda, J.↗
NASA NTRS2012

A Distributive, Non-Destructive, Real-Time Approach to Snowpack Monitoring

This invention is designed to ascertain the snow water equivalence (SWE) of snowpacks with better spatial and temporal resolutions than present techniques. The approach is ground-based, as opposed to some techniques that are air-based. In addition, the approach is compact, non-destructive, and can be communicated with remotely, and thus can be deployed in areas not possible with current methods. Presently there are two principal ground-based techniques for obtaining SWE measurements. The first is manual snow core measurements of the snowpack. This approach is labor-intensive, destructive, and has poor temporal resolution. The second approach is to deploy a large (e.g., 3x3 m) snowpillow, which requires significant infrastructure, is potentially hazardous [uses a approximately equal to 200-gallon (approximately equal to 760-L) antifreeze-filled bladder], and requires deployment in a large, flat area. High deployment costs necessitate few installations, thus yielding poor spatial resolution of data. Both approaches have limited usefulness in complex and/or avalanche-prone terrains. This approach is compact, non-destructive to the snowpack, provides high temporal resolution data, and due to potential low cost, can be deployed with high spatial resolution. The invention consists of three primary components: a robust wireless network and computing platform designed for harsh climates, new SWE sensing strategies, and algorithms for smart sampling, data logging, and SWE computation.

Frolik, Jeff↗
DOE OSTI2026

Describing Point Defect Topology in 2D Energy Materials Through Computer Vision

Point defects such as vacancies and impurity atoms strongly impact the performance of 2D materials. Traditional efforts often rely on manual detection, a process that is time-intensive, prone to human error, and challenging to scale. Here we leverage machine learning (ML) methods to identify and quantify vacancies within 2D transition metal carbides (Ti3C2, MXenes), aiming to expedite detection while improving accuracy. MXenes exhibit valuable defect-defined electrochemical properties, but we currently lack statistical understanding of defect topology needed to fully harness these materials. Here we employ a convolutional neural network for semantic segmentation of experimental MXene images, opening an opportunity to conduct a rigorous statistical study on defect hierarchy while investigating local relaxation in the lattice. We show how the integration of ML can yield fundamental insight into point defects, providing a powerful tool that will play an increasingly crucial role in the future of materials science. ML is often not just a matter of straightforward application, and pretrained models proved ineffective in this case. Instead, we trained our own neural network (NN) and applied data augmentation techniques and fine-tuning to the training dataset. Since labeled microscopy data is often scarce, we developed training data from a previously published wide-frame MXene image, using customized Gaussian fitting to locate atomic positions. Our trained model was then applied to a large dataset of experimental images, enabling a statistical study of defect configurations across three samples prepared with different HF etchant concentrations (5%, 9.1%, and 12.5%), as shown in Fig. 1. This also allowed us to investigate local strain around vacancies, though we find that we are limited by the precision of measurements using high-angle annular dark field (HAADF) images, as shown in Fig. 2. This study demonstrates how ML enables large-scale, quantitative analysis of atomic defects - an otherwise infeasible task with traditional methods. While our NN was specialized for Ti3C2 MXenes, the pipeline we developed provides a foundation for future ML models tailored to other materials. Ultimately, we envision embedding the NN onto the microscope to give real-time feedback to the user. To make this a reality, continued work is necessary to fully understand the NN's capabilities and limitations. This study gets one step closer to our goals of automated experimentation moving away from traditional methods of manual labeling. As ML capabilities advance, we hope to continue adapting and applying these techniques in microscopy.

2D materials↗
DOE OSTI2025

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗
DOE OSTI2015

LOOKING FOR GRANULATION AND PERIODICITY IMPRINTS IN THE SUNSPOT TIME SERIES

The sunspot activity is the end result of the cyclic destruction and regeneration of magnetic fields by the dynamo action. We propose a new method to analyze the daily sunspot areas data recorded since 1874. By computing the power spectral density of daily data series using the Mexican hat wavelet, we found a power spectrum with a well-defined shape, characterized by three features. The first term is the 22 yr solar magnetic cycle, estimated in our work to be 18.43 yr. The second term is related to the daily volatility of sunspots. This term is most likely produced by the turbulent motions linked to the solar granulation. The last term corresponds to a periodic source associated with the solar magnetic activity, for which the maximum power spectral density occurs at 22.67 days. This value is part of the 22–27 day periodicity region that shows an above-average intensity in the power spectra. The origin of this 22.67 day periodic process is not clearly identified, and there is a possibility that it can be produced by convective flows inside the star. The study clearly shows a north–south asymmetry. The 18.43 yr periodical source is correlated between the two hemispheres, but the 22.67 day one is not correlated. It is shown that toward the large timescales an excess occurs in the northern hemisphere, especially near the previous two periodic sources. To further investigate the 22.67 day periodicity, we made a Lomb–Scargle spectral analysis. The study suggests that this periodicity is distinct from others found nearby.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Statistical and Spectral Analysis of Wind Characteristics Relevant to Wind Energy Assessment Using Tower Measurements in Complex Terrain

The main objective of the study was to investigate spatial and temporal characteristics of the wind speed and direction in complex terrain that are relevant to wind energy assessment and development, as well as to wind energy system operation, management, and grid integration. Wind data from five tall meteorological towers located in Western Nevada, USA, operated from August 2003 to March 2008, used in the analysis. The multiannual average wind speeds did not show significant increased trend with increasing elevation, while the turbulence intensity slowly decreased with an increase were the average wind speed. The wind speed and direction were modeled using the Weibull and the von Mises distribution functions. The correlations show a strong coherence between the wind speed and direction with slowly decreasing amplitude of the multiday periodicity with increasing lag periods. The spectral analysis shows significant annual periodicity with similar characteristics at all locations. The relatively high correlations between the towers and small range of the computed turbulence intensity indicate that wind variability is dominated by the regional synoptic processes. Knowledge and information about daily, seasonal, and annual wind periodicities are very important for wind energy resource assessment, wind power plant operation, management, and grid integration.

Belu, Radian↗