Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical techniques”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Structure–property relationships in atomic-scale junctions: Histograms and beyond

Over the past 10 years, there has been tremendous progress in the measurement, modeling and understanding of structure–function relationships in single molecule junctions. Numerous research groups have addressed significant scientific questions, directed both to conductance phenomena at the single molecule level and to the fundamental chemistry that controls junction functionality. Many different functionalities have been demonstrated, including single-molecule diodes, optically and mechanically activated switches, and, significantly, physical phenomena with no classical analogues, such as those based on quantum interference effects. Experimental techniques for reliable and reproducible single molecule junction formation and characterization have led to this progress. In particular, the scanning tunneling microscope based break-junction (STM-BJ) technique has enabled rapid, sequential measurement of large numbers of nanoscale junctions allowing a statistical analysis to readily distinguish reproducible characteristics. Furthermore, harnessing fundamental link chemistry has provided the necessary chemical control over junction formation, enabling measurements that revealed clear relationships between molecular structure and conductance characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal Experiments for Fractured Rock Characterization: Theoretical Analysis and Inverse Modeling

Abstract Field‐scale properties of fractured rocks play a crucial role in many subsurface applications, yet methodologies for identification of the statistical parameters of a discrete fracture network (DFN) are scarce. We present an inversion technique to infer two such parameters, fracture density and fractal dimension, from cross‐borehole thermal experiments data. It is based on a particle‐based heat‐transfer model, whose evaluation is accelerated with a deep neural network (DNN) surrogate that is integrated into a grid search. The DNN is trained on a small number of the heat‐transfer model runs and predicts the cumulative density function of the thermal field. The latter is used to compute fine posterior distributions of the (to be estimated) parameters. Our synthetic experiments reveal that fracture density is well constrained by data, while fractal dimension is harder to determine. Adding nonuniform prior information related to the DFN connectivity improves the inference of this parameter.

Zhou, Zitong↗
DOE OSTI2014

CHARA/MIRC observations of two M supergiants in Perseus OB1: Temperature, bayesian modeling, and compressed sensing imaging

Two red supergiants (RSGs) of the Per OB1 association, RS Per and T Per, have been observed in the H band using the Michigan Infra-Red Combiner (MIRC) instrument at the CHARA array. The data show clear evidence of a departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened disks based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of RSGs. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3%-5% of the stellar flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Deep inference of simulated strong lenses in ground-based surveys

The large number of strong lenses discoverable in future astronomical surveys will likely enhance the value of strong gravitational lensing as a cosmic probe of dark energy and dark matter. However, leveraging the increased statistical power of such large samples will require further development of automated lens modeling techniques. We show that deep learning and simulation-based inference (SBI) methods produce informative and reliable estimates of parameter posteriors for strong lensing systems in ground-based surveys. We present the examination and comparison of two approaches to lens parameter estimation for strong galaxy-galaxy lenses — Neural Posterior Estimation (NPE) and Bayesian Neural Networks (BNNs). We perform inference on 1-, 5-, and 12-parameter lens models for ground-based imaging data that mimics the Dark Energy Survey (DES). We find that NPE outperforms BNNs, producing posterior distributions that are more accurate, precise, and well-calibrated for most parameters. For the 12-parameter NPE model, the calibration is consistently within <10% of optimal calibration for all parameters, while the BNN is rarely within 20% of optimal calibration for any of the parameters. Similarly, residuals for most of the parameters are smaller (by up to an order of magnitude) with the NPE model than the BNN model. This work takes important steps in the systematic comparison of methods for different levels of model complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Isovector and isoscalar tensor charges of the nucleon from lattice QCD

We present results for the iso-vector and flavor diagonal tensor charges $g$$^{u−d}_{T}$, $g$$^{u}_{T}$, $g$$^{d}_{T}$, and $g$$^{s}_{T}$ needed to probe novel tensor interactions at the TeV scale in neutron and nuclear β-decays and the contribution of the quark electric dipole moment (EDM) to the neutron EDM. The lattice QCD calculations were done using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. These ensembles span three lattice spacings a ≈ 0.06, 0.09 and 0.12 fm and three quark masses corresponding to the pion masses M π ≈ 130, 220 and 310 MeV. Using estimates from these ensembles, we quantify all systematic uncertainties and perform a simultaneous extrapolation in the lattice spacing, volume and light quark masses for the connected contributions. The final estimates of the connected nucleon (proton) tensor charge for the iso-vector combination is $g$$^{u−d}_{T}$ = 1.020(76) in the MS scheme at 2 GeV. The additional disconnected quark loop contributions needed for the flavor-diagonal matrix elements are calculated using a stochastic estimator employing the truncated solver method with the all-mode-averaging technique. We find that the size of the disconnected contribution is smaller than the statistical error in the connected contribution. This allows us to bound the disconnected contribution and include it as an additional uncertainty in the flavor-diagonal charges. After a continuum extrapolation, we find $g$$^{u}_{T}$ = 0.774(66), $g$$^{d}_{T}$ = −0.233(28) and $g$$^{u+d}_{T}$ = 0.541(67). The strangeness tensor charge, that can make a significant contribution to the neutron EDM due to the large ratio ms/mu,d, is $g$$^{s}_{T}$ = 0.008(9) in the continuum limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Nuclear Deformation at Finite Temperature

Deformation, a key concept in our understanding of heavy nuclei, is based on a mean-field description that breaks the rotational invariance of the nuclear many-body Hamiltonian. Here, we present a method to analyze nuclear deformations at finite temperature in a framework that preserves rotational invariance. The auxiliary-field Monte-Carlo method is used to generate the statistical ensemble and calculate the probability distribution associated with the quadrupole operator. Applying the technique to nuclei in the rare-earth region, we identify model-independent signatures of deformation and find that deformation effects persist to higher temperatures than the spherical-to-deformed shape phase-transition temperature of mean-field theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
DOE OSTI2016

SU-C-207A-05: Feature Based Water Equivalent Path Length (WEPL) Determination for Proton Radiography by the Technique of Time Resolved Dose Measurement

Purpose: Studies show that WEPL can be determined from modulated dose rate functions (DRF). However, the previous calibration method based on statistics of the DRF is sensitive to energy mixing of protons due to scattering through different materials (termed as range mixing here), causing inaccuracies in the determination of WEPL. This study intends to explore time-domain features of the DRF to reduce the effect of range mixing in proton radiography (pRG) by this technique. Methods: An amorphous silicon flat panel (PaxScan™ 4030CB, Varian Medical Systems, Inc., Palo Alto, CA) was placed behind phantoms to measure DRFs from a proton beam modulated by a specially designed modulator wheel. The performance of two methods, the previously used method based on the root mean square (RMS) and the new approach based on time-domain features of the DRF, are compared for retrieving WEPL and RSP from pRG of a Gammex phantom. Results: Calibration by T{sub 80} (the time point for 80% of the major peak) was more robust to range mixing and produced WEPL with improved accuracy. The error of RSP was reduced from 8.2% to 1.7% for lung equivalent material, with the mean error for all other materials reduced from 1.2% to 0.7%. The mean error of the full width at half maximum (FWHM) of retrieved inserts was decreased from 25.85% to 5.89% for the RMS and T{sub 80} method respectively. Monte Carlo simulations in simplified cases also demonstrated that the T{sub 80} method is less sensitive to range mixing than the RMS method. Conclusion: WEPL images have been retrieved based on single flat panel measured DRFs, with inaccuracies reduced by exploiting time-domain features as the calibration parameter. The T{sub 80} method is validated to be less sensitive to range mixing and can thus retrieve the WEPL values in proximity of interfaces with improved numerical and spatial accuracy for proton radiography.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Quantifying EGS Reservoir Complexity with an Integrated Geophysical Approach-Improved Resolution Ambient Seismic Noise Interferometry

The primary goal of the project was met, as improved resolution (from kilometers to tens/hundreds of meters and to a depth of at least 3 km) was added to the seismic interferometry-estimated parameters relative to the work reported in AltaRock Energy Inc. (2014) for the Dixie Valley Geothermal Wellfield. A procedure for initial fault identification, using seismic interferometry, was developed. The project team further improved, and applied a new, cost-effective, and non-invasive exploration method, developed and validated since 2010, during two previous experiments, at Dixie Valley and Soda Lake, Nevada. The method, initially designed for hydrothermal exploration, was adapted to EGS exploration conditions, which require deeper crustal information. Regions of elevated energy for both the 2-7 Hz and 5-15 Hz frequency bands are located above intersections of multiple faults. Thus, an unprecedented, innovative seismic field deployment and ambient seismic noise and signal processing were successfully used to estimate crustal structure and to identify fault locations at depths up to 3 km, at a cost per kilometer of survey that in some cases may be as much as ten times less than the cost of a controlled-source survey. Based on valid and useful seismic data collected and analyzed during Phase I, in Phase II the team has further quantified statistically the advantage of using seismic analysis, in combination with other geological and geophysical techniques, for EGS target identification. A power spectral density analysis of the passive-seismic empirical Green’s functions reveals significantly higher power recorded over the exploration area at the power plant, and in the drilling area, near the BB station 10201, especially on the horizontal components. In this study, the most important relationship observed during the prior report, namely the one between velocity and temperature was explored. The areas that show anomalously high predicted temperature in the geostatistics occur where multiple faults intersect. Not surprisingly, such areas span between the Dixie Valley range front fault and the basin-ward Piedmont faults, in the vicinity of numerous wells. These locations also fall in the area of high EGS favorability from the prior 2014 study. The anomalies also coincide with the high power spectral density values. The result of altering all favorability criteria for the new data reached a maximum of 1/8th of a point and were not large enough to change any favorability value reported in the 2014 study. The new data reinforce the results obtained from the previous 2014 study. The current high-resolution ambient seismic noise project provided excellent results in defining faults in the wellfield. As such the identification of known faults by the project demonstrated the utility of this method for identifying faults at an exploration stage in geothermal development. Such fault identification, while extremely useful in exploration investigation, has not added sufficient new data for the generation of significantly different DVGW favorability maps in the current project. We believe the cost savings of our passive-seismic technique for EGS, over common industry active-source survey techniques, may have a minimum of 30%, and a maximum of 90%. There are many significant non-cost benefits of passive seismic surveys. Our passive-source method provides an environmentally-friendly seismic survey with minimal to negligible impacts. The economic benefit is greater EGS project feasibility.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Main Bearing Replacement and Damage - A Field Data Study on 15 Gigawatts of Wind Energy Capacity

This study seeks to establish a comprehensive baseline of knowledge for the replacement and damage of main bearings in wind turbines. The purpose of this report is to provide a high-level summary of the data set, methodology, and results of this work. Full technical details and an extended analysis will be made available in a future publication. We collected data on main bearing replacements and reported damage from industrial partners based in Europe and the United States. In total, we obtained data for 167 wind power plants, with a combined capacity of 15.3 gigawatts (GW). Most of the data set was comprised of land-based, three-point mount, spherical roller bearings. Within this data set were 689 instances of main bearing replacement. Analysis was undertaken in two parts: first, a statistical analysis of the main bearing time-to-replacement data using survival analysis techniques; second, quantitative and qualitative analyses of the obtained damage information. Our results showed that 10% of a fixed main bearing population would be expected to have been replaced by 10.5 years. This is close to half of the 20-year design value. Fitted parametric distributions then indicated that by year 20, between 22% and 25% of main bearings are expected to have been replaced. Analysis of the damage reports revealed spalling to be the main type of damage listed. The additional presence of surface damage in the collected data indicates that at least part of the spalling cases are likely due to surface-initiated rolling contact fatigue. At this stage is not clear what proportion of spalling cases result from "wear induced", surface-initiated and subsurface-initiated rolling contact fatigue. While this work provides important insights into the current state of main bearing replacements and damage, many questions remain. An ongoing and expanding data collection and analysis effort focused on wind turbine main bearings is therefore recommended.

17 WIND ENERGY↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

Makovec, Alajos↗
DOE OSTI2018

Quantitative Constraints on the Reionization History from the IGM Damping Wing Signature in Two Quasars at z > 7

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral universe, and enables measurements of the hydrogen neutral fraction x HI (z) at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Lyα. Here, we combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two quasars with the highest redshifts known, ULAS J1120+0641 and ULAS J1342+0928, and measure volume-averaged neutral fractions {x HI }(Z = 7.09) = 0.48 $+0.26\atop{-0.26}$and {x HI }(Z = 7.54) = 0.60 $+0.20\atop{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of 10 3 ≤ t q ≤ 10 8 yr.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Study of on-board compression of earth resources data

The current literature on image bandwidth compression was surveyed and those methods relevant to compression of multispectral imagery were selected. Typical satellite multispectral data was then analyzed statistically and the results used to select a smaller set of candidate bandwidth compression techniques particularly relevant to earth resources data. These were compared using both theoretical analysis and simulation, under various criteria of optimality such as mean square error (MSE), signal-to-noise ratio, classification accuracy, and computational complexity. By concatenating some of the most promising techniques, three multispectral data compression systems were synthesized which appear well suited to current and future NASA earth resources applications. The performance of these three recommended systems was then examined in detail by all of the above criteria. Finally, merits and deficiencies were summarized and a number of recommendations for future NASA activities in data compression proposed.

Habibi, A.↗
NASA NTRS1978

Ocean waves

Ocean wave data can be obtained from such active microwave probe techniques as monostatic HF and VHF, bistatic HF, HF synthetic aperture radar, altimeters, satellite and airborne synthetic aperture radar, carrier wave or pulsed dual-frequency radars, and coastal surveillance radar. Approaches to texture analysis of ocean wave imagery are discussed, with attention given to transform techniques or spatial frequency analysis, and the analysis of second-order gray level statistics. In addition, recommendations are made for further work on the modulation of short gravity waves by longer waves as a function of wind speed and wave direction, and the derivation of transfer functions for the ocean response of dual-frequency radars.

Bartsch, N.↗
NASA NTRS1978

Pulse compression inside an actively AM mode-locked Nd:YAG laser using a liquid Kerr cell

A mode-locked Nd:YAG laser capable of generating pulses substantially shorter than the present limit of 30-50 ps, while retaining the stability of active mode locking, would be a useful development. This study presents a simplified theoretical analysis and a preliminary experimental verification of the use of an optical Kerr cell as a passive lossless pulse compressor inside a Q-switched and actively mode-locked YAG laser. The theoretical model indicates possible compression to as short as about 5 ps in a typical low-power Nd:YAG laser. The preliminary experimental results show substantial pulse compression via the predicted mechanism, though still not below about 30 ps. The technique appears capable, however, of pulse compression to substantially shorter values without the instabilities, statistical fluctuations, or energy losses associated with saturable-absorber mode locking.

Duerinckx, A. J.↗