Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 1,441 records · Page 80

Intrepid MCMC: Metropolis-Hastings with exploration

In engineering examples, one often encounters the need to sample from unnormalized distributions with complex shapes that may also be implicitly defined through a physical or numerical simulation model, making it computationally expensive to evaluate the associated density function. For such cases, MCMC has proven to be an invaluable tool. Random-walk Metropolis Methods (also known as Metropolis-Hastings (MH)), in particular, are highly popular for their simplicity, flexibility, and ease of implementation. However, most MH algorithms suffer from significant limitations when attempting to sample from distributions with multiple modes (particularly disconnected ones). Here, in this paper, we present Intrepid MCMC - a novel MH scheme that utilizes a simple coordinate transformation to significantly improve the mode-finding ability and convergence rate to the target distribution of random-walk Markov chains while retaining most of the simplicity of the vanilla MH paradigm. Through multiple examples, we showcase the improvement in the performance of Intrepid MCMC over vanilla MH for a wide variety of target distribution shapes. We also provide an analysis of the mixing behavior of the Intrepid Markov chain, as well as the efficiency of our algorithm for increasing dimensions. A thorough discussion is presented on the practical implementation of the Intrepid MCMC algorithm. Finally, its utility is highlighted through a Bayesian parameter inference problem for a two-degree-of-freedom oscillator under free vibration.

97 - MATHEMATICS AND COMPUTING↗

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗

Visibility related to backscatter at 1.54 micron

The lidar process was shown to have the necessary potential to fulfill the need for a remote measurement of visibility. Visibility can be inferred from a lidar return optical extinction. The wavelength 1.54 micron was chosen, being near the visible wavelength region and having a high eye safety threshol, 200,000 times higher than 1.06 micron; 1.54 is the erbium laser wavelength. This research utilized 105 measured height profiles of natural droplet size distributions data, taken in clouds, fog, and haze. These profiles were examined to determine the completeness of the droplet counting data. It was found that the particle spectrometer data were incomplete in the very light ford and haze so this portion of the data was eliminated. Utilizing the Mie theory, these droplet size distribution profiles were converted to backscatter at 1.54 micron and extinction in the visible region, 0.55 micron. Using Koschmeider's relationship, the extinction profiles were converted to visibility. The visibility and backscatter profiles were compared to develop a relationship between visibility and backscatter at 1.54 micron.

Barber, T. L.↗

On convection-induced viscosity in accretion disks in cataclysmic variables

The vertical structure of accretion disks in cataclysmic variable systems is examined on the assumption that the convection is solely responsible for inducing the viscous heating. The viscous heating is distributed over the entire vertical structure so that purely radiative zones are able to carry a radial mass flux. The convective flux is computed by two different methods in separate calculations and the results from each method are compared. The values of the viscosity obtained are compared with those values inferred from the observed time scales of dwarf nova eruptions. Sequences of equilibrium solutions for constant radius and varying mass accretion rate are examined, the relations between viscosity and surface density are discussed in the context of recent models for the eruptions of dwarf novae.

Cannizzo, John K.↗

Derivation of Delaware Bay tidal parameters from Space Shuttle photography

The tide-related parameters of the Delaware Bay are derived from Space Shuttle time-series photographs. The water areas in the bay are measured from interpretation maps of the photographs with a CALCOMP 9100 digitizer and ERDAS Image Processing System. The corresponding tidal levels are calculated using the exposure time annotated on the photographs. From these data, an approximate function relating the water area to the tidal level at a reference point is determined. Based on the function, the water areas of the Delaware Bay at mean high water (MHW) and mean low water (MLW), below 0 m, and for the tidal zone are inferred. With MHW and MLW areas and the mean tidal range, we calculate the tidal influx of the Delaware Bay, which is 2.76 x 10 exp 9 cu m. The velocity of flood tide at the bay mouth is determined using the tidal flux and an integral of the velocity distribution function at the cross section between Cape Henlopen and Cape May. The result is 132 cm/s, which compares well with the data on tidal current charts.

Zheng, Quanan↗

Recent deformation rates on Venus

Constraints on the recent geological evolution of Venus may be provided by quantitative estimates of the rates of the principal resurfacing processes, volcanism and tectonism. This paper focuses on the latter, using impact craters as strain indicators. The total postimpact tectonic strain lies in the range 0.5-6.5%, which defines a recent mean strain rate of 10(exp -18)-10(exp -17)/s when divided by the mean surface age. Interpretation of the cratering record as one of pure production requires a decline in resurfacing rates at about 500 Ma (catastrophic resurfacing model). If distributed tectonic resurfacing contributed strongly before that time, as suggested by the widespread occurrence of tessera as inliers, the mean global strain rate must have been at least approximately 10(exp -15)/s, which is also typical of terrestrial active margins. Numerical calculations of the response of the lithosphere to inferred mantle convective forces were performed to test the hypothesis that a decrease in surface strain rate by at least two orders of magnitude could be caused by a steady decline in heat flow over the last billion years. Parameterized convection models predict that the mean global thermal gradient decreases by only about 5 K/km over this time; even with the exponential dependence of viscosity upon temperature, the surface strain rate drops by little more than one order of magnitude. Strongly unsteady cooling and very low thermal gradients today are necessary to satisfy the catastrophic model. An alternative, uniformitarian resurfacing hypothesis holds that Venus is resurfaced in quasi-random 'patches' several hundred kilometers in size that occur in response to changing mantle convection patterns.

Grimm, Robert E.↗

The Influence of Planetary Wave Transport on Arctic Ozone as Observed by POAM III

Three recent Arctic winters have been unusually warm or cold compared to the 1979-2001 mean, resulting in either a warm, disturbed vortex or a cold, quiet vortex. Interannual differences n ozone transport in the lower stratospheric Arctic vortex are investigated using data from the Polar Ozone and Aerosol Measurement (POAM) III satellite instrument obtained between May 1998 and April 2001. Time series of probability distribution functions (PDFs) of the FOAM data are used to identify seasonal and interannual variations in the transport processes controlling polar ozone. A major warming occurred in December, 1998, that caused the middle stratospheric vortex to be displaced from the pole for nearly a month and dissipate. Ozone transport is inferred from the PDFs and supported by potential vorticity (PV) analyses. The vortex reformed in January, 1999, filled with high O3 air from lower latitudes, which then cooled and descended. As a result, by the end of winter 1999, ozone at 500K in the vortex was 0.5-1.0 ppm higher than in the cold winter of 2000. The winter of 2000-1 had frequent wave disturbances, beginning with a fairly large event in November. However, at the end of this wave activity, the high PV core of the vortex was still intact, and by February, 03 in the vortex looked the same as in 2000. This study demonstrates that interannual variability in the 03-PV relationship can be caused by transport, independent of 03 loss by polar stratospheric clouds (PSCs). Ozone in the vortex in a very disturbed winter does not necessarily represent'pre-PSC' O3 levels in other years.

Strahan, Susan↗

Selective Amnesia using Contrastive Subnet Erasure for Class Level Unlearning in Vision Models

We study concept-level forgetting in pretrained vision models: removing an entire semantic category so the system no longer recognizes that object in unseen images and contexts, rather than merely forgetting specific training examples. Prior work either applies blunt global projections or fine-tunes parameters, which can introduce collateral damage to unrelated features, add compute, and become unstable as forgetting strength increases. We introduce Contrastive Subnet Erasure (CSE), a training-free, encoder-centric edit that targets a compact set of channels most responsible for the class and attenuates them in a calibrated manner. The modification is algebraically folded into the subsequent layer, yielding no inference-time overhead and leaving task heads unchanged. To evaluate whether forgetting generalizes beyond the data used to specify the class, we introduce a cross dataset protocol in which the class is defined on a source dataset and performance is measured on a disjoint target dataset drawn from a different distribution with no shared images. This setup tests whether the model still fails to recognize the object when it looks different or appears in new scenes, and it helps avoid overfitting to patterns in the source dataset. Across CIFAR 10, CIFAR 100, and ImageNet under this protocol, CSE achieves stronger forgetting of the target class while better preserving non target utility than existing baselines in both single class and multi class settings. Overall, CSE provides a simple, stable, and deployment-ready mechanism for class-level unlearning in vision.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

A genetic algorithm approach for reconstructing spectral content from filtered x-ray diode array spectrometers

Filtered diode array spectrometers are routinely employed to infer the temporal evolution of spectral power from x-ray sources, but uniquely extracting spectral content from a finite set of broad, spectrally overlapping channel spectral sensitivities is decidedly nontrivial in these under-determined systems. Here, we present the use of genetic algorithms to reconstruct a probabilistic spectral intensity distribution and compare to the traditional approach most commonly found in the literature. Unlike many of the previously published models, spectral reconstructions from this approach are neither limited by basis functional forms nor do they require a priori spectral knowledge. While the original intent of such measurements was to diagnose the temporal evolution of spectral power from quasi-blackbody radiation sources—where the exact details of spectral content were not thought to be crucial—we demonstrate that this new technique can greatly enhance the utility of the diagnostic by providing more physical spectra and improved robustness to hardware configuration for even strongly non-Planckian distributions.

47 OTHER INSTRUMENTATION↗

Whole-genome sequencing distinguishes the two most common giant kelp ecomorphs

Abstract Giant kelp, Macrocystis pyrifera, exists as distinct morphological variants—or “ecomorphs”—in different populations, yet the mechanism for this variation is uncertain, and environmental drivers for either adaptive or plastic phenotypes have not been identified. The ecomorphs Macrocystis “pyrifera” and M. “integrifolia” are distributed throughout temperate waters of North and South America with almost no geographic overlap and exhibit an incongruous, non-mirrored, distribution across the equator. This study evaluates the degree of genetic divergence between M. “pyrifera” and M. “integrifolia” across 18 populations in Chile and California using whole-genome sequencing and single-nucleotide polymorphism markers. Our results based on a principal component analysis, admixture clustering by genetic similarity, and phylogenetic inference demonstrate that M. “pyrifera” and M. “integrifolia” are genetically distinguishable. Analyses reveal separation by Northern and Southern Hemispheres and between morphs within hemispheres, suggesting that the convergent “integrifolia” morphology arose separately in each hemisphere. This is the first study to use whole-genome sequencing to understand genetic divergence in giant kelp ecomorphs, identifying 83 potential genes under selection and providing novel insights about Macrocystis evolution that were not evident with previous genetic techniques. Future studies are needed to uncover the environmental forces driving local adaptation and presumed convergent evolution of these morphs.

Environmental Sciences & Ecology↗

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗

High latitude minor ion enhancements: A clue for studies of magnetosphere-atmosphere coupling

Unexpectedly abrupt and pronounced distributions of the thermal molecular ions NO(+), O2(+) and N2(+) were observed at mid and high latitudes by the OGO-6 ion mass spectrometer. These minor ions may reach concentration levels exceeding 1000 ions/cu cm at altitudes as great as 1000 km, suggestive of scale heights well in excess of those inferred from low and mid-latitude measurements, under relatively undisturbed conditions. The high latitude ion enhancements were observed to be narrowly defined in time and space, with molecular ion concentrations changing by as much as an order of magnitude between successive orbits.

Taylor, H. A., Jr.↗

Interpretation of spectral optical thickness observations of the El Chichon eruption cloud

It is pointed out that the recent El Chichon eruption has provided a unique opportunity to observe a short-term climate change event. In October 1982 and in May 1983 observation experiments were undertaken to determine the spectral optical thickness of the El Chichon eruption cloud during airborne latitudinal surveys. One of two basic objectives to the observation experiment was the conduction of direct measurements of the optical thickness for radiation and climate studies. In connection with the second objective, inferences of the physical characteristic of the volcanic cloud were to be obtained on the basis of an application of remote sensing techniques. The present paper is concerned with the second objective, taking into account questions regarding the feasibility to conduct a meaningful analysis of the variation of the aerosol size distribution with latitude on the basis of the accuracy of the measured data.

Spinhirne, J. D.↗

Cratering and obliteration history of the south polar region of Mars

Results are presented from a comprehensive analysis of crater size frequency distributions compiled from Viking Orbiter images of Martian south polar terrains. Depositional histories for the various terrain units are modelled based on deviations of cumulative crater size frequency plots from an assumed production function. Stratigraphic and morphologic data obtained from visual examination of Viking images are also used, primarily to corroborate inferences concerning ages and obliteration histories derived from the crater data.

Plaut, J. J.↗

The cratering record in the inner solar system: Implications for earth

Internal and external processes have reworked the Earth's surface throughout its history. In particular, the effect of meteorite impacts on the early history of the earth is lost due to fluvial, aeolian, volcanic and plate tectonic action. The cratering record on other inner solar system bodies often provides the only clue to the relative cratering rates and intensities that the earth has experienced throughout its history. Of the five major bodies within the inner solar system, Mercury, Mars, and the Moon retain scars of an early episode of high impact rates. The heavily cratered regions on Mercury, Mars, and the Moon show crater size-frequency distribution curves similar in shape and crater density, whereas the lightly cratered plains on the Moon and Mars show distribution curves which, although similar to each other, are statistically different in shape and density from the more heavily cratered units. The similarities among crater size-frequency distribution curves for the Moon, Mercury, and Mars suggest that the entire inner solar system was subjected to the two populations of impacting objects but Earth and Venus have lost their record of heavy bombardment impactors. Thus, based on the cratering record on the Moon, Mercury, and Mars, it can be inferred that the Earth experienced a period of high crater rates and basin formation prior to about 3.8 BY ago. Recent studies have linked mass extinctions to large terrestrial impacts, so life forms were unable to establish themselves until impact rates decreased substantially and terrestrial conditions became more benign. The possible periodicity of mass extinctions has led to the theory of fluctuating impact rates due to comet showers in the post heavy bombardment period. The active erosional environment on the Earth complicates attempts to verify these showers by erasing geological evidence of older impact craters. The estimated size of the impactor purportedly responsible for the Cretaceous-Tertiary mass extinctions is 10 km in diameter. Thus impactors greater than or equal to the size postulated for K-T impactor are rare within the inner solar system since the end of heavy bombardment.

Barlow, N. G.↗

Global evolution of interplanetary sector structure, coronal holes, and solar wind streams during 1976-1993: Stackplot displays based on solar magnetic observations

We use potential field calculations and solar magnetic observations during 1976 to 1993 to infer the evolution of interplanetary sector structure, coronal holes, and solar wind streams at heliographic latitudes ranging from 80 deg S to 80 deg N. The results are presented in the form of stackplots, which show long-lived patterns that rotate quasi-rigidly at rates determined by the photospheric distribution of nonaxisymmetric magnetic flux. The fastest wind streams and their coronal hole sources form slowly rotating patterns near the poles just after sunspot maximum but migrate to lower latitudes and tend to rotate at near-equatorial rates as sunspot activity declines.

Wang, Y. M.↗