Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 1,369 records · Page 76

iTOUGH2: A multiphysics simulation-optimization framework for analyzing subsurface systems

We report iTOUGH2 is a simulation-optimization framework for the TOUGH suite of nonisothermal multiphase flow models and related simulators of geophysical, geochemical, and geomechanical processes. After appropriate parameterization of subsurface structures and their properties, iTOUGH2 runs simulations for multiple parameter sets and analyzes the resulting output for parameter estimation through automatic model calibration, local and global sensitivity analyses, data-worth analyses, and uncertainty propagation analyses. Development of iTOUGH2 is driven by scientific challenges and user needs, with new capabilities continually added to both the forward simulator and the optimization framework. This paper provides a summary description of methods and features implemented in iTOUGH2, and discusses the usefulness and limitations of an integrated simulation-optimization workflow in support of the characterization and analysis of complex multiphysics subsurface systems.

58 GEOSCIENCES↗

A Nonempirical Relative Permeability Model for Hydrate-Bearing Sediments

There are currently two types of relative permeability models that are used to model gas production from hydrate-bearing sediments: fully empirical parameter-fitting models [such as the University of Tokyo model (Masuda et al. 1997) and the Brooks and Corey model (Brooks and Corey 1964)] and partially empirical models [such as the Kozeny and Carman model (Wyllie and Gardner 1958) and capillary-tube-based models that assume only a single phase]. This study proposes an analytical model to estimate relative permeability of gas and water in a hydrate-bearing porous medium without curve fitting or use of any empirical parameters. The model is derived by conserving the momentum balance with the steady-state form of the Navier-Stokes equation for gas/water flow in a hydrate-bearing porous medium. The model is validated against a number of laboratory studies and is shown to perform better than most empirical models over a full range of experimental data. The proposed model is an analytical function of rock properties (average pore size and shape, porosity, irreducible water saturation, and saturation of hydrate), fluid properties (gas/water saturations and viscosities), and the hydrate-growth pattern [pore filling (PF), wall coating (WC), and a combination of PF and WC]. Here, the benefits of the proposed model include sensitivity analysis of relevant physical parameters on relative permeability and estimation of rock parameters (such as porosity, pore size, and residual water saturation) using inverse modeling. The model can also be used to estimate two-phase permeability in a permeable medium without hydrates.

58 GEOSCIENCES↗

Correlations of the Elemental Compositions of Primary Coal Tar and Char

We report that simulations of large-scale coal combustors rely on accurate submodels to describe the chemical and physical changes in coal during reaction. Typically, researchers use simplified empirical submodels tuned to experimental data to reduce the computational complexity. When data are not readily available, researchers use simplifying assumptions, which can create inaccuracies and biases in a large simulation. One such simplifying assumption in coal research is how to describe the elemental composition of primary pyrolysis products. This paper explores several different empirical model forms to predict the dry, ash-free fractions of C, H, O, N, and S in both the char and the tar, using variables such as parent coal composition, reaction conditions (temperatures and particle residence times), and key coal structural parameters derived from NMR measurements. These model forms were correlated to existing data using a wide range of experimental data using a cross-validation procedure. Since coal structural values can be expensive to measure, several correlations from the literature were used to estimate these values based on information from the proximate and ultimate analyses of the parent coal, including a new correlation for the coal aromaticity. These model forms were tested against a set of measured elemental compositions of tar and char to find the best fit to use in the cross-validation process. Lastly, the best empirical models are presented that predict elemental composition of the coal char and tar after devolatilization.

01 COAL, LIGNITE, AND PEAT↗

Optical and UV spectroscopy of the black hole binary candidate LMC X-1

Both further optical spectroscopy of the binary star identified with LMC X-1, obtained between 1983 and 1985, and a series of IUE UV spectra taken during a 5 day interval in 1984 are presented. The optical data are used to refine the orbital period to 4.2288 days, and improved orbital parameters are derived. The velocity of the optical emission lines is antiphased with the absorption lines and has twice the velocity amplitude. These new results support the estimates of the masses in the system given earlier. The most probable component masses are approximately 20 solar masses for the primary and near 6 solar masses (for the x-ray star), suggesting the the latter may be a black hole. The UV spectra show very weak, low-velocity stellar-wind lines. It is suggested that much of the surrounding medium is highly ionized by the X-ray flux. The 'nonwind' UV spectral lines and the UV continuum temperature are consistent with the optical data, indicating a late O type star of M(bol) = -8.5. There is a weak modulation of absorption-line strengths with orbital phase, suggestive of a lack of axisymmetry in the X-irradiation of the primary star and indicative of a fairly low orbital inclination.

Hutchings, J. B.↗

Improved canopy reflectance modeling and scene inference through improved understanding of scene pattern

The Li-Strahler reflectance model, driven by LANDSAT Thematic Mapper (TM) data, provided regional estimates of tree size and density within 20 percent of sampled values in two bioclimatic zones in West Africa. This model exploits tree geometry in an inversion technique to predict average tree size and density from reflectance data using a few simple parameters measured in the field (spatial pattern, shape, and size distribution of trees) and in the imagery (spectral signatures of scene components). Trees are treated as simply shaped objects, and multispectral reflectance of a pixel is assumed to be related only to the proportions of tree crown, shadow, and understory in the pixel. These, in turn, are a direct function of the number and size of trees, the solar illumination angle, and the spectral signatures of crown, shadow and understory. Given the variance in reflectance from pixel to pixel within a homogeneous area of woodland, caused by the variation in the number and size of trees, the model can be inverted to give estimates of average tree size and density. Because the inversion is sensitive to correct determination of component signatures, predictions are not accurate for small areas.

Franklin, Janet↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

Variation of Ice Microphysical Properties With Temperature and Humidity at Tops of Convective Clouds

A better understanding of the many interacting processes governing the evaluation of ice in natural clouds is required to improve the representation of ice clouds in global circulation models. Recent studies suggest a dominant role of vapor growth processes in determining the temperature dependence of cloud top ice sizes and shapes. Using airborne cloud remote sensing along with reanalysis data, here we show that observed cloud top ice effective radii and estimated normalized growth rates at cloud top highly correlate with an approximately linear relationship, which is consistent with a conceptual model also presented. Furthermore, significant differences in crystal shape characteristics and scattering asymmetry parameters are found between sub- and super-saturated cloud tops over ocean, although not over land. These results provide valuable observational targets for studying ice formation and evolution processes using models, while also helping interpretation of satellite observations of ice microphysical properties at cloud tops.

ice↗

CONUS and Global Carbon Consumption Database for Wildland Fire

Fire plays a significant role on both national and global scales, profoundly impacting landscapes shaped by human activity as well as wildlands. Even though fire can be devastating, wildland fire is a natural and integral force on our landscapes. Fires can also serve to reduce fuels to mitigate wildfire risk and maintain healthy ecosystem functions. However, the smoke produced by fires, regardless of their size or purpose, can pose adverse effects on human health when inhaled downwind. Understanding the influence of smoke on air quality and human well-being necessitates the quantification of emissions that fires release into the atmosphere. In response to this need, we have developed a carbon consumption database for the Continental United States (CONUS) and are developing a global carbon consumption database, both of which are directly link to distinct fuels within various fire danger categories. Fuel Characteristic Classification System (FCCS) fuels are used to parameterize biomass, and consumption is broken down into five Fire Danger categories (Low, Moderate, High, Very High, Extreme), for both ‘new’ and ‘residual’ burning scenarios. Residual burned area is defined as burning in areas that have recently burned. We implemented the FCCS 30-meter CONUS fuelbed dataset during the 2019 Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) campaign to estimate daily carbon emissions. Our emissions estimates were rigorously compared against in-situ measurements of CO2, CO, and black carbon aerosols, revealing a robust agreement between the two datasets. The 300-meter global database builds upon the foundations of the Pettinari, M. Lucrecia (2015) Global Fuelbed database, a global fuel map with standardized FCCS biomass parameters. These products can serve as valuable tools when used in conjunction with burned area data to rapidly and accurately estimate carbon consumed and released into the atmosphere.

FCCS↗

Multistage variable probability forest volume inventory

An inventory scheme based on the use of computer processed LANDSAT MSS data was developed. Output from the inventory scheme provides an estimate of the standing net saw timber volume of a major timber species on a selected forested area of the Navajo Nation. Such estimates are based on the values of parameters currently used for scaled sawlog conversion to mill output. The multistage variable probability sampling appears capable of producing estimates which compare favorably with those produced using conventional techniques. In addition, the reduction in time, manpower, and overall costs lend it to numerous applications.

Anderson, J. E.↗

DREISS: Using State-Space Models to Infer the Dynamics of Gene Expression Driven by External and Internal Regulatory Networks

Gene expression is controlled by the combinatorial effects of regulatory factors from different biological subsystems such as general transcription factors (TFs), cellular growth factors and microRNAs. A subsystem’s gene expression may be controlled by its internal regulatory factors, exclusively, or by external subsystems, or by both. It is thus useful to distinguish the degree to which a subsystem is regulated internally or externally–e.g., how non-conserved, species-specific TFs affect the expression of conserved, cross-species genes during evolution. We developed a computational method (DREISS, dreiss.gerteinlab.org) for analyzing the Dynamics of gene expression driven by Regulatory networks, both External and Internal based on State Space models. Given a subsystem, the “state” and “control” in the model refer to its own (internal) and another subsystem’s (external) gene expression levels. The state at a given time is determined by the state and control at a previous time. Because typical time-series data do not have enough samples to fully estimate the model’s parameters, DREISS uses dimensionality reduction, and identifies canonical temporal expression trajectories (e.g., degradation, growth and oscillation) representing the regulatory effects emanating from various subsystems. To demonstrate capabilities of DREISS, we study the regulatory effects of evolutionarily conserved vs. divergent TFs across distant species. In particular, we applied DREISS to the time-series gene expression datasets of C. elegans and D. melanogaster during their embryonic development. We analyzed the expression dynamics of the conserved, orthologous genes (orthologs), seeing the degree to which these can be accounted for by orthologous (internal) versus species-specific (external) TFs. We found that between two species, the orthologs have matched, internally driven expression patterns but very different externally driven ones. This is particularly true for genes with evolutionarily ancient functions (e.g. the ribosomal proteins), in contrast to those with more recently evolved functions (e.g., cell-cell communication). This suggests that despite striking morphological differences, some fundamental embryonic-developmental processes are still controlled by ancient regulatory systems.

Biochemistry & Molecular Biology↗

On the estimation of boundary layer heights: a machine learning approach

Abstract. The planetary boundary layer height (zi) is a key parameter used in atmospheric models for estimating the exchange of heat, momentum, and moisture between the surface and the free troposphere. Near-surface atmospheric and subsurface properties (such as soil temperature, relative humidity, etc.) are known to have an impact on zi. Nevertheless, precise relationships between these surface properties and zi are less well known and not easily discernible from the multi-year dataset. Machine learning approaches, such as random forest (RF), which use a multi-regression framework, help to decipher some of the physical processes linking surface-based characteristics to zi. In this study, a 4-year dataset from 2016 to 2019 at the Southern Great Plains site is used to develop and test a machine learning framework for estimating zi. Parameters derived from Doppler lidars are used in combination with over 20 different surface meteorological measurements as inputs to a RF model. The model is trained using radiosonde-derived zi values spanning the period from 2016 through 2018 and then evaluated using data from 2019. Results from 2019 showed significantly better agreement with the radiosonde compared to estimates derived from a thresholding technique using Doppler lidars only. Noteworthy improvements in daytime zi estimates were observed using the RF model, with a 50 % improvement in mean absolute error and an R2 of greater than 85 % compared to the Tucker method zi. We also explore the effect of zi uncertainty on convective velocity scaling and present preliminary comparisons between the RF model and zi estimates derived from atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Design criteria and eigensequence plots for satellite computed tomography

The use of the degrees of freedom for signal is proposed as a design criteria for comparing different designs for satellite and other measuring systems. It is also proposed that certain eigensequence plots be examined at the design stage along with appropriate estimates of the parameter lambda playing the role of noise to signal ratio. The degrees of freedom for signal and the eigensequence plots may be determined using prior information in the spectral domain which is presently available along with a description of the system, and simulated data for estimating lambda. This work extends the 1972 work of Weinreb and Crosby.

Wahba, G.↗

Design criteria and eigensequence plots for satellite-computed tomography

The use of the degrees of freedom for signal is proposed as a design criteria for comparing different designs for satellite and other measuring systems. It is also proposed that certain eigensequence plots be examined at the design stage along with appropriate estimates of the parameter lambda playing the role of noise to signal ratio. The degrees of freedom for signal and the eigensequence plots may be determined using prior information in the spectral domain which is presently available along with a description of the system, and simulated data for estimating lambda. This work extends the 1972 work of Weinreb and Crosby.

Wahba, G.↗

Optical Navigation for Autonomous Approach of Small Unknown Bodies

State of the practice in navigation around small celestial bodies heavily relies on ground sup- port and human skill, in particular, for perception-based operations such as optical navigation and mapping. This leads to longer duration and more complex mission operations and sub- sequently higher cost. Furthermore, it imposes limitations for certain missions such as fast fly-bys or multi-agent operations. In this work, we present an autonomous navigation strat- egy suitable for approaching small unexplored bodies. During the approach, we estimate the body’s physical properties as well as the spacecraft’s relative trajectory and associated un- certainties. The autonomous navigation strategy, which is solely based on optical measure- ments, begins as soon as the body becomes resolved in the navigation camera and terminates at the start of proximity operations, when the spacecraft makes its first trajectory correction to stay in the vicinity of the body. Our strategy uses multiple image-processing algorithms: light-curve analysis for estimating the target body’s rotation rate, Shape-from-Silhouette for reconstructing the 3D shape and estimating its rotation pole, and feature tracking tailored to Small-Body images for estimating relative navigation parameters. We used the Mission Analysis, Operations, and Navigation Toolkit Environment (MONTE) developed by the Jet Propulsion Laboratory to evaluate the feasibility of this multi-phase navigation strategy using simulated images of an approach trajectory. We used the Rosetta mission data to generate photorealistic images to characterise the performance of this approach. This work is based on the assumptions that the spacecraft attitude is known, the body is a principal-axis rotator, a-priori estimates of ephemerides and scale are available, and the body is observed from a zero sun phase only during initial approach. Preliminary results show orbit determination performance that is on par with the human navigation from the Rosetta mission; albeit with a 1% bias in spacecraft-target radial distance estimate. The bias error is likely due to the robustness and accuracy of the visual tracking under dynamic lighting conditions and per- spective changes, which decrease accuracy.

Villa, Jacopo↗

Bayes estimation on parameters of the single-class classifier

Normal procedures used for designing a Bayes classifier to classify wheat as the major crop of interest require not only training samples of wheat but also those of nonwheat. Therefore, ground truth must be available for the class of interest plus all confusion classes. The single-class Bayes classifier classifies data into the class of interest or the class 'other' but requires training samples only from the class of interest. This paper will present a procedure for Bayes estimation on the mean vector, covariance matrix, and a priori probability of the single-class classifier using labeled samples from the class of interest and unlabeled samples drawn from the mixture density function.

Lin, G. C.↗

Uncertainty Analysis of Inertial Model Attitude Sensor Calibration and Application with a Recommended New Calibration Method

Statistical tools, previously developed for nonlinear least-squares estimation of multivariate sensor calibration parameters and the associated calibration uncertainty analysis, have been applied to single- and multiple-axis inertial model attitude sensors used in wind tunnel testing to measure angle of attack and roll angle. The analysis provides confidence and prediction intervals of calibrated sensor measurement uncertainty as functions of applied input pitch and roll angles. A comparative performance study of various experimental designs for inertial sensor calibration is presented along with corroborating experimental data. The importance of replicated calibrations over extended time periods has been emphasized; replication provides independent estimates of calibration precision and bias uncertainties, statistical tests for calibration or modeling bias uncertainty, and statistical tests for sensor parameter drift over time. A set of recommendations for a new standardized model attitude sensor calibration method and usage procedures is included. The statistical information provided by these procedures is necessary for the uncertainty analysis of aerospace test results now required by users of industrial wind tunnel test facilities.

Tripp, John S.↗

Analysis of Particle Image Velocimetry (PIV) Data for Acoustic Velocity Measurements

Acoustic velocity measurements were taken using Particle Image Velocimetry (PIV) in a Normal Incidence Tube configuration at various frequency, phase, and amplitude levels. This report presents the results of the PIV analysis and data reduction portions of the test and details the processing that was done. Estimates of lower measurement sensitivity levels were determined based on PIV image quality, correlation, and noise level parameters used in the test. Comparison of measurements with linear acoustic theory are presented. The onset of nonlinear, harmonic frequency acoustic levels were also studied for various decibel and frequency levels ranging from 90 to 132 dB and 500 to 3000 Hz, respectively.

Blackshire, James L.↗

Method and apparatus for autonomous, in-receiver prediction of GNSS ephemerides

Methods and apparatus for autonomous in-receiver prediction of orbit and clock states of Global Navigation Satellite Systems (GNSS) are described. Only the GNSS broadcast message is used, without need for periodic externally-communicated information. Earth orientation information is extracted from the GNSS broadcast ephemeris. With the accurate estimation of the Earth orientation parameters it is possible to propagate the best-fit GNSS orbits forward in time in an inertial reference frame. Using the estimated Earth orientation parameters, the predicted orbits are then transformed into Earth-Centered-Earth-Fixed (ECEF) coordinates to be used to assist the GNSS receiver in the acquisition of the signals. GNSS satellite clock states are also extracted from the broadcast ephemeris and a parameterized model of clock behavior is fit to that data. The estimated modeled clocks are then propagated forward in time to enable, together with the predicted orbits, quicker GNSS signal acquisition.

Bar-Sever, Yoaz E.↗