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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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 NTRS2019

Predicting the X-Ray Spectra of Stellar-mass Black Holes from Simulations

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a harm3d GRMHD simulation are post-processed by an improved version of the Monte Carlo radiation transport code pandurata (in the corona) and the Feautrier solver PTRANSX (in the disk), with xstar subroutines. The codes are run in a sequential, iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer, including features, such as the Fe Kα emission line and corresponding K-edge absorption trough, due to disk reprocessing of coronal power. For the example cases we consider—a non-rotating 10 M(sub ⊙) black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington—we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1–3 keV), steep power laws extending to high energy (Γ = 2.7–4.5), and prominent, asymmetric Fe Kα emission lines with equivalent widths in the range 40–400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry—no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., is needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

Kinch, Brooks E.↗
NASA NTRS2000

On the Origin and Evolution of Stellar Chromospheres, Coronae and Winds

This grant was awarded by NASA to The University of Alabama in Huntsville (UAH) to construct state-of-the-art, theoretical, two-component, chromospheric models for single stars of different spectral types and different evolutionary status. In our proposal, we suggested to use these models to predict the level of the "basal flux", the observed range of variation of chromospheric activity for a given spectral type, and the decrease of this activity with stellar age. In addition, for red giants and supergiants, we also proposed to construct self-consistent, purely theoretical wind models, and used these models to investigate the origin of "dividing lines" in the H-R diagram. In the following, we describe our completed work. We have accomplished the first main goal of our proposal by constructing first purely theoretical, time-dependent and two-component models of stellar chromospheres.1 The models require specifying only three basic stellar parameters, namely, the effective temperature, gravity and rotation rate, and they take into account non-magnetic and magnetic regions in stellar chromospheres. The non-magnetic regions are heated by acoustic waves generated by the turbulent convection in the stellar subphotospheric layers. The magnetic regions are identified with magnetic flux tubes uniformly distributed over the entire stellar surface and they are heated by longitudinal tube waves generated by turbulent motions in the subphotospheric and photospheric layers. The coverage of stellar surface by magnetic regions (the so-called filling factor) is estimated for a given rotation rate from an observational relationship. The constructed models are time-dependent and are based on the energy balance between the amount of mechanical energy supplied by waves and radiative losses in strong Ca II and Mg II emission lines. To calculate the amount of wave energy in the non-magnetic regions, we have used the Lighthill-Stein theory for sound generation.

Musielak, Z. E.↗
NASA NTRS2014

The Hot Horizontal-Branch Stars in NGC288 - Effects of Diffusion and Stratification on Their Atmospheric Parameters*

Context. NGC288 is a globular cluster with a well developed blue horizontal branch covering the so-called u-jump which indicates the onset of diffusion. It is therefore well suited to study the effects of diffusion in blue horizontal branch (HB) stars. Aims. We compare observed abundances to predictions from stellar evolution models calculated with diffusion and from stratified atmospheric models. We verify the effect of using stratified model spectra to derive atmospheric parameters. In addition we investigate the nature of the overluminous blue HB stars around the u-jump. Methods. We define a new photometric index sz from uvby measurements that is gravity sensitive between 8 000K and 12 000 K. Using medium-resolution spectra and Stroemgren photometry we determine atmospheric parameters (Teff, logg) and abundances for the blue HB stars. We use both homogeneous and stratified model spectra for our spectroscopic analyses. Results. The atmospheric parameters and masses of the hot HB stars in NGC288 show a behaviour seen also in other clusters for temperatures between 9 000K and 14 000 K. Outside this temperature range, however, they follow rather the results found for such stars in (omega)Cen. The abundances derived from our observations are for most elements (except He and P) within the abundance range expected from evolutionary models that include the effects of atomic diffusion and assume a surface mixed mass of 10(exp −7) M. The abundances predicted by stratified model atmospheres are generally significantly more extreme than observed, except for Mg. The use of stratified model spectra to determine effective temperatures, surface gravities and masses moves the hotter stars to a closer agreement with canonical evolutionary predictions. Conclusions. Our results show definite promise towards solving the long-standing issue of surface gravity and mass discrepancies for hot HB stars, but there is still much work needed to arrive at a self-consistent solution.

horizontal branch (HB)↗
NASA NTRS1994

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

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

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

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

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

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

The structure of a three-dimensional turbulent boundary layer

The three-dimensional turbulent boundary layer is shown to have a self-consistent two-layer asymptotic structure in the limit of large Reynolds number. In a streamline coordinate system, the streamwise velocity distribution is similar to that in two-dimensional flows, having a defect-function form in the outer layer which is adjusted to zero at the wall through an inner wall layer. An asymptotic expansion accurate to two orders is required for the cross-stream velocity which is shown to exhibit a logarithmic form in the overlap region. The inner wall-layer flow is collateral to leading order but the influence of the pressure gradient, at large but finite Reynolds numbers, is not negligible and can cause substantial skewing of the velocity profile near the wall. Conditions under which the boundary layer achieves self-similarity and the governing set of ordinary differential equations for the outer layer are derived. The calculated solution of these equations is matched asymptotically to an inner wall-layer solution and the composite profiles so formed describe the flow throughout the entire boundary layer. The effects of Reynolds number and cross-stream pressure gradient on the crossstream velocity profile are discussed and it is shown that the location of the maximum cross-stream velocity is within the overlap region.

Degani, A. T.↗
NASA NTRS2000

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗
NASA NTRS2006

Ion Thruster Discharge Performance Per Magnetic Field Topography

DC-ION is a detailed computational model for predicting the plasma characteristics of rain-cusp ion thrusters. The advanced magnetic field meshing algorithm used by DC-ION allows precise treatment of the secondary electron flow. This capability allows self-consistent estimates of plasma potential that improves the overall consistency of the results of the discharge model described in Reference [refJPC05mod1]. Plasma potential estimates allow the model to predict the onset of plasma instabilities, and important shortcoming of the previous model for optimizing the design of discharge chambers. A magnetic field mesh simplifies the plasma flow calculations, for both the ions and the secondary electrons, and significantly reduces numerical diffusion that can occur with meshes not aligned with the magnetic field. Comparing the results of this model to experimental data shows that the behavior of the primary electrons, and the precise manner of their confinement, dictates the fundamental efficiency of ring-cusp. This correlation is evident in simulations of the conventionally sized NSTAR thruster (30 cm diameter) and the miniature MiXI thruster (3 cm diameter).

plasma processes↗
NASA NTRS2011

Design and Impacts of Land-Biogenic-Atmosphere Coupling in the NASA-Unified WRF (NU-WRF) Modeling System

Land-Atmosphere coupling is typically designed and implemented independently for physical (e.g. water and energy) and chemical (e.g. biogenic emissions and surface depositions)-based models and applications. Differences in scale, data requirements, and physics thus limit the ability of Earth System models to be fully coupled in a consistent manner. In order for the physical-chemical-biological coupling to be complete, treatment of the land in terms of surface classification, condition, fluxes, and emissions must be considered simultaneously and coherently across all components. In this study, we investigate a coupling strategy for the NASA-Unified Weather Research and Forecasting (NU-WRF) model that incorporates the traditionally disparate fluxes of water and energy through NASA's LIS (Land Information System) and biogenic emissions through BEIS (Biogenic Emissions Inventory System) and MEGAN (Model of Emissions of Gases and Aerosols from Nature) into the atmosphere. In doing so, inconsistencies across model inputs and parameter data are resolved such that the emissions from a particular plant species are consistent with the heat and moisture fluxes calculated for that land cover type. In turn, the response of the atmospheric turbulence and mixing in the planetary boundary layer (PBL) acts on the identical surface type, fluxes, and emissions for each. In addition, the coupling of dust emission within the NU-WRF system is performed in order to ensure consistency and to maximize the benefit of high-resolution land representation in LIS. The impacts of those self-consistent components on' the simulation of atmospheric aerosols are then evaluated through the WRF-Chem-GOCART (Goddard Chemistry Aerosol Radiation and Transport) model. Overall, this ambitious project highlights the current difficulties and future potential of fully coupled. components. in Earth System models, and underscores the importance of the iLEAPS community in supporting improved knowledge of processes and innovative approaches for models and observations.

Tan, Qian↗