Search NASA⌕ Search

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.

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

At least 271 records · Page 15

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.↗

Cyclotron resonant scattering in the spectra of gamma-ray bursts

Fits of theoretical spectra from Monte Carlo radiation-transfer calculations to dips at approximately 20 and 40 keV in a spectrum of the gamma-ray burst source GB 880 205 give best-fit values and 68 percent-confidence intervals for the magnetic field of (1.71 + or - 0.07) x 10 to the 12th G, the electron density of (1.2 + or - 0.6) x 10 to the 21st electrons/cm-squared, and the cosine of the viewing angle relative to the field of 0.31 + or - 0.05. The dips observed at approximately 20 keV in the spectra are interpreted as cyclotron resonant scattering, in which electrons undergo radiative 0 to 1 to 0 Landau transitions initiated by photons near the first harmonic. Physical self-consistency fixes the temperature, and the equilibrium temperature equals 5.3 + 0.3 or - 0.2 keV. These results suggest that this gamma-ray burst and many others which exhibit a low-energy dip originate from strongly magnetic neutron stars and are galactic in origin.

Wang, J. C. L.↗

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.↗

One-dimensional cloud fluid model for propagating star formation

The aim of this project was to study the propagation of star formation (SF) with a self-consistent deterministic model for the interstellar gas. The questions of under what conditions does star formation propagate in this model and what are the mechanisms of the propagation are explored. Here, researchers used the deterministic Oort-type cloud fluid model of Scalo and Struck-Marcell (1984, also see the review of Struck-Marcell, Scalo and Appleton 1987). This cloud fluid approach includes simple models for the effects of cloud collisional coalescence or disruption, collisional energy dissipation, and cloud disruption and acceleration as the result of young star winds, HII regions and supernovae. An extensive one-zone parameter study is presented in Struck-Marcell and Scalo (1987). To answer the questions above, researchers carried out one-dimensional calculations for an annulus within a galactic disk, like the so-called solar neighborhood of the galactic chemical evolution. In the calculations the left-hand boundary is set equal to the right hand boundary. The calculation is obviously idealized; however, it is computationally convenient to study the first order effects of propagating star formation. The annulus was treated as if it were at rest, i.e., in the local rotating frame. This assumption may remove some interesting effects of a supersonic gas flow, but was necessary to maintain a numerical stability in the annulus. The results on the one-dimensional propagation of SF in the Oort cloud fluid model follow: (1) SF is propagated by means of hydrodynamic waves, which can be generated by external forces or by the pressure generated by local bursts. SF is not effectively propagated via diffusion or variation in cloud interaction rates without corresponding density and velocity changes. (2) The propagation and long-range effects of SF depend on how close the gas density is to the critical threshold value, i.e., on the susceptibility of the medium.

Titus, Timothy N.↗

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.↗

Localization effects in radiationally disordered high-temperature superconductors: Theoretical interpretation

Theoretical interpretation of recent experiments on radiationally disordered high-temperature superconductors is presented, based on the concepts of mutual interplay of Anderson localization and superconductivity. Microscopic derivation of Ginzburg-Landau coefficients for the quasi-two-dimensional system in the vicinity of localization transition is given in the framework of the self-consistent theory of localization. The 'minimal metallic conductivity' for the quasi-two-dimensional case is enhanced due to a small overlap of electronic states on the nearest neighbor conducting planes. This leads to a stronger influence of localization effects than in ordinary (three-dimensional) superconductors. From this point of view even the initial samples of high-temperature superconductors are already very close to Anderson transition. Anomalies of H(c2) are also analyzed, explaining the upward curvature of H(c2)(T) and apparent independence of dH(c2)/dT (T = T(sub c)) on the degree of disorder as due to localization effects. Researchers discuss the possible reasons of fast T(sub c) degradation due to the enhanced Coulomb effects caused by the disorder induced decrease of localization length. The appearance and growth of localized magnetic moments is also discussed. The disorder dependence of localization length calculated from the experimental data on conductivity correlates reasonably with the theoretical criterion for suppression of superconductivity in the system with localized electronic states.

Goshchitskii, B. N.↗

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.↗

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.↗

Asymptotic Giant Branch stars as a source of short-lived radioactive nuclei in the solar nebula

We carried out a theoretical evaluation of the contribution of Asymptotic Giant Branch (AGB) stars to some short-lived (10(exp 6) less than or equal to Tau-bar less than or equal to 2 x 10(exp 7) yr) isotopes in the Interstellar Medium (ISM) and in the early solar system using stellar model calculations for thermally pulsing evolutionary phases of low-mass stars. The yields of s-process nuclei in the convective He-shell for different neutron exposures tau(sub 0) were obtained, and AGB stars were shown to produce several radioactive nuclei (especially Pd-107, Pb-205, Fe-60, Zr-93, Tc-99, Cs-135, and Hf-182) in diferent amounts. Assuming either contamination of the solar nebula from a single AGB star or models for continuous injection and mixing from many stars into the ISM, we calculate the ratios of radioactive to stable nuclei at the epoch of the Sun's formation. The dilution factor between the AGB ejecta and the early solar system matter is obtained by matching the observed Pd-107/Pd-108 and depends on the value of tau(sub 0). It is found that small masses M(sub He) of He-shell material (10(exp -4)-10(exp -7) solar mass) enriched in s-process nuclei are sufficient to contaminate 1 solar mass of the ISM to produce the Pd-107 found in the early solar system. Predictions are made for all of the other radioactive isotopes. The optimal model to explain several observed radioactive species at different states of the proto-solar nebula involves a single AGB star with a low neutron exposure (tau(sub 0) = 0.03 mbarn(sup -1)) which contaminated the cloud with a dilution factor of M(sub He)/solar mass approximately 1.5 x 10(exp -4). This will also contribute newly synthesized stable s-process nuclei in the amount of approximately 10(exp -4) of their abundances already present in the proto-solar cloud. Variations in the degree of homogenization (approximately 30%) of the injected material may account for some of the small general isotopic anomalies found in meteorites. It is also found that Fe-60 is produced in small but significant quantities that may be sufficient to explain the observations if the time elapsed delta from the contamination of the ISM to the formation of protoplanetary bodies is not higher than delta = 5 x 10(exp 6) yr. If delta is longer, up to 10 x 10(exp 6) yr, this would require the single AGB star to experience enhanced neutron densities (n(sub n) approximately 3 x 10(exp 9)n/cu cm) in the s-processing zone in order to compensate for the branching at Fe-59. The alternative model of long-term continuous ejection of matter from many AGB stars does not appear to match the observations. We also estimate the Al-26 production from the H-shell and find that the Al-26 abundance in the early solar system may be readily explained in a self-consistent manner. Moreover, Al-26 from AGB stars may contribute substantially to the galactic Al-26 gamma-source, while no significant gamma-flux from Co-60 (deriving from Fe-60 decay) is to be expected.

Wasserburg, G. J.↗

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗

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.↗

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.↗

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.↗

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.↗

On the tidal interaction between protoplanets and the primordial solar nebula. II - Self-consistent nonlinear interaction

A method to analyze the full nonlinear response and physical processes associated with the tidal interaction between a binary system and a thin disk in the steady state is presented. Using this approach, density wave propagation, induced by tidal interaction, may be studied for a wide range of sound speeds and viscosities. The effect of self-gravity may also be incorporated. The results of several calculations relevant to the tidal interaction between a protoplanet and the primordial solar nebula are also presented.

Lin, D. N. C.↗

Renormalization group analysis of turbulence

The objective is to understand and extend a recent theory of turbulence based on dynamic renormalization group (RNG) techniques. The application of RNG methods to hydrodynamic turbulence was explored most extensively by Yakhot and Orszag (1986). An eddy viscosity was calculated which was consistent with the Kolmogorov inertial range by systematic elimination of the small scales in the flow. Further, assumed smallness of the nonlinear terms in the redefined equations for the large scales results in predictions for important flow constants such as the Kolmogorov constant. It is emphasized that no adjustable parameters are needed. The parameterization of the small scales in a self-consistent manner has important implications for sub-grid modeling.

Smith, Leslie M.↗

Asymptotic structure and similarity solutions for three-dimensional turbulent boundary layers

The asymptotic structure of the three-dimensional turbulent boundary layer is investigated in the limit of large Reynolds numbers. A self-consistent, but relatively complex, two-layer structure exists and the simplest situation, corresponding to a plane of symmetry, is considered in this paper as a first step. The adjustment of the streamwise velocity to relative rest, through an outer defect layer and then an inner wall layer, is similar to that in two-dimensional flow. The adjustment of the cross-streamwise velocity is more complicated and it is shown that two terms in the expansion are required to obtain useful results, and in particular to obtain the velocity skew angle at the wall near the symmetry plane. The conditions under which self-similarity is achieved near a plane of symmetry are investigated. A set of ordinary differential equations is developed which describe the streamwise and cross-streamwise velocities near a plane of symmetry in a self-similar flow through two orders of magnitude. Calculated numerical solutions of these equations yield trends which are consistent with experimental observations.

Degani, A. T.↗

The three-dimensional turbulent boundary layer near a plane of symmetry

The asymptotic structure of the three-dimensional turbulent boundary layer near a plane of symmetry is considered in the limit of large Reynolds number. A self-consistent two-layer structure is shown to exist wherein the streamwise velocity is brought to rest through an outer defect layer and an inner wall layer in a manner similar to that in two-dimensional boundary layers. The cross-stream velocity distribution is more complex and two terms in the asymptotic expansion are required to yield a complete profile which is shown to exhibit a logarithmic region. The flow in the inner wall layer is demonstrated to be collateral to leading order; pressure-gradient effects are formally of higher order but can cause the velocity profile to skew substantially near the wall at the large but finite Reynolds numbers encountered in practice. The governing set of ordinary differential equations describing a self-similar flow is derived. The calculated numerical solutions of these equations are matched asymptotically to an inner wall-layer solution and the results show trends that are consistent with experimental observations.

Degani, A. T.↗