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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Crashworthiness of light aircraft fuselage structures: A numerical and experimental investigation

The dynamic behavior of aircraft fuselage structures subject to various impact conditions was investigated. An analytical model was developed based on a self-consistent finite element (CFE) formulation utilizing shell, curved beam, and stringer type elements. Equations of motion were formulated and linearized (i.e., for small displacements), although material nonlinearity was retained to treat local plastic deformation. The equations were solved using the implicit Newmark-Beta method with a frontal solver routine. Stiffened aluminum fuselage models were also tested in free flight using the UTIAS pendulum crash test facility. Data were obtained on dynamic strains, g-loads, and transient deformations (using high speed photography in the latter case) during the impact process. Correlations between tests and predicted results are presented, together with computer graphics, based on the CFE model. These results include level and oblique angle impacts as well as the free-flight crash test. Comparisons with a hybrid, lumped mass finite element computer model demonstrate that the CFE formulation provides the test overall agreement with impact test data for comparable computing costs.

Nanyaro, A. P.↗
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 NTRS1976

Computation on collisionless steady-state plasma flow past a charged disk

A computer method is presented using the 'inside-out' approach, for predicting the structure of the disturbed zone near a moving body in space. The approach uses fewer simplifying assumptions than other available methods, and is applicable to large ranges of the values of body and plasma parameters. Two major advances concerning 3-dimensional bodies are that thermal motions of ions as well as of electrons are treated realistically by following their trajectories in the electric field, and the technique for achieving self-consistency is promising for very large bodies. Three sample solutions were obtained for a disk-shaped body, charged negatively to a potential 4kT/e. With ion Mach number 4, and equal ion and electron temperatures, the wakes of a relatively small body (radius 5 Debye lengths) and a relatively large body (radius 100 Debye lengths) both begin to fill up between 2 and 3 body radii downstream. For the large body there is in addition a potential well (about 6kT/e deep) behind the body. Increasing the ion Mach number to 8 for the large body causes the potential well to become wider and longer but not deeper. For the large body, the quasineutrality assumption is validated outside of a cone-shaped region in the very near wake. For the large as well as the small body, the disturbed zone behind the body extends transversely no more than 2 or 3 body radii, a result of significance for the design of spacecraft boom instrumentation.

Parker, L. W.↗
NASA NTRS1980

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

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

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

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

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

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

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

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

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

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

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

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

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