Abstract:Purpose The purpose of this paper is to provide a framework for the implementation of an adjoint sensitivity formulation for least‐squares partial differential equations constrained optimization problems exploiting a multiphysics finite elements package. The estimation of the diffusion coefficient in a Poisson‐type diffusion equation is used as an example. Design/methodology/approach The authors derive the adjoint formulation in a continuous setting allowing to attribute to the direct and adjoint states the role of different fields to be solved for. They are one‐way coupled through the mismatch between measured and direct states acting as a source term in the adjoint equation. Having solved for the direct and adjoint state, the sensitivity of the cost function with respect to the design variables can then be obtained by a suitable post‐processing procedure. This sensitivity can then be used to efficiently solve the least‐squares problem. Findings The authors derived the adjoint formulation in a continuous setting allowing the direct and adjoint states to be attributed the role of different fields to be solved. They are one‐way coupled through the mismatch between measured and direct states acting as a source term in the adjoint equation. It is found that, having solved for the direct and adjoint state, the sensitivity of the cost function with respect to the design variables can then be obtained by a suitable post‐processing procedure. Research limitations/implications This paper implies that modern multiphysics finite elements packages provide a flexible and extendable software environment for the experimentation with different adjoint formulations. Such tools are therefore expected to become increasingly important in solving notoriously difficult partial differential equation (PDE)‐constrained least‐squares problems. The framework also provides the possibility of experimentation with different regularization techniques (total variation and multiscale techniques for instance) to handle the ill‐posedness of the problem. Originality/value In this paper the adjoint sensitivity computation is casted as a multiphysics problem allowing for a flexible and extendable implementation.

Enabling hyper-differential sensitivity analysis for ill-posed inverse problems

Hyper-Differential Sensitivity Analysis for Inverse Problems Constrained by Partial Differential Equations

Hyper-Differential Sensitivity Analysis of Uncertain Parameters in PDE-Constrained Optimization

Hyper-differential sensitivity analysis with respect to model discrepancy: Posterior Optimal Solution Sampling

Hyper-differential sensitivity analysis with respect to model discrepancy: Optimal solution updating

Hyper-differential sensitivity analysis with respect to model discrepancy: Mathematics and computation

Sensitivity Analysis of the Information Gain in Infinite-Dimensional Bayesian Linear Inverse Problems

Variance-based sensitivity of Bayesian inverse problems to the prior distribution

Surrogate-Based Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Structural Error.

Generalized Bayes Approach to Inverse Problems with Model Misspecification

Bayesian hybrid-kernel machine-learning-assisted sensitivity analysis and sensitivity-relevant inverse modeling for groundwater DNAPL contamination

Adjoint sensitivity in PDE constrained least squares problems as a multiphysics problem

A data-driven and model-based accelerated Hamiltonian Monte Carlo method for Bayesian elliptic inverse problems

Generalized sparsity-promoting solvers for Bayesian inverse problems: Versatile sparsifying transforms and unknown noise variances

Objective-sensitive principal component analysis for high-dimensional inverse problems

Derivative-Based Global Sensitivity Analysis for Models with High-Dimensional Inputs and Functional Outputs

A greedy sensor selection algorithm for hyperparameterized linear Bayesian inverse problems with correlated noise models

Goal-Oriented Sensitivity Analysis of Hyperparameters in Deep Learning

HJ-sampler: A Bayesian sampler for inverse problems of a stochastic process by leveraging Hamilton-Jacobi PDEs and score-based generative models

A-optimal encoding weights for nonlinear inverse problems, with application to the Helmholtz inverse problem

Adaptive Hyperbolic-cross-space Mapped Jacobi Method on Unbounded Domains with Applications to Solving Multidimensional Spatiotemporal Integrodifferential Equations