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 307 records · Page 17

Equation of State of CO 2 Shock Compressed to 1 TPa

Equation-of-state (pressure, density, temperature, internal energy) and reflectivity measurements on shock-compressed CO 2 at and above the insulating-to-conducting transition reveal new insight into the chemistry of simple molecular systems in the warm-dense-matter regime. CO 2 samples were precompressed in diamond-anvil cells to tune the initial densities from 1.35 g/cm 3 (liquid) to 1.74 g/cm 3 (solid) at room temperature and were then shock compressed up to 1 TPa and 93 000 K. Lastly, variation in initial density was leveraged to infer thermodynamic derivatives including specific heat and Gruneisen coefficient, exposing a complex bonded and moderately ionized state at the most extreme conditions studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance Study of CXL Memory Topology

This paper presents a comprehensive evaluation of the performance impact of various Compute Express Link (CXL) memory topologies, with a particular emphasis on CXL switches, in the context of High- Performance Computing (HPC) and Large Language Model (LLM) inference workloads. Our study unveils significant performance variations across different topologies, demonstrating that certain configurations yield superior performance for specific workloads. These findings underscore the critical importance of tailored topol- ogy selection in optimizing system performance. Additionally, we address the inherent challenges associated with integrating CXL switches, including overhead considerations and routing complex- ities. Our research highlights the necessity for thorough evalua- tion methodologies to fully leverage CXL technology’s potential in contemporary computing environments. These insights provide valuable guidance for system architects and data center operators in designing and optimizing CXL-based infrastructures for diverse workload requirements.

CXL, memory, Artificial Intelligence (AI), HPC↗

Spatiotemporal Evolution of Hanle and Zeeman Synthetic Polarization in a Chromospheric Spectral Line

Due to the quick evolution of the solar chromosphere, its magnetic field cannot be inferred reliably without accounting for the temporal variations of its polarized light. This has been broadly overlooked in the modeling and interpretation of the polarization, due to technical problems (e.g., lack of temporal resolution or of time-dependent MHD solar models) and/or because many polarization measurements can apparently be explained without dynamics. Here, we show that the temporal evolution is critical for explaining the spectral-line scattering polarization because of its sensitivity to rapidly varying physical quantities and the possibility of signal cancellations and attenuation during extended time integration. For studying the combined effect of time-varying magnetic fields and kinematics, we solved the 1.5D non-LTE problem of the second kind in time-dependent 3D R-MHD solar models and synthesized the Hanle and Zeeman polarization in forward scattering for the chromospheric λ 4227 line. We find that the quiet-Sun polarization amplitudes depend on the periodicity and spectral coherence of the signal enhancements produced by kinematics, but that substantially larger linear polarization signals should exist all over the solar disk for short integration times. The spectral morphology of the polarization is discussed as a combination of Hanle, Zeeman, partial redistribution and dynamic effects. We give physical references for observations by degrading and characterizing our slit time series in different spatiotemporal resolutions. The implications of our results for the interpretation of the second solar spectrum and for the investigation of the solar atmospheric heatings are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Coincident observations of ionospheric troughs and the equatorial plasmapause

Electron-density observations made in the topside ionosphere by the Ariel 4 and Isis 2 satellites are examined in conjunction with results obtained by Explorer 45 when it traversed the near-equatorial plasmapause with one hour (both UT and MLT) of the Ariel and Isis traversals of the same L coordinate. Both dusk and night observations are analyzed, and an attempt is made to show that depressions in ionospheric electron density occur in the vicinity of the plasmapause field line. It is concluded that the electron distributions observed in the electron-density troughs at 550 km near dusk by Ariel and at 1400 km near midnight by Isis do not always parallel variations in the light-ion distribution inferred from the Explorer plasmapause traversals and that there appears to be no specific feature of the main ionospheric trough which can be used to identify the plasmapause field line except in a statistical sense.

Grebowsky, J. M.↗

Chemistry of Apollo 12 mare basalts - Magma types and fractionation processes

Major and trace element data for a large suite of petrographically diverse Apollo 12 mare basalts are presented, and magma types sampled at the Apollo 12 site are characterized. The data confirm earlier classifications of the basalts into olivine, pigeonite, ilmenite, and feldspathic basalts. The olivine and pigeonite basalts are shown to be comagnetic and related by olivine fractionation. The other types, which differ in trace element and isotopic characteristics, are derived from different sources within the lunar interior. The spatial relations between the main basalt types are discussed in terms of local cratering events, and it is suggested that the younger ilmenite basalts overlie the olivine-pigeonite basalts. The role of olivine-dominated near-surface crystal fractionation in causing chemical variation is examined, and a relation between inferred cooling rate and the position of a sample in the fractionation sequence is determined.

Rhodes, J. M.↗

High-altitude atomic nitrogen densities

Theoretical calculations of the seasonal and diurnal variations of atomic nitrogen are compared with measurements made by the open source neutral mass spectrometer on the AE-C satellite. With the simultaneous measurements of molecular nitrogen and atomic oxygen densities as input, model calculations of odd nitrogen densities predict the same trends in atomic nitrogen as those observed. From these comparisons it is inferred that horizontal transport significantly reduces the diurnal variation of atomic nitrogen. Estimates are given of the sensitivity of atomic nitrogen densities to variations in the photoelectron flux, the neutral temperatures, and the neutral winds.

Oran, E. S.↗

Spatial dependences in the distant solar wind: Pioneers 10 and 11

Pioneer 10, 11 observations of the solar wind and magnetic field between 1 and 20 AU are reviewed. Spatial dependences, which are emphasized, must be inferred in the presence of large temporal variations including solar cycle effects. The separation of spatial and temporal dependences is achieved principally through the use of multipoint observations including baseline measurements at 1 AU. Measurements of the solar wind parameters (radial speed, flux, proton temperature) and of the magnetic field magnitude and components are compared with two theories, the Parker theory which assumes radial, azimuthally symmetric flow and the Goldstein-Jokipii theory which includes effects associated with stream-stream interactions. The observed radial gradients in the proton density and velocity and the magnetic field are consistent with the Parker model. A qualitative dependence of field magnitude on heliomagnetic latitude, i.e., referred to the observed location of the heliospheric current sheet, was derived. The field strength was found to decrease with distance from the current sheet.

Smith, E. J.↗

Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities

Here we propose a high dimensional Bayesian inference framework for learning heterogeneous dynamics of a COVID-19 model, with a specific application to the dynamics and severity of COVID-19 inside and outside long-term care (LTC) facilities. We develop a heterogeneous compartmental model that accounts for the heterogeneity of the time-varying spread and severity of COVID-19 inside and outside LTC facilities, which is characterized by time-dependent stochastic processes and time-independent parameters in ~ 1500 dimensions after discretization. To infer these parameters, we use reported data on the number of confirmed, hospitalized, and deceased cases with suitable post-processing in both a deterministic inversion approach with appropriate regularization as a first step, followed by Bayesian inversion with proper prior distributions. To address the curse of dimensionality and the ill-posedness of the high-dimensional inference problem, we propose use of a dimension-independent projected Stein variational gradient descent method, and demonstrate the intrinsic low-dimensionality of the inverse problem. We present inference results with quantified uncertainties for both New Jersey and Texas, which experienced different epidemic phases and patterns. Moreover, we also present forecasting and validation results based on the empirical posterior samples of our inference for the future trajectory of COVID-19.

42 ENGINEERING↗

Brittleness and Bayesian Inference

In a recent report, “Bayesian Brittleness: Why no Bayesian model is ‘good enough,’” (arXiv, 1304.6772v1) the authors prove rigorously that under a seemingly natural class of priors Π, the posterior expectation of a function Φ(θ) can achieve essentially any value that Φ(θ) alone can achieve, so that Bayesian inference appears to have no robustness whatsoever. We explain this puzzling result, and show why it does not imply a breakdown of Bayesian inference. The problem is that the priors leading to the extreme results depend on the observed data, which means they are not valid priors. The corresponding posteriors are also invalid, so the extreme variation in these posteriors does not imply that Bayesian inference is brittle. The report, however, provides a new framework which may be used for any assumed class of priors Π, and which we explain in detail. This framework may be a useful new tool for deriving quantitative robustness results.

97 MATHEMATICS AND COMPUTING↗

Time spectral analysis of midlatitude disturbances in the Martian atmosphere

Time-spectral and cross-spectral analyses are used to investigate synoptic pressure variations in the Viking-2 pressure, wind, and temperature data for selected portions of the Mars fall, winter, and spring seasons. Estimates of the phase relationships between the highly coherent pressure, wind, and temperature oscillations are obtained, and found to be very similar to those expected for baroclinic waves and to those obtained from terrestrial surface data. Phase speeds and zonal wavenumbers are inferred by interpreting the pressure and meridional wind variations in terms of eastward-traveling quasi-geostrophic waves.

Barnes, J. R.↗

Persistent Sub-Yearly Chromospheric Variations in Lower Main-Sequence Stars: Tau Booe and alpha Com

The recent discoveries of extrasolar planetary systems around lower main-sequence stars such as tau Booe (HD 120136) has prompted further investigation into their stellar activity. A cursory analysis of tau Booe for cyclic chromospheric activity, based on its 30-yr record of Ca 2 H and K fluxes obtained as part of the HK Project from Mount Wilson Observatory, finds an intermediate, sub-yearly period (approximately 117 d) in chromospheric activity in addition to, and separate from, both its rotation (3.3 d) and long-term variability. As a persistent subyearly period in surface magnetic activity is unprecedented, we investigate this apparent anomaly further by examining chromospheric activity levels of other stars with similar mass, searching for variability in chromospheric activity with periods of less than one year, but longer than measured or predicted rotation. An examination of the time series of 40 mid-to-late F dwarfs yielded one other star for further analysis: alpha Com (HD 114378, P approximately 132 d). The variations for these two stars were checked for persistence and coherence. Based on these determinations, we eliminate the possibilities of rotation, long-term activity cycle, and the evolution of active regions as the cause of this variation in both stars. In particular, for tau Booe we infer that the phenomenon may be chromospheric in origin; however, beyond this, it is difficult to identify anything further regarding the cause of the activity variations, or even whether the observed modulation in the two stars have the same origin.

Maulik, Davesh↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

Variational Bayesian strategies for high-dimensional, stochastic design problems

This paper is concerned with a lesser-studied problem in the context of model-based, uncertainty quantification (UQ), that of optimization/design/control under uncertainty. The solution of such problems is hindered not only by the usual difficulties encountered in UQ tasks (e.g. the high computational cost of each forward simulation, the large number of random variables) but also by the need to solve a nonlinear optimization problem involving large numbers of design variables and potentially constraints. We propose a framework that is suitable for a class of such problems and is based on the idea of recasting them as probabilistic inference tasks. To that end, we propose a Variational Bayesian (VB) formulation and an iterative VB–Expectation-Maximization scheme that is capable of identifying a local maximum as well as a low-dimensional set of directions in the design space, along which, the objective exhibits the largest sensitivity. We demonstrate the validity of the proposed approach in the context of two numerical examples involving thousands of random and design variables. In all cases considered the cost of the computations in terms of calls to the forward model was of the order of 100 or less. The accuracy of the approximations provided is assessed by information-theoretic metrics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗

Curie isotherm surfaces inferred from high-altitude magnetic anomaly data

Two-dimensional Curie depth models along two sections across the western United States are derived using an equivalent layer magnetization model derived from Magsat data, and the results are used to constrain finite element thermal models developed by Mayhew and Majer (1980). Regional heat flow variations predicted by the models compare favorably with those inferred from direct measurements. The methodology is applied to upward-continued aeromagnetic data, and the results are found to be in good agreement with previous Curie depth estimates and average measured heat flow.

Mayhew, M. A.↗

TIME EVOLUTION OF PLASMA PARAMETERS DURING THE RISE OF A SOLAR PROMINENCE INSTABILITY

We present high-spatial resolution spectropolarimetric observations of a quiescent hedgerow prominence taken in the He I 1083.0 nm triplet. The observation consisted of a time series in sit-and-stare mode of ∼36 minutes duration. The spectrograph's slit crossed the prominence body and we recorded the time evolution of individual vertical threads. Eventually, we observed the development of a dark Rayleigh-Taylor plume that propagated upward with a velocity, projected onto the plane of the sky, of 17 km s{sup –1}. Interestingly, the plume apex collided with the prominence threads pushing them aside. We inferred Doppler shifts, Doppler widths, and magnetic field strength variations by interpreting the He I Stokes profiles with the HAZEL code. The Doppler shifts show that clusters of threads move coherently while individual threads have oscillatory patterns. Regarding the plume we found strong redshifts (∼9-12 km s{sup –1}) and large Doppler widths (∼10 km s{sup –1}) at the plume apex when it passed through the prominence body and before it disintegrated. We associate the redshifts with perspective effects while the Doppler widths are more likely due to an increase in the local temperature. No local variations of the magnetic field strength associated with the passage of the plume were found; this leads us to conclude that the plumes are no more magnetized than the surroundings. Finally, we found that some of the threads' oscillations are locally damped, what allowed us to apply prominence seismology techniques to infer additional prominence physical parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Geophysical inferences of thermal-chemical structures in the lower mantle

Lateral variations of the temperature field in the lower mantle have been reconstructed using new results in mineral physics and seismic tomographic data. We show that, with the application of high-pressure experimental values of thermal expansivity and of sound velocities, the slow seismic anomalies in the lower mantle under the Pacific and Africa can be converted into realistic-looking plume structures with large dimensions of 0(1000 km). The outer fringes of the plumes have an excess temperature of around 400 K. In the core of the plumes are found tonguelike structures with extremely high thermal anomalies. These values can exceed 1200 K and are too high to be explained on the basis of thermal anomalies alone. We suggest that these major plumes in the deep mantle may be driven by both thermal and chemical buoyancies or that enhanced conductive heat-transfer may be important there.

Yuen, D. A.↗