Simplification of second order systems of ordinary differential equations with a turning point.
Second order systems of linear differential equations eliminating all nonessential difficulties associated with initial value problem
SEARCH · Search NASA
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.
Second order systems of linear differential equations eliminating all nonessential difficulties associated with initial value problem
Accelerating convergence of one step methods for numerical solution of ordinary differential equations, discussing various initial value problems
Relativistic solution by Landau procedure of initial value problem for one dimensional electron plasma wave coexisting with charged immobile background
Surface wave terms found in long wavelength limit during solution of initial value problem for semiinfinite hot plasma
Unstable electron plasma initial value problem solved by numerically integrating Vlasov equation in one dimension
Book on linear and ordinary celestial mechanics covering perturbed two body motion, numerical methods, canonical theory and initial value problems
Predictor-corrector schemes are discussed for the numerical integration of initial value problems by computer. Descloux formulas, which are considered modified Nordsieck formulas are studied because of their desirable stability and truncation properties along with the standard Nordsieck formulas. The stability of Nordsieck and modified Nordsieck algorithms are computed, and the results are graphed. Test problems are included as examples for comparing the efficiency of the algorithms.
Near resonance charge exchange between ions of a large ionization potential gas such as helium or neon and vapors of metals such as zinc, cadmium, selenium, or tellurium has produced laser action in the metal ion gas. The possibility of obtaining population inversions in near resonant charge exchange systems (Xe-Ca, Xe-Mg, Xe-Sr, Xe-Ba, Ar-Mg, N-Ca) was investigated. The analysis is an initial value problem that utilizes rate equations for the densities of relevant levels of the laser gas (Ca, Ba, Mg, or Sr) and an electron energy equation. Electron excitation rates are calculated using the Bohr-Thomson approximation for the cross section. Approximations to experimental values of the electron ionization cross section and the ion-atom charge exchange cross section are used. Preliminary results have been obtained for the Ca-Xe system and show that it is possible to obtain gains greater than 10 to the 14th power/m with inversion times up to 8x10 to the minus 7th power second. A possible charge exchange laser system using a MPD arc plasma accelerator is also described.
The distortion of the electron velocity distribution caused by a large amplitude Landau wave is determined analytically for the initial-value problem. The resulting stability of electrostatic perturbations impressed on the evolving plasma is studied. Narrow sidebands of the applied frequency experience consecutive growths of large magnitude during the early stages of the nonlinear wave-particle interaction. The significance of the derived results to both wave propagation experiments and triggered VLF emissions in the magnetosphere is discussed.
Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.
The calculation of the outer inviscid flow about a rectangular wing moving at supersonic speeds is reported. The inviscid equations of motion governing the flow generated by the wing form a set of hyperbolic differential equations. The flow field about the rectangular wing is separated into three regions consisting of the forebody, the afterbody, and the wing wake. Solutions for the forebody are obtained using conical flow techniques while the afterbody and the wing wake regions are treated as initial value problems. The numerical solutions are compared in the two dimensional regions with known exact solutions.
Numerical methods of integration of the equations of motion of a controlled satellite under the influence of gravity-gradient torque are considered. The results of computer experimentation using a number of Runge-Kutta, multi-step, and extrapolation methods for the numerical integration of this differential system are presented, and particularly efficient methods are noted. A large bibliography of numerical methods for initial value problems for ordinary differential equations is presented, and a compilation of Runge-Kutta and multistep formulas is given. Less common numerical integration techniques from the literature are noted for further consideration.
Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.
The protection of composite turbine fan blades against impact forces has prompted the study of dynamic stresses in composites due to transient loads. The mathematical model treats the laminated plate as an equivalent anisotropic material. The use of Mindlin's approximate theory of crystal plates results in five two-dimensional stress waves. Three of the waves are flexural and two involve in-plane extensional strains. The initial value problem due to a transient distributed transverse force on the plate is solved using Laplace and Fourier transforms. A fast computer program for inverting the two-dimensional Fourier transform is used. Stress contours for various stresses and times after application of load are obtained for a graphite fiber-epoxy matrix composite plate. Results indicate that the points of maximum stress travel along the fiber directions.
A Runge-Kutta algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative. The algorithm includes the Fehlberg step control procedure.
For various applications in fluid dynamics, it is assumed that the total temperature is constant. Therefore, the energy equation can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed. It is shown that the system is strictly hyperbolic and well posed for the initial value problems. Boundary conditions are described such that the linearized system is well posed. The Hopscotch method is investigated and numerical results are presented.
The theory of Burdet's focal elements (1969) is outlined. The differential equations are presented, and the initial value problem is described together with the transformation to rectangular coordinates and classical elements. The focal elements are well defined for zero eccentricity and inclination and can be adopted for the computation of elliptic, parabolic, and hyperbolic motion. For the numerical integration of near-geostationary orbits, a comparison of the efficiency is made between focal elements, Kustaanheimo-Stiefel (1971) theory, and rectangular coordinates. For this class of orbits, a higher accuracy has been obtained by integrating elements than integrating rectangular coordinates.
Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.