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.
Method of solving Navier-Stokes equations of incompressible flow in general nonorthogonal curvilinear coordinates incorporates fractional-step and finite-volume approaches. Developed in continuing effort to achieve accuracy without need for excessive computation time in numerical simulation of time-dependent, three-dimensional flows bounded by surfaces of complicated shape.
We have applied and compared the efficiency and accuracy of two commonly used numerical methods for the solution of Navier-Stokes equations. The artificial compressibility method augments the continuity equation with a transient pressure term and allows one to solve the modified equations as a coupled system. Due to its implicit nature, one can have the luxury of taking a large temporal integration step at the expense of higher memory requirement and larger operation counts per step. Meanwhile, the fractional step method splits the Navier-Stokes equations into a sequence of differential operators and integrates them in multiple steps. The memory requirement and operation count per time step are low, however, the restriction on the size of time marching step is more severe. To explore the strengths and weaknesses of these two methods, we used them for the computation of a two-dimensional driven cavity flow with Reynolds number of 100 and 1000, respectively. Three grid sizes, 41 x 41, 81 x 81, and 161 x 161 were used. The computations were considered after the L2-norm of the change of the dependent variables in two consecutive time steps has fallen below 10(exp -5).
An enhancement of the crack diffusion model of brittle fracture is proposed. In the model, certain probabilities are evaluated as averages over the set of all virtual crack paths, and the latter are modeled as realization of a Wiener process. However, commonly observed fracture profiles are smoother. In fact, one can show that their fractal dimension is smaller than the value 3/2 for a Wiener process. Accordingly, it is proposed to represent fracture profiles by graphs of fractional integrals of a Wiener process. This allows one to employ additional fractographic information obtained from fracture profiles, e.g. their fractal dimension, in a predictive formalism, evaluating probabilities of failure.
The lubrication problem of infinitely long slider bearings with a mixture of fluid and particulate solid at small volume fraction level is studied. Closed-form analytical solutions for pressure and shear stress are obtained for a class of solid aggregates. The results reduce to those of pure fluid in the limiting case. A parametric study of the bearing performance with particulate solid is presented.
A bulk PBL model with a simple internal vertical structure and a simple second-order closure has been developed for use as a PBL parameterization in a large-scale model. A method to determine the value of sigma, the fractional area covered by rising motion, is given. Large-eddy simulations of Moeng (1984, 1986) are used to validate various aspects of the parameterization.
The efficiency of the INS3D-fractional step (FS) code has been significantly enhanced by accelerating the Poisson solver with a multigrid (MG) method. This is apparently the first implementation of the MG method for solving the 3D discrete Poisson-like equation obtained in the derivation of FS solution methods for the incompressible Navier-Stokes equations in generalized nonorthogonal coordinate systems using a staggered arrangement of the variables. The MG solver is insensitive to the geometry, even in cases of highly nonorthogonal clustered meshes, as well as to the free parameters of MG methods.
A preliminary presentation is given of a model which explains the picritic lunar mare green glass compositions as composites or magmas derived by polybaric fractional melting. The model accommodates changes in source composition and solidus temperature. As in the Klein and Langmuir model, the final compositions of the melts signify only average depths of melting. The onset of melting is deeper, final segregation is shallower.
The order reduction problem of z-transfer functions is solved by using the multipoint Jordan continued-fraction expansion (MJCFE) technique. An efficient algorithm that does not require the use of complex algebra is presented for obtaining an MJCFE from a stable z-transfer function with expansion points selected from the unit circle and/or the positive real axis of the z-plane. The reduced-order models are exactly the multipoint Pade approximants of the original system and, therefore, they match the (weighted) time-moments of the impulse response and preserve the frequency responses of the system at some characteristic frequencies, such as gain crossover frequency, phase crossover frequency, bandwidth, etc.
Ruthenium (Ru) and iridium (Ir) are the least mobile platinum group elements (PGE's) within the Cretaceous-Tertiary (K-T) boundary clay (BC). The Ru/Ir ratio is, therefore, the most useful PGE interelement ratio for distinguishing terrestrial and extraterrestrial contributions to the BC. The Ru/Ir ratio of marine K-T sections (1.77 +/- 0.53) is statistically different from that of the continental sections (0.93 +/- 0.28). The marine Ru/Ir ratios are chondritic (C1 = 1.48 +/- 0.09), but the continental ratios are not. We discovered an inverse correlation of shocked quartz size (or distance from the impact site) and Ru/Ir ratio. This correlation may arise from the difference in Ru and Ir vaporization temperature and/or fractionation during condensation from the ejecta cloud. Postsedimentary alteration, remobilization, or terrestrial PGE input may be responsible for the Ru/Ir ratio variations within the groups of marine and continental sites studied. The marine ratios could also be attained if approximately 15 percent of the boundary metals were contributed by Deccan Trap emissions. However, volcanic emissions could not have been the principal source of the PGE's in the BC because mantle PGE ratios and abundances are inconsistent with those measured in the clay. The Ru/Ir values for pristine Tertiary mantle xenoliths (2.6 +/- 0.48), picrites (4.1 +/- 1.8), and Deccan Trap basalt (3.42 +/- 1.96) are all statistically distinct from those measured in the K-T BC.
Tables and figures are presented which show local site observations of cloud fractions, the number of cloud layers, direction of movement, and precipitation data collected during the FIRE (First ISCCP Regional Experiment) Phase 2 Cirrus Intensive Field Observations (IFO) conducted in Coffeyville, Kansas during November and December, 1991. Selected data are also presented at the times of the TIROS Operational Vertical Sounder (TOVS) satellite overpass. Several major scientific projects have used surface-based observations of clouds to compare directly with those being observed from satellites. Characterizing the physical properties of clouds is extremely useful in obtaining a more accurate analysis of the effect of clouds and their movements on weather and climate. It is the purpose of this paper to report data collected during the FIRE Phase 2 IFO experiment and to provide a brief history of such a surface-based system and the technical information required for recording local cloud parameters.
A chemical model for simulating the sources of the lunar mare basalts was developed by considering a modified mafic cumulate source formed during the combined equilibrium and fractional crystallization of a lunar magma ocean (LMO). The parameters which influence the initial LMO and its subsequent crystallization are examined, and both trace and major elements are modeled. It is shown that major elements tightly constrain the composition of mare basalt sources and the pathways to their creation. The ability of this LMO model to generate viable mare basalt source regions was tested through a case study involving the high-Ti basalts.
Results are presented from the analyses of small ballistically dispersed melt samples (occurring in the form of aerodynamically shaped spheres, dumbbells, and teardrops) from the Wabar Crater (Saudi Arabia) and of melts from the Wabar and Nejed meteorites. On the basis of the data obtained and of models for crater forming processes it is concluded that the ballistically-dispersed glasses were formed from material at shallow levels in the target zone where temperatures were higher and where intimate physical contact with the disrupting impactor was achieved. Siderophile element fractionation was completed early, while the impactor was undergoing decompression and prior to mixing with the target glasses.
The fraction of photosynthetically active radiation intercepted by vegetation, F(sub ipar) is an important parameter for modeling the interactions between the land-surface and atmosphere and for estimating vegetation biomass productivity. This study was, therefore, an integral part of FIFE. The specific purpose of this experiment was to find out how well definitive measurements of F(sub ipar) on the ground relate to near-ground and satellite based spectral reflectance measurements. Concurrent measurements of F(sub ipar) and ground, helicopter, and satellite based reflectance measurements were taken at thirteen tall-grass prairie sites within the FIFE experimental area. The sites were subjected to various combinations of burning and grazing managements. The ground and helicopter based reflectance measurements were taken on the same day or few days from the time of the overpass of LANDSAT and SPOT satellites. Ground-based reflectance measurements and sun photometer readings taken at the times of the satellite overpasses were used to correct for atmospheric attenuation. Hand-held radiometer spectral indices were strongly correlated with helicopter and satellite based values (r = 0.94 for helicopter, 0.93 for LANDSAT Thematic Mapper, and 0.86 for SPOT). However, the ground, helicopter, and satellite based normalized difference spectral vegetation indices showed low sensitivity to changes in F(sub ipar). Reflectance measurements were only moderately well correlated with measurements of F(sub ipar) (r = 0.82 for hand-held radiometer, 0.84 for helicopter measurements, and 0.75 for the LANDSAT Thematic Mapper and SPOT). Improved spectral indices which can compensate for site differences are needed in order to monitor F(sub ipar) more reliably.
The collisional ejection process for hydrogen on Venus is reanalyzed. Improved values for the efficiency of H and D escape as a function of the ionospheric temperature are reported. It is proposed that the reduction of the hydrogen flux for collisional ejection be reduced from 8 to 3.5 x 10 exp 6/sq cm/s, and a revised D/H fractional factor of 0.47 due to collisional ejection is suggested. The resulting deuterium flux is 3.1 x 10 exp 4/sq cm/s, roughly six times the flux due to charge exchange, making collisional ejection the dominant escape mechanism for deuterium on Venus.
This paper presents the use of matrix fraction descriptions for least-squares curve fitting of the frequency spectra to compute two matrix polynomials. The matrix polynomials are intermediate step to obtain a linearized representation of the experimental transfer function. Two approaches are presented: first, the matrix polynomials are identified using an estimated transfer function; second, the matrix polynomials are identified directly from the cross/auto spectra of the input and output signals. A set of Markov parameters are computed from the polynomials and subsequently realization theory is used to recover a minimum order state space model. Unevenly spaced frequency response functions may be used. Results from a simple numerical example and an experiment are discussed to highlight some of the important aspect of the algorithm.
Results of an investigation of the crystallization history of the asteroidal core that produced nearly 200 iron meteorites of group IIIAB are presented. The accuracy of the published distribution of coefficients between solid and liquid metal for seven elements and their dependence on the concentrations of S and P are critically assessed to develop a nonideal fractional crystallization model that reproduces all of the major features of the overall chemical trends in group IIIAB. It is suggested that the apparent distribution coefficient of S was was much higher than its equilibrium value (not greater than 0.01) and that it increased during crystallization. In the present preferred model the apparent distribution coefficient for S increases from 0.6 to 0.8, while the S concentration of the liquid increases from 6 to 13 wt pct. It is inferred that light S-rich liquid accumulates preferentially at the top of the core and in structural traps formed by the advancing solid.
Mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs) are derived by partial melting of the upper mantle and are marked by systematic excesses of thorium-230 activity relative to the activity of its parent, uranium-238. Experimental measurements of the distribution of thorium and uranium between the melt and solid residue show that, of the major phases in the upper mantle, only garnet will retain uranium over thorium. This sense of fractionation, which is opposite to that caused by clinopyroxene-melt partitioning, is consistent with the thorium-230 excesses observed in young oceanic basalts. Thus, both MORBs and OIBs must begin partial melting in the garnet stability field or below about 70 kilometers. A calculation shows that the thorium-230-uranium-238 disequilibrium in MORBs can be attributed to dynamic partial melting beginning at 80 kilometers with a melt porosity of 0.2 percent or more. This result requires that melting beneath ridges occurs in a wide region and that the magma rises to the surface at a velocity of at least 0.9 meter per year.
It was shown that in diamonds of Efremovka CV3 the noble gases with normal isotopic compositions are fractionated in different degree while the correlation of isotopic anomalous components is nearly constant. Some data for noble gases in DE-4 sample of Efremovka chondrite are considered. In contrast to DE-2 sample the DE-4 was treated except conc. HClO4, 220 C in addition with mixture of conc. H2SO4+H3PO4 (1:1), 220 C, twice. Noble gases analysis were performed in Germany at Max Plank Institute fur Chemie. Noble gases were released by oxidation of samples at stepped heating from 420 C to 810 C and by pyrolysis at 580, 590, and 680 C.