Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse analysis”

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 37 records · Page 2

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗

Adsorption of ethane, ethene, and ethyne on reconstructed Fe 3 O 4 (001)

The adsorption of the C2 hydrocarbons including ethane, ethene, and ethyne are studied on magnetite Fe 3 O 4 (001) by a combination of molecular beam dosing, temperature programmed desorption, and X-ray photoelectron spectroscopy. The ethane desorption profile has a single temperature invariant peak at 100 K while ethene and ethyne exhibit additional peaks at ~120 K and ~135 K. An inversion analysis is used to extract coverage dependent desorption energies as well as coverage-averaged prefactors for each molecule. Ethene and ethyne exhibit moderate coverage dependent desorption energies decreasing from ~80 to ~30 kJ/mol at saturation while ethane shows a relatively coverage invariant desorption energy ~28 kJ/mol. The desorption energies of the C2 hydrocarbons increase in order of increasing bond order. This is likely due to the enhanced interaction of the ethene and ethyne p system to the coordinatively unsaturated octahedral Fe sites (Fe oct ) on the ¬oxide surface. The saturation coverages of each C2 hydrocarbon at 82 K were also determined. Furthermore, these match well with the 2-dimensional spacing determined from their liquid densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adsorption and reaction of methanol on Fe 3 O 4 (001)

The interaction of methanol with iron oxide surfaces is of interest due to its potential in hydrogen storage, and from a fundamental perspective as a chemical probe of the reactivity. We present here a study examining the adsorption and reaction of methanol on magnetite Fe 3 O 4 (001) at cryogenic temperatures using a combination of temperature programmed desorption, x-ray photoelectron spectroscopy, and scanning tunneling microscopy. We report the methanol desorption profile from Fe 3 O 4 (001) is complex, exhibiting peaks at 145, 175, 238, and 273 K, corresponding to the desorption of intact methanol, as well as peaks at 350 and 500 K due to the reaction of methoxy intermediates. The saturation of a monolayer of methanol corresponds to ~5 molecules/unit cell (u.c.), which is slightly higher than the number of surface octahedral iron atoms of 4/u.c. We probe the kinetics and thermodynamics of the desorption of molecular methanol using inversion analysis. The deconvolution of the complex desorption profile into individual peaks allows for calculations of both the desorption energy and the prefactor of each feature. The initial 0.7 methanol/u.c. reacts to form a methoxy and hydroxy intermediates by 180 K, which remains on the surface above room temperature after intact methanol has desorbed. The methoxy species react via one of two channels, a recombination reaction with surface hydroxyls to form additional methanol at ~350 K, and a disproportionation reaction to form methanol and formaldehyde at ~500 K. Only 20% of the methoxy species undergo the disproportionation reaction, with most of them reacting via the 350 K pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Protorheology

We all instinctively poke, bounce, scoop, and observe materials to understand rheological properties quickly. Yet, these observations are rarely analyzed quantitatively. To address this, here we introduce the paradigm of protorheology: approximate quantitative inference from simple observations. Several case studies demonstrate how protorheology is an inclusive entry to rheology for a broad range of practitioners and strengthens the confidence and interpretation of accurate laboratory measurements. We survey a range of creative tests according to which rheological phenomenon is revealed. Some new working equations are derived, and all working equations are summarized for convenient reference and comparison across different methods. Finally, this establishes a framework to enable increased use of photos, videos, and quantitative inference and to support the increasing interest in digital image analysis, inverse methods, and high-throughput characterization being applied to rheological properties.

42 ENGINEERING↗

A Multi-Scale Computational Platform for Predictive Modeling of Corrosion in Al-Steel Joints (Final Report)

The research team proposed to develop innovative multi-scale models to predict corrosion and the resulting mechanical performances in aluminum-steel joints. The methods of joining considered are resistance spot welding, self-piercing riveting, and rivet-welding, all suitable for mass production applications. The multi-scale models integrate high throughput first-principle calculations based on density functional theory (DFT), high throughput calculation of phase diagrams (CALPHAD) modeling, and finite element method (FEM) simulations. These models are to be validated through laboratory experiments. Furthermore, the models are available as open source so as to enable scientists and engineers in the community to adapt and contribute to the development and application. The approaches rely on the research team’s extensive experience on the prediction of properties of individual phases at finite temperatures and variable compositions through DFT calculations, and our broad expertise on dissimilar material joining and their corrosion. The proposed computational framework enables high throughput computations for improved predictions of corrosion and the associated mechanical performance in dissimilar material joints, resulting in significant reduction in computational time needed by the current state-of-the-art methods. With the participation of researchers from three universities, an auto manufacturer, two manufacturing technology/equipment suppliers, and a software developer/vendor, the interdisciplinary research team applies the technical development on both phase-based modeling and laboratory experiments into the automobile body joining processes for validation and technology demonstration. The global cost of corrosion was estimated at about 3.4% of the global GDP in 2013. By using available corrosion control practices, it is estimated a saving between 15-35% of the cost of corrosion. In the U.S., more than $276 billion is spent repairing corrosion damage. Prediction of the corrosion and its impact on performance of the dissimilar material joints is critical for reducing the massive number of the current corrosion-based recalls for automobiles. Thus, the project goal is to develop models to enable predictive maintenance and end-of-life planning of multi-metal joints with risk of corrosion under different conditions such as exposure to high temperatures in summer and salt solutions in winter, quantified through its pH. An academia-industry consortium led by the University of Michigan and including Pennsylvania State University, University of Illinois Urbana-Champaign, University of Georgia, General Motors Company, Livermore Software Technology Corporation, and Optimal Process Technologies, LLC. created multi-scale models for prediction of corrosion in aluminum-steel joint structures such of them used in vehicle subassemblies – chassis and transmission systems. Starting from the first principle calculations, the team developed mathematical and data-driven models to predict the metallic components, which are formed during joining of two metals, for example aluminum and steel - a lightweight multilateral system which is currently used in more than 60% car bodies. These models were used for simulating chemical reactions that are happening when the joining metallic components are exposed to high temperatures and different pH values. The team was able to predict how the corrosion installs on the metallic components and how they lead to a sudden failure of components in cars. Newly developed machine learning algorithms combining Science, Technology, Engineering and Math disciplines, advanced finite element simulation and experimental validations have been integrated in a platform for prediction of the corrosion evolution and prediction the failure of joints under mechanical loadings and fatigue. Moreover, based on machine learning and inverse analysis, the team proposed solutions for designing new metallic alloys less susceptible to corrosion when joining multi-material assembles. An average of 4% error compared with experiments was achieved for the most common joints that are used in vehicle subassemblies.

36 MATERIALS SCIENCE↗

Three-Dimensional Imaging Lidar for Characterizing Particle Fields and Organisms in the Mesopelagic Zone

The ocean’s mesopelagic zone is largely uncharacterized despite its vital role in sustaining ocean ecosystems. The composition, cycling, and fate of particle fields in the mesopelagic lacks an integrative multi-scale understanding of organism migration patterns, distribution, and diversity. This problem is addressed by combining complementary technologies with overlapping size spectra, including profiler mounted optical scattering sensors, profiler, and ship mounted acoustic devices, and a custom Unobtrusive Multi-Static Lidar Imager (UMSLI). This unique sensor suite can observe distributions of particles including organisms over a six order of magnitude dynamic size range, from microns to meters. Overlapping size ranges between different methods allows for cross-validation. This work focuses on the lidar imaging measurements and optical backscattering and attenuation, covering a combined particle size range of 0.1 mm to several cm. Particles at the small end of this range are sized using an existing backscattering time series inversion method after Briggs et al. (2013). Larger particles are resolved with UMSLI over an expanding volume using three-dimensional photo-realistic laser serial imaging. UMSLI’s image rectifying ability over time allows for derivation of particle concentration, size, and spatial distribution. Technical details on the development and post-processing methods for the novel UMSLI system are provided. Image resolved particle size distributions (PSDs) revealed a size shift from smaller to larger particles (>0.5 mm) as indicated by flatter slopes from dawn (slope = 2.6) to dusk (slope = 3.0). PSD trends are supported by an optical backscatter and transmissometer time series inversion analysis. Size shifts in the particle field are largely attributed to aggregation effects. Images support evidence of temporal variation between dusk and dawn stations through statistical analysis of particle concentrations for particle sizes 0.50–5.41 mm. Spatial analysis of the particle field revealed a dominantly uniform distributed marine snow background. The importance and potential of integrated approaches to studying particle and organism dynamics in ocean environments are discussed.

54 ENVIRONMENTAL SCIENCES↗

Mars Upper Atmospheric Temperature and Atomic Oxygen Density Derived from the O i 130.4 nm Emission Observed by NASA’s MAVEN Mission

The abundance of atomic oxygen is a major unknown for the Mars upper atmosphere, and systematic measurements of this key species is a primary objective of NASA’s Mars Atmosphere and Volatile Evolution (MAVEN) mission. Moreover, the Martian dayside temperature structure and its variability has been a subject of considerable debate in the past. In this study, an inversion analysis of the periapsis limb scans of the O i 130.4 nm emission observed by MAVEN is performed to quantify the neutral temperature and atomic oxygen abundance in the Mars upper atmosphere. Our analysis yields unprecedented resolution of the spatial variation of Mars atomic oxygen, showing that its exobase density decreases with the solar zenith angle (SZA) from ∼0° to 60° by about an order of magnitude, from ∼10{sup 8} to ∼10{sup 7} cm{sup −3} near perihelion (2015 April) and from ∼10{sup 7} to ∼10{sup 6} cm{sup −3} near aphelion (2017 November). For a given SZA the exobase densities near perihelion are about an order of magnitude higher than those near aphelion. The [O]/[CO{sub 2}] mixing ratio at 130 km altitude can vary from ∼0.1% to 4.6%, which also decreases with increasing SZA. Moreover, it is shown that the Mars exobase temperature does not vary significantly with SZA from ∼0° to ∼60°, which is estimated to be ∼239.2 ± 27.6 K near perihelion and ∼162.4 ± 19.9 K near aphelion, with a difference of 76.8 K between these two seasonal extremes.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation and (Inverse Modified) Error Analysis for Implicitly Templated ODE-Nets

We focus on learning unknown dynamics from data using ODE-nets templated on implicit numerical initial value problem solvers. First, we perform inverse modified error analysis of the ODE-nets using unrolled implicit schemes for ease of interpretation. It is shown that training an ODE-net using an unrolled implicit scheme returns a close approximation of an inverse modified differential equation (IMDE). In addition, we establish a theoretical basis for hyperparameter selection when training such ODE-nets, whereas current strategies usually treat numerical integration of ODE-nets as a black box. We thus formulate an adaptive algorithm which monitors the level of error and adapts the number of (unrolled) implicit solution iterations during the training process, so that the error of the unrolled approximation is less than the current learning loss. This helps accelerate training while maintaining accuracy. Several numerical experiments are performed to demonstrate the advantages of the proposed algorithm compared to nonadaptive unrollings and validate the theoretical analysis. Here, we also note that this approach naturally allows for incorporating partially known physical terms in the equations, giving rise to what is termed “gray box” identification.

ODE-nets↗

DNA Sequence Analysis of an Inversion Hot Spot in Lobeliaceae Plastomes

The evolution of plastid genomes (plastomes) in land plants is typically conservative, with extensive structural rearrangements present in only a few groups. Early Southern blot analysis identified two Lobelia species that minimally required deletion of the plastid gene accD and five inversions to account for their plastome arrangement relative to the ancestral organization. Sixty alternative 5-step inversion scenarios could account for the observed arrangement, but only one scenario was consistent with the criterion of ‘common cause’ attributable to a putative rearrangement hot spot at the accD deletion-site. Plastome sequencing demonstrated that this previously hypothesized inversion order is historically accurate. Detailed reconstructions of the ancestral plastome organization before and after each inversion are presented herein. Stem-loop and disruption-rescue models were evaluated for each inversion. One inversion has an obvious stem-loop basis, but the other four inversions were primarily caused by serial insertion of foreign (extra-plastid) DNA bearing large open-reading frames that disrupted plastome organization at the accD deletion-site, and complete plastomes were rescued by seemingly arbitrary ligation or fortuitous recombination at the other inversion endpoint. Transposed copies of DNA segments from elsewhere in the plastome are frequently inserted at inversion junctions, and four junctions are consistent with the stem-loop ligation model.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

Inverse prediction of PuO2 processing conditions using Bayesian seemingly unrelated regression with functional data

Over the past decade, a variety of innovative methodologies have been developed to better characterize the relationships between processing conditions and the physical, morphological, and chemical features of special nuclear material (SNM). Different processing conditions generate SNM products with different features, which are known as “signatures” because they are indicative of the processing conditions used to produce the material. These signatures can potentially allow a forensic analyst to determine which processes were used to produce the SNM and make inferences about where the material originated. This article investigates a statistical technique for relating processing conditions to the morphological features of PuO 2 particles. We develop a Bayesian implementation of seemingly unrelated regression (SUR) to inverse-predict unknown PuO 2 processing conditions from known PuO 2 features. Model results from simulated data demonstrate the usefulness of the technique. Applied to empirical data from a bench-scale experiment specifically designed with inverse prediction in mind, our model successfully predicts nitric acid concentration, while results for Pu concentration and precipitation temperature were equivalent to a simple mean model. Our technique compliments other recent methodologies developed for forensic analysis of nuclear material and can be generalized across the field of chemometrics for application to other materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Revisiting the Temporal Leontief Inverse: New Insights on the Analysis of Regional Technological Economic Change

The current availability of longer series of input-output tables, as well as the release of global input-output databases, has fostered a growing literature analyzing changes in the economic structure and their drivers. In this paper, we take advantage of these time-series by proposing a methodology designed to trace the contribution of different drivers of the change in interindustrial relationships over time. Based on the Temporal Leontief Inverse (TLI), the Extended TLI (ETLI) decomposes the economy-wide effects of changes in direct interindustrial links between years, isolating the impact of different determinants of economic (environmental, energy, etc.) spillovers according to the interests of the researcher. For example, one can explore how the multipliers of a particular industry were affected by changes in technology of other sectors and in the own sector; by changes in trade patterns in specific countries; by indirect changes in intraregional production chains in foreign nations; etc. The ETLI is illustrated by uncovering certain hidden effects not captured in a previous application of the original TLI to the Chicago region between 1980-1997.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗