SEARCH · Engineering Papers
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.
Vacuum polarization in the field of a high-Z finite-size nucleus is examined, and the polarization charge density in coordinate space of order alpha(Z alpha)exp n not less than 3 is calculated. Energy-level shifts of K- and L-shell electrons in hydrogenlike systems are given.
An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented which provides high flowfield resolution when applied to two-dimensional and three-dimensional flow problems. The present scheme locally adjusts the grid density using gradients of selected flow variables from either finite difference or finite volume calculations, and it can specify user-prescribed grid stretching such that control of the grid spacing can be maintained in areas of known flowfield behavior. A robust and efficient grid smoothing technique has been incorporated into the adaptive grid routine.
This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.
Scintillations (high frequency variations) observed in the radio signal during the occultation of Voyager 1 by Titan (Hinson and Tyler, 1983) provide information concerning neutral atmospheric density fluctuations on scales on hundreds of meters to a few kilometers. Those seen at altitudes higher than 25 km above the surface were interpreted by Hinson and Tyler as being caused by linear, freely propagating (energy-conserving) gravity waves, but this interpretation was found to be inconsistent with the scintillation data below the 25-km altitude level. Here an attempt is made to interpret the entire scintillation profile between the surface and the 90-km altitude level in terms of gravity waves generated at the surface. Numerical calculations of the density fluctuations caused by two-dimensional, nonhydrostatic, finite-amplitude gravity waves propagating vertically through Titan's atmosphere are performed to produce synthetic scintillation profiles for comparison with the observations. The numerical model accurately treats the effects of wave transience, nonlinearity, and breakdown due to convective instability in the overturned part of the wave. The high-altitude scintillation data were accurately recovered with a freely propagating wave solution, confirming the analytic model of Hinson and Tyler. It is found that the low-altitude scintillation data can be fit by a model where a component of the gravity waves becomes convectively unstable and breaks near the 15 km level. The large-scale structure of the observed scintillation profile in the entire altitude range between 5 and 85 km can be simulated by a model where the freely propagating and breaking waves are forced at the surface simultaneously. Further analysis of the Voyager 1 Titan low-altitude scintillation data, using inversion theory appropriate for strong scattering, could potentially remove some of the ambiguities remaining in this analysis and allow a better determination of the strength and source of the waves.
The Monte Carlo simulation method in conjunction with the finite element large deflection modal formulation are used to estimate fatigue life of aircraft panels subjected to stationary Gaussian band-limited white-noise excitations. Ten loading cases varying from 106 dB to 160 dB OASPL with bandwidth 1024 Hz are considered. For each load case, response statistics are obtained from an ensemble of 10 response time histories. The finite element nonlinear modal procedure yields time histories, probability density functions (PDF), power spectral densities and higher statistical moments of the maximum deflection and stress/strain. The method of moments of PSD with Dirlik's approach is employed to estimate the panel fatigue life.
The inductive formation of current sheets in a conical theta pinch FARAD (Faraday Accelerator with Radio-frequency Assisted Discharge) thruster is investigated experimentally with time-integrated photography. The goal is to help in understanding the mechanisms and conditions controlling the strength and extent of the current sheet, which are two indices important for FARAD as a propulsion concept. The profiles of these two indices along the inside walls of the conical acceleration coil are assumed to be related to the profiles of the strength and extent of the luminosity pattern derived from photographs of the discharge. The variations of these profiles as a function of uniform back-fill neutral pressure (with no background magnetic field and all parameters held constant) provided the first clues on the nature and qualitative dependencies of current sheet formation. It was found that there is an optimal pressure for which both indices reach a maximum and that the rate of change in these indices with pressure differs on either side of this optimal pressure. This allowed the inference that current sheet formation follows a Townsend-like breakdown mechanism modified by the existence of a finite pressure-dependent radio-frequency-generated electron density background. The observation that the effective location of the luminosity pattern favors the exit-half of the conical coil is explained as the result of the tendency of the inductive discharge circuit to operate near its minimal self-inductance. Movement of the peak in the luminosity pattern towards the upstream side of the cone with increasing pressure is believed to result from the need of the circuit to compensate for the increase in background plasma resistivity due to increasing pressure.
A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.
Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.
The approach solves for both Reynolds and Favre averaged quantities and calculates the scalar pdf. Turbulent models used to close the governing equations are formulated to account for complex mixing and variable density effects. In addition, turbulent mass diffusivities are not assumed to be in constant proportion to turbulent momentum diffusivities. The governing equations are solved by a combination of finite-difference technique and Monte-Carlo simulation. Some preliminary results on simple variable density shear flows are presented. The differences between these results and those obtained using conventional models are discussed.
A quantile-based approach to statistical analysis and probability modeling of data is presented which formulates statistical inference problems as functional inference problems in which the parameters to be estimated are density functions. Density estimators can be non-parametric (computed independently of model identified) or parametric-select (approximated by finite parametric models that can provide standard models whose fit can be tested). Exponential models and autoregressive models are approximating densities which can be justified as maximum entropy for respectively the entropy of a probability density and the entropy of a quantile density. Applications of these ideas are outlined to the problems of modeling: (1) univariate data; (2) bivariate data and tests for independence; and (3) two samples and likelihood ratios. It is proposed that bi-information estimation of a density function can be developed by analogy to the problem of identification of regression models.
Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.
Ion bombardment of a growing film during deposition is necessary in many instances to ensure a fully dense coating, particularly for hard coatings. Until the recent advent of unbalanced magnetron (UBM) cathodes, reactive sputtering had not been able to achieve the same degree of ion bombardment as other physical vapor deposition processes. The amount of ion bombardment of the substrate depends on the plasma density at the substrate, and in a UBM system the amount of bombardment will depend on the degree of unbalance of the cathode. In multi-cathode systems, the magnetic fields between the cathodes must be linked to confine the fast electrons that collide with the gas atoms. Any break in this linkage results in electrons being lost and a low plasma density. Modeling of the magnetic fields in a UBM cathode using a finite element analysis program has provided great insight into the interaction between the magnetic fields in multi-cathode systems. Large multi-cathode systems will require very strong magnets or many cathodes in order to maintain the magnetic field strength needed to achieve a high plasma density. Electromagnets offer the possibility of independent control of the plasma density. Such a system would be a large-scale version of an ion beam enhanced deposition (IBED) system, but, for the UBM system where the plasma would completely surround the substrate, the acronym IBED might now stand for Ion Blanket Enhanced Deposition.
A general expression for the ponderomotive force of electromagnetic field with slowly varying amplitude in a warm plasma is obtained in the collisionless two-fluid model. Compared to the existing expression obtained for a cold and stationary plasma, this result contains two additional terms for the general case. One of them is involved with the density gradient and the wavenumber derivative of the dielectric tensor of a plasma. Thus, for an inhomogeneous plasma, it accounts for the effects of spatial dispersion, hence of finite temperature of the plasma. The other term is proportional to the time derivative of the density and the frequency derivative of the dielectric tensor, and is nonzero unless the plasma is stationary or nondispersive. Inclusion of these new terms in the ponderomotive force will enable self-consistent analysis of the general behavior of nonlinear waves with finite wavelengths in a warm plasma.
The analysis of plasma beam and shear problems in magnetic fields is usually based on a hydromagnetic fluid model. In a low-density collisionless plasma, however, the kinetic effects of the plasma, such as finite Larmor radius effects, are not yet clearly understood. In this paper, the kinetic equilibria of plasma shears in a uniform and fixed magnetic field, with full ion motion, are discussed by solving the Vlasov equation with a given electric field and drift velocity. In this model, the ion density profile through the plasma shear layer is quite different from the one predicted by a hydromagnetic model. As a result of a complicated ion gyromotion through the shear layer, single- and double-humped ion density profiles are obtained. The dependence on the temperature and the strength of the shear will be discussed. The results show a significant difference between positive and negative shears.
A study was performed to examine the influence of varying mesh density on an LS-DYNA simulation of a rectangular-shaped foam projectile impacting the space shuttle leading edge Panel 6. The shuttle leading-edge panels are fabricated of reinforced carbon-carbon (RCC) material. During the study, nine cases were executed with all possible combinations of coarse, baseline, and fine meshes of the foam and panel. For each simulation, the same material properties and impact conditions were specified and only the mesh density was varied. In the baseline model, the shell elements representing the RCC panel are approximately 0.2-in. on edge, whereas the foam elements are about 0.5-in. on edge. The element nominal edge-length for the baseline panel was halved to create a fine panel (0.1-in. edge length) mesh and doubled to create a coarse panel (0.4-in. edge length) mesh. In addition, the element nominal edge-length of the baseline foam projectile was halved (0.25-in. edge length) to create a fine foam mesh and doubled (1.0-in. edge length) to create a coarse foam mesh. The initial impact velocity of the foam was 775 ft/s. The simulations were executed in LS-DYNA for 6 ms of simulation time. Contour plots of resultant panel displacement and effective stress in the foam were compared at four discrete time intervals. Also, time-history responses of internal and kinetic energy of the panel, kinetic and hourglass energy of the foam, and resultant contact force were plotted to determine the influence of mesh density.
The structure of the drift instability was examined in several density regimes. Let sub e be the total electron mean free path, k sub z the wave-vector component along the magnetic field, and the ratio of perpendicular ion diffusion to parallel electron streaming rates. At low densities (k sub z lambda 1) the drift mode is isothermal and should be treated kineticly. In the finite heat conduction regime square root of m/M k sub z Lambda sub 1) the drift instability threshold is reduced at low densities and increased at high densities as compared to the isothermal threshold. Finally, in the energy transfer limit (k sub z kambda sub e square root of m/M) the drift instability behaves adiabatically in a fully ionized plasma and isothermally in a partially ionized plasma for an ion-neutral to Coulomb collision frequency ratio.
This report summarizes work that has been performed to date on the strength of a cryotank insulation system using Rohacell foam and TUFI-coated AETB-12 ceramic tiles directly bonded to a simulated graphite-epoxy tank wall. Testing utilized a custom specimen design which consists of a long tensile specimen with eccentric loading to induce curvature similar to the curvature expected due to 'pillowing' of the tank when pressurized. A finite element model was constructed to predict the specific element strains in the test article, and to assist with design of the test specimen to meet the specific goals of curvature and laminate strain. The results indicate that the heat treated 3.25-pcf density Rohacell foam does not provide sufficient strength for the induced stresses due to curvature and stress concentration at the RTV bondline to the TUFI tile. The test was repeated using higher density non-heat treated Rohacell foam (6.9 pcf) without foam failure. The finite element model was shown to predict specimen behavior, and validation of the model was successful. It is pertinent to mention that the analyses described herein accurately predicted the failure of the heat treated foams and based on this analysis method it is expected that the untreated 3.25 pcf Rohacell foam will be successful.
Density matrices of symmetry projected single determinant wave functions for finite groups, considering many particle system