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.
This report describes the development of an enstrophy model capable of predicting turbulence separation and its application to two airfoils at various angles of attack and Mach numbers. In addition, a two equation kappa-xi model with a tensor eddy viscosity was developed. Plans call for this model to be used in calculating three dimensional turbulent flows.
A numerical investigation was conducted to assess the accuracy of two turbulence models when computing non-axisymmetric nozzle-afterbody flows with propulsive jets. Navier-Stokes solutions were obtained for a Convergent-divergent non-axisymmetric nozzle-afterbody and its associated jet exhaust plume at free-stream Mach numbers of 0.600 and 0.938 at an angle of attack of 0 deg. The Reynolds number based on model length was approximately 20 x 10(exp 6). Turbulent dissipation was modeled by the algebraic Baldwin-Lomax turbulence model with the Degani-Schiff modification and by the standard Jones-Launder kappa-epsilon turbulence model. At flow conditions without strong shocks and with little or no separation, both turbulence models predicted the pressures on the surfaces of the nozzle very well. When strong shocks and massive separation existed, both turbulence models were unable to predict the flow accurately. Mixing of the jet exhaust plume and the external flow was underpredicted. The differences in drag coefficients for the two turbulence models illustrate that substantial development is still required for computing very complex flows before nozzle performance can be predicted accurately for all external flow conditions.
A thin-layer Navier-Stokes code, CFL3D, was utilized to compute the flow over a high-lift multi-element airfoil. This study was conducted to improve the prediction of high-lift flowfields using various turbulence models and improved glidding techniques. An overset Chimera grid system is used to model the three element airfoil geometry. The effects of wind tunnel wall modeling, changes to the grid density and distribution, and embedded grids are discussed. Computed pressure and lift coefficients using Spalart-Allmaras, Baldwin-Barth, and Menter's kappa-omega - Shear Stress Transport (SST) turbulence models are compared with experimental data. The ability of CFL3D to predict the effects on lift coefficient due to changes in Reynolds number changes is also discussed.
In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.
We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.
A general volume rendering technique is described that efficiently produces images of excellent quality from data defined over irregular grids having a wide variety of formats. Rendering is done in software, eliminating the need for special graphics hardware, as well as any artifacts associated with graphics hardware. Images of volumes with about one million cells can be produced in one to several minutes on a workstation with a 150 MHz processor. A significant advantage of this method for applications such as computational fluid dynamics is that it can process multiple intersecting grids. Such grids present problems for most current volume rendering techniques. Also, the wide range of cell sizes (by a factor of 10,000 or more), which is typical of such applications, does not present difficulties, as it does for many techniques. A spatial hierarchical organization makes it possible to access data from a restricted region efficiently. The tree has greater depth in regions of greater detail, determined by the number of cells in the region. It also makes it possible to render useful 'preview' images very quickly (about one second for one-million-cell grids) by displaying each region associated with a tree node as one cell. Previews show enough detail to navigate effectively in very large data sets. The algorithmic techniques include use of a kappa-d tree, with prefix-order partitioning of triangles, to reduce the number of primitives that must be processed for one rendering, coarse-grain parallelism for a shared-memory MIMD architecture, a new perspective transformation that achieves greater numerical accuracy, and a scanline algorithm with depth sorting and a new clipping technique.
The photochemistry of 32 neutral and 21 ion species in Triton's atmosphere is considered. Parent species N2, CH4, and CO (with a mixing ratio of 3 x 10(exp -4) in our basic model) sublime from the ice with rates of 40, 208, and 0.3 g/sq cm/b.y., respectively. Chemistry below 50 km is driven mostly by photolysis of methane by the solar and interstellar medium Lyman-alpha photons, producing hydrocarbons C2H4, C2H6, and C2H2 which form haze particles with precipitation rates of 135, 28, and 1.3 g/sq cm/b.y., respectively. Some processes are discussed which increase the production of HCN (by an order of magnitude to a value of 29 g/sq cm/b.y.) and involve indirect photolysis of N2 by neutrals. Reanalysis of the measured methane profiles gives an eddy diffusion coefficient K = 4 x 10(exp 3)sq cm/s above the tropopause and a more accurate methane number density near the surface, (3.1 +/- 0.8)x IO(exp 11)/cu cm. Chemistry above 200 km is driven by the solar EUV radiation (lambda less than 1000 A) and by precipitation of magnetospheric electrons with a total energy input of 10(exp 8) W (based on thermal balance calculations). The most abundant photochemical species are N, H2, H, 0, and C. They escape with the total rates of 7.7 x 10(exp 24)/ s, 4.5 x 10(exp 25)/s, 2.4 x 10(exp 25)/s, 4.4 x 10(exp 22)/s, and 1.1 x 10(exp 24), respectively. Atomic species are transported to a region of 50-200 km and drive the chemistry there. Ionospheric chemistry explains the formation of an E region at 150-240 km with HCO(+) as a major ion, and of an F region above 240 km with a peak at 320 km and C(+) as a major ion. The ionosphere above 500 km consists of almost equal densities of C(+) and N(+) ions. The model profiles agree with the measured atomic nitrogen and electron density profiles. A number of other models with varying rate coefficients of some reactions, differing properties of the haze particles (chemically passive or active), etc., were developed. These models show that there are four basic unknown values which have strong impacts on the composition and structure of the atmosphere and ionosphere. These values and their plausible ranges are the CO mixing ratio f(sub co) = 10(exp -4) - 10(exp -3), the magnetospheric electron energy input (1 +/- 0.5) x 10(exp 8) W, the rate coefficient of charge-exchange reaction N2(+) + C(kappa) = 10(exp -11) - 10(exp -10)cu cm/s, and the ion escape velocity upsilon(sub i) approx. equals 150 cm/s.
A fibrous rare earth selective emitter is approximated as an infinitely long cylinder. The spectral emittance, epsilon(lambda), is obtained by solving the radiative transfer equations with appropriate boundary conditions and uniform temperature. For optical depths, Kappa(R) = alpha(lambda)R, where alpha(lambda) is the extinction coefficient and R is the cylinder radius, greater than 1 the spectral emittance is nearly at its maximum value. There is an optimum cylinder radius, R(opt), for maximum emitter efficiency, eta(E). Values for R(opt) are strongly dependent on the number of emission bands of the material. The optimum radius decreases slowly with increasing emitter temperature, while the maximum efficiency and useful radiated power increase rapidly with increasing temperature.
The phase lag by which the earth's body tide follows the tidal potential is estimated for the principal lunar semidiurnal tide M(sub 2). The estimate results from combining recent tidal solutions from satellite tracking data and from Topex/Poseidon satellite altimeter data. Each data type is sensitive to the body-tide lag: gravitationally for the tracking data, geometrically for the altimetry. Allowance is made for the lunar atmospheric tide. For the tidal potential Love number kappa(sub 2) we obtain a lag epsilon of 0.20 deg +/- 0.05 deg, implying an effective body-tide Q of 280 and body-tide energy dissipation of 110 +/- 25 gigawatts.
FClO has been proposed as an intermediate in reactions involving ClF, Cl2O, and ClF3O, and it has been suggested as a molecule of atmospheric interest. It has been prepared in situ by the hydrolysis of ClF3. The pure rotational spectrum of FClO has been studied by conventional millimeter wave techniques and by microwave Fourier transform spectroscopy. Selected transitions were searched for using predictions based on an analysis of the nu(sub 1) band. FClO is an asymmetric prolate top, kappa = -0.8950 for F(35)ClO, with a rather small dipole component of 0.093 (4) D along the a-axis and a larger one of 1.93 (5) D along the b-axis. Transitions with 0 <= J <= 54 and 0 <= K(sub a) <= 18 were observed. Cl hyperfine splitting was generally observable throughout the spectrum with F-19 spin-rotation splitting observable as well in the microwave region. Structural parameters, harmonic force constants, and nuclear magnetic shielding parameters were derived and will be compared with data of related molecules, such as ClF3, ClF, FClO2, and FClO3. High resolution infrared spectra were taken in the regions of the FCl stretching mode and bending mode around 600 and 310/cm, respectively. A preliminary analysis indicates that the FCl stretch, near 596.86/cm for F(35)ClO, is in resonance with the dark overtone of delta near 617/cm. A brief progress report will be given.
The rotational spectra of OIO in its ground vibrational and first excited bending states have been observed for the first time. OIO was formed initially from the products of a microwave discharge in O2 passing over molecular iodine and later with greater yield in a DC discharge through a mixture of O2 and I2 vapor. OIO is an asymmetric prolate rotor (kappa = -0.690) with a (sup 2)B(sub 1) electronic ground state. Over 550 ground state transitions and over 160 transitions of the excited bending state have been included in the fits. The resulting parameters are well determined and will be compared to those recently published for OBrO and OClO. These will be interpreted in terms of the molecular geometry, harmonic force field, and electronic structure.
The Raman spectra of singled-wall carbon nanotubes (SWNT) produced by laser and are process were studied between 5 and 500 kappa. The line width vs. temperature dependence of the low-frequency Raman bands between 150 and 200/ cm deviates from that expected for phonon decay through phonon-phonon scattering mechanism. The experimental results and their analysis provided convincing evidence that each of the low-frequency Raman lines is a superposition of several narrower Raman lines corresponding to tubes of nearly the same diameter. The application of Raman spectroscopy to probe the distribution of SWNT by both diameter and chirality is discussed.
This paper presents diagnostic data obtained from the plume of a graphite composite target during carbon nanotube production by the double-pulse laser oven method. The insitu emission spectrum (300 to 650 nm) is recorded at different locations upstream of the target and at different delay times from the lasers (IR and green). Spectral features are identified as emissions from C2 (Swan System: a (sup 3)pi(sub g) - delta (sup 3)pi(sub u) and C3 (Comet Head System: A (sup 1)pi(sub u) - chi (sup 1)sigma(sub u) (sup +). Experimental spectra are compared with computed spectra to estimate vibrational temperatures of excited state C2 in the range of 2500 to 4000 kappa The temporal evolution of the 510 nm band of C2 is monitored for two target positions in various locations which shows confinement of the plume in the inner tube and increase in plume velocity with temperature. The excitation spectra of C2 are obtained by using a dye laser to pump the (0,1) transition of the Swan System and collecting the Laser Induced Fluorescence signal from C2 These are used to obtain "ground-state" rotational and vibrational temperatures which are close to the oven temperature. Images of the plume are also collected and are compared with the spectral measurements.
A supersonic coaxial jet facility is designed and experimental data are acquired suitable for the validation of CFD codes employed in the analysis of high-speed air-breathing engines. The center jet is of a light gas, the coflow jet is of air, and the mixing layer between them is compressible. The jet flow field is characterized using schlieren imaging, surveys with pitot, total temperature and gas sampling probes, and RELIEF velocimetry. VULCAN, a structured grid CFD code, is used to solve for the nozzle and jet flow, and the results are compared to the experiment for several variations of the kappa - omega turbulence model
Fires associated with tropical deforestation, land conversion, and land use greatly contribute to emissions as well as the depletion of carbon and nutrient pools. The objective of this research was to compare change detection techniques for identifying deforestation and cattle pasture formation during a period of early colonization and agricultural expansion in the vicinity of Jamari, Rond6nia. Multi-date Landsat Thematic Mapper (TM) data between 1984 and 1992 was examined in a 94 370-ha area of active deforestation to map land cover change. The Tasseled Cap (TC) transformation was used to enhance the contrast between forest, cleared areas, and regrowth. TC images were stacked into a composite multi-date TC and used in a principal components (PC) transformation to identify change components. In addition, consecutive TC image pairs were differenced and stacked into a composite multi-date differenced image. A maximum likelihood classification of each image composite was compared for identification of land cover change. The multi-date TC composite classification had the best accuracy of 78.1% (kappa). By 1984, only 5% of the study area had been cleared, but by 1992, 11% of the area had been deforested, primarily for pasture and 7% lost due to hydroelectric dam flooding. Finally, discrimination of pasture versus cultivation was improved due to the ability to detect land under sustained clearing opened to land exhibiting regrowth with infrequent clearing.
Major results of the S3-3 Langmuir sweep study are published. Studies show statistics and average density and temperature variation on auroral field lines up to 8000 km altitude. Alfven wave papers were published. Our model of Alfven wave propagation on auroral field lines was successfully extended to handle varying density and magnetic field for the inertial mode. The study showed that Alfven wave can create time-dispersed electron signatures. A study was undertaken to extend Langmuir sweep I-V curves to handle the case of an kappa electron distribution as well as Maxwellian. The manuscript is in preparation. Participated in International Space Science Institute study of Alfvenic structures which resulted in a group review paper. The proposed work was to develop an extended model of Alfven wave propagation along auroral field lines to study electron acceleration. As part of this work, a major task was to characterize density and temperature along auroral field lines by using spacecraft Langmuir sweep data. The work that was completed under this funding was successful at both tasks. Three papers have been published as part of this work and a fourth manuscript is in preparation.
Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.
The measurement methodologies of astrophysics experiments reflect the enormous variation of the astrophysical radiation itself. The diverse nature of the astrophysical radiation, e.g. cosmic rays, electromagnetic radiation, and neutrinos, is further complicated by the enormous span in energy, from the 1.95 Kappa relic neutrino background to cosmic rays with energy greater than 10(exp 20)eV. The measurement of gravity waves and search for dark matter constituents are also of astrophysical interest. Thus, the experimental techniques employed to determine the energy of the incident particles are strongly dependent upon the specific particles and energy range to be measured. This paper summarizes some of the calorimetric methodologies and measurements planned by future astrophysics experiments. A focus will be placed on the measurement of higher energy astrophysical radiation. Specifically, future cosmic ray, gamma ray, and neutrino experiments will be discussed.