SEARCH · Search NASA

Results for “Self-consistent calculations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Computer models of the spacecraft wake

Until recently, computations of space plasma flow over a spacecraft have been unstable for ratios of spacecraft dimension to Debye length typical of the low Earth orbit environment. Calculations are presented of the spacecraft/environment interaction based on two computer codes, MACH and POLAR. MACH, an inside-out particle tracking code, was developed for the purpose of validating the physics of POLAR in regimes where these are no comprehensive theoretical or experimental results. While the spacecraft which can be treated by MACH are restricted to simple geometries, the methodology is more fundamental than POLAR. MACH generates self-consistent solutions within the context of quasisteady Vlasov plasma flow and achieves Debye ratios previously unobtainable. POLAR uses a three-dimensional finite-element representation of the vehicle in a staggered mesh. The plasma sheath is modeled by outside-in particle tracking. Solutions for the plasma flow, wake and vehicle charging are obtained by Vlasov-Poisson iteration; charge stabilization techniques make the results virtually insensitive to the Debye ratio. POLAR reproduces the Laframboise static plasma solutions for sperical probes and fits the Makita-Kuriki probe data for spheres in a flowing plasma in regions where comparisons are valid. POLAR and MACH solutions for the particle and electrostatic potential structure of the wake of a charged disk in a low-altitude flow are shown for Mach numbers 4, 5, and 8. New features of the solutions include ion focussing in the wake and a definitive determination of the sheath edge in the wake which shows that the sheath is not an equipotential.

Rubin, A. G.↗

NASA NTRS1981

Hot relativistic winds and the Crab Nebula

Efforts to formulate a self-consistent model of pulsar magnetospheres which links the particle source near the pulsar to the outflowing relativistic wind and couples the wind to the surrounding nebula are reviewed. The use of a relativistic MHD wind is recommended to account for global photon emission and the invisibility of the method of plasma transport. Consideration of a magnetic monopole relativistic wind due to an axially symmetric aligned rotator is combined with calculations of the initial velocity of the wind to show that the flow velocity in such a model will never exceed Mach 1. Extending the solution to the case of a hot relativistic wind at supersonic speeds is noted to yield results consistent with observations of the Crab Nebula

Fujimura, F. S.↗

NASA NTRS1979

Energy deposition in the earth's atmosphere due to impact of solar activity-generated disturbances

Energy deposition in and dynamic responses of the terrestrial atmosphere to solar flare-generated shocks and other physical processes - such as particle precipitation and local heating - are investigated self-consistently in the context of hydrodynamics, the problem being treated as an initial boundary-value problem. It is extremely difficult to construct a general model for the line solar activity-magnetosphere-atmosphere; however, a limited model for this link is possible. The paper describes such a model, and presents some results on energy deposition into the earth's atmosphere due to solar activity-generated disturbances. Results from the present calculations are presented and discussed.

Wu, S. T.↗

NASA NTRS1976

Accretion onto magnetized neutron stars - Structure and interchange instability of a model magnetosphere

A self-consistent model is analyzed for the spherical infall of weakly magnetized plasma into the magnetosphere of a slowly rotating, strongly magnetized neutron star. It is shown that spherical infall is probably a good approximation for X-ray sources which accrete from a stellar wind. The location of the standoff shock which halts the hypersonic infall is estimated along with the emission from the shocked layer. The location of the equilibrium magnetopause and the structure of the magnetic field within it are calculated; it is found that the magnetic poles are true cusps and that the entry of gas due to equilibrium flow across a cusp is almost certainly dominated by the interchange instability near the magnetic equator. The energy principle is applied to derive necessary conditions for the occurrence of this instability. The results indicate that the strong magnetic-pressure gradient stabilizes the gas unless moderately strong radiative cooling takes place and that the cooled plasma enters the magnetosphere as long filaments capable of moving between field lines. The rate at which the equilibrium magnetopause can 'absorb' mass and momentum is derived, the validity of the approximations employed is discussed, and the likely evolution of the sinking filaments is outlined to show that the spatial distribution of the plasma is determined mainly by the dynamics and thermodynamics of the filaments rather than the magnetic-field structure.

Arons, J.↗

NASA NTRS1983

Crater ejecta scaling laws - Fundamental forms based on dimensional analysis

Self-consistent scaling laws are developed for meteoroid impact crater ejecta. Attention is given to the ejection velocity of material as a function of the impact point, the volume of ejecta with a threshold velocity, and the thickness of ejecta deposit in terms of the distance from the impact. Use is made of recently developed equations for energy and momentum coupling in cratering events. Consideration is given to scaling of laboratory trials up to real-world events and formulations are developed for calculating the ejection velocities and ejecta blanket profiles in the gravity and strength regimes of crater formation. It is concluded that, in the gravity regime, the thickness of an ejecta blanket is the same in all directions if the thickness and range are expressed in terms of the crater radius. In the strength regime, however, the ejecta velocities are independent of crater size, thereby allowing for asymmetric ejecta blankets. Controlled experiments are recommended for the gravity/strength transition.

Housen, K. R.↗

NASA NTRS1977

The continuous UV flux of Alpha Lyrae - Non-LTE results

Non-LTE calculations for the ultraviolet C I and Si I continuous opacity show that LTE results overestimate the importance of these sources of opacity and underestimate the emergent flux in Alpha Lyr. The largest errors occur between 1100 and 1160 A, where the predicted flux in non-LTE is as much as 50 times larger than in LTE, in reasonable accord with Copernicus observations. The discrepancy between LTE models and observations has been interpreted to result from the existence of a chromosphere. Until a self-consistent non-LTE model atmosphere becomes available, such an interpretation is premature.

Snijders, M. A. J.↗

NASA NTRS1977

Boundary-value problem of configurations with compressible free vortex flow

A self-consistent version of the compressible boundary-value problem for configurations with leading-edge vortex separation is formulated, based on the assumption that the compressible flow field is controlled by the linearized potential equation. The stream surface boundary condition and the zero pressure jump condition of the compressible free vortex flows are analyzed; application of the Goethert rule permits the compressible nonlinear boundary-value problem for the subsonic flow domain to be transformed into an equivalent nonlinear incompressible problem. The compressibility corrections developed are used in numerical calculations of subsonic leading-edge vortex flows about planar wing geometries. The sample calculations, employing an inviscid flow model in which the wing and vortex sheets are represented by piecewise continuous quadratic doublet sheet distributions, are applicable to high subsonic Mach numbers.

Brune, G. W.↗

NASA NTRS1981

Theory of dense hydrogen

Equations of state for molecular and metallic hydrogen are calculated in an accurate and comparable manner by using the density-functional method applied to static lattices. By isolating the electronic energies in this way, attention can be drawn to the relative importance of the protonic degrees of freedom. From the equation-of-state results it can be concluded that a remnant molecular pairing is preferred in a band-overlap metallic state, complete dissociation occurring only at very high densities. An accurate determination of the pressure needed to achieve complete dissociation is shown to require a self-consistent treatment of electron and proton degrees of freedom.

Chakravarty, S.↗

NASA NTRS1978

Effect of temperature coupling on ozone depletion prediction

The effects of chlorine perturbations on both the temperature and the ozone distribution in the stratosphere have been studied using a simplified radiative-photochemical model. The model solves the hydrostatic equation for total density in a self-consistent manner as the temperature is changed. Radiative coupling is found to have a significant effect on both the thermal structure and the ozone distribution, particularly in the 35-50-km region. By increasing the ClX mixing ratio by 5.0 ppbv, the temperature in this region is decreased by 5 to 10 K with a slight increase below 30 km. The local ozone depletion around 40 km due to added ClX is smaller compared with the estimate made by keeping the temperature fixed to the ambient condition. However, the integrated effect of radiative coupling is to increase the calculated column ozone depletion by 15% to 25% in this model.

Chandra, S.↗

NASA NTRS1974

Temperature and density structure of thermal proton flows

Thermal proton flows along magnetic field lines are an important feature of magnetosphere-ionosphere coupling. In this paper we report the results of a theoretical study of the thermal structure of such flows. The adopted steady state model is based upon O+, H+, and electrons with self-consistent solutions for the separate O+, H+, and electron temperatures, the O+ and H+ densities and the H+ drift velocity. Through investigation of a number of parameters affecting the model, it is shown that Joule heating arising from the flow of H+ through O+ preferentially heats H+, so that the H+ temperature is substantially greater than the O+ temperature. Low O+ densities characteristic of the trough region appear to give high H+/O+ temperature ratios. Typical O+ densities characteristic of polar wind flow regions give moderate H+/O3 temperature ratios. The Mach number of H+ outflow is substantially reduced in the present models in comparison with the older fixed temperature calculations.

Banks, P. M.↗

NASA NTRS1978

On the stability of self-consistent large amplitude waves in a cold plasma. I - Transverse circularly polarized waves in the absence of a large scale magnetic field

It is demonstrated that a self-consistent circularly polarized wave in an otherwise field-free homogeneous cold plasma is unstable to small amplitude perturbations. For either an electron-positron plasma or an electron-proton plasma the instability rate is at least about the order of the effective plasma frequency when the bulk flow speed is zero. For finite bulk flow speeds of the plasma, it is shown that the electron-positron plasma is unstable, again with a growth rate of the order of the effective plasma frequency; it is also shown that the electron-proton plasma is unstable (at least at small wave numbers, k) with a growth rate proportional to k. The calculated instability rates are conservative, for other modes not investigated here may be more unstable. The results of these calculations bear directly on the understanding of plasma systems thought to be driven by large amplitude waves.

Lee, M. A.↗

NASA NTRS1981

Simulation of the current-driven electrostatic ion cyclotron instability

Results are presented of a self-consistent particle simulation of the growth and saturation of electrostatic ion cyclotron (EIC) waves along auroral field lines. The driving mechanism for the waves is an initial drifting Maxwellian distribution for the electrons. The magnetic field is taken to lie in the simulation plane, so that modes such as the ion acoustic instability and the EIC instability can be excited, thus making it possible to study the behavior of EIC waves in the presence of ion acoustic turbulence. A dc electric field, which tends to keep the current constant along the field lines, is included. It is shown that the EIC instability grows to a moderate level before saturating due to plateau formation. The perpendicular ion temperature increases by typically 40%, with an increase arising from bulk heating only. The anomalous resistivity, associated with the EIC instability is determined and compared with previous calculations and with that due to the ion acoustic instability. Implications of the results for observation of ion cyclotron waves on auroral field lines are discussed.

Pritchett, P. L.↗

NASA NTRS1977

Departures of the electron energy distribution from a Maxwellian in hydrogen. I - Formulation and solution of the electron kinetic equation. II - Consequences

The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.

Shoub, E. C.↗

NASA NTRS1986

Nonlinear analysis of a relativistic beam-plasma cyclotron instability

A self-consistent set of nonlinear and relativistic wave-particle equations are derived for a magnetized beam-plasma system interacting with electromagnetic cyclotron waves. In particular, the high-frequency cyclotron mode interacting with a streaming and gyrating electron beam within a background plasma is considered in some detail. This interaction mode may possibly find application as a high-power source of coherent short-wavelength radiation for laboratory devices. The background plasma, although passive, plays a central role in this mechanism by modifying the dielectric properties in which the magnetized electron beam propagates. For a particular choice of the transverse beam velocity (i.e., the speed of light divided by the relativistic mass factor), the interaction frequency equals the nonrelativistic electron cyclotron frequency times the relativistic mass factor. For this choice of transverse beam velocity the detrimental effects of a longitudinal beam velocity spread is virtually removed. Power conversion efficiencies in excess of 18 percent are both analytically calculated and obtained through numerical simulations of the wave-particle equations. The quality of the electron beam, degree of energy and pitch angle spread, and its effect on the beam-plasma cyclotron instability is studied.

Sprangle, P.↗

NASA NTRS1978

Structure of current sheets in magnetic holes at 1 AU

Current density profiles in several types of interplanetary magnetic holes were calculated assuming that the currents flow in planar sheets and that the magnetic field varies only in the direction normal to the sheet. The planarity was verified in four holes which were observed by two suitably spaced spacecraft. The structure of the current sheets ranges from very simple in some holes to very complex in others. The observed structures are found to be qualitatively consistent with models based on self-consistent solutions of Vlasov's equation and Maxwell's equations. Examples of complex, irregular magnetic holes are also presented, and they are shown to contain multiple, current sheets in which currents flow parallel to one another.

Fitzenreiter, R. J.↗

NASA NTRS1979

Negative streamer development in FEP teflon

A computational model is developed which describes the evolution and propagation of an ionizing front (negative streamer) in solid materials. The ionization front consists of drifting avalanching electrons moving self-consistently under the influence of their own space-charge field together with an applied external field. The required input information for the model consists of the functional dependence of the macroscopic transport coefficients on the local electric field, the initial conditions for beginning the calculation, and the strength of the applied field. A computational approach for specifying the transport coefficients and initional conditions is also described. The approach has been implemented by constructing three computer codes which sequentially interface, beginning with single electron scattering, and ending with streamer development. Computational results are presented for model calculations in Teflon. The overall model is perceived to provide a picture of the initiation phase of a propagating discharge in electron-irradiated dielectrics.

Beers, B. L.↗

NASA NTRS1982

Quantum-chemical modeling of smectite clays

A self-consistent charge extended Hueckel program is used in modeling isomorphic substitution of Al(3+) by Na(+), K(+), Mg(2+), Fe(2+), and Fe(3+) in the octahedral layer of a dioctahedral smectite clay, such as montmorillonite. Upon comparison of the energies involved in the isomorphic substitution, it is found that the order for successful substitution is as follows: Al(3+), Fe(3+), Mg(2+), Fe(2+), Na(+), which is equivalent to Ca(2+), and then K(+). This ordering is found to be consistent with experimental observation. The calculations also make it possible to determine the possible penetration of metal ions into the clay's 2:1 crystalline layer. For the cases studied, this type of penetration can occur at elevated temperatures into regions where isomorphic substitution has occurred with metal ions that bear a formal charge of less than 3+. The computed behavior of the electronic structure in the presence of isomorphic substitution is found to be similar to behavior associated with semiconductors.

Aronowitz, S.↗

NASA NTRS1986

Production of flickering aurora and field-aligned electron flux by electromagnetic ion cyclotron waves

Recent observations have suggested that flickering aurora is produced by a modulation of the field-aligned component of the electron flux within an auroral arc. It is proposed that a portion of the field-aligned electrons are of ionospheric origin and that these electrons are accelerated and their flux modulated by electromagnetic ion cyclotron waves that occur below the main acceleration region on auroral arc field lines. A model of the electromagnetic ion cyclotron wave shows that the parallel phase velocity of the wave increase as the wave propagates toward the ionosphere. A test particle calculation shows that ionospheric electrons trapped or reflected by the wave are accelerated to energies of several keV and that their flux is modulated at the wave frequency. The relative amplitudes of the model wave electric fields are consistent with the observations of small-scale low-frequency ionospheric and magnetospheric electric fields near auroral arcs of approximately 10 mV/m and 100 mV/m, respectively. The large-amplitude ion cyclotron waves also produce a ponderomotive force and a self-consistent ambipolar electric field. Energy considerations show that the downward energy flux in the electromagnetic ion cyclotron wave can be several percent of the total downward auroral electron energy flux.

Temerin, M.↗