Errors in hybrid computers
Method is described for reduction of error components in numerical integration, sampling with zero hold order, and execution time delay.
SEARCH · Engineering Papers
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.
Method is described for reduction of error components in numerical integration, sampling with zero hold order, and execution time delay.
Modified program has been developed using improved variation of Encke method which avoids accumulation of round-off errors and avoids numerical ambiguities arising from near-circular orbits of low inclination. Variety of interplanetary trajectory problems can be computed with maximum accuracy and efficiency.
Computer software designed to perform interval arithmetic is described. An interval is defined as the set of all real numbers between two given numbers including or excluding one or both endpoints. Interval arithmetic consists of the various elementary arithmetic operations defined on the set of all intervals, such as interval addition, subtraction, union, etc. One of the main applications of interval arithmetic is in the area of error analysis of computer calculations. For example, it has been used sucessfully to compute bounds on sounding errors in the solution of linear algebraic systems, error bounds in numerical solutions of ordinary differential equations, as well as integral equations and boundary value problems. The described software enables users to implement algorithms of the type described in references efficiently on the IBM 360 system.
A refined elastic analysis of bonded joints which accounts for transverse shear deformation and transverse normal stress was developed to obtain the stresses and displacements in the adherends and in the bond. The displacements were expanded in terms of polynomials in the thicknesswise coordinate; the coefficients of these polynomials were functions of the axial coordinate. The stress distribution was obtained in terms of these coefficients by using strain-displacement and stress-strain relations. The governing differential equations were obtained by integrating the equations of equilibrium, and were solved. The boundary conditions (interface or support) were satisfied to complete the analysis. Single-lap, flush, and double-lap joints were analyzed, along with the effects of adhesive properties, plate thicknesses, material properties, and plate taper on maximum peel and shear stresses in the bond. The results obtained by using the thin-beam analysis available in the literature were compared with the results obtained by using the refined analysis. In general, thin-beam analysis yielded reasonably accurate results, but in certain cases the errors were high. Numerical investigations showed that the maximum peel and shear stresses in the bond can be reduced by (1) using a combination of flexible and stiff bonds, (2) using stiffer lap plates, and (3) tapering the plates.
Interferometric observations taken at 7.8 GHz (gamma approximately = 3.8 cm) with five pairings of antennae of 25 extragalactic radio sources between April, 1972 and May, 1973 are reported. These sources exhibit a broad variety of fine structure from very simple to complex. The total flux and the correlated flux of some of the sources underwent large changes in a few weeks, while the structure and total power of others remained constant during the entire period of observation. Some aspects of the data processing and a discussion of errors are presented. Numerous figures are provided and explained. The individual radio sources are described in detail.
The Heat Capacity Mapping Mission (HCMM), scheduled for launch in 1978, will be three-axis stabilized relative to the earth in a 600-kilometer altitude, polar orbit. The autonomous attitude control system consists of three torquing coils and a momentum wheel driven in response to error signals computed from data received from an infrared horizon sensor and a magnetometer. This paper presents a simple model of the attitude dynamics and derives the equations that determine the stability of the system during both attitude acquisition (acquisition-mode) and mission operations (mission-mode). Modifications to the proposed mission-mode control laws which speed the system's response to transient attitude errors and reduce the steady-state attitude errors are suggested. Numerical simulations are performed to validate the results obtained with the simple model.
Nine simplified analytic models based on average properties accurately predict heat rates for silicon/germanium thermoelectric generators. Solutions from simplified models were compared with those obtained using sophisticated numerical analysis. Maximum errors in calculated heat rate range from about 4 percent to about 0.2 percent. Models also used to calculate power delivered to load and thermodynamic efficiency.
The effects of truncation error on the numerical solution of transonic flows using the full potential equation are studied. The effects of adapting grid point distributions to various solution aspects including shock waves is also discussed. A conclusion is that a rapid change of grid spacing is damaging to the accuracy of the flow solution. Therefore, in a solution adaptive grid application an optimal grid is obtained as a tradeoff between the amount of grid refinement and the rate of grid stretching.
A complete integral formulation is developed for radiative transfer in emitting and absorbing media contained between infinitely long, concentric cylinders with black surfaces. The resulting flux equation is used to state both the gray and nongray radiative equilibrium problems that are solved numerically with controlled error. The results are compared with those from other approaches, both approximate and geometrically limiting, to check and put into perspective the present and prior work. The nongray solution is unique for concentric cylindrical media. The independent parameters of optical thickness and radius ratio are varied to study the emissive power (temperature) distribution, total heat flux, and the required computational time.
The computation of orbits can be done more efficiently by the use of any of several new formulations of the perturbed two body problem which consider the total energy of the orbital system as one of the dependent variables. The total energy is the osculating two body energy plus the potential energy due to perturbing masses. The use of the total energy as the dependent variable instead of the two body energy is a relatively new idea. The advantage of using total energy arises from the fact that the more perturbing potential energy that is accounted for in the total energy variable, the more nearly constant is the total energy. In fact, except for dissipative forces such as drag, the only reason for the total energy not being constant is the rotation or revolution of the perturbing mass. This near constancy of the total energy has the effect of inhibiting error growth during numerical solution. This paper will present the results of an application of total energy formulation to the problem of the precise computation of orbits.
The authors apply a general method of bounding the event error probability of trellis-coded modulation schemes to fading channels and use the effective length and the minimum-squared-product distance to replace the minimum-free-squared-Euclidean distance as code design parameters for Rayleigh and Rician fading channels with a substantial multipath component. They present 8-PSK trellis codes specifically constructed for fading channels that outperform equivalent codes designed for the additive white Gaussian noise channel when v is greater than or equal to 5. For quasiregular trellis codes there exists an efficient algorithm for evaluating event error probability, and numerical results on Pe which demonstrate the importance of the effective length as a code design parameter for fading channels with or without side information have been obtained. This is consistent with the case for binary signaling, where the Hamming distance remains the best code design parameter for fading channels. The authors show that the use of Reed-Solomon block codes with expanded signal sets becomes interesting only for large values of E(s)/N(0), where they begin to outperform trellis codes.
A complete numerical simulation and error analysis was performed for the proposed experiment with the objective of establishing an accurate assessment of the feasibility and the potential accuracy of the measurement of the Lense-Thirring precession. Consideration was given to identifying the error sources which limit the accuracy of the experiment and proposing procedures for eliminating or reducing the effect of these errors. Analytic investigations were conducted to study the effects of major error sources with the objective of providing error bounds on the experiment. The analysis of realistic simulated data is used to demonstrate that satellite laser ranging of two Lageos satellites, orbiting with supplemental inclinations, collected for a period of 3 years or more, can be used to verify the Lense-Thirring precession. A comprehensive covariance analysis for the solution was also developed.
An optimal least squares finite element method is proposed for two dimensional and three dimensional elliptic problems and its advantages are discussed over the mixed Galerkin method and the usual least squares finite element method. In the usual least squares finite element method, the second order equation (-Delta x (Delta u) + u = f) is recast as a first order system (-Delta x p + u = f, Delta u - p = 0). The error analysis and numerical experiment show that, in this usual least squares finite element method, the rate of convergence for flux p is one order lower than optimal. In order to get an optimal least squares method, the irrotationality Delta x p = 0 should be included in the first order system.
Designer and design team productivity improves with skill, experience, and the tools available. The design process involves numerous trials and errors, analyses, refinements, and addition of details. Computerized tools have greatly speeded the analysis, and now new theories and methods, emerging under the label Artificial Intelligence (AI), are being used to automate skill and experience. These tools improve designer productivity by capturing experience, emulating recognized skillful designers, and making the essence of complex programs easier to grasp. This paper outlines the aircraft design process in today's technology and business climate, presenting some of the challenges ahead and some of the promising AI methods for meeting these challenges.
Designer and design team productivity improves with skill, experience, and the tools available. The design process involves numerous trials and errors, analyses, refinements, and addition of details. Computerized tools have greatly speeded the analysis, and now new theories and methods, emerging under the label Artificial Intelligence (AI), are being used to automate skill and experience. These tools improve designer productivity by capturing experience, emulating recognized skillful designers, and making the essence of complex programs easier to grasp. This paper outlines the aircraft design process in today's technology and business climate, presenting some of the challenges ahead and some of the promising AI methods for meeting those challenges.
A Navier-Stokes procedure to calculate the flow about an airfoil at incidence was developed. The parabolized equations are solved in the streamline coordinates generated for an arbitrary airfoil shape using conformal mapping. A modified k-epsilon turbulence model is applied in the entire domain, but the eddy viscosity in the laminar region is suppressed artificially to simulate the region correctly. The procedure was applied to airfoils at various angles of attack, and the results are quite satisfactory for both laminar and turbulent flows. It is shown that the present choice of the coordinate system reduces the error due to numerical diffusion, and that the lift is accurately predicted for a wide range of incidence.
This paper describes the use of a combined model and parameter identification approach, based on modal analysis and artificial intelligence (AI) techniques, for identifying damage or flaws in a rotating truss structure incorporating embedded piezoceramic sensors. This smart structure example is representative of a class of structures commonly found in aerospace systems and next generation space structures. Artificial intelligence techniques of classification, heuristic search, and an object-oriented knowledge base are used in an AI-based model identification approach. A finite model space is classified into a search tree, over which a variant of best-first search is used to identify the model whose stored response most closely matches that of the input. Newly-encountered models can be incorporated into the model space. This adaptativeness demonstrates the potential for learning control. Following this output-error model identification, numerical parameter identification is used to further refine the identified model. Given the rotating truss example in this paper, noisy data corresponding to various damage configurations are input to both this approach and a conventional parameter identification method. The combination of the AI-based model identification with parameter identification is shown to lead to smaller parameter corrections than required by the use of parameter identification alone.
This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.