Search NASASearch

SEARCH · Search NASA

Results for “Gradient estimation”

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

Multigrid properties of upwind-biased data reconstructions

The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.

Warren, Gary P.

Advanced Observation Operators for GPS Radio Occultation: Validation of the 2D Ray Tracing Approach with CHAMP and SAC-C bending angle and refractivity data - Part 1

Global Positioning System (GPS) Radio Occultations (RO) bending angles and refractivity data characterize mostly the vertical structure of the Earth's atmosphere. We answer the question whether proper simulation of GPS RO data for data assimilation can be obtained with one-dimensional vertical operators, or if accounting also for horizontal atmospheric structures via ray-tracing makes a positive difference when compared with real data. We present a detailed implementation of a geometrical optics multi-plane two-dimensional (2D) ray-tracing as an observation operator to simulate GPS RO bending angles and refractivities within the Finite Volume Data Assimilation System (FVDAS). Comparisons of the outputs of that 2D observation operator with those of simpler ID observation operators are used to generate estimates of errors induced by neglecting tangent point drift (TPD) and horizontal gradients (HG). These error estimates are then confronted with errors estimates derived using 6335 real CHAMP and SAC-C occultations. The agreement for TPD-induced (HG-induced) errors is remarkably positive at altitudes 10-30 km (below 10 km). Comparisons in bending angles O - B STD of the outputs of the multi-plane 2D ray-tracer with those of a vertical Abel transform show reductions of about 8% of the usual O - B bending angle STD due to TPD in the stratosphere (3% due to HG, in the troposphere only). In terms of refractivity, the O - B STD reductions are about 1520% for TPD and 3-5% for HG in the same regions. These reductions are obtained using either 6-hour forecasts or analyses as backgrounds, and using Geometrical Optics (GO) or Canonical Transform (CT) data.

Poli, P.

Kurtosis Approach Nonlinear Blind Source Separation

In this paper, we introduce a new algorithm for blind source signal separation for post-nonlinear mixtures. The mixtures are assumed to be linearly mixed from unknown sources first and then distorted by memoryless nonlinear functions. The nonlinear functions are assumed to be smooth and can be approximated by polynomials. Both the coefficients of the unknown mixing matrix and the coefficients of the approximated polynomials are estimated by the gradient descent method conditional on the higher order statistical requirements. The results of simulation experiments presented in this paper demonstrate the validity and usefulness of our approach for nonlinear blind source signal separation Keywords: Independent Component Analysis, Kurtosis, Higher order statistics.

Duong, Vu A.

Kurtosis Approach for Nonlinear Blind Source Separation

In this paper, we introduce a new algorithm for blind source signal separation for post-nonlinear mixtures. The mixtures are assumed to be linearly mixed from unknown sources first and then distorted by memoryless nonlinear functions. The nonlinear functions are assumed to be smooth and can be approximated by polynomials. Both the coefficients of the unknown mixing matrix and the coefficients of the approximated polynomials are estimated by the gradient descent method conditional on the higher order statistical requirements. The results of simulation experiments presented in this paper demonstrate the validity and usefulness of our approach for nonlinear blind source signal separation.

kurtosis

Influence of Exercise Modality on Cerebral-Ocular Hemodynamics and Pressures

Background: Moderate and high intensity aerobic or resistance exercise has clearly identified benefits for cardiac, muscle, and bone health. However, the impact of such exercise - either as a mitigating or an exacerbating factor - on the development of the visual impairment and intracranial pressure syndrome (VIIP) is unknown. Accordingly, our aim was to characterize the effect of an acute bout of resistance (RE), moderate-intensity continuous (CE), and high-intensity interval exercise (IE) during a cephalad fluid shift on cerebral-ocular hemodynamics and pressures. Methods: 10 male subjects (36 plus or minus 9 years) completed 4 testing days in a 15 degree head-down tilt (HDT): (1) assessment of maximum volume of O (sub 2), (2) RE session (4 sets of 12 repetition maximum leg press exercise), (3) CE session (30 minutes of cycling at 60 percent maximum volume of O (sub 2)), and (4) IE session (4 by 4-minute intervals of exercise at 85 percent maximum volume of O (sub 2) with 3-minute active rest periods). During each session, blood flow (Vivid-e, GE Healthcare) in extracranial arteries (common carotid artery, CCA; internal carotid artery, ICA; external carotid artery, ECA and vertebral artery, VA), and mean blood flow velocity in middle cerebral artery (MCA), internal jugular pressure (IJP; VeinPress), and intraocular pressure (IOP; Icare PRO) were measured at rest, at the end of each resistance or interval set, and every 5 minutes during continuous exercise. Translaminar pressure gradient (TLPG) was estimated by subtracting IJP from IOP. Results: There were no differences across days in pre-exercise resting blood flows or pressures. IOP decreased slightly from HDT rest (20.2 plus or minus 2.3 millimeters of mercury) to exercise (RE: 19.2 plus or minus 2.8 millimeters of mercury; CE: 18.9 plus or minus 3.2 millimeters of mercury; IE: 20.1 plus or minus 2.8 millimeters of mercury), while IJP decreased during CE (31.6 plus or minus 9.5 millimeters of mercury) and RE (32.0 plus or minus 8.1 millimeters of mercury), and increased during IE (35.1 plus or minus 9.5 millimeters of mercury) from HDT rest (33.3 plus or minus 6.5 millimeters of mercury). Estimated TLPG was increased during IE only. Compared to RE and CE, IE resulted in the greatest increase in MCA blood flow velocity and extracranial artery blood flow. Conclusions: These preliminary results suggest that high-intensity IE acutely increases cerebral blood flow, IJP, and TLPG. Alterations in TLPG is one mechanism that may contribute to optic nerve sheath edema in astronauts. Accordingly, acutely raising IOP and/or orbital pressure during exercise could optimize cerebral-ocular pressures during spaceflight.

Scott, J.

Scientific and Technical Assistance for the Deployment of a Flexible Airborne Spectrometer System During C-MAPExp and COMEX: COMEX Version - 1.9T

The COMEX (CO2 and MEthane eXperiment) campaign supports the mission definition of CarbonSat and HyspIRI (Hyperspectral Infrared Imager) by providing representative airborne remote sensing data ‐ MAMAP (Methane Airborne MAPper) for CarbonSat; the Airborne Visual InfraRed Imaging Spectrometer (Classic & Next Generation) AVIRIS‐C/AVIRIS‐NG for HyspIRI ‐ as well as ground-based and airborne in‐situ data. The objectives of the COMEX campaign activities are (see Campaign Implementation Plan (RD‐4)): 1. Investigate spatial/spectral resolution trade‐offs for CH4 anomaly detection and flux inversion by comparison of MAMAP‐derived emission estimates with AVIRIS/AVIRIS‐NG derived data. 2. Evaluate sun‐glint observation geometry on CH4 retrievals for marine sources. 3. Characterize the effect of Surface Spectral Reflectance (SSR) heterogeneity on trace gas retrievals of CO2 and CH4 for medium and low‐resolution spectrometry. 4. Identify benefits from joint SWIR/TIR (ShortWave InfraRed/Thermal InfraRed ) data for trace gas detection and retrieval by comparison of MAMAP and AVIRIS/AVIRIS‐NG NIR/SWIR data with MAKO (Aerospace Corp.)TIR data. The ability to derive emission source strength for a range of strong emitting targets by remote sensing will be evaluated from combined AVIRIS‐NG and MAMAP data, adding significant value to the HyspIRI campaign AVIRIS‐NG dataset. The data will be used to quantify anomalies in atmospheric CO2 and CH4 from strong local greenhouse gas sources e.g. localized industrial complexes, landfills, etc. and to derive CO2 and CH4 emissions estimates from atmospheric gradient measurements. The original campaign concept was developed by University of Bremen and BRI. The COMEX campaign is funded bilaterally by NASA and ESA (European Space Agency). Whereas NASA funds the US part of the project via a contract with Dr. Ira Leifer, BRI (Bubbleology Research International), the contribution of MAMAP to the COMEX campaign is funded by ESA within the COMEX‐E project and NASA with respect to a 50 percent contribution to the flight-related costs of flying MAMAP on a US aircraft. The Data Acquisition Report (RD‐9) describes the instrumentation used, the measurements made by the team during the COMEX campaign in May/June 2014 and August/September 2014 in California, and an initial assessment of the data quality.

Point Sources

Motion-Contrast Sensitivity: Visibility of Motion Gradients of Various Spatial Frequencies

The purpose of our experiments was to estimate basic sensitivity to motion gradients and to evaluate the evidence for second-order integration and differentiation of motion signals. We measured sensitivity to spatially sinusoidal contrast modulation between two oppositely moving bandpass-filtered noise images. The motion-contrast sensitivity function, defined as the inverse of threshold modulation amplitude as a function of modulation spatial frequency, was bandpass in shape with declines at both highest and lowest frequencies. The functions for three noise spatial frequencies had approximately the same shape when modulation frequency was expressed as a fraction of noise frequency. We compared the data with a model in which linear motion filters, whose outputs are squared or rectified, are followed by a second stage of excitatory or inhibitory pooling. The data are consistent with a model in which (1) all excitatory pooling occurs at the linear stage and (2) the second stage contains a large inhibitory pooling area, with a radius approximately eight times that of the linear receptive field.

Watson, Andrew B.

Simulation studies of wide and medium field of view earth radiation data analysis

A parameter estimation technique is presented to estimate the radiative flux distribution over the earth from radiometer measurements at satellite altitude. The technique analyzes measurements from a wide field of view (WFOV), horizon to horizon, nadir pointing sensor with a mathematical technique to derive the radiative flux estimates at the top of the atmosphere for resolution elements smaller than the sensor field of view. A computer simulation of the data analysis technique is presented for both earth-emitted and reflected radiation. Zonal resolutions are considered as well as the global integration of plane flux. An estimate of the equator-to-pole gradient is obtained from the zonal estimates. Sensitivity studies of the derived flux distribution to directional model errors are also presented. In addition to the WFOV results, medium field of view results are presented.

Green, R. N.

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.

Using SMAP Level-4 Soil Moisture to Constrain MOD16 Evapotranspiration over the Contiguous USA

Evapotranspiration (ET) is a key hydrologic variable linking the Earth’s water, carbon and energy cycles. At large spatial scales, remote sensing-based (RS) models are often used to quantify ET. Despite the large number of RS ET models available, few include soil moisture as a key environmental input, which can degrade model accuracy and utility. Here, we use model assimilation enhanced soil moisture estimates from the NASA SMAP (Soil Moisture Active Passive) mission as a water supply control in the MOD16 ET algorithm framework. SMAP-derived daily surface (0-5cm depth) and root zone (0-1m depth) soil moisture are used with MODIS (Moderate Resolution Imaging Spectroradiometer) vegetation observations, and 4km gridded regional surface meteorology (Gridmet) as primary inputs for estimating daily ET and underlying model soil and stomatal conductance terms. We calibrated the model environmental response parameters using tower eddy covariance ET observations representing major North American biomes. The model ET results were validated using a holdout set of tower observations spanning a large regional climate gradient. The updated ET estimates outperform the baseline MOD16 product across all tower validation sites (RMSE = 0.758 vs 1.108 mm day-1; R2 = 0.68 vs0.45, respectively). Smaller relative improvements were obtained using a recalibrated model with 4km Gridmet meteorology, but no soil moisture control (RMSE = 0.813 mm day-1; R2 = 0.66), indicating that these changes are essential for the improved model performance. The soil moisture-constrained model improvements and relative benefits from the SMAP observations are greater in arid climates, consistent with stronger soil moisture control on ET in water-limited regions. The use of SMAP soil moisture as an additional model constraint improves MOD16 regional performance and provides a new framework for investigating both soil and atmosphere controls on ET.

SMAP

Differential measurement of cosmic-ray gradient with respect to interplanetary current sheet

Simultaneous magnetic field and charged particle measurements from the Voyager pacecraft at heliographic latitude separations from 10 deg to 21 deg are used to determine the latitude gradient of the galactic cosmic ray flux with respect to the interplanetary current sheet. By comparing the ratio of cosmic ray flux at Voyager 1 to that at Voyager 2 during periods when both spacecraft are first north and then south of the interplanetary current sheet, we find an estimate of the latitudinal gradient with respect to the current sheet of approximately -0.15 + or- 0.05%/deg under restricted interplanetary conditions.

Christon, S. P.

Differential measurement of cosmic-ray gradient with respect to interplanetary current sheet

Simultaneous magnetic field and charged particle measurements from the Voyager spacecraft at heliographic latitude separations from 10 deg. to 21 deg. are used to determine the latitude gradient of the galactic cosmic ray flux with respect to the interplanetary current sheet. By comparing the ratio of cosmic ray flux at Voyager 1 to that a Voyager 2 during periods when both spacecraft are first nort and then south of the interplanetary current sheet, we find an estimate of the latitudinal gradient with respect to the current sheet of approximately -0.15 + or 0.05% deg under restricted interplanetary conditions.

Christon, S. P.

Metamorphic P-T paths and Precambrian crustal growth in East Antarctica

The metamorphic constraints on crustal thicknesses in Archean and post-Archean terranes are summarized along with possible implications for tectonic processes. It is important to recognize that P-T estimates represent perturbed conditions and should not be used to estimate steady state geothermal gradients or crustal thicknesses. The example is cited of the Dora Maira complex in the French Alps, where crustal rocks record conditions of 35 kbar and 800 C, implying their subduction to depths of 100 km or more, followed by subsequent uplift to the surface. Therefore such P-T estimates tell more about processes than crustal thicknesses. More importantly, according to the author, are determinations of P-T paths, particularly coupled with age measurements, because these may provide constraints on how and when perturbed conditions relax back to steady state conditions. P-T paths are illustrated that should be expected from specific tectonic processes, including Tibetan style collision, with and without subsequent extension, rifting of thin or thickened crust, and magmatic accretion. Growth of new crust, associated with magmatic accretion, for example, could possibly be monitored with these P-T paths.

Harley, S. L.

Studies using IMP, Voyager and Pioneer cosmic ray data to determine the size of the heliosphere

The purpose of this project was to use the cosmic ray data from the IMP, Voyager and Pioneer spacecraft in the heliosphere out to approximately 65 AU to estimate the size of the heliosphere. We used several techniques to develop a consistent picture of the size of the heliosphere. The first method used a response function approach which determined the intensity as a function of time by scaling the modulation effect as they move outward and eventually reach the boundary of the heliosphere. In this model the effects of transient cosmic ray disturbances is included. A second approach using the perturbation method in which drifts are considered as a perturbation to the standard diffusion-convection modulation models was not fully developed. In a third approach the location of the modulation boundary beyond the termination shock was estimated using observations of the intensity and radial gradients between Voyager 2 and Pioneer 10 along with new estimates of the interstellar intensity of more than 70 MeV galactic cosmic rays. Using this method we found that for 7 years, from 1983 to 1990, the modulation boundary remained constant at 83 +/- 5 AU. We infer from these studies that a modulation boundary can be estimated only by extrapolating the observed radial gradients when the solar magnetic field polarity is such that cosmic-ray particles are drifting in the heliosphere inward toward the Earth along the neutral sheet. The boundary distance is larger than the estimates of the location of the termination shock at 67 +/- 5 AU using the same method. Two other studies partially supported by this grant are attached. The first deals with the recovery period of the greater than 70 MeV cosmic rays in the outer heliosphere from 1992-1995. In the second paper we compare the rigidity dependence of the 11-year cosmic ray variation at the Earth in two cycles of opposite solar magnetic field polarity.

Lockwood, John A.

Arctic Snow Depth and Sea Ice Thickness From ICESat‐2 and CryoSat‐2 Freeboards: A First Examination

We present a first examination of Arctic sea ice snow depth estimates from differencing satellite lidar (ICESat‐2) and radar (CryoSat‐2) freeboards. These estimates cover the period between 14 October 2018 and the end of April 2019. Snow depth is related to freeboard differences by the refractive index/bulk density of the snow layer—the only free parameter in the approach. Area‐averaged snow depth ranges from 9 cm (on first‐year ice: 5 cm, multiyear ice: 14 cm) in late October to 19 cm (first‐year ice: 17 cm, multiyear ice: 27 cm) in April; on average, this snow is thinner over FYI. Spatial patterns and gradients of snow depth estimates compare well with reconstructions using snowfall from ERA‐Interim and ERA5, although snowfall from ERA5 is systematically higher. For all months, the results suggest that ~50% of the total freeboard is comprised of snow. Retrievals are within a few centimeters of snow depth data acquired by Operation IceBridge in April 2019. Sources of uncertainties associated with this freeboard‐differencing approach are discussed. Further, sea ice thicknesses calculated using the retrieved snow depth and a modified climatology are contrasted. Comparatively, the snow depth and calculated ice thickness using a modified climatology are higher by ~5 cm and 0.33 m, although these differences are not uniform throughout the season. Snow accumulation was slower between October and December but increased between December and January, unlike the modified climatology, which exhibited a monotonic accumulation for all months. Future opportunities for assessment and improvement of these estimates are discussed.

Arctic Ocean

Using PEACE Data from the four CLUSTER Spacecraft to Measure Compressibility, Vorticity, and the Taylor Microscale in the Magnetosheath and Plasma Sheet

We present determinations of compressibility and vorticity in the magnetosheath and plasma sheet using moments from the four PEACE thermal electron instruments on CLUSTER. The methodology used assumes a linear variation of the moments throughout the volume defined by the four satellites, which allows spatially independent estimates of the divergence, curl, and gradient. Once the vorticity has been computed, it is possible to estimate directly the Taylor microscale. We have shown previously that the technique works well in the solar wind. Because the background flow speed in the magnetosheath and plasma sheet is usually less than the Alfven speed, the Taylor frozen-in-flow approximation cannot be used. Consequently, this four spacecraft approach is the only viable method for obtaining the wave number properties of the ambient fluctuations. Our results using electron velocity moments will be compared with previous work using magnetometer data from the FGM experiment on Cluster.

Goldstein, Melvyn L.

Internal temperature gradient of alloy semiconductor melts from interrupted growth experiments

Interrupted growth experiments on Hg(0.8)Cd(0.2)Se alloys were performed in an attempt to estimate the actual internal temperature gradient in the melt during directional solidification. The results have been analyzed in terms of a theoretical model which assumes an axial composition profile characteristic of diffusion controlled solute redistribution during growth. A comparison of the calculated and measured values suggests a reduction in the applied (empty furnace) temperature gradient by about a factor of three. Values of the interface segregation coefficient (k) determined from the phase diagram, predict values of solute concentration build up in the solid during growth interruptions which are inconsistent with measured results. It appears that the k values used tend to underestimate the actual k values for higher alloy compositions and overestimate k for the lower alloy compositions.

Andrews, R. N.

Estimate of the distance to the solar wind termination shock from gradients of anomalous cosmic ray oxygen

The radial gradient of anomalous cosmic ray oxygen measured instantaneously between Voyager 2 and Pioneer 10 during 1985-1988 is correlated with the inferred tilt of the heliospheric neutral sheet. This is consistent with a simple model in which the radial gradient is related to the length of the neutral sheet between the two spacecraft. With this model we show that the radial gradient and the tilt of the neutral sheet near the solar wind termination shock can be inferred from the Voyager and Pioneer observations. By comparing the time history of the inferred tilt with that derived from solar observations, we estimate that the termination shock was at 62 +5/-7 AU at solar minimum in 1987. At solar maximum the shock should be located at about 90 AU due to the increased pressure of the solar wind.

Cummings, A. C.