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,279 records · Page 71

Uncertainty Quantification of CFD Data Generated for a Model Scramjet Isolator Flowfield

Computational fluid dynamics is now considered to be an indispensable tool for the design and development of scramjet engine components. Unfortunately, the quantification of uncertainties is rarely addressed with anything other than sensitivity studies, so the degree of confidence associated with the numerical results remains exclusively with the subject matter expert that generated them. This practice must be replaced with a formal uncertainty quantification process for computational fluid dynamics to play an expanded role in the system design, development, and flight certification process. Given the limitations of current hypersonic ground test facilities, this expanded role is believed to be a requirement by some in the hypersonics community if scramjet engines are to be given serious consideration as a viable propulsion system. The present effort describes a simple, relatively low cost, nonintrusive approach to uncertainty quantification that includes the basic ingredients required to handle both aleatoric (random) and epistemic (lack of knowledge) sources of uncertainty. The nonintrusive nature of the approach allows the computational fluid dynamicist to perform the uncertainty quantification with the flow solver treated as a "black box". Moreover, a large fraction of the process can be automated, allowing the uncertainty assessment to be readily adapted into the engineering design and development workflow. In the present work, the approach is applied to a model scramjet isolator problem where the desire is to validate turbulence closure models in the presence of uncertainty. In this context, the relevant uncertainty sources are determined and accounted for to allow the analyst to delineate turbulence model-form errors from other sources of uncertainty associated with the simulation of the facility flow.

Baurle, R. A.↗

A Multi-Satellite Assimilation and Modeling Platform to Construct a Global and Complete view of the Hydrologic Cycle

What captivates me about the research field of satellite remote sensing of the earth system science is the opportunity to help human beings by addressing challenging questions related to the current and future availability of natural resources for domestic, agricultural, and industrial needs. For example, a continuously varying climate poses a threat to the water resource, not only in developing countries, but also in economically well-developed regions, creating the urgent need to characterize and predict its spatial and temporal availability. Specifically, my curiosity is driven by the challenge of merging cutting-edge space technology and earth observations (i.e., remote sensing) with state-of-the-art models for the purpose of improving our scientific knowledge about the variability and the change of the hydrologic cycle.The only practical way to observe the land surface processes (e.g., hydrologic cycle) on continental to global scales is via satellite remote sensing. Though remote sensing can make spatially comprehensive measurements of various components of the land surface system, it cannot provide direct information on the entire system, and the measurements represent only a snapshot in time. Physical based models may be used to continuously predict the temporal and spatial earth processes, but these predictions are often poor, due to model initialization, parameter and meteorological forcing errors, and inadequate model physics and/or resolution. Thus, an attractive prospect is to combine the strengths of land surface models and observations (and minimize the weaknesses of both) to provide a superior land surface state estimate. This is the goal of land surface data assimilation. Data assimilation provides a better estimate of the environmental states than either models or observations could individually do. The broad concept of data assimilation can be applied to various disciplines such as hydrology, ecology, environmental hazards, agriculture and economy.My research is targeted to integrating multiple satellite information having multi-sensor and multi-resolution assimilation within today's state-of-the-art hydrologic models. Multi-sensor and multi-resolution assimilation techniques represent a necessary milestone in future earth science applications because they offer the capability of comprehensively integrating disparate types of earth observations via deep-learning techniques. A multi-satellite assimilation and modeling platform will ultimately provide a robust and complete dataset to more fully understand earth's dynamics. In fact, the multi-satellite assimilation platforms build a comprehensive description of the earth processes that is geared toward accurately representing the complexity of natural and anthropogenic interactions in land surface processes. Most importantly, these platforms are suitable as decision support tools for applications across different aspects of earth system science.My recent work has focused on recent multi-sensor data assimilation techniques targeted at improving snow, soil moisture, groundwater, and terrestrial water storage hydrological states.

Girotto, Manuela↗

Modeling and Simulation Techniques for the NASA SLS Service Module Panel Separation Event; from Loosely-Coupled Euler to Fully-Coupled 6-DOF, Time-Accurate, Navier-Stokes Methodologies

An aerodynamic database has been generated for use by the Orion Multi-Purpose Crew Vehicle (MPCV) Program to analyze Service Module (SM) panel jettison from the NASA SLS vehicle. The database is a combination of CFD data for the panel aerodynamic coefficients, and MATLAB code written to query the CFD data. The Cart3D inviscid CFD flow solver was used to generate the panel aerodynamic coefficients for static panel orientations and free stream conditions that can occur during the jettison event. The MATLAB code performs the multivariate interpolation to obtain aerodynamic coefficients. The MATLAB code uses input for SM panel parameters and returns the SM panel aerodynamic force and moment coefficients for use with a Six-Degree-of-Freedom (6-DOF) motion solver to model the jettison event. This paper examines the accuracy of the sequential-static database approach by modeling the panel jettison event with a fully-coupled, time-dependent, viscous, moving-body CFD simulation. The fully-coupled simulation is obtained using the Loci/Chem unstructured Navier-Stokes CFD solver. The results show that the fully-coupled approach agrees well with the loosely-coupled database/6-DOF approach, indicating that unsteady effects are minimal for the panel jettison event. These results suggest that the database/6-DOF approach is sufficient. In addition, this paper presents the development of an uncertainty model for use in Monte Carlo analysis of the panel jettison event. Here viscous CFD simulations are obtained with Loci/Chem and compared to the inviscid CFD forces and moments. An uncertainty model based on model-form error and numerical error is presented.

Hall, Leslie H.↗

Progress in Atmospheric Carbon Monitoring Using NASA's GEOS Model and Data from the OCO and GOSAT Missions

NASA's Global Modeling and Assimilation Office (GMAO) produces a variety of carbon products based the synthesis of satellite remote sensing data and outputs of the Goddard Earth Observing System (GEOS). This includes bottom-up surface fluxes due to fossil fuel emissions, biomass burning, terrestrial biospheric exchange, and ocean exchangeconstrained by measurements of nighttime lights, fire radiative power, normalized difference vegetation index, and ocean color. These fluxes are the basis of top-down estimates of carbon concentrations and fluxes. In particular, the GMAO system processes retrievals of column carbon dioxide (XCO2) from GOSAT and OCO-2 to produce a high-resolution, long-term global analysis of CO2 in three dimensions every 6 hours. Here, we discuss the potential applications of such products for satellite intercomparison and evaluation against independent, non-coincident data. We also highlight the ability to provide monthly global atmospheric growth rates inferred from the assimilated CO2 concentration product. Finally, we discuss the challenges facing such products including bias correction and the estimation and analysis of model transport errors.

Weir, B.↗

Orbit Determination Covariance Analyses for the Parker Solar Probe Mission

This paper details pre-launch navigation covariance analyses for the Parker Solar Probe mission. Baseline models and error assumptions are outlined. The results demonstrate how navigation will satisfy requirements and are used to define operational plans. A few sensitivities are identified and the accompanying investigations are described. Predicted state uncertainty results show that most requirements are met with substantial margin. Moreover, navigation sensitivities may be accommodated operationally and this has been incorporated into project planning. Detailed results are presented only for select launch dates, however twenty unique trajectories (one per launch opportunity) have been assessed.

Jones, Drew Ryan↗

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗

Application of Compressive Sensing for Gravitational Microlensing Events

Compressive sensing is a novel mathematical technique for simultaneous data acquisition and compression. This technique is particularly apt for time-series photometric measurements. In this work we apply compressive sensing to time-series photometric measurements specifically obtained due to gravitational microlensing events. We show through simulation modelling the error sensitivity for detecting microlensing event parameters. Particularly, we show the relation of the amount of error and its impact on the microlensing parameters of interest. We derive statistical error bounds to apply those as a baseline for analyzing the effectiveness of compressive sensing application. Our results conclude that for single and binary microlensing events we can obtain error less than 1% over a 3-pixel radius of the center of the microlensing star by using 25% Nyquist rate measurements.

Compressive Sensing↗

CropManage Application for Vineyard Irrigation Decision-Support

CropManage is a free web-application developed by U.C. Cooperative Extension to support evapotranspiration based irrigation scheduling and nutrient management for major specialty crops. Prescribed phenology curves are used to develop daily estimates of canopy cover within a given field, based on days since planting (annual crops) or budbreak (trees, vines). These curves are modulated by a MaxCan parameter representing seasonal maximum canopy cover. Crop development observations can be used to adjust for such factors as weather anomalies or non-standard agronomic practice, as needed. Canopy cover is converted to crop coefficient and combined with reference evapotranspiration to derive daily water consumption. Guidance on crop water requirement is then conveyed to users in terms of system runtime issued on-demand for a given date, largely based on total evapotranspiration since last irrigation event. In this study, CropManage was adapted to vineyards by adding modules accounting for early-season soil moisture depletion and cover crop presence. A crop stress parameter was added to accommodate deficit irrigation practice, allowing the user to specify percentage departure from full water requirement along with start/stop dates. An initial verification exercise was performed on three winegrape vineyards located in California’s Central Coast (2020), North Coast (2020) and Central Valley (2019). Daily crop evapotranspiration was monitored by eddy-covariance fluxtowers. MaxCan was measured by ground and satellite observation. Stress regime was specified by grower practice where available, otherwise stress levels were inferred from applied water records. Mean absolute error and mean bias error of modeled cumulative evapotranspiration were computed with respect to the eddy covariance measurements collected throughout the growing season. Results indicate the modified CropManage water management module performs reasonably well for winegrape. Additional effort is planned to modify the nutrient module for vineyard use.

CropManage↗

Remote Radiation Sensing Using Aerial and Ground Platforms

Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗

Mathematical Models for Doppler Measurements

Error analysis increases precision of navigation. Report presents improved mathematical models of analysis of Doppler measurements and measurement errors of spacecraft navigation. To take advantage of potential navigational accuracy of Doppler measurements, precise equations relate measured cycle count to position and velocity. Drifts and random variations in transmitter and receiver oscillator frequencies taken into account. Mathematical models also adapted to aircraft navigation, radar, sonar, lidar, and interferometry.

Lear, William M.↗

Prediction of pilot reserve attention capacity during air-to-air target tracking

Reserve attention capacity of a pilot was calculated using a pilot model that allocates exclusive model attention according to the ranking of task urgency functions whose variables are tracking error and error rate. The modeled task consisted of tracking a maneuvering target aircraft both vertically and horizontally, and when possible, performing a diverting side task which was simulated by the precise positioning of an electrical stylus and modeled as a task of constant urgency in the attention allocation algorithm. The urgency of the single loop vertical task is simply the magnitude of the vertical tracking error, while the multiloop horizontal task requires a nonlinear urgency measure of error and error rate terms. Comparison of model results with flight simulation data verified the computed model statistics of tracking error of both axes, lateral and longitudinal stick amplitude and rate, and side task episodes. Full data for the simulation tracking statistics as well as the explicit equations and structure of the urgency function multiaxis pilot model are presented.

Onstott, E. D.↗

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗

Background Error Covariance Estimation using Information from a Single Model Trajectory with Application to Ocean Data Assimilation into the GEOS-5 Coupled Model

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory. SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Data Assimilation↗

Spin glasses and error-correcting codes

In this article, we study a model for error-correcting codes that comes from spin glass theory and leads to both new codes and a new decoding technique. Using the theory of spin glasses, it has been proven that a simple construction yields a family of binary codes whose performance asymptotically approaches the Shannon bound for the Gaussian channel. The limit is approached as the number of information bits per codeword approaches infinity while the rate of the code approaches zero. Thus, the codes rapidly become impractical. We present simulation results that show the performance of a few manageable examples of these codes. In the correspondence that exists between spin glasses and error-correcting codes, the concept of a thermal average leads to a method of decoding that differs from the standard method of finding the most likely information sequence for a given received codeword. Whereas the standard method corresponds to calculating the thermal average at temperature zero, calculating the thermal average at a certain optimum temperature results instead in the sequence of most likely information bits. Since linear block codes and convolutional codes can be viewed as examples of spin glasses, this new decoding method can be used to decode these codes in a way that minimizes the bit error rate instead of the codeword error rate. We present simulation results that show a small improvement in bit error rate by using the thermal average technique.

Belongie, M. L.↗

Introduction to CAUSES: Description of Weather and Climate Models and Their Near‐Surface Temperature Errors in 5 Day Hindcasts near the Southern Great Plains

We introduce the Clouds Above the United States and Errors at the Surface (CAUSES) project with its aim of better understanding the physical processes leading to warm screen temperature biases over the American Midwest in many numerical models. In this first of four companion papers, 11 different models, from nine institutes, perform a series of 5 day hindcasts, each initialized from reanalyses. After describing the common experimental protocol and detailing each model configuration, a gridded temperature data set is derived from observations and used to show that all the models have a warm bias over parts of the Midwest. Additionally, a strong diurnal cycle in the screen temperature bias is found in most models. In some models the bias is largest around midday, while in others it is largest during the night. At the Department of Energy Atmospheric Radiation Measurement Southern Great Plains (SGP) site, the model biases are shown to extend several kilometers into the atmosphere. Finally, to provide context for the companion papers, in which observations from the SGP site are used to evaluate the different processes contributing to errors there, it is shown that there are numerous locations across the Midwest where the diurnal cycle of the error is highly correlated with the diurnal cycle of the error at SGP. This suggests that conclusions drawn from detailed evaluation of models using instruments located at SGP will be representative of errors that are prevalent over a larger spatial scale.

Morcrette, C. J.↗