10,539 research outputs found
Adaptive Reconstruction for Electrical Impedance Tomography with a Piecewise Constant Conductivity
In this work we propose and analyze a numerical method for electrical
impedance tomography of recovering a piecewise constant conductivity from
boundary voltage measurements. It is based on standard Tikhonov regularization
with a Modica-Mortola penalty functional and adaptive mesh refinement using
suitable a posteriori error estimators of residual type that involve the state,
adjoint and variational inequality in the necessary optimality condition and a
separate marking strategy. We prove the convergence of the adaptive algorithm
in the following sense: the sequence of discrete solutions contains a
subsequence convergent to a solution of the continuous necessary optimality
system. Several numerical examples are presented to illustrate the convergence
behavior of the algorithm.Comment: 26 pages, 12 figure
Gradient-Based Estimation of Uncertain Parameters for Elliptic Partial Differential Equations
This paper addresses the estimation of uncertain distributed diffusion
coefficients in elliptic systems based on noisy measurements of the model
output. We formulate the parameter identification problem as an infinite
dimensional constrained optimization problem for which we establish existence
of minimizers as well as first order necessary conditions. A spectral
approximation of the uncertain observations allows us to estimate the infinite
dimensional problem by a smooth, albeit high dimensional, deterministic
optimization problem, the so-called finite noise problem in the space of
functions with bounded mixed derivatives. We prove convergence of finite noise
minimizers to the appropriate infinite dimensional ones, and devise a
stochastic augmented Lagrangian method for locating these numerically. Lastly,
we illustrate our method with three numerical examples
An optimal adaptive wavelet method for First Order System Least Squares
In this paper, it is shown that any well-posed 2nd order PDE can be
reformulated as a well-posed first order least squares system. This system will
be solved by an adaptive wavelet solver in optimal computational complexity.
The applications that are considered are second order elliptic PDEs with
general inhomogeneous boundary conditions, and the stationary Navier-Stokes
equations.Comment: 40 page
Progressive construction of a parametric reduced-order model for PDE-constrained optimization
An adaptive approach to using reduced-order models as surrogates in
PDE-constrained optimization is introduced that breaks the traditional
offline-online framework of model order reduction. A sequence of optimization
problems constrained by a given Reduced-Order Model (ROM) is defined with the
goal of converging to the solution of a given PDE-constrained optimization
problem. For each reduced optimization problem, the constraining ROM is trained
from sampling the High-Dimensional Model (HDM) at the solution of some of the
previous problems in the sequence. The reduced optimization problems are
equipped with a nonlinear trust-region based on a residual error indicator to
keep the optimization trajectory in a region of the parameter space where the
ROM is accurate. A technique for incorporating sensitivities into a
Reduced-Order Basis (ROB) is also presented, along with a methodology for
computing sensitivities of the reduced-order model that minimizes the distance
to the corresponding HDM sensitivity, in a suitable norm. The proposed reduced
optimization framework is applied to subsonic aerodynamic shape optimization
and shown to reduce the number of queries to the HDM by a factor of 4-5,
compared to the optimization problem solved using only the HDM, with errors in
the optimal solution far less than 0.1%
- …