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

Abner J. Salgado

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.OCmath.AP

Catalog footprint

What is connected

30works
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

30 published item(s)

preprint2022arXiv

Convergent, with rates, methods for normalized infinity Laplace, and related, equations

We propose a monotone, and consistent numerical scheme for the approximation of the Dirichlet problem for the normalized Infinity Laplacian, which could be related to the family of so--called two--scale methods. We show that this method is convergent, and prove rates of convergence. These rates depend not only on the regularity of the solution, but also on whether or not the right hand side vanishes. Some extensions to this approach, like obstacle problems and symmetric Finsler norms are also considered.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Diagonally implicit Runge-Kutta schemes: Discrete energy-balance laws and compactness properties

We study diagonally implicit Runge-Kutta (DIRK) schemes when applied to abstract evolution problems that fit into the Gelfand-triple framework. We introduce novel stability notions that are well-suited to this setting and provide simple, necessary and sufficient, conditions to verify that a DIRK scheme is stable in our sense and in Bochner-type norms. We use several popular DIRK schemes in order to illustrate cases that satisfy the required structural stability properties and cases that do not. In addition, under some mild structural conditions on the problem we can guarantee compactness of families of discrete solutions with respect to time discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Preconditioned accelerated gradient descent methods for locally Lipschitz smooth objectives with applications to the solution of nonlinear PDEs

We develop a theoretical foundation for the application of Nesterov's accelerated gradient descent method (AGD) to the approximation of solutions of a wide class of partial differential equations (PDEs). This is achieved by proving the existence of an invariant set and exponential convergence rates when its preconditioned version (PAGD) is applied to minimize locally Lipschitz smooth, strongly convex objective functionals. We introduce a second-order ordinary differential equation (ODE) with a preconditioner built-in and show that PAGD is an explicit time-discretization of this ODE, which requires a natural time step restriction for energy stability. At the continuous time level, we show an exponential convergence of the ODE solution to its steady state using a simple energy argument. At the discrete level, assuming the aforementioned step size restriction, the existence of an invariant set is proved and a matching exponential rate of convergence of the PAGD scheme is derived by mimicking the energy argument and the convergence at the continuous level. Applications of the PAGD method to numerical PDEs are demonstrated with certain nonlinear elliptic PDEs using pseudo-spectral methods for spatial discretization, and several numerical experiments are conducted. The results confirm the global geometric and mesh size-independent convergence of the PAGD method, with an accelerated rate that is improved over the preconditioned gradient descent (PGD) method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Time fractional gradient flows: Theory and numerics

We develop the theory of fractional gradient flows: an evolution aimed at the minimization of a convex, l.s.c.~energy, with memory effects. This memory is characterized by the fact that the negative of the (sub)gradient of the energy equals the so-called Caputo derivative of the state. We introduce the notion of energy solutions, for which we provide existence, uniqueness and certain regularizing effects. We also consider Lipschitz perturbations of this energy. For these problems we provide an a posteriori error estimate and show its reliability. This estimate depends only on the problem data, and imposes no constraints between consecutive time-steps. On the basis of this estimate we provide an a priori error analysis that makes no assumptions on the smoothness of the solution.

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 diffuse interface model for two-phase ferrofluid flows

We develop a model describing the behavior of two-phase ferrofluid flows using phase field-techniques and present an energy-stable numerical scheme for it. For a simplified, yet physically realistic, version of this model and the corresponding numerical scheme we prove, in addition to stability, convergence and as by-product existence of solutions. With a series of numerical experiments we illustrate the potential of these simple models and their ability to capture basic phenomenological features of ferrofluids such as the Rosensweig instability.

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

Numerical analysis of strongly nonlinear PDEs

We review the construction and analysis of numerical methods for strongly nonlinear PDEs, with an emphasis on convex and nonconvex fully nonlinear equations and the convergence to viscosity solutions. We begin by describing a fundamental result in this area which states that stable, consistent, and monotone schemes converge as the discretization parameter tends to zero. We review methodologies to construct finite difference, finite element, and semi-Lagrangian schemes that satisfy these criteria, and, in addition, discuss some rather novel tools that have paved the way to derive rates of convergence within this framework.

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
preprint2016arXiv

Stability analysis of pressure correction schemes for the Navier-Stokes equations with traction boundary conditions

We present a stability analysis for two different rotational pressure correction schemes with open and traction boundary conditions. First, we provide a stability analysis for a rotational version of the grad-div stabilized scheme of [A. Bonito, J.-L. Guermond, and S. Lee. Modified pressure-correction projection methods: Open boundary and variable time stepping. In Numerical Mathematics and Advanced Applications - ENUMATH 2013, volume 103 of Lecture Notes in Computational Science and Engineering, pages 623-631. Springer, 2015]. This scheme turns out to be unconditionally stable, provided the stabilization parameter is suitably chosen. We also establish a conditional stability result for the boundary correction scheme presented in [E. Bansch. A finite element pressure correction scheme for the Navier-Stokes equations with traction boundary condition. Comput. Methods Appl. Mech. Engrg., 279:198-211, 2014]. These results are shown by employing the equivalence between stabilized gauge Uzawa methods and rotational pressure correction schemes with traction boundary conditions.

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 Isaacs equation

We propose and analyze a two-scale finite element method for the Isaacs equation. The fine scale is given by the mesh size $h$ whereas the coarse scale $\varepsilon$ is dictated by an integro-differential approximation of the partial differential equation. We show that the method satisfies the discrete maximum principle provided that the mesh is weakly acute. This, in conjunction with weak operator consistency of the finite element method, allows us to establish convergence of the numerical solution to the viscosity solution as $\varepsilon, h\to0$, and $\varepsilon \gtrsim h^{1/2}|\log h|$. In addition, using a discrete Alexandrov Bakelman Pucci estimate we deduce rates of convergence, under suitable smoothness assumptions on the exact solution.

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
preprint2015arXiv

The equations of ferrohydrodynamics: modeling and numerical methods

We discuss the equations describing the motion of ferrofluids subject to an external magnetic field. We concentrate on the model proposed by R. Rosensweig, provide an appropriate definition for the effective magnetizing field, and explain the simplifications behind this definition. We show that this system is formally energy stable, and devise a numerical scheme that mimics the same stability estimate. We prove that solutions of the numerical scheme always exist and, under further simplifying assumptions, that the discrete solutions converge. We also discuss alternative formulations proposed in pre-existing work, primarily involving a regularization of the magnetization equation and supply boundary conditions which lead to an energy stable system. We present a series of numerical experiments which illustrate the potential of the scheme in the context of real applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A PDE approach to fractional diffusion: a posteriori error analysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence rates for the classical, thin and fractional elliptic obstacle problems

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.

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
preprint2013arXiv

Convergence analysis of a fractional time-stepping technique for incompressible fluids with microstructure

We present and analyze a fully discrete fractional time stepping technique for the solution of the micropolar Navier Stokes equations, which is a system of equations that describes the evolution of an incompressible fluid whose material particles possess both translational and rotational degrees of freedom. The proposed scheme uncouples the computation of the linear and angular velocity and the pressure. It is unconditionally stable and delivers optimal convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The micropolar Navier-Stokes equations: A priori error analysis

The unsteady Micropolar Navier-Stokes Equations (MNSE) are a system of parabolic partial differential equations coupling linear velocity and pressure with angular velocity: material particles have both translational and rotational degrees of freedom. We propose and analyze a first order semi-implicit fully-discrete scheme for the MNSE, which decouples the computation of the linear and angular velocities, is unconditionally stable and delivers optimal convergence rates under assumptions analogous to those used for the Navier-Stokes equations. With the help of our scheme we explore some qualitative properties of the MNSE related to ferrofluid manipulation and pumping. Finally, we propose a second order scheme and show that it is almost unconditionally stable.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Total Variation Diminishing Interpolation Operator and Applications

We construct, on continuous $Q_1$ finite elements over Cartesian meshes, an interpolation operator that does not increase the total variation. The operator is stable in $L^1$ and exhibits second order approximation properties. With the help of it we provide improved error estimates for discrete minimizers of the total variation denoising problem and for total variation flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discrete Total Variation Flows Without Regularization

We propose and analyze an algorithm for the solution of the $L^2$-subgradient flow of the total variation functional. The algorithm involves no regularization, thus the numerical solution preserves the main features that motivate practitioners to consider this type of energy. We propose an iterative scheme for the solution of the arising problems, show that the iterations converge, and develop a stopping criterion for them. We present numerical experiments which illustrate the power of the method, explore the solution behavior, and compare with regularized flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Diffuse Interface Model for Electrowetting with Moving Contact Lines

We introduce a diffuse interface model for the phenomenon of electrowetting on dielectric and present an analysis of the arising system of equations. Moreover, we study discretization techniques for the problem. The model takes into account different material parameters on each phase and incorporates the most important physical processes, such as incompressibility, electrostatics and dynamic contact lines; necessary to properly reflect the relevant phenomena. The arising nonlinear system couples the variable density incompressible Navier-Stokes equations for velocity and pressure with a Cahn-Hilliard type equation for the phase variable and chemical potential, a convection diffusion equation for the electric charges and a Poisson equation for the electric potential. Numerical experiments are presented, which illustrate the wide range of effects the model is able to capture, such as splitting and coalescence of droplets.

0 citations0 reviewsNo code linkedOpen access