This content will become publicly available on October 31, 2023
 Award ID(s):
 1835860
 Publication Date:
 NSFPAR ID:
 10377502
 Journal Name:
 Inverse Problems
 Volume:
 38
 Issue:
 12
 Page Range or eLocationID:
 125006
 ISSN:
 02665611
 Sponsoring Org:
 National Science Foundation
More Like this

The deformations of several slender structures at nanoscale are conceivably sensitive to their nonhomogenous elasticity. Owing to their small scale, it is not feasible to discern their elasticity parameter fields accurately using observations from physical experiments. Molecular dynamics simulations can provide an alternative or additional source of data. However, the challenges still lie in developing computationally efficient and robust methods to solve inverse problems to infer the elasticity parameter field from the deformations. In this paper, we formulate an inverse problem governed by a linear elastic model in a Bayesian inference framework. To make the problem tractable, we use a Gaussian approximation of the posterior probability distribution that results from the Bayesian solution of the inverse problem of inferring Young’s modulus parameter fields from available data. The performance of the computational framework is demonstrated using two representative loading scenarios, one involving cantilever bending and the other involving stretching of a helical rod (an intrinsically curved structure). The results show that smoothly varying parameter fields can be reconstructed satisfactorily from noisy data. We also quantify the uncertainty in the inferred parameters and discuss the effect of the quality of the data on the reconstructions.

For the inverse problem in physical models, one measures the solution and infers the model parameters using information from the collected data. Oftentimes, these data are inadequate and render the inverse problem illposed. We study the illposedness in the context of optical imaging, which is a medical imaging technique that uses light to probe (bio)tissue structure. Depending on the intensity of the light, the forward problem can be described by different types of equations. Highenergy light scatters very little, and one uses the radiative transfer equation (RTE) as the model; lowenergy light scatters frequently, so the diffusion equation (DE) suffices to be a good approximation. A multiscale approximation links the hyperbolictype RTE with the parabolictype DE. The inverse problems for the two equations have a multiscale passage as well, so one expects that as the energy of the photons diminishes, the inverse problem changes from well to illposed. We study this stability deterioration using the Bayesian inference. In particular, we use the Kullback–Leibler divergence between the prior distribution and the posterior distribution based on the RTE to prove that the information gain from the measurement vanishes as the energy of the photons decreases, so that the inverse problem is illposedmore »

Nonlinear statespace models are powerful tools to describe dynamical structures in complex time series. In a streaming setting where data are processed one sample at a time, simultaneous inference of the state and its nonlinear dynamics has posed significant challenges in practice. We develop a novel online learning framework, leveraging variational inference and sequential Monte Carlo, which enables flexible and accurate Bayesian joint filtering. Our method provides an approximation of the filtering posterior which can be made arbitrarily close to the true filtering distribution for a wide class of dynamics models and observation models. Specifically, the proposed framework can efficiently approximate a posterior over the dynamics using sparse Gaussian processes, allowing for an interpretable model of the latent dynamics. Constant time complexity per sample makes our approach amenable to online learning scenarios and suitable for realtime applications.

SUMMARY We introduce a new finiteelement (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an inplane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjointbased optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posteriormore »

We investigate solution methods for largescale inverse problems governed by partial differential equations (PDEs) via Bayesian inference. The Bayesian framework provides a statistical setting to infer uncertain parameters from noisy measurements. To quantify posterior uncertainty, we adopt Markov Chain Monte Carlo (MCMC) approaches for generating samples. To increase the efficiency of these approaches in highdimension, we make use of local information about gradient and Hessian of the target potential, also via Hamiltonian Monte Carlo (HMC). Our target application is inferring the field of soil permeability processing observations of pore pressure, using a nonlinear PDE poromechanics model for predicting pressure from permeability. We compare the performance of different sampling approaches in this and other settings. We also investigate the effect of dimensionality and nongaussianity of distributions on the performance of different sampling methods.