Search NASA⌕ Search

SEARCH · Search NASA

Results for “Estimation method”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quasar Identification Using Multivariate Probability Density Estimated from Nonparametric Conditional Probabilities

Nonparametric estimation for a probability density function that describes multivariate data has typically been addressed by kernel density estimation (KDE). A novel density estimator recently developed by Farmer and Jacobs offers an alternative high-throughput automated approach to univariate nonparametric density estimation based on maximum entropy and order statistics, improving accuracy over univariate KDE. This article presents an extension of the single variable case to multiple variables. The univariate estimator is used to recursively calculate a product array of one-dimensional conditional probabilities. In combination with interpolation methods, a complete joint probability density estimate is generated for multiple variables. Good accuracy and speed performance in synthetic data are demonstrated by a numerical study using known distributions over a range of sample sizes from 100 to 10 6 for two to six variables. Performance in terms of speed and accuracy is compared to KDE. The multivariate density estimate developed here tends to perform better as the number of samples and/or variables increases. As an example application, measurements are analyzed over five filters of photometric data from the Sloan Digital Sky Survey Data Release 17. The multivariate estimation is used to form the basis for a binary classifier that distinguishes quasars from galaxies and stars with up to 94% accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Efficient calibration for imperfect computer models

Many computer models contain unknown parameters which need to be estimated using physical observations. Furthermore, the calibration method based on Gaussian process models may lead to unreasonable estimate for imperfect computer models. In this work, we extend their study to calibration problems with stochastic physical data. We propose a novel method, called the L 2 calibration, and show its semiparametric efficiency. The conventional method of the ordinary least squares is also studied. Theoretical analysis shows that it is consistent but not efficient. Here, numerical examples show that the proposed method outperforms the existing ones.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Spillover as a cause of bias in baseline evaluation methods for demand response programs

Prior research has shown load reduction estimates from residential event-driven demand response programs (e.g., Critical Peak Pricing)using X of the highest Y days with a weather adjustment method are the best performing within the class of currently used baseline methods. However, they are still biased relative to estimates produced from randomized control trials (RCTs), the unbiased “gold standard” evaluation method. In this paper we identify underlying factors that cause some of the bias found in one commonly used baseline method. Rather than simply quantifying bias, our research provides a deeper understanding of what causes the bias and thus can be used to develop more accurate methods that are not subject to these underlying factors. Previous studies have compared various baseline methods relative to each other; however, because all baseline methods are biased, it is impossible to determine the true bias that exists in them. We have access to a unique RCT dataset: the Sacramento Municipal Utility District's study of critical peak pricing. Our analysis of 23 event days over two summers allows us to identify the true bias on load reduction estimates by using the RCT estimates as the unbiased gold standard against which we compare the estimates from the baseline methods. We found that spillover of energy reductions, from hours targeted by a program onto other hours, is one underlying factor that is a major cause of bias in baseline methods. We discuss alternative baseline methods that may not be subject to this same bias.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

On a framework for generating PoD curves assisted by numerical simulations

The Probability of Detection (PoD) curve method has emerged as an important tool for the assessment of the performance of NDE techniques, a topic of particular interest to the nuclear industry where inspection qualification is very important. The conventional experimental means of generating PoD curves though, can be expensive, requiring large data sets (covering defects and test conditions), and equipment and operator time. Several methods of achieving faster estimates for PoD curves using physics-based modelling have been developed to address this problem. Numerical modelling techniques are also attractive, especially given the ever-increasing computational power available to scientists today. Here we develop procedures for obtaining PoD curves, assisted by numerical simulation and based on Bayesian statistics. Numerical simulations are performed using Finite Element analysis for factors that are assumed to be independent, random and normally distributed. PoD curves so generated are compared with experiments on austenitic stainless steel (SS) plates with artificially created notches. We examine issues affecting the PoD curve generation process including codes, standards, distribution of defect parameters and the choice of the noise threshold. We also study the assumption of normal distribution for signal response parameters and consider strategies for dealing with data that may be more complex or sparse to justify this. These topics are addressed and illustrated through the example case of generation of PoD curves for pulse-echo ultrasonic inspection of vertical surface-breaking cracks in SS plates.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Global carbon budget 2014

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO 2 emissions from fossil fuel combustion and cement production ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The mean ocean CO 2 sink ( S OCEAN ) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S OCEAN is evaluated with data products based on surveys of ocean CO 2 measurements. The global residual terrestrial CO 2 sink ( S LAND ) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO 2 , and land-cover-change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ;, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2004–2013), E FF was 8.9 ± 0.4 GtC yr⁻¹, E LUC 0.9 ± 0.5 GtC yr⁻¹, G ATM 4.3 ± 0.1 GtC yr⁻¹, S OCEAN 2.6 ± 0.5 GtC yr⁻¹, and S LAND 2.9 ± 0.8 GtC yr⁻¹. For year 2013 alone, E FF grew to 9.9 ± 0.5 GtC yr⁻¹, 2.3% above 2012, continuing the growth trend in these emissions, E LUC was 0.9 ± 0.5 GtC yr⁻¹, G ATM was 5.4 ± 0.2 GtC yr⁻¹, S OCEAN was 2.9 ± 0.5 GtC yr⁻¹, and S LAND was 2.5 ± 0.9 GtC yr⁻¹. G ATM was high in 2013, reflecting a steady increase in E FF and smaller and opposite changes between S OCEAN and S LAND compared to the past decade (2004–2013). The global atmospheric CO 2 concentration reached 395.31 ± 0.10 ppm averaged over 2013. We estimate that E FF will increase by 2.5% (1.3–3.5%) to 10.1 ± 0.6 GtC in 2014 (37.0 ± 2.2 GtCO 2 yr⁻¹), 65% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the global economy. From this projection of E FF and assumed constant E LUC for 2014, cumulative emissions of CO 2 will reach about 545 ± 55 GtC (2000 ± 200 GtCO 2 ) for 1870–2014, about 75% from E FF and 25% from E LUC . This paper documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this living data set (Le Quéré et al., 2013, 2014). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center ( doi:10.3334/CDIAC/GCP_2014 ).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Managing negative values is reservoir inflow computation: A case study

Reservoir inflow is conventionally estimated using the water balance method, which involves the reservoir release and the change in storage during the period considered. As a result, the estimated inflow may sometimes be negative as the errors involved in each input variable build-up to the output. In our study, the fleet data was provided by the Tennessee Valley Authority (TVA) for their Norris Hydropower facility. Unlike the flow release data, which was readily accessible, the change in storage had to be calculated using the reservoir elevation and volume relationship. The original inflow estimates produced a wide range of negative values with large outliers, making it difficult to visualize the current trends. This paper describes a methodology to remove the negative values encountered during the inflow computation, and the results were analyzed by correlating with the nearby streamflow gaging stations.

Shibu, Asha↗
DOE OSTI2020

Development of a Comparison Framework for Evaluating Environmental Contours of Extreme Sea States

Environmental contours of extreme sea states are often utilized for the purposes of reliability-based offshore design. Many methods have been proposed to estimate environmental contours of extreme sea states, including, but not limited to, the traditional inverse first-order reliability method (I-FORM) and subsequent modifications, copula methods, and Monte Carlo methods. These methods differ in terms of both the methodology selected for defining the joint distribution of sea state parameters and in the method used to construct the environmental contour from the joint distribution. It is often difficult to compare the results of proposed methods to determine which method should be used for a particular application or geographical region. The comparison of the predictions from various contour methods at a single site and across many sites is important to making environmental contours of extreme sea states useful in practice. The goal of this paper is to develop a comparison framework for evaluating methods for developing environmental contours of extreme sea states. This paper develops generalized metrics for comparing the performance of contour methods to one another across a collection of study sites, and applies these metrics and methods to develop conclusions about trends in the wave resource across geographic locations, as demonstrated for a pilot dataset. These proposed metrics and methods are intended to judge the environmental contours themselves relative to other contour methods, and are thus agnostic to a specific device, structure, or field of application. The metrics developed and applied in this paper include measures of predictive accuracy, physical validity, and aggregated temporal performance that can be used to both assess contour methods and provide recommendations for the use of certain methods in various geographical regions. The application and aggregation of the metrics proposed in this paper outline a comparison framework for environmental contour methods that can be applied to support design analysis workflows for offshore structures. This comparison framework could be extended in future work to include additional metrics of interest, potentially including those to address issues pertinent to a specific application area or analysis discipline, such as metrics related to structural response across contour methods or additional physics-based metrics based on wave dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

TH-AB-202-12: The First Clinical Implementation of a Real-Time Six Degree of Freedom Tracking System During Radiation Therapy

Purpose: In current practice, imaging is typically performed prior to treatment; the cancer target motion during treatment is unknown. We present the first clinical implementation of real time Kilovoltage Intrafraction Monitoring (KIM) system which tracks the cancer target translational and rotational motions during treatment. Methods: KIM technology: KIM estimates the 3D position of the target tumour based on segmented 2D positions of the three implanted fiducials in each of the kV images (125 kV, 10 mA at 11 fps) taken continuously during the treatment arc. The 2D-3D target estimation is based on a probability distribution function, obtained during pre-treatment CBCT. Rotations about each axis with the centroid of the markers as the pivot were calculated using the iterative closest point algorithm in real time. Patient: A patient with prostate adenocarcinoma undergoing stereotactic body radiotherapy (SBRT) with 36.25 Gy delivered in 5 fractions (Varian TrueBeam, 6X, VMAT) was enrolled in the study. The trial complies with Australian ethical and regulatory standards. Results: Of the 5 fractions of treatment the patient received, KIM was utilised successfully in 4 fractions with 3 couch shifts due to large persistent prostate movements (>2mm for more than 5 seconds). KIM translational accuracy and precision in comparison with post treatment kV-MV triangulation are 0.28±0.59 mm, −0.19±0.25 mm and 0.23±0.69 mm in the Left-Right, Superior-Inferior and Anterior-Posterior directions, respectively. KIM rotational accuracy as compared with triangulation is: 0.429°±2.22°, −0.44°±4.7° and 0.06°±1.08° in the roll, pitch and yaw direction, respectively. Conclusion: The first six degree of freedom KIM system was successfully implemented clinically. The presented KIM system has sub-millimeters accuracy and precision in all three translational axes, and less than 1° of mean error in all three rotational axes. Acknowledgement: This work is supported by Cancer Australia grants APPXXX, APPYYY.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Bayesian Model Averaging for Ensemble-Based Estimates of Solvation Free Energies

This paper applies the Bayesian Model Averaging (BMA) statistical ensemble technique to estimate small molecule solvation free energies. There is a wide range methods for predicting solvation free energies, ranging from empirical statistical models to {\it ab initio} quantum mechanical approaches. Each of these methods are based on a set of conceptual assumptions that can affect a method's predictive accuracy and transferability. Using an iterative statistical process, we have selected and combined solvation energy estimates using an ensemble of 17 diverse methods from the SAMPL4 blind prediction study to form a single, aggregated solvation energy estimate. The ensemble design process evaluates the statistical information in each individual method as well as the performance of the aggregate estimate obtained from the ensemble as a whole. Methods that possess minimal or redundant information are pruned from the ensemble and the evaluation process repeats until aggregate predictive performance can no longer be improved. We show that this process results in a final aggregate estimate that outperforms all individual methods by reducing estimate errors by as much as 91% to 1.2 kcal/mol accuracy. We also compare our iterative refinement approach to other statistical ensemble approaches and demonstrate that this iterative process reduces estimate errors by as much as 61%. This work provides a new approach for accurate solvation free energy prediction and lays the foundation for future work on aggregate models that can balance computational cost with predictive accuracy.

solvation, Bayesian model averaging, Electrostatic↗
DOE OSTI2025

Hierarchical Speed Planner for Automated Vehicles: A Framework for Lagrangian Variable Speed Limit in Mixed-Autonomy Traffic

Here, this article presents a novel hierarchical speed planning framework for variable speed limits in mixed-autonomy traffic environments, leveraging server-side macroscopic control and vehicle-side microscopic execution. The framework integrates real-time traffic state estimation (TSE) and reinforcement learning (RL)-based control to mitigate congestion and improve traffic flow. A TSE enhancement module combines macroscopic data from sources like INRIX with high-resolution observations from connected autonomous vehicles (CAVs), enabling predictive modeling to address latency and noise. The target speed design module employs kernel smoothing and a buffer zone strategy to optimize traffic density and flow around bottlenecks. The proposed system was validated in the largest open-road test to date with 100 CAVs, demonstrating an overall 8% traffic density decrease, with a specific decrease of 7% upstream, 10% downstream, and a 52% decrease during the congestion formation phase at bottlenecks.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗
DOE OSTI2022

Transportable absolute yields of underground nuclear explosions. Application to the North Korean nuclear tests

The method developed and presented here improves our nuclear explosion monitoring capability by providing a new independent and transportable method to determine the explosive yields of seismically recorded events associated with foreign weapons development. Here we report new independent absolute yield and depth of burial estimates of the North Korean nuclear tests which are commensurate with previously determined source parameters. The key advantage of this method over other methods which estimate absolute explosive yields is that this method does not require any a priori calibration, and can be immediately applied to any region of interest. This method uses the source information obtained from the spectral ratios of envelopes of measured seismic coda waves to invert for model-derived yield and depth of burial. The only requirements of this method are that the region of interest must contain several (> 3) seismic sources which have been well recorded (SNR > 2) by a set of shared stations. For regions with only a few events (~ 3) and/or extremely band-limited observations for which the depth of burial is unknown, an earthquake with a known magnitude may be required to obtain an accurate and robust yield estimate.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗
DOE OSTI2015

Global Carbon Budget 2015

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We also discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, while emissions from land-use change ( E LUC ), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. Moreover, the mean ocean CO 2 sink ( S OCEAN ) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S OCEAN is evaluated with data products based on surveys of ocean CO 2 measurements. The global residual terrestrial CO 2 sink ( S LAND ) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO 2 , and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), E FF was 9.0 ± 0.5 GtC yr -1 , E LUC was 0.9 ± 0.5 GtC yr -1 , G ATM was 4.4 ± 0.1 GtC yr -1 , S OCEAN was 2.6 ± 0.5 GtC yr -1 , and S LAND was 3.0 ± 0.8 GtC yr -1 . For the year 2014 alone, E FF grew to 9.8 ± 0.5 GtC yr -1 , 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr -1 that took place during 2005–2014. Also, for 2014, E LUC was 1.1 ± 0.5 GtC yr -1 , G ATM was 3.9 ± 0.2 GtC yr -1 , S OCEAN was 2.9 ± 0.5 GtC yr -1 , and S LAND was 4.1 ± 0.9 GtC yr -1 . G ATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger S LAND for that year. The global atmospheric CO 2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in E FF will be near or slightly below zero, with a projection of -0.6 [range of -1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of E FF and assumed constant E LUC for 2015, cumulative emissions of CO 2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO 2 ) for 1870–2015, about 75 % from E FF and 25 % from E LUC . Finally, this living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center ( doi:10.3334/CDIAC/GCP_2015 ).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Uncertainty-aware Continuous Implicit Neural Representations for Remote Sensing Object Counting

Many existing object counting methods rely on density map estimation (DME) of the discrete grid representation by decoding extracted image semantic features from designed convolutional neural networks (CNNs). Relying on discrete density maps not only leads to information loss dependent on the original image resolution, but also has a scalability issue when analyzing high-resolution images with cubically increasing memory complexity. Furthermore, none of the existing methods can offer reliable uncertainty quantification (UQ) for the derived count estimates. To overcome these limitations, we design UNcertainty-aware, hypernetwork-based Implicit neural representations for Counting (UNIC) to assign probabilities and the corresponding counting confidence over continuous spatial coordinates. We derive a sampling-based Bayesian counting loss function and develop the corresponding model training algorithm. UNIC outperforms existing methods on the Remote Sensing Object Counting (RSOC) dataset with reliable UQ and improved interpretability of the derived count estimates. Our code is available at https://github.com/SiyuanXu-tamu/UNIC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spherical Indentation and Implementation of S 3 /P for yield stress determination of brittle materials

Here, a mathematically transparent and robust experimental method has been developed to estimate the yield stress of brittle materials through the analysis of depth-sensing spherical indentation. Employing Hertzian contact mechanics, an elastically invariant ratio based on the simple equation S 3 / P = 6RE r 2 , (where S and P are contact stiffness and indentation load, respectively) has been derived that enables more accurate and confident determination of the transition from elastic to inelastic deformation; a transition that the yield stress dictates and represents. Using two diamond spheres with radii of 3.2 and 8.6 μm, the indentation test method and analyses are applied to two vitreous silicates: Corning's HPFS 7980® fused silica and Vitro's Starphire® soda lime silicate. The estimated yield strengths are 8.15 GPa ± 2.5% for the fused silica and 6.1 GPa ± 3.3% for the soda-lime silicate, and both were independent of indenter radius. Verification of this new experimental method is demonstrated with an as-drawn titanium by showing equivalence of measured yield stress by its spherical indentation and that from uniaxial compression testing. This method will enable easier and more confident estimation of yield stress in brittle materials - a property that historically has been elusive to measure for these materials using common laboratory mechanical test methods.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Energetic electron precipitation associated with pulsating aurora: EISCAT and Van Allen Probe observations

Pulsating auroras show quasi-periodic intensity modulations caused by the precipitation of energetic electrons of the order of tens of keV. It is expected theoretically that not only these electrons but also subrelativistic/relativistic electrons precipitate simultaneously into the ionosphere owing to whistler mode wave-particle interactions. The height-resolved electron density profile was observed with the European Incoherent Scatter (EISCAT) Tromsø VHF radar on 17 November 2012. Electron density enhancements were clearly identified at altitudes >68 km in association with the pulsating aurora, suggesting precipitation of electrons with a broadband energy range from ~10 keV up to at least 200 keV. The riometer and network of subionospheric radio wave observations also showed the energetic electron precipitations during this period. During this period, the footprint of the Van Allen Probe-A satellite was very close to Tromsø and the satellite observed rising tone emissions of the lower band chorus (LBC) waves near the equatorial plane. Considering the observed LBC waves and electrons, we conducted a computer simulation of the wave-particle interactions. This showed simultaneous precipitation of electrons at both tens of keV and a few hundred keV, which is consistent with the energy spectrum estimated by the inversion method using the EISCAT observations. This result revealed that electrons with a wide energy range simultaneously precipitate into the ionosphere in association with the pulsating aurora, providing the evidence that pulsating auroras are caused by whistler chorus waves. We suggest that scattering by propagating whistler simultaneously causes both the precipitations of subrelativistic electrons and the pulsating aurora.

58 GEOSCIENCES↗
DOE OSTI2017

Turbulent mixing and removal of ozone within an Amazon rainforest canopy

Simultaneous profiles of turbulence statistics and mean ozone mixing ratio are used to establish a relation between eddy diffusivity and ozone mixing within the Amazon forest. A one–dimensional diffusion model is proposed and used to infer mixing time scales from the eddy diffusivity profiles. Data and model results indicate that during daytime conditions, the upper (lower) half of the canopy is well (partially) mixed most of the time and that most of the vertical extent of the forest can be mixed in less than an hour. During nighttime, most of the canopy is predominantly poorly mixed, except for periods with bursts of intermittent turbulence. Even though turbulence is faster than chemistry during daytime, both processes have comparable time scales in the lower canopy layers during nighttime conditions. Nonchemical loss time scales (associated with stomatal uptake and dry deposition) for the entire forest are comparable to turbulent mixing time scale in the lower canopy during the day and in the entire canopy during the night, indicating a tight coupling between turbulent transport and dry deposition and stomatal uptake processes. Because of the significant time of day and height variability of the turbulent mixing time scale inside the canopy, it is important to take it into account when studying chemical and biophysical processes happening in the forest environment. In conclusion, the method proposed here to estimate turbulent mixing time scales is a reliable alternative to currently used models, especially for situations in which the vertical distribution of the time scale is relevant.

54 ENVIRONMENTAL SCIENCES↗