Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

The rate parameters for coupled vibration-dissociation in a generalized SSH approximation

A theoretical study of vibrational excitations and dissociations of nitrogen undergoing a nonequilibrium relaxation process upon heating and cooling is reported. The rate coefficients for collisional induced vibrational transitions and transitions from a bound vibrational state into a dissociative state have been calculated using an extension of the theory originally proposed by Schwarz (SSH) et al. (1952). High-lying vibrational states and dissociative states were explicitly included but rotational energy transfer was neglected. The transition probabilities calculated from the SSH theory were fed into the master equation, which was integrated numerically to determine the population distribution of the vibrational states as well as bulk thermodynamic properties. The results show that: (1) the transition rates have a minimum near the middle of the bound vibrational levels, causing a bottleneck in the vibrational relaxation and dissociation rates; (2) high vibrational states are always in equilibrium with the dissociative state; (3) for the heating case, only the low vibrational states relax according to the Landau-Teller theory; (4) for the cooling case, vibrational relaxation cannot be described by a rate equation; (5) Park's (1985, 1988) two-temperature model is approximately valid; and (6) the average vibrational energy removed in dissociation is about 30 percent of the dissociation energy.

Sharma, Surendra P.↗
NASA NTRS1988

Developments in the simulation of a geopotential research mission

An essential element of any satellite system that will be used to recover information about high degree and order terms in the geopotential model of the earth is one or more low altitude (about 160 km) satellites equipped with a drag compensation mechanism. To study the effects of various error sources and to test the new theoretical and numerical techniques that will be developed for such a mission, two simulated scenarios have been used with a geopotential model complete to degree and order 360 and Encke's method to numerically integrate the equations of motion. The two scenarios include the low-low dual satellite system with integrated, one-way Doppler measurements and the single satellite system with gradiometer measurements. The simulations include a reference orbit which is assumed to be available from conventional tracking systems.

Schutz, B. E.↗
NASA NTRS1988

The nonlinear development of Goertler vortices in growing boundary layers

The development of Goertler vortices in boundary layers over curved walls in the nonlinear regime is investigated. The growth of the boundary layer makes a parallel flow analysis impossible except in the high wavenumber regime so in general the instability equations must be integrated numerically. Here the spanwise dependence of the basic flow is described using Fourier series expansion while the normal and streamwise variations are taken into account using finite differences. The calculations suggest that a given disturbance imposed at some position along the wall will eventually reach a local equilibrium state essentially independent of the initial conditions. In fact, the equilibrium state reached is qualitatively similar to the large amplitude high wave-number solution described asymptotically by Hall (1982). In general, it is found that the nonlinear interactions are dominated by a mean field type of interaction between the mean flow and the fundamental. Thus, even though higher harmonics of the fundamental are necessarily generated, most of the disturbance energy is confined to the mean flow correction and the fundamental. A major result of the calculations is finding that the downstream velocity field develops a strongly inflectional character as the flow moves downstream. The latter result suggests that the major effect of Goertler vortices on boundary layers of practical importance might be to make them highly receptive to rapidly growing Rayleigh modes of instability.

Hall, Philip↗
NASA NTRS1989

Comprehensive Silicon-Solar-Cell Program

Comprehensive Silicon Solar Cell Computer Modeling (SICELL) program simulates silicon solar cell. Predicts device parameters as efficiency, voltage-vs.-current characteristic fill factor, and temperature coefficients of parameters. Technique used similar to numerical-integration methods, but commonality described by use of mesh-point field. Validation studies show accuracies of simulations range from 0.08 percent to 3.6 percent for 27 experimental data points over temperature range of 300K to 421K. Results obtained by use of 10 mesh points in n- and p-type regions and for two iterations. SICELL runs interactively on VAX computer under VMS and written in VAX/VMS FORTRAN 77.

Lamorte, Michael F.↗
NASA NTRS1989

Theoretical investigation of EM wave generation and radiation in the ULF, ELF, and VLF bands by the electrodynamic orbiting tether

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.

Estes, Robert D.↗
NASA NTRS1989

Navier-Stokes analysis of cold scramjet-afterbody flows

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.

Baysal, Oktay↗
NASA NTRS1989

A review of some problems in global-local stress analysis

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.

Nelson, Richard B.↗
NASA NTRS1989

Solitosynthesis: Cosmological evolution of non-topological solitons

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.

Griest, Kim↗
NASA NTRS1989

The effect of spatial and temporal wet-troposphere fluctuations on connected element interferometry

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.

Edwards, C. D.↗
NASA NTRS1988

Satellite ephemerides for the Voyager Neptune encounter

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.

Jacobson, Robert A.↗
NASA NTRS1988

The dynamics of coorbital satellite systems

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.

Salo, H.↗
NASA NTRS1988

High resolution absorption cross-sections and band oscillator strengths of the Schumann-Runge absorption bands of isotopic oxygen, (O-18)2, at 79 K

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.

Yoshino, K.↗
NASA NTRS1988

The nonlinear breakup of a magnetic layer - Instability to interchange modes

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.

Cattaneo, F.↗
NASA NTRS1988

Computer modeling of test particle acceleration at oblique shocks

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.

Decker, Robert B.↗
NASA NTRS1989

Numerical investigation of chemically reacting flows in ramjet dump combustors

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.

Hsieh, Kwang-Chung↗
NASA NTRS1989

A solar system dust ring with earth as its shepherd

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.

Jackson, A. A.↗
NASA NTRS1989

Thermal-structural finite element analysis using linear flux formulation

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.

Pandey, Ajay K.↗