Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sufficient 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 199 records · Page 11

Verification of the Uniformly-Ordered Binary Decision Algorithm in Correlated-Benchmark Whisper Calculations

Whisper is a nuclear criticality safety code package that aids analysts in validation exercises by computing upper subcritical limits (USL) for applications of interest. To obtain statistically meaningful, significant, and conservative USLs, the analyst must ensure that Whisper selects a sufficient number of benchmarks that are neutronically similar to the application. Many of the available benchmarks are correlated but are currently treated as independent, leading to an artificially small sample size, as their individual information contributions will be overestimated. To aid the analyst in obtaining a sufficient sample size, prior work [2] demonstrated application of the Uniformly-Ordered Binary Decision (UOBD) algorithm in adjusting benchmark weights to account for benchmark correlations. This work provides verification of the Whisper implementation and considers the impact of updated benchmark correlations compared to those available previously. We demonstrate that the UOBD algorithm performs as expected with an analytic example. With HEU-SOL-THERM-001 cases 1 through 10 as the applications, we compare the USLs computed with benchmark correlations available in the Whisper 1.1 release only to those computed with additional benchmark correlations from DICE 2023 and demonstrate substantive differences.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Challenges for the statistical gravitational-wave method to measure the Hubble constant

Gravitational waves (GW) can be employed as standard sirens that will soon measure the Hubble constant with sufficient precision to weigh in on the ~5σ Hubble tension. Most GW sources will have no identified electromagnetic counterpart, leading to uncertainty in the redshift of the source, and in turn a degeneracy between host galaxy distance, redshift, and H 0 . In the case where no electromagnetic counterparts are identified, it has been proposed that a statistical canvassing of candidate GW hosts, found in a large galaxy survey for example, can be used to accurately constrain the Hubble constant. We study and simulate this "galaxy voting" method to compute H 0 . Here, we find that the Hubble constant posterior is in general biased relative to the true value even when making optimistic assumptions about the statistical properties of the sample. Using the MICECAT light-cone catalog, we find that the bias in the H0 posteriors depends on the realization of the underlying galaxy sample and the precision of the GW source distance measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗

Implementation and Assessment of Advanced Analog Vector-Matrix Processor

This paper discusses the design and implementation of an analog optical vecto-rmatrix coprocessor with a throughput of 128 Mops for a personal computer. Vector matrix calculations are inherently parallel, providing a promising domain for the use of optical calculators. However, to date, digital optical systems have proven too cumbersome to replace electronics, and analog processors have not demonstrated sufficient accuracy in large scale systems. The goal of the work described in this paper is to demonstrate a viable optical coprocessor for linear operations. The analog optical processor presented has been integrated with a personal computer to provide full functionality and is the first demonstration of an optical linear algebra processor with a throughput greater than 100 Mops. The optical vector matrix processor consists of a laser diode source, an acoustooptical modulator array to input the vector information, a liquid crystal spatial light modulator to input the matrix information, an avalanche photodiode array to read out the result vector of the vector matrix multiplication, as well as transport optics and the electronics necessary to drive the optical modulators and interface to the computer. The intent of this research is to provide a low cost, highly energy efficient coprocessor for linear operations. Measurements of the analog accuracy of the processor performing 128 Mops are presented along with an assessment of the implications for future systems. A range of noise sources, including cross-talk, source amplitude fluctuations, shot noise at the detector, and non-linearities of the optoelectronic components are measured and compared to determine the most significant source of error. The possibilities for reducing these sources of error are discussed. Also, the total error is compared with that expected from a statistical analysis of the individual components and their relation to the vector-matrix operation. The sufficiency of the measured accuracy of the processor is compared with that required for a range of typical problems. Calculations resolving alloy concentrations from spectral plume data of rocket engines are implemented on the optical processor, demonstrating its sufficiency for this problem. We also show how this technology can be easily extended to a 100 x 100 10 MHz (200 Cops) processor.

Gary, Charles K.↗

Paleomagnetism of the Mesozoic in Alaska

Over 400 oriented cores of Permian, Triassic, Jurassic, and Cretaceous sedimentary and igneous rocks were collected from 34 sites at 10 areas throughout southern Alaska. After magnetic cleaning in successively higher alternating fields 179 samples were considered to be stable and to give statistically consistent results within each site and age group. Due to the lack of a sufficient number of stable samples, the results from Permian, Triassic, and Cretaceous rocks were inconclusive. The nine remaining Jurassic sites represent 100 samples from three general areas in southern Alaska. The southern Alaskan Jurassic paleomagnetic pole is significantly different from the North American Jurassic pole. This suggests that since the Jurassic, southern Alaska must have moved approximately 18 degrees north and rotated 52 degrees clockwise to reach its present position. Tectonic interpretation of these results give a possible explanation for many of the geologic features observed in southern Alaska.

Packer, D. R.↗

Service Life Extension of the Propulsion System of Long-Term Manned Orbital Stations

One of the critical non-replaceable systems of a long-term manned orbital station is the propulsion system. Since the propulsion system operates beginning with the launch of station elements into orbit, its service life determines the service life of the station overall. Weighing almost a million pounds, the International Space Station (ISS) is about four times as large as the Russian space station Mir and about five times as large as the U.S. Skylab. Constructed over a span of more than a decade with the help of over 100 space flights, elements and modules of the ISS provide more research space than any spacecraft ever built. Originally envisaged for a service life of fifteen years, this Earth orbiting laboratory has been in orbit since 1998. Some elements that have been launched later in the assembly sequence were not yet built when the first elements were placed in orbit. Hence, some of the early modules that were launched at the inception of the program were already nearing the end of their design life when the ISS was finally ready and operational. To maximize the return on global investments on ISS, it is essential for the valuable research on ISS to continue as long as the station can be sustained safely in orbit. This paper describes the work performed to extend the service life of the ISS propulsion system. A system comprises of many components with varying failure rates. Reliability of a system is the probability that it will perform its intended function under encountered operating conditions, for a specified period of time. As we are interested in finding out how reliable a system would be in the future, reliability expressed as a function of time provides valuable insight. In a hypothetical bathtub shaped failure rate curve, the failure rate, defined as the number of failures per unit time that a currently healthy component will suffer in a given future time interval, decreases during infant-mortality period, stays nearly constant during the service life and increases at the end when the design service life ends and wear-out phase begins. However, the component failure rates do not remain constant over the entire cycle life. The failure rate depends on various factors such as design complexity, current age of the component, operating conditions, severity of environmental stress factors, etc. Development, qualification and acceptance test processes provide rigorous screening of components to weed out imperfections that might otherwise cause infant mortality failures. If sufficient samples are tested to failure, the failure time versus failure quantity can be analyzed statistically to develop a failure probability distribution function (PDF), a statistical model of the probability of failure versus time. Driven by cost and schedule constraints however, spacecraft components are generally not tested in large numbers. Uncertainties in failure rate and remaining life estimates increase when fewer units are tested. To account for this, spacecraft operators prefer to limit useful operations to a period shorter than the maximum demonstrated service life of the weakest component. Running each component to its failure to determine the maximum possible service life of a system can become overly expensive and impractical. Spacecraft operators therefore, specify the required service life and an acceptable factor of safety (FOS). The designers use these requirements to limit the life test duration. Midway through the design life, when benefits justify additional investments, supplementary life test may be performed to demonstrate the capability to safely extend the service life of the system. An innovative approach is required to evaluate the entire system, without having to go through an elaborate test program of propulsion system elements. Evaluating every component through a brute force test program would be a cost prohibitive and time consuming endeavor. ISS propulsion system components were designed and built decades ago. There are no representative ground test articles for some of the components. A 'test everything' approach would require manufacturing new test articles. The paper outlines some of the techniques used for selective testing, by way of cherry picking candidate components based on failure mode effects analysis, system level impacts, hazard analysis, etc. The type of testing required for extending the service life depends on the design and criticality of the component, failure modes and failure mechanisms, life cycle margin provided by the original certification, operational and environmental stresses encountered, etc. When specific failure mechanism being considered and the underlying relationship of that mode to the stresses provided in the test can be correlated by supporting analysis, time and effort required for conducting life extension testing can be significantly reduced. Exposure to corrosive propellants over long periods of time, for instance, lead to specific failure mechanisms in several components used in the propulsion system. Using Arrhenius model, which is tied to chemically dependent failure mechanisms such as corrosion or chemical reactions, it is possible to subject carefully selected test articles to accelerated life test. Arrhenius model reflects the proportional relationship between time to failure of a component and the exponential of the inverse of absolute temperature acting on the component. The acceleration factor is used to perform tests at higher stresses that allow direct correlation between the times to failure at a high test temperature to the temperatures to be expected in actual use. As long as the temperatures are such that new failure mechanisms are not introduced, this becomes a very useful method for testing to failure a relatively small sample of items for a much shorter amount of time. In this article, based on the example of the propulsion system of the first ISS module Zarya, theoretical approaches and practical activities of extending the service life of the propulsion system are reviewed with the goal of determining the maximum duration of its safe operation.

Kamath, Ulhas↗

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗

Statistical significance test for transition matrices of atmospheric Markov chains

Low-frequency variability of large-scale atmospheric dynamics can be represented schematically by a Markov chain of multiple flow regimes. This Markov chain contains useful information for the long-range forecaster, provided that the statistical significance of the associated transition matrix can be reliably tested. Monte Carlo simulation yields a very reliable significance test for the elements of this matrix. The results of this test agree with previously used empirical formulae when each cluster of maps identified as a distinct flow regime is sufficiently large and when they all contain a comparable number of maps. Monte Carlo simulation provides a more reliable way to test the statistical significance of transitions to and from small clusters. It can determine the most likely transitions, as well as the most unlikely ones, with a prescribed level of statistical significance.

Vautard, Robert↗

A Data-Driven Framework for Power System Event Type Identification via Safe Semi-Supervised Techniques

Herein this paper investigates the use of phasor measurement unit (PMU) data with deep learning techniques to construct real-time event identification models for transmission networks. Increasing penetration of distributed energy resources represents a great opportunity to achieve decarbonization, as well as challenges in systematic situational awareness. When high-resolution PMU data and sufficient manually recorded event labels are available, the power event identification problem is defined as a statistical classification problem that can be solved by numerous cutting-edge classifiers. However, in real grids, collecting tremendous high-quality event labels is quite expensive. Utilities frequently have a large number of event records without in-depth details (i.e., unlabeled events). To bridge this gap, we propose a novel semi-supervised learning-based method to improve the performance of event classifiers trained with a limited number of labeled events by exploiting the information from massive unlabeled events. In other words, compared to existing data-driven methods, our method requires only a small portion of labeled data to achieve a similar level of accuracy. Meanwhile, this work discusses and addresses the performance degradation caused by class distribution mismatch between the training set and the real applications. Based on the proposed safe learning mechanism, our model does not directly use all unlabeled events during model training, but selectively uses them through a comprehensive evaluation procedure. Numerical studies on a sizable PMU dataset have been used to validate the performance of the proposed method.

42 ENGINEERING↗

In situ measurement of the cosmic microwave background temperature at a distance of 7.5 kiloparsecs

Absorption of orthoformaldehyde (H2CO) toward the giant H II region W51A (G49.5 - 0.4) has been used as a remote probe of the cosmic microwave background (CMB) at a distance of 7.5 kpc. VLA observations of the 6 cm and 2 cm wavelength transitions provided sufficient resolution and sensitivity to resolve condensations within the foregound clouds. Solutions to the equations of statistical equilibrium within each condensation in the context of a large velocity gradient model yielded an estimate for the CMB temperature at 2.1 mm wavelength of 3.2 + or - 0.9 K. The uncertainty is dominated by modeling of collisional pumping by neutral hydrogen molecules (H2).

Kogut, A.↗

Using Global Simulations of the Magnetosphere for Multi-Satellite Mission Planning and Data Analysis

We use global simulations of Earth's magnetosphere to assess the scientific return from a multi-satellite mission in the magnetosphere. We examine 4 different scenarios with 20, 40, 80, and 160 satellites, respectively. The satellite orbits are randomized with perigee distances ranging from 2 to 5R(sub E), apogee distances between 10 and 50 R(sub E), and within +/-5R(sub E) of the geocentric solar ecliptic (GSE) equator. For each of these satellite configurations we examine the expected observations during a typical substorm by using time traces obtained from a global simulation at the satellite positions. The 160 satellite configuration yields sufficient information to distinguish between different substorm models without any temporal/spatial ambiguities. An 80 satellite configuration still provides sufficient information for this task, however for fewer events with good satellite conjunctions and with less statistical certainty. For constellations with fewer than 40 satellites time-space ambiguities are likely to remain in the observation. However, any multi-satellite constellation would be a quantum leap in magnetospheric research because of the unprecedented coverage of other regions, because it would enable new measurement techniques that are unique to multi-satellite missions, and because it would enable the use of data assimilation techniques in global models for the first time.

Raeder, Joachim↗

Wind Retrievals under Rain for Passive Satellite Microwave Radiometers and its Applications to Hurricane Tracking

We have developed an algorithm that retrieves wind speed under rain using C-hand and X-band channels of passive microwave satellite radiometers. The spectral difference of the brightness temperature signals due to wind or rain allows to find channel combinations that are sufficiently sensitive to wind speed but little or not sensitive to rain. We &ve trained a statistical algorithm that applies under hurricane conditions and is able to measure wind speeds in hurricanes to an estimated accuracy of about 2 m/s. We have also developed a global algorithm, that is less accurate but can be applied under all conditions. Its estimated accuracy is between 2 and 5 mls, depending on wind speed and rain rate. We also extend the wind speed region in our model for the wind induced sea surface emissivity from currently 20 m/s to 40 mls. The data indicate that the signal starts to saturate above 30 mls. Finally, we make an assessment of the performance of wind direction retrievals from polarimetric radiometers as function of wind speed and rain rate

Meissner, Thomas↗

Spatial and Temporal Variability of Vertical Velocity under Shallow Cumulus

Vertical velocity distribution below cloud is one of the key determinants of cloud life cycle, but observations of this variable are extremely sparse in space. Doppler lidar retrievals and large-eddy simulations at the U.S. Department of Energy’s Atmospheric Radiation Measurement User Facility Southern Great Plains site are used to determine whether vertical velocity statistics from temporally dense profiles at a single location can be substituted for spatial vertical velocity statistics. We show that even a small number (five) of widely distributed [ O (1°) latitude/longitude spacing] lidars is sufficient sampling to reconstruct domainwide spatial vertical velocity variance, but not higher moments of the vertical velocity distribution. Spatial and temporal vertical velocity variances in the Doppler lidar observations are nearly interchangeable as long as the spatial variance is temporally averaged and the temporal variance is averaged across lidars. This is true even though the dominant spatial scales of vertical velocity variability are ≲ 3 km, more than an order of magnitude smaller than the spacing between the lidars. Further, in the limit where the temporal variance does not vary across a spatial domain (e.g., if the meteorological and surface forcing of the atmospheric turbulence is homogeneous across the domain) and the domain-mean vertical velocity is zero, the commonly available retrieval of temporal vertical velocity variance at one site is equivalent to the spatial variance over the domain. We use an updraft parcel model to show that substituting temporal for spatial vertical velocity statistics will have a relatively minor effect on cloud droplet number concentrations.

54 ENVIRONMENTAL SCIENCES↗

Clustering and runaway merging in a primordial black hole dominated universe

If primordial black holes (PBHs) are present in the early Universe, their contribution to the energy budget grows relative to that of radiation and generically becomes dominant unless the initial abundance is exponentially small. This black hole domination scenario is largely unconstrained for PBHs with masses ≲ 10 9 g, which evaporate prior to big bang nucleosynthesis. However, if the era of PBH domination is sufficiently long, the PBHs form clusters and can merge appreciably within these objects. We calculate the population statistics of these clusters within the Press-Schechter formalism and find that, for a wide range of PBH masses and Hubble rates at the onset of PBH domination, the mergers within PBH clusters can exhibit runaway behavior, where the majority of the cluster will eventually form a single black hole with a mass much greater than the original PBH mass. Finally, these mergers can dramatically alter the PBH mass distribution and leave behind merged relic black holes that evaporate after big bang nucleosynthesis and yield various observational signatures, excluding parameter choices previously thought to be viable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Correcting the PFNS for more consistent fission modeling

For FY20, we had a deliverable to write a report detailing efforts to simultaneously evaluate both the prompt neutron multiplicity, $\overline{ν}$, and the prompt neutron fission neutron spectrum, PFNS, using CGMF. CGMF is the LANL-developed fission fragment decay code that consistently evaporates prompt neutrons and γ rays using the Hauser-Feshbach statistical theory of compound reactions. The decay begins by constructing the initial conditions of the fission fragments, then decaying each one from the excited state by neutrons and γ rays, conserving energy, momentum, spin, and parity in each step of the emission. The initial conditions of the fragments, along with the multiplicity, energy, and direction of each emitted neutron or γ ray, are recorded, allowing for the full reconstruction of the fission event. These event histories allow us to reconstruct average quantities, as well as correlations between observables, that can be compared with experimental or evaluated data. In that initial report, although there was already a favorable comparison between $\overline{ν}$ from CGMF, experiment, and the current ENDF/B-VIII.0 evaluation, we showed that there was still significant work to be done to improve the PFNS from CGMF. Historically, the PFNS is calculated too soft by Hauser-Feshbach fission models, and CGMF is no exception. The incorrect shape presents a significant challenge in fission modeling, including for our understanding of the fission process and for our ability to consistently calculate and predict a variety of prompt fission observables (such as fission fragment initial conditions, neutron and γ-ray multiplicities and energies, and the correlations between all observables). In our companion report, we detail our success in using CGMF to evaluate $\overline{ν}$. Although not included in the optimization explicitly, we also keep the initial conditions of the fission fragments physical, along with reproducing reasonably well the neutron multiplicity distribution. As we would expect from the sensitivities calculations from, the average neutron energies change very little from the $\overline{ν}$ optimization along with the PFNS (as will be shown in Sec. 2.6). The conclusion was that the global and statistical models would have to be investigated instead of just the fission fragment initial conditions (as is sufficient for $\overline{ν}$). This report details those efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nationwide forestry applications program. Ten-Ecosystem Study (TES) site 8, Grays Harbor County, Washington

The author has identified the following significant results. Level 2 forest features (softwood, hardwood, clear-cut, and water) can be classified with an overall accuracy of 71.6 percent plus or minus 6.7 percent at the 90 percent confidence level for the particular data and conditions existing at the time of the study. Signatures derived from training fields taken from only 10 percent of the site are not sufficient to adequately classify the site. The level 3 softwood age group classification appears reasonable, although no statistical evaluation was performed.

Prill, J. C.↗

Physical conditions in the cool parts of prominences

A spectral analysis of the cool parts of seven prominences reported by Shin-Huei (1961) and Yakovkin and Zel'dina (1964) is presented. The analysis requires only measurements of the Na and Sr(+) resonance lines and the middle H Balmer lines, all of which are optically thin, and uses the Na and Sr(+) statistical equilibrium solutions of Landman (1983). It is shown that the Na/H integrated intensity ratios alone are sufficient to determine the degree of H ionization. The major result of this investigation is that the degree of H ionization in the cool prominence parts is very low and the total gas pressure very high.

Landman, D. A.↗

Remarks on cometary astrometry

A statistical examination of the astrometric observations made during the past 20 years indicates that there are generally sufficient data on bright and well-publicized comets. On the other hand, astrometric coverage of fainter objects has recently become a severe problem. In spite of an increase in the number of participating observatories, fewer observatories are making significant contributions. Automated methods have generally improved the quality of observations. The timely communication of observations is also discussed.

Marsden, B. G.↗