Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sparse approximation”

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 91 records · Page 5

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.↗

Space Weathering of Olivine: Samples, Experiments and Modeling

Olivine is a major constituent of chondritic bodies and its response to space weathering processes likely dominates the optical properties of asteroid regoliths (e.g. S- and many C-type asteroids). Analyses of olivine in returned samples and laboratory experiments provide details and insights regarding the mechanisms and rates of space weathering. Analyses of olivine grains from lunar soils and asteroid Itokawa reveal that they display solar wind damaged rims that are typically not amorphized despite long surface exposure ages, which are inferred from solar flare track densities (up to 10 (sup 7 y)). The olivine damaged rim width rapidly approaches approximately 120 nm in approximately 10 (sup 6 y) and then reaches steady-state with longer exposure times. The damaged rims are nanocrystalline with high dislocation densities, but crystalline order exists up to the outermost exposed surface. Sparse nanophase Fe metal inclusions occur in the damaged rims and are believed to be produced during irradiation through preferential sputtering of oxygen from the rims. The observed space weathering effects in lunar and Itokawa olivine grains are difficult to reconcile with laboratory irradiation studies and our numerical models that indicate that olivine surfaces should readily blister and amorphize on relatively short time scales (less than 10 (sup 3 y)). These results suggest that it is not just the ion fluence alone, but other variable, the ion flux that controls the type and extent of irradiation damage that develops in olivine. This flux dependence argues for caution in extrapolating between high flux laboratory experiments and the natural case. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind processing of olivine.

Keller, L. P.↗

The stellar population in the core of M15

The inner core of the globular cluster M15 within approximately 2 sec of the geometrical center has been explored with high-resolution images taken through several broad-band UV filter with the Faint Object Camera (FOC) on board the Hubble Space Telescope (HST). Approximately 210 stars in this region down to a 5 sigma detection limit of m(sub 220) = 21.5 were reliably identified and located on a UV - U color magnitude diagram for the first time. A majority of stars (about 70% of the total) observed this way lie above the expected main-sequence turn-off of this cluster and below the sparsely populated horizontal branch. The extension of the main sequence above the turn-off separates this population in two roughly equal components situated to the right and left of this line. Most of the former must be classical blue stragglers while the rest belong to a new, as yet unidentified, population of very blue stars. Possibilities include, but are not restricted to, well-mixed single stars, subdwarfs, and helium white dwarfs. Similar objects are also found just outside the core out to approximately 6 sec from the center, but the brighter, presumably more massive ones, are sharply confined to the core itself. The measured excess of bright blue stars and the relative deficiency of bright red giants in the core are consistent with the blue inward color gradient measured from the ground and imply that dynamical evolution can significantly affect the stellar population in the very dense central regions of a high-concentration globular cluster like M15.

De Marchi, Guido↗

Assessing Disaggregated SMAP Soil Moisture Products in the United States

A soil moisture (SM) disaggregation algorithm based on thermal inertia (TI) theory was implemented to downscale the Soil Moisture Active Passive (SMAP) Enhanced product(SPL2SMPE) from 9 km to 1 km over the continental United States. The algorithm applies land surface temperature and normalized difference vegetation index from Moderate Resolution Imaging Spectroradiometer (MODIS) at higher spatial resolution to estimate relative soil wetness within a coarse SMAP grid -this MODIS-derived relative wetness is then used to produce the downscaled SMAP SM. Results from the algorithm were evaluated in terms of their spatiotemporal coverage and accuracy using in situ measurements from SMAP Core Validation Sites (CVS), the US Department of Agriculture Soil Climate Analysis Network (USDA-SCAN), and the National Oceanic and Atmospheric Administration Climate Reference Network(NOAA-CRN). Results were also compared with the baselineSPL2SMPE and the SMAP/Sentinel-1 (SPL2SMAPS) 1kmproduct. Overall, the unbiased root mean square error (ubRMSE)of the disaggregated SM at the CVS using the TI approach is approximately 0.04 m3=m3, which is the SMAP mission requirement for the baseline products. The TI approach out performs the SMAP/Sentinel SL2SMAPS 1km product by approximately0.02 m3=m3. Over the agriculture/crop areas from SCAN and CRN sparse network stations, the TI approach exhibits better ubRMSE compared to SPL2SMPE andSPL2SMAPS byabout0.01 and 0.02 m3=m3, indicating its advantage in these areas. However, a drawback of this approach is that there are data gaps due to cloud cover as optical sensors cannot have a clear view of the land surface.

Soil moisture↗

Understanding Olivine and Pyroxene in Chondritic IDPs

There is a widespread depletion of Cr from the FeO-rich olivine in Wild-2, relative to type IIA chondrules(Frank et al., 2014)-the degree of this depletion is akin to that in unequilibrated chondrites such as Krymka (LL3.2) or Rainbow (CO3.2).Since Cr is highly mobile under even mild thermal metamorphism(as low as 200ºC), it is a sensitive indicator of such an event, showing greater depletion in smaller grains (i.e., matrix vs. chondrules) (Grossman and Brearley, 2005). Grossman and Brearley (2005) found that the distribution of Cr in FeO-rich olivine systematically changes as metamorphism increases between type 3.0 and type 3.2. Thus, Wild-2 olivine appears to carry evidence of mild thermal metamorphism. It is critical to determine whether there is any evidence for this from samples of other comets. For this we can examine an hydrous chondritic IDPs. Unfortunately, the data for olivine and pyroxene in IDPs is sparse, mainly decades old and lacks useful information on minor elements. We have begun to collect new compositional data for olivine and pyroxene in chondritic interplanetary dust particles (filling a major gap in our knowledge, as approximately half of these grains derive from comets), to determine whether the results obtained for Wild-2 are typical for comets in general, or whether Wild-2 (or its accretion components) has experienced an atypical geological history. A determination that a comet, or its accretion components, had experienced significant thermal metamorphism would greatly alter models of early solar system history.

Michael Zolensky↗

Low-order design and high-order simulation of active closed-loop control for aerospace structures under construction

Partially constructed/assembled structures in space are complicated enough but their dynamics will also be operating in closed-loop with feedback controllers. The dynamics of such structures are modeled by large-scale finite element models. The model dimension L is extremely large (approximately 10,000) while the numbers of actuators (M) and sensors (P) are small. The model parameters M(sub m) mass matrix, D(sub o) damping matrix, and K(sub o) stiffness matrix, are all symmetric and sparse (banded). Thus simulation of open-loop structure models of very large dimension can be accomplished by special integration techniques for sparse matrices. The problem of simulation of closed-loop control of such structures is complicated by the addition of controllers. Simulation of closed-loop controlled structures is an essential part of the controller design and evaluation process. Current research in the following areas is presented: high-order simulation of actively controlled aerospace structures; low-order controller design and SCI compensation for unmodeled dynamics; prediction of closed-loop stability using asymptotic eigenvalue series; and flexible robot manipulator control experiment.

Balas, Mark J.↗

Mitigating the Impacts of Measurement Error in the Quesst Mission Community Noise Study

Beginning in 2025, the NASA Quesst mission will conduct a series of community response tests involving flyovers of the X-59 aircraft at select localities across the United States. Several waves of a longitudinal survey will be administered over approximately one month of testing in order to capture perceptual responses to low-amplitude sonic booms, or “sonic thumps”. Simultaneously, noise exposure levels will be estimated by fusing model-based predictions with measurements taken from a sparse network of monitors in the region. As one of the aims of the study is to produce a dose-response curve, a regression model relating perceptual response to noise exposure levels, it is important to acknowledge the potential attenuation bias that results from measurement error in the estimated noise exposure levels. In this presentation we review and compare several methods for dealing with measurement error in generalized linear mixed models. The methods are demonstrated on simulated data and real data collected during past NASA risk reduction studies.

measurement error↗

Zeolite Formation and Weathering Processes in Dry Valleys of Antartica: Martian Analogs

Terrestrial weathering processes in cold-desert climates such as the Dry Valleys of Antarctica may provide an excellent analog to chemical weathering and diagenesis of soils on Mars. Detailed studies of soil development and the chemical and mineralogical alterations occurring within soil columns in Wright Valley, Antarctica show incredible complexity in the upper meter of soil. Previous workers noted the ice-free Dry Valleys are the best terrestrial approximations to contemporary Mars. Images returned from the Pathfinder and Spirit landers show similarities to surfaces observed within the Dry Valleys. Similarities to Mars that exist in these valleys are: mean temperatures always below freezing (-20 C), no rainfall, sparse snowfall-rapidly removed by sublimation, desiccating winds, diurnal freeze-thaw cycles (even during daylight hours), low humidity, oxidative environment, relatively high solar radiation and low magnetic fields . The Dry Valley soils contain irregular distributions and low abundances of soil microorganisms that are somewhat unusual on Earth. Physical processes-such as sand abrasion-are dominant mechanisms of rock weathering in Antarctica. However, chemical weathering is also an important process even in such extreme climates. For example, ionic migration occurs even in frozen soils along liquid films on individual soil particles. It has also been shown that water with liquid-like properties is present in soils at temperatures on the order of approx.-80 C and it has been observed that the percentage of oxidized iron increases with increasing soil age and enrichments in oxidized iron occurs toward the surface. The presence of evaporates is evident and appear similar to "evaporite sites" within the Pathfinder and Spirit sites. Evaporites indicate ionic migration and chemical activity even in the permanently frozen zone. The presence of evaporates indicates that chemical weathering of rocks and possibly soils has been active. Authogenic zeolites have been identified within the soil columns because they are fragile; i.e. they are euhedral, unabraded, and unfractured, strongly suggesting in situ formation. Their presence in Antarctic samples is another indication that diagenic processes are active in cold-desert environments. The presence of zeolites, and other clays along with halites, sulfates, carbonates, and hydrates are to be expected within the soil columns on Mars at the Gusev and Isidis Planitia regions. The presence of such water-bearing minerals beneath the surface supplies one of the requirements to support biological activity on Mars.

Gibson, E. K., Jr.↗

A Note on Substructuring Preconditioning for Nonconforming Finite Element Approximations of Second Order Elliptic Problems

In this paper an algebraic substructuring preconditioner is considered for nonconforming finite element approximations of second order elliptic problems in 3D domains with a piecewise constant diffusion coefficient. Using a substructuring idea and a block Gauss elimination, part of the unknowns is eliminated and the Schur complement obtained is preconditioned by a spectrally equivalent very sparse matrix. In the case of quasiuniform tetrahedral mesh an appropriate algebraic multigrid solver can be used to solve the problem with this matrix. Explicit estimates of condition numbers and implementation algorithms are established for the constructed preconditioner. It is shown that the condition number of the preconditioned matrix does not depend on either the mesh step size or the jump of the coefficient. Finally, numerical experiments are presented to illustrate the theory being developed.

Maliassov, Serguei↗

Towards a Multi-Variable Parametric Cost Model for Ground and Space Telescopes

Parametric cost models can be used by designers and project managers to perform relative cost comparisons between major architectural cost drivers and allow high-level design trades; enable cost-benefit analysis for technology development investment; and, provide a basis for estimating total project cost between related concepts. This paper hypothesizes a single model, based on published models and engineering intuition, for both ground and space telescopes: OTA Cost approximately (X) D(exp (1.75 +/- 0.05)) lambda(exp(-0.5 +/- 0.25) T(exp -0.25) e (exp (-0.04)Y). Specific findings include: space telescopes cost 50X to 100X more ground telescopes; diameter is the most important CER; cost is reduced by approximately 50% every 20 years (presumably because of technology advance and process improvements); and, for space telescopes, cost associated with wavelength performance is balanced by cost associated with operating temperature. Finally, duplication only reduces cost for the manufacture of identical systems (i.e. multiple aperture sparse arrays or interferometers). And, while duplication does reduce the cost of manufacturing the mirrors of segmented primary mirror, this cost savings does not appear to manifest itself in the final primary mirror assembly (presumably because the structure for a segmented mirror is more complicated than for a monolithic mirror).

Space Mirror Technology↗

Results of a statistical approach to rainfall estimation using Nimbus 5 6.7 micrometers and 11.5 micrometers THIR data

Nimbus 5 6.7 mm and 11.5 mm temperature humidity infrared radiometer (THIR) data were used in a simple multiple regression scheme to test the feasibility of using these data to estimate hourly rainfall. Throughout the test area (85 W to 105 W and 45 N to 30 N) subareas (8 deg x 6 deg) were chosen from which point to point and areal statistics were obtained. Four subsets of data were used. The first consisted of only those surface stations indicating precipitation whose latitude and longitude coincided with the THIR grid points. A second used surface stations 0.1 degree from the THIR grid points. The third was a combination of subsets one and two. A reciprocal distance weighting scheme was used to derive precipitation values in data sparse areas. A fourth subset was made using these data combined with the data from subsets one and two. Point estimates resulted in negative correlations between estimated and grid derived "surface" precipitation. One degree areal estimates showed a slight improvement with a correlation coefficient of approximately 0.11. Single regression areal estimates resulted in correlations of approximately 0.11 and 0.20 for the 6.7 mm and 11.5 mm data respectively. These poor results were attributed to problems which are inherent in the satellite data (location errors, short temporal span of data, wavelength of sensors, etc.) and the lack of sufficient surface data to better verify the satellite estimate.

Ormsby, J. P.↗

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.↗

Two-dimensional radiative transfer in cloudy atmospheres - The spherical harmonic spatial grid method

A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels of approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres.

Evans, K. F.↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

Identification of a precambrian rift through Missouri by digital image processing of geophysical and geological data

A newly discovered feature in the midcontinent - a gravity low that begins at a break in the midcontinent gravity high in SE Nebraska, extends across Missouri in a NW-SE direction, and intersects the Mississippi Valley graben to form the Pascola arch - is discussed. The anomaly varies from 120 to 160 km in width, extends approximately 700 km, and is best expressed in southern Missouri, where it has a Bouguer amplitude of about -34 mGal. It is noted that the magnitude of the anomaly cannot be explained on the basis of a thickened section of Paleozoic sedimentary rock. The gravity data and the sparse seismic refraction data for the region are found to be consistent with an increased crustal thickness beneath the gravity low. It is thought that the gravity anomaly is probably the present expression of a failed arm of a rifting event, perhaps one associated with the spreading that led to or preceded formation of the granite and rhyolite terrain of southern Missouri.

Guinness, E. A.↗

Treatment of cloudiness in regional scale - Numerical weather prediction models

A numerical weather prediction model is considered in which cloud water and precipitation are explicitly predicted as functions of three space dimensions and time. Approximations in the thermodynamic and continuity equations are discussed in detail, as well as the cumulus parameterization and the cloud microphysics parameterization. An example of a real-data forecast indicates that very accurate quantitative precipitation and cloud cover forecasts are possible using initial observations from only the relatively sparse rawinsonde network in the U.S.

Kreitzberg, C. W.↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗