45,055 research outputs found
Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method
This paper addresses optimization problems constrained by partial
differential equations with uncertain coefficients. In particular, the robust
control problem and the average control problem are considered for a tracking
type cost functional with an additional penalty on the variance of the state.
The expressions for the gradient and Hessian corresponding to either problem
contain expected value operators. Due to the large number of uncertainties
considered in our model, we suggest to evaluate these expectations using a
multilevel Monte Carlo (MLMC) method. Under mild assumptions, it is shown that
this results in the gradient and Hessian corresponding to the MLMC estimator of
the original cost functional. Furthermore, we show that the use of certain
correlated samples yields a reduction in the total number of samples required.
Two optimization methods are investigated: the nonlinear conjugate gradient
method and the Newton method. For both, a specific algorithm is provided that
dynamically decides which and how many samples should be taken in each
iteration. The cost of the optimization up to some specified tolerance
is shown to be proportional to the cost of a gradient evaluation with requested
root mean square error . The algorithms are tested on a model elliptic
diffusion problem with lognormal diffusion coefficient. An additional nonlinear
term is also considered.Comment: This work was presented at the IMG 2016 conference (Dec 5 - Dec 9,
2016), at the Copper Mountain conference (Mar 26 - Mar 30, 2017), and at the
FrontUQ conference (Sept 5 - Sept 8, 2017
Optimization Methods for Inverse Problems
Optimization plays an important role in solving many inverse problems.
Indeed, the task of inversion often either involves or is fully cast as a
solution of an optimization problem. In this light, the mere non-linear,
non-convex, and large-scale nature of many of these inversions gives rise to
some very challenging optimization problems. The inverse problem community has
long been developing various techniques for solving such optimization tasks.
However, other, seemingly disjoint communities, such as that of machine
learning, have developed, almost in parallel, interesting alternative methods
which might have stayed under the radar of the inverse problem community. In
this survey, we aim to change that. In doing so, we first discuss current
state-of-the-art optimization methods widely used in inverse problems. We then
survey recent related advances in addressing similar challenges in problems
faced by the machine learning community, and discuss their potential advantages
for solving inverse problems. By highlighting the similarities among the
optimization challenges faced by the inverse problem and the machine learning
communities, we hope that this survey can serve as a bridge in bringing
together these two communities and encourage cross fertilization of ideas.Comment: 13 page
A posteriori modeling error estimates in the optimization of two-scale elastic composite materials
The a posteriori analysis of the discretization error and the modeling error
is studied for a compliance cost functional in the context of the optimization
of composite elastic materials and a two-scale linearized elasticity model. A
mechanically simple, parametrized microscopic supporting structure is chosen
and the parameters describing the structure are determined minimizing the
compliance objective. An a posteriori error estimate is derived which includes
the modeling error caused by the replacement of a nested laminate
microstructure by this considerably simpler microstructure. Indeed, nested
laminates are known to realize the minimal compliance and provide a benchmark
for the quality of the microstructures. To estimate the local difference in the
compliance functional the dual weighted residual approach is used. Different
numerical experiments show that the resulting adaptive scheme leads to simple
parametrized microscopic supporting structures that can compete with the
optimal nested laminate construction. The derived a posteriori error indicators
allow to verify that the suggested simplified microstructures achieve the
optimal value of the compliance up to a few percent. Furthermore, it is shown
how discretization error and modeling error can be balanced by choosing an
optimal level of grid refinement. Our two scale results with a single scale
microstructure can provide guidance towards the design of a producible
macroscopic fine scale pattern
- …