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 1,009 records · Page 56

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Theoretical Determination of Size Effects in Zeolite-Catalyzed Alcohol Dehydration

In the upgrading of biomass pyrolysis vapors to hydrocarbons, dehydration accomplishes a primary objective of removing oxygen, and acidic zeolites represent promising catalysts for the dehydration reaction. Here, we utilized density functional theory calculations to estimate adsorption energetics and intrinsic kinetics of alcohol dehydration over H-ZSM-5, H-BEA, and H-AEL zeolites. The ONIOM (our Own N-layered Integrated molecular Orbital and molecular Mechanics) calculations of adsorption energies were observed to be inconsistent when benchmarked against QM (Quantum Mechanical)/Hartree–Fock and periodic boundary condition calculations. However, reaction coordinate calculations of adsorbed species and transition states were consistent across all levels considered. Comparison of ethanol, isopropanol (IPA), and tert-amyl alcohol (TAA) over these three zeolites allowed for a detailed examination of how confinement impacts on reaction mechanisms and kinetics. The TAA, seen to proceed via a carbocationic mechanism, was found to have the lowest activation barrier, followed by IPA and then ethanol, both of which dehydrate via a concerted mechanism. Barriers in H-BEA were consistently found to be lower than in H-ZSM-5 and H-AEL, attributed to late transition states and either elevated strain or inaccurately estimating long-range electrostatic interactions in H-AEL, respectively. Molecular dynamics simulations revealed that the diffusivity of these three alcohols in H-ZSM-5 were significantly overestimated by Knudsen diffusion, which will complicate experimental efforts to develop a kinetic model for catalytic fast pyrolysis.

09 BIOMASS FUELS↗

( Hubble Space Telescope Hubble Space Telescope) proper motion (HSTPROMO) catalogs of galactic globular clusters. IV. Kinematic profiles and average masses of blue straggler stars

We make use of the Hubble Space Telescope proper-motion catalogs derived by Bellini et al. to produce the first radial velocity dispersion profiles σ(R) for blue straggler stars (BSSs) in Galactic globular clusters (GCs), as well as the first dynamical estimates for the average mass of the entire BSS population. We show that BSSs typically have lower velocity dispersions than stars with mass equal to the main-sequence turnoff mass, as one would expect for a more massive population of stars. Since GCs are expected to experience some degree of energy equipartition, we use the relation σ∝M{sup −η}, where η is related to the degree of energy equipartition, along with our velocity dispersion profiles to estimate BSS masses. We estimate η as a function of cluster relaxation from recent Monte Carlo cluster simulations by Bianchini et al. and then derive an average mass ratio M{sub BSS}/M{sub MSTO}=1.50±0.14 and an average mass M{sub BSS}=1.22±0.12 M {sub ⊙} from 598 BSSs across 19 GCs. The final error bars include any systematic errors that are random between different clusters, but not any potential biases inherent to our methodology. Our results are in good agreement with the average mass of M{sub BSS}=1.22±0.06 M {sub ⊙} for the 35 BSSs in Galactic GCs in the literature with properties that have allowed individual mass determination.

79 ASTRONOMY AND ASTROPHYSICS↗

Ruprecht 147: A Paradigm of Dissolving Star Cluster

We employed recent Gaia/DR2 data to investigate the dynamical status of the nearby (300 pc), old (2.5 Gyr) open cluster Ruprecht 147. We found prominent leading and trailing tails of stars along the cluster orbit, which demonstrates that Ruprecht 147 is losing stars at fast pace. Star counts indicate the cluster has a core radius of 33.3 arcmin and a tidal radius of 137.5 arcmin. The cluster also possesses an extended corona, which cannot be reproduced by a simple King model. We computed the present-day cluster mass using its luminosity and mass function, and derived an estimate of 234 ± 52 M {sub ⊙}. We also estimated the cluster original mass using available recipes extracted from N-body simulations obtaining a mass at birth of 50,000 ± 6500 M {sub ⊙}. Therefore dynamical mass loss, mostly caused by tidal interaction with the Milky Way, reduced the cluster mass by about 99%. We then conclude that Ruprecht 147 is rapidly dissolving into the general Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗

A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field

We have conducted a study of star formation in the outer Galaxy from 65°< l < 265° in the region observed by the GLIMPSE360 program. This Spitzer warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5 μm. We combine the IRAC, Wide-field Infrared Survey Explorer (WISE), and Two Micron All Sky Survey catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 young stellar objects (YSOs) across the field spanning >180° in Galactic longitude. Using the DBSCAN method on the combined catalog, we identify 618 clusters or aggregations of YSOs having five or more members. We identify 10,476 class I, 29,604 class II, and 7325 anemic class II/class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found that 100 of the clusters identified have previously measured distances in the WISE H ii survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with <3σ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters; dividing the clusters by galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M {sub ⊙} for R {sub Gal} < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M {sub ⊙} for R {sub Gal} > 11.5 kpc. The slope of the combined IMF was found to be Γ = 1.92 ± 0.42 above 3 M {sub ⊙}. These values are consistent with each other within the uncertainties and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗