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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Ground-state energy and relativistic corrections for positronium hydride

Variational calculations of the ground state of positronium hydride (HPs) are reported, including various expectation values, electron-positron annihilation rates, and leading relativistic corrections to the total and dissociation energies. The calculations have been performed using a basis set of 4000 thoroughly optimized explicitly correlated Gaussian basis functions. The relative accuracy of the variational energy upper bound is estimated to be of the order of 2x10{sup -10}, which is a significant improvement over previous nonrelativistic results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

A stochastic method for stand-alone photovoltaic system sizing

Photovoltaic systems utilize solar energy to generate electrical energy to meet load demands. Optimal sizing of these systems includes the characterization of solar radiation. Solar radiation at the Earth's surface has random characteristics and has been the focus of various academic studies. The objective of this study was to stochastically analyze parameters involved in the sizing of photovoltaic generators and develop a methodology for sizing of stand-alone photovoltaic systems. Energy storage for isolated systems and solar radiation were analyzed stochastically due to their random behavior. For the development of the methodology proposed stochastic analysis were studied including the Markov chain and beta probability density function. The obtained results were compared with those for sizing of stand-alone using from the Sandia method (deterministic), in which the stochastic model presented more reliable values. Both models present advantages and disadvantages; however, the stochastic one is more complex and provides more reliable and realistic results. (author)

14 SOLAR ENERGY↗
DOE OSTI2012

Linearized Functional Minimization for Inverse Modeling

Heterogeneous aquifers typically consist of multiple lithofacies, whose spatial arrangement significantly affects flow and transport. The estimation of these lithofacies is complicated by the scarcity of data and by the lack of a clear correlation between identifiable geologic indicators and attributes. We introduce a new inverse-modeling approach to estimate both the spatial extent of hydrofacies and their properties from sparse measurements of hydraulic conductivity and hydraulic head. Our approach is to minimize a functional defined on the vectors of values of hydraulic conductivity and hydraulic head fields defined on regular grids at a user-determined resolution. This functional is constructed to (i) enforce the relationship between conductivity and heads provided by the groundwater flow equation, (ii) penalize deviations of the reconstructed fields from measurements where they are available, and (iii) penalize reconstructed fields that are not piece-wise smooth. We develop an iterative solver for this functional that exploits a local linearization of the mapping from conductivity to head. This approach provides a computationally efficient algorithm that rapidly converges to a solution. A series of numerical experiments demonstrates the robustness of our approach.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2011

Computational technique for plasma parameters determination using Langmuir probe data

In the present work, we consider a new numerical method for processing the experimental information on the electron energy distribution function obtained with a Langmuir probe in a low-pressure plasma. This method offers the possibility to establish the temperature and concentration of the electrons for different forms of the distribution function. Some specific difficulties of the previous methods used to do such estimations are surpassed using the method proposed in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

GEOMETRIC AND DYNAMICAL MODELS OF REVERBERATION MAPPING DATA

We present a general method to analyze reverberation (or echo) mapping data that simultaneously provides estimates for the black hole mass and for the geometry and dynamics of the broad-line region (BLR) in active galactic nuclei (AGNs). While previous methods yield a typical scale size of the BLR or a reconstruction of the transfer function, our method directly infers the spatial and velocity distribution of the BLR from the data, from which a transfer function can be easily derived. Previous echo mapping analysis requires an independent estimate of a scaling factor known as the virial coefficient to infer the mass of the black hole, but this is not needed in our more direct approach. We use the formalism of Bayesian probability theory and implement a Markov Chain Monte Carlo algorithm to obtain estimates and uncertainties for the parameters of our BLR models. Fitting of models to the data requires knowledge of the continuum flux at all times, not just the measured times. We use Gaussian Processes to interpolate and extrapolate the continuum light curve data in a fully consistent probabilistic manner, taking the associated errors into account. We illustrate our method using simple models of BLR geometry and dynamics and show that we can recover the parameter values of our test systems with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping observing campaign. With a geometry model we can recover the mean radius of the BLR to within {approx}0.1 dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within {approx}0.05 dex random uncertainty, for simulated data with a line profile average signal-to-noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗
DOE OSTI2022

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗
DOE OSTI2016

On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The stellar mass assembly of low-redshift, massive, central galaxies in SDSS and the TNG300 simulation

ABSTRACT The stellar mass assembly of galaxies can be affected by both secular and environmental processes. In this study, for the first time, we investigate the stellar mass assembly of $\sim 90\, 000$ low-redshift, central galaxies selected from SDSS group catalogues ($M_{\rm Stellar}\gtrsim 10^{9.5}\, \mathrm{M}_{\odot }$, $M_{\rm Halo}\gtrsim 10^{12}\, \mathrm{M}_{\odot }$) as a function of both stellar mass and halo mass. We use estimates of the times at which 10, 50, and 90 per cent of the stellar mass were assembled from photometric spectral energy distribution fitting, allowing a more complete investigation than single stellar ages alone. We consider trends in both stellar mass and halo mass simultaneously, finding dependences of all assembly times on both. We find that galaxies with higher stellar masses (at constant halo mass) have on average older lookback times, similar to previous studies of galaxy assembly. We also find that galaxies at higher halo mass (at constant stellar mass) have younger lookback times, possibly due to a larger reservoir of gas for star formation. An exception to this is a subsample with high stellar-to-halo mass ratios, which are likely massive, field spirals. We compare these observed trends to those predicted by the TNG300 simulation, finding good agreement overall as a function of either stellar mass or halo mass. However, some differences in the assembly times (of up to ∼3 Gyr) appear when considering both stellar mass and halo mass simultaneously, noticeably at intermediate stellar masses (MStellar ∼ 1011 M⊙). These discrepancies are possibly linked to the quenched fraction of galaxies and the kinetic mode active galactic nucleus feedback implemented in TNG300.

Jackson, Thomas M.↗
DOE OSTI2022

Sensitivity Analysis of Fluid–Fluid Interfacial Area, Phase Saturation and Phase Connectivity on Relative Permeability Estimation Using Machine Learning Algorithms

Recent studies have shown that relative permeability can be modeled as a state function which is independent of flow direction and dependent upon phase saturation (S), phase connectivity (X), and fluid–fluid interfacial area (A). This study evaluates the impact of each of the three state parameters (S, X, and A) in the estimation of relative permeability. The relative importance of the three state parameters in four separate quadrants of S-X-A space was evaluated using a machine learning algorithm (out-of-bag predictor importance method). The results show that relative permeability is sensitive to all the three parameters, S, X, and A, with varying magnitudes in each of the four quadrants at a constant value of wettability. We observe that the wetting-phase relative permeability is most sensitive to saturation, while the non-wetting phase is most sensitive to phase connectivity. Although the least important, fluid–fluid interfacial area is still important to make the relative permeability a more exact state function.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗
DOE OSTI2020

Future projections and uncertainty assessment of precipitation extremes in the Korean peninsula from the CMIP5 ensemble

Abstract Projections of changes in extreme climate are sometimes predicted by multimodel ensemble methods that combine forecasts from individual simulation models using weighted averaging. One method to assign weight to each model is the Bayesian model averaging (BMA) in which posterior probability is used. For the cases of extreme climate, the generalized extreme value distribution (GEVD) is typically used. We applied the approach of GEV‐embedded BMA to a series of 35 years of the annual maximum daily precipitation data (both historical data and data gathered from simulation experiments for future periods) over the Korean peninsula as simulated by the models in the Coupled Model Intercomparison Project Phase 5 (CMIP5). Simulation data under two Representative Concentration Pathway (RCP) scenarios, namely RCP4.5 and RCP8.5, were used. Observed data and 17 CMIP5 models for 12 gird cells in Korea have been examined to predict future changes in precipitation extremes. A simple regional frequency analysis of pooling observations from three stations in each cell was employed to reduce the estimation variance and local fluctuations. A bias correction technique using the regression‐type transfer function was applied to these simulation data. Return levels spanning over 20 and 50 years, as well as the return periods relative to the reference years (1971–2005), were estimated for two future periods, namely Period 1 (2021–2050) and Period 2 (2066–2095). From these analyses, relative increase observed in the spatially averaged 20‐year (50‐year) return level was approximately 23% (16%) and 45% (36%) in the RCP4.5 and RCP8.5 experiments, respectively, by the end of the 21st century. We concluded that extreme rainfalls will likely occur two times and four times more frequently in the RCP4.5 and RCP8.5 scenarios, respectively, as compared to in the reference years by the end of the 21st century.

Lee, Youngsaeng↗
DOE OSTI2018

Global land model development: time to shift from a plant functional type to a plant functional trait approach (Final Report)

This project advanced global land models by shifting from the current plant functional type approach to one that better utilizes what is known about the importance and variability of plant traits, within a framework of simultaneously improving fundamental physiological relations that are at the core of model carbon cycling algorithms. A primary goal for Earth system modeling is to make accurate predictions of the future trajectory of the climate system, based on a mechanistic understanding of processes regulating fluxes of mass and energy among system components. Land plays an important role in modifying the Earth’s mass and energy balance, as a critical link in the global cycling of carbon, among others. Land surface models have developed to include mechanistic representations of vegetation physiology, carbon and nutrient dynamics in plants and soils, how they might respond to changing climate and chemistry, and how those changes might feedback to influence changes in atmospheric greenhouse gases themselves. Existing models represent the global distribution of vegetation types using the Plant Functional Type concept. Plant Functional Types are classes of plant species with similar evolutionary and life history with presumably similar responses to environmental conditions like CO 2 , water and nutrient availability. Fixed properties for each Plant Functional Type are specified through a collection of physiological parameters, or traits. These traits, mostly physiological in nature (e.g., leaf nitrogen and longevity) are used in model algorithms to estimate ecosystem properties and/or drive calculated process rates. In most models, 5 to 15 functional types represent terrestrial vegetation; in essence, they assume there are a total of only 5 to 15 different kinds of plants on the entire globe. This assumption of constant plant traits captured within the functional type concept has serious limitations, as a single set of traits does not reflect trait variation observed within and between species and communities. While this simplification was necessary decades past, substantial improvement is now possible. Rather than assigning a small number of constant parameter values to all grid cells in a model, we developed procedures that predict a frequency distribution of values for any given grid cell. Thus, by incorporating the mean and variance model fidelity has been improved. The trait-based approach has improved land modeling by (1) incorporating patterns and heterogeneity of traits into model parameterization, thus evolving away from a framework that considers large areas of vegetation to have near identical trait values and (2) allowing for improved treatment of physiological responses to environmental drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Electron Cloud Trapping in Recycler Combined Function Dipole Magnets

Electron cloud can lead to a fast instability in intense proton and positron beams in circular accelerators. In the Fermilab Recycler the electron cloud is confined within its combined function magnets. We show that the field of combined function magnets traps the electron cloud, present the results of analytical estimates of trapping, and compare them to numerical simulations of electron cloud formation. The electron cloud is located at the beam center and up to 1% of the particles can be trapped by the magnetic field. Since the process of electron cloud build-up is exponential, once trapped this amount of electrons significantly increases the density of the cloud on the next revolution. In a Recycler combined function dipole this multiturn accumulation allows the electron cloud reaching final intensities orders of magnitude greater than in a pure dipole. The multi-turn build-up can be stopped by injection of a clearing bunch of 1010 p at any position in the ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

A Method for Estimation and Suppression of Quasi-Periodic Noise

Electromagnetic interference, particularly from sources such as light fixtures and power supplies, negatively affects the function and performance of electronic systems. This noise, often quasi-periodic in nature, contains harmonic structure that can be used for noise estimation and suppression. This work proposes a novel method for tracking and filtering such noise signals that reduces the comb filter bandwidth required for noise estimation and suppression by tracking the instantaneous frequency of the fundamental component.

42 ENGINEERING↗
DOE OSTI2014

Vapor liquid solid-hydride vapor phase epitaxy (VLS-HVPE) growth of ultra-long defect-free GaAs nanowires: Ab initio simulations supporting center nucleation

High aspect ratio, rod-like and single crystal phase GaAs nanowires (NWs) were grown by gold catalyst-assisted hydride vapor phase epitaxy (HVPE). High resolution transmission electron microscopy and micro-Raman spectroscopy revealed polytypism-free zinc blende (ZB) NWs over lengths of several tens of micrometers for a mean diameter of 50 nm. Micro-photoluminescence studies of individual NWs showed linewidths smaller than those reported elsewhere which is consistent with the crystalline quality of the NWs. HVPE makes use of chloride growth precursors GaCl of which high decomposition frequency after adsorption onto the liquid droplet catalysts, favors a direct and rapid introduction of the Ga atoms from the vapor phase into the droplets. High influxes of Ga and As species then yield high axial growth rate of more than 100 μm/h. The diffusion of the Ga atoms in the liquid droplet towards the interface between the liquid and the solid nanowire was investigated by using density functional theory calculations. The diffusion coefficient of Ga atoms was estimated to be 3 × 10{sup −9} m{sup 2}/s. The fast diffusion of Ga in the droplet favors nucleation at the liquid-solid line interface at the center of the NW. This is further evidence, provided by an alternative epitaxial method with respect to metal-organic vapor phase epitaxy and molecular beam epitaxy, of the current assumption which states that this type of nucleation should always lead to the formation of the ZB cubic phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

WE-AB-202-06: Correlating Lung CT HU with Transformation-Based and Xe-CT Derived Ventilation

Purpose: Regional lung ventilation is useful to reduce radiation-induced function damage during lung cancer radiation therapy. Recently a new direct HU (Hounsfield unit)-based method was proposed to estimate the ventilation potential without image registration. The purpose of this study is to examine if there is a functional dependence between HU values and transformation-based or Xe-CT derived ventilation. Methods: 4DCT images acquired from 13 patients prior to radiation therapy and 4 mechanically ventilated sheep subjects which also have associated Xe-CT images were used for this analysis. Transformation-based ventilation was computed using Jacobian determinant of the transformation field between peak-exhale and peak-inhale 4DCT images. Both transformation and Xe-CT derived ventilation was computed for each HU bin. Color scatter plot and cumulative histogram were used to compare and validate the direct HU-based method. Results: There was little change of the center and shape of the HU histograms between free breathing CT and 4DCT average, with or without smoothing, and between the repeated 4DCT scans. HU of −750 and −630 were found to have the greatest transformation-based ventilation for human and sheep subjects, respectively. Maximum Xe-CT derived ventilation was found to locate at HU of −600 in sheep subjects. The curve between Xe-CT ventilation and HU was noisy for tissue above HU −400, possibly due to less intensity change of Xe gas during wash-out and wash-in phases. Conclusion: Both transformation-based and Xe-CT ventilation demonstrated that lung tissues with HU values in the range of (-750, −600) HU have the maximum ventilation potential. The correlation between HU and ventilation suggests that HU might be used to help guide the ventilation calculation and make it more robust to noise and image registration errors. Research support from NIH grants CA166703 and CA166119 and a gift from Roger Koch.

60 APPLIED LIFE SCIENCES↗