SEARCH · Engineering Papers

Results for “numerical integration”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Discrete distributed strain sensing of intelligent structures

Techniques are developed for the design of discrete highly distributed sensor systems for use in intelligent structures. First the functional requirements for such a system are presented. Discrete spatially averaging strain sensors are then identified as satisfying the functional requirements. A variety of spatial weightings for spatially averaging sensors are examined, and their wave number characteristics are determined. Preferable spatial weightings are identified. Several numerical integration rules used to integrate such sensors in order to determine the global deflection of the structure are discussed. A numerical simulation is conducted using point and rectangular sensors mounted on a cantilevered beam under static loading. Gage factor and sensor position uncertainties are incorporated to assess the absolute error and standard deviation of the error in the estimated tip displacement found by numerically integrating the sensor outputs. An experiment is carried out using a statically loaded cantilevered beam with five point sensors. It is found that in most cases the actual experimental error is within one standard deviation of the absolute error as found in the numerical simulation.

Anderson, Mark S.↗

NASA NTRS1978

The determination of gravity anomalies from geoid heights using the inverse Stokes' formula, Fourier transforms, and least squares collocation

A numerical method for the determination of gravity anomalies from geoid heights is described using the inverse Stokes formula. This discrete form of the inverse Stokes formula applies a numerical integration over the azimuth and an integration over a cubic interpolatory spline function which approximates the step function obtained from the numerical integration. The main disadvantage of the procedure is the lack of a reliable error measure. The method was applied on geoid heights derived from GEOS-3 altimeter measurements in the calibration area of the GEOS-3 satellite.

Rummel, R.↗

NASA NTRS1984

Comparison of free-piston Stirling engine model predictions with RE1000 engine test data

Predictions of a free-piston Stirling engine model are compared with RE1000 engine test data taken at NASA-Lewis Research Center. The model validation and the engine testing are being done under a joint interagency agreement between the Department of Energy's Oak Ridge National Laboratory and NASA-Lewis. A kinematic code developed at Lewis was upgraded to permit simulation of free-piston engine performance; it was further upgraded and modified at Lewis and is currently being validated. The model predicts engine performance by numerical integration of equations for each control volume in the working space. Piston motions are determined by numerical integration of the force balance on each piston or can be specified as Fourier series. In addition, the model Fourier analyzes the various piston forces to permit the construction of phasor force diagrams. The paper compares predicted and experimental values of power and efficiency and shows phasor force diagrams for the RE1000 engine displacer and piston. Further development plans for the model are also discussed.

Tew, R. C., Jr.↗

NASA NTRS1984

Comparison of free-piston Stirling engine model predictions with RE1000 engine test data

Tew, R. C.↗

NASA NTRS1992

On the thermal conductivity due to collisions between relativistic degenerate electrons

A numerical integration is undertaken of the integral expression for the thermal conductivity due to collisions of relativistic degenerate electrons and compared to the limiting-case analytic formula. The integration is designed to encompass all temperature/electron-plasma-frequency temperature ratios between the melting temperature and the Fermi temperature. High accuracy fits are demonstrated by interpolating the values of the integrals of the function and by using an asymptotic function by Urpin and Yakovlev (1980). The numerically integrated expression compares favorably to the limiting-case analytic asymptotic formula by Urpin and Yakovlev, and the results are valid for temperatures above and below the electron-plasma-frequency temperature. The present techniques can be used in stellar opacity calculations and in the study of the propagation of deflagration fronts of compact stellar remnants.

Timmes, F. X.↗

NASA NTRS1992

Comparison Of Two Viscous-Flow Computer Codes

Two viscous-flow computer codes compared by applying them to five test cases of steady-state transonic viscous flows about transonic airfoils. Two codes were: FLOMG, which solves Navier-Stokes differential equations of flow by implementing explicit, Runge-Kutta, finite-volume, multigrid numerical-integration algorithm; and ARC2D, which implements implicit, finite-difference, approximate-factorization, eigenvector-diagonalization numerical-integration algorithm.

Maksymiuk, C. M.↗

NASA NTRS1974

Special perturbations employing osculating reference states

The concept of employing osculating reference position and velocity vectors in the numerical integration of the equations of motion of a satellite is examined. The choice of the reference point is shown to have a significant effect upon numerical efficiency and the class of trajectories described by the differential equations of motion. For example, when the position and velocity vectors on the osculating orbit at a fixed reference time are chosen, a universal formulation is yielded. For elliptical orbits, however, this formulation is unattractive for numerical integration purposes due to Poisson terms (mixed secular) appearing in the equations of motion. Other choices for the reference point eliminate this problem but usually at the expense of universality. A number of these formulations, including a universal one, are considered here. Comparisons of the numerical characteristics of these techniques with those of the Encke method are presented.

Born, G. H.↗

NASA NTRS1973

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 1: Technical discussion

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical solution to the problem on N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The derivation of the second-order solution is summarized by showing the essential steps, some in functional form. The general asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-earth, and interplanetary solutions. The results show that the accuracies of the asymptotic solutions range from an order of magnitude better than conic approximations to that of numerical integration itself. Also, since no iterations are required, the asymptotic boundary value solutions are obtained in a fraction of the time required for comparable numerically integrated solutions. The subject of minimizing the second-order error is discussed, and recommendations made for further work directed toward achieving a uniform accuracy in all applications.

Lancaster, J. E.↗

NASA NTRSDate not supplied

Derivatives Across Adaptive-Step Integration Accurately and Efficiently

The numerical integration of ordinary differential equations (ODEs) is common in trajectory optimization. Previous studies have shown that the use of adaptive-step Runge-Kutta integration schemes can introduce errors in the derivatives due to non-differentiable behavior within the integration scheme itself. The authors previously demonstrated the use of an implicit adaptation to the Runge-Kutta methods to mitigate this issue, but at a significant computational cost. In this work, the authors demonstrate an explicit adaptation to the Runge-Kutta methods to address the issue that does not incur a significant computational cost.

Kaushik Ponnapalli↗

NASA NTRS1986

Analysis and application of frozen orbits for the TOPEX mission

This paper identifies and analyzes frozen orbits for use by the Topography Experiment (TOPEX) mission. Frozen orbits are characterized by almost no long-term change in eccentricity or argument of periapsis. The standard method of frozen orbit prediction is shown to be inadequate for TOPEX inclinations due to the effect of higher degree zonal harmonic gravity terms. A method is described from which long-term motion in mean eccentricity, argument of periapsis, and inclination can be predicted without numerical integration and is used to locate frozen orbits. A zonal gravity field of degree 13 is shown to be necessary and sufficient for TOPEX frozen orbit prediction. Results are verified by numerical integration methods. Frozen orbits are available for most TOPEX orbits under consideration and short-term altitude variations and maximum altitude rates are specified. Alternatives to the frozen orbit which may result in lower altitude variations and rates are briefly examined.

Smith, J. C., Jr.↗

NASA NTRS1990

Trees, bialgebras and intrinsic numerical algorithms

Preliminary work about intrinsic numerical integrators evolving on groups is described. Fix a finite dimensional Lie group G; let g denote its Lie algebra, and let Y(sub 1),...,Y(sub N) denote a basis of g. A class of numerical algorithms is presented that approximate solutions to differential equations evolving on G of the form: dot-x(t) = F(x(t)), x(0) = p is an element of G. The algorithms depend upon constants c(sub i) and c(sub ij), for i = 1,...,k and j is less than i. The algorithms have the property that if the algorithm starts on the group, then it remains on the group. In addition, they also have the property that if G is the abelian group R(N), then the algorithm becomes the classical Runge-Kutta algorithm. The Cayley algebra generated by labeled, ordered trees is used to generate the equations that the coefficients c(sub i) and c(sub ij) must satisfy in order for the algorithm to yield an rth order numerical integrator and to analyze the resulting algorithms.

Crouch, Peter↗

NASA NTRS1972

Numerical analysis of the asymptotic two-point boundary value solution for N-body trajectories.

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical boundary value solution applicable to a broad class of trajectory problems. In addition, the earlier first-order solutions have been extended to second-order to determine if improved accuracy is possible. Comparisons between the asymptotic solution and numerical integration for several lunar and interplanetary trajectories show that the asymptotic solution is generally quite accurate. Also, since no iterations are required, a solution to the boundary value problem is obtained in a fraction of the time required for numerically integrated solutions.

Lancaster, J. E.↗

NASA NTRS2004

Analyzing Dynamics of Cooperating Spacecraft

A software library has been developed to enable high-fidelity computational simulation of the dynamics of multiple spacecraft distributed over a region of outer space and acting with a common purpose. All of the modeling capabilities afforded by this software are available independently in other, separate software systems, but have not previously been brought together in a single system. A user can choose among several dynamical models, many high-fidelity environment models, and several numerical-integration schemes. The user can select whether to use models that assume weak coupling between spacecraft, or strong coupling in the case of feedback control or tethering of spacecraft to each other. For weak coupling, spacecraft orbits are propagated independently, and are synchronized in time by controlling the step size of the integration. For strong coupling, the orbits are integrated simultaneously. Among the integration schemes that the user can choose are Runge-Kutta Verner, Prince-Dormand, Adams-Bashforth-Moulton, and Bulirsh- Stoer. Comparisons of performance are included for both the weak- and strongcoupling dynamical models for all of the numerical integrators.

Hughes, Stephen P.↗

NASA NTRS1985

Optical analysis of parabolic dish concentrators for solar dynamic power systems in space

An optical analysis of a parabolic solar collection system operating in Earth orbit was performed using ray tracing techniques. The analysis included the effects of: (1) solar limb darkening, (2) parametric variation of mirror surface error, (3) parametric variation of mirror rim angle, and (4) parametric variation of alignment and pointing error. This ray tracing technique used numerical integration to combine the effects of rays emanating from different parts of the sun at different intensities with the effects of normally distributed mirror-surface errors to compute the angular intensity distribution of rays leaving the mirror surface. A second numerical integration was then performed over the surface of the parabolic mirror to compute the radial distribution of brightness at the mirror focus. Major results of the analysis included: (1) solar energy can be collected at high temperatures with high efficiency, (2) higher absorber temperatures can be achieved at lower efficiencies, or higher efficiencies can be achieved at lower temperatures, and (3) collection efficiency is near its maximum level across a broad plateau of rim angles from 40 deg to 70 deg.

Jefferies, K. S.↗

NASA NTRS1976

Earth-moon system: Dynamics and parameter estimation; numerical considerations and program documentation

Major activities included coding and verifying equations of motion for the earth-moon system. Some attention was also given to numerical integration methods and parameter estimation methods. Existing analytical theories such as Brown's lunar theory, Eckhardt's theory for lunar rotation, and Newcomb's theory for the rotation of the earth were coded and verified. These theories serve as checks for the numerical integration. Laser ranging data for the period January 1969 - December 1975 was collected and stored on tape. The main goal of this research is the development of software to enable physical parameters of the earth-moon system to be estimated making use of data available from the Lunar Laser Ranging Experiment and the Very Long Base Interferometry experiment of project Apollo. A more specific goal is to develop software for the estimation of certain physical parameters of the moon such as inertia ratios, and the third and fourth harmonic gravity coefficients.

Breedlove, W. J., Jr.↗

NASA NTRS1992

A doubly averaging method for third body perturbations in planet equator coordinates

The first order doubly averaged potential due to third-body gravity is derived in any arbitrary coordinates. The equations of motion are nonsingular at zero eccentricity. The derivation uses a recursive method which allows easy expansion to higher order terms. Instead of using analytical quadrature to obtain the doubly averaged potential, the method presented in this paper simply eliminates the mean anomaly of the perturbed and perturbing bodies by inspection of the recursive formulation. The derivatives of the orbital elements can be numerically integrated rapidly. When a planet equator coordinate system is used, they can be added directly to the derivatives due to gravity harmonics without any coordinate transformation. The method is applied to various high altitude missions. The results are compared with a high precision numerical integration method and are found to provide excellent agreement.

Hwok, Johnny H.↗

NASA NTRS2007

Microphysical Timescales in Clouds and their Application in Cloud-Resolving Modeling

Independent prognostic variables in cloud-resolving modeling are chosen on the basis of the analysis of microphysical timescales in clouds versus a time step for numerical integration. Two of them are the moist entropy and the total mixing ratio of airborne water with no contributions from precipitating particles. As a result, temperature can be diagnosed easily from those prognostic variables, and cloud microphysics be separated (or modularized) from moist thermodynamics. Numerical comparison experiments show that those prognostic variables can work well while a large time step (e.g., 10 s) is used for numerical integration.

Zeng, Xiping↗

NASA NTRS1991

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗