NAApr 4, 2017
Sparse optimal control for fractional diffusionEnrique Otárola, Abner J. Salgado
We consider an optimal control problem that entails the minimization of a nondifferentiable cost functional, fractional diffusion as state equation and constraints on the control variable. We provide existence, uniqueness and regularity results together with first order optimality conditions. In order to propose a solution technique, we realize fractional diffusion as the Dirichlet-to-Neumann map for a nonuniformly elliptic operator and consider an equivalent optimal control problem with a nonuniformly elliptic equation as state equation. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We propose a fully discrete scheme: piecewise constant functions for the control variable and first--degree tensor product finite elements for the state variable. We derive a priori error estimates for the control and state variables which are quasi--optimal with respect to degrees of freedom.
NAJan 11, 2019
Adaptive finite element methods for sparse PDE-constrained optimizationAlejandro Allendes, Francisco Fuica, Enrique Otárola
We propose and analyze reliable and efficient a posteriori error estimators for an optimal control problem that involves a nondifferentiable cost functional, the Poisson problem as state equation and control constraints. To approximate the solution to the state and adjoint equations we consider a piecewise linear finite element method whereas three different strategies are used to approximate the control variable: piecewise constant discretization, piecewise linear discretization and the so-called variational discretization approach. For the first two aforementioned solution techniques we devise an error estimator that can be decomposed as the sum of four contributions: two contributions that account for the discretization of the control variable and the associated subgradient, and two contributions related to the discretization of the state and adjoint equations. The error estimator for the variational discretization approach is decomposed only in two contributions that are related to the discretization of the state and adjoint equations. On the basis of the devised a posteriori error estimators, we design simple adaptive strategies that yield optimal rates of convergence for the numerical examples that we perform.
NAOct 4, 2018
Adaptive finite element methods for the pointwise tracking optimal control problem of the Stokes equationsAlejandro Allendes, Francisco Fuica, Enrique Otárola et al.
We propose and analyze a reliable and efficient a posteriori error estimator for the pointwise tracking optimal control problem of the Stokes equations. This linear-quadratic optimal control problem entails the minimization of a cost functional that involves point evaluations of the velocity field that solves the state equations. This leads to an adjoint problem with a linear combination of Dirac measures as a forcing term and whose solution exhibits reduced regularity properties. We also consider constraints on the control variable. The proposed a posteriori error estimator can be decomposed as the sum of four contributions: three contributions related to the discretization of the state and adjoint equations, and another contribution that accounts for the discretization of the control variable. On the basis of the devised a posteriori error estimator, we design a simple adaptive strategy that illustrates our theory and exhibits a competitive performance.
NANov 26, 2014
A PDE approach to fractional diffusion: a posteriori error analysisLong Chen, Ricardo H. Nochetto, Enrique Otárola et al.
We derive a computable a posteriori error estimator for the $α$-harmonic extension problem, which localizes the fractional powers of elliptic operators supplemented with Dirichlet boundary conditions. Our a posteriori error estimator relies on the solution of small discrete problems on anisotropic cylindrical stars. It exhibits built-in flux equilibration and is equivalent to the energy error up to data oscillation, under suitable assumptions. We design a simple adaptive algorithm and present numerical experiments which reveal a competitive performance.
NAOct 8, 2014
Convergence rates for the classical, thin and fractional elliptic obstacle problemsRicardo H. Nochetto, Enrique Otárola, Abner J. Salgado
We review the finite element approximation of the classical obstacle problem in energy and max-norms and derive error estimates for both the solution and the free boundary. On the basis of recent regularity results we present an optimal error analysis for the thin obstacle problem. Finally, we discuss the localization of the obstacle problem for the fractional Laplacian and prove quasi-optimal convergence rates.