Engineering Papers⌕ Search

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.

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

At least 1,621 records · Page 90

Benchmarking performance changes in the simulation of extratropical modes of variability across CMIP generations

In this work, we evaluate extratropical modes of variability in the three most recent phases of the Coupled Model Intercomparison Project (CMIP3, 5, and 6) to gauge improvement of climate models over time. A suite of high-level metrics is employed to objectively evaluate how well climate models simulate the qobserved Northern Annular Mode (NAM), North Atlantic Oscillation (NAO), Pacific North America pattern (PNA), Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO), North Pacific Oscillation (NPO), and North Pacific Gyre Oscillation (NPGO). We apply a Common Basis Function (CBF) approach that projects model anomalies onto observed empirical orthogonal functions (EOFs), together with the traditional EOF approach, to CMIP Historical and AMIP models. We find simulated spatial patterns of those modes have been significantly improved in the newer models, although the skill improvement is sensitive to the mode and season considered. We identify some potential contributions to the pattern improvement of certain modes (e.g., the Southern Hemisphere jet and High-top vertical coordinate), however the performance changes are likely attributed to gradual improvement of the base climate and multiple relevant processes. Less performance improvement is evident in the mode amplitude of these modes and systematic overestimation of the mode amplitude in spring remains in the newer climate models. We find that the post-dominant season amplitude errors in atmospheric modes are not limited to coupled runs but are often already evident in AMIP simulations. This suggests that rectifying the egregious post-dominant season amplitude errors found in many models can be addressed in an atmospheric-only framework making it more tractable to address in the model development process.

54 ENVIRONMENTAL SCIENCES↗

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

Assessment of recent resolution and parameterization changes in the GLA fourth order GCM

The Goddard Laboratory for Atmospheres' fourth-order GCM is under evaluation for the impact on model integrations of enhanced horizontal and vertical resolution, as well as the effects of such novel parameterization schemes as that of gravity-wave-drag, the Arakawa-Schubert (1974) cumulus parameterization, and an explicitly-resolved planetary boundary layer. While the doubling of the GMC's horizontal resolution to 2 deg in latitude and 2.5 deg in longitude has improved the model's predictive skill for 6-7 day forecasts, systematic errors associated with the model's climate drift lead to a deterioration in predictions for longer forecasts.

Helfand, H. M.↗

A New Global Regression Analysis Method for the Prediction of Wind Tunnel Model Weight Corrections

A new global regression analysis method is discussed that predicts wind tunnel model weight corrections for strain-gage balance loads during a wind tunnel test. The method determines corrections by combining "wind-on" model attitude measurements with least squares estimates of the model weight and center of gravity coordinates that are obtained from "wind-off" data points. The method treats the least squares fit of the model weight separate from the fit of the center of gravity coordinates. Therefore, it performs two fits of "wind- off" data points and uses the least squares estimator of the model weight as an input for the fit of the center of gravity coordinates. Explicit equations for the least squares estimators of the weight and center of gravity coordinates are derived that simplify the implementation of the method in the data system software of a wind tunnel. In addition, recommendations for sets of "wind-off" data points are made that take typical model support system constraints into account. Explicit equations of the confidence intervals on the model weight and center of gravity coordinates and two different error analyses of the model weight prediction are also discussed in the appendices of the paper.

strain-gage balance↗

Mass Inferencing Model Creation And Deployment To Lunar Excavation Robot, RASSOR

NASA’s Regolith Advanced Surface Systems Operations Robot (RASSOR) Excavator is a teleoperated mobile robotic platform with a unique space regolith excavation capability. This research project developed functionality for inferencing regolith mass ingested during RASSOR operation, enhancing RASSOR’s ability to successfully complete ISRU missions. Radio wave propagation time to the Moon and back is ~2.56 seconds. Though teleoperation is possible with this delay, autonomous capability that enables RASSOR to plan and execute excavation missions intelligently and efficiently is preferred. To teleoperate or run autonomously, it is crucial for the quantity of regolith mass ingested by RASSOR to be available as a system state for efficient operation (e.g. knowledge of whether drums are full informs the task of highest priority, whether it be continuing to dig, or returning to a processing plant to offload regolith). A configurable data reduction and analysis pipeline was created to allow for straightforward incorporation of new data, such as that from lunar excavation, to improve model performance in new environments. Four distinct modeling approaches were employed in developing a mass inferencing approach that could work on RASSOR. All four models take in system states and output a mass prediction for each set of the robot’s bucket drums. Initial results from deployment to RASSOR and testing in a simulated lunar environment show that the models have <10% mean error during robot operation. Future work includes refinement of a model that estimates regolith mass in real-time during excavation as well as further testing of the developed models on the hardware.

ROS↗

Examples of Nonconservatism in the CARE 3 Program

This paper presents parameter regions in the CARE 3 (Computer-Aided Reliability Estimation version 3) computer program where the program overestimates the reliability of a modeled system without warning the user. Five simple models of fault-tolerant computer systems are analyzed; and, the parameter regions where reliability is overestimated are given. The source of the error in the reliability estimates for models which incorporate transient fault occurrences was not readily apparent. However, the source of much of the error for models with permanent and intermittent faults can be attributed to the choice of values for the run-time parameters of the program.

Dotson, Kelly J.↗

Anomaly detection for MPC forecast in Fleet of Water Heaters

Among residential devices, water heaters consume 20% of home energy use in the United States. Water heaters possess the capability to store energy within their reservoirs, enabling the ability to decouple energy use from hot water use. This capability can be used to reduce energy usage and costs while also supporting grid services. This requires accurate forecasting of the parameters of the water heater such as upper and lower temperatures. In this study, we analyzed the performance and behavior of a water heater model used in the real-world to predict a control mechanism that is implemented in a smart residential neighborhood. The model forecasts are accurate in most cases but not all. In such scenarios, error correction of the model is necessary to further improve model predictive control accuracy. Anomaly detection is the first step of error correction. This study complements existing research by grouping time series data into two clusters one with anomalies and another without anomalies. To achieve this task, we explored and compared multiple unsupervised machine learning algorithms to perform clustering. Among these algorithms, Ward clustering has the lowest running time and identified the highest number of anomalies for the upper temperature limit. The proposed approach is tested based on the data collected in a neighborhood with 46 townhomes located in Atlanta, GA.

Lebakula, Viswadeep↗

Internal Model-Based Active Damping Strategy for a Back-to-Back Modular Multilevel Converter System for Advanced Grid Support: Preprint

The proposed work focuses on the possibility of achieving the reduction in the size of the interfacing filter for the grid connection for a back-to-back modular multilevel converter system with the use of third order LCL filters. To achieve the same attenuation, it is possible to reduce the size of the interfacing filters by the usage of third order filters. However, this kind of filtering comes with the limitation of having sustained oscillations due to lack of damping especially during transient changes. Therefore, in this work, a methodology has been proposed based on the principle of internal model to cater for this unwanted oscillations. The third order system is modeled inside the microcontroller with the damping enabled. The error between the output from the model and the actual are compared and the error is utilized to accomplish the active damping strategy. The proposed architecture is verified via computer simulations based on MATLAB/Simulink domain and various case study results along with their discussion is presented in this paper.

dq control↗

The Separation and Hα Contrasts of Massive Accreting Planets in the Gaps of Transitional Disks: Predicted Hα Protoplanet Yields for Adaptive Optics Surveys

We present a massive accreting gap planet model that ensures large gaps in transitional disks are kept dust free by the scattering action of three coplanar quasi-circular planets in a 1:2:4 mean motion resonance (MMR). This model uses the constraint of the observed gap size, and the dust-free nature of the gap, to determine within ∼10% the possible orbits for three massive planets in an MMR. Calculated orbits are consistent with the observed orbits and Hα emission (the brightest line to observe these planets) for LkCa 15 b, PDS 70 b, and PDS 70 c within observational errors. Moreover, the model suggests that the scarcity of detected Hα planets is likely a selection effect of the current limitations of non-coronagraphic, low (<10%) Strehl, Hα imaging with adaptive optics (AO) systems used in past Hα surveys. We predict that as higher Strehl AO systems (with high-performance custom coronagraphs; like the 6.5 m Magellan Telescope MagAO-X system) are utilized at Hα, the number of detected gap planets will substantially increase by more than tenfold. For example, we show that >25 ± 5 new Hα “gap planets” are potentially discoverable by a survey of the best 19 transitional disks with MagAO-X. Detections of these accreting protoplanets will significantly improve our understanding of planet formation, planet growth and accretion, solar system architectures, and planet–disk interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Terrain Relative Navigation for Guided Descent on Titan

Titan’s dense atmosphere, low gravity, and high winds at high altitudes create descent times of >90 minutes with standard entry/descent/landing (EDL) architectures and result in large unguided landing ellipses, with 99% values of 110x110 km and 149x72 km in recent Titan lander proposals. Enabling precision landing on Titan could increase science return for the types of missions proposed to date and make additional types of landing sites accessible, opening up new possibilities for science investigations. Precision landing on Titan has unique challenges, because the hazy atmosphere makes it difficult to see the surface and because it requires guided descent with divert ranges that are one to two orders of magnitude larger than needed for other target bodies, i.e. up to on the order of 100 km. It is conceivable that such a divert capability could be provided economically by a parafoil or other steerable aerodynamic decelerator deployed several 10s of km above the surface. The long descent times lead to large inertial navigation errors, hence a need for terrain relative navigation (TRN). This would require a TRN capability that can operate at such altitudes, despite challenges of seeing the surface sufficiently clearly and of depending on map products that are two orders of magnitude lower in spatial resolution than those for Mars and airless bodies. We then develop algorithms for map matching and feature tracking with descent images and test these with synthetic images created from Cassini/Huygens data sets and our radiative transfer model. We also introduce new possibilities for TRN based on the potential to discriminate some specific types of terrain onboard in descent imagery, such as lake vs adjacent ground and dune vs interdune. We use sensor measurement noise models in simulations of state estimation with an extended Kalman filter that includes coordinates of a set of tracked features in the state vector. Case studies were done for two notional landing sites, one in a site with only dry ground and one in a Titan lake district. In both cases, the filter error model shows 3 position error at touchdown on the order of 2 km. More work is needed to validate these results with higher fidelity camera models and larger data sets, but this is very promising.

Matthies, Larry↗

Comparative Study on the Machine Learning-Based Prediction of Adsorption Energies for Ring and Chain Species on Metal Catalyst Surfaces

Computation of adsorption and transition state energies for a large number of surface intermediates for numerous active site models pose significant computational overhead in computational screening of catalysts. Machine learning (ML) techniques can be used to predict part of these energies. To predict the energies, ML models need to be fed appropriate metal and species descriptors. For complex surface chemistries, the structures of the intermediate species can vary greatly. In this paper, working with the hydrodeoxygenation of succinic acid on six different metal surfaces, we have studied the effect of linear and non-linear ML models used along with pen-and-paper based species descriptors and two categories of metal descriptors on two different categories of intermediate species: chain and ring. More specifically, our computations include the prediction of chain species when trained on only chain species and also when trained on both chain and ring species. Similar computations were performed for predictions of ring species. In each case, results of linear ML models were compared with kernel based non-linear models. Our results indicate that ring species data does not improve the prediction of chain species. Similarly, chain species data does not improve the prediction of ring species. The use of non-linear ML models, however, did help to minimize the prediction errors compared to the linear models. Furthermore, the study also shows that electronic or adsorption energy based metal descriptors along with bond count based species fingerprints can achieve a mean absolute error (MAE) of less than 0.2 eV for complex chain molecules when used with an appropriate machine learning model.

Adsorption↗

Dynamic Modeling of Near Isothermal Compressor for Transcritical Carbon Dioxide Cycle

Compressors are the major energy consumption components in vapor compression systems, drawing much research effort in reducing carbon emissions. The isothermal compressor integrates the compressor chamber and gas cooler to achieve near isothermal compression, reaching up to 30% energy reduction compared to the traditional isentropic compression work. This paper presents a detailed isothermal compressor model combined with a generalized liquid piston model to account for the carbon dioxide (CO2) isothermal compression process. The model is established based on MATLAB environment. The model uses the real experimental data as boundary conditions and initial settings, which also considers the CO2 solubility in liquid piston (mineral oil) for designing, optimizing and customizing the compression chambers. The validation was carried out with experimental data using a prototype with 3.5 kW capacity. The results have demonstrated the accuracy of the dynamic model (6.2% relative error for chamber pressure and 0.5 K deviation for chamber temperature), which provide a guideline for designing and customizing the isothermal compression cycle.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution and seismology of Procyon

We have calculated an array of stellar models for Procyon A that are based on the best physics available to us, including the latest opacities and nuclear cross sections. The array of models spans the error space centered on Procyon's mass, chemical composition, effective temperature, and luminosity. We find that with OPAL opacities no convective overshoot is needed at the edge of the convective core to match Procyon's position in the H-R diagram. We discuss the discrepancy between the astrophysical mass and the astrometric mass, which is reduced by the more up-to-date physics in the models, and describe how oscillation data can be used to distinguish among the possible models and help resolve the discrepancy. We have calculated oscillation frequencies for the l = 0, 1, 2, and 3 p-modes and oscillation periods for the l = 1, 2, and 3 g-modes for several of these models.

Guenther, D. B.↗

Accelerating multilevel Markov Chain Monte Carlo using machine learning models

Here, this work presents an efficient approach for accelerating multilevel Markov Chain Monte Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning models. While conventional techniques for large-scale Bayesian inference often substitute computationally expensive high-fidelity models with machine learning models, thereby introducing approximation errors, our approach offers a computationally efficient alternative by augmenting high-fidelity models with low-fidelity ones within a hierarchical framework. The multilevel approach utilizes the low-fidelity machine learning model (MLM) for inexpensive evaluation of proposed samples thereby improving the acceptance of samples by the high-fidelity model. The hierarchy in our multilevel algorithm is derived from geometric multigrid hierarchy. We utilize an MLM to accelerate the coarse level sampling. Training machine learning model for the coarsest level significantly reduces the computational cost associated with generating training data and training the model. We present an MCMC algorithm to accelerate the coarsest level sampling using MLM and account for the approximation error introduced. We provide theoretical proofs of detailed balance and demonstrate that our multilevel approach constitutes a consistent MCMC algorithm. Additionally, we derive the expression for cost reduction due to machine learning model to facilitate cost analysis of the hierarchical sampling algorithm. Our technique is demonstrated on a standard benchmark inference problem in groundwater flow, where we estimate the probability density of a quantity of interest using a four-level MCMC algorithm. Our proposed algorithm accelerates multilevel sampling by a factor of two while achieving similar accuracy compared to sampling using the standard multilevel algorithm.

97 MATHEMATICS AND COMPUTING↗

TPSAS-NF1676L-19238-DND

The 24-hour cycle of solar insolation drives diurnal cycles in Earth system processes critical to climate: radiation, convection, turbulence, and cloud processes. Cloud properties, temperature, radiative fluxes, and precipitation exhibit robust diurnal cycles as a result. The presence of these robust diurnal cycles fundamentally changes the time mean and variability of the TOA energy budget. Therefore, there is a need to understand the diurnal cycle contributions to the TOA energy budget in both observations and climate models. The first objective of this study is to quantify the diurnal cycle impact on the time mean and variability of the TOA energy budget in observations defined as the difference between TOA fluxes computed with diurnally uniform and diurnally varying cloud properties. The required observational input is observed from the combined Terra+Aqua Cloud and Earth's Radiant Energy System (CERES) data products ranging from June 2002 through October 2012. The second objective of this study is to evaluate the impact of known errors in the model diurnal cycle simulation on the time mean and variability of the TOA energy budget. The results indicate that (1) the diurnal cycle impacts on the regional time mean Tropical longwave and shortwave fluxes range from 1-3 W m 2 and 5 25 W m 2, respectively; (2) the diurnal cycle contributions to TOA flux variability exceeds 50% in land convective regions (e.g., central South America and central Africa);(3) the diurnal cycle contributions to the TOA shortwave and longwave flux time mean and variability in the CanAM4 GCM are too small in land convective regions due to errors in the diurnal distribution of convective clouds grid box errors exceed 100%.

Patrick C Taylor↗

A three-parameter approximation for radiative transfer in nonhomogeneous atmospheres - Application to the O3 9.6-micron band

A three-parameter scaling approximation is developed in conjunction with random models for the computation of transmittance along nonhomogeneous atmospheres. In addition to the scaled path length and pressure, a third parameter is introduced to account for intermediate absorption. With reference to the results derived from line-by-line calculations, the accuracy of this scaling approximation is assessed in different atmospheric conditions. For the O3 9.6-micron band where the Curtis-Godson and other conventional scaling methods usually produce large errors, the present model can achieve relative accuracies within about 1 percent in transmittance and the downward surface flux. Accuracy of about 0.15 K/d can be obtained in cooling rate calculations. The present three-parameter approximation is well suited for radiative transfer parameterizations in which the demand for high accuracy is required in the transmittance and flux calculations.

Fu, Qiang↗

Massively Parallel Assimilation of TOGA/TAO and Topex/Poseidon Measurements into a Quasi Isopycnal Ocean General Circulation Model Using an Ensemble Kalman Filter

A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Keppenne, Christian L.↗