SEARCH · Search NASA

Results for “numerical error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Investigation of ODE integrators using interactive graphics

In the present paper, two FORTRAN programs are described which use an interactive graphic terminal to generate accuracy and stability plots for a given multistep integrator of ordinary differential equations (ODE). One program treats the fixed-stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution. It also generates the linear stability region. The other program generates an analog with an accuracy plot. Both systems are capable of computing method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗

NASA NTRS1980

Errors in finite-difference computations on curvilinear coordinate systems

Curvilinear coordinate systems were used extensively to solve partial differential equations on arbitrary regions. An analysis of truncation error in the computation of derivatives revealed why numerical results may be erroneous. A more accurate method of computing derivatives is presented.

Mastin, C. W.↗

NASA NTRS1979

Multiclass Bayes error estimation by a feature space sampling technique

A general Gaussian M-class N-feature classification problem is defined. An algorithm is developed that requires the class statistics as its only input and computes the minimum probability of error through use of a combined analytical and numerical integration over a sequence simplifying transformations of the feature space. The results are compared with those obtained by conventional techniques applied to a 2-class 4-feature discrimination problem with results previously reported and 4-class 4-feature multispectral scanner Landsat data classified by training and testing of the available data.

Mobasseri, B. G.↗

NASA NTRS1980

Inventory estimation on the massively parallel processor

This paper describes algorithms for efficiently computing inventory estimates from satellite based images. The algorithms incorporate a one dimensional feature extraction which optimizes the pairwise sum of Fisher distances. Biases are eliminated with a premultiplication by the inverse of the analytically derived error matrix. The technique is demonstrated with a numerical example using statistics obtained from an actual Landsat scene. Attention was given to implementation of the Massively Parallel processor (MPP). A timing analysis demonstrates that the inventory estimation can be performed an order of magnitude faster on the MPP than on a conventional serial machine.

Argentiero, P. D.↗

NASA NTRS1981

An optimal numerical filter for wide-field-of-view measurements of earth-emitted radiation

A technique is described in which all data points along an arc of the orbit may be used in an optimal numerical filter for wide-field-of-view measurements of earth emitted radiation. The statistical filter design is derived whereby the filter is required to give a minimum variance estimate of the radiative exitance at discrete points along the ground track of the satellite. An equation for the optimal numerical filter is given by minimizing the estimate error variance equation with respect to the filter weights, resulting in a discrete form of the Wiener-Hopf equation. Finally, variances of the errors in the radiant exitance can be computed along the ground track and in the cross track directions.

Smith, G. L.↗

NASA NTRS1982

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗

NASA NTRS1982

Determination of normal modes from measured complex modes

A technique is presented for computing a set of normal modes from a set of measured complex modes. The number of elements in the modal vectors, which is equal to the number of measurements, can be larger than the number of modes under consideration. It is also shown that the practice of normal mode approximation to complex modes can lead to very large errors when the modes are too complex. A numerical example and a simulated experiment illustrate the concepts discussed and support the theory presented.

Ibrahim, S. R.↗

NASA NTRS1985

Computation of Normal Modes from Identified Complex Modes

A technique is presented to compute a set of normal modes from a set of measured (damped) complex modes. The number of elements in the modal vectors, which is equal to the number of measurements, can be larger than the number of modes under consideration. It is also shown in this paper that the practice of normal mode approximation to complex modes can lead to considerably large errors when the modes are too complex. A numerical example and a simulated experiment are presented to illustrate the concepts discussed and to support the theory presented.

Ibrahim, S. R.↗

NASA NTRS1985

Spectrum Orbit Utilization Program Documentation: SOUP5 Version 3.8 User's Manual, Volume 2, Appendices a Through G

The appendixes of the user manual are presented. Input forms which may be used to prepare data for the SOUP5V3.4 of the R2BCSAT-83 data base are given. The IBM job control language which can be used to run the SOUP5 system from a magnetic tape is described. Copies of a run using the delivered tape and IBM OS/MVS Job Control Language card deck are illustrated. Numerical limits on scenario data requests are listed. Error handling, error messages and editing procedures are also listed. Instructions as to how to enter a protection ratio template are given. And relation between PARC prameter, channelization, channel families, and interference categories are also listed.

Davidson, J.↗

NASA NTRS1985

Grid-free simulation of diffusion using random walk methods

The simulation of the diffusion of a continuum field by the random walk (RW) displacement of a set of particles is considered. Elements of the gradients of the diffusive concentration are transported by computational particles. It is demonstrated that, by the use of concentration gradients in the RW process, statistical errors are reduced and each realization of the numerical solution is a representation of the exact solution. The algorithm is grid-free, and the computational elements move to follow the gradients; hence, the algorithm is self-adaptive, and uniform resolution is achieved for all times.

Ghoniem, A. F.↗

NASA NTRS1986

Sequential ranging integration times in the presence of CW interference in the ranging channel

The Deep Space Network (DSN), managed by the Jet Propulsion Laboratory for NASA, is used primarily for communication with interplanetary spacecraft. The high sensitivity required to achieve planetary communications makes the DSN very susceptible to radio-frequency interference (RFI). In this paper, an analytical model is presented of the performance degradation of the DSN sequential ranging subsystem in the presence of downlink CW interference in the ranging channel. A trade-off between the ranging component integration times and the ranging signal-to-noise ratio to achieve a desired level of range measurement accuracy and the probability of error in the code components is also presented. Numerical results presented illustrate the required trade-offs under various interference conditions.

Mathur, Ashok↗

NASA NTRS1991

A non-linear mathematical model for dynamic analysis of spur gears including shaft and bearing dynamics

A six-degree-of-freedom nonlinear semi-definite model with time varying mesh stiffness has been developed for the dynamic analysis of spur gears. The model includes a spur gear pair, two shafts, two inertias representing load and prime mover, and bearings. As the shaft and bearing dynamics have also been considered in the model, the effect of lateral-torsional vibration coupling on the dynamics of gears can be studied. In the nonlinear model developed several factors such as time varying mesh stiffness and damping, separation of teeth, backlash, single- and double-sided impacts, various gear errors and profile modifications have been considered. The dynamic response to internal excitation has been calculated by using the 'static transmission error method' developed. The software prepared (DYTEM) employs the digital simulation technique for the solution, and is capable of calculating dynamic tooth and mesh forces, dynamic factors for pinion and gear, dynamic transmission error, dynamic bearing forces and torsions of shafts. Numerical examples are given in order to demonstrate the effect of shaft and bearing dynamics on gear dynamics.

Ozguven, H. Nevzat↗

NASA NTRS1992

Implementation of Vigneron's streamwise pressure gradient approximation in the PNS equations

For single sweep parabolized Navier-Stokes solvers, the streamwise pressure gradient must be modified in the subsonic region to eliminate numerical instabilities. The accuracy of this modification on the solution of the parabolized Navier-Stokes equations with Vigneron's technique is shown to depend on how the numerical approximation of the pressure gradient is formed. A simple test case of supersonic laminar flow over a flat plate is computed with two different numerical methods for solving the PNS equations. Significant errors in the temperature profile and skin friction coefficient are demonstrated using a fully conservative differencing treatment of Vigneron's splitting for the pressure gradient typically used in parabolized Navier-Stokes solvers. The physical reason for this error is discussed. An alternate formulation is demonstrated which minimizes these errors.

Morrison, Joseph H.↗

NASA NTRS1991

Antenna pattern control using impedence surfaces

During this research period, September 16, 1990 to March 15, 1991, a design method for selecting a low-loss impedance material coating for a horn antenna pattern control has been developed. This method and the stepped waveguide technique can be employed to accurately compute the electromagnetic wave phenomenon inside the transition region of the horn antenna, with or without the impedance surfaces, from the feed to the radiating aperture. For moment method solutions of the electric and magnetic current distributions on the radiating aperture and the outer surface of the horn antenna, triangular surface-patch modes are introduced to replace the sinusoidal surface-patch modes as expansion and testing functions to provide a more physical expansion of the current distributions. In the synthesis problem, a numerical optimization process is formulated to minimize the error function between the desired waveguide modes and the modes provided by the horn transition with impedance surfaces. Since the modes generated by the horn transition with impedance surface are computed by analytical techniques, the computational error involved in the synthesis of the antenna pattern is minimum. Therefore, the instability problem can be avoided. A preliminary implementation of the techniques has demonstrated that the developed theory of the horn antenna pattern control using the impedance surfaces is realizable.

Balanis, Constantine A.↗

NASA NTRS1991

IMP: A performance code

IMP (Integrated Mission Program) is a simulation language and code used to model present and future Earth, Moon, or Mars missions. The profile is user controlled through selection from a large menu of events and maneuvers. A Fehlberg 7/13 Runge-Kutta integrator with error and step size control is used to numerically integrate the differential equations of motion (DEQ) of three spacecraft, a main, a target, and an observer. Through selection, the DEQ's include guided thrust, oblate gravity, atmosphere drag, solar pressure, and Moon gravity effects. Guide parameters for thrust events and performance parameters of velocity changes (Delta-V) and propellant usage (maximum of five systems) are developed as needed. Print, plot, summary, and debug files are output.

Dauro, Vincent A., Sr.↗

NASA NTRS1994

Fracture mechanics life analytical methods verification testing

The objective was to evaluate NASCRAC (trademark) version 2.0, a second generation fracture analysis code, for verification and validity. NASCRAC was evaluated using a combination of comparisons to the literature, closed-form solutions, numerical analyses, and tests. Several limitations and minor errors were detected. Additionally, a number of major flaws were discovered. These major flaws were generally due to application of a specific method or theory, not due to programming logic. Results are presented for the following program capabilities: K versus a, J versus a, crack opening area, life calculation due to fatigue crack growth, tolerable crack size, proof test logic, tearing instability, creep crack growth, crack transitioning, crack retardation due to overloads, and elastic-plastic stress redistribution. It is concluded that the code is an acceptable fracture tool for K solutions of simplified geometries, for a limited number of J and crack opening area solutions, and for fatigue crack propagation with the Paris equation and constant amplitude loads when the Paris equation is applicable.

Favenesi, J. A.↗

NASA NTRS1995

Covolume solutions of three dimensional div-curl equations

Delaunay-Voronoi mesh systems provide a generalization of the classical rectangular staggered meshes to unstructured meshes. It is shown how such 'covolume' discretizations may be applied to div-curl systems in three dimensions. Error estimates are proved and confirmed by a numerical illustration.

Nicolaides, R. A.↗

NASA NTRS1994

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗