6,650 research outputs found
Sinc-Galerkin estimation of diffusivity in parabolic problems
A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise
A residual based snapshot location strategy for POD in distributed optimal control of linear parabolic equations
In this paper we study the approximation of a distributed optimal control
problem for linear para\-bolic PDEs with model order reduction based on Proper
Orthogonal Decomposition (POD-MOR). POD-MOR is a Galerkin approach where the
basis functions are obtained upon information contained in time snapshots of
the parabolic PDE related to given input data. In the present work we show that
for POD-MOR in optimal control of parabolic equations it is important to have
knowledge about the controlled system at the right time instances. For the
determination of the time instances (snapshot locations) we propose an
a-posteriori error control concept which is based on a reformulation of the
optimality system of the underlying optimal control problem as a second order
in time and fourth order in space elliptic system which is approximated by a
space-time finite element method. Finally, we present numerical tests to
illustrate our approach and to show the effectiveness of the method in
comparison to existing approaches
Improving Newton's method performance by parametrization: the case of Richards equation
The nonlinear systems obtained by discretizing degenerate parabolic equations
may be hard to solve, especially with Newton's method. In this paper, we apply
to Richards equation a strategy that consists in defining a new primary unknown
for the continuous equation in order to stabilize Newton's method by
parametrizing the graph linking the pressure and the saturation. The resulting
form of Richards equation is then discretized thanks to a monotone Finite
Volume scheme. We prove the well-posedness of the numerical scheme. Then we
show under appropriate non-degeneracy conditions on the parametrization that
Newton\^as method converges locally and quadratically. Finally, we provide
numerical evidences of the efficiency of our approach
A Certified Trust Region Reduced Basis Approach to PDE-Constrained Optimization
Parameter optimization problems constrained by partial differential equations (PDEs) appear in many science and engineering applications. Solving these optimization problems may require a prohibitively large number of computationally expensive PDE solves, especially if the dimension of the design space is large. It is therefore advantageous to replace expensive high-dimensional PDE solvers (e.g., finite element) with lower-dimensional surrogate models. In this paper, the reduced basis (RB) model reduction method is used in conjunction with a trust region optimization framework to accelerate PDE-constrained parameter optimization. Novel a posteriori error bounds on the RB cost and cost gradient for quadratic cost functionals (e.g., least squares) are presented and used to guarantee convergence to the optimum of the high-fidelity model. The proposed certified RB trust region approach uses high-fidelity solves to update the RB model only if the approximation is no longer sufficiently accurate, reducing the number of full-fidelity solves required. We consider problems governed by elliptic and parabolic PDEs and present numerical results for a thermal fin model problem in which we are able to reduce the number of full solves necessary for the optimization by up to 86%. Key words: model reduction, optimization, trust region methods, partial differential equations, reduced basis methods, error bounds, parametrized systemsFulbright U.S. Student ProgramNational Science Foundation (U.S.). Graduate Research Fellowship ProgramHertz FoundationUnited States. Department of Energy. Office of Advanced Scientific Computing Research (Award DEFG02-08ER2585)United States. Department of Energy. Office of Advanced Scientific Computing Research (Award DE-SC0009297
Optimal control of Allen-Cahn systems
Optimization problems governed by Allen-Cahn systems including elastic
effects are formulated and first-order necessary optimality conditions are
presented. Smooth as well as obstacle potentials are considered, where the
latter leads to an MPEC. Numerically, for smooth potential the problem is
solved efficiently by the Trust-Region-Newton-Steihaug-cg method. In case of an
obstacle potential first numerical results are presented
- …