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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗
NASA NTRS1981

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗
NASA NTRS1985

Proximity Indicator for Remote Manipulator

Display indicates pitch, yaw, and distance of remote manipulator with respect to object to be grasped. Displays numerical values and bargraph simulation of position and orientation of hand. When errors in position and orientation are small enough to ensure successful grasp, unit alerts operator with audible and visible signals.

Bejczy, A. K.↗
NASA NTRS1987

Experiments on thermoacoustic convection heat transfer in gravity and zero-gravity environments

The results of an experimental study of thermoacoustic convection (TAC) heat transfer in gravity and zero-gravity environments are presented. The experimental apparatus consisted of a cylinder containing air as the compressible fluid. The enclosed air was heated electrically at the top surface which consisted of a thin high-resistance steel foil connected to a power source. Thermocouples were used to measure the transient temperature of the air on the axis of the cylinder and the heated surface in the both zero-gravity and gravity environments. The zero-gravity tests were performed in the Zero-Gravity Drop Tower Facility of NASA-Lewis Research Center. The experimental results were corrected for the error due to radiation absorption by the thermocouples. A conduction-only numerical heat transfer model was developed to compute the transient air temperature in the cylindrical geometry. The results were compared to the experimental data to determine the significance of the thermoacoustic convection heat transfer mechanism. It is observed that the rate of heat transfer to the air measured during the experiments is consistently higher than that obtained by the conduction-only solution indicating a significant presence of the TAC heat transfer. Further experiments are planned to measure directly (1) the radiative heat transfer contribution to the rise in the air temperature, and (2) the air pressure oscillations within the cylinder that are responsible for the convective heat transfer mode.

Parang, Masood↗
NASA NTRS1988

A multivariate variational objective analysis-assimilation method. Part 2: Case study results with and without satellite data

The variational multivariate assimilation method described in a companion paper by Achtemeier and Ochs is applied to conventional and conventional plus satellite data. Ground-based and space-based meteorological data are weighted according to the respective measurement errors and blended into a data set that is a solution of numerical forms of the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation for a dry atmosphere. The analyses serve first, to evaluate the accuracy of the model, and second to contrast the analyses with and without satellite data. Evaluation criteria measure the extent to which: (1) the assimilated fields satisfy the dynamical constraints, (2) the assimilated fields depart from the observations, and (3) the assimilated fields are judged to be realistic through pattern analysis. The last criterion requires that the signs, magnitudes, and patterns of the hypersensitive vertical velocity and local tendencies of the horizontal velocity components be physically consistent with respect to the larger scale weather systems.

Achtemeier, Gary L.↗
NASA NTRS1989

The preliminary SOL (Sizing and Optimization Language) reference manual

The Sizing and Optimization Language, SOL, a high-level special-purpose computer language has been developed to expedite application of numerical optimization to design problems and to make the process less error-prone. This document is a reference manual for those wishing to write SOL programs. SOL is presently available for DEC VAX/VMS systems. A SOL package is available which includes the SOL compiler and runtime library routines. An overview of SOL appears in NASA TM 100565.

Lucas, Stephen H.↗
NASA NTRS1988

Development of a satellite SAR image spectra and altimeter wave height data assimilation system for ERS-1

The applicability of ERS-1 wind and wave data for wave models was studied using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR data. A series of global wave hindcasts is made for the surface stress and surface wind fields by assimilation of scatterometer data for the full 96-day SEASAT and also for two wind field analyses for shorter periods by assimilation with the higher resolution ECMWF T63 model and by subjective analysis methods. It is found that wave models respond very sensitively to inconsistencies in wind field analyses and therefore provide a valuable data validation tool. Comparisons between SEASAT SAR image spectra and theoretical SAR spectra derived from the hindcast wave spectra by Monte Carlo simulations yield good overall agreement for 32 cases representing a wide variety of wave conditions. It is concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra with confidence for wave studies if supported by realistic wave models and theoretical computations of the strongly nonlinear mapping of the wave spectrum into the SAR image spectrum. A closed nonlinear integral expression for this spectral mapping relation is derived which avoids the inherent statistical errors of Monte Carlo computations and may prove to be more efficient numerically.

Hasselmann, Klaus↗
NASA NTRS1992

Spurious Solutions Of Nonlinear Differential Equations

Report utilizes nonlinear-dynamics approach to investigate possible sources of errors and slow convergence and non-convergence of steady-state numerical solutions when using time-dependent approach for problems containing nonlinear source terms. Emphasizes implications for development of algorithms in CFD and computational sciences in general. Main fundamental conclusion of study is that qualitative features of nonlinear differential equations cannot be adequately represented by finite-difference method and vice versa.

Yee, H. C.↗
NASA NTRS1992

Adaptive Grids For 3-D Parabolized Navier-Stokes Computations

Computational grids adjusted iteratively to reduce errors. Procedure for iterative adjustment of computational grids developed for use in numerical solution of three-dimensional parabolized Navier-Strokes equations of flow. Refines grid in regions of high gradients in initial computed flow so on subsequent iterations, one can satisfy competing requirements to produce solution capturing pertinent features as shock waves and boundary layers, accurate resolution of which typically requires fine meshes; and prevent unnecessary refinement of mesh elsewhere, preventing undue increase in amount of computation.

Harvey, A. D.↗
NASA NTRS1991

Finite element structural redesign by large admissible perturbations

In structural redesign, two structural states are involved; the baseline (known) State S1 with unacceptable performance, and the objective (unknown) State S2 with given performance specifications. The difference between the two states in performance and design variables may be as high as 100 percent or more depending on the scale of the structure. A Perturbation Approach to Redesign (PAR) is presented to relate any two structural states S1 and S2 that are modeled by the same finite element model and represented by different values of the design variables. General perturbation equations are derived expressing implicitly the natural frequencies, dynamic modes, static deflections, static stresses, Euler buckling loads, and buckling modes of the objective S2 in terms of its performance specifications, and S1 data and Finite Element Analysis (FEA) results. Large Admissible Perturbation (LEAP) algorithms are implemented in code RESTRUCT to define the objective S2 incrementally without trial and error by postprocessing FEA results of S1 with no additional FEAs. Systematic numerical applications in redesign of a 10 element 48 degree of freedom (dof) beam, a 104 element 192 dof offshore tower, a 64 element 216 dof plate, and a 144 element 896 dof cylindrical shell show the accuracy, efficiency, and potential of PAR to find an objective state that may differ 100 percent from the baseline design.

Bernitsas, Michael M.↗
NASA NTRS1992

Sampling and Representativeness for a Spaceborne Wind Lidar

The capability of an Earth-orbiting lidar to produce a reliable windfield sampling with an error of 1-3 ms(exp -1) and that is relevant to numerical forecasting and climate studies is discussed. The spatial and temporal resolutions range from 100 km to 1000 km and 3 hr to 1 month respectively. In this respect cloud obstruction is of great concern, as it can prevent the lidar probing whole parts of the atmosphere, for both short and long periods of time. A worldwide analysis of cloudiness from either visual observations made from the ground or else from the Stratospheric Aerosol and Gas Experiment (SAGE) instrument, has shown that the average cloud cover is 60 percent. This analysis is only valid at large scales (time and space) compatible with the currently operating satellites. However, a lidar footprint is only about a hundred meters or less, and it was observed from ground based lidars that, even when the cloudiness is 100 percent, some lidar shots get through. A complete analysis of the probability of probing through clouds by a spaceborne lidar is presently beyond our capacities. Hence, no small scale cloudiness data set is yet available representing all kinds of meteorological situations at all latitudes.

Lieutaud, F.↗
NASA NTRS1992

Model reduction using new optimal Routh approximant technique

An optimal Routh approximant of a single-input single-output dynamic system is a reduced-order transfer function of which the denominator is obtained by the Routh approximation method while the numerator is determined by minimizing a time-response integral-squared-error (ISE) criterion. In this paper, a new elegant approach is presented for obtaining the optimal Routh approximants for linear time-invariant continuous-time systems. The approach is based on the Routh canonical expansion, which is a finite-term orthogonal series of rational basis functions, and minimization of the ISE criterion. A procedure for combining the above approach with the bilinear transformation is also presented in order to obtain the optimal bilinear Routh approximants of linear time-invariant discrete-time systems. The proposed technique is simple in formulation and is amenable to practical implementation.

Hwang, Chyi↗
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 NTRS1993

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗
NASA NTRS1995

Accurate Evaluation of Quantum Integrals

Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schrodinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.

Galant, D. C.↗
NASA NTRS2004

Multipulse PPM on memoryless channels

memoryless channels. First, we derive the maximum likelihood decision rule and an exact expression for the symbol error rate for n (greater-than or equal) 1, avoiding a numerically unstable aspect of = 1 formula of [GK76] and generalizing the n = 2 result of [SV03]. Next, we compare the capacity of multipulse PPM to that of conventional single-pulse PPM when averaged power, peak, and bandwidth constraints are simultaneously imposed. On the basis, we demonstrate that multipulse PPM does not produce appreciable gains over conventional PPM except at high average power.

optical communications↗
NASA NTRS2016

Solar Probe Plus: Unique Navigation Modeling Challenges

The Solar Probe Plus (SPP) mission is preparing to launch in 2018, and will directly investigate the outer atmosphere of our star. At 9.86 solar radii, SPP must operate in an unexplored regime. The environment and aspects of the mission design present some unique challenges for navigation, particularly in terms of modeling the dynamics. Non-gravitational force models, unique to this mission, are given with analytical expressions. For each of these models (and error sources), a maximum bound on the force perturbation magnitude is quantified numerically. Additionally, the effect of charged particles on radiometric observables is discussed, along with methods being employed to pre-process the measurements. This survey is an overview of unique modeling employed by SPP navigation, but also a reference for future missions traveling near the Sun.

Jones, Drew Ryan↗
DOE OSTI2020

Comparison of plug flow and multi-node stratified tank modeling approaches regarding computational efficiency and accuracy

Residential water heaters contain water stratified by temperature-driven density differences. This implies that a water tank can reach a state in which the top and bottom sections have different temperatures, unless mixing happens. A high degree of thermal stratification can improve the efficiency of some water heaters, by saving the amount of energy required for the heat-up process. Studies of stratification became popular in the 1970s and it remains an active research topic today. The research has led to the development of different models and techniques to better predict and define a stratified tanks behavior. By comparing these models and techniques used previously to describe thermal stratification, the phenomenon could be better understood, exploited, and used to increase efficiency and thermal energy capacity in modern water tanks. From the existing models, we found the one-dimensional standard plug-flow and a multi node model to be appropriate for analyzing the processes of the heat up and cool-down in a water tank. These two models are based on energy balances. This work involved comparing the accuracy and computational effort needed to implement these models. To assess accuracy, we compared both types of existing models to experimental data (also collected in this work) which included a heat up process using an external heat pump. This external process included a layering process that has an eddy diffusivity at five times the rate of thermal diffusion. For this project, we implemented the models in MATLAB, the multi-paradigm numerical computing environment. We quantified model accuracy using the root mean squared error between modeled data and experimental data for six measured tank temperatures. Comparing the accuracy and the computational time taken to run the simulation provides a method to contrast the performance of each model and a way to rate it. The multi node model was run using from 6 to 96 spatial nodes; the plug flow model was run using 1 to 0.001 º C temperature bin sizes. Additionally, timesteps were varied from 4 to 236 s. The results quantify the tradeoff between accuracy and computational time, providing guidance for simulations to intelligently select the best model type and simulation parameters. This research can be used to validate the pre-existing models and possibly improve the modern water tank.

Bulnes, Fernando Karg↗