Search NASA⌕ Search

SEARCH · Search NASA

Results for “Self-consistent methods”

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 73 records · Page 4

Modeling of infrared flux spectra from disk-shaped interstellar dust clouds

Although there is growing observational and theoretical evidence for many disk-shaped objects of astrophysical interest, spherical geometry is assumed in most radiative transfer models. Recently we generalized the quasi-diffusion method developed by Leung (1975, 1976) for spherical geometry to solve the problem of scattering, absorption, and reemission by dust grains in a medium of disk geometry. The method is applicable to a variety of astronomical sources whose dynamics are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, bipolar-flow molecular clouds, accretion disks and disk galaxies. Using this technique and realistic grain opacities we construct theoretical models to determine self-consistently the dust temperature distribution and infrared emission from disk-shaped, quiescent dark globules heated externally by the intersteller radiation field. The effects of the following parameters on the temperature structure and the emergent spectrum are studied: grain type (graphite and silicate), optical depth, density inhomogeneity, and degree of disk flattening. The disk models are characterized by an aspect ratio and an optical depth at 0.55 microns. For inhomogeneous models, a gaussian density distribution is assumed such that the ratio of central to surface density is 100. To study the effects of source geometry, we also compare results for models with spherical and disk geometry. In this case both models have the same radius and central optical depth, the disk models having a 1:1 aspect ratio. While the dust temperature distributions in the two cases are very similar, the emergent flux for the disk model depends sensitively on the viewing angle. For clouds which are unresolved, one would expect, since the thermal emission is isotropic in the neighborhood of an emitting grain, and since the emission (in the far infrared) is optically thin, that the emergent flux spectrum should be characteristic only of the dust temperature, and independent of viewing angle. However, because of the lack of complete symmetry and the resulting radiation anisotropy, this is found not to be the case for the disk models. Angle is that it implies a large uncertainty in estimating the radiation dust mass.

Spagna, G. F., Jr.↗

Stellar convection. III - Convection at large Rayleigh numbers

The results of a numerical, three-dimensional study of convection in a self-gravitating sphere of Boussinesq fluid with a Rayleigh number of 10-billion and a Prandtl number of 1 are presented. The velocity and temperature are computed by using spectral methods in the horizontal and finite-differencing in the radial directions. An eddy viscosity and diffusivity are needed to model the subresolution flow. For Rayleigh numbers much less than 10-billion the flows do not have well-defined inertial ranges, and an eddy viscosity and diffusivity cannot be assigned in this self-consistent manner. By computing the energy spectra as well as the detailed energy budgets as a function of wavenumber, it is shown that for Rs = 10-billion there is an inertial range for the modes corresponding to spherical harmonics with l greater than 6.

Marcus, P. S.↗

Relativistic scattered wave calculations on UF6

Self-consistent Dirac-Slater multiple scattering calculations are presented for UF6. The results are compared critically to other relativistic calculations, showing that the results of all molecular orbital calculations are in qualitative agreement, as measured by energy levels, population analyses, and spin-orbit splittings. A detailed comparison is made to the relativistic X alpha(RX alpha) method of Wood and Boring, which also uses multiple scattering theory, but incorporates relativistic effects in a more approximate fashion. For the most part, the RX alpha results are in agreement with the present results.

Case, D. A.↗

Evaluation and optimization of silicon sheet solar cells

This paper describes the results and procedures to evaluate and improve the efficiency of solar cells made from various unconventional silicon sheets. The performance parameters included photovoltaic characteristics, spectral response, dark I-V characteristics, and diffusion length. The evaluation techniques used provided accurate and reliable information on sheet performance, and self-consistent results were obtained from the various measurement techniques used. Minority carrier diffusion length (L) was shown to be the ultimate limiting factor for the sheet cell performance (efficiency) and other back-up measurements confirmed this L-dependence. Limited efforts were made to identify defects which influence cell performance, and to use some improved process methods to increase cell efficiency.

Yoo, H.↗

Effective core potentials for the cadmium and mercury atoms

Ab initio effective core potentials have been obtained for the cadmium and mercury atoms by the methods of Kahn et al. (1976). Both two and twelve valence electron representations of Cd and Hg were tested for various atom state-configurations by comparison with all-electron calculations. The generated potentials were used to obtain the equilibrium bond distances and molecular binding energies for the dichloride and dimethyl compounds of both atoms from single and optimum-double configuration self-consistent field calculations.

Basch, H.↗

A model study of diurnal behavior of the ionosphere and the protonosphere coupling

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.

Marubashi, K.↗

MCSCF wave functions for excited states of polar molecules - Application to BeO

A previously reported multi-configuration self-consistent field (MCSCF) algorithm based on the generalized Brillouin theorem is extended in order to treat the excited states of polar molecules. In particular, the algorithm takes into account the proper treatment of nonorthogonality in the space of single excitations and invokes, when necessary, a constrained optimization procedure to prevent the variational collapse of excited states. In addition, a configuration selection scheme (suitable for use in conjunction with extended configuration interaction methods) is proposed for the MCSCF procedure. The algorithm is used to study the low-lying singlet states of BeO, a system which has not previously been studied using an MCSCF procedure. MCSCF wave functions are obtained for three 1 Sigma + and two 1 Pi states. The 1 Sigma + results are juxtaposed with comparable results for MgO in order to assess the generality of the description presented here.

Bauschlicher, C. W., Jr.↗

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

Self-consistent inclusion of space-charge in the traveling wave tube

It is shown how the complete field of the electron beam may be incorporated into the transmission line model theory of the traveling wave tube (TWT). The fact that the longitudinal component of the field due to the bunched beam is not used when formulating the beam-to-circuit coupling equation is not well-known. The fundamental partial differential equation for the traveling wave field is developed and compared with the older (now standard) one. The equation can be solved numerically using the same algorithms, but now the coefficients can be updated continuously as the calculation proceeds down the tube. The coefficients in the older equations are primarily derived from preliminary measurements and some trial and error. The newer coefficients can be found by a recursive method, since each has a well defined physical interpretation and can be calculated once a reasonable first trial solution is postulated. The results of the new expression were compared with those of the older forms, as well as to a field theory model to show the ease in which a reasonable fit to the field prediction is obtained. A complete summary of the existing transmission line modeling of the TWT is given to explain the somewhat vague ideas and techniques in the general area of drifting carrier-traveling circuit wave interactions. The basic assumptions and inconsistencies of the existing theory and areas of confusion in the general literature are examined and hopefully cleared up.

Freeman, Jon C.↗

The Third General Scientific Assembly of the International Association of Geomagnetism and Aeronomy - Special sessions of auroral processes

Methods of timing magnetic substorms, the rapid fluctuations of aurorae, electromagnetic and electrostatic instabilities observed on the field lines of aurorae, the auroral microstructure, and the relationship of currents, electric field and particle precipitation to auroral form are discussed. Attention is given to such topics as D-perturbations as an indicator of substorm onset, the role of the magnetotail in substorms, spectral information derived from imaging data on aurorae, terrestrial kilometric radiation, and the importance of the mirror force in self-consistent models of particle fluxes, currents and potentials on auroral field lines.

Russell, C. T.↗

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

Quantitative simulation of a magnetospheric substorm. I - Model logic and overview

Results of a comprehensive computer simulation of the behavior of the earth's inner magnetosphere during a substorm-type event are reported. It is pointed out that the computer model self-consistently computes electric fields, currents, and plasma distributions and velocities in the inner-magnetosphere/ionosphere system; parallel electric fields and ionospheric neutral winds, however, are not included. The basic equations of the model are derived, and the inputs are described. An overview of the results is also given. The first appendix contains derivations of general, useful laws of bounce-averaged gradient, curvature, and E x B drifts in a plasma with isotropic pitch angle distribution. The second appendix gives a description of the numerical method used in the simulation.

Harel, M.↗

Evolution of magnetic field and atmospheric response. I - Three-dimensional formulation by the method of projected characteristics. II - Formulation of proper boundary equations

The method described as the method of nearcharacteristics by Nakagawa (1980) is renamed the method of projected characteristics. Making full use of properties of the projected characteristics, a new and simpler formulation is developed. As a result, the formulation for the examination of the general three-dimensional problems is presented. It is noted that since in practice numerical solutions must be obtained, the final formulation is given in the form of difference equations. The possibility of including effects of viscous and ohmic dissipations in the formulation is considered, and the physical interpretation is discussed. A systematic manner is then presented for deriving physically self-consistent, time-dependent boundary equations for MHD initial boundary problems. It is demonstrated that the full use of the compatibility equations (differential equations relating variations at two spatial locations and times) is required in determining the time-dependent boundary conditions. In order to provide a clear physical picture as an example, the evolution of axisymmetric global magnetic field by photospheric differential rotation is considered.

Nakagawa, Y.↗

N-body simulations of disks

The methods used in large-scale n-body simulations are discussed. However, the present review concentrates on results already obtained in n-body simulations using systems containing up to 200,000 simulation stars. Results are presented which show that the stability criterion developed for flattened systems applies only to truly axisymmetric instabilities. Purely stellar disks acquire rather large velocity dispersions, generally two or more times the velocity dispersion required by Toomre (1964) for axisymmetric stability. In computer simulations, the bar-forming instability can be prevented only by comparatively large velocity dispersions. However, simulations including the effects of the galactic halo and core as a fixed background field show that bar formation can be prevented for fixed halo components as large or larger than the self-consistent disk component. Experiments performed to determine the collisional relaxation time for large-scale gravitational n-body calculations show that these models are indeed 'collisionless'.

Hohl, F.↗

A modified Monte Carlo model for the ionospheric heating rates

A Monte Carlo method is adopted as a basis for the derivation of the photoelectron heat input into the ionospheric plasma. This approach is modified in an attempt to minimize the computation time. The heat input distributions are computed for arbitrarily small source elements that are spaced at distances apart corresponding to the photoelectron dissipation range. By means of a nonlinear interpolation procedure their individual heating rate distributions are utilized to produce synthetic ones that fill the gaps between the Monte Carlo generated distributions. By varying these gaps and the corresponding number of Monte Carlo runs the accuracy of the results is tested to verify the validity of this procedure. It is concluded that this model can reduce the computation time by more than a factor of three, thus improving the feasibility of including Monte Carlo calculations in self-consistent ionosphere models.

Mayr, H. G.↗

Flare energetics

In this investigation of flare energetics, researchers sought to establish a comprehensive and self-consistent picture of the sources and transport of energy within a flare. To achieve this goal, they chose five flares in 1980 that were well observed with instruments on the Solar Maximum Mission, and with other space-borne and ground-based instruments. The events were chosen to represent various types of flares. Details of the observations available for them and the corresponding physical parameters derived from these data are presented. The flares were studied from two perspectives, the impulsive and gradual phases, and then the results were compared to obtain the overall picture of the energics of these flares. The role that modeling can play in estimating the total energy of a flare when the observationally determined parameters are used as the input to a numerical model is discussed. Finally, a critique of the current understanding of flare energetics and the methods used to determine various energetics terms is outlined, and possible future directions of research in this area are suggested.

Wu, S. T.↗

A modified Monte Carlo model for the ionospheric heating rates.

A Monte Carlo method is adopted as a basis for the derivation of the photoelectron-heat input into the ionospheric plasma. Since a great number of Monte Carlo runs are required normally for the computation of the heating rates, this approach is modified in an attempt to minimize the computation time. The heat-input distributions are computed for arbitrarily small source elements that are spaced apart at distances corresponding to the photoelectron dissipation range. By means of a nonlinear interpolation procedure their individual heating-rate distributions are utilized to produce synthetic ones that fill the gaps between the Monte Carlo generated distributions. By varying these gaps and the corresponding number of Monte Carlo runs the accuracy of the results is tested to verify the validity of this procedure. It is concluded that this model can reduce the computation time by as much as an order of magnitude, thus improving the feasibility of including Monte Carlo calculations in self-consistent ionosphere models.

Mayr, H. G.↗

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗