BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Enrique Otarola

Enrique Otarola appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.OCNumerical Analysismath.AP

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

A posteriori error estimates for an optimal control problem with a bilinear state equation

We propose and analyze a posteriori error estimators for an optimal control problem that involves an elliptic partial differential equation as state equation and a control variable that enters the state equation as a coefficient; pointwise constraints on the control variable are considered as well. We consider two different strategies to approximate optimal variables: a fully discrete scheme in which the admissible control set is discretized with piecewise constant functions and a semi-discrete scheme where the admissible control set is not discretized; the latter scheme being based on the so-called variational discretization approach. We design, for each solution technique, an a posteriori error estimator and show, in two and three dimensional Lipschitz polygonal/polyhedral domains (not necessarily convex), that the proposed error estimator is reliable and efficient. We design, based on the devised estimators, adaptive strategies that deliver optimal experimental rates of convergence for the performed numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Optimal Control Problem for the Navier--Stokes Equations with Point Sources

We analyze, in two dimensions, an optimal control problem for the Navier--Stokes equations where the control variable corresponds to the amplitude of forces modeled as point sources; control constraints are also considered. This particular setting leads to solutions to the state equation exhibiting reduced regularity properties. We operate under the framework of Muckenhoupt weights, Muckenhoupt-weighted Sobolev spaces, and the corresponding weighted norm inequalities and derive the existence of optimal solutions and first and, necessary and sufficient, second order optimality conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A posteriori error estimates for a distributed optimal control problem of the stationary Navier-Stokes equations

In two and three dimensional Lipschitz, but not necessarily convex, polytopal domains, we propose and analyze a posteriori error estimators for an optimal control problem involving the stationary Navier--Stokes equations; control constraints are also considered. We devise two strategies of discretization: a semi discrete scheme where the control variable is not discretized -- the so-called variational discrezation approach -- and a fully discrete scheme where the control is discretized with piecewise quadratic functions. For each solution solution technique, we design an a posteriori error estimator that can be decomposed as the sum of contributions related to the discretization of the state and adjoint equations and, additionally, the discretization of the control variable for when the fully discrete scheme is considered. We prove that the devised error estimators are reliable and also explore local efficiency estimates. Numerical experiments reveal a competitive performance of adaptive loops based on the devised a posteriori error estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Error estimates for optimal control problems involving the Stokes system and Dirac measures

The aim of this work is to derive a priori error estimates for finite element discretizations of control--constrained optimal control problems that involve the Stokes system and Dirac measures. The first 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. The second problem involves a control variable that corresponds to the amplitude of forces modeled as point sources. This leads to a solution of the state equations with reduced regularity properties. For each problem, we propose a finite element solution technique and derive a priori error estimates. Finally, we present numerical experiments, in two and three dimensions, that illustrate our theoretical developments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal control of a parabolic fractional PDE: analysis and discretization

We consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(Ω))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions and derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The stationary Boussinesq problem under singular forcing

In Lipschitz two and three dimensional domains, we study the existence for the so--called Boussinesq model of thermally driven convection under singular forcing. By singular we mean that the heat source is allowed to belong to $H^{-1}(\varpi,Ω)$, where $\varpi$ is a weight in the Muckenhoupt class $A_2$ that is regular near the boundary. We propose a finite element scheme and, under the assumption that the domain is convex and $\varpi^{-1} \in A_1$, show its convergence. In the case that the thermal diffusion and viscosity are constants, we propose an a posteriori error estimator and show its reliability and local efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A piecewise linear FEM for an optimal control problem of fractional operators: error analysis on curved domains

We propose and analyze a new discretization technique for a linear-quadratic optimal control problem involving the fractional powers of a symmetric and uniformly elliptic second oder operator; control constraints are considered. Since these fractional operators can be realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic equation, we recast our problem as a nonuniformly elliptic optimal control problem. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We propose a fully discrete scheme that is based on piecewise linear functions on quasi-uniform meshes to approximate the optimal control and first-degree tensor product functions on anisotropic meshes for the optimal state variable. We provide an a priori error analysis that relies on derived Holder and Sobolev regularity estimates for the optimal variables and error estimates for an scheme that approximates fractional diffusion on curved domains; the latter being an extension of previous available results. The analysis is valid in any dimension. We conclude by presenting some numerical experiments that validate the derived error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive finite element methods for an optimal control problem involving Dirac measures

The purpose of this work is the design and analysis of a reliable and efficient a posteriori error estimator for the so-called pointwise tracking optimal control problem. This linear-quadratic optimal control problem entails the minimization of a cost functional that involves point evaluations of the state, thus leading to an adjoint problem with Dirac measures on the right hand side; control constraints are also considered. The proposed error estimator relies on a posteriori error estimates in the maximum norm for the state and in Muckenhoupt weighted Sobolev spaces for the adjoint state. We present an analysis that is valid for two and three-dimensional domains. We conclude by presenting several numerical experiments which reveal the competitive performance of adaptive methods based on the devised error estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An a posteriori error analysis for an optimal control problem involving the fractional Laplacian

In a previous work, we introduced a discretization scheme for a constrained optimal control problem involving the fractional Laplacian. For such a control problem, we derived optimal a priori error estimates that demand the convexity of the domain and some compatibility conditions on the data. To relax such restrictions, in this paper, we introduce and analyze an efficient and, under certain assumptions, reliable a posteriori error estimator. We realize the fractional Laplacian as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi--infinite cylinder in one more spatial dimension. This extra dimension further motivates the design of an posteriori error indicator. The latter is defined as the sum of three contributions, which come from the discretization of the state and adjoint equations and the control variable. The indicator for the state and adjoint equations relies on an anisotropic error estimator in Muckenhoupt weighted Sobolev spaces. The analysis is valid in any dimension. On the basis of the devised a posteriori error estimator, we design a simple adaptive strategy that exhibits optimal experimental rates of convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimization with respect to order in a fractional diffusion model: analysis, approximation and algorithmic aspects

We consider an identification problem, where the state $u$ is governed by a fractional elliptic equation and the unknown variable corresponds to the order $s \in (0,1)$ of the underlying operator. We study the existence of an optimal pair $(\bar s, \bar u)$ and provide sufficient conditions for its local uniqueness. We develop semi-discrete and fully discrete algorithms to approximate the solutions to our identification problem and provide a convergence analysis. We present numerical illustrations that confirm and extend our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Some applications of weighted norm inequalities to the error analysis of PDE constrained optimization problems

The purpose of this work is to illustrate how the theory of Muckenhoupt weights, Muckenhoupt weighted Sobolev spaces and the corresponding weighted norm inequalities can be used in the analysis and discretization of PDE constrained optimization problems. We consider: a linear quadratic constrained optimization problem where the state solves a nonuniformly elliptic equation; a problem where the cost involves pointwise observations of the state and one where the state has singular sources, e.g. point masses. For all three examples we propose and analyze numerical schemes and provide error estimates in two and three dimensions. While some of these problems might have been considered before in the literature, our approach allows for a simpler, Hilbert space-based, analysis and discretization and further generalizations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A FEM for an optimal control problem of fractional powers of elliptic operators

We study solution techniques for a linear-quadratic optimal control problem involving fractional powers of elliptic operators. These fractional operators can be realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi-infinite cylinder in one more spatial dimension. Thus, we consider an equivalent formulation with a nonuniformly elliptic operator as state equation. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We discretize the proposed truncated state equation using first degree tensor product finite elements on anisotropic meshes. For the control problem we analyze two approaches: one that is semi-discrete based on the so-called variational approach, where the control is not discretized, and the other one is fully discrete via the discretization of the control by piecewise constant functions. For both approaches, we derive a priori error estimates with respect to the degrees of freedom. Numerical experiments validate the derived error estimates and reveal a competitive performance of anisotropic over quasi-uniform refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A fractional space-time optimal control problem: analysis and discretization\

We study a linear-quadratic optimal control problem involving a parabolic equation with fractional diffusion and Caputo fractional time derivative of orders $s \in (0,1)$ and $γ\in (0,1]$, respectively. The spatial fractional diffusion is realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic operator. Thus, we consider an equivalent formulation with a quasi-stationary elliptic problem with a dynamic boundary condition as state equation. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We consider a fully-discrete scheme: piecewise constant functions for the control and, for the state, first-degree tensor product finite elements in space and a finite difference discretization in time. We show convergence of this scheme and, for $s \in (0,1)$ and $γ= 1$, we derive a priori error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A PDE Approach to Numerical Fractional Diffusion

Fractional diffusion has become a fundamental tool for the modeling of multiscale and heterogeneous phenomena. However, due to its nonlocal nature, its accurate numerical approximation is delicate. We survey our research program on the design and analysis of efficient solution techniques for problems involving fractional powers of elliptic operators. Starting from a localization PDE result for these operators, we develop local techniques for their solution: a priori and a posteriori error analyses, adaptivity and multilevel methods. We show the flexibility of our approach by proposing and analyzing local solution techniques for a space-time fractional parabolic equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A PDE approach to space-time fractional parabolic problems

We study solution techniques for parabolic equations with fractional diffusion and Caputo fractional time derivative, the latter being discretized and analyzed in a general Hilbert space setting. The spatial fractional diffusion is realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi-infinite cylinder in one more spatial dimension. We write our evolution problem as a quasi-stationary elliptic problem with a dynamic boundary condition. We propose and analyze an implicit fully-discrete scheme: first-degree tensor product finite elements in space and an implicit finite difference discretization in time. We prove stability and error estimates for this scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite element approximation of the parabolic fractional obstacle problem

We study a discretization technique for the parabolic fractional obstacle problem in bounded domains. The fractional Laplacian is realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic equation posed on a semi-infinite cylinder, which recasts our problem as a quasi-stationary elliptic variational inequality with a dynamic boundary condition. The rapid decay of the solution suggests a truncation that is suitable for numerical approximation. We discretize the truncation with a backward Euler scheme in time and, for space, we use first-degree tensor product finite elements. We present an error analysis based on different smoothness assumptions

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multilevel methods for nonuniformly elliptic operators

We develop and analyze multilevel methods for nonuniformly elliptic operators whose ellipticity holds in a weighted Sobolev space with an $A_2$--Muckenhoupt weight. Using the so-called Xu-Zikatanov (XZ) identity, we derive a nearly uniform convergence result, under the assumption that the underlying mesh is quasi-uniform. We also consider the so-called $α$-harmonic extension to localize fractional powers of elliptic operators. Motivated by the scheme proposed in [R.H. Nochetto, E. Otarola, and A.J. Salgado. A PDE approach to fractional diffusion in general domains: a priori error analysis. arXiv:1302.0698, 2013] we present a multilevel method with line smoothers and obtain a nearly uniform convergence result on anisotropic meshes. Numerical experiments reveal a competitive performance of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Piecewise polynomial interpolation in Muckenhoupt weighted Sobolev spaces and applications

We develop a constructive piecewise polynomial approximation theory in weighted Sobolev spaces with Muckenhoupt weights for any polynomial degree. The main ingredients to derive optimal error estimates for an averaged Taylor polynomial are a suitable weighted Poincare inequality, a cancellation property and a simple induction argument. We also construct a quasi-interpolation operator, built on local averages over stars, which is well defined for functions in $L^1$. We derive optimal error estimates for any polynomial degree on simplicial shape regular meshes. On rectangular meshes, these estimates are valid under the condition that neighboring elements have comparable size, which yields optimal anisotropic error estimates over $n$-rectangular domains. The interpolation theory extends to cases when the error and function regularity require different weights. We conclude with three applications: nonuniform elliptic boundary value problems, elliptic problems with singular sources, and fractional powers of elliptic operators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A PDE approach to fractional diffusion in general domains: a priori error analysis

The purpose of this work is the study of solution techniques for problems involving fractional powers of symmetric coercive elliptic operators in a bounded domain with Dirichlet boundary conditions. These operators can be realized as the Dirichlet to Neumann map for a degenerate/singular elliptic problem posed on a semi-infinite cylinder, which we analyze in the framework of weighted Sobolev spaces. Motivated by the rapid decay of the solution of this problem, we propose a truncation that is suitable for numerical approximation. We discretize this truncation using first degree tensor product finite elements. We derive a priori error estimates in weighted Sobolev spaces. The estimates exhibit optimal regularity but suboptimal order for quasi-uniform meshes. For anisotropic meshes, instead, they are quasi-optimal in both order and regularity. We present numerical experiments to illustrate the method's performance.

0 citations0 reviewsNo code linkedOpen access