Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Requirements for facilities and measurement techniques to support CFD development for hypersonic aircraft

The design of a hypersonic aircraft poses unique challenges to the engineering community. Problems with duplicating flight conditions in ground based facilities have made performance predictions risky. Computational fluid dynamics (CFD) has been proposed as an additional means of providing design data. At the present time, CFD codes are being validated based on sparse experimental data and then used to predict performance at flight conditions with generally unknown levels of uncertainty. This paper will discuss the facility and measurement techniques that are required to support CFD development for the design of hypersonic aircraft. Illustrations are given of recent success in combining experimental and direct numerical simulation in CFD model development and validation for hypersonic perfect gas flows.

Sellers, William L., III↗
NASA NTRS1992

LAI inversion from optical reflectance using a neural network trained with a multiple scattering model

The inversion of the leaf area index (LAI) canopy parameter from optical spectral reflectance measurements is obtained using a backpropagation artificial neural network trained using input-output pairs generated by a multiple scattering reflectance model. The problem of LAI estimation over sparse canopies (LAI < 1.0) with varying soil reflectance backgrounds is particularly difficult. Standard multiple regression methods applied to canopies within a single homogeneous soil type yield good results but perform unacceptably when applied across soil boundaries, resulting in absolute percentage errors of >1000 percent for low LAI. Minimization methods applied to merit functions constructed from differences between measured reflectances and predicted reflectances using multiple-scattering models are unacceptably sensitive to a good initial guess for the desired parameter. In contrast, the neural network reported generally yields absolute percentage errors of <30 percent when weighting coefficients trained on one soil type were applied to predicted canopy reflectance at a different soil background.

Smith, James A.↗
NASA NTRS1992

Estimation of effective aerodynamic roughness with altimeter measurements

A new method is presented for estimating the aerodynamic roughness length of heterogeneous land surfaces and complex landscapes using elevation measurements performed with an airborne laser altimeter and the Seasat radar altimeter. Land surface structure is characterized at increasing length scales by considering three basic landscape elements: (1) partial to complete canopies of herbaceous vegetation; (2) sparse obstacles (e.g., shrubs and trees); and (3) local relief. Measured parameters of land surface geometry are combined to obtain an effective aerodynamic roughness length which parameterizes the total atmosphere-land surface stress.

Menenti, M.↗
NASA NTRS1993

A three-dimensional pressure flux-split RNS application to sub/supersonic flow in inlets and ducts

The reduced Navier-Stokes (RNS) formulation is combined with a pressure based flux split procedure for the computational analysis of three dimensional flow. Only lowest order diffusion, that is required to satisfy the no-slip boundary condition is retained in this approximation. The governing equations are appropriately differenced such that the physical boundary conditions are closely coupled and combined to form a closed discrete system. No numerical or characteristic type boundary conditions are required. A global pressure relaxation procedure, i.e. multi-sweep PNS, is considered. A sparse matrix direct solver is applied to the crossplane solution. At the outflow boundary only the pressure or pressure gradient is prescribed. The resulting formulation has been applied to the solution of three-dimensional subsonic/supersonic internal flows in inlets and ducts.

Khosla, P. K.↗
NASA NTRS1993

The North Atlantic circulation in the early 1980s - An estimate from inversion of a finite-difference model

A finite-difference model of the North Atlantic is constructed for the purpose of making an estimate of the circulation through an inverse calculation. The data base is eclectic, and includes hydrography, oxygen, nutrients, current meter and float records, atmospheric momentum, heat and water vapor transfers, as well as estimates of certain integral fluxes. Owing to the available hydrographic database, the model resolution is restricted to 1 deg at best, and is much coarser in many aspects. This limited resolution is a major obstacle to accurate estimates of climatological fluxes. In its final form, there are about 9000 constraints in 29,000 formal unknowns plus 9000 noise unknowns. The system is solved as a tapered least-squares system by a sparse conjugate gradient algorithm. With the exception of a few float velocities, all constraints are found to be consistent within error estimates. The model produces estimates of large-scale fluxes and flux divergences for all conventional properties including heat and nutrients as well as carbon dioxide and alkalinity. Meridional fluxes of carbon are found to be indistinguishable from zero, whereas the North Atlantic tends to export nutrients to the south, but carry heat to the north. Traditional oceanographic depictions of the circulation through combination of nonsynoptic data into steady models may have reached their useful limit in the present calculation, as the conflicts between the data and physical requirements become quantitatively apparent.

Martel, Francoise↗
NASA NTRS1985

Theoretical studies of the marine sulfur cycle

Several reduced sulfur compounds are produced by marine organisms and then enter the atmosphere, where they are oxidized and ultimately returned to the ocean or the land. The oceanic dimethyl sulfide (DMS) flux, in particular, represents a significant fraction of the annual global sulfur input to the atmosphere. In the atmosphere, this gas is converted to sulfur dioxide (SO2), methane sulfonic acid, and other organic acids which are relatively stable and about which little is known. SO2 is a short lived gas which, in turn, is converted to sulfuric acid and other sulfate compounds which contribute significantly to acid rain. Because of the complexity of the sulfur system, it is not well understood even in the unperturbed atmosphere. However, a number of new observations and experiments have led to a significant increase in the understanding of this system. A number of one dimensional model experiments were conducted on the gas phase part of the marine sulfur cycle. The results indicate the measured concentration of DMS and the amplitude of its diurnal cycle are in agreement with estimates of its global flux. It was also found that DMS can make a large contribution to the background SO2 concentration in the free troposphere. Estimates of CS2 concentrations in the atmosphere are inconsistent with estimated fluxes; however, measured reaction rates are consistent with the observed steep tropospheric gradient in CS2. Observations of CS2 are extremely sparse. Further study is planned.

Toon, Owen B.↗
NASA NTRS1993

Issues and approach to develop validated analysis tools for hypersonic flows: One perspective

Critical issues concerning the modeling of low density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools, and the activity in the NASA Ames Research Center's Aerothermodynamics Branch is described. Inherent in the process is a strong synergism between ground test and real gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flowfield simulation codes are discussed. These models were partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions is sparse and reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high enthalpy flow facilities, such as shock tubes and ballistic ranges.

Deiwert, George S.↗
NASA NTRS1991

Discrete Kalman filtering equations of second-order form for control-structure interaction simulations

A second-order form of discrete Kalman filtering equations is proposed as a candidate state estimator for efficient simulations of control-structure interactions in coupled physical coordinate configurations as opposed to decoupled modal coordinates. The resulting matrix equation of the present state estimator consists of the same symmetric, sparse N x N coupled matrices of the governing structural dynamics equations as opposed to unsymmetric 2N x 2N state space-based estimators. Thus, in addition to substantial computational efficiency improvement, the present estimator can be applied to control-structure design optimization for which the physical coordinates associated with the mass, damping and stiffness matrices of the structure are needed instead of modal coordinates.

Park, K. C.↗
NASA NTRS1993

Th-230 - U-238 series disequilibrium of the Olkaria basalts Gregory Rift Valley, Kenya

U-Th disequilibrium analyses of the Naivasha basalts show a very small (U-238/Th-230) ratios which are lower than any previously analyzed basalts. The broadly positive internal isochron trend from one sample indicates that the basalts may have source heterogeneities, this is supported by earlier work. The Naivasha complex comprises a bimodal suite of basalts and rhyolites. The basalts are divided into two stratigraphic groups each of a transitional nature. The early basalt series (EBS) which were erupted prior to the Group 1 comendites and, the late basalt series (LBS) which erupted temporally between the Broad Acres and the Ololbutot centers. The basalts represent a very small percentage of the overall eruptive volume of material at Naivasha (less than 2 percent). The analyzed samples come from four stratigraphic units in close proximity around Ndabibi, Hell's Gate and Akira areas. The earliest units occur as vesicular flows from the Ndabibi plain. These basalts are olivine-plagioclase phyric with the associated hawaiites being sparsely plagioclase phyric. An absolute age of 0.5Ma was estimated for these basalts. The next youngest basalts flows occur as younger tuft cones in the Ndabibi area and are mainly olivine-plagioclase-clinopyroxcene phyric with one purely plagioclase phyric sample. The final phase of activity at Ndabibi resulted in much younger tuft cones consisting of air fall ashes and lapilli tufts. Many of these contain resorbed plagioclase phenocrysts with sample number 120c also being clinopyroxene phyric. The isotopic evidence for the basalt formation is summarized.

Black, S.↗
NASA NTRS1993

A structural origin for the cantaloupe terrain of Triton

Cantaloupe terrain is unique to Triton. It is Triton's oldest terrain and includes about 250,000 km sq. region displaying sparsely cratered, closely spaced, nearly circular dimples about 30-40 km across. This terrain is found on no other planet because, only on Triton the final major global thermal pulse (1) caused completed (or nearly) interior melting resulting in a cooling history where large thermal stresses shattered and contorted a thin, weak lithosphere, and (2) occurred after heavy bombardment so that the surface features were preserved. The cantaloupe terrain is composed of intersecting sets of structures (folds and/or faults) that have developed as a result of global compression generated by volumetric changes associated with cooling of Triton's interior. Further, it is proposed that these structures developed after the period of heavy bombardment, and resulted from the last major global thermal epoch in Triton's unique history (either caused by tidal or radio metric heating). Initially, as the body cooled and the structures formed, their surface topography was most likely modified by thermal relaxation of the warm surface ices. In other bodies like Mercury, thermal stresses generated from global cooling and contraction have resulted in widely spaced thrust faults, whereas on Triton, thermal stresses produced more closely-spaced folds and faults sets. This difference in structural style is probably due to differences in lithospheric properties (thickness, strength, etc.), the magnitude of stress (directly dependent on the thermal history), and when the structures formed, relative to the period of heavy bombardment.

Boyce, Joseph M.↗
NASA NTRS1993

Lunar Scout Two spacecraft gravity experiment

Measurement of the gravity field of the Moon has a high science priority because of its implications for the internal structure and thermal history of the Moon, and it has a high priority for future exploration activities because of the influence of lunar gravity on spacecraft navigation and orbit maintenance. The current state of knowledge in the lunar gravity field (and the uncertainty in the knowledge) is based primarily on data accumulated from the Lunar Orbiter and Apollo programs. Data are sparse and emphasize the equatorial band (+/- 30 deg) on the near side of the Moon. There are no tracking data on the far side and only the Lunar Orbiter 5 provides a small amount of high inclination data. A host of gravity models developed from different combinations of tracking data have large discrepancies in their predictions of spacecraft motion and orbit lifetimes. There are also large disagreements in the Mercator projections of the gravity acceleration from each model, especially on the far side, where the contours tend to have no obvious relationship with the local topography. The science and engineering requirements for global gravity field mapping will be satisfied with continuous radio metric tracking of Lunar Scout 1 in a low polar orbit using the Deep Space Network and Lunar Scout 2 in a high elliptical orbit. The gravity field of the Moon will be mapped during the Scout Program using a two spacecraft concept. In the two spacecraft concept, one spacecraft is placed in a high altitude eccentric orbit while the second spacecraft is in a low altitude polar orbit. The gravity experiment requires a radio frequency that will permit two-way Doppler tracking between the spacecraft and the Deep Space Network (DSN). Both spacecraft carry NASA standard transponder systems for data transmission to Earth as well as for tracking and orbit determination. Data sufficient to produce a gravity field map could be acquired within one month with this system.

Cheng, Andrew F.↗
NASA NTRS1993

A path-oriented matrix-based knowledge representation system

Experience has shown that designing a good representation is often the key to turning hard problems into simple ones. Most AI (Artificial Intelligence) search/representation techniques are oriented toward an infinite domain of objects and arbitrary relations among them. In reality much of what needs to be represented in AI can be expressed using a finite domain and unary or binary predicates. Well-known vector- and matrix-based representations can efficiently represent finite domains and unary/binary predicates, and allow effective extraction of path information by generalized transitive closure/path matrix computations. In order to avoid space limitations a set of abstract sparse matrix data types was developed along with a set of operations on them. This representation forms the basis of an intelligent information system for representing and manipulating relational data.

Feyock, Stefan↗
NASA NTRS1992

Preliminary results of polarization signatures for glacial moraines in the Mono Basin, Eastern Sierra Nevada

The valleys of the Mono Basin contain several sets of lateral and terminal moraines representing multiple stages of glaciation. The semi-arid climate with slow weathering rates preserved sequences of nested younger moraines within older ones. There is a well established relative chronology and recently exposure dating provided a new set of numerical dates. The moraines span the late Wisconsin (11-25 ka) to the Illinoian (130-190 ka) glaciations. The Mono Basin area was used as a 'calibration site' to establish remote dating techniques for eventual transfer to the more inaccessible but geomorphically and climatically similar moraines of the South American Andes Mountains. Planned polarimetric synthetic aperture radar (SAR) imagery acquired by JPL AIRSAR (South American Campaign) and SIR-C (Andes super-site) are analyzed to establish chronologies of previously undated moraine sequences in a study of Pleistocene climatic change in the Southern Hemisphere. The dry climate and sparse vegetation is also favorable for correlation of ground surface roughness with radar polarization signature. The slow weathering processes acting over thousands of years reduce the size, frequency, and angularity of surface boulders while increasing soil development on the moraines. Field observations based on this hypothesis result in relative ages consistent with those inferred from nested position within the valley. Younger moraines, therefore, will appear rougher than the older smoother moraines at scales measurable at AIRSAR wavelengths. Previously documented effects of ground surface roughness on polarization signatures suggest that analysis of moraine polarization signatures can be useful for relative dating. The technique may be extended to predict numerical ages. The data set reported were acquired on 8 Sep. 1989 with the JPL Airborne SAR (AIRSAR) collecting polarimetric imagery at C- (5.6 cm), L- (24 cm), and P-band (68 cm) with a flight-line parallel to the strike of the mountains. Phase calibration was performed on the analyzed scene by setting the co-phase of a smooth lake to zero as described. Absolute amplitude calibration was not possible because corner reflectors were not deployed.

Forster, Richard R.↗
NASA NTRS1993

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗
NASA NTRS1993

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗
NASA NTRS1993

Relationship of attenuation in a vegetation canopy to physical parameters of the canopy

A discrete scatter model is employed to compute the radiometric response (i.e. emissivity) of a layer of vegetation over a homogeneous ground. This was done to gain insight into empirical formulas for the emissivity which have recently appeared in the literature and which indicate that the attenuation through the canopy is proportional to the water content of the vegetation and inversely proportional to wavelength raised to a power around unity. The analytical result assumes that the vegetation can be modeled by a sparse layer of discrete, randomly oriented particles (leaves, stalks, etc.). The attenuation is given by the effective wave number of the layer obtained from the solution for the mean wave using the effective field approximation. By using the Ulaby-El Rayes formula to relate the dielectric constant of the vegetation to its water content, it can be shown that the attenuation is proportional to water content. The analytical form offers insight into the dependence of the empirical parameters on other variables of the canopy, including plant geometry (i.e. shape and orientation of the leaves and stalks of which the vegetation is comprised), frequency of the measurement and even the physical temperature of the vegetation. Solutions are presented for some special cases including layers consisting of cylinders (stalks) and disks (leaves).

Karam, M. A.↗
NASA NTRS1993

U-Pb isotopic ages and characteristics of ancient (greater than 4.0 Ga) lunar highland rocks

A review of the present (and certainly sparse) U-Pb isotopic database for ancient (greater than 4.0 Ga) lunar highland rocks indicates that both early-forming anorthosites and high-Mg suite rocks were formed from sources with a wide variety of U-238/Pb-204(mu) values as early as 4.44 Ga. The U-Pb data from Apollo 16 anorthosites 67075, 60025, and 62337 indicate source-mu values between approximately 55 and approximately 350; whereas values for high-Mg suite rocks are tentatively considered to be over 500, similar to KREEP sources. At present, I interpret the data to indicate that either early cumulates were formed slowly in a magma environment that allowed isotopic reequilibration or mixing over a long period of crystallization, or these cumulates formed from rapidly-evolving (both geochemically and isotopically) residual liquids. It is possible that an early bombardment period played a role in the mixing process.

Premo, W. R.↗
NASA NTRS1993

Computation of signal delays in RC networks

A model for signal delay computation in RC networks is presented. The strength of the paradigm is its generality and simplicity. The definition of delay is applicable to RC meshes with potential resistive attenuating paths to ground. The algorithms can also be applied to undriven circuits (static charge sharing) and circuits with initial charge. To compute the delays, each node in the network is explored locally to derive a system of sparse linear equations. The solutions of the system are delay values based on the Elmore time constant at each point in the circuit.

Hidalgo, Juan Carlos↗