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.
This article describes the most important features of the elliptic restricted three-body problem. The methods of numerical integration with recurrent power series are developed for both the equations of motion and the variational equations. The conditions for the existence of periodic orbits and families of periodic orbits are also outlined in detail, and finally illustrated with a family of symmetric periodic orbits.
A numerical study of finite-strength, isentropic pressure waves transverse to the axis of a circular cylinder was made for the radial resonant mode. The waves occur in a gas otherwise at rest, filling the cylinder. A method of characteristics was used for the numerical solution. For small but finite amplitudes, calculations indicate the existence of waves of permanent potential form. For larger amplitudes, a shock is indicated to occur. The critical value of the initial amplitude parameter in the power series is found to be 0.06 to 0.08, under various types of initial conditions.
On a plane earth, the ionospheric or tropospheric range error depends only on the total refractivity content or zeroth moment of the refracting layer and the elevation angle. On a spherical earth, however, the dependence is more complex; so for more accurate results it has been necessary to resort to complex ray-tracing calculations. A simple, high-accuracy alternative to the ray-tracing calculation is presented. By appropriate expansion of the angular dependence in the ray-tracing integral in a power series in height, an expression is obtained for the range error in terms of a simple function of elevation angle, E, at the expansion height and of the mth moment of the refractivity, N, distribution about the expansion height. The rapidity of convergence is heavily dependent on the choice of expansion height. For expansion heights in the neighborhood of the centroid of the layer (300-490 km), the expansion to N = 2 (three terms) gives results accurate to about 0.4% at E = 10 deg. As an analytic tool, the expansion affords some insight on the influence of layer shape on range errors in special problems.
The evolution of a disturbance in the presence of a uniform mean velocity gradient is calculated by a power-series solution of the incompressible Navier-Stokes equations. Terms through those in time cubed are retained in the solution. For the initial condition a three-dimensional cosine distribution with two harmonic terms is assumed. The nonlinear interaction of these harmonic terms produces new harmonics which in turn interact. For large velocity gradients the energy of the disturbance can grow with time. The results shed some light on the maintenance or growth of turbulence in a shear flow.
This paper deals with elastic stress intensity factors associated with stress fields near partially debonded interfaces in finite, bi-material plates. Bi-material stress intensity factors K1 and K2 are not associated, respectively, with symmetric and skew-symmetric loading and are functions of an arbitrary length parameter. Any fracture condition involving bi-material stress intensity factors must be independent of the arbitrary length parameter and must be consistent with experiment. We define a new parameter called the interfacial stress intensity factor which satisfies the former condition. A truncated Williams power series stress function for stresses in a cracked bi-material plate, as modified by Sih and Rice, is used with a modified boundary collocation procedure to generate interfacial stress intensity factors for finite bi-material plates.
Standard thermoelectric data have been determined for four thermocouple combinations used below the ice point: commercial types E, K, and T and the special combination of KP vs Au-0.07 at.% Fe. Power series coefficients necessary to generate E = f(T) data for these combinations are given. In addition, extensive tests were carried out to ascertain the variability between materials from different manufacturers. Typical wires processed for low-temperature usage often have slightly thermoelectric properties (up to 1%) from those, nominally identical, wires made for high-temperature applications. Generally, type E is the most satisfactory standardized commercial combination. The Au-0.07 at.% Fe material has been examined in particular detail, because it is not yet an ISA standard material. It has excellent sensitivity in the liquid helium/liquid hydrogen temperature range and a nearly linear thermovoltage at higher temperatures.
Quantitative models of the external magnetospheric field were derived by making least-squares fits to magnetic field measurements from four IMP satellites. The data were fit to a power series expansion in the solar magnetic coordinates and the solar wind-dipole tilt angle, and thus the models contain the effects of seasonal north-south asymmetries. The expansions are divergence-free, but unlike the usual scalar potential expansions, the models contain a nonzero curl representing currents distributed within the magnetosphere. Characteristics of four models are presented, representing different degrees of magnetic disturbance as determined by the range of Kp values. The latitude at the earth separating open polar cap field lines from field lines closing on the dayside is about 5 deg lower than that determined by previous theoretically-derived models. At times of high Kp, additional high latitude field lines are drawn back into the tail.
Description of the theoretical spectroscopic properties of the 2p pi/sub u/ and 3d sigma/sub g/ excited states of the H2/+/ hydrogen molecular ion. Numerical integration of the Schrodinger equation is used to determine vibration-rotation eigenvalues. Dunham power series expansions are used to determine the equilibrium separation, potential coefficients, and spectroscopic constants. The eigenvalues are used to determine delta-G, Bv, Dv, and Hv.
A number of hybrid methods for solving Cauchy problems are described on the basis of an evaluation of advantages of single and multiple-point numerical integration methods. The selection criterion is the principle of minimizing computer time. The methods discussed include the Nordsieck method, the Bulirsch-Stoer extrapolation method, and the method of recursive Taylor-Steffensen power series.
Theoretical vibration frequencies and mode shapes are obtained for skewed plates and helicoidal shells with a cantilever boundary. Using Hamilton's law of varying action, a power series solution is developed to obtain converged numerical results for the five lowest frequencies. Effects of geometrical variables such as aspect ratio, sweep angle and shell radius to thickness ratio are investigated. Accuracy of the solution method is substantiated by comparison with existing skewed plate spherical cap, and conical shell results.
A computer algorithm has been developed to determine the blunt-body flowfields supporting symmetric parabolic and paraboloidal shock waves at infinite free-stream Mach number. Solutions are expressed in an analytic form as high-order power series, in the coordinate normal to the shock, whose coefficients can be determined exactly. Analytic continuation is provided by the use of Pade approximations. Test cases provide solutions of very high accuracy. In the axisymmetric case for gamma equals 715 the solution has been found far downstream, where it agrees with the modified blast-wave results. For plane flow, on the other hand, a limit line appears within the shock layer, a short distance past the sonic line, suggesting the presence of an imbedded shock. Local solutions in the downstream limit are discussed.
The model is derived by making least squares fits to magnetic field measurements from four Imp satellites. It includes four sets of coefficients, representing different degrees of magnetic disturbance as determined by the range of Kp values. The data are fit to a power series expansion in the solar magnetic coordinates and the solar wind-dipole tilt angle, and thus the effects of seasonal north-south asymmetries are contained. The expansion is divergence-free, but unlike the usual scalar potential expansion, the model contains a nonzero curl representing currents distributed within the magnetosphere. The latitude at the earth separating open polar cap field lines from field lines closing on the day side is about 5 deg lower than that determined by previous theoretically derived models. At times of high Kp, additional high-latitude field lines extend back into the tail. Near solstice, the separation latitude can be as low as 75 deg in the winter hemisphere. The average northward component of the external field is much smaller than that predicted by theoretical models; this finding indicates the important effects of distributed currents in the magnetosphere.
Within the context of the restricted problem of three bodies the effects caused by varying the mass ratio of the primaries and the eccentricity of their orbits, upon periodic orbits of the infinitesimal mass which are numerical continuations of circular orbits in the ordinary problem of two bodies are shown. A recursive power series technique is used to numerically integrate the equations of motion as well as the first variational equations in order to generate a two parameter family of perodic orbits and identify the linear stability characteristics. Seven such families are investigated with equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the orbits of the primaries in a range 0.0 to 0.6. Stable orbits are associated with large distances of the infinitesimal mass from the perturbing primary, nearly circular motion of the primaries, and small mass ratios of the primaries. Unstable orbits for the infinitesimal mass are associated with small distances from the perturbing primary, highly elliptic orbits of the primaries and large mass ratios.
Numerical instabilities often arise in the use of high-ordered collocation approximations for numerically solving parabolic partial differential equations. These problems may be reduced by formulations involving evaluation of collocation polynomials rather than combination of the polynomials into a power series. As an illustration, two formulations using shifted Legendre polynomials of order 26 and less are compared.
The general equation for the velocity potential of quasi-one-dimensional acoustic wave motion in a variable area, finite duct with one-dimensional flow is derived by using a perturbation technique. The nonlinear second-order partial differential equation is linearized and then solved, by either a power series expansion method or the Runge-Kutta fourth-order method, for harmonic time dependence. The boundary condition taken at the duct mouth is that of matching the impedance of the duct sound field to that of the radiation field at the duct opening. Three axial Mach number variations along the duct axis are considered and the results obtained are compared with those for the case of constant Mach number, to determine the influence of the axial velocity gradient on sound propagation. The effect of flow on the radiation impedance is also considered.
A third-order solution was developed for the motions of artificial satellites moving in the gravitational field of the earth, whose potential includes the second-, third-, and fourth-order zonal harmonics. Third-order periodic perturbations with fourth-order secular perturbations were derived by the Hori perturbation method. All quantities were expanded into power series of the eccentricity, but the solution was obtained so as to be closed with respect to the inclination. A comparison with the results of numerical integration of the equations of motion indicates that the solution can predict the position of a close-earth, small-eccentricity satellite with an accuracy of better than one cm over a period of one month.
Basic information about the computer code STAGS (Structural Analysis of General Shells) is presented to describe to potential users the scope of the code and the solution procedures that are incorporated. Primarily, STAGS is intended for analysis of shell structures, although it has been extended to more complex shell configurations through the inclusion of springs and beam elements. The formulation is based on a variational approach in combination with local two dimensional power series representations of the displacement components. The computer code includes options for analysis of linear or nonlinear static stress, stability, vibrations, and transient response. Material as well as geometric nonlinearities are included. A few examples of applications of the code are presented for further illustration of its scope.
The steady, three-dimensional, potential flow (neglecting surface tension) of a slender stream of fluid in the presence of gravity is considered and a description is presented of the first part of an asymptotic theory of such a flow. The theory will provide the ingredients for the construction of the solution to many flow problems involving slender streams. The flows and their boundaries are represented as asymptotic power series in the slenderness ratio of the stream. Attention is given to the formulation of the outer expansion problems, first approximations for a jet, a first approximation for pipe flow, and a first approximation for channel flow.