Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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 955 records · Page 53

Statistical and Scaling Properties of Remotely-Sensed Soil Moisture in Two Contrasting Domains in the North American Monsoon Region

Characterizing soil moisture (theta) variability is important for inferring high-resolution information from coarse estimates provided by remote sensors. In this study, we analyze the spatial variability and scale invariance of high-resolution theta estimates collected in two contrasting semiarid areas, Arizona (AZ) and Sonora (SON), during the Soil Moisture Experiment - North American Monsoon in 2004 (SMEX04- NAME). Results reveal that as the mean theta condition ( ) becomes drier, the spatial standard deviation becomes smaller in both domains. The coefficient of variation of theta decreases with in SON, but does not display a clear tendency with in AZ. We also found the presence of scale invariance and multifractality in the range of support scales from 51.2 km to 0.8 km for all soil moisture fields in the two regions. The multifractal properties of theta are clearly linked to in SON, while the relation is affected by more dispersion in AZ. We argue this is due to differences in the dynamic (rainfall) and static (vegetation) controls on theta in the two domains.

Mascaro, Giuseppe↗

Magnetospheric electric fields deduced from drifting whistler paths.

Results of a study in which the amplitude of the E-W component E sub w of the convection electric field in the nightside magnetosphere has been inferred from the observed cross-L motions of whistler ducts within the plasmasphere, and several ducts distributed over 1 to 2 earth radii in L space and over plus or minus deg 15 around the longitude of the Eights, Antarctica, whistler station have been tracked simultaneously. The method appears capable of resolving fluctuations in E sub w with period T equal to approximately 15 min and rms amplitude as low as 0.05 mV/m. For variations with T greater than 1 hour the method has a sensitivity of the order of 0.01 mV/m. Three case studies are presented, two of which illustrate convection activity associated with relatively isolated substorms. In these two cases E sub w reversed from westward to eastward for a period following the decay of substorm bay activity. In the third case the substorm bay activity was prolonged, and E sub w remained westward and at enhanced levels until local dawn. Evidence was found that, at least in a limited longitudinal sector, perturbing substorm E sub w fields can penetrate deep within the plasmasphere. In two of the case studies comparisons of E sub w and the interplanetary magnetic-field theta component show evidence of a possible relation based on brief (less than or equal to 1 hour) southward excursions but not on long preceding southward events.

Carpenter, D. L.↗

Millimeter continuum observations of Galactic center giant molecular cloud cores

Results are presented of observations of 1.3- and 0.8-mm continuum emission toward the cores of three Galactic center molecular clouds with ongoing massive star formation, Sagittarius B2, C, and D, which were made in order to study possible variations in the high-mass star formation rate per unit mass between the Galactic center and the disk. The luminosity-to-mass ratio, based on the mass estimates derived from the millimeter continuum emission, is used as a tracer of the high-mass star formation rate in GMC cores. The magnitude of errors involved in using millimeter continuum emission for determining the core mass is estimated through radiative transfer modeling. It is inferred from the present millimeter data, along with previously published far-infrared data, that the Sgr C and D cores are very similar in terms of mean dust optical depth and temperature. The luminosity-to-mass ratios derived for the Sgr C and D cores are found to be consistent with those of typical disk GMC cores with comparable far-infrared luminosities.

Lis, D. C.↗

Antarctic Rebound and the Time-Dependence of the Earth's Shape

Great strides have been made during the past 30 years in refining models of the last global glaciation. The refinements draw upon a vastly expanded relative sea level and sedimentary core record. Furthermore, we now possess a sharpened understanding of the mechanisms that drive climate changes associated with deglaciation. Some 15 years ago, using only 5.5 years of ranging data, analyses of the drift in LAGEOS I node acceleration was used to infer that postglacial rebound was responsible for a secular change in the Earth's ellipsoidal shape (Yoder et al., .1983]. Today there exists a wealth of geodynamics satellite orbit data that constrain the secular time-dependence of the Earth's shape and low order gravity field, which includes mass redistribution from present-day glacier and great ice sheet imbalance and from postglacial rebound. We have shown that an unambiguous determination of the secular variation in the Earth's pear shaped harmonic (l = 3, m = 0) might provide information that bears on the present-day mass balance of Antarctica. This issue is revisited in light of new constraints on glacial loading during the late-Pleistocene and Holocene. An especially critical issue for the interpretation of secular odd degree zonal harmonics, l = 3 to 7, is the timing and magnitude of the deglaciation of Antarctica from Last Glacial Maximum. We explore ways in which the recovery of secular variation in both zonal and non-zonal harmonics for l = 2 through 7 can improve constraints on both rebound and present-day ice sheet balance.

Ivins, Erik R.↗

Image of Fomalhaut Dust Ring at 350 Microns: The Relative Column Density Map Shows Pericenter-Apocenter Asymmetry

We have imaged the circumstellar disk of Fomalhaut at 350 mm wavelength, using SHARC II (Submillimeter High Angular Resolution Camera II) at the Caltech Submillimeter Observatory. The spatial resolution of the raw images (9") has been enhanced by a factor of 3 using the HiRes deconvolution procedure. We find that at this wavelength and signal-to-noise ratio (approx.12), the observed morphology is that of a simple inclined ring (i approx. 70 deg), with little or no other apparent structure--this is the first observation that shows clearly the ring morphology of the disk. We have combined our 350 mm data with Spitzer Space Telescope images at 24, 70, and 160 mm in order to estimate the two-dimensional spatial variation of relative column density ("tau map") using our DISKFIT procedure. The tau map is based on the following physical assumptions: (1) the wavelength variation of opacity is the same throughout the disk, (2) the radial variation of dust temperature is dictated by the energy balance of individual grains in the stellar radiation field, and (3) the vertical scale height of the disk follows a power-law radial variation. The results confirm the ringlike morphology but also show that the geometric center is displaced from the star by about 8 AU and that the ring has an apocentric enhancement of approximately 14% in integrated column density. If we interpret the displacement in terms of elliptical orbital motion due to gravitational perturbation by an unseen planet, then the implied forced eccentricity is ~0.06; dynamical modeling then predicts an apocentric density enhancement consistent with that inferred from the tau map.

Fomalhaut↗

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗

Physical characteristics of cometary dust from optical studies

Observations of the sunlight scattered and thermal emission from cometary dust, which may be used to infer the physical properties of the dust grains, are reviewed. Consideration is given to the observed wavelength dependence of the scattered light from cometary coma and tails, the average scattering function of the dust grains, the average grain Bond albedo, the polarization of the scattered light, and grain temperatures deduced from thermal infrared emission. The thermal properties of dust grains are illustrated for models based on magnetite or olivine grain materials, with consideration given to the variation of thermal properties with particle radius and heliocentric distance. Comparison of the models with observations indicates that a disordered or amorphous olivine composition can give a reasonable fit to the data for appropriate grain sizes and temperatures. The observations acquired are noted to indicate an optically important particle size of 1 micron, with silicate particles not larger than a few microns usually present although pure silicate grains can not be responsible for the thermal emission, and the cometary dust grains are most likely not spherical. Further observations needed in the infrared are indicated.

Hanner, M. S.↗

Information theory lateral density distribution for Earth inferred from global gravity field

Information Theory Inference, better known as the Maximum Entropy Method, was used to infer the lateral density distribution inside the Earth. The approach assumed that the Earth consists of indistinguishable Maxwell-Boltzmann particles populating infinitesimal volume elements, and followed the standard methods of statistical mechanics (maximizing the entropy function). The GEM 10B spherical harmonic gravity field coefficients, complete to degree and order 36, were used as constraints on the lateral density distribution. The spherically symmetric part of the density distribution was assumed to be known. The lateral density variation was assumed to be small compared to the spherically symmetric part. The resulting information theory density distribution for the cases of no crust removed, 30 km of compensated crust removed, and 30 km of uncompensated crust removed all gave broad density anomalies extending deep into the mantle, but with the density contrasts being the greatest towards the surface (typically + or 0.004 g cm 3 in the first two cases and + or - 0.04 g cm 3 in the third). None of the density distributions resemble classical organized convection cells. The information theory approach may have use in choosing Standard Earth Models, but, the inclusion of seismic data into the approach appears difficult.

Rubincam, D. P.↗

Short-term variability of nitrogen dioxide in the winter stratosphere

Measurements of limb radiance from the Solar Mesosphere Explorer (SME) satellite are used to infer the NO2 density at the 10- and 16-mbar pressure levels from January 1 to March 31 of 1982. A photochemical-dynamical model is developed using the presently accepted chemistry of the NO(x) family. The dynamical model produces isentropic trajectories which simulate the history of air parcels. From the trajectories the photochemical model calculates NO2 densities, which are compared to those observed by SME. Although the model generally reproduces the spatial and temporal variations rather well, some disagreement was noted for conditions of exceptionally low temperatures. Further analysis indicates that the temperature sensitivity of the N2O5 photolysis cross sections may be overestimated at low temperatures.

Zawodny, J. M.↗

The nature of cometary dust as determined from infrared observations

The infrared measurements of comets, the compositional information available from interplanetary dust particles (IDPs), and the recent results of flybys to Comet Halley can help in restricting the nature and composition of cometary dust models (c.f., Proceedings of the 20th ESLAB Symposium on Exploration of Halley's Comet, 1986). Researchers tried to incorporate some of these results into a coherent model to account for the observed cometary infrared emission. The presence of 10 and 3.4 micron features in Comet Halley (c.f. Bregman et al. 1987; Wickramasinghe and Allen 1986) indicated the presence of at least two components in the grain material, namely silicates and some form of amorphous carbon. These two components could reside in separate grains or may be parts of composite particles. Both these cases have been considered (see Krishna Swamy el a. 1988a, 1988b). In the absence of refractive index data for cometary analogs, the authors used the optical constants of olivine-rich lunar material 12009.48 (Perry et al. 1972) for the infrared region and that of alpha:C-H film for amorphous carbon (angus et al. 1986). For the visible region, a value of m = 1.38-0.39i was used for the silicates, and values published by Arakawa et al. (1985) were used for the amorphous carbon. These materials should give a representative behavior of the expected results. The model results were compared to observational data. The strength of the 3.4 micron and 10 micron features relative to the adjacent continuum, as well as the slope of the continuum between 2500 and 1250 cm(exp -1) (4 to 8 microns), were used as criteria for comparison. Model calculations with alpha approx. equals -3.5, and also the size distribution function inferred for Comet Halley, with a mass fraction (X) of silicate to amorphous carbon grains of about 40 to 1 can fit the data. A good match is obtained for the infrared spectra of Comets Halley and West from a 40 to 1 mixture of silicate and amorphous carbon grains with a a(exp -3.5) size distribution function. The results are consistent with compositional constraints provided by interplanetary dust particles (IPDs) and Halley flyby data. The variation of grain temperature with heliocentric distance appears to account for the major changes observed in cometary spectra.

Swamy, K. S. Krishna↗

Dark material in the polar layered deposits on Mars

The Martian polar layered deposits probably record periodic variations in the deposition of dust and ice caused by climate changes over the past 10 to 100 million years. The polar layered deposits appear to be the source of dark, saltating material that has been distributed over the surface of Mars, but the mechanisms by which this material is incorporated and eroded from the layered deposits is unknown. These mechanisms must be understood before the processes that formed and modified the layered deposits can be inferred and related to Martian climate changes. In summary, weathering of the Martian layered deposits by sublimation of water ice may account for the geologic relationships observed in the Martian polar regions. The non-volatile components of the layered deposits appears to consist mainly of bright red dust, with small amounts of dark dust or sand. Alternatively, dark sand-sized basaltic particles may occasionally be transported onto the layered deposits, forming thin layers. Once eroded, particles of either type may saltate to form the dark dunes found in both polar regions.

Herkenhoff, Ken↗

Current sheet Formation in a Conical Theta Pinch Faraday Accelerator with Radio-Frequency Assisted Discharge

The inductive formation of current sheets in a conical theta pinch FARAD (Faraday Accelerator with Radio-frequency Assisted Discharge) thruster is investigated experimentally with time-integrated photography. The goal is to help in understanding the mechanisms and conditions controlling the strength and extent of the current sheet, which are two indices important for FARAD as a propulsion concept. The profiles of these two indices along the inside walls of the conical acceleration coil are assumed to be related to the profiles of the strength and extent of the luminosity pattern derived from photographs of the discharge. The variations of these profiles as a function of uniform back-fill neutral pressure (with no background magnetic field and all parameters held constant) provided the first clues on the nature and qualitative dependencies of current sheet formation. It was found that there is an optimal pressure for which both indices reach a maximum and that the rate of change in these indices with pressure differs on either side of this optimal pressure. This allowed the inference that current sheet formation follows a Townsend-like breakdown mechanism modified by the existence of a finite pressure-dependent radio-frequency-generated electron density background. The observation that the effective location of the luminosity pattern favors the exit-half of the conical coil is explained as the result of the tendency of the inductive discharge circuit to operate near its minimal self-inductance. Movement of the peak in the luminosity pattern towards the upstream side of the cone with increasing pressure is believed to result from the need of the circuit to compensate for the increase in background plasma resistivity due to increasing pressure.

Hallock, Ashley K.↗

Combined Heat Transfer in High-Porosity High-Temperature Fibrous Insulations: Theory and Experimental Validation

Combined radiation and conduction heat transfer through various high-temperature, high-porosity, unbonded (loose) fibrous insulations was modeled based on first principles. The diffusion approximation was used for modeling the radiation component of heat transfer in the optically thick insulations. The relevant parameters needed for the heat transfer model were derived from experimental data. Semi-empirical formulations were used to model the solid conduction contribution of heat transfer in fibrous insulations with the relevant parameters inferred from thermal conductivity measurements at cryogenic temperatures in a vacuum. The specific extinction coefficient for radiation heat transfer was obtained from high-temperature steady-state thermal measurements with large temperature gradients maintained across the sample thickness in a vacuum. Standard gas conduction modeling was used in the heat transfer formulation. This heat transfer modeling methodology was applied to silica, two types of alumina, and a zirconia-based fibrous insulation, and to a variation of opacified fibrous insulation (OFI). OFI is a class of insulations manufactured by embedding efficient ceramic opacifiers in various unbonded fibrous insulations to significantly attenuate the radiation component of heat transfer. The heat transfer modeling methodology was validated by comparison with more rigorous analytical solutions and with standard thermal conductivity measurements. The validated heat transfer model is applicable to various densities of these high-porosity insulations as long as the fiber properties are the same (index of refraction, size distribution, orientation, and length). Furthermore, the heat transfer data for these insulations can be obtained at any static pressure in any working gas environment without the need to perform tests in various gases at various pressures.

Daryabeigi, Kamran↗

Ionospheric plasma cloud dynamics

Measurements of the thermospheric neutral wind and ionospheric drift made at Eglin AFB, Florida and Kwajalein Atoll are discussed. The neutral wind measurements at Eglin had little variation over a period of four years for moderate magnetic activity (Kp 4); the ionospheric drifts are small. Evidence is presented that indicates that increased magnetic activity has a significant effect on the neutral wind magnitude and direction at this midlatitude station. The neutral wind at dusk near the equator is generally small although in one case out of seven it was significantly larger. It is described how observations of large barium releases can be used to infer the degree of electrodynamic coupling of ion clouds to the background ionosphere. Evidence is presented that indicates that large barium releases are coupled to the conjugate ionosphere at midlatitudes.

Source record↗

On the determination of the Hermaean magnetic moment - A critical review

The determination of Mercury's magnetic moment from the spatially and temporally limited observations obtained by the Mariner 10 mission is dependent upon the assumed nature of both the intrinsic planetary magnetic field and that of the magnetospheric current systems. In this paper the methods that have been used for this purpose are reviewed. The results that have been obtained are then considered in comparison with the constraints placed on the planetary field by the solar wind dynamic pressures and stand-off distances inferred from the Mariner 10 magnetic field data by Slavin and Holtzer (1979) which are consistent with an effective dipole moment of 6 + or - 2 x 10 to the 22nd G cu cm. It is found that the models which fit the observations with substantial quadrupole and octupole moments are not consistent with the magnetospheric boundary conditions, presumably owing to incorrect assumptions regarding the magnetopause position, incorrect assumptions regarding solar wind dynamic pressure, and/or averaging over temporal variations in the Mariner 10 data.

Slavin, J. A.↗

Directional emittance corrections for thermal infrared imaging

A simple measurement technique for measuring the variation of directional emittance of surfaces at various temperatures using commercially available radiometric IR imaging systems was developed and tested. This technique provided the integrated value of directional emittance over the spectral bandwidth of the IR imaging system. The directional emittance of flat black lacquer and red stycast, an epoxy resin, measured using this technique were in good agreement with the predictions of the electromagnetic theory. The data were also in good agreement with directional emittance data inferred from directional reflectance measurements made on a spectrophotometer.

Daryabeigi, Kamran↗

Metallographic Cooling Rate of IVA Irons Revisited

There is long standing problem reconciling the chemical evidence that the IVA iron meteorites formed in a core with the diverse cooling rates reported by several researchers. This large inferred range of cooling rates suggests that the IVA irons were distributed at different depths in a parent body with a complex structure when the Widmanstatten pattern formed. On the other hand, some researchers argued that the diverse cooling rates in group IVA result from inaccurate model parameters such as phase diagram, interdiffusion coefficients, and kamacite nucleation and growth mechanisms. In addition, the measured cooling rates may not apply for the same cooling temperature ranges, and the variation in the crystallographic orientations of the Widmanstatten plates on the analysis surface may result in inaccurate measurements of widths needed for the computer simulation models. We have revaluated the major parameters in computer model developed by Hopfe and Goldstein and measured cooling rates for the IVA irons. Such data are useful in evaluating whether these meteorites were part of a single core of a parent body during the formation of the Widmanstatten pattern.

Yang, J.↗

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

The triple Langmuir probe is an electrostatic probe in which three probe tips collect current when inserted into a plasma. The triple probe differs from a simple single Langmuir probe in the nature of the voltage applied to the probe tips. In the single probe, a swept voltage is applied to the probe tip to acquire a waveform showing the collected current as a function of applied voltage (I-V curve). In a triple probe three probe tips are electrically coupled to each other with constant voltages applied between each of the tips. The voltages are selected such that they would represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Under the assumption of a Maxwellian plasma, one can determine the time-varying plasma temperature T(sub e)(t) and number density n(sub e)(t) from the applied voltage levels and the time-histories of the collected currents. In the present paper we examine the theory of triple probe operation, specifically focusing on the assumption of a Maxwellian plasma. Triple probe measurements have been widely employed for a number of pulsed and timevarying plasmas, including pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. While the equilibrium assumption may be justified for some applications, it is unlikely that it is fully justifiable for all pulsed and time-varying plasmas or for all times during the pulse of a plasma device. To examine a simple non-equilibrium plasma case, we return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. We compare the distribution function that was assumed a priori with the distribution function one would infer when applying standard triple probe theory to analyze the collected currents. For the assumed class of non-Maxwellian distribution functions this serves to illustrate the effect a non-Maxwellian plasma would have on results interpreted using the equilibrium triple probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗