Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data 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,639 records · Page 91

Uncertainties for two-dimensional models of solar rotation from helioseismic eigenfrequency splitting

Observed solar p-mode frequency splittings can be used to estimate angular velocity as a function of position in the solar interior. Formal uncertainties of such estimates depend on the method of estimation (e.g., least-squares), the distribution of errors in the observations, and the parameterization imposed on the angular velocity. We obtain lower bounds on the uncertainties that do not depend on the method of estimation; the bounds depend on an assumed parameterization, but the fact that they are lower bounds for the 'true' uncertainty does not. Ninety-five percent confidence intervals for estimates of the angular velocity from 1986 Big Bear Solar Observatory (BBSO) data, based on a 3659 element tensor-product cubic-spline parameterization, are everywhere wider than 120 nHz, and exceed 60,000 nHz near the core. When compared with estimates of the solar rotation, these bounds reveal that useful inferences based on pointwise estimates of the angular velocity using 1986 BBSO splitting data are not feasible over most of the Sun's volume. The discouraging size of the uncertainties is due principally to the fact that helioseismic measurements are insensitive to changes in the angular velocity at individual points, so estimates of point values based on splittings are extremely uncertain. Functionals that measure distributed 'smooth' properties are, in general, better constrained than estimates of the rotation at a point. For example, the uncertainties in estimated differences of average rotation between adjacent blocks of about 0.001 solar volumes across the base of the convective zone are much smaller, and one of several estimated differences we compute appears significant at the 95% level.

Genovese, Christopher R.↗

Microwave Sensing of Ocean Winds Including the Use of SAR Imagery

This paper presents aircraft measurements of near-surface atmospheric boundary layer wind roll signatures and radar-derived sea surface roughness. These data are completely coincident in space and time and this unique feature supports attempts to definitively link SAR backscatter signatures to boundary layer roll impacts. The open-ocean data were collected at an altitude of 20 m from NOAA's Long-EZ aircraft using its turbulence probe and down-looking Ka-band radar scatterometer. Several flight legs of 20-30 km were flown with a heading across the wind direction, which is also roughly perpendicular to the roll vortices. We find remarkable correlation between measured modulations in the along-wind component of wind speed and radar backscatter for the spatial scale of 1 to 1.5 km. Close agreement between normalized modulation amplitudes suggests the radar-inferred surface slope variance is changing linearly with wind speed. These data were collected within 30 minutes of a RADARSAT SAR overpass where apparent boundary layer impacts of the same orientation and spatial dimension are prevalent in the SAR backscatter image. Quantitative comparison between modulations in the aircraft and satellite radar data will be discussed.

Vandemark, Douglas↗

Analysis of X-ray observations of the 15 June 1973 flare in active region NOAA 131

Observations and analyses of the 1B/M3 flare of 15 June, 1973 in active region NOAA 131 (McMath 12379) are presented. The X-ray observations, consisting of broadband photographs and proportional counter data from the Skylab/ATM NASA-MSFC/Aerospace S-056 experiment, are used to infer temperatures, emission measures, and densities for the flaring plasma. The peak temperature from the spatially resolved photographs is 25,000,000 K, while the temperature from the full-disk proportional counter data is approximately 15,000,000 K. The density is 3 times 10 to the 10th/cu cm. The X-ray flare emission appears to come primarily from two low-lying curvilinear features lying perpendicular to and centered on the line where the photospheric longitudinal magnetic field is zero. Similarities in the preflare and postflare X-ray emission patterns indicate that no large-scale relaxation of the coronal magnetic configuration was observed. Also discussed are H-alpha and magnetic field observations of the flare and the active region. Finally, results of numerical calculations, including thermal conduction, radiative loss, and chromospheric evaporation, are in qualitative agreement with the decay phase observations.

Krall, K. R.↗

Degree of Ice Particle Surface Roughness Inferred from Polarimetric Observations

The degree of surface roughness of ice particles within thick, cold ice clouds is inferred from multidirectional, multi-spectral satellite polarimetric observations over oceans, assuming a column-aggregate particle habit. An improved roughness inference scheme is employed that provides a more noise-resilient roughness estimate than the conventional best-fit approach. The improvements include the introduction of a quantitative roughness parameter based on empirical orthogonal function analysis and proper treatment of polarization due to atmospheric scattering above clouds. A global 1-month data sample supports the use of a severely roughened ice habit to simulate the polarized reflectivity associated with ice clouds over ocean. The density distribution of the roughness parameter inferred from the global 1- month data sample and further analyses of a few case studies demonstrate the significant variability of ice cloud single-scattering properties. However, the present theoretical results do not agree with observations in the tropics. In the extra-tropics, the roughness parameter is inferred but 74% of the sample is out of the expected parameter range. Potential improvements are discussed to enhance the depiction of the natural variability on a global scale.

surface↗

A biology-informed similarity metric for simulated patches of human cell membrane

Complex scientific inquiries rely increasingly upon large and autonomous multiscale simulation campaigns, which fundamentally require similarity metrics to quantify ‘sufficient’ changes among data and/or configurations. However, subject matter experts are often unable to articulate similarity precisely or in terms of well-formulated definitions, especially when new hypotheses are to be explored, making it challenging to design a meaningful metric. Furthermore, the key to practical usefulness of such metrics to enable autonomous simulations lies in in situ inference, which requires generalization to possibly substantial distributional shifts in unseen, future data. Here, we address these challenges in a cancer biology application and develop a meaningful similarity metric for ‘patches’—regions of simulated human cell membrane that express interactions between certain proteins of interest and relevant lipids. In the absence of well-defined conditions for similarity, we leverage several biology-informed notions about data and the underlying simulations to impose inductive biases on our metric learning framework, resulting in a suitable similarity metric that also generalizes well to significant distributional shifts encountered during the deployment. We combine these intuitions to organize the learned embedding space in a multiscale manner, which makes the metric robust to incomplete and even contradictory intuitions. Our approach delivers a metric that not only performs well on the conditions used for its development and other relevant criteria, but also learns key spatiotemporal relationships without ever being exposed to any such information during training.

97 MATHEMATICS AND COMPUTING↗

Utilizing Calibrated GPS Reflected Signals to Estimate Soil Reflectivity and Dielectric Constant: Results from SMEX02

Extensive reflected GPS data was collected using a GPS reflectometer installed on an HC130 aircraft during the Soil Moisture Experiment 2002 (SMEX02) near Ames, Iowa. At the same time, widespread surface truth data was acquired in the form of point soil moisture profiles, areal sampling of near-surface soil moisture, total green biomass and precipitation history, among others. Previously, there have been no reported efforts to calibrate reflected GPS data sets acquired over land. This paper reports the results of two approaches to calibration of the data that yield consistent results. It is shown that estimating the strength of the reflected signals by either (1) assuming an approximately specular surface reflection or (2) inferring the surface slope probability density and associated normalization constants give essentially the same results for the conditions encountered in SMEX02. The corrected data is converted to surface reflectivity and then to dielectric constant as a test of the calibration approaches. Utilizing the extensive in-situ soil moisture related data this paper also presents the results of comparing the GPS-inferred relative dielectric constant with the Wang-Schmugge model frequently used to relate volume moisture content to dielectric constant. It is shown that the calibrated GPS reflectivity estimates follow the expected dependence of permittivity with volume moisture, but with the following qualification: The soil moisture value governing the reflectivity appears to come from only the top 1-2 centimeters of soil, a result consistent with results found for other microwave techniques operating at L-band. Nevertheless, the experimentally derived dielectric constant is generally lower than predicted. Possible explanations are presented to explain this result.

Katzberg, Stephen J.↗

Estimating the R-curve from residual strength data

A method is presented for estimating the crack-extension resistance curve (R-curve) from residual-strength (maximum load against original crack length) data for precracked fracture specimens. The method allows additional information to be inferred from simple test results, and that information can be used to estimate the failure loads of more complicated structures of the same material and thickness. The fundamentals of the R-curve concept are reviewed first. Then the analytical basis for the estimation method is presented. The estimation method has been verified in two ways. Data from the literature (involving several materials and different types of specimens) are used to show that the estimated R-curve is in good agreement with the measured R-curve. A recent predictive blind round-robin program offers a more crucial test. When the actual failure loads are disclosed, the predictions are found to be in good agreement.

Orange, T. W.↗

Precipitation Characteristics of ISCCP Cloud Regimes for Improving Model Hydrological Budgets

The key in unraveling relationships between precipitation and atmospheric circulations is their common linkage to clouds. Clouds can be described in a variety of ways and several approaches can be adopted to examine their connections to precipitation. We claim that when cloud regimes (aka weather states) from the International Satellite Cloud Climatology Project (ISCCP) are used to conditionally sample/sort and average precipitation data, useful insight and GCM-appropriate diagnostics on the origins and distribution of precipitation can be obtained. The ISCCP cloud regimes are mesoscale (2.5 ) cloud mixtures determined by cluster analysis on joint histograms of cloud optical thickness and cloud top pressure inferred from geostationary and polar orbiter satellite passive retrievals. The ISCCP cloud regime data are combined with GPCP IDD merged surface precipitation data and/or higher temporal and spatial resolution TRMM Multi-Satellite Precipitation Analysis (TMPA) data. The analysis is performed separately for three geographical zones, tropics, and northern/southern midlatitudes (for GPCP; only the tropics can be examined with TMPA data). Our presentation aspires to provide answers to the following questions: (l) What is the mean and variability of surface precipitation produced by each cloud regime at the time of regime occurrence? (2) What is the relative contribution of each cloud regime to the total precipitation within its geographical zone? (3) What is the geographical distribution of precipitation corresponding to particular cloud regime? (4) To what extent are the cloud regimes distinct in terms of their precipitation characteristics and is the regime ordering in terms of convective strength consistent with the observed precipitation intensity?

Lee, D.↗

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Coronal structures deduced from photospheric magnetic field and He I lambda 10830 observations

The National Solar Observatory synoptic program provides an extensive and unique data base of high-resolution full-disk observations of the line-of-sight photospheric magnetic fields and of the He I lambda 10830 equivalent width. These data have been taken nearly daily for more than 21 years since 1974 and provide the opportunity to investigate the behavior of the magnetic fields in the photosphere and those inferred for the corona spanning on the time scales of a day to that of a solar cycle. The intensity of structures observed in He I lambda 10830 are strongly modulated by overlying coronal radiation; areas with low coronal emission are generally brighter in He I lambda 10830, while areas with high coronal emission are darker. For this reason, He I lambda 10830 was selected in the mid-1970's as way to identify and monitor coronal holes, magnetic fields with an open configuration, and the sources of high-speed solar wind streams. The He I lambda 10830 spectroheliograms also show a wide variety of other structures from small-scale, short-lived dark points (less than 30 arc-sec, hours) to the large-scale, long-lived two 'ribbon' flare events that follow the filament eruptions (1000 arc-sec, days). Such structures provide clues about the connections and changes in the large-scale coronal magnetic fields that are rooted in concentrations of magnetic network and active regions in the photosphere. In this paper, what observations of the photospheric magnetic field and He I lambda 10830 can tell us about the short- and long-term evolution of the coronal magnetic fields will be discussed, focussing on the quiet Sun and coronal holes. These data and what we infer from them will be compared with direct observations of the coronal structure from the Yohkoh Soft X-ray Telescope.

Harvey, Karen L.↗

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗

Systematic characterization of human gut microbiome-secreted molecules by integrated multi-omics

The human gut microbiome produces a complex mixture of biomolecules that interact with human physiology and play essential roles in health and disease. Crosstalk between micro-organisms and host cells is enabled by different direct contacts, but also by the export of molecules through secretion systems and extracellular vesicles. The resulting molecular network, comprised of various biomolecular moieties, has so far eluded systematic study. Here we present a methodological framework, optimized for the extraction of the microbiome-derived, extracellular biomolecular complement, including nucleic acids, (poly)peptides, and metabolites, from flash-frozen stool samples of healthy human individuals. Our method allows simultaneous isolation of individual biomolecular fractions from the same original stool sample, followed by specialized omic analyses. The resulting multi-omics data enable coherent data integration for the systematic characterization of this molecular complex. Our results demonstrate the distinctiveness of the different extracellular biomolecular fractions, both in terms of their taxonomic and functional composition. This highlights the challenge of inferring the extracellular biomolecular complement of the gut microbiome based on single-omic data. The developed methodological framework provides the foundation for systematically investigating mechanistic links between microbiome-secreted molecules, including those that are typically vesicle-associated, and their impact on host physiology in health and disease.

59 BASIC BIOLOGICAL SCIENCES↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Furthermore, our examples can be reproduced with the code publicly available at Zenodo.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quartet-based inference is statistically consistent under the unified duplication-loss-coalescence model

Abstract Motivation The classic multispecies coalescent (MSC) model provides the means for theoretical justification of incomplete lineage sorting-aware species tree inference methods. This has motivated an extensive body of work on phylogenetic methods that are statistically consistent under MSC. One such particularly popular method is ASTRAL, a quartet-based species tree inference method. Novel studies suggest that ASTRAL also performs well when given multi-locus gene trees in simulation studies. Further, Legried et al. recently demonstrated that ASTRAL is statistically consistent under the gene duplication and loss model (GDL). GDL is prevalent in evolutionary histories and is the first core process in the powerful duplication-loss-coalescence evolutionary model (DLCoal) by Rasmussen and Kellis. Results In this work, we prove that ASTRAL is statistically consistent under the general DLCoal model. Therefore, our result supports the empirical evidence from the simulation-based studies. More broadly, we prove that the quartet-based inference approach is statistically consistent under DLCoal. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗