Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function Estimation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Calculations of electronic excitation by protons and alpha particles in silicon

This paper presents an ab initio method of calculating the valence electron contribution to the stopping and interaction cross sections of swift ions in matter using linear response calculations of the dielectric function in crystals. We estimate the contribution to the stopping and cross sections by repeated electron capture and stripping of the projectile ion by the medium. We also predict the electronic stopping of protons and alpha particles in silicon and compare to measured values. The predictions generally agree with the data at ion speeds greater than 1 atomic unit. We find that electron capture and stripping is at most a small contribution to the stopping for both ions studied.

electron excitation, silicon, nuclear interactions↗
DOE OSTI2014

Lattice thermal conductivity of Bi, Sb, and Bi-Sb alloy from first principles

Using first principles, we calculate the lattice thermal conductivity of Bi, Sb, and Bi-Sb alloys, which are of great importance for thermoelectric and thermomagnetic cooling applications. Our calculation reveals that the ninth-neighbor harmonic and anharmonic force constants are significant; accordingly, they largely affect the lattice thermal conductivity. Several features of the thermal transport in these materials are studied: (1) the relative contributions from phonons and electrons to the total thermal conductivity as a function of temperature are estimated by comparing the calculated lattice thermal conductivity to the measured total thermal conductivity, (2) the anisotropy of the lattice thermal conductivity is calculated and compared to that of the electronic contribution in Bi, and (3) the phonon mean free path distributions, which are useful for developing nanostructures to reduce the lattice thermal conductivity, are calculated. In conclusion, the phonon mean free paths are found to range from 10 to 100 nm for Bi at 100 K.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Uncertainty-quantification-enabled inversion of nuclear responses

Nuclear quantum many-body methods rely on integral transform techniques to infer properties of electroweak response functions from ground-state expectation values. Retrieving the energy dependence of these responses is highly nontrivial, especially for quantum Monte Carlo methods, as it requires inverting the Laplace transform, a notoriously ill-posed problem. Here, in this work, we propose an artificial neural network architecture suitable for accurate response function reconstruction with precise estimation of the uncertainty of the inversion. We demonstrate the capabilities of this new architecture benchmarking it against maximum entropy and previously developed neural network methods designed for a similar task, paying particular attention to its robustness noise in the Euclidean

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

An iterative reconstruction algorithm for pole figure inversion using total variation regularization

This paper describes a new discrete method for inverting X-ray pole figures to estimate the orientation distribution function (ODF). The method employs the equal-volume `cubochoric' representation for uniform discretization of orientation space, SO(3). The forward-projection model is combined with an anisotropic total variation term to iteratively determine the ODF. Here, the efficacy of the new method is evaluated with both model and experimental data and compared with existing discrete and series expansion methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Generalized Graph Laplacian Based Anomaly Detection for Spatiotemporal MicroPMU Data

This letter develops a novel anomaly detection method using the generalized graph Laplacian (GGL) matrix to visualize the spatiotemporal relationship of distribution-level phasor measurement unit (uPMU) data. The uPMU data in a specific time horizon is segregated into multiple segments. An optimization problem formulated as a Lagrangian function is utilized to estimate the GGL matrix. During the iterative process, an optimal update is constituted as a quadratic program (QP) problem. To perform the uPMU-based spatiotemporal analysis, normalized diagonal elements of GGL matrix are proposed as a quantitative metric. The effectiveness of the developed method is demonstrated through real-world uPMU measurements gathered from test feeders in Riverside, CA.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Real-time deep-learning inversion of seismic full waveform data for CO 2 saturation and uncertainty in geological carbon storage monitoring

Deep-learning inversion has recently drawn attention in geological carbon storage research due to its potential of imaging and monitoring carbon storage in real time, significantly improving efficiency and safety of carbon storage operations. We present a deep-learning full waveform inversion method that after the neural network has been trained can image CO 2 saturation and its uncertainty in real time. Our deep-learning inversion method is based on the U-Net architecture with the neural network trained on pairs of synthetic seismic data and CO 2 saturation models. Accordingly, our training establishes a mapping relationship between seismic data and CO 2 saturation models and once fully trained directly estimates CO 2 saturation as a function of subsurface location. We further quantify uncertainties of CO 2 saturation estimates using the Monte Carlo dropout method and a bootstrap aggregating method. For this proof-of-concept study, the CO 2 training models and data are derived from the Kimberlina 1.2 model, a hypothetical 3D geological carbon storage model that is constructed based on various geological and hydrological data from the Southern San Joaquin Basin, California. We perform deep-learning inversion experiments using noise-free and noisy training and test data sets and compare the results. Our modelling experiments show that (1) the deep-learning inversion can estimate 2D distributions of CO 2 fairly well even in the presence of Gaussian random noise and (2) both CO 2 saturation imaging and uncertainty quantification can be done in real time. Our results suggest that the deep-learning inversion method can serve as a robust real-time monitoring tool for geological carbon storage and/or other time-varying reservoir/aquifer properties that result from injection, extraction, and/or other subsurface transport phenomena.

58 GEOSCIENCES↗
DOE OSTI2025

Tidal–hydrological dynamics of water temperature across freshwater forested wetlands on the northeastern Pacific coast

Abstract Tidal freshwater forests were once extensive across temperate coastlines, but loss and fragmentation have made estimation of their ecosystem functions challenging. We measured water temperature for 2 years in three Sitka spruce tidal forests, a restoration site, and an adjacent emergent marsh on the Columbia River, Washington, United States. We assessed spatial variability of water temperature within sites including the effects of hydrology, differences among bay and tributary tidal forests, and differences between the tidal forests and the mainstem Columbia, the restoration site, and the emergent marsh. The tidal forests nearest to the bay had lower interior water temperatures than their channel confluences by up to 2.5°C (weekly median temperature) and 2.0°C (weekly maximum temperature), with most cooling occurring during the low‐flow months of July–September. Tributary sites had maximum temperatures up to 1.9°C cooler than bay sites and 4.2°C cooler than the mainstem. Temperatures in the two bay sites decreased by −0.16°C/100 m and −0.07°C/100 m, on average. The restoration site had the smallest within‐site temperature gradient. Differences in maximum temperatures were greatest when tidal range was low, while higher tidal ranges were associated with warmer and more variable site interiors relative to their confluences. These results suggest that water temperatures in these tidal forests can provide temperature refugia for cold water biota including salmon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Two Methods for Retrieving UV Index for All Cloud Conditions from Sky Imager Products or Total SW Radiation Measurements

In the present study, we assess the cloud effects on UV Index (UVI) and total solar radiation (TR) as a function of cloud cover estimations and sunny conditions (from sky imaging products) as well as of solar zenith angle (SZA). These analyses are undertaken for a southern-hemisphere mid-latitude site where a 10-years dataset is available. It is confirmed that clouds reduce TR more than UV, in particular for obscured Sun conditions, low cloud fraction (< 60%) and large SZA (> 60º). Similarly, clouds enhance TR more than UV, mainly for visible Sun conditions, large cloud fraction and large SZA. Two methods to estimate UVI are developed: 1) from sky imaging cloud cover and sunny conditions, and 2) from TR measurements. Both methods may be used in practical operational applications, although Method 2 shows overall the best performance, since TR allows accounting for cloud optical properties. The mean absolute differences of Method 2 estimations with respect to measured values are 0.17 UVI units (for 1-minute data) and 0.79 Standard Erythemal Dose (SED) units (for daily integrations). Method 1 shows less accurate results but it is still suitable to estimate UVI: mean absolute differences are 0.37 UVI units and 1.6 SED.

clouds↗
DOE OSTI2026

Characterizing degraded visual environments with scanning lenses to improve contrast recovery

Atmospheric aerosols, such as fog and smoke, degrade optical signals via the scattering and absorption of light, making it difficult to recover information from the surrounding environment. Here, we present a technique that measures the angular distribution of light in scattering environments to help recover these degraded optical signals. This is done by passing a collimated beam through the aerosol environment and collecting the scattered and unscattered light with an f-theta scan lens. The lens transforms the angular distribution of the scattered light into the linear domain, mapping the collected light onto a detector. The measured angular scattering distribution is then used to estimate the point spread function of the environment at an arbitrary stand-off distance. This point spread function can be used to deconvolve images blurred by the aerosol environment in which the angular scattering distribution was initially measured. With this approach, we demonstrate improvements in the effective resolution of deblurred images up to 26% over a 120-cm stand-off distance. Experimentation was performed in the Sandia National Laboratories tabletop fog chamber, a platform that generates aerosol environments analogous to real-world fogs. We intend to leverage this technique in future work to develop a low-SWaP system for the in situ characterization of aerosol environments.

Aerosols↗
DOE OSTI2018

Reproducing heavy-ion fusion cross sections at extreme sub-barrier energies with a simple formula: Analysis of fusion hindrance with a Gaussian barrier distribution

Heavy-ion fusion hindrance occurs at extreme sub-barrier energies. This behavior is well reproduced with a simple cross section formula, which was developed by Siwek-Wilczynska et al., based on a single-Gaussian distribution of fusion barrier heights, before the discovery of the hindrance phenomenon. This expression has not yet been widely used and referenced in the literature. An analysis by using this simple formula is presented for 29 systems, from 16 O + 18 O to 64 Ni + 124 Sn , all being measured down to less than 10μb. The agreement with the data is even better than the ones from sophisticated Coupled-channels calculations. This simple expression also applies to fusion reactions in lighter systems. Here, the three parameters contained in this formula vary in a relatively smooth fashion over the whole mass range, and can be used to extrapolate cross sections or to obtain an estimate of the excitation function for systems which have not been measured. Extensions and restrictions of this method are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A multiresolution adaptive wavelet method for nonlinear partial differential equations

We report the multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electron-deuteron deep-inelastic scattering with spectator tagging at EIC

We summarize the phenomenological description of electron-deuteron deep-inelastic scattering (DIS) with detection of a nucleon in the nuclear fragmentation region, eD → e' + N + X ("spectator tagging"), at collision energies accessible with an Electron-Ion Collider (EIC), √s eD ~ several 10 GeV. This includes (a) the kinematic variables, differential cross section, conditional structure functions, and phase space elements; (b) a cross section model based on the light-front impulse approximation and an analytic parametrization of the deuteron wave function; (c) counting rate estimates for simulations of cross section measurements and on-shell extrapolation in the recoil nucleon momentum; (d) a FORTRAN 77 code package for computing the cross section and performing counting rate estimates. The document summarizes the definitions and results without derivation and is intended as a reference for programming and simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Estimating Carbon Flux and Storage: Constraint of the Community Land Model Using Observations at Different Temporal Scales

Terrestrial ecosystems deliver significant global services through their role in the global carbon cycle, energy balance and feedbacks to the climate system. Model projections generally agree under current ‘equilibrium’ conditions, but, despite extensive empirical research, current ESMs disagree on the long-term effects of climatic change and disturbance on the carbon cycle in part because of knowledge gaps and in part because empirical information is not readily incorporated into models. One way to reduce uncertainty in ESMs is to force them to agree with our best estimates of how ecosystems function. Major activities: We 1) developed and tested phenology and biomass data products, 2) tested which datasets were most suitable for data assimilation, 3) developed the data assimilation approaches at a single site and then expanded this to the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Study Results and Plans for Future Improvement of the NSLS-II FOFB System

Currently, the fast orbit feedback (FOFB) system correct the transverse orbit in both planes via air coil magnets installed over low conductivity sections of vacuum chamber (stainless steel bellows) with a fast acquisition (FA) that samples at 10kHz. Bandwidth of a few kHz has been achieved (power-supply + corrector + vacuum chamber) for the purpose that the FOFB bandwidth limitation is not anymore limited by the correctors. The bandwidths of the FOFB system in both planes have been experimentally characterized in the open-loop condition via sinwave excitation method. They are significantly lower than the measurements (2.5kHz in x) done by Weixing Chen et al when the control system was largely bypassed, only from PSI to fast correctors, final to BPM. The estimated dI/dt as a function of the frequency from the beam noise spectrum has been compared to the slew rate of the power supply. The most recent study evidences: after the optimization of fast corrector compensation circuit, the power supply system along shows the 3dB bandwidth of 5.6kHz in the horizontal plane and 3.5kHz in the vertical plane. Preliminary beam studies indicate the potential increase of the system bandwidth and the gain within the FOFB bandwidth. Finally, the plan for combing the FOFB and SOFB systems to correct down to DC has been discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗
DOE OSTI2024

Temporal Explosion Source Processes of Declared Nuclear Tests in the Democratic People’s Republic of Korea

In this work we highlight a preliminary temporal source analysis of the six declared Democratic People's Republic of Korea (DPRK) nuclear tests. We use regional seismic data to estimate relative source time functions (RSTFs) via iterative time-domain deconvolution (Ammon, 2006; Pippin, 2022) of vertical-component ground motions recorded within 2000 km of the source region. Since RSTFs are ideally independent of site and propagation effects, their amplitude spectrum is equivalent to the source spectral ratio, but they also retain phase information. We compare observed RSTFs (in the time and frequency domains) with synthetic RSTFs derived from the Mueller & Murphy (1971) explosion source model. The resolution of these time functions varies, however, we generally obtain high-quality results within the limitations of the recording broadband instrumentation. The results indicate that this method effectively preserves source time-history information that can be used for temporal analysis of remote nuclear explosions. This preliminary analysis is intended to assess the viability of using time-domain deconvolution methods for extracting temporal source information.

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

A Physical Model for Understanding the Activation of MoS 2 Basal-plane Sulfur Atoms for the Hydrogen Evolution Reaction

All the DFT calculations are done with the Vienna Ab Initio Simulation Package (VASP) using the projector augmented wave method. The Bayesian error estimation exchange-correlation functionals (BEEF) with van der Waals interactions are employed. A plane-wave cutoff energy of 400 eV is used together with PAW-PBE potentials where semi core p states are treated as valence. All the calculations allow for spin-polarization. The structures are relaxed until the force is converged to < 0.01 eV/Å. The lattice parameter of MoS 2 unit cell, optimized with this functional, is 3.19 Å. A (4×4) supercell is used to model all the transition metal doped MoS 2 systems studied here, including those with S vacancies. For calculations in the initial dopant structure exploration, the Brillouin zone is sampled with a 3x3x1 Monkhorst-Pack k-point mesh. A 6×6×1 Monkhorst-Pack k-point mesh is used for the H binding energy and density of states calculations. In all calculations, the vacuum layer is set as 15 Å to eliminate periodic interaction perpendicular to the basal plane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Globally optimal interferometry with lossy twin Fock probes

Parity or quadratic spin (e.g., J z 2 ) readouts of a Mach–Zehnder (MZ) interferometer probed with a twin Fock (TF) input state allow saturating the optimal sensitivity attainable among all mode-separable states with a fixed total number of particles but only when the interferometer phase θ is near zero. When more general Dicke state probes are used, the parity readout saturates the quantum Fisher information (QFI) at θ = 0, whereas better-than-standard quantum limit performance of the J z 2 readout is restricted to an o ( N ) occupation imbalance. We show that a method of moments readout of two quadratic spin observables J z 2 and J + 2 + J − 2 is globally optimal for Dicke state probes; i.e., the error saturates the QFI for all θ . In the lossy setting, we derive the time-inhomogeneous Markov process describing the effect of particle loss on TF states, showing that the method of moments readout of four at-most-quadratic spin observables is sufficient for globally optimal estimation of θ when two or more particles are lost. The analysis culminates in a numerical calculation of the QFI matrix for distributed MZ interferometry on the four-mode state | N 4 , N 4 , N 4 , N 4 〉 and its lossy counterparts, showing that an advantage for the estimation of any linear function of the local MZ phases θ 1 and θ 2 (compared to independent probing of the MZ phases by two copies of | N 4 , N 4 〉 ) appears when more than one particle is lost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗