SEARCH · Engineering Papers

Results for “Parameter estimation using data”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Convective Parameters Derived from Radiosonde Data (SONDEPARAM) Value-Added Product Report

Radiosondes provide fundamental observations of the vertical profile of atmospheric state (pressure, temperature, humidity, and winds), with important implications for subsequent studies on environmental controls on cloud conditions. Within convective cloud environments, there is an increasing demand for additional value-added products (VAPs) to facilitate the use of U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility radiosonde data sets. Such VAPs should provide quick and reliable estimates for several standard radiosonde parameters or quantities of interest using common assumptions, as well as open, flexible code for visualization and user interaction. The Convective Parameters Derived from Radiosonde Data (SONDEPARAM) VAP will apply several robust algorithms used in Wang et al. (2020) for the calculation of useful radiosonde convective cloud parameters, including the convective available potential energy (CAPE), convective inhibition (CIN), and other convective parameters, for several different assumptions regarding the initial parcel characteristics (i.e., surface-based, most unstable, mixed layer). These ARM VAP codes are developed in open, flexible Python formats, with the intention that these parameters/calculations will be incorporated into traditional ARM quick-look radiosonde plotting, yet associated with user-available codes for ease in user reproduction and assumption modification.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

THE BIVARIATE SIZE-LUMINOSITY RELATIONS FOR LYMAN BREAK GALAXIES AT z {approx} 4-5

We study the bivariate size-luminosity distribution of Lyman break galaxies (LBGs) selected at redshifts around 4 and 5 in GOODS and the HUDF fields. We model the size-luminosity distribution as a combination of log-normal distribution (in size) and Schechter function (in luminosity), therefore it enables a more detailed study of the selection effects. We perform extensive simulations to quantify the dropout-selection completenesses and measurement biases and uncertainties in two-dimensional size and magnitude bins, and transform the theoretical size-luminosity distribution to the expected distribution for the observed data. Using maximum-likelihood estimator, we find that the Schechter function parameters for B {sub 435}-dropouts and are consistent with the values in the literature, but the size distributions are wider than expected from the angular momentum distribution of the underlying dark matter halos. The slope of the size-luminosity (RL) relation is similar to those found for local disk galaxies, but considerably shallower than local early-type galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Updated universal relations for tidal deformabilities of neutron stars from phenomenological equations of state

Equation of state (EOS) insensitive relations, so-called universal relations, between the neutron star (NS) compactness, its multipolar tidal deformability coefficients, and between the tidal parameters for binary systems are essential to break degeneracies in gravitational wave data analysis. Here, we validate and recalibrate these universal relations using a large set of almost 2 million phenomenological EOSs that are consistent with current observations. In doing so, we extend universal relations to a larger region of the EOS parameter space, most notably to softer EOSs and larger compactnesses. We show that waveform models that neglect higher-than-leading-order tidal deformations of the NSs accumulate as much as 3.5 radians of dephasing from 20Hz to merger. Further, we also perform a full Bayesian parameter estimation of the GW170817 data, and we compare the NS radius constraints produced using universal relations from the literature and the updated fits we propose here. We find that the new fits yield a NS radius that is smaller by about 500 meters. This difference is less than the statistical uncertainty on the radius at the signal-to-noise-ratio of GW170817, but it is significantly larger than the precision anticipated for next-generation detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Systematic validation of non-equilibrium thermochemical models using Bayesian inference

The validation process proposed by Babuška et al. [1] is applied to thermochemical models describing post-shock flow conditions. In this validation approach, experimental data is involved only in the calibration of the models, and the decision process is based on quantities of interest (QoIs) predicted on scenarios that are not necessarily amenable experimentally. Moreover, uncertainties present in the experimental data, as well as those resulting from an incomplete physical model description, are propagated to the QoIs. We investigate four commonly used thermochemical models: a one-temperature model (which assumes thermal equilibrium among all inner modes), and two-temperature models developed by Macheret et al. [2], Marrone and Treanor [3], and Park [4]. Up to 16 uncertain parameters are estimated using Bayesian updating based on the latest absolute volumetric radiance data collected at the Electric Arc Shock Tube (EAST) installed inside the NASA Ames Research Center. Following the solution of the inverse problems, the forward problems are solved in order to predict the radiative heat flux, QoI, and examine the validity of these models. Our results show that all four models are invalid, but for different reasons: the one-temperature model simply fails to reproduce the data while the two-temperature models exhibit unacceptably large uncertainties in the QoI predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗

DOE OSTI2016

Reconstruction of neuronal input through modeling single-neuron dynamics and computations

Mathematical models provide a mathematical description of neuron activity, which can better understand and quantify neural computations and corresponding biophysical mechanisms evoked by stimulus. In this paper, based on the output spike train evoked by the acupuncture mechanical stimulus, we present two different levels of models to describe the input-output system to achieve the reconstruction of neuronal input. The reconstruction process is divided into two steps: First, considering the neuronal spiking event as a Gamma stochastic process. The scale parameter and the shape parameter of Gamma process are, respectively, defined as two spiking characteristics, which are estimated by a state-space method. Then, leaky integrate-and-fire (LIF) model is used to mimic the response system and the estimated spiking characteristics are transformed into two temporal input parameters of LIF model, through two conversion formulas. We test this reconstruction method by three different groups of simulation data. All three groups of estimates reconstruct input parameters with fairly high accuracy. We then use this reconstruction method to estimate the non-measurable acupuncture input parameters. Results show that under three different frequencies of acupuncture stimulus conditions, estimated input parameters have an obvious difference. The higher the frequency of the acupuncture stimulus is, the higher the accuracy of reconstruction is.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion

This study examined the impacts of reservoir properties on CO2 migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO2 monitoring data such as saturation distribution. The injection reservoir was assumed to be located 1400-1500 m below the ground surface such that CO2 remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO2 injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended up to 8000 m from the injection well. The CO2 migration was simulated using the PNNL-developed simulator STOMP-CO2e (the water-salt-CO2 module). We adopted a nonlinear parameter estimation and optimization modeling software package, PEST, for automated reservoir parameter estimation. We explored the effects of data quality, data worth, and data redundancy on the detectability of reservoir parameters using CO2 saturation monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use of CO2 saturation monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy.

Carbon sequestration↗

DOE OSTI2013

Anisotropic parameter estimation using velocity variation with offset analysis

Seismic anisotropy is defined as velocity dependent upon angle or offset. Knowledge about anisotropy effect on seismic data is important in amplitude analysis, stacking process and time to depth conversion. Due to this anisotropic effect, reflector can not be flattened using single velocity based on hyperbolic moveout equation. Therefore, after normal moveout correction, there will still be residual moveout that relates to velocity information. This research aims to obtain anisotropic parameters, ε and δ, using two proposed methods. The first method is called velocity variation with offset (VVO) which is based on simplification of weak anisotropy equation. In VVO method, velocity at each offset is calculated and plotted to obtain vertical velocity and parameter δ. The second method is inversion method using linear approach where vertical velocity, δ, and ε is estimated simultaneously. Both methods are tested on synthetic models using ray-tracing forward modelling. Results show that δ value can be estimated appropriately using both methods. Meanwhile, inversion based method give better estimation for obtaining ε value. This study shows that estimation on anisotropic parameters rely on the accuracy of normal moveout velocity, residual moveout and offset to angle transformation.

58 GEOSCIENCES↗

DOE OSTI2024

Blueprints for Training Information Bottlenecks for Collider Analyses

Dimensionality reduction is a crucial aspect of data analysis in high energy physics, even if accompanied by information loss. Several methods, including histogram- and kernel-based analyses, are only computationally feasible for low-dimensional data. Furthermore, simulation models used in HEP can often only be validated for low-dimensional data. We provide several blueprints for using machine learning to create low-dimensional data representations (continuous event variables and discrete classification labels) for use in signal discovery and parameter estimation tasks. We also describe how to design the learned representation to facilitate a) searches with unknown model parameters and b) validation of simulation models in data control regions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

How deep is your soil? Quantifying and spatially analyzing understudied deep soil in the United States

Deep soil is largely understudied and important in understanding biogeochemical processes in soil. Here, understudied soil is defined as the difference between soil studied to a known depth and the estimated bedrock depth. To understand more about deep soil, the understudied soil in the US was quantified and spatially analyzed using soil survey data and model estimates of bedrock depth. An equation was derived to find understudied soil using the dataset parameters “max lower depth studied”, “depth to bedrock”, and “likelihood of bedrock in the top 200 cm”. The survey data and bedrock model revealed that soil has been studied to an average depth of 1-2 meters, and the average depth to bedrock is 20 meters. Soil data density in the soil surveys was greatest in the West Coast, Midwest, and areas historically managed for agricultural, while the non-contiguous US and interior West were underrepresented. The soil had been studied deeper than the estimated soil depth in 455 out of 56,889 observation points concentrated in Alaska, California, Texas, Florida, Puerto Rico, and the US Virgin Islands. To understand the diversity and any taxonomic bias of the global soil data available, soil order was compared to US-based National Resource Conservation Service percentages and it was found that Oxisols, Alfisols, Ultisols, Andisols, and Histosols were overrepresented while Gelisols, Aridisols, Vertisols, Entisols, and Spodosols are underrepresented. Soil depth is important in exploring the complexity of biogeochemical processes that take place in soil.

Bedrock↗

DOE OSTI2016

Parameter Estimator for Engineering Systems

This software model generates vehicle parameter sets for downstream modeling applications using real-world data. Value to users is provided through improvements to modeling activities as a result of decreases in time taken during model development and validation. Currently configured for transportation modeling activities, the software provides initial estimates for multiple fundamental vehicle parameters including mass, coefficient of drag, and tire rolling resistance using real-world drive and duty cycle information as inputs.

Source record↗

DOE OSTI2012

MODELING EXTRAGALACTIC FOREGROUNDS AND SECONDARIES FOR UNBIASED ESTIMATION OF COSMOLOGICAL PARAMETERS FROM PRIMARY COSMIC MICROWAVE BACKGROUND ANISOTROPY

Using the latest physical modeling and constrained by the most recent data, we develop a phenomenological parameterized model of the contributions to intensity and polarization maps at millimeter wavelengths from external galaxies and Sunyaev-Zeldovich effects. We find such modeling to be necessary for estimation of cosmological parameters from Planck data. For example, ignoring the clustering of the infrared background would result in a bias in n{sub s} of 7{sigma} in the context of an eight-parameter cosmological model. We show that the simultaneous marginalization over a full foreground model can eliminate such biases, while increasing the statistical uncertainty in cosmological parameters by less than 20%. The small increases in uncertainty can be significantly reduced with the inclusion of higher-resolution ground-based data. The multi-frequency analysis we employ involves modeling 46 total power spectra and marginalization over 17 foreground parameters. We show that we can also reduce the data to a best estimate of the cosmic microwave background power spectra, with just two principal components (with constrained amplitudes) describing residual foreground contamination.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

Perfusion Estimated With Rapid Dynamic Contrast-Enhanced Magnetic Resonance Imaging Correlates Inversely With Vascular Endothelial Growth Factor Expression and Pimonidazole Staining in Head-and-Neck Cancer: A Pilot Study

Purpose: To analyze, in a pilot study, rapidly acquired dynamic contrast-enhanced (DCE)-MRI data with a general two-compartment exchange tracer kinetic model and correlate parameters obtained with measurements of hypoxia and vascular endothelial growth factor (VEGF) expression in patients with squamous cell carcinoma of the head and neck. Methods and Materials: Eight patients were scanned before surgery. The DCE-MRI data were acquired with 1.5-s temporal resolution and analyzed using the two-compartment exchange tracer kinetic model to obtain estimates of parameters including perfusion and permeability surface area. Twelve to 16 h before surgery, patients received an intravenous injection of pimonidazole. Samples taken during surgery were used to determine the level of pimonidazole staining using immunohistochemistry and VEGF expression using quantitative real-time polymerase chain reaction. Correlations between the biological and imaging data were examined. Results: Of the seven tumors fully analyzed, those that were poorly perfused tended to have high levels of pimonidazole staining (r = -0.79, p = 0.03) and VEGF expression (r = -0.82, p = 0.02). Tumors with low permeability surface area also tended to have high levels of hypoxia (r = -0.75, p = 0.05). Hypoxic tumors also expressed higher levels of VEGF (r = 0.82, p = 0.02). Conclusions: Estimates of perfusion obtained with rapid DCE-MRI data in patients with head-and-neck cancer correlate inversely with pimonidazole staining and VEGF expression.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2016

Quantifying shallow subsurface water and heat dynamics using coupled hydrological-thermal-geophysical inversion

Improving our ability to estimate the parameters that control water and heat fluxes in the shallow subsurface is particularly important due to their strong control on recharge, evaporation and biogeochemical processes. The objectives of this study are to develop and test a new inversion scheme to simultaneously estimate subsurface hydrological, thermal and petrophysical parameters using hydrological, thermal and electrical resistivity tomography (ERT) data. The inversion scheme – which is based on a nonisothermal, multiphase hydrological model – provides the desired subsurface property estimates in high spatiotemporal resolution. A particularly novel aspect of the inversion scheme is the explicit incorporation of the dependence of the subsurface electrical resistivity on both moisture and temperature. The scheme was applied to synthetic case studies, as well as to real datasets that were autonomously collected at a biogeochemical field study site in Rifle, Colorado. At the Rifle site, the coupled hydrological-thermal-geophysical inversion approach well predicted the matric potential, temperature and apparent resistivity with the Nash–Sutcliffe efficiency criterion greater than 0.92. Synthetic studies found that neglecting the subsurface temperature variability, and its effect on the electrical resistivity in the hydrogeophysical inversion, may lead to an incorrect estimation of the hydrological parameters. The approach is expected to be especially useful for the increasing number of studies that are taking advantage of autonomously collected ERT and soil measurements to explore complex terrestrial system dynamics.

58 GEOSCIENCES↗

DOE OSTI2016

PEES (Parameter Estimator for Engineering Systems) [SWR-16-38]

Miller, Eric↗

DOE OSTI2014

EXONEST: Bayesian model selection applied to the detection and characterization of exoplanets via photometric variations

EXONEST is an algorithm dedicated to detecting and characterizing the photometric signatures of exoplanets, which include reflection and thermal emission, Doppler boosting, and ellipsoidal variations. Using Bayesian inference, we can test between competing models that describe the data as well as estimate model parameters. We demonstrate this approach by testing circular versus eccentric planetary orbital models, as well as testing for the presence or absence of four photometric effects. In addition to using Bayesian model selection, a unique aspect of EXONEST is the potential capability to distinguish between reflective and thermal contributions to the light curve. A case study is presented using Kepler data recorded from the transiting planet KOI-13b. By considering only the nontransiting portions of the light curve, we demonstrate that it is possible to estimate the photometrically relevant model parameters of KOI-13b. Furthermore, Bayesian model testing confirms that the orbit of KOI-13b has a detectable eccentricity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗

DOE OSTI2014

SU-E-I-95: Personalized Radiography Technical Parameters for Each Patient and Exam

Purpose: To determine exact electrical parameters (kV, mAs) a radiological technologist shall use taking account the exam and patient's structure, with guarantee of minimum dose and adequate quality image. Methods: A patient's absorbed dose equation was developed by means of Entrance Skin Dose (ESD), irradiated area and patient width for specific anatomy. ESD is calculated from a developed equation, where entrance surface air-KERMA and backscatter factor are included, with air-to-skin coefficient conversion. We developed specific Lambert-Beer attenuation equations derived from mass energy-absorption coefficients data for skin, fat, and muscle and bone as one tissue. Anatomy tissue thickness distribution at central X-ray location in anteroposterior incidence for hand and chest, was estimate by discounting constant skin and bone thickness from patient measured width, assuming the result as muscle and fat. A clinical research at a big hospital were executed when real parameters (kV, mAs, filtration, ripple) used by technologists were combined with the image quality and patient's data: anatomy width, height and weight. A correlation among the best images acquired and electrical parameters used were confronted with patient's data and dose estimation. The best combinations were used as gold standards. Results: For each anatomy, two equations were developed to calculate voltage (kV) and exposure (mAs) to reproduce and interpolate the gold standards. Patient is measured and data are input into equations, giving radiological technologists the right set of electrical parameters for that specific exam. Conclusion: This work indicates that radiological technologist can personalize the exact electrical parameters for each patient exam, instead of using standard values. It also guarantee that patients under or over-sized measures will receive the right dose for the best image. It will stop wrong empiric adjusts technologists do when examining a non-standard patient and reduce probability of radiography retaken because of over or under exposition.

60 APPLIED LIFE SCIENCES↗