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.
A procedure to calculate all second-order polarizability tensor elements in static fields has been developed. The calculations are based on semiempirical Hamiltonians (MNDO) that include shaped electric fields. The technique has been applied to mono-, di-, and trisubstituted benzenes incorporating nitro, methyl, and primary and secondary amino groups.
The strain-induced splitting of the heavy-hole (hh) and light-hole (lh) valence bands for 4-microns thick GaAs/Si is examined on a microscopic scale using linear polarized-cathodoluminescence imaging and spectroscopy. The energies and intensities of the hh- and lh-exciton luminescence are quantitatively analyzed to determine spatial variations in the stress tensor. The results indicate that regions near and far from the microcracks are primarily subject to uniaxial and biaxial tensile stresses, respectively. The transition region where biaxial stress gradually converts to uniaxial stress is analyzed, and reveals a mixing of lh and hh characters in the strain-split bands.
The fine-scale structure of turbulence in a fully developed turbulent duct flow is examined by considering the 3D velocity derivative field obtained from direct numerical simulations at two relatively small Reynolds numbers. The magnitudes of all mean-square derivatives (normalized by wall variables) increase with the Reynolds number, the increase being largest at the wall. These magnitudes are not consistent with the assumption of local isotropy except perhaps near the duct center-line. When the assumption of local isotropy is relaxed to one of local axisymmetry, or invariance with respect to rotation about a coordinate axis (here chosen in the streamwise direction), satisfactory agreement is indicated by the data outside the wall region. Support for axisymmetry is demonstrated by anisotropy invariant maps of the dissipation and vorticity tensors.
The turbulent mixing noise of a supersonic jet is calculated for a round convergent-divergent nozzle at the design pressure ratio. Aerodynamic computations are performed using the PARC code with a k-epsilon turbulence model. Lighthill's acoustic analogy combined with Ribner's assumption is adopted. The acoustics solution is based upon the methodology followed by GE in the MGB code. The source correlation function is expressed as a linear combination of second-order tensors. Assuming separable second-order correlations and incorporating Batchelor's isotropic turbulence model, the source term was calculated from the kinetic energy of turbulence. A Gaussian distribution for the time-delay of correlation was introduced. The computational fluid dynamics (CFD) solution was used to obtain the source strength as well as the characteristic time-delay of correlation. The effect of sound/flow interaction was incorporated using the high frequency asymptotic solution to Lilley's equation for axisymmetric geometries. Acoustic results include sound pressure level directivity and spectra at different polar angles. The aerodynamic and acoustic results demonstrate favorable agreement with experimental data.
A generalized curvilinear coordinate Taylor weak statement implicit finite element algorithm is developed for the two-dimensional and axisymmetric compressible Navier-Stokes equations for ideal and reacting gases. For accurate hypersonic simulation, air is modeled as a mixture of five perfect gases, i.e., molecular and atomic oxygen and nitrogen as well as nitric oxide. The associated pressure is then determined via Newton solution of the classical chemical equilibrium equation system. The directional semidiscretization is achieved using an optimal metric data Galerkin finite element weak statement, on a developed 'companion conservation law system', permitting classical test and trial space definitions. Utilizing an implicit Runge-Kutta scheme, the terminal algorithm is then nonlinearly stable, and second-order accurate in space and time on arbitrary curvilinear coordinates. Subsequently, a matrix tensor product factorization procedure permits an efficient numerical linear algebra handling for large Courant numbers. For ideal- and real-gas hypersonic flows, the algorithm generates essentially nonoscillatory numerical solutions in the presence of strong detached shocks and boundary layer-inviscid flow interactions.
The equilibrium structure of homogeneous turbulent shear flow is investigated from a theoretical standpoint. Existing turbulence models, in apparent agreement with physical and numerical experiments, predict an unbounded exponential time growth of the turbulent kinetic energy and dissipation rate; only the anisotropy tensor and turbulent time scale reach a structural equilibrium. It is shown that if a residual vortex stretching term is maintained in the dissipation rate transport equation, then there can exist equilibrium solutions, with bounded energy states, where the turbulence production is balanced by its dissipation. Illustrative calculations are presented for a k-epsilon model modified to account for net vortex stretching.
The effect of including the Braginskii viscous stress tensor in magnetohydrodynamics is remarked upon. It is shown that semiquantitative agreement with a recently observed anisotropy in the turbulent solar wind spectrum can be achieved in this way. The modifications of the dynamical equations are simple enough to permit their inclusion in numerical codes. The effects of large 'ion parallel viscosities' also may be significant for plasmas in quite different regimes than the solar wind.
The time-dependent Navier-Stokes equations are numerically integrated for two-dimensional incompressible viscous flow in a shear-driven square cavity. Using a time-splitting method and finite differences on a staggered mesh, the momentum and pressure equations are directly solved by a tensor product method where one finite difference direction is diagonalized by eigenvalue decomposition. The effects of increasing Reynolds number are studied and the developing boundary layer is captured by using a finely clustered mesh. At Re = 30,000 the flow is in a continuously developing unsteady regime. Power spectrum plots indicate that the unsteady flow oscillates with one fundamental frequency and exhibits some characteristics of transition between laminar and turbulent states.
Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.
Methods employed to derive recent kinetic theories for rapid noncomminuting granular flows are extended to homogeneous flows in which a fraction of the repeated collisions produce tiny fractures on the particles' peripheries and gradually reduce their effective diameters. The theory consists of balance equations for mass, momentum, and energy, as well as constitutive relations for the presence tensor and collisional rates of mass and energy lost. The work of Richman and Chou (1989) is improved by incorporating into the constitutive theory the critical impact energy below which no mass loss occurs in a binary collision. The theory is applied to granular shear flows and, for fixed shear rates, predicts the time variations of the solid fraction granular temperature, and induced stresses, as well as their extreme sensitivities to small changes in the critical impact energy.
Results of three-dimensional MHD simulations of magnetotail dynamics with anisotropic pressure are presented. The pressure tensor is assumed to be gyrotropic, satisfying a modified double adiabatic approximation including Ohmic heating. It is found that these constraints tend to stabilize the tail with respect to the resistive tearing instability. Including an increasing level of additional isotropization results in a gradual transition to the fast instability found in an isotropic model, which is discussed for comparison. Possible isotropization mechanisms yielding magnetotail instability are discussed also. The stabilizing anisotropies are found primarily in the plasma sheet boundary region. This fact indicates a possibly important role of isotropization mechanisms operating in this region in destabilizing the magnetotail and initiating a tearing instability without the necessity of changes in fluctuation levels in the neutral sheet.
Report describes experiments in which linearly polarized cathodoluminesce scanning electron microscopy, in both imaging and spectroscopic forms, used to analyze spatial variations in stress tensors in specimens of GaAs/Si. Specimens fabricated by atmospheric-pressure metallo-organic chemical vapor deposition of GaAs to thickness of 4 micrometers on Si substrates. Combination of mismatch between Si and GaAs crystal lattices and differential thermal contraction of materials from deposition temperature (700 degrees C) to room temperature gives rise to stresses in GaAs layer. Experiments conducted to contribute to understanding of effects of these stresses and ramifications thereof for fabrication of electronic devices in GaAs/Si system.
New formulation for Nystrom integrators offers increased precision and speed of computation. Formulation takes advantage of capabilities of modern electronic computers. Independent and dependent variables and various partial and total derivatives expressed via compact tensor notation. In this approach, typical set of constraint equations derived in about 1 day. Solution of constraint equations speeded by SEARCH computer program. Is general-purpose optimization program, which maximizes or minimizes cost function defined by user. Used successfully, to make least-squares fit of about 4,000 equations in 12 unknowns.
The flow field inside a whirling annular seal was measured using a 3-D Laser Doppler Anemometer (LDA) system. The seal investigated has a clearance of 1.27 mm, a length of 37.3 mm, and is mounted on a drive shaft with a 50 percent eccentricity ratio. This results in the rotor whirling at the same speed as the shaft rotation (whirl ratio = 1.0). The seal is operated at Reynolds number of 12,000 and a Taylor number of 6,300 (3,600 rpm). The 3-D LDA system is equipped with a rotary encoding system which is used to produce phase averaged measurements of the entire mean velocity vector field and Reynolds stress tensor field from 0.13 mm upstream to 0.13 mm downstream of the seal. The mean velocity field reveals a highly three dimensional flow field with large radial velocities near the inlet of the seal as well as a recirculation zone on the rotor surface. The location of maximum mean axial velocity migrates from the pressure side of the rotor at the inlet to the suction side at turbulence kinetic energy. However, turbulence production and dissipation attain equilibrium fairly quickly with remaining relatively constant over the last half of the seal.
The mean velocity and Reynolds stress tensor throughout a whirling annular seal are presented. The data was collected with a three dimensional laser Doppler velocimeter using phase averaging. Two axial flow conditions (Re = 12,000 and 24,000) were studied at one shaft speed (Ta = 6,600). The eccentricity and whirl ratios were 50 and 100 percent, respectively. There is a region of high axial momentum in this region is higher in the low Reynolds number case due to an axial recirculation zone that occurs on the suction side of the rotor at the inlet. The recirculation zone does not occur in the high Reynolds number case. At both Reynolds numbers, there is a recirculation zone on the rotor surface in the pressure side of the inlet. This recirculation zone extends from 20 to 200 degrees rotor zenith in the tangential direction, and is one third of a clearance wide radially. The high Reynolds number recirculation zone is 1.5 mean clearances long, while the low Reynolds number zone extends 2 mean clearances downstream. When compared to previous studies, it is apparent that the tangential momentum is no greater for a seal with whirl than for one without if other parameters are constant. Areas of high tangential momentum occur in the clearance where the axial momentum is low. Average exit plane tangential velocities in the high Reynolds number case are 1.5 times greater than those in the other flow case. These results are in general agreement with predictions made by other investigators.
The shape of seven large domes on the plains of Venus, with volumes between 100 and 1000 cu km, is compared with that of an axisymmetric gravity current spreading over a rigid horizontal surface. Both the altimetric profiles and the horizontal projection of the line of intersection of domes on the SAR images agree well with the theoretical similarity solution for a newtonian fluid, but not with the shape calculated for a rigid-plastic rheology, nor with that for a static model with a strong skin. As a viscous current spreads, it generates an isotropic strain rate tensor whose magnitude is independent of radius. Such a flow can account for the randomly oriented cracks that are uniformly distributed on the surface of the domes. The stress induced by the flow in the plains material below is obtained, and is probably large enough to produce the short radial cracks in the surface of the plains beyond the domes. The viscosity of the domes can be estimated from their thermal time constants if spreading is possible only when the fluid is hot, and lies between 10(exp 14) and 10(exp 17) Pa s. Laboratory experiments show that such viscosities correspond to temperatures of 610 - 690 C in dry rhyolitic magmas. These temperatures agree with laboratory measurements of the solidus temperature of wet rhyolite. These results show that the development of the domes can be understood using simple fluid dynamical ideas, and that the magmas involved can be produced by wet melting at depths below 10 km, followed by eruption and degassing.
A new and comprehensive dielectric tensor characterization instrument is presented for characterization of magneto-optical recording media and non-magnetic thin films. Random and systematic errors of the system are studied. A series of TbFe, TbFeCo, and Co/Pt samples with different composition and thicknesses are characterized for their optical and magneto-optical properties. The optical properties of several non-magnetic films are also measured.