Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 469 records · Page 26

LSPRAY: Lagrangian Spray Solver for Applications With Parallel Computing and Unstructured Gas-Phase Flow Solvers

Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, Manthena S.↗

Serial segment method for measuring remnant

For tracking tasks where a sum of sine waves forcing function is used it is often difficult and/or expensive to obtain the pilot's remnant in the vicinity of the sine waves. For the case where each sine wave has at least four times an integer number of cycle per run length, this paper illustrates the serial segments method for measuring remnant power spectral density in a frequency band centered on each sine wave. This method can be implemented on digital, hybrid, or analog Fourier coefficient analyzers, and is particularly advantageous on the latter since properties of Fourier coefficients are exploited to yield both a remnant measure and an improved estimate of the correlated component.

Magdaleno, R. E.↗

Cost comparison of competing local distribution systems for communication satellite traffic

The boundaries of market areas which favor various means for distributing communications satellite traffic are considered. The distribution methods considered are: control Earth station with cable access, rooftop Earth stations, Earth station with radio access, and various combinations of these methods. The least cost system for a hypothetical region described by number of users and the average cable access mileage is discussed. The region is characterized by a function which expresses the distribution of users. The results indicate that the least cost distribution is central Earth station with cable access for medium to high density areas of a region combined with rooftop Earth stations or (for higher volumes) radio access for remote users.

Dopfel, F. E.↗

Far-infrared investigation of the Taurus star-forming region using the IRAS database

The Infrared Astronomical Satellite (IRAS) has given us the first completely unbiased sky-survey in the far-infrared with wavebands centered at 12, 25, 60 and 100 microns. The Taurus-Auriga complex was selected as the first molecular cloud to be investigated in this study. The Taurus clouds were defined as lying between 04h and 05h in R.A. and +16 to +31 degrees in Dec., then the IRAS point-source catalogue was searched for sources with good or moderate quality fluxes in all three of the shortest IRAS bands. The sources which were selected in this way were then classified into subgroups according to their IRAS colors. Taurus is generally believed to be an area of low-mass star formation, having no luminous O-B associations within or near to the cloud complex. Once field stars, galaxies and planetary nebulae had been removed from the sample only the molecular cloud cores, T Tauri stars and a few emission-line A and B stars remained. The great majority of these objects are pre-main sequence in nature and, as stated by Chester (1985), main sequence stars without excess far-infrared emission would only be seen in Taurus if their spectral types were earlier than about A5 and then not 25 microns. By choosing our sample in this way we are naturally selecting the hotter and thus more evolved sources. To counteract this, the molecular cloud core-criterion was applied to soruces with good or moderate quality flux at 25, 60 and 100 microns, increasing the core sample by about one third. The candidate protostar B335 is only detected by IRAS at 60 and 100 microns while Taurus is heavily contaminated by cirrus at 100 microns. This means that detection at 25 microns is also required with those at 60 and 100 microns to avoid confusing a ridge of cirrus with a genuine protostar. The far-infrared luminosity function of these sources is then calculated and converted to the visual band by a standard method to compare with the field star luminosity function of Miller and Scalo. The eventual aim of this work is to obtain far-infrared luminosity functions for a number of molecular clouds which are known to be forming low-mass stars and to investigate how the slope is affected by changes in the density and turbulence of material.

Hughes, J. D.↗

Time-periodic variations in stratospheric ozone from satellite observations

The time-periodic variations in stratospheric zonal mean ozone number density are examined with emphasis on annual and semi-annual oscillations. The data base for this analysis is the ozone observation from the SAGE II satellite instrument. A multiple linear regression method is adopted for the analysis. The results show that the amplitudes and phases of the time-periodic ozone variations are functions of altitude and latitude.

Wang, Pi-Huan↗

Electrodeposited CuInSe2 Thin Film Junctions

We have investigated thin films and junctions based on copper indium diselenide (CIS) which have been grown by electrochemical deposition. CIS is a leading candidate for use in polycrystalline thin film photovoltaic solar cells. Electrodeposition is a cost-effective method for producing thin-film CIS. We have produced both p and n type CIS thin films from the same aqueous solution by simply varying the deposition potential. A CIS pn junction was deposited using a step-function potential. Stoichiometry of the single layer films was determined by energy dispersive spectroscopy. Carrier densities of these films increased with deviation from stoichiometry, as determined by the capacitance versus voltage dependence of Schottky contacts. Optical bandgaps for the single layer films as determined by transmission spectroscopy were also found to increase with deviation from stoichiometry. Rectifying current versus voltage characteristics were demonstrated for the Schottky barriers and for the pn junction.

Raffaelle, R. P.↗

Electrodeposited CulnSe2 Thin Film Junctions

We have investigated thin films and junctions based on copper indium diselenide (CIS) which have been grown by electrochemical deposition. CIS is a leading candidate for use in polycrystalline thin film photovoltaic solar cells. Electrodeposition is a cost-effective method for producing thin-film CIS. We have produced both p and n type CIS thin films from the same aqueous solution by simply varying the deposition potential. A CIS pn junction was deposited using a step-function potential. Stoichiometry of the single layer films was determined by energy dispersive spectroscopy. Carrier densities of these films increased with deviation from stoichiometry, as determined by the capacitance versus voltage dependence of Schottky contacts. Optical bandgaps for the single layer films as determined by transmission spectroscopy were also found to increase with deviation from stoichiometry. Rectifying current versus voltage characteristics were demonstrated for the Schottky barriers and for the pn junction.

Raffaelle, R. P.↗

A Hybrid Numerical Method for Turbulent Mixing Layers

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental data.

Georgiadis, Nicholas J.↗

Modeling Growth of Nanostructures in Plasmas

As semiconductor circuits shrink to CDs below 0.1 nm, it is becoming increasingly critical to replace and/or enhance existing technology with nanoscale structures, such as nanowires for interconnects. Nanowires grown in plasmas are strongly dependent on processing conditions, such as gas composition and substrate temperature. Growth occurs at specific sites, or step-edges, with the bulk growth rate of the nanowires determined from the equation of motion of the nucleating crystalline steps. Traditional front-tracking algorithms, such as string-based or level set methods, suffer either from numerical complications in higher spatial dimensions, or from difficulties in incorporating surface-intense physical and chemical phenomena. Phase field models have the robustness of the level set method, combined with the ability to implement surface-specific chemistry that is required to model crystal growth, although they do not necessarily directly solve for the advancing front location. We have adopted a phase field approach and will present results of the adatom density and step-growth location in time as a function of processing conditions, such as temperature and plasma gas composition.

Hwang, Helen H.↗

The Plasma Environment at Mercury

Mercury is the least explored terrestrial planet, and the one subjected to the highest flux of solar radiation in the heliosphere. Its highly dynamic, miniature magnetosphere contains ions from the exosphere and solar wind, and at times may allow solar wind ions to directly impact the planet's surface. Together these features create a plasma environment that shares many features with, but is nonetheless very different from, that of Earth. The first in situ measurements of plasma ions in the Mercury space environment were made only recently, by the Fast Imaging Plasma Spectrometer (FIPS) during the MESSENGER spacecraft's three flybys of the planet in 2008-2009 as the probe was en route to insertion into orbit about Mercury earlier this year. Here. we present analysis of flyby and early orbital mission data with novel techniques that address the particular challenges inherent in these measurements. First. spacecraft structures and sensor orientation limit the FIPS field of view and allow only partial sampling of velocity distribution functions. We use a software model of FIPS sampling in velocity space to explore these effects and recover bulk parameters under certain assumptions. Second, the low densities found in the Mercury magnetosphere result in a relatively low signal-to-noise ratio for many ions. To address this issue, we apply a kernel density spread function to guide removal of background counts according to a background-signature probability map. We then assign individual counts to particular ion species with a time-of-flight forward model, taking into account energy losses in the carbon foil and other physical behavior of ions within the instrument. Using these methods, we have derived bulk plasma properties and heavy ion composition and evaluated them in the context of the Mercury magnetosphere.

Raines, James M.↗

Surface photometry of spiral galaxies. I - Theoretical color variation and surface brightness across spiral arms

A method for the construction of color variations and surface brightness across a galactic spiral arm, using the concepts and results of the density-wave theory, is developed. The color variation measured by the index Q varies with the ratio of the light of young stars to that of the old stars, the initial mass function, and the pattern speed. The fact that the variation of brightness and color across a spiral arm is relatively smooth and symmetrical and not sharp as in the galactic shock picture is due to the smearing-out effects of star formation inside a cloud in the regions behind the galactic shock, and the tendency of new stars to fall back to small galactocentric distances after their formation because of their lower angular velocity at birth. Both effects arise naturally in the context of the density-wave theory.

Yuan, C.↗

Theory of inhomogeneous fluids

A theory of nonuniform liquids is presented which is based on the Yvon-Born-Green equation for the one-particle density and the Ornstein-Zernike equation of an inhomogeneous system. The necessary closure is effected by exploiting the solution of a modified hypernetted-chain equation and making a local approximation on the (highly universal) bridge function. The theory is successfully applied to model problems that have also been studied by direct-simulation methods. A brief generalization to quantum liquids is also given.

Nieminen, R. M.↗

A mathematical characterization of vegetation effect on microwave remote sensing from the Earth

In passive microwave remote sensing of the earth, a theoretical model that utilizes the radiative transfer equations was developed to account for the volume scattering effects of the vegetation canopy. Vegetation canopies such as alfalfa, sorghum, and corn are simulated by a layer of ellipsoidal scatterers and cylindrical structures. The ellipsoidal scatterers represent the leaves of vegetation and are randomly positioned and oriented. The orientation of ellipsoids is characterized by a probability density function of Eulerian angles of rotation. The cylindrical structures represent the stalks of vegetation and their radii are assumed to be much smaller than their lengths. The underlying soil is represented by a half-space medium with a homogeneous permittivity and uniform temperature profile. The radiative transfer quations are solved by a numerical method using a Gaussian quadrature formula to compute both the vertical and horizontal polarized brightness temperature as a function of observation angle. The theory was applied to the interpretation of experimental data obtained from sorghum covered fields near College Station, Texas.

Choe, Y.↗

An APL program for the distribution of energy deposition by charged particles passing through thin absorbers

An APL program which numerically evaluates the probability density function (PDF) for the energy deposited in a thin absorber by a charged particle is proposed, with application to the construction, pointing, and control of spacecraft. With this program, the PDF of the restricted energy loss distribution of Watts (1973) is derived, and Vavilov's (1957) distribution is obtained by proper parameter selection. The method is demonstrated with the example of the effect of charged particle induced radiation on the Hubble Space Telescope (HST) pointing accuracy. A Monte Carlo study simulates the photon noise caused by charged particles passing through the photomultiplier tube window, and the stochastic variation of energy loss is introduced into the simulation by generating random energy losses from a power law distribution. The program eliminates annoying loop procedures, and model parameter sensitivity can be studied using the graphical output.

Howell, L. W.↗

Equivalence of physically based statistical fracture theories for reliability analysis of ceramics in multiaxial loading

The present comparison of the Batdorf (1974) flaw density and orientation distribution approach with Evans' (1978) elemental strength approach, with a view to identities in fracture criteria and distribution functions, notes that despite their fundamental differences in multiaxial loading fracture probabilities, the two approaches yield identical predictions. Lamon's (1988) assertion to the contrary, in light of different theoretical predictions by the two methods for the case of alumina disks loaded in flexure, is demonstrated to be in error.

Chao, Luen-Yuan↗

New method for estimating low-earth-orbit collision probabilities

An unconventional but general method is described for estimating the probability of collision between an earth-orbiting spacecraft and orbital debris. This method uses a Monte Caralo simulation of the orbital motion of the target spacecraft and each discrete debris object to generate an empirical set of distances, each distance representing the separation between the spacecraft and the nearest debris object at random times. Using concepts from the asymptotic theory of extreme order statistics, an analytical density function is fitted to this set of minimum distances. From this function, it is possible to generate realistic collision estimates for the spacecraft.

Vedder, John D.↗

Measuring ionospheric electron density using the plasma frequency probe

During the past decade, the plasma frequency probe (PFP) has evolved into an accurate, proven method of measuring electron density in the ionosphere above about 90 km. The instrument uses an electrically short antenna mounted on a sounding rocket that is immersed in the plasma and notes the frequency where the antenna impedance is large and nonreactive. This frequency is closely related to the plasma frequency, which is a direct function of free electron concentration. The probe uses phase-locked loop technology to follow a changing electron density. Several sections of the plasma frequency probe circuitry are unique, especially the voltage-controlled oscillator that uses both an electronically tuned capacitor and inductor to give the wide tuning range needed for electron density measurements. The results from two recent sounding rocket flights (Thunderstorm II and CRIT II) under vastly different plasma conditions demonstrate the capabilities of the PFP and show the importance of in situ electron density measurements of understanding plasma processes.

Jensen, Mark D.↗

Bayesian Inference in Satellite Gravity Inversion

To solve a geophysical inverse problem means applying measurements to determine the parameters of the selected model. The inverse problem is formulated as the Bayesian inference. The Gaussian probability density functions are applied in the Bayes's equation. The CHAMP satellite gravity data are determined at the altitude of 400 kilometer altitude over the South part of the Pannonian basin. The model of interpretation is the right vertical cylinder. The parameters of the model are obtained from the minimum problem solved by the Simplex method.

Kis, K. I.↗