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.
Program uses Navier-Stokes and continuity equations for incompressible, viscous fluid as the basic equations governing reorientation flow dynamics. Program can simulate curved as well as straight-walled boundaries; has ability to calculate both free-surface and confined flows; and can be used in either cylindrical or plane geometry.
Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.
Data obtained with the quadrupole mass spectrometer aboard the OGO 6 satellite show the effects of significant surface interaction processes, including nearly complete recombination of incoming atomic oxygen on the walls of the instrument antechamber plus absorption and desorption of oxygen and carbon monoxide. The observed data are fit by solving the time-dependent continuity equations accounting for production and loss of atomic oxygen, molecular oxygen, and (in the case of mass 28) carbon monoxide. The surface parameters that best fit the data are selected and applied to the determination of ambient densities and their estimated errors.
Using temperature fields derived from the Mariner 9 infrared spectroscopy experiment, the Martian atmospheric tidal pressure and wind fields are calculated. Temperature as a function of local time, latitude, and atmospheric pressure level is obtained by secular and longitudinal averaging of the data. The resulting temperature field is approximated by a spherical harmonic expansion, retaining one symmetric and one asymmetric term for wavenumber zero and wavenumber one. Vertical averaging of the linearized momentum and continuity equations results in an inhomogeneous tidal equation for surface pressure fluctuations with the driving function related to the temperature field through the geopotential function and the hydrostatic equation. Solutions of the tidal equation show a diurnal fractional pressure amplitude approximately equal to one half of the vertically averaged diurnal fractional temperature amplitude.
Examination of the perturbed continuity equation when the perturbations are the result of an internal atmospheric gravity wave in the E region. The transient response of the ionization is interpreted as the gradient instability and the values of the vertical and horizontal wave numbers that will induce it are plotted for various heights. Only in the presence of westward directed electric fields, which are believed to occur only at night, will the gravity waves induce the gradient instability. Approximate analytic expressions are obtained for the permitted wave numbers as well as for the instability growth times. In the course of this analysis it is shown that in the D region all irregularities, even those that are field-aligned, will tend to move with the ion velocity.
The simultaneous time-dependent continuity equations for O(+), H(+) and He(+) in the low latitude F-region are solved. Account is taken of E x B drift, a meridional neutral wind, and ion-ion and ion-neutral drag. The calculated profiles of O(+) and H(+) concentrations at 1630 LT are in fair agreement with the observations of Hanson et al. The He(+) field-aligned velocity is almost matched to the O(+) field-aligned velocity and, above the chemical equilibrium region and around the He(+) peak, the He(+) concentration is determined largely by production and transport. There is disagreement between the theoretical vertical He(+) profile and the profile observed by Hanson et al. Satisfactory agreement is obtained with Taylor's satellite results at fixed height for O(+) and H(+). It is found that the He(+) concentration is greater in the winter hemisphere than in the summer hemisphere, even if the neutral helium distribution is symmetrical about the Equator. The He(+) results are consistent with Taylor's results.
Data obtained with the quadrupole mass spectrometer aboard the Ogo 6 satellite show the effects of significant surface interaction processes, including nearly complete recombination of incoming atomic oxygen on the walls of the instrument antechamber plus adsorption and desorption of oxygen and carbon monoxide. The observed data are fit by solving the time-dependent continuity equations accounting for production and loss of atomic oxygen, molecular oxygen, and (in the case of mass 28) carbon monoxide. The surface parameters that best fit the data are selected and applied to the determination of ambient densities and their estimated errors.
An explanation was developed for the formation, near midnight at midlatitudes, of a broad electron density layer extending approximately from 120 to 180 km and usually referred to as the intermediate E layer. The responsible mechanism is believed to be the converging vertical ion drifts resulting from winds of the solar semidiurnal tide. Numerical solutions of the continuity equation appropriate to the intermediate layer is described for particular models of ion drift, diffusion coefficents, and ionization production. Analysis of rocket observations of the layer show that the ionization rate is highly correlated with the planetary geomagnetic index, K sub p. Particle flux measurements support the idea that energetic electrons are the principal source of this ionization. A semiconductor spectrometer experiment for investigation of the particle flux, spectrum, and angular properties was designed and successfully flown on a Nike Apache rocket. A detailed description of the theory, design, and calibration of the experiment and some preliminary results presented.
We have solved the coupled momentum and continuity equations for NO(+), O2(+), and O(+) ions in the E- and F-regions of the ionosphere. This theoretical model has enabled us to examine the relative importance of various processes that affect molecular ion densities. We find that transport processes are not important during the day; the molecular ions are in chemical equilibrium at all altitudes. At night, however, both diffusion and vertical drifts induced by winds or electric fields are important in determining molecular ion densities below about 200 km. Nitric oxide plays an important role in determining the NO(+) to O2(+) ratio in the E-region, particularly at night. Nocturnal sources of ionization are required to maintain the E-region through the night. Vertical velocities induced by expansion and contraction of the neutral atmosphere are too small to affect ion densities at any altitude.
Construction of a model neutral and ion exosphere for a planet weakly interacting with the solar wind. The model is constructed in general terms and is then specialized to possible neutral and ion exospheres for the planet Mercury. The neutral exosphere model allows for density and temperature variations and for rotation at the exobase. The ion exosphere is produced by ionization of the neutral exosphere in the solar wind, and its density distribution is obtained by solving the continuity equation in the drift approximation. Applying to Mercury a surface temperature distribution inferred from infrared data and a vanishing bound neutral flux at the base, He and He(+) density distributions are found. When the He atmosphere of Mercury is due entirely to the surface bombardment by solar wind He(2+), the resulting He(+) density is found to vary from 0.15 to 0.001 per cu cm over the range from 1.5 to 5 planetocentric radii on the dayside. These densities are found to be detectable by typical solar-wind plasma instruments.
The determination of the relaxation time of a temperature perturbation in the stratosphere must take into account the effects of the absorption of solar energy by ozone, while the ozone density itself is dependent on temperature. A photochemical model, consisting of continuity equations for each of 29 constituents with reaction rates and adjustment times, is used to obtain the vertical distribution of ozone. The altitude range from 35 to 60 km in non-Arctic regions is shown to be in approximate joint radiative-photochemical equilibrium. The temperature and ozone distributions are thus well buffered. Modest changes in ozone and temperature were calculated on the basis of this model from large changes in cooling rates and reaction rates, and the results are shown to be more in line with actual observations. The importance of estimating the mixing ratios for NO and H2O is emphasized.
Using temperature fields derived from the Mariner 9 infrared spectroscopy experiment, the Martian atmospheric tidal pressure and wind fields are calculated. Temperature as a function of local time, latitude, and atmospheric pressure level is obtained by secular and longitudinal averaging of the data. The resulting temperature field is approximated by a spherical harmonic expansion, retaining one symmetric and one asymmetric term each for wavenumber zero and wavenumber one. Vertical averaging of the linearized momentum and continuity equations results in an inhomogeneous tidal equation for surface pressure fluctuations with the driving function related to the temperature field through the geopotential function and the hydrostatic equation. Solutions of the tidal equation show a diurnal fractional pressure amplitude approximately equal to one-half the vertically averaged diurnal fractional temperature amplitude.
The time-dependent continuity equations, including diffusion, were solved for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to stable auroral red (SAR) arcs. The results show that molecular nitrogen is excited vibrationally to the degree that the rate constant for the ionospheric loss process, O(+) + N2 yields NO(+) + N, is increased by as much as a factor of 7.6 at F2 region altitudes. It was found that deviations from the energetically equivalent Boltzmann distribution were large, causing the rate constant to be as much as 1.6 times the rate constant calculated for the Boltzmann distribution. These results indicate that SAR arc intensities as small as 58 R can produce noticeable increases in the ionosphere ion-atom interchange reaction rate and hence in the rate of loss of ionospheric electrons. It is suggested that the observed decrease of electron density in the F2 region in SAR arcs can probably be explained by enhanced reaction rates for ion-atom interchange between O(+) and N2 caused by vibrational excitation of molecular nitrogen by electron impact.
An analysis of the oscillatory fluid flow in the vicinity of a circular orifice with a steady grazing flow is presented. The study is similar to that of Hersh and Rogers but with the addition of the grazing flow. Starting from the momentum and continuity equations, a considerably simplified system of partial differential equations is developed with the assumption that the flow can be described by an oscillatory motion superimposed upon the known steady flow. The equations are seen to be linear in the region where the grazing flow effects are dominant, and a solution and the resulting orifice impedance are presented for this region. The nonlinearity appears to be unimportant for the usual conditions found in aircraft noise suppressors. Some preliminary conclusions of the study are that orifice resistance is directly proportional to grazing flow velocity (known previously from experimental data) and that the orifice inductive (mass reactance) end correction is not a function of grazing flow. This latter conclusion is contrary to the widely held notion that grazing flow removes the effect of the orifice inductive end correction. This conclusion also implies that the experimentally observed total inductance reduction with grazing flow might be in the flow within the orifice rather than in the end correction.
A general device-analysis program has been utilized to study the efficiency of silicon solar cells. The analysis is applied to specific geometries of both n(+)-p and n(+)-p-p(+) solar cells, and involves a numerical solution of the basic transport and continuity equations. This approach allows solutions that are free of typical limiting assumptions involved in solving the device equations as well as solutions relating to lifetime, mobility variations, and diffused-region profiles. The analysis includes available empirical information on diffusion length, mobility, and lifetime as a function of doping, as well as a Gaussian profile for the diffused region. Results are presented which illustrate the limitations of efficiency as a function of doping. It is found that the maximum efficiencies for both types of cell converge at lower resistivities to around 16% with air-mass-zero radiation and a single-layer absorbing-SiO antireflecting film. It is also found that the minority-carrier lifetime, both in the n(+) surface and p-type bulk regions, presents serious limitations to conversion efficiency, particularly in the low-resistivity cells.
The continuity equation for cosmic ray propagation is used to derive a set of linear equations interrelating the fluxes of multiply charged nuclei as observed at any particular part of the galaxy. The derivation leads to model independent definitions for cosmic ray storage time, mean density of target nuclei and effective mass traversed. The set of equations form a common framework for comparisons of theories and observations. As an illustration, it is shown that there exists a large class of propagation models which give the same result as the exponential path length model. The formalism is shown to accommodate dynamic as well as equilibrium models of production and propagation.
A new method to study dynamic behavior of the ionosphere-protonosphere coupling within the plasmasphere is developed and used to calculate ion distributions above 500 km using observed electron densities at 500 km. The method is based on the relation between the total magnetic flux tube content of H(+) above some reference height (e.g., 3000 km) and the H(+) flux at that height, which is uniquely determined by the coupled momentum and continuity equations for the quasi-steady state. The O(+) profile is perturbed from diffusive equilibrium only by ion drag with H(+). The time dependency is taken into account by applying self-consistent boundary conditions which are obtained from the solution of the time-dependent equation for the tube content of H(+). Calculations are carried out for two models. The first model deals with the ion dynamics in a tube of force which rotates with the earth. In the second model, effects of the cross-L plasma drifts are considered. Both models produce diurnal variations of H(+) and O(+) densities near 1000 km similar to observational results.
New flux-corrected transport algorithms are described for solving generalized continuity equations. These techniques were developed by requiring that the finite difference formulas used ensure positivity for an initially positive convected quantity. Thus FCT is particularly valuable for fluid-like problems with strong gradients or shocks. Repeated application of the same subroutine to mass, momentum, and energy conservation equations gives a simple solution of the coupled time-dependent equations of ideal compressible fluid dynamics without introducing an artificial viscosity. FCT algorithms span Eulerian, sliding-rezone, and Lagrangian finite difference grids in several coordinate systems. The latest FCT techniques are fully vectorized for parallel/pipeline processing.