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

Francisco-Javier Sayas

Francisco-Javier Sayas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysis

Catalog footprint

What is connected

29works
3topics
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

29 published item(s)

preprint2020arXiv

A note on devising HDG+ projections on polyhedral elements

In this note, we propose a simple way of constructing HDG+ projections on polyhedral elements. The projections enable us to analyze the Lehrenfeld-Schöberl HDG (HDG+) methods in a very concise manner, and make many existing analysis techniques of standard HDG methods reusable for HDG+. The novelty here is an alternative way of constructing the projections without using $M$-decompositions as a middle step. This extends our previous results [S. Du and F.-J. Sayas, SpringerBriefs in Mathematics (2019)] (elliptic problems) and [S. Du and F.-J. Sayas, Math.~Comp.~89~(2020), 1745-1782] (elasticity) to polyhedral meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unified error analysis of HDG methods for the static Maxwell equations

We propose a framework that allows us to analyze different variants of HDG methods for the static Maxwell equations using one simple analysis. It reduces all the work to the construction of projections that best fit the structures of the approximation spaces. As applications, we analyze four variants of HDG methods (denoted by B, H, B+, H+), where two of them are known (variants H, B+) and the other two are new (variants H+, B). Under certain regularity assumption, we show that all the four variants are optimally convergent and that variants B+ and H+ achieve superconvergence without post-processing. For the two known variants, we prove their optimal convergence under weaker requirements of the meshes and the stabilization functions thanks to the new analysis techniques being introduced. For solution with low-regularity, we give an analysis to these methods and investigate the effect of different stabilization functions on the convergence. At the end, we provide numerical experiments to support the analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An overlapping decomposition framework for wave propagation in heterogeneous and unbounded media: Formulation, analysis, algorithm, and simulation

A natural medium for wave propagation comprises a coupled bounded heterogeneous region and an unbounded homogeneous free-space. Frequency-domain wave propagation models in the medium, such as the variable coefficient Helmholtz equation, include a faraway decay radiation condition (RC). It is desirable to develop algorithms that incorporate the full physics of the heterogeneous and unbounded medium wave propagation model, and avoid an approximation of the RC. In this work we first present and analyze an overlapping decomposition framework that is equivalent to the full-space heterogeneous-homogenous continuous model, governed by the Helmholtz equation with a spatially dependent refractive index and the RC. Our novel overlapping framework allows the user to choose two free boundaries, and gain the advantage of applying established high-order finite and boundary element methods (FEM and BEM) to simulate an equivalent coupled model. The coupled model comprises auxiliary interior bounded heterogeneous and exterior unbounded homogeneous Helmholtz problems. A smooth boundary can be chosen for simulating the exterior problem using a spectrally accurate BEM, and a simple boundary can be used to setup a high-order FEM for the interior problem. Thanks to the spectral accuracy of the exterior computational model, the resulting coupled system in the overlapping region is relatively very small. Using the decomposed equivalent framework, we develop a novel overlapping FEM-BEM algorithm for simulating the acoustic or electromagnetic wave propagation in two dimensions. Our FEM-BEM algorithm for the full-space model incorporates the RC exactly. Numerical experiments demonstrate the efficiency of the FEM-BEM approach for simulating smooth and non-smooth wave fields, with the latter induced by a complex heterogeneous medium and a discontinuous refractive index.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A fully discrete BEM-FEM scheme for transient acoustic waves

We study a symmetric BEM-FEM coupling scheme for the scattering of transient acoustic waves by bounded inhomogeneous anisotropic obstacles in a homogeneous field. An incident wave in free space interacts with the obstacles and produces a combination of transmission and scattering. The transmitted part of the wave is discretized in space by finite elements while the scattered wave is reduced to two fields defined on the boundary of the obstacles and is discretized in space with boundary elements. We choose a coupling formulation that leads to a symmetric system of integro-differential equations. The retarded boundary integral equations are discretized in time by Convolution Quadrature, and the interior field is discretized in time with the trapezoidal rule. We show that the scattering problem generates a C_0 group of isometries in a Hilbert space, and use associated estimates to derive stability and convergence results. We provide numerical experiments and simulations to validate our results and demonstrate the flexibility of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

HDG methods for elastodynamics

We derive and analyze a hybridizable discontinuous Galerkin (HDG) method for approximating weak solutions to the equations of time-harmonic linear elasticity on a bounded Lipschitz domain in three dimensions. The real symmetry of the stress tensor is strongly enforced and its coefficients as well as those of the displacement vector field are approximated simultaneously at optimal convergence with respect to the choice of approximating spaces, wavenumber, and mesh size. Sufficient conditions are given so that the system is indeed transferable onto a global hybrid variable that, for larger polynomial degrees, may be approximated via a smaller-dimensional space than the original variables. We construct several variants of this method and discuss their advantages and disadvantages, and give a systematic approach to the error analysis for these methods. We touch briefly on the application of this error analysis to the time-dependent problem, and finally, we examine two different implementations of the method over various polynomial degrees and numerically demonstrate the convergence properties proven herein.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Boundary integral solvers for an evolutionary exterior Stokes problem

This paper proposes and analyzes a full discretization of the exterior transient Stokes problem with Dirichlet boundary conditions. The method is based on a single layer boundary integral representation, using Galerkin semidiscretization in the space variables and multistep Convolution Quadrature in time. Convergence estimates are based on a Laplace domain analysis, which translates into a detailed study of the exterior Brinkman problem. Some numerical experiments are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

New developments on the coupling of mixed-FEM and BEM for the three-dimensional Stokes problem

In this paper we consider the three dimensional exterior Stokes problem and study the solvability of the corresponding continuous and discrete formulations that arise from the coupling of a dual-mixed variational formulation with the boundary integral equation method. More precisely, after employing the incompressibility condition to eliminate the pressure, we consider the resulting velocity-stress-vorticity approach with different kind of boundary conditions on an annular bounded domain, and couple the underlying equations with either one or two boundary integral equations arising from the application of the usual and normal traces to the Green representation formula in the exterior unbounded region. As a result, we obtain saddle point operator equations, which are then analyzed by the well-known Babuska-Brezzi theory. We prove the well-posedness of the continuous formulations, identifying previously the space of solutions of the associated homogeneous problem, and specify explicit hypotheses to be satisfied by the finite element and boundary element subspaces in order to guarantee the stability of the respective Galerkin schemes. In particular, following a similar analysis given recently for the Laplacian, we are able to extend the classical Johnson & Nédélec procedure to the present case, without assuming any restrictive smoothness requirement on the coupling boundary, but only Lipschitz-continuity. In addition, and differently from known approaches for the elasticity problem, we are also able to extend the Costabel & Han coupling procedure to the 3D Stokes problem by providing a direct proof of the required coerciveness property, that is without argueing by contradiction, and by using the natural norm of each space instead of mesh-dependent norms. Finally, we briefly describe concrete examples of discrete spaces satisfying the aforementioned hypotheses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Costabel-Stephan system of Boundary Integral Equations in the Time Domain

In this paper we formulate a transmission problem for the transient acoustic wave equation as a system of retarded boundary integral equations. We then analyse a fully discrete method using a general Galerkin semidiscretization-in-space and Convolution Quadrature in time. All proofs are developed using recent techniques based on the theory of evolution equations. Some numerical experiments are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time-Dependent Fluid-Structure Interaction

The problem of determining the manner in which an incoming acoustic wave is scattered by an elastic body immersed in a fluid is one of central importance in detecting and identifying submerged objects. The problem is generally referred to as a fluid-structure interaction and is mathematically formulated as a time-dependent transmission problem. In this paper, we consider a typical fluid-structure interaction problem by using a coupling procedure which reduces the problem to a nonlocal initial-boundary problem in the elastic body with a system of integral equations on the interface between the domains occupied by the elastic body and the fluid. We analyze this nonlocal problem by the Lubich approach via the Laplace transform, an essential feature of which is that it works directly on data in the time domain rather than in the transformed domain. Our results may serve as a mathematical foundation for treating time-dependent fluid-structure interaction problems by convolution quadrature coupling of FEM and BEM.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fully discrete Kirchhoff formulas with CQ-BEM

In this paper we propose and analyze a fully discrete method for a direct boundary integral formulation of the scattering of a transient acoustic wave by a sound-soft obstacle. The method uses Galerkin-BEM in the space variables and three different choices of time-stepping strategies based on Convolution Quadrature. The numerical analysis of the method is carried out directly in the time domain, not reverting to Laplace transform techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Matlab tools for HDG in three dimensions

In this paper we provide some Matlab tools for efficient vectorized coding of the Hybridizable Discontinuous Galerkin for linear variable coefficient reaction-diffusion problems in polyhedral domains. The resulting tools are modular and include enhanced structures to deal with convection-diffusion problems, plus several projections and a superconvergent postprocess of the solution. Loops over the elements are exclusively local and, as such, have been parallelized.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetric coupling of LDG-FEM and DG-BEM

We analyze a discontinuous Galerkin FEM-BEM scheme for a second order elliptic transmission problem posed in the three-dimensional space. The symmetric variational formulation is discretized by nonconforming Raviart-Thomas finite elements on a general partition of the interior domain coupled with discontinuous boundary elements on an independent quasi-uniform mesh of the transmission interface. We prove (almost) quasi-optimal convergence of the method and confirm the theory by a numerical experiment. In addition, we consider the case when continuous rather than discontinuous boundary elements are used.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational views of stokeslets and stresslets

In this paper we present a self-contained variational theory of the layer potentials for the Stokes problem on Lipschitz boundaries. We use these weak definitions to show how to prove the main theorems about the associated Calderón projector. Finally, we relate these variational definitions to the integral forms. Instead of working these relations from scratch, we show some formulas parametrizing the Stokes layer potentials in terms of those for the Lamé and Laplace operators. While all the results in this paper are well known for smooth domains, and most might be known for non-smooth domains, the approach is novel a gives a solid structure to the theory of Stokes layer potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A decoupled preconditioning technique for a mixed Stokes-Darcy model

We propose an efficient iterative method to solve the mixed Stokes-Dracy model for coupling fluid and porous media flow. The weak formulation of this problem leads to a coupled, indefinite, ill-conditioned and symmetric linear system of equations. We apply a decoupled preconditioning technique requiring only good solvers for the local mixed-Darcy and Stokes subproblems. We prove that the method is asymptotically optimal and confirm, with numerical experiments, that the performance of the preconditioners does not deteriorate on arbitrarily fine meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A fully discrete Calderon Calculus for two dimensional time harmonic waves

In this paper, we present a fully discretized Calderón Calculus for the two dimensional Helmholtz equation. This full discretization can be understood as highly non-conforming Petrov-Galerkin methods, based on two staggered grids of mesh size $h$, Dirac delta distributions substituting acoustic charge densities and piecewise constant functions for approximating acoustic dipole densities. The resulting numerical schemes from this calculus are all of order $h^2$ provided that the continuous equations are well posed. We finish by presenting some numerical experiments illustrating the performance of this discrete calculus.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Nystrom method for the two dimensional Helmholtz hypersingular equation

In this paper we propose and analyze a class of simple Nyström discretizations of the hypersingular integral equation for the Helmholtz problem on domains of the plane with smooth parametrizable boundary. The method depends on a parameter (related to the staggering of two underlying grids) and we show that two choices of this parameter produce convergent methods of order two, while all other stable methods provide methods of order one. Convergence is shown for the density (in uniform norm) and for the potential postprocessing of the solution. Some numerical experiments are given to illustrate the performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence analysis of a high-order Nystrom integral-equation method for surface scattering problems

In this paper we present a convergence analysis for the Nystrom method proposed in [Jour. Comput. Phys. 169 pp. 2921-2934, 2001] for the solution of the combined boundary integral equation formulations of sound-soft acoustic scattering problems in three-dimensional space. This fast and efficient scheme combines FFT techniques and a polar change of variables that cancels out the kernel singularity. We establish the stability of the algorithms in the $L^2$ norm and we derive convergence estimates in both the $L^2$ and $L^\infty$ norms. In particular, our analysis establishes theoretically the previously observed super-algebraic convergence of the method in cases in which the right-hand side is smooth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coupling of HDG with a double-layer potential BEM

In this paper we propose and analyze a new coupling procedure for the Hybridizable Discontinuous Galerkin Method with Galerkin Boundary Element Methods based on a double layer potential representation of the exterior component of the solution of a transmission problem. We show a discrete uniform coercivity estimate for the non-symmetric bilinear form and prove optimal convergence estimates for all the variables, as well as superconvergence for some of the discrete fields. Some numerical experiments support the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong coupling of finite element methods for the Stokes-Darcy problem

The aim of this paper is to propose a systematic way to obtain convergent finite element schemes for the Darcy-Stokes flow problem by combining well-known mixed finite elements that are separately convergent for Darcy and Stokes problems. In the approach in which the Darcy problem is set in its natural $\mathbf{H}(\text{div})$ formulation and the Stokes problem is expressed in velocity-pressure form, the transmission condition ensuring global mass conservation becomes essential. As opposed to the strategy that handles weakly this transmission condition through a Lagrange multiplier, we impose here this restriction exactly in the space of global velocity field. Our analysis of the Galerkin discretization of the resulting problem reveals that, if the mixed finite element space used in the Darcy domain admits an $\mathbf{H}(\text{div})$-stable discrete lifting of the normal trace, then it can be combined with any stable Stokes mixed finite element of the same order to deliver a stable global method with quasi-optimal convergence rate. Finally, we present a series of numerical tests confirming our theoretical convergence estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of a non-symmetric coupling of Interior Penalty DG and BEM

We analyze a non-symmetric coupling of interior penalty discontinuous Galerkin and boundary element methods in two and three dimensions. Main results are discrete coercivity of the method, and thus unique solvability, and quasi-optimal convergence. The proof of coercivity is based on a localized variant of the variational technique from [F.-J. Sayas, The validity of Johnson-Nédeléc's BEM-FEM coupling on polygonal interfaces, {\em SIAM J. Numer. Anal.}, 47(5):3451--3463, 2009]. This localization gives rise to terms which are carefully analyzed in fractional order Sobolev spaces, and by using scaling arguments for rigid transformations. Numerical evidence of the proven convergence properties has been published previously.

0 citations0 reviewsNo code linkedOpen access