Engineering Papers⌕ Search

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.

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

At least 829 records · Page 46

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (T(sub g)) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing longwave radiation (CLR) is sensitive to upper level moisture (q(sub l)) over wet regions and (T(sub g)) over dry regions The clear sky window flux from 800 to 1200/cm (RadWn) is sensitive to low level moisture (q(sub t)) and T(sub g). Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub l) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in T(sub g), q(sub h) and q(sub l). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing Delta(T(sub g)) and broad layer Delta(q(sub l) from .500 hPa to surface and Delta(q(sub h)) from 200 to .300 hPa provides a good approximation to the full radiative transfer calculations. typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for Delta(T(sub g)) to about 20% for Delta(q(sub l)). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 min). We note that these errors, as well as a cold bias in the T(sub g). have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSM/I total precipitable water.

Wu, Man Li C.↗

Size_resolved_ice_nucleating_particle_INP_concentrations_from_the_INPOP_campaign

Ground-based immersion mode ice nucleation measurements were conducted on samples collected at the US DOE ARM AMF3 in Oliktok Point, AK from 1 Mar to 31 May 2017. Samples were collected using a Davis Rotating-drum Unit for Monitoring (DRUM; DRUMAir, LLC) cascading impactor in 4 stages (i.e., size bins). Instrument information on the model used, the DA400, can be found at the DRUMAir website (http://www.drumair.com/). Stages A, B, C, and D, had particle size cuts at 2.96, 1.21, 0.34, and 0.15 µm in particle diameter. Strips of sample substrate (PFA, perfluoroalkoxy; McMaster) were adhered to each disc in each stage and coated with petrolatum to enable particles to stick on impact. The discs move slowly over time, such that aerosol loading is streaked onto the PFA. Every 24 hours, a blank is created on the PFA to separate daily samples. The pump typically pulled 25 lpm at the inlet with all the discs in place in each stage. Daily samples were collected for each of the four stages, equaling 36,000 total L of air per sample. Discs rotated for approximately 24 days before the sampling per strip was complete and changed to a new disc. Discs were stored frozen after collection. Following collection and transport from AMF3 to PI Creamean's laboratory in Boulder, CO, strips were removed from the discs and cut up to separate daily samples. Each daily portion of the strip was stored in separate sterile bags and stored in a standard commercial freezer until analysis (approximately 4 – 5 months after collection). Drop freezing assays (DFAs) were conducted on one daily sample per week of the study for all four stages, with daily samples analyzed from 22 - 29 May 2017. Immediately after removal from the freezer, 2 mL of ultrapure water was added with a single-use sterile syringe directly to the bags, then mixed using a vortex mixer for 2 hours at 500 rpm to resuspend particles from the PFA into the ultrapure water. Following sample preparation, another single-use syringe was used to create 100 aliquots on a 3-inch diameter copper plate, then covered with a transparent plastic dome. The drops were approximately the same volume (2.5 µL each), though careful inspection. The plate was cooled at approximately 3 – 10 K min-1 from room temperature until all drops froze on the plate. Temperature was measured by an Omega thermocouple probes and meter. Drop freezing was visually detected, but recorded through monitoring software to provide the time frozen, channel of the meter used, temperature, and cooling rate for each drop. Not all 100 drops were always detected; the total number of rows in each data file equals the number of drops recorded, which typically was > 80 %. Each sample was tested three times (i.e., a new set of 100 drops was created for each test). Data provided are files containing these parameters. From the freezing temperatures recorded, the fraction frozen can be calculated as can the estimated INP concentration per L of air using the equation from Vali (1971): [INPs] = −(ln(f)/Valiquot)*(Vsuspension/Vair) where f is the proportion of droplets not frozen, Valiquot is the volume of each aliquot (i.e., 2.5e-6 L), Vsuspension is the volume of the suspension (i.e., 0.002 L), and Vair is the volume of air per sample (i.e., 36000 L). Vali, G.: Quantitative Evaluation of Experimental Results an the Heterogeneous Freezing Nucleation of Supercooled Liquids, Journal of the Atmospheric Sciences, 28, 402-409, 1971.

4-stage DRUM aerosol collector (DA400),Aerosol imp↗

SAR terrain classifier and mapper of biophysical attributes

In preparation for the launch of SIR-C/X-SAR and design studies for future orbital SAR, a program has made considerable progress in the development of an SAR terrain classifier and algorithms for quantification of biophysical attributes. The goal of this program is to produce a generalized software package for terrain classification and estimation of biophysical attributes and to make this package available to the larger scientific community. The basic elements of the SAR (Synthetic Aperture Radar) terrain classifier are outlined. An SAR image is calibrated with respect to known system and processor gains and external targets (if available). A Level 1 classifier operates on the data to differentiate: urban features, surfaces and tall and short vegetation. Level 2 classifiers further subdivide these classes on the basis of structure. Finally, biophysical and geophysical inversions are applied to each class to estimate attributes of interest. The process used to develop the classifiers and inversions is shown. Radar scattering models developed from theory and from empirical data obtained by truck-mounted polarimeters and the JPL AirSAR are validated. The validated models are used in sensitivity studies to understand the roles of various scattering sources (i.e., surface trunk, branches, etc.) in determining net backscatter. Model simulations of sigma (sup o) as functions of the wave parameters (lambda, polarization and angle of incidence) and the geophysical and biophysical attributes are used to develop robust classifiers. The classifiers are validated using available AirSAR data sets. Specific estimators are developed for each class on the basis of the scattering models and empirical data sets. The candidate algorithms are tested with the AirSAR data sets. The attributes of interest include: total above ground biomass, woody biomass, soil moisture and soil roughness.

Ulaby, Fawwaz T.↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

58 GEOSCIENCES↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗

Greybox Thermal Parameter Identification of Electric Machine Stators

The parameters of electric machine thermal equivalent circuit networks are difficult to predict due to material and manufacturing uncertainties. In this paper, a Greybox system identification approach is used to identify parameters of electric machine stator lumped parameter thermal networks (LPTNs). LPTNs provide a low order, computationally efficient, dynamic model of temperatures at specific locations. Second and third order LPTN model structures are defined as state space equations with stator thermal parameters to be identified. To test the Greybox electric machine stator thermal system identification, five stator motorette prototypes were constructed with controlled variations in slot fill and slot liner thickness. The variation in the motorette thermal parameters and thermal time constants are detected using the Greybox identification. Special attention is given to the impact of sampling rate and Greybox data record length on parameter estimation accuracy.

33 ADVANCED PROPULSION SYSTEMS↗

Discovering the Unknowns: A First Step

This article aims at discovering the unknown variables in the system through data analysis. The main idea is to use the time of data collection as a surrogate variable and try to identify the unknown variables by modeling gradual and sudden changes in the data. We use Gaussian process modeling and a sparse representation of the sudden changes to efficiently estimate the large number of parameters in the proposed statistical model. The method is tested on a realistic dataset generated using a one-dimensional implementation of a Magnetized Liner Inertial Fusion (MagLIF) simulation model, and encouraging results are obtained.

42 ENGINEERING↗

Mechanistic within-host mathematical model of inhalational anthrax

We present a mathematical model of the dynamics of Bacillus anthracis bacteria within the lymph nodes and blood of a host, following inhalation of an initial dose of spores. We also incorporate the dynamics of protective antigen, which is the binding component of the anthrax toxin produced by the bacteria. The model offers a mechanistic description of the early infection dynamics of inhalational anthrax, while its stochastic nature allows us to study the probabilities of different outcomes (for example, how likely it is that the infection will be cleared for a given inhaled dose of spores) in order to explain dose-response data for inhalational anthrax. The model is calibrated via a Bayesian approach, using in vivo data from New Zealand white rabbit and guinea pig infection studies, enabling within-host parameters to be estimated. We also leverage incubation-period data from the Sverdlovsk 1979 anthrax outbreak to show that the model can accurately describe human time-to-symptoms data under reasonable parameter regimes. Finally, we derive a simple approximate formula for the probability of symptom onset before time t, assuming that the number of inhaled spores has a Poisson distribution.

59 BASIC BIOLOGICAL SCIENCES↗

Impact of Water Use by Utility-Scale Solar on Groundwater Resources of the Chuckwalla Basin, CA: Final Modeling Report

The DOE and BLM identified 285,000 acres of desert land in the Chuckwalla valley in the western U.S., for solar energy development. In addition to several approved solar projects, a pumped storage project was recently proposed to pump nearly 8000 acre-ft-yr of groundwater to store and stabilize solar energy output. This study aims at providing estimates of the amount of naturally-occurring recharge, and to estimate the impact of the pumping on the water table. To better provide the locations and intensity of natural recharge, this study employs an integrated, physically-based hydrologic model, PAWS+CLM, to calculate recharge. Then, the simulated recharge is used in a parameter estimation package to calibrate spatially-distributed K field. This design is to incorporate all available observational data, including soil moisture monitoring stations, groundwater head, and estimates of groundwater conductivity, to constrain the modeling. To address the uncertainty of the soil parameters, an ensemble of simulations are conducted, and the resulting recharges are either rejected or accepted based on calibrated groundwater head and local variation of the K field. The results indicate that the natural total inflow to the study domain is between 7107 and 12,772 afy. During the initial-fill phase of pumped storage project, the total outflow exceeds the upper bound estimate of the inflow. If the initial-fill is annualized to 20 years, the average pumping is more than the lower bound of inflows. The results indicate after adding the pumped storage project, the system will nearing, if not exceeding, its maximum renewable pumping capacity. The accepted recharges lead to a drawdown range of 24 to 45 ft for an assumed specific yield of 0.05. However, the drawdown is sensitive to this parameter, whereas there is insufficient data to adequately constrain this parameter.

14 SOLAR ENERGY↗

Machine-Learning-Assisted Hybrid Earth System Modelling

We propose a novel hybrid modeling approach that combines a physics-based, numerical model of the Earth system with a machine learning model. We also envision to integrate the hybrid model training with a data assimilation cycle, so that new observations update both the state estimate used as initial conditions and the machine learning model parameters.

58 GEOSCIENCES↗

Xanthos – A Global Hydrologic Model

Xanthos is an open-source hydrologic model, written in Python, designed to quantify and analyse global water availability. Xanthos simulates historical and future global water availability on a monthly time step at a spatial resolution of 0.5 geographic degrees. Xanthos was designed to be extensible and used by scientists that study global water supply and work with the Global Change Assessment Model (GCAM). Xanthos uses a user-defined configuration file to specify model inputs, outputs and parameters. Xanthos has been tested using actual global data sets and the model is able to provide historical observations and future estimates of renewable freshwater resources in the form of total runoff.

58 GEOSCIENCES↗

The effects of correlation on goodness of fit

Autocorrelated normal random variates were generated via computer and the effects of various levels of correlation on goodness of fit problems were studied. The results are useful in determining the distribution or estimating the parameters of populations that correlate observations such as wind speeds and temperature. The model used to generate the autocorrelated data is an autoregressive process of order 1. The Kolmogorov-Smirnov and chi-square statistics are used in the analysis. It was observed in the simulation that high positive correlations tend to shift the sample mean away from the population mean and negative correlations tend to shift the sample mean towards the population mean. In many cases, it was observed that positive and negative correlations tend to decrease the standard deviation. However, since this did not occur in all cases, no definite conclusion can be made regarding the standard deviation. Since the autoregressive process is a linear transformation, it is not surprising that normality was preserved. However, a possible extension of this problem could be to generate non-normal data and observe how the distribution is affected by correlation. Another extension might utilize another model such as autoregressive of order k or a moving average process of order k.

Kitchens, L. J.↗

The application of remote sensing to the development and formulation of hydrologic planning models: Executive summary

Methods for the reduction of remotely sensed data and its application in hydrologic land use assessment, surface water inventory, and soil property studies are presented. LANDSAT data is used to provide quantitative parameters and coefficients to construct watershed transfer functions for a hydrologic planning model aimed at estimating peak outflow from rainfall inputs.

Castruccio, P. A.↗

Infinite-dimensional approach to system identification of Space Control Laboratory Experiment (SCOLE)

The identification of a unique set of system parameters in large space structures poses a significant new problem in control technology. Presented is an infinite-dimensional identification scheme to determine system parameters in large flexible structures in space. The method retains the distributed nature of the structure throughout the development of the algorithm and a finite-element approximation is used only to implement the algorithm. This approach eliminates many problems associated with model truncation used in other methods of identification. The identification is formulated in Hilbert space and an optimal control technique is used to minimize weighted least squares of error between the actual and the model data. A variational approach is used to solve the problem. A costate equation, gradients of parameter variations and conditions for optimal estimates are obtained. Computer simulation studies are conducted using a shuttle-attached antenna configuration, more popularly known as the Space Control Laboratory Experiment (SCOLE) as an example. Numerical results show a close match between the estimated and true values of the parameters.

Hossain, S. A.↗

Statistical analysis of censored motion sickness latency data using the two-parameter Weibull distribution

The suitability of the two-parameter Weibull distribution for describing highly censored cat motion sickness latency data was evaluated by estimating the parameters with the maximum likelihood method and testing for goodness of fit with the Kolmogorov-Smirnov statistic. A procedure for determining confidence levels and testing for significance of the difference between Weibull parameters is described. Computer programs for these procedures may be obtained from an archival source.

Park, Won J.↗

The Carina dwarf spheroidal galaxy - How dark is it?

Precise radial velocities obtained with a photon-counting echelle spectrograph for a sample of 17 red giants in the Carina dwarf spheroidal galaxy are presented. The calculation of the systemic velocity and central velocity dispersion of Carina is described, the existing data constraining the structural parameters of Carina are reviewed, and an estimate of the central surface brightness of the galaxy is derived. These data are used to estimate the central mass density of Carina, as well as central and global mass-to-light ratios. It is concluded that the inferred mass densities and mass-density limits for all acceptable models imply the presence of a significant DM component in Carina. DM properties of all well-studied dSph systems are summarized and compared.

Mateo, Mario↗

SOFIA Closed- and Open-Door Aerodynamic Analyses

Work to evaluate the aerodynamic characteristics and the cavity acoustic environment of the SOFIA (Stratospheric Observatory for Infrared Astronomy) airplane has been completed. The airplane has been evaluated in its closed-door configuration, as well as several open-door configurations. Work performed included: acoustic analysis tool development, cavity acoustic evaluation, stability and control parameter estimation, air data calibration, and external flow evaluation. Qualitative airflow data were obtained during the closed- and open-door flights using tufts on the aft portion of the fuselage. Video was taken from a chase plane. This video was analyzed for various flight conditions, and general flow descriptions of the aft fuselage of the 747SP were developed for the different closed and open door configurations.

Cumming, Stephen↗

Autonomous Modal Identification of the Space Shuttle Tail Rudder

Autonomous modal identification automates the calculation of natural vibration frequencies, damping, and mode shapes of a structure from experimental data. This technology complements damage detection techniques that use continuous or periodic monitoring of vibration characteristics. The approach shown in the paper incorporates the Eigensystem Realization Algorithm (ERA) as a data analysis engine and an autonomous supervisor to condense multiple estimates of modal parameters using ERA's Consistent-Mode Indicator and correlation of mode shapes. The procedure was applied to free-decay responses of a Space Shuttle tail rudder and successfully identified the seven modes of the structure below 250 Hz. The final modal parameters are a condensed set of results for 87 individual ERA cases requiring approximately five minutes of CPU time on a DEC Alpha computer.

Pappa, Richard S.↗