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 325 records · Page 18

Model Evaluation and Hindcasting: An Experiment with an Integrated Assessment Model

Integrated assessment models have been extensively used for analyzing long term energy and greenhouse emissions trajectories and have influenced key policies on this subject. Though admittedly these models are focused on the long term trajectories, how well these models are able to capture historical dynamics is an open question. In a first experiment of its kind, we present a framework for evaluation of such integrated assessment models. We use Global Change Assessment Model for this zero order experiment, and focus on the building sector results for USA. We calibrate the model for 1990 and run it forward up to 2095 in five year time steps. This gives us results for 1995, 2000, 2005 and 2010 which we compare to observed historical data at both fuel level and service level. We focus on bringing out the key insights for the wider process of model evaluation through our experiment with GCAM. We begin with highlighting that creation of an evaluation dataset and identification of key evaluation metric is the foremost challenge in the evaluation process. Our analysis highlights that estimation of functional form of the relationship between energy service demand, which is an unobserved variable, and its drivers is a significant challenge in the absence of adequate historical data for both the dependent and driver variables. Historical data availability for key metrics is a serious limiting factor in the process of evaluation. Interestingly, service level data against which such models need to be evaluated are itself a result of models. Thus for energy services, the best we can do is compare our model results with other model results rather than observed and measured data. We show that long term models, by the nature of their construction, will most likely underestimate the rapid growth in some services observed in a short time span. Also, we learn that modeling saturated energy services like space heating is easier than modeling unsaturated services like space cooling and understanding that how far a service is from its saturation level is a key question which we probably don’t have an answer to. Finally and most importantly, even if long term models partially miss the short term dynamics, the long term insights provides by these models is fairly robust. We conclude by highlighting that our work is the first step in the much wider process of integrated assessment model evaluation and will hence have its own limitations. Future evaluation research work should build upon this zero order experiment for improving our modeling of human and coupled earth systems.

building energy↗

Mechanism-based constitutive modeling of ZEK100 magnesium alloy with crystal plasticity and in-situ HEXRD experiment

The constitutive behavior of a hexagonal close-packed (HCP) polycrystalline ZEK100 magnesium alloy was investigated here using combined high energy X-ray diffraction (HEXRD) from a synchrotron source and crystal plasticity modeling approach. The in-situ tensile test data coupled with the HEXRD enabled the tracking of the lattice strain evolution during deformation. The microscopic behavior represented by lattice strain and the macroscopic behavior represented by stress-strain curves were then used together as objective function to estimate the critical resolved shear stress (CRSS) and hardening parameters of available slip and deformation twin systems in the ZEK100 alloy. An enhanced predominant twinning reorientation (ePTR) scheme was proposed in the current work, and the ePTR parameters were determined for the first time by the use of basal plane peak intensity along loading direction measured from HEXRD. Two crystal plasticity models, the computationally efficient elastic-plastic self-consistent (EPSC) and crystal plasticity finite element (CPFE) models, were developed incorporating the deformation twinning for the HCP-structured metals. The determined constitutive parameters were further validated by comparing the predicted deformation texture with the measured one. The work provides a useful and computationally-efficient modeling scheme to understand the slip/twin induced deformation behaviors of the ZEK100 alloy in micro- and macro-scales.

36 MATERIALS SCIENCE↗

Experimental and Computational Evaluation of Lipidomic In-Source Fragmentation as a Result of Postionization with Matrix-Assisted Laser Desorption/Ionization

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can provide spatially resolved molecular information about a sample. Recently, a postionization approach (MALDI-2) has been commercially integrated with MALDI-MSI, allowing for bettered sensitivity and consequent improved spatial resolution. While advantages of MALDI-2 have previously been established, we demonstrate here statistically increased in-source fragmentation (ISF) results from postionization with a commercial instrument. Via lipid standard analyses, known MALDI ISF pathways (e.g., loss of trimethylamine) were statistically increased in MALDI-2 compared to MALDI-1 (65–172% increase in fragmentation). Gas phase molecular modeling with density functional theory estimated that the most-weighted virtual orbitals to excite within lipids involve ester and phosphate bonds. Protonated lipid excitation energies are furthermore red-shifted compared to those of other adduct types [e.g., 254 nm for protonated PC(16:0/18:1)] and approach the MALDI-2 laser energy (266 nm). Analysis of rat brain homogenate detected statistically more positive-ion mode peaks with MALDI-2 (1090) than that with MALDI-1 (719), where Kernel density estimations showed that the majority of this enhancement occurs with low m/z ions (i.e., m/z 75–500). Taken together with the lipid standard data, these observations may indicate ISF due to postionization. Finally, while artifact contributions from matrix blanks were also noted, both experimental and computational data sets suggest that the overall extent of ISF is statistically increased in MALDI-2 compared to MALDI-1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatiotemporal and Statistical Mapping of Transition Metal Equilibria in Alkaline Media

Transition metal dissolution and redeposition (D/R) kinetics in alkaline media play a critical role in various chemical and electrochemical processes. Competitive reaction kinetics between different transition metals can modulate individual metal behavior in these processes. To date, these phenomena have remained largely unmeasured, and even when captured, they are difficult to statistically characterize due to their dynamic nature, simultaneous occurrence, and spatially heterogeneous nature. Here, in this study, we develop a statistical analysis framework based on in situ and operando X-ray fluorescence microscopy (XFM) to investigate the relative D/R kinetics of multiple transition metals in alkaline media. By employing statistical analysis, we quantify the spatial distribution of D/R species and assess the rate at which the system reaches equilibrium under varying reaction conditions. We show that pH does not simply change the rate of dissolution and redeposition, but reorganizes the cross-element kinetic correlations among Ni, Fe, and Mn and accelerates the spatial equilibration of D/R events, as quantified through correlation analysis, reaction-rate estimation, probability function distributions, and texture-based monitoring statistics. Additionally, we demonstrate how modifying the solvent environment can influence D/R kinetics, providing a pathway for tuning materials synthesis and process optimization. Our study offers valuable insights into the complex interplay between different transition metals and provides a reliable statistical framework for spatial analysis of diverse imaging data sets, enabling deeper extraction of latent information across multiple modalities.

36 MATERIALS SCIENCE↗

Stochastic simulation of electron avalanches on supercomputers

In the paper, we present a three-dimensional parallel Monte Carlo algorithm named ELSHOW which is developed for simulation of electron avalanches in gases. Parallel implementation of the ELSHOW was made on supercomputers with different architectures (massive parallel and hybrid ones). Using the ELSHOW, calculations of such integral characteristics as the number of particles in an avalanche, the coefficient of impact ionization, the drift velocity, and the others were made. Also, special precise computations were made to select an appropriate size of the time step using the technique of dependent statistical tests. Particularly, the algorithm consists of special methods of distribution modeling, a lexicographic implementation scheme for “branching” of trajectories, justified estimation of functionals. A comparison of the obtained results for nitrogen with previously published theoretical and experimental data was made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Importance sampling large deviations in nonequilibrium steady states. I

Large deviation functions hold information on the stability and response of systems driven into nonequilibrium steady states and in such a way are similar to free energies for systems at equilibrium. As with equilibrium free energies, evaluating large deviation functions numerically for all but the simplest systems is difficult because by construction they depend on exponentially rare events. In this first paper of a series, we assess different trajectory-based sampling methods capable of computing large deviation functions of time integrated observables within nonequilibrium steady states. We highlight some convergence criteria and best practices using a number of different models, including a biased Brownian walker, a driven lattice gas, and a model of self-assembly. We show how two popular methods for sampling trajectory ensembles, transition path sampling and diffusion Monte Carlo, suffer from exponentially diverging correlations in trajectory space as a function of the bias parameter when estimating large deviation functions. Improving the efficiencies of these algorithms requires introducing guiding functions for the trajectories.

74 ATOMIC AND MOLECULAR PHYSICS↗

The small numbers of large Kuiper Belt objects

We explore the brightness distribution of the largest and brightest (m(R) < 22) Kuiper Belt Objects (KBOs). We construct a luminosity function of the dynamically excited or hot Kuiper Belt (orbits with inclinations >5°) from the very brightest to m(R) = 23. We find for m(R) ≲ 23, a single slope appears to describe the luminosity function. We estimate that ∼12 KBOs brighter than m(R) ∼ 19.5 are present in the Kuiper Belt today. With nine bodies already discovered this suggests that the inventory of bright KBOs is nearly complete.

79 ASTRONOMY AND ASTROPHYSICS↗

Revealing asymmetries in the HD 181327 debris disk: A recent massive collision or interstellar medium warping

New multi-roll coronagraphic images of the HD 181327 debris disk obtained using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope reveal the debris ring in its entirety at high signal-to-noise ratio and unprecedented spatial resolution. We present and apply a new multi-roll image processing routine to identify and further remove quasi-static point-spread function-subtraction residuals and quantify systematic uncertainties. We also use a new iterative image deprojection technique to constrain the true disk geometry and aggressively remove any surface brightness asymmetries that can be explained without invoking dust density enhancements/deficits. The measured empirical scattering phase function for the disk is more forward scattering than previously thought and is not well-fit by a Henyey-Greenstein function. The empirical scattering phase function varies with stellocentric distance, consistent with the expected radiation pressured-induced size segregation exterior to the belt. Within the belt, the empirical scattering phase function contradicts unperturbed debris ring models, suggesting the presence of an unseen planet. The radial profile of the flux density is degenerate with a radially varying scattering phase function; therefore estimates of the ring's true width and edge slope may be highly uncertain. We detect large scale asymmetries in the disk, consistent with either the recent catastrophic disruption of a body with mass >1% the mass of Pluto, or disk warping due to strong interactions with the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

DESI mock challenge: constructing DESI galaxy catalogues based on FastPM simulations

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. $k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$ or $s\gt 10\, \mathrm{Mpc}\, h^{-1}$. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10 per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Load Profile Inpainting for Missing Load Data Restoration and Baseline Estimation

This paper introduces a Generative Adversarial Nets (GAN) based, Load Profile Inpainting Network (Load-PIN) for restoring missing load data segments and estimating the baseline for a demand response event. The inputs are time series load data before and after the inpainting period together with explanatory variables (e.g., weather data). Here, we propose a Generator structure consisting of a coarse network and a fine-tuning network. The coarse network provides an initial estimation of the data segment in the inpainting period. The fine-tuning network consists of self-attention blocks and gated convolution layers for adjusting the initial estimations. Loss functions are specially designed for the fine-tuning and the discriminator networks to enhance both the point-to-point accuracy and realisticness of the results. We test the Load-PIN on three real-world data sets for two applications: patching missing data and deriving baselines of conservation voltage reduction (CVR) events. We benchmark the performance of Load-PIN with five existing deep-learning methods. Our simulation results show that, compared with the state-of-the-art methods, Load-PIN can handle varying-length missing data events and achieve 15-30% accuracy improvement.

14 SOLAR ENERGY↗

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Next, using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing.

58 GEOSCIENCES↗

SU-G-JeP2-02: A Unifying Multi-Atlas Approach to Electron Density Mapping Using Multi-Parametric MRI for Radiation Treatment Planning

Purpose: MRI has a number of advantages over CT as a primary modality for radiation treatment planning (RTP). However, one key bottleneck problem still remains, which is the lack of electron density information in MRI. In the work, a reliable method to map electron density is developed by leveraging the differential contrast of multi-parametric MRI. Methods: We propose a probabilistic Bayesian approach for electron density mapping based on T1 and T2-weighted MRI, using multiple patients as atlases. For each voxel, we compute two conditional probabilities: (1) electron density given its image intensity on T1 and T2-weighted MR images, and (2) electron density given its geometric location in a reference anatomy. The two sources of information (image intensity and spatial location) are combined into a unifying posterior probability density function using the Bayesian formalism. The mean value of the posterior probability density function provides the estimated electron density. Results: We evaluated the method on 10 head and neck patients and performed leave-one-out cross validation (9 patients as atlases and remaining 1 as test). The proposed method significantly reduced the errors in electron density estimation, with a mean absolute HU error of 138, compared with 193 for the T1-weighted intensity approach and 261 without density correction. For bone detection (HU>200), the proposed method had an accuracy of 84% and a sensitivity of 73% at specificity of 90% (AUC = 87%). In comparison, the AUC for bone detection is 73% and 50% using the intensity approach and without density correction, respectively. Conclusion: The proposed unifying method provides accurate electron density estimation and bone detection based on multi-parametric MRI of the head with highly heterogeneous anatomy. This could allow for accurate dose calculation and reference image generation for patient setup in MRI-based radiation treatment planning.

60 APPLIED LIFE SCIENCES↗

k -d Darts: Sampling by k -dimensional flat searches

In this work, we introduced k-d darts as a particular type of higherdimensional sampling. We described a k-d dart framework for hyperplanes of general dimension k, and then demonstrated efficiency and accuracy over three applications using k = 1, and accuracy for one application using k = 1. In particular, darts produce accurate estimates of the volume of an object regardless of the dimension, orientation, and aspect ratio of the object. Axis-aligned darts are universally preferable to unaligned ones for sampling square domains, and we expect this to extend to hyperrectangles, such as bounding boxes. Darts also produce accurate mean estimates for function integration.

97 MATHEMATICS AND COMPUTING↗

Feedback-Based, System-Level Properties of Vertebrate-Microbial Interactions

Improved characterization of infectious disease dynamics is required. To that end, three-dimensional (3D) data analysis of feedback-like processes may be considered. To detect infectious disease data patterns, a systems biology (SB) and evolutionary biology (EB) approach was evaluated, which utilizes leukocyte data structures designed to diminish data variability and enhance discrimination. Using data collected from one avian and two mammalian (human and bovine) species infected with viral, parasite, or bacterial agents (both sensitive and resistant to antimicrobials), four data structures were explored: (i) counts or percentages of a single leukocyte type, such as lymphocytes, neutrophils, or macrophages (the classic approach), and three levels of the SB/ EB approach, which assessed (ii) 2D, (iii) 3D, and (iv) multi-dimensional (rotating 3D) host-microbial interactions. In all studies, no classic data structure discriminated disease-positive (D+, or observations in which a microbe was isolated) from disease-negative (D–, or microbial-negative) groups: D+ and D– data distributions overlapped. In contrast, multi-dimensional analysis of indicators designed to possess desirable features, such as a single line of observations, displayed a continuous, circular data structure, whose abrupt inflections facilitated partitioning into subsets statistically significantly different from one another. In all studies, the 3D, SB/EB approach distinguished three (steady, positive, and negative) feedback phases, in which D– data characterized the steady state phase, and D+ data were found in the positive and negative phases. In humans, spatial patterns revealed false-negative observations and three malaria-positive data classes. In both humans and bovines, methicillin-resistant Staphylococcus aureus (MRSA) infections were discriminated from non-MRSA infections. More information can be extracted, from the same data, provided that data are structured, their 3D relationships are considered, and well-conserved (feedback-like) functions are estimated. Patterns emerging from such structures may distinguish well-conserved from recently developed host-microbial interactions. Applications include diagnosis, error detection, and modeling.

59 BASIC BIOLOGICAL SCIENCES↗

An integrative computational approach for prioritization of genomic variants

An essential step in the discovery of molecular mechanisms contributing to disease phenotypes and efficient experimental planning is the development of weighted hypotheses that estimate the functional effects of sequence variants discovered by high-throughput genomics. With the increasing specialization of the bioinformatics resources, creating analytical workflows that seamlessly integrate data and bioinformatics tools developed by multiple groups becomes inevitable. Here we present a case study of a use of the distributed analytical environment integrating four complementary specialized resources, namely the Lynx platform, VISTA RViewer, the Developmental Brain Disorders Database (DBDB), and the RaptorX server, for the identification of high-confidence candidate genes contributing to pathogenesis of spina bifida. The analysis resulted in prediction and validation of deleterious mutations in the SLC19A placental transporter in mothers of the affected children that causes narrowing of the outlet channel and therefore leads to the reduced folate permeation rate. The described approach also enabled correct identification of several genes, previously shown to contribute to pathogenesis of spina bifida, and suggestion of additional genes for experimental validations. The study demonstrates that the seamless integration of bioinformatics resources enables fast and efficient prioritization and characterization of genomic factors and molecular networks contributing to the phenotypes of interest.

59 BASIC BIOLOGICAL SCIENCES↗

1999 Puget Sound Household Travel Survey

The 1999 Puget Sound Household Travel Survey was conducted between July and November. NuStats Research and Consulting conducted the survey on behalf of the Puget Sound Regional Council. The purpose of the study was to provide data to continue developing and refining the Regional Travel Demand Forecasting Model, as well as to provide a better understanding of travel behavior in the Puget Sound region. The study area consists of King, Kitsap, Pierce, and Snohomish counties. The resultant dataset will be used to fulfill the model's functions of estimating trip generation and distribution, mode choice, and assignments. The study had household members 16 years old or older keep track of travel for a 48-hour period. A total of 9,028 households were recruited to participate in the study. Of these, 6,000 households (66.5%) completed travel diaries, and the information was retrieved from all household members regardless of age. An “attitude” survey about transportation and land use issues was also mailed to household members 16 years old or older.

1Hz data↗

1998/99 Thurston County Household Travel Study

The survey was conducted under the auspices of the Thurston Regional Planning Council, and it was funded through a state grant awarded to Intercity Transit of Olympia, Washington. Data collection was from September 1998 through March 1999. The purpose of the study was to provide data for the continuing development and refinement of the Regional Travel Demand Forecasting Model, as well as to provide a better understanding of travel behavior in the southern Puget Sound region of Washington. The resultant data set will be used to fulfill the model's functions of estimating trip generation and distribution, mode choice, and assignments. Participating households were assigned specific “travel days” to record their travel over a 48-hour period. A total of 2,465 households were recruited to participate in the study. Of these, 1,537 households completed travel diaries, and the information was retrieved from 3,653 household members regardless of age. Households member made 25,278 total trips during their 48-hour diary period.

1Hz data↗

SHASTA-HELP

This tool contains functionality for estimating the storage potential of pure and blended natural gas-hydrogen mixtures in various subsurface formations and is intended to be used for pre-characterization or site screening purposes. Detailed storage analysis should be conducted by full reservoir simulation models.

Hydrogen↗