SEARCH · Engineering Papers

Results for “Model error”

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 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Extracting galactic structure parameters from multivariated density estimation

Multivariate statistical analysis, including includes cluster analysis (unsupervised classification), discriminant analysis (supervised classification) and principle component analysis (dimensionlity reduction method), and nonparameter density estimation have been successfully used to search for meaningful associations in the 5-dimensional space of observables between observed points and the sets of simulated points generated from a synthetic approach of galaxy modelling. These methodologies can be applied as the new tools to obtain information about hidden structure otherwise unrecognizable, and place important constraints on the space distribution of various stellar populations in the Milky Way. In this paper, we concentrate on illustrating how to use nonparameter density estimation to substitute for the true densities in both of the simulating sample and real sample in the five-dimensional space. In order to fit model predicted densities to reality, we derive a set of equations which include n lines (where n is the total number of observed points) and m (where m: the numbers of predefined groups) unknown parameters. A least-square estimation will allow us to determine the density law of different groups and components in the Galaxy. The output from our software, which can be used in many research fields, will also give out the systematic error between the model and the observation by a Bayes rule.

Chen, B.↗

DOE OSTI2022

Sensitivity of Cloud-Radiative Effects to Cloud Fraction Parametrizations in Tropical, Mid-Latitude and Arctic Kilometre-Scale Simulations

The Regional Atmosphere (RA) configuration of the Met Office Unified Model currently requires different cloud fraction parametrizations (CFP) for tropical and midlatitude simulations. To explore the scope for unification of these two RA configurations, this paper presents a detailed evaluation of simulations over tropical, midlatitude and arctic domains, with two different diagnostic CFPs, a prognostic CFP, and no CFP at all. Furthermore, a novel, hybrid approach was used that treats liquid cloud diagnostically and ice cloud prognostically. Using observations from three U.S. Department of Energy Atmospheric Radiation Measurement supersites, it is shown that none of these CFPs stands out as superior over all domains. Over the frequently overcast Arctic, the all-or-nothing approach best captures the cloud radiative properties. Conversely, CFPs are of benefit in regions with frequent partial cloudiness, such as the midlatitudes and the tropics. However, their improved cloud radiative properties often hide an error compensation. All models underestimate overcast, low-base cloud with small water paths in convective environments. In addition, mid-latitude overcast, low-base, optically thick clouds in the morning, possibly associated with overnight convection, are frequently too broken. Diagnostic schemes compensate for these errors by producing spurious, scattered afternoon cloud, which could be due to a correct cloud response to too eager convective initiation. Winter clouds over the midlatitudes are improved when liquid cloud is represented diagnostically with a bimodal saturation-departure PDF, without error compensation. While it is difficult to unify the RA across the globe around a single CFP scheme, the newly proposed Hybrid scheme performs reasonably well for cloud cover across all regions. It also exhibits SW biases that are smaller than most other configurations and is less affected by excessive liquid water paths and compensating errors than fully diagnostic schemes. Surface precipitation is fairly insensitive to the CFP in the simulations shown here.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2017

Program to Optimize Simulated Trajectories II (POST2) Surrogate Models for Mars Ascent Vehicle (MAV) Performance Assessment

The primary purpose of the multiPOST tool is to enable the execution of much larger sets of vehicle cases to allow for broader trade space exploration. However, this exploration is not achieved solely with the increased case throughput. The multiPOST tool is applied to carry out a Design of Experiments (DOE), which is a set of cases that have been structured to capture a maximum amount of information about the design space with minimal computational effort. The results of the DOE are then used to fit a surrogate model, ultimately enabling parametric design space exploration. The approach used for the MAV study includes both DOE and surrogate modeling. First, the primary design considerations for the vehicle were used to develop the variables and ranges for the multiPOST DOE. The final set of DOE variables were carefully selected in order to capture the desired vehicle trades and take into account any special considerations for surrogate modeling. Next, the DOE sets were executed through multiPOST. Following successful completion of the DOE cases, a manual verification trial was performed. The trial involved randomly selecting cases from the DOE set and running them by hand. The results from the human analyst's run and multiPOST were then compared to ensure that the automated runs were being executed properly. Completion of the verification trials was then followed by surrogate model fitting. After fits to the multiPOST data were successfully created, the surrogate models were used as a stand-in for POST2 to carry out the desired MAV trades. Using the surrogate models in lieu of POST2 allowed for visualization of vehicle sensitivities to the input variables as well as rapid evaluation of vehicle performance. Although the models introduce some error into the output of the trade study, they were very effective at identifying areas of interest within the trade space for further refinement by human analysts. The next section will cover all of the ground rules and assumptions associated with DOE setup and multiPOST execution. Section 3.1 gives the final DOE variables and ranges, while section 3.2 addresses the POST2 specific assumptions. The results of the verification trials are given in section 4. Section 5 gives the surrogate model fitting results, including the goodness-of-fit metrics for each fit. Finally, the MAV specific results are discussed in section 6.

Zwack, M. R.↗

NASA NTRS1999

Comparing Parameter Estimation Techniques for an Electrical Power Transformer Oil Temperature Prediction Model

This paper examines various sources of error in MIT's improved top oil temperature rise over ambient temperature model and estimation process. The sources of error are the current parameter estimation technique, quantization noise, and post-processing of the transformer data. Results from this paper will show that an output error parameter estimation technique should be selected to replace the current least squares estimation technique. The output error technique obtained accurate predictions of transformer behavior, revealed the best error covariance, obtained consistent parameter estimates, and provided for valid and sensible parameters. This paper will also show that the output error technique should be used to minimize errors attributed to post-processing (decimation) of the transformer data. Models used in this paper are validated using data from a large transformer in service.

Morris, A. Terry↗

NASA NTRS1993

The North Atlantic circulation in the early 1980s - An estimate from inversion of a finite-difference model

A finite-difference model of the North Atlantic is constructed for the purpose of making an estimate of the circulation through an inverse calculation. The data base is eclectic, and includes hydrography, oxygen, nutrients, current meter and float records, atmospheric momentum, heat and water vapor transfers, as well as estimates of certain integral fluxes. Owing to the available hydrographic database, the model resolution is restricted to 1 deg at best, and is much coarser in many aspects. This limited resolution is a major obstacle to accurate estimates of climatological fluxes. In its final form, there are about 9000 constraints in 29,000 formal unknowns plus 9000 noise unknowns. The system is solved as a tapered least-squares system by a sparse conjugate gradient algorithm. With the exception of a few float velocities, all constraints are found to be consistent within error estimates. The model produces estimates of large-scale fluxes and flux divergences for all conventional properties including heat and nutrients as well as carbon dioxide and alkalinity. Meridional fluxes of carbon are found to be indistinguishable from zero, whereas the North Atlantic tends to export nutrients to the south, but carry heat to the north. Traditional oceanographic depictions of the circulation through combination of nonsynoptic data into steady models may have reached their useful limit in the present calculation, as the conflicts between the data and physical requirements become quantitatively apparent.

Martel, Francoise↗

NASA NTRS2018

Preliminary Jitter Stability Results for the Large UV/Optical/Infrared (LUVOIR) Surveyor Concept Using a Non-Contact Vibration Isolation and Precision Pointing System

The need for high payload dynamic stability and ultra-stable mechanical systems is an overarching technology need for large space telescopes such as the Large Ultraviolet / Optical / Infrared (LUVOIR) Surveyor concept. The LUVOIR concept includes a 15-meter-diameter segmented-aperture telescope with a suite of serviceable instruments operating over a range of wavelengths between 100nm to 2.5 um. Wavefront error (WFE) stability of less than 10 picometers RMS of uncorrected system WFE per wavefront control step represents a drastic performance improvement over current space-based telescopes being fielded. Through the utilization of an isolation architecture that involves no mechanical contact between the telescope and the host spacecraft structure, a system design is realized that maximizes the telescope dynamic stability performance without driving stringent technology requirements on spacecraft structure, sensors or actuators. Through analysis of the LUVOIR finite element model and linear optical model, the wavefront error and Line-Of-Sight (LOS) jitter performance is discussed in this paper when using the Vibration Isolation and Precision Pointing System (VIPPS) being developed cooperatively with Lockheed Martin in addition to a multi-loop control architecture. The multi-loop control architecture consists of the spacecraft Attitude Control System (ACS), VIPPS, and a Fast Steering Mirror on the instrument. While the baseline attitude control device for LUVOIR is a set of Control Moment Gyroscopes (CMGs), Reaction Wheel Assembly (RWA) disturbance contribution to wavefront error stability and LOS stability are presented to give preliminary results in this paper. CMG disturbance will be explored in further work to be completed.

Sacks, Lia W.↗

NASA NTRS1991

Progress in the development of PDF turbulence models for combustion

A combined Monte Carlo-computational fluid dynamic (CFD) algorithm was developed recently at Lewis Research Center (LeRC) for turbulent reacting flows. In this algorithm, conventional CFD schemes are employed to obtain the velocity field and other velocity related turbulent quantities, and a Monte Carlo scheme is used to solve the evolution equation for the probability density function (pdf) of species mass fraction and temperature. In combustion computations, the predictions of chemical reaction rates (the source terms in the species conservation equation) are poor if conventional turbulence modles are used. The main difficulty lies in the fact that the reaction rate is highly nonlinear, and the use of averaged temperature produces excessively large errors. Moment closure models for the source terms have attained only limited success. The probability density function (pdf) method seems to be the only alternative at the present time that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus may be the only viable approach for more accurate turbulent combustion calculations. Assumed pdf's are useful in simple problems; however, for more general combustion problems, the solution of an evolution equation for the pdf is necessary.

Hsu, Andrew T.↗

NASA NTRS1996

Utilization of STRAT Data

One of the goals of the STRAT mission was to improve our understanding of transport in the stratosphere. We have used our coupled 2-D model to analyse transport between the tropics and midlatitudes using aircraft measurements of (CO2) and (N2O) made during the SPADE and ASHOE/MAESA campaigns (Boering et al. 1994) for model validation. We have further used the model to estimate errors in the experimental determination of the age of stratospheric air that arise from nonlinear increases of real tracers such as (CO2) and (SF6). The interpretation of measurements made during STRAT is still in progress. Most of the work performed under this contract has been described in the second year progress report, submitted in April of 1996. In the following, the research is summarized and some additional work that has become relevant in the meantime is described. We have developed an interactive model of the dynamics, chemistry and radiation of the stratosphere with partial support from this contract. The dynamics module integrates the primitive equations on a sphere. Lower boundary zonal winds and eddy geopotential heights are specified from observations. Chemistry and tracer transport are done in two dimensions. The chemistry model is a family model. Long lived species are transported and short lived species are determined in intervals of 5 to 10 days by assuming photochemical steady state. Full diurnal integrations with the fast species are performed periodically and diurnal coefficients, derived from the explicit integration are used to determine the steady state solution. Heating rates in the stratosphere are calculated from model ozone and the model temperatures and circulation fields are used in the chemistry transport module.

Source record↗

DOE OSTI2020

A data-driven model for thermodynamic properties of a steam generator under cycling operation

The varying electricity demand from coal power plants due to the intermittent nature of renewable sources leads to load-follow and on/off operations referred to as cycling. Cycling causes transients of properties such as pressure and temperature within various components of the steam generation system.These transients cause increased damage because of fatigue and creep-fatigue interactions shortening the life of components. An algorithm is developed to identify cycling operations from the gross power data. The data-driven model based on artificial neural networks (ANN) is developed using 10 years data from Coal Creek Station power plant located in North Dakota, USA to estimate properties of the steam generator components during cycling operations. Furthermore, the uniqueness of this model is the ability to predict component properties for the cycling as well as base-load operations and is reported for the first time. The ANN model estimates the component properties, for a given gross power profile and initial conditions, as they vary during cycling operations. As a representative example, the ANN estimates are presented for the superheater outlet pressure, reheater inlet temperature, and flue gas temperature at the air heater inlet. The changes in these variables as a function of the gross power over the time duration are compared with measurements to assess the predictive capability of the model. Mean square errors of 4.49E-04 for superheater outlet pressure, 1.62E-03 for reheater inlet temperature, and 4.14E-04 for flue gas temperature at the air heater inlet were observed.

01 COAL, LIGNITE, AND PEAT↗

NASA NTRS1987

Modeling and control of beam-like structures

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for beam-like continua. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

NASA NTRS1994

Comprehensive stellar population models and the disentanglement of age and metallicity effects

The construction of detailed models for intermediate and old stellar populations is described. Input parameters include metallicity (-2 less than (Fe/H) less than 0.5), single-burst age (between 1.5 and 17 Gyr), and initial mass function (IMF) exponent. Quantities output include broadband magnitudes, spectral energy distributions, surface brightness fluctuation magnitudes, and a suite of 21 absorption feature indices. The models are checked against a wide variety of available observations. Examinations of model output yield the following conclusions. (1) If the percentage change delta age/delta Z approximately equals 3/2 for two populations, they will appear almost identical in most indices. A few indices break this degeneracy by being either more abundance sensitive (Fe4668, Fe5015, Fe5709, and Fe5782) or more age sensitive (G4300, H beta, and presumably higher order Balmer lines) than usual. (2) Present uncertainties in stellar evolution are of the same magnitude as the effects of IMF and Y in the indices studied. (3) Changes in abundance ratios (like (Mg/Fe)) are predicted to be readily apparent in the spectra of old stellar populations. (4) The I-band flux of a stellar population is predicted to be nearly independent of metallicity and only modestly sensitive to age. The I band is therefore recommended for standard candle work or studies of M/L in galaxies. Other conclusions stem from this work. (1) Intercomparison of models and observations of two TiO indices seem to indicate variation of the (V/Ti) ratio among galaxies, but it is not clear how this observation ties into the standard picture of chemical enrichment. (2) Current estimates of (Fe/H) for the most metal-rich globulars that are based on integrated indices are probably slightly too high. (3) Colors of population models from different authors exhibit a substantial range. At solar metallicity and 13 Gyr, this range corresponds to an age error of roughly +/- 7 Gyr. Model colors from different authors applied in a differential sense have smaller uncertainties. (4) In the present models the dominant error for colors is probably the transformation from stellar atmospheric parameters to stellar colors. (5) Stellar B - V is difficult to model, and current spreads among different authors can reach 0.2 mag. (6) If known defects in the stellar flux library are corrected, the population model colors of this work in passbands redder than U would be accurate to roughly 0.03 mag in an absolute sense. These corrections are not made in the tables of model output.

Worthey, Guy↗

DOE OSTI2020

Tensor Basis Gaussian Process Models of Hyperelastic Materials

In this study, we develop Gaussian process regression (GPR) models of isotropic hyperelastic material behavior. First, we consider the direct approach of modeling the components of the Cauchy stress tensor as a function of the components of the Finger stretch tensor in a Gaussian process. We then consider an improvement on this approach that embeds rotational invariance of the stress-stretch constitutive relation in the GPR representation. This approach requires fewer training examples and achieves higher accuracy while maintaining invariance to rotations exactly. Finally, we consider an approach that recovers the strain-energy density function and derives the stress tensor from this potential. Although the error of this model for predicting the stress tensor is higher, the strain-energy density is recovered with high accuracy from limited training data. The approaches presented here are examples of physics-informed machine learning. They go beyond purely data-driven approaches by embedding the physical system constraints directly into the Gaussian process representation of materials models.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Accuracy of hourly energy predictions for demand flexibility applications

Decarbonization goals in the United States electricity sector are increasing the levels of renewable energy generation in the electricity supply system, and are driving increased attention to building electrification, which will increase the magnitude and shift the timing of the electricity system peak. These changes are motivating new approaches to coordinate building electricity demand with low-carbon renewable generation, elevating the importance of demand flexibility (DF) in buildings and the need to quantify the temporal impacts of DF. In this paper, we first characterize the hourly predictive accuracy of six commonly used baseline models in an application context of quantifying building-level load shift. Our analysis revealed insights such as hours of the day (afternoons), periods of the week (weekends), and seasons (summer) that were predicted with more accuracy than other time periods. In addition, the analysis showed tendencies toward overprediction or underprediction of load. Secondly, we provide the first published investigation of baseline erosion from repeated dispatch of building load shifting. We observed that as the baseline period is pushed back further from the prediction day, the distribution of errors across baseline model predictions increases, with notable inflection points near the three-week erosion point for two of the three models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2011

Orion Handling Qualities During ISS Rendezvous and Docking

The Orion spacecraft was designed to rendezvous with multiple vehicles in low earth orbit (LEO) and beyond. To perform the required rendezvous and docking task, Orion must provide enough control authority to perform coarse translational maneuvers while maintaining precision to perform the delicate docking corrections. While Orion has autonomous docking capabilities, it is expected that final approach and docking operations with the International Space Station (ISS) will initially be performed in a manual mode. A series of evaluations was conducted by NASA and Lockheed Martin at the Johnson Space Center to determine the handling qualities (HQ) of the Orion spacecraft during different docking and rendezvous conditions using the Cooper-Harper scale. This paper will address the specifics of the handling qualities methodology, vehicle configuration, scenarios flown, data collection tools, and subject ratings and comments. The initial Orion HQ assessment examined Orion docking to the ISS. This scenario demonstrates the Translational Hand Controller (THC) handling qualities of Orion. During this initial assessment, two different scenarios were evaluated. The first was a nominal docking approach to a stable ISS, with Orion initializing with relative position dispersions and a closing rate of approximately 0.1 ft/sec. The second docking scenario was identical to the first, except the attitude motion of the ISS was modeled to simulate a stress case ( 1 degree deadband per axis and 0.01 deg/sec rate deadband per axis). For both scenarios, subjects started each run on final approach at a docking port-to-port range of 20 ft. Subjects used the THC in pulse mode with cues from the docking camera image, window views, and range and range rate data displayed on the Orion display units. As in the actual design, the attitude of the Orion vehicle was held by the automated flight control system at 0.5 degree deadband per axis. Several error sources were modeled including Reaction Control System (RCS) jet angular and position misalignment, RCS thrust magnitude uncertainty, RCS jet force direction uncertainty due to self plume impingement, and Orion center of mass uncertainty.

Hart, Jeremy J.↗

NASA NTRS1988

An alternative to Guyan reduction of finite-element models

Structural modeling is a key part of structural system identification for large space structures. Finite-element structural models are commonly used in practice because of their general applicability and availability. The initial models generated by using a standard computer program such as NASTRAN, ANSYS, SUPERB, STARDYNE, STRUDL, etc., generally contain tens of thousands of degrees of freedom. The models must be reduced for purposes of identification. Not only does the magnitude of the identification effort grow exponentially as a function of the number of degrees of freedom, but numerical procedures may also break down because of accumulated round-off errors. Guyan reduction is usually applied after a static condensation. Misapplication of Guyan reduction can lead to serious modeling errors. It is quite unfortunate and disappointing, since the accuracy of the original detailed finite-element model one tries very hard to achieve is lost by the reduction. First, why and how Guyan reduction always causes loss of accuracy is examined. An alternative approach is then introduced. The alternative can be thought of as an improvement of Guyan reduction, the Rayleigh-Ritz method, and in particular the recent algorithm of Wilson, Yuan, and Dickens. Unlike Guyan reduction, the use of the alternative does not need any special insight, experience, or skill for partitioning the structural degrees of freedom. In addition to model condensation, this alternative approach can also be used for predicting analytically, quickly, and economically, what are those structural modes that are excitable by a force actuator at a given trial location. That is, in the excitation of the structural modes for identification, it can be used for guiding the placement of the force actuators.

Lin, Jiguan Gene↗

NASA NTRS1985

A comparative study of spherical and flat-Earth geopotential modeling at satellite elevations

Flat-Earth modeling is a desirable alternative to the complex spherical-Earth modeling process. These methods were compared using 2 1/2 dimensional flat-earth and spherical modeling to compute gravity and scalar magnetic anomalies along profiles perpendicular to the strike of variably dimensioned rectangular prisms at altitudes of 150, 300, and 450 km. Comparison was achieved with percent error computations (spherical-flat/spherical) at critical anomaly points. At the peak gravity anomaly value, errors are less than + or - 5% for all prisms. At 1/2 and 1/10 of the peak, errors are generally less than 10% and 40% respectively, increasing to these values with longer and wider prisms at higher altitudes. For magnetics, the errors at critical anomaly points are less than -10% for all prisms, attaining these magnitudes with longer and wider prisms at higher altitudes. In general, in both gravity and magnetic modeling, errors increase greatly for prisms wider than 500 km, although gravity modeling is more sensitive than magnetic modeling to spherical-Earth effects. Preliminary modeling of both satellite gravity and magnetic anomalies using flat-Earth assumptions is justified considering the errors caused by uncertainties in isolating anomalies.

Parrott, M. H.↗