Design and implementation of models for the double precision trajectory program /DPTRAJ/
Development and numerical integration of differential equations of spacecraft motion used in Double Precision Trajectory Program
Spier, G. W.↗
SEARCH · Engineering Papers
Results for “numerical differentiation”
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.
Collisionless sheath near an electrode in the presence of a magnetic field
Magnetic field effects on collisionless plasma sheath near planar electrode, solving electric and magnetic potentials nonlinear differential equations by numerical scheme
Hassan, H. A.↗
Numerical stability in problems of linear algebra.
Mathematical problems are introduced as mappings from the space of input data to that of the desired output information. Then a numerical process is defined as a prescribed recurrence of elementary operations creating the mapping of the underlying mathematical problem. The ratio of the error committed by executing the operations of the numerical process (the roundoff errors) to the error introduced by perturbations of the input data (initial error) gives rise to the concept of lambda-stability. As examples, several processes are analyzed from this point of view, including, especially, old and new processes for solving systems of linear algebraic equations with tridiagonal matrices. In particular, it is shown how such a priori information can be utilized as, for instance, a knowledge of the row sums of the matrix. Information of this type is frequently available where the system arises in connection with the numerical solution of differential equations.
Babuska, I.↗
A generalized theory for the design of contraction cones and other low speed ducts
A generalization of the Tsien method of contraction cone design is described. The design velocity distribution is expressed in such a form that the required high order derivatives can be obtained by recursion rather than by numerical or analytic differentiation. The method is applicable to the design of diffusers and converging-diverging ducts as well as contraction cones. The computer program is described and a FORTRAN listing of the program is provided.
Barger, R. L.↗
Computer program grade for design and analysis of graded-porosity heat-pipe wicks
A computer program for numerical solution of differential equations that describe heat pipes with graded-porosity fibrous wicks is discussed. A mathematical problem is provided with a summary of the input and output steps used to solve it. The program is also applied to the analysis of a typical heat pipe.
Eninger, J. E.↗
Arbitrarily Curved and Twisted Space Beams
A derivation of the equations which govern the deformation of an arbitrarily curved and twisted space beam is presented. These equations differ from those of the classical theory in that (1) extensional effects are included; (2) the strain-displacement relations are derived; and (3) the expressions for the stress resultants are developed from the strain displacement relations. It is shown that the torsional stress resultant obtained by the classical approach is basically incorrect except when the cross-section is circular. The governing equations are given in the form of first-order differential equations. A numerical algorithm is given for obtaining the natural vibration characteristics and example problems are presented.
Hunter, W. F.↗
Special Perturbations Using Back-Correction Methods of Numerical Integration
A new class of linear multistep methods for numerical integration of differential equations is reported that permits satellite computation solutions to be corrected at certain points in the past as the integration advances in time. Algorithms have been developed for the solution of both first- and second-order differential equations. The back correction method appears to be more efficient than classical methods when dominant and perturbing forces can be separated.
Feagin, T.↗
A note on the generation of phase plane plots on a digital computer
A technique is presented for generating phase plane plots on a digital computer which circumvents the difficulties associated with more traditional methods of numerical solving nonlinear differential equations. In particular, the nonlinear differential equation of operation is formulated.
Simon, M. K.↗
Dynamic crack propagation and arrest in orthotropic DCB fiber composite specimens
An orthotropic double cantilever beam (DCB) model is used to study dynamic crack propagation and arrest in 90 deg unidirectional Hercules AS/3501-6 graphite fiber epoxy composites. The dynamic fracture toughness of the composite is determined from tests performed on the long-strip specimen and DCB crack arrest experiments are conducted. By using the dynamic fracture toughness in a finite-difference solution of the DCB governing partial differential equations, a numerical solution of the crack propagation and arrest events is computed. Excellent agreement between the experimental and numerical crack arrest results are obtained.
Williams, J. H., Jr.↗
Air-to-air combat analysis - Review of differential-gaming approaches
The problem of evaluating the combat performance of fighter/attack aircraft is discussed, and the mathematical nature of the problem is examined. The following approaches to air combat analysis are reviewed: (1) differential-turning differential game and (2) coplanar differential game. Selected numerical examples of these approaches are presented. The relative advantages and disadvantages of each are analyzed, and it is concluded that air combat analysis is an extremely difficult mathematical problem and that no one method of approach is best for all purposes. The paper concludes with a discussion of how the two approaches might be used in a complementary manner.
Ardema, M. D.↗
Conference on Complex Turbulent Flows: Comparison of Computation and Experiment, Stanford University, Stanford, CA, September 14-18, 1981, Proceedings. Volume 2 - Taxonomies, reporters' summaries, evaluation, and conclusions
Computational techniques for simulating turbulent flows were explored, together with the results of experimental investigations. Particular attention was devoted to the possibility of defining a universal closure model, applicable for all turbulence situations; however, conclusions were drawn that zonal models, describing localized structures, were the most promising techniques to date. The taxonomy of turbulent flows was summarized, as were algebraic, differential, integral, and partial differential methods for numerical depiction of turbulent flows. Numerous comparisons of theoretically predicted and experimentally obtained data for wall pressure distributions, velocity profiles, turbulent kinetic energy profiles, Reynolds shear stress profiles, and flows around transonic airfoils were presented. Simplifying techniques for reducing the necessary computational time for modeling complex flowfields were surveyed, together with the industrial requirements and applications of computational fluid dynamics techniques.
Kline, S. J.↗
The Control of Float Zone Interfaces by the Use of Selected Boundary Conditions
Since the float zone process involves two solid-melt interfaces, possible gas interfaces, heat and mass transfers, various driving forces and complex heating sources, an analysis of the entire process would be very complex. For an initial investigation, a more feasible approach is to examine each component of the process separately. The three principal components are: (1) the shapes of the melt and solid-melt interfaces, (2) the heat and mass transfers, and (3) the heating and cooling sources. This study combined facets of all three components. The purpose of this effort was to study and compute the surface boundary conditions required to give flat (float zone) FZ solid-melt interfaces. This study was undertaken in two phases. The first phase was to investigate the solid zones surface boundary conditions required for flat solid-melt interfaces when given the melt zone surface boundary conditions. The second phase complemented the first and was to investigate the melt zone surface boundary conditions required for flat solid-melt interfaces if given the solid zones surface boundary conditions. Dual integral transform methods were used in both phases; in addition, the use of various numerical methods for differential equations and linear systems of equations were required.
Foster, L. M.↗
Positioning Rotors In Turbine Flowmeters
Lengths of wakes roughly proportional to thickness of vanes. Mathematical model simplifies analysis of effects of flow-straightening vanes in turbine flowmeter. Yields numerical solution of differential equations of flow for quick examination of efforts of thicknesses of vanes and rate of flow on extent of wake behind vanes. From examination, minimum distance at which flowmeter rotor placed behind vanes determined.
Lynch, Edward D.↗
Goertler vortices in supersonic boundary layers
The problem of Goertler vortices in compressible boundary layers over concave walls is studied by solving parabolic partial differential equations. The numerical scheme used is based upon a fourth-order accurate compact difference scheme. The results show that both the velocity and temperature fluctuations may lead to a Goertler vortex. The vortex growth rates determined from the present method are found to differ somewhat from those given by normal mode solution. The effects of cooling, heating and pressure gradient are also studied. At a Mach number of 3.5, cooling has a small destabilizing effect. The effect of adverse pressure gradient is also found to be destabilizing on a Mach 8 boundary layer.
Spall, R. E.↗
GPS-aided gravimetry at 30 km altitude from a balloon-borne platform
A balloon-borne experiment, flown at 30 km altitude over New Mexico, was used to test dynamic differential Global Positioning System (GPS) tracking in support of gravimetry at high-altitudes. The experiment package contained a gravimeter (Vibrating String Accelerometer), a full complement of inertial instruments, a TI-4100 GPS receiver and a radar transponder. The flight was supported by two GPS receivers on the ground near the flight path. From the 8 hour flight, about a forty minute period was selected for analysis. Differential GPS phase measurements were used to estimate changes in position over the sample time interval, or average velocity. In addition to average velocity, differential positions and numerical averages of acceleration were obtained in three components. Gravitational acceleration was estimated by correcting for accelerations due to translational motion, ignoring all rotational effects.
Lazarewicz, Andrew R.↗
A new second-order integration algorithm for simulating mechanical dynamic systems
A new integration algorithm which has the simplicity of Euler integration but exhibits second-order accuracy is described. In fixed-step numerical integration of differential equations for mechanical dynamic systems the method represents displacement and acceleration variables at integer step times and velocity variables at half-integer step times. Asymptotic accuracy of the algorithm is twice that of trapezoidal integration and ten times that of second-order Adams-Bashforth integration. The algorithm is also compatible with real-time inputs when used for a real-time simulation. It can be used to produce simulation outputs at double the integration frame rate, i.e., at both half-integer and integer frame times, even though it requires only one evaluation of state-variable derivatives per integration step. The new algorithm is shown to be especially effective in the simulation of lightly-damped structural modes. Both time-domain and frequency-domain accuracy comparisons with traditional integration methods are presented. Stability of the new algorithm is also examined.
Howe, R. M.↗
Radiation transport around axisymmetric blunt body vehicles using a modified differential approximation
A moment method for computing 3D radiative transport in axisymmetric thermochemical nonequilibrium flows is developed. The method uses the P-1 approximation to reduce the governing system of integro-differential equations to a coupled set of partial differential equations. The numerical solution of these equations for realistic variations of the radiation properties is discussed. Representative results from the method are shown and compared to tangent slab calculations. The agreement between the transport methods is found to be about 10 percent in the stagnation region, with the difference increasing along the flank of the vehicle.
Hartung, Lin C.↗
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.↗