Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical errors”

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,153 records · Page 64

Numerical Studies of the On-Line Reactivity Monitoring Techniques in Subcritical Systems: Revisit of the Area-Ratio Method, the Neutron Source-Modulation Method, and the Flux-to-Current Ratio Method

In the subcritical systems, the subcriticality level must be monitored continuously to ensure the safe operation of the systems. Techniques such as the neutron source-modulation method and the flux-to-current ratio method have been proposed for on-line reactivity measurements of the subcritical systems. However, the accuracy of these methods is unknown in the actual applications. In this report, the fundamental concepts of these techniques are visited, with numerical simulation techniques developed to validate these on-line reactivity monitoring methodologies. Specifically, the Monte Carlo numerical simulation results show that the neutron source-modulation method is not accurate in determining the reactivity of the subcritical systems. The flux-to-current ratio method can provide on-line reactivity measurements, and its spatial effects have been studied and well understood through the numerical simulations of the fuel loading of the neutron source facility of Ukraine. In addition, a modified version of the area-ratio method is proposed for the on-line measurements of the absolute reactivity of the subcritical systems. Monte Carlo numerical simulations showed that this method can provide the absolute reactivity of the subcritical systems within ~5% error, if the neutron detectors are placed close to the fuel assembly interface with the reflector zone or in the reflector zones away from the neutron source.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Error Correction for Tandem Data-Transmission Paths

Mathematical analysis for digital data transmission calculates optimum number of binary error-correcting repeaters to install in given number of wideband channel links. Asymptotic results compared to computed numerical results.

Posner, E. C.↗

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗

Three-dimensional Navier-Stokes calculations of multiple interacting vortex rings

Results from a finite-difference Navier-Stokes code for three-dimensional, unsteady, vortical flows in unbounded domains are presented and analyzed in this paper. The vortical flows presented are representative of vortex rings and other closed vortical tubes or structures in fluid mechanics. Such structures are important elements in fluid flows such as jets, atmospheric turbulence, and the far-field wakes of aircraft, and studies of their interaction may aid in an understanding of complex fluid flows. The paper demonstrates that computational methods can be used as a viable alternative or supplement to experimental techniques for studying the physics of vortex flows. The separate visualization of vortex stretching, convection, and diffusion is presented in this paper for a single elliptical vortex ring.The calculations employ a truncated series expansion technique to simulate the unbounded nature of the fluid flow with a finite computational domain, which is a more accurate technique than the conventional freestream boundary specification. The numerical divergence of the three-dimensional vorticity field is considered as a useful estimate of truncation error, and the use of a kinetic energy decay law as a calculation check is demonstrated. Results from the Navier-Stokes code are presented for the unsteady motion of two and four vortex rings along parallel axes, and the results agree qualitatively with experimental flow visualization.

Chamberlain, J. P.↗

Unsteady lifting-line theory as a singular-perturbation problem

Unsteady lifting-line theory is developed for a wing of large aspect ratio oscillating at low frequency in inviscid incompressible flow. The wing is assumed to have a rigid chord but a flexible span. Use of the method of matched asymptotic expansions reduces the problem from a singular integral equation to quadrature. The pressure field and airloads, for a prescribed wing shape and motion, are obtained in closed form as expansions in inverse aspect ratio. A rigorous definition of unsteady induced downwash is also obtained. Numerical calculations are presented for an elliptic wing in pitch and heave; compared with numerical lifting-surface theory, computation time is reduced significantly. The present work also identifies and resolves errors in the unsteady lifting line theory of James (1975), and points out a limitation in that of Van Holten (1975, 1976, 1977).

Ahmadi, A.↗

Mean-square error bounds for reduced-error linear state estimators

The mean-square error of reduced-order linear state estimators for continuous-time linear systems is investigated. Lower and upper bounds on the minimal mean-square error are presented. The bounds are readily computable at each time-point and at steady state from the solutions to the Ricatti and the Lyapunov equations. The usefulness of the error bounds for the analysis and design of reduced-order estimators is illustrated by a practical numerical example.

Baram, Y.↗

Mean-square-error bounds for reduced-order linear state estimators

The mean-square error of reduced-order linear state estimators for continuous-time linear systems is investigated. Lower and upper bounds on the minimal mean-square error are presented. The bounds are readily computable at each time-point and at steady state from the solutions to the Ricatti and the Liapunov equations. The usefulness of the error bounds for the analysis and design of reduced-order estimators is illustrated by a practical numerical example.

Baram, Y.↗

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗

High-resolution processing of long-pulse-lidar data

The purpose of this work is to demonstrate the performance of the Fourier-deconvolution technique developed to improve the resolution of the long-pulse coherent lidar power profile, and to take into account the multiplicative fluctuations. The behavior of the error due to the power fluctuations and a method to reduce it are analyzed theoretically and simulated numerically. A processing of real data obtained from the National Oceanic and Atmospheric Administration (NOAA) ground-based Doppler lidar is also presented. Similar problems were investigated using a different approach based on introducing a correction function to the lidar equation.

Gurdev, L. L.↗

Practical aspects of spatially high accurate methods

The computational qualities of high order spatially accurate methods for the finite volume solution of the Euler equations are presented. Two dimensional essentially non-oscillatory (ENO), k-exact, and 'dimension by dimension' ENO reconstruction operators are discussed and compared in terms of reconstruction and solution accuracy, computational cost and oscillatory behavior in supersonic flows with shocks. Inherent steady state convergence difficulties are demonstrated for adaptive stencil algorithms. An exact solution to the heat equation is used to determine reconstruction error, and the computational intensity is reflected in operation counts. Standard MUSCL differencing is included for comparison. Numerical experiments presented include the Ringleb flow for numerical accuracy and a shock reflection problem. A vortex-shock interaction demonstrates the ability of the ENO scheme to excel in simulating unsteady high-frequency flow physics.

Godfrey, Andrew G.↗

Damage detection in elastic structures using vibratory residual forces and weighted sensitivity

A methodology is presented for detecting structural damage in elastic structures by nondestructive means. Measured modal test data along with a correlated analytical structural model are used to locate potentially damaged regions using residual modal force vectors and to conduct a weighted sensitivity analysis to assess the extent of mass and/or stiffness variations, where damage is characterized as a stiffness reduction. The current approach is unique among other approaches in that it accounts for (1) variations in system mass, system stiffness, and mass center (locations), (2) perturbations of both the natural frequencies and modal vectors, and (3) statistical confidence factors for the structural parameters and potential experimental instrumentation error. Moreover, this procedure can be used with either full or reduced models. A wide variety of numerical examples are presented that show that the current method provides a precise indication of both the location and the extent of structural damage.

Ricles, J. M.↗

A study of the torque equilibrium of an autogiro rotor

Two improvements have been made in the method developed in NACA Reports nos. 487 and 591 for the estimation of the inflow velocity required to overcome a given decelerating torque in an autogiro rotor. At low tip-speed ratios, where the assumptions necessary for the analytical integrations of the earlier papers are valid, the expressions therein derived are greatly simplified by combining and eliminating terms with a view of minimizing the numerical computations required. At high tip-speed ratios, by means of charts based on graphical integrations, errors inherent in the assumptions associated with the analytical method are largely eliminated. The suggested method of estimating the inflow velocity presupposes a knowledge of the decelerating torque acting on the rotor; all available full-scale experimental information on this subject is included.

Bailey, F J , Jr↗

Multigrid Approach to Incompressible Viscous Cavity Flows

Two-dimensional incompressible viscous driven-cavity flows are computed for Reynolds numbers on the range 100-20,000 using a loosely coupled, implicit, second-order centrally-different scheme. Mesh sequencing and three-level V-cycle multigrid error smoothing are incorporated into the symmetric Gauss-Seidel time-integration algorithm. Parametrics on the numerical parameters are performed, achieving reductions in solution times by more than 60 percent with the full multigrid approach. Details of the circulation patterns are investigated in cavities of 2-to-1, 1-to-1, and 1-to-2 depth to width ratios.

Wood, William A.↗

Mesoscale Assimilation of TMI Rainfall Data with 4DVAR: Sensitivity Studies

Sensitivity studies are performed on the assimilation of TRMM (Tropical Rainfall Measurement Mission) Microwave Imager (TMI) derived rainfall data into a mesoscale model using a four-dimensional variational data assimilation (4DVAR) technique. A series of numerical experiments is conducted to evaluate the impact of TMI rainfall data on the numerical simulation of Hurricane Bonnie (1998). The results indicate that rainfall data assimilation is sensitive to the error characteristics of the data and the inclusion of physics in the adjoint and forward models. In addition, assimilating the rainfall data alone is helpful for producing a more realistic eye and rain bands in the hurricane but does not ensure improvements in hurricane intensity forecasts. Further study indicated that it is necessary to incorporate TMI rainfall data together with other types of data such as wind data into the model, in which case the inclusion of the rainfall data further improves the intensity forecast of the hurricane. This implies that proper constraints may be needed for rainfall assimilation.

Tao, Wei-Kuo↗

PlanWorks: A Debugging Environment for Constraint Based Planning Systems

Numerous planning and scheduling systems employ underlying constraint reasoning systems. Debugging such systems involves the search for errors in model rules, constraint reasoning algorithms, search heuristics, and the problem instance (initial state and goals). In order to effectively find such problems, users must see why each state or action is in a plan by tracking causal chains back to part of the initial problem instance. They must be able to visualize complex relationships among many different entities and distinguish between those entities easily. For example, a variable can be in the scope of several constraints, as well as part of a state or activity in a plan; the activity can arise as a consequence of another activity and a model rule. Finally, they must be able to track each logical inference made during planning. We have developed PlanWorks, a comprehensive system for debugging constraint-based planning and scheduling systems. PlanWorks assumes a strong transaction model of the entire planning process, including adding and removing parts of the constraint network, variable assignment, and constraint propagation. A planner logs all transactions to a relational database that is tailored to support queries for of specialized views to display different forms of data (e.g. constraints, activities, resources, and causal links). PlanWorks was specifically developed for the Extensible Universal Remote Operations Planning Architecture (EUROPA(sub 2)) developed at NASA, but the underlying principles behind PlanWorks make it useful for many constraint-based planning systems. The paper is organized as follows. We first describe some fundamentals of EUROPA(sub 2). We then describe PlanWorks' principal components. We then discuss each component in detail, and then describe inter-component navigation features. We close with a discussion of how PlanWorks is used to find model flaws.

Daley, Patrick↗

Design of Multifunctional Hierarchical Space Structures

We describe a system for the design of space structures with tunable structural properties based on the discrete assembly of modular lattice elements. These lattice elements can be constructed into larger beam-like elements, which can then be assembled into large scale truss structures. These discrete lattice elements are reversibly assembled with mechanical fasteners, which allows them to be arbitrarily reconfigured into various application-specific designs. In order to assess the validity of this approach, we design two space structures with similar geometry but widely different structural requirements: an aerobrake, driven by strength requirements, and a precision segmented reflector, driven by stiffness and accuracy requirements. We will show agreement between simplified numerical models based on hierarchical assembly and analytical solutions. We will also present an assessment of the error budget resulting from the assembly of discrete structures. Lastly, we will address launch vehicle packing efficiency issues for transporting these structures to lower earth orbit.

Aeronautics↗

Integrated Multi-Satellite Evaluation for the Global Precipitation Measurement: Impact of Precipitation Types on Spaceborne Precipitation Estimation

Integrated multi-sensor assessment is proposed as a novel approach to advance satellite precipitation validation in order to provide users and algorithm developers with an assessment adequately coping with the varying performances of merged satellite precipitation estimates. Gridded precipitation rates retrieved from space sensors with quasi-global coverage feed numerous applications ranging from water budget studies to forecasting natural hazards caused by extreme events. Characterizing the error structure of satellite precipitation products is recognized as a major issue for the usefulness of these estimates. The Global Precipitation Measurement (GPM) mission aims at unifying precipitation measurements from a constellation of low-earth orbiting (LEO) sensors with various capabilities to detect, classify and quantify precipitation. They are used in combination with geostationary observations to provide gridded precipitation accumulations. The GPM Core Observatory satellite serves as a calibration reference for consistent precipitation retrieval algorithms across the constellation. The propagation of QPE uncertainty from LEO active/passive microwave (PMW) precipitation estimates to gridded QPE is addressed in this study, by focusing on the impact of precipitation typology on QPE from the Level-2 GPM Core Observatory Dual-frequency Precipitation Radar (DPR) to the Microwave Imager (GMI) to Level-3 IMERG precipitation over the Conterminous U.S. A high-resolution surface precipitation used as a consistent reference across scales is derived from the ground radar-based Multi-Radar/Multi-Sensor. While the error structure of the DPR, GMI and subsequent IMERG is complex because of the interaction of various error factors, systematic biases related to precipitation typology are consistently quantified across products. These biases display similar features across Level-2 and Level-3, highlighting the need to better resolve precipitation typology from space and the room for improvement in global-scale precipitation estimates. The integrated analysis and framework proposed herein applies more generally to precipitation estimates from sensors and error sources affecting low-earth orbiting satellites and derived gridded products.

Rainfall↗

MRAC with Time-Varying Reference Model from Convex Combination

An extension of model reference adaptive control is presented that accommodates use of a time-varying reference model. Specifically, the reference model is taken to be a time-varying convex combination of two linear, time-invariant models. The design is intended to act as a way to smoothly transition between two different reference models without resorting to a scheduled switch. It also provides the ability to use an interpolated reference model when the plant is operating between design points. The time variation of the combination must satisfy some requirements to ensure stability but is otherwise user choice. Subject to these requirements, bounded tracking error behavior is demonstrated via Lyapunov stability analysis for the single-input, single-output, output feedback case. Tracking error convergence is asymptotic when time variation ceases. The proposed design is demonstrated in simulation of a numerical model.

Hashemi, Kelley E.↗