Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “testing statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

Testing the thermal Sunyaev-Zel’dovich power spectrum of a halo model using hydrodynamical simulations

Statistical properties of large-scale cosmological structures serve as powerful tools for constraining the cosmological properties of our Universe. Tracing the gas pressure, the thermal Sunyaev-Zel’dovich (tSZ) effect is a biased probe of mass distribution and, hence, can be used to test the physics of feedback or cosmological models. Therefore, it is crucial to develop robust modelling of hot gas pressure for applications to tSZ surveys. Since gas collapses into bound structures, it is expected that most of the tSZ signal is within halos produced by cosmic accretion shocks. Hence, simple empirical halo models can be used to predict the tSZ power spectra. In this study, we employed the HMx halo model to compare the tSZ power spectra with those of several hydrodynamical simulations: the Horizon suite and the Magneticum simulation. We examine various contributions to the tSZ power spectrum across different redshifts, including the one- and two-halo term decomposition, the amount of bound gas, the importance of different masses, and the electron pressure profiles. Our comparison of the tSZ power spectrum reveals discrepancies between the halo model and cosmological simulations that increase with redshift. We find a 20% to 50% difference between the measured and predicted tSZ angular power spectrum over the multipole range ℓ = 10 3 − 10 4 . Our analysis reveals that these differences are driven by the excess of power in the predicted two-halo term at low k and in the one-halo term at high k . At higher redshifts ( z ∼ 3), simulations indicate that more power comes from outside the virial radius than from inside, suggesting a limitation in the applicability of the halo model. We also observe differences in the pressure profiles, despite the fair level of agreement on the tSZ power spectrum at low redshift with the default calibration of the halo model. In conclusion, our study suggests that the properties of the halo model need to be carefully controlled against real or mock data to be proven useful for cosmological purposes.

Ayçoberry, Emma (ORCID:0000000292351195)↗

Gender Consideration in Experiment Design for Airbrake in Prebreathe

If gender is a confounder of the decompression sickness (DCS) or venous gas emboli (VGE) outcomes of a proposed air break in oxygen prebreathe (PB) project, then decisions about the final experiment design must be made. We evaluated if the incidence of DCS and VGE from tests in altitude chambers over 20 years were different between men and women after resting and exercise prebreathe protocols. Nitrogen washout during PB is our primary risk mitigation strategy to prevent subsequent DCS and VGE in subjects. Bubbles in the pulmonary artery (venous blood) were detected from the precordial position using Doppler ultrasound bubble detectors. The subjects were monitored for VGE for four min at about 15 min intervals for the duration of the altitude exposure, with maximum bubble grade assigned a Spencer Grade of IV. There was no difference in DCS incidence between men and women in either PB protocol. The incidence of VGE and Grade IV VGE is statistically lower in women compared to men after resting PB. Even when 10 tests were compared with Mantel-Haenszel 2 where both men (n = 168) and women (n = 92) appeared, the p-value for VGE incidence was still significant at 0.03. The incidence of VGE and Grade IV VGE is not statistically lower in women compared to men after exercise PB. Even when six tests were compared with Mantel-Haenszel x2 where both men (n = 165) and women (n = 49) appeared, the p-value for VGE incidence was still not significant at 0.90. Our goal is to understand the risk of brief air breaks during PB without other confounding variables invalidating our conclusions. The cost to additionally account for the confounding role of gender on VGE outcome after resting PB is judged excessive. Our decision is to only evaluate air breaks in the exercise PB protocol. So there is no restriction to recruiting women as test subjects.

Conkin, Johnny↗

Effect of Aluminum Titanate (Al 2 TiO 5 ) Doping on the Mechanical Performance of Solid Oxide Fuel Cell Ni-YSZ Anode

The mechanical behavior of un-doped and 1–10 wt.% aluminum titanate doped NiO-YSZ anodes was evaluated in the oxidized and reduced state of the material. Sample bars 25 × 5 × 2 mm were fabricated and tested in a three-point bending apparatus and statistical analyzes of the collected data were performed. A remarkable enhancement of the flexural strength was found both for the reduced and oxidized state of the material when compared to the un-doped samples. In both cases, the development of a secondary phase was observed proportional to the doping amount of aluminum titanate. Morphological analysis along with preliminary phase identification and suspected mechanisms for strength enhancement are presented and discussed.

25 ENERGY STORAGE↗

A formation water-based nutrient recipe for potentially increasing methane release from coal in situ

Biogasifying coal to methane represents an environmentally benign way to utilize the abundant and inexpensive coal resource. To increase methane yield from coal, numerous studies have investigated the approach of biostimulation through finding the best nutrient solutions to enhance microbial activities. Toward this end, however, almost all studies have adopted laboratory made medium that is tap water- or deionized and distilled water- based. As a matter of fact, this water is dramatically different from formation water in coal basins. Thus, in order to enhance methane release from coal in situ, this study aimed to design a formation water-based recipe. To accomplish this objective, the chemical and microbial compositions of the formation water collected from the San Juan basin were analyzed first. Equipped with this fundamental knowledge, a screening test was conducted to evaluate nine parameters to identify statistically significant ones affecting methane yield from coal. For those critical parameters, the optimal value for each was determined through response surface methodology. Finally, the predicted results by the models were verified by an experimental study adopting all optimum conditions. Here this study demonstrated that microbes capable of converting coal to methane were present at the San Juan basin and the developed recipe increased methane yield 24.3-fold compared to those without.

01 COAL, LIGNITE, AND PEAT↗

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗

Deep neural network improves the estimation of polygenic risk scores for breast cancer

Polygenic risk scores (PRS) estimate the genetic risk of an individual for a complex disease based on many genetic variants across the whole genome. Here, we compared a series of computational models for estimation of breast cancer PRS. A deep neural network (DNN) was found to outperform alternative machine learning techniques and established statistical algorithms, including BLUP, BayesA, and LDpred. In the test cohort with 50% prevalence, the Area Under the receiver operating characteristic Curve (AUC) were 67.4% for DNN, 64.2% for BLUP, 64.5% for BayesA, and 62.4% for LDpred. BLUP, BayesA, and LPpred all generated PRS that followed a normal distribution in the case population. However, the PRS generated by DNN in the case population followed a bimodal distribution composed of two normal distributions with distinctly different means. This suggests that DNN was able to separate the case population into a high-genetic-risk case subpopulation with an average PRS significantly higher than the control population and a normal-genetic-risk case subpopulation with an average PRS similar to the control population. This allowed DNN to achieve 18.8% recall at 90% precision in the test cohort with 50% prevalence, which can be extrapolated to 65.4% recall at 20% precision in a general population with 12% prevalence. Interpretation of the DNN model identified salient variants that were assigned insignificant p values by association studies, but were important for DNN prediction. These variants may be associated with the phenotype through nonlinear relationships.

59 BASIC BIOLOGICAL SCIENCES↗

Bayesian stability and force modeling for uncertain machining processes

Accurately simulating machining operations requires knowledge of the cutting force model and system frequency response. However, this data is collected using specialized instruments in an ex-situ manner. Bayesian statistical methods instead learn the system parameters using cutting test data, but to date, these approaches have only considered milling stability. This paper presents a physics-based Bayesian framework which incorporates both spindle power and milling stability. Initial probabilistic descriptions of the system parameters are propagated through a set of physics functions to form probabilistic predictions about the milling process. The system parameters are then updated using automatically selected cutting tests to reduce parameter uncertainty and identify more productive cutting conditions, where spindle power measurements are used to learn the cutting force model. The framework is demonstrated through both numerical and experimental case studies. Results show that the approach accurately identifies both the system natural frequency and cutting force model.

42 ENGINEERING↗

Generalized framework for likelihood-based field-level inference of growth rate from velocity and density fields

Measuring the growth rate of large-scale structures ( f ) as a function of redshift has the potential to break degeneracies between modified gravity and dark energy models, when combined with expansion-rate probes. Direct estimates of peculiar velocities of galaxies have attracted interest as a means of estimating fσ 8 . In particular, field-level methods can be used to fit the field nuisance parameter along with cosmological parameters simultaneously. This article aims to provide the community with a unified framework for the theoretical modeling of the likelihood-based field-level inference by performing fast field covariance calculations for velocity and density fields. Our purpose is to lay the foundations for a nonlinear extension of the likelihood-based method at the field level. We have developed a generalized framework, implemented in the dedicated software flip to perform a likelihood-based inference of fσ 8 . We derived a new field covariance model, which includes wide-angle corrections. We also included the models previously described in the literature inside our framework. We compared their performance against ours, and we validated our model by comparing it with the two-point statistics of a recent N-body simulation. The tests we performed have allowed us to validate our software and determine the appropriate wavenumber range to integrate our covariance model and its validity in terms of separation. Our framework allows for a wider wavenumber coverage to be used in our calculations than in previous works, which is particularly interesting for nonlinear model extensions. Finally, our generalized framework allows us to efficiently perform a survey geometry-dependent Fisher forecast of the fσ 8 parameter. We show that the Fisher forecast method we developed gives an error bar that is 30% closer to a full likelihood-based estimation than a standard volume Fisher forecast.

Ravoux, Corentin↗

A hybrid approach to excited-state-specific variational Monte Carlo and doubly excited states

We extend our hybrid linear-method/accelerated-descent variational Monte Carlo optimization approach to excited states and investigate its efficacy in double excitations. In addition to showing a superior statistical efficiency when compared to the linear method, our tests on small molecules show good energetic agreement with benchmark methods. We also demonstrate the ability to treat double excitations in systems that are too large for a full treatment by using selected configuration interaction methods via an application to 4-aminobenzonitrile. Finally, we investigate the stability of state-specific variance optimization against collapse to other states’ variance minima and find that symmetry, Ansatz quality, and sample size all have roles to play in achieving stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

QUIJOTE scientific results – IV. A northern sky survey in intensity and polarization at 10–20 GHz with the multifrequency instrument

ABSTRACT We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17, and 19 GHz, and covering approximately 29 000 deg2, including most of the northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE multifrequency instrument (MFI), and have angular resolutions of around 1°, and sensitivities in polarization within the range 35–40 µK per 1° beam, being a factor ∼2–4 worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterize the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5 per cent. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A, and 3C274), which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11 GHz and WMAP 23 GHz map peaks at β = −3.09 with a standard deviation of 0.14. The measured BB/EE ratio at scales of ℓ = 80 is 0.26 ± 0.07 for a Galactic cut |b| > 10°. We find a positive TE correlation for 11 GHz at large angular scales (ℓ ≲ 50), while the EB and TB signals are consistent with zero in the multipole range 30 ≲ ℓ ≲ 150. The maps discussed in this paper are publicly available.

Astronomy & Astrophysics↗

A unified catalogue-level reanalysis of stage-III cosmic shear surveys

ABSTRACT Cosmological parameter constraints from recent galaxy imaging surveys are reaching percent-level accuracy on the effective amplitude of the lensing signal, S8. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce subpercent level measurements of cosmological parameters, providing a milestone test of the ΛCDM model. To supply guidance to the upcoming LSST analysis, it is important to understand thoroughly the results from different recent galaxy imaging surveys and assess their consistencies. In this work, we perform a unified catalogue-level reanalysis of three cosmic shear data sets: the first year data from the Dark Energy Survey (DES-Y1), the 1000 deg2 data set from the Kilo-Degree Survey (KiDS-1000), and the first year data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-Y1). We utilize a pipeline developed and rigorously tested by the LSST Dark Energy Science Collaboration to perform the reanalysis and assess the robustness of the results to analysis choices. We find the S8 constraint to be robust to two different small-scale modelling approaches, and varying choices of cosmological priors. Our unified analysis allows the consistency of the surveys to be rigorously tested, and we find the three surveys to be statistically consistent. Due to the partially overlapping footprint, we model the cross-covariance between KiDS-1000 and HSC-Y1 approximately when combining all three data sets, resulting in a 1.6–1.9 per cent constraint on S8 given different assumptions on the cross-covariance.

79 ASTRONOMY AND ASTROPHYSICS↗

Impacts of floating-point non-associativity on reproducibility for HPC and deep learning applications

Run to run variability in parallel programs caused by floating-point non-associativity has been known to significantly affect reproducibility in iterative algorithms, due to accumulating errors. Non-reproducibility can critically affect the efficiency and effectiveness of correctness testing for stochastic programs. Recently, the sensitivity of deep learning training and inference pipelines to floating-point non-associativity has been found to sometimes be extreme. It can prevent certification for commercial applications, accurate assessment of robustness and sensitivity, and bug detection. New approaches in scientific computing applications have coupled deep learning models with high-performance computing, leading to an aggravation of debugging and testing challenges. Here we perform an investigation of the statistical properties of floating-point non-associativity within modern parallel programming models, and analyze performance and productivity impacts of replacing atomic operations with deterministic alternatives on GPUs. We examine the recently-added deterministic options in PyTorch within the context of GPU deployment for deep learning, uncovering and quantifying the impacts of input parameters triggering run to run variability and reporting on the reliability and completeness of the documentation. Finally, we evaluate the strategy of exploiting automatic determinism that could be provided by deterministic hardware, using the Groq LPUTM accelerator for inference portions of the deep learning pipeline. We demonstrate the benefits that a hardware-based strategy can provide within reproducibility and correctness efforts.

Shanmugavelu, Sanjif↗

High Pressure Hydrogen Pressure Relief Devices: Accelerated Life Testing and Application Best Practices

Pressure relief devices (PRDs ) are used to protect high pressure systems from burst failure caused by overpressurization. Codes and standards require the use of PRDs for the safe design of many pressurized systems. These systems require high reliability due to the risks associated with a burst failure. Hydrogen service can increase the risk of PRD failure due to material property degradation caused by hydrogen attack. The National Renewable Energy Laboratory (NREL) has conducted an accelerated life test on a conventional spring loaded PRD. Based on previous failures in the field, the nozzles specific to these PRDs are of particular interest. A nozzle in a PRD is a small part that directs the flow of fluid toward the sealing surface to maintain the open state of the valve once the spring force is overcome. The nozzle in this specific PRD is subjected to the full tensile force of the fluid pressure. These nozzles are made from 440C material, which is a type of hardened steel that is commonly chosen for high pressure applications because of its high strength properties. In a hydrogen environment, however, 440C is considered a worst case material since hydrogen attack results in a loss of almost all ductility and thus 440C is prone to fatigue and material failure. Accordingly, 440C is not recommended for hydrogen service. Conducting an accelerated life test on a PRD with 440C material provides information on necessary and sufficient conditions required to produce crack initiation and failure. The accelerated life test also provides information on other PRD failure modes that are somewhat statistically random in nature.

33 ADVANCED PROPULSION SYSTEMS↗

A Non-hydrostatic Variable Resolution Atmospheric Model in ACME

The Accelerated Climate Modeling for Energy (ACME) Project is designed to perform coupled climate simulations at high resolution on current and next generation DOE leadership class supercomputers. Prior to this project, the ACME atmospheric model employed the primitive equations of motion, which constrained its maximum resolution and limited the physical phenomena that it could resolve. Our team removed this constraint by implementing a non-hydrostatic (NH) atmospheric model in the ACME Framework. Collaboration with SciDAC FASTMath institutes facilitated improvements to key technologies needed for the non-hydrostatic model. This model was tested in a variable resolution framework, which enables high resolution simulation regions at reasonable computational costs. The correct performance of the NH model was verified with standard dynamical-core tests, and by demonstrating that it reproduced expected climate statistics in the hydrostatic regime. The completed model was applied to evaluate the limitations of existing turbulent mountain stress and gravity wave drag parameterizations, in order to demonstrate the power of the new nonhydrostatic model to model new phenomena and enhance the quality of present day 1/4 degree simulations.

Hall, David Matthew↗

Modeling the Functional Forms of Grid Disturbances

This report introduces a functional form that may be used to quantitatively predict the impacts of new grid tools and changing system qualities on the likelihoods, durations, and depths of various grid disturbances. Each disturbance scenario is modeled to have three component stages—avoidance, reactance and recovery, which together parametrically estimate one disturbance’s impacts. The modeled scenario is then placed and replicated within an analysis period to represent the likelihood or frequency of the scenario and its consequent impacts. Whereas analysts have struggled to define and apply metrics for grid resilience, the functional form introduced by this report shares units of measurement with accepted grid-status measures (e.g., numbers of customers currently experiencing a service outage). Furthermore, the integrated and averaged functional form over an analysis period provides a meaningful normalized performance metric (e.g., customer outage minutes per year) that is ultimately independent of the duration of the period of. The approach may be applied similarly regardless of the severity or frequency of the disturbances that are being analyzed. Because metrics can be chosen to be identical in both the hypothetical future and the actual historical past, the historical past eventually becomes the test of the future predictions, at least in a statistical sense. The authors originally developed this approach to facilitate analysis of the effects of transactive energy (TE) systems effects on electric power grid resilience. TE systems invite energy suppliers and consumers to actively collaborate toward the discovery of, and their responses to, the locational value of energy. The findings from this process are often embodied as energy prices, the dynamics of which indicate the locational value of energy and can further represent important grid service needs. While some academic papers claim to quantify the value of a specific TE system design toward grid resilience, the answer, in general, has been elusive. Not only do multiple and conflicting definitions of resilience and reliability exist, but countless TE systems are being invented. We conclude the following: (1) The ideal analysis should harmonize rather than differentiate resilience and reliability. Therefore, this report uses the more general term disturbance whenever the overloaded terms resilience and reliability can be avoided. (2) The effectiveness of TE systems must be mapped to underlying qualities of a TE system, thereby avoiding presumptions that every TE design offers similar advantages. The authors seek to evaluate the parametric effects of TE system qualities (e.g., spatial granularity, granularity of time steps, length of future prediction horizon) on avoiding, reacting to, and recovering from grid disturbances. Furthermore, any advantages (or disadvantages) must be fairly compared with the many alternative tools, systems, and strategies that might offer comparable benefits.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Improving the Parameterization of Cloud and Rain Microphysics in E3SM using Novel Observationally-Constrained Bayesian Approach (Final Technical Report)

In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.

54 ENVIRONMENTAL SCIENCES↗

Charge States and FIP Bias of the Solar Wind from Coronal Holes, Active Regions, and Quiet Sun

Connecting in situ measured solar-wind plasma properties with typical regions on the Sun can provide an effective constraint and test to various solar wind models. We examine the statistical characteristics of the solar wind with an origin in different types of source regions. We find that the speed distribution of coronal-hole (CH) wind is bimodal with the slow wind peaking at ∼400 km s{sup −1} and the fast at ∼600 km s{sup −1}. An anti-correlation between the solar wind speeds and the O{sup 7+}/O{sup 6+} ion ratio remains valid in all three types of solar wind as well during the three studied solar cycle activity phases, i.e., solar maximum, decline, and minimum. The N{sub Fe}/N{sub O} range and its average values all decrease with the increasing solar wind speed in different types of solar wind. The N{sub Fe}/N{sub O} range (0.06–0.40, first ionization potential (FIP) bias range 1–7) for active region wind is wider than for CH wind (0.06–0.20, FIP bias range 1–3), while the minimum value of N{sub Fe}/N{sub O} (∼ 0.06) does not change with the variation of speed, and it is similar for all source regions. The two-peak distribution of CH wind and the anti-correlation between the speed and O{sup 7+}/O{sup 6+} in all three types of solar wind can be explained qualitatively by both the wave-turbulence-driven and reconnection-loop-opening (RLO) models, whereas the distribution features of N{sub Fe}/N{sub O} in different source regions of solar wind can be explained more reasonably by the RLO models.

79 ASTRONOMY AND ASTROPHYSICS↗