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.
Results for trajectory analysis, fuel analysis, and analysis of aerobraking maneuvers are presented. Calculations are included for thickness and mass of the propellant tanks; determination of weight (avionics module); numerical integration program; hydrostatic pressure due to acceleration; and propellant tank pressurization.
The problem of electromagnetic wave generation by an electrodynamic tethered satellite system is important both for the ordinary operation of such systems and for their possible application as orbiting transmitters. The tether's ionospheric circuit closure problem is closely linked with the propagation of charge-carrying electromagnetic wave packets away from the tethered system. Work is reported which represents a step towards a solution to the problem that takes into account the effects of boundaries and of vertical variations in plasma density, collision frequencies, and ion species. The theory of Alfen wave packet generation by an electrodynamic tethered system in an infinite plasma medium is reviewed, and brief summary of previous work on the problem is given. The consequences of the presence of the boundaries and the vertical nonuniformity are then examined. One of the most significant new features to emerge when ion-neutral collisions are taken into account is the coupling of the Alfven waves to the fast magnetosonic wave. This latter wave is important, as it may be confined by vertical variations in the Alfven speed to a sort of leaky ionospheric wave guide, the resonances of which could be of great importance to the signal received on the Earth's surface. The infinite medium solution for this case where the (uniform) geomagnetic field makes an arbitrary angle with the vertical is taken as the incident wave-packet. Even without a full solution, a number of conclusions can be drawn, the most important of which may be that the electromagnetic field associated with the operation of a steady-current tethered system will probably be too weak to detect on the Earth's surface, even for large tethered currents. This is due to the total reflection of the incident wave at the atmospheric boundary and the inability of a steady-current tethered system to excite the ionospheric wave-guide. An outline of the approach to the numerical problem is given. The use of numerical integrations and boundary conditions consistent with a conducting Earth is proposed to obtain the solution for the horizontal electromagnetic field components at the boundary of the ionosphere with the atmospheric cavity.
The progress of two efforts in coding solutions of Navier-Stokes equations is summarized. The first effort concerns a 3-D space marching parabolized Navier-Stokes (PNS) code being modified to compute the supersonic mixing flow through an internal/external expansion nozzle with multicomponent gases. The 3-D PNS equations, coupled with a set of species continuity equations, are solved using an implicit finite difference scheme. The completed work is summarized and includes code modifications for four chemical species, computing the flow upstream of the upper cowl for a theoretical air mixture, developing an initial plane solution for the inner nozzle region, and computing the flow inside the nozzle for both a N2/O2 mixture and a Freon-12/Ar mixture, and plotting density-pressure contours for the inner nozzle region. The second effort concerns a full Navier-Stokes code. The species continuity equations account for the diffusion of multiple gases. This 3-D explicit afterbody code has the ability to use high order numerical integration schemes such as the 4th order MacCormack, and the Gottlieb-MacCormack schemes. Changes to the work are listed and include, but are not limited to: (1) internal/external flow capability; (2) new treatments of the cowl wall boundary conditions and relaxed computations around the cowl region and cowl tip; (3) the entering of the thermodynamic and transport properties of Freon-12, Ar, O, and N; (4) modification to the Baldwin-Lomax turbulence model to account for turbulent eddies generated by cowl walls inside and external to the nozzle; and (5) adopting a relaxation formula to account for the turbulence in the mixing shear layer.
The various types of local-global finite-element problems point out the need to develop a new generation of software. First, this new software needs to have a complete analysis capability, encompassing linear and nonlinear analysis of 1-, 2-, and 3-dimensional finite-element models, as well as mixed dimensional models. The software must be capable of treating static and dynamic (vibration and transient response) problems, including the stability effects of initial stress, and the software should be able to treat both elastic and elasto-plastic materials. The software should carry a set of optional diagnostics to assist the program user during model generation in order to help avoid obvious structural modeling errors. In addition, the program software should be well documented so the user has a complete technical reference for each type of element contained in the program library, including information on such topics as the type of numerical integration, use of underintegration, and inclusion of incompatible modes, etc. Some packaged information should also be available to assist the user in building mixed-dimensional models. An important advancement in finite-element software should be in the development of program modularity, so that the user can select from a menu various basic operations in matrix structural analysis.
The thermal creation, fusion, evaporation, and destruction of non-topological solitons (NTS) after a phase transition in the early universe is considered. By defining and following NTS statistical equilibrium and departures from it, and depending on particle physics parameters, one of three possible scenarios occurs. If reaction rates are high enough, a period of equilibrium occurs and relic abundances are determined by the freeze-out temperature. Equilibrium first drives most NTS's into their constituents (free phi particles) and then causes rapid fusion into large NTS's. If freeze-out occurs during the first phase, the NTS's are almost entirely destroyed, while if it occurs during the second phase, solitosynthesis occurs and NTS's may be cosmically relevant. For slow reaction rates the NTS's are born frozen out and have the abundance determined by the phase transition. Analytic approximations for determining the abundances are developed, and tested by numerically integrating a reaction network in an expanding universe. Unfortunately, for most of the parameter space considered, solito-destruction/evaporation occurs.
Numerical integrations of the structure function of tropospheric inhomogeneities have been performed to assess the impact of water vapor fluctuations on connected element interferometry (CEI). The expectation value of the RMS troposphere error for a differential spacecraft-quasar observation is derived by integrating the spatial refractivity structure function along raypaths to both the spacecraft and quasar from two spatially separated sites. Correlations between the tropospheric conditions at the two sites, which can become significant for short baseline observations, are fully accounted for in this calculation. Temporal effects are treated by assuming a frozen-flow model in which a fixed spatial distribution blows over both sites. Two nominal observation scenarios are considered, along with variations to study the dependence of the resultant differential troposphere errors on baseline length, observation time, source separation angle, and elevation. Consecutive differential observations are found to be almost completely uncorrelated, implying that averaging many repeated differential observations can quickly reduce the troposphere error.
This paper presents the results of the latest fits of both analytical theory and numerically integrated Neptunian satellite orbits to Earth-based astrometric observations. Ephemerides based on the integrated orbits will be used by the Voyager project for pre-encounter planning and analysis until late 1988 when the final pre-encounter ephemerides will be produced. As a by-product of the orbit fits, new estimates of the Neptune mass, the second zonal harmonic of Neptune, and the pole orientation of Neptune are obtained. The theory and integrated orbits are compared with each other and with orbits obtained by previous investigators.
The dynamical behavior of N coorbital satellites moving with the same average mean motion around a primary has been studied both analytically and by numerical integrations for N from 2 to 9 satellites. Simplified dynamical equations are used to determine the different stationary configurations and their local stability against infinitesimal perturbations. The ring of equally spaced identical satellites is found to be locally unstable for N of six or less, while for N from two to eight there exists another stable compact stationary configuration with separations of 60 deg or less between adjacent satellites. For N of seven or more the equally spaced configuration becomes locally stable, and for N of nine or more it is the only stationary configuration. The motion becomes chaotic for large-amplitude perturbations. The chaotic motion fills a restricted region of phase space whose outer boundary is determined by the maximum velocity curve.
Cross-sections of (O-18)2 at 79 K have been obtained from photoabsorption measurements at various pressures throughout the wavelength region 177.8-197.8 nm with a 6.65 m photoelectric scanning spectrometer equipped with a 2400 lines/mm grating and having an instrumental width (FWHM) of 0.0013 nm. The measured absorption cross-sections of the Schumann-Runge bands (14,0) through (2,0) are, with the exception of the (12,0) band, independent of the instrumental width. The measured cross-sections are presented graphically here and are available at wavenumber intervals of about 0.1/cm as numerical compilations stored on magnetic tape. Band oscillator strengths of those bands have been determined by direct numerical integration of the measured absolute cross-sections and are in excellent agreement with these theoretically calculated values.
Motivated by considerations of the solar toroidal magnetic field, the behavior of a layer of uniform magnetic field embedded in a convectively stable atmosphere is studied. Since the field can support extra mass, such a configuration is top-heavy and thus instabilities of the Rayleigh-Taylor type can occur. For both static and rotating basic states, the evolution of the interchange modes (no bending of the field lines) is followed by integrating numerically the nonlinear compressible MHD equations. The initial Rayleigh-Taylor instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The subsequent motions are determined primarily by the strong interactions between vortices which are responsible for the rapid disruption of the magnetic layer.
The present evaluation of the basic techniques and illustrative results of charged particle-modeling numerical codes suitable for particle acceleration at oblique, fast-mode collisionless shocks emphasizes the treatment of ions as test particles, calculating particle dynamics through numerical integration along exact phase-space orbits. Attention is given to the acceleration of particles at planar, infinitessimally thin shocks, as well as to plasma simulations in which low-energy ions are injected and accelerated at quasi-perpendicular shocks with internal structure.
The time-dependent Navier-Stokes equations, including second-order turbulence model, are numerically integrated by using four-stage Runge-Kutta scheme to predict the steady-state supersonic flow structures in ramjet dump combustors. The formulation is derived for reacting flows with finite-rate chemistry. In the present study, it is firstly attempted to assess the accuracy of existing high-order turbulence model in supersonic flows. The comparison shows reasonable agreement between calculated and measured data in terms of velocity distributions. It is indicated that a modified constant C-mu for calculating turbulent eddy viscosity is needed in the supersonic flow regime and the adaptive meshing is preferred to capture the recirculation zone. In the reacting flow calculation, the results from a test case of hydrogen and air combustion at premixed conditon show that the rearward facing step is able to increase flow residence time and stabilize the flame in supersonic flows.
Numerical integrations are used here to show that small dust grains can be temporarily captured into exterior orbit-orbit resonances with the earth, lasting from less than 10,000 years to more than 100,000 years. Grains with radii of 30-100 microns, orbiting in planes less than 10 deg from the plane of the solar system and with orbital eccentricities of less than 0.3, are captured most easily. It is argued that there should be an approximately toroidal cloud of particles, derived mostly from the asteroid belt, trapped into a variety of these exterior resonances. The cloud is mostly beyond the earth's orbit, but includes it.
A linear flux approach is developed for a finite element thermal-structural analysis of steady-state thermal and structural problems. The element fluxes are assumed to vary linearly in the same form as the element unknown variables, and the finite element matrices are evaluated in closed form. Since numerical integration is avoided, significant computational time saving is achieved. Solution accuracy and computational speed improvements are demonstrated by solving several two-and three-dimensional thermal-structural examples.
Vast ethereal rings composed of micron-sized particles have been detected about the three giant planets so far visited by Voyager. An orbital evolution mechanism for gravitationally dominated circumplanetary dust is identified here that may temporarily overwhelm other processes. Evolution rates that agree with numerical integrations are derived, and it is suggested how this process could affect the location of the faint rings evident today.
The nonlinear spatial evolution, from the source to the atmosphere, of the auroral electron beam and the beam-generated electrostatic whistler noise was studied, calculating changes in beam parameters from equations for the conservation of total particle and wave energy and momentum flux density. Wave power fluxes were calculated by numerically integrating the wave kinetic equations, and the levels of beam-generated noise were determined by using thermal levels of Cerenkov radiation as a source. It was found that beam parameters evolve on ionospheric scale lengths, and their positive slope feature in velocity space is maintained over altitudes measured in thousands of kilometers of altitude, even though they can generate wave energy density fluxes sufficient to modify the ionospheric density profile.
A method is presented for calculating the electromagnetic wave field on the earth's surface associated with the operation of an electrodynamic tethered satellite system of constant or slowly varying current in the upper ionosphere. The wave field at the ionospheric boundary and on the earth's surface is obtained by numerical integration. The results suggest that the ionospheric waves do not propagate into the atmosphere and that the image of the Alfven wings from a steady-current tether should be greatly broadened on the earth's surface.
The dynamical properties and spectral characteristics of the ionospheric interchange instability in the collisional and inertial regimes are investigated numerically and analytically, using a two-dimensional fluid model. This model, which describes the collisional and inertial dynamics of the F region ionosphere, is similar to the models of Ott (1978) and Kintner and Seyler (1985). By numerically integrating the model equations, the nonlinear evolution of three interchange unstable configurations in the collisional and inertial regimes are followed, and their dynamical and spectral properties are studied. The solutions indicate that the dynamical properties of a collision-dominated regime are distinctly different from those of an inertia-dominated regime. The qualitative characteristics of a plasma bubble in both regime-types are described.