Source author record

Carmen Rodrigo

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

7 published item(s)

preprint2020arXiv

A two level method for isogeometric discretizations

Isogeometric Analysis (IGA) is a computational technique for the numerical approximation of partial differential equations (PDEs). This technique is based on the use of spline-type basis functions, that are able to hold a global smoothness and allow to exactly capture a wide set of common geometries. The current rise of this approach has encouraged the search of fast solvers for isogeometric discretizations and nowadays this topic is full of interest. In this framework, a desired property of the solvers is the robustness with respect to both the polinomial degree $p$ and the mesh size $h$. For this task, in this paper we propose a two-level method such that a discretization of order $p$ is considered in the first level whereas the second level consists of a linear or quadratic discretization. On the first level, we suggest to apply one single iteration of a multiplicative Schwarz method. The choice of the block-size of such an iteration depends on the spline degree $p$, and is supported by a local Fourier analysis (LFA). At the second level one is free to apply any given strategy to solve the problem exactly. However, it is also possible to get an approximation of the solution at this level by using an $h-$multigrid method. The resulting solver is efficient and robust with respect to the spline degree $p$. Finally, some numerical experiments are given in order to demonstrate the good performance of the proposed solver.

preprint2020arXiv

Robust preconditioners for a new stabilized discretization of the poroelastic equations

In this paper, we present block preconditioners for a stabilized discretization of the poroelastic equations developed in [45]. The discretization is proved to be well-posed with respect to the physical and discretization parameters, and thus provides a framework to develop preconditioners that are robust with respect to such parameters as well. We construct both norm-equivalent (diagonal) and field-of-value-equivalent (triangular) preconditioners for both the stabilized discretization and a perturbation of the stabilized discretization that leads to a smaller overall problem after static condensation. Numerical tests for both two- and three-dimensional problems confirm the robustness of the block preconditioners with respect to the physical and discretization parameters.

preprint2020arXiv

Using hierarchical matrices in the solution of the time-fractional heat equation by multigrid waveform relaxation

This work deals with the efficient numerical solution of the time-fractional heat equation discretized on non-uniform temporal meshes. Non-uniform grids are essential to capture the singularities of "typical" solutions of time-fractional problems. We propose an efficient space-time multigrid method based on the waveform relaxation technique, which accounts for the nonlocal character of the fractional differential operator. To maintain an optimal complexity, which can be obtained for the case of uniform grids, we approximate the coefficient matrix corresponding to the temporal discretization by its hierarchical matrix (${\cal H}$-matrix) representation. In particular, the proposed method has a computational cost of ${\cal O}(k N M \log(M))$, where $M$ is the number of time steps, $N$ is the number of spatial grid points, and $k$ is a parameter which controls the accuracy of the ${\cal H}$-matrix approximation. The efficiency and the good convergence of the algorithm, which can be theoretically justified by a semi-algebraic mode analysis, are demonstrated through numerical experiments in both one- and two-dimensional spaces.

preprint2016arXiv

A Nonconforming Finite Element Method for the Biot's Consolidation Model in Poroelasticity

A stable finite element scheme that avoids pressure oscillations for a three-field Biot's model in poroelasticity is considered. The involved variables are the displacements, fluid flux (Darcy velocity), and the pore pressure, and they are discretized by using the lowest possible approximation order: Crouzeix-Raviart finite elements for the displacements, lowest order Raviart-Thomas-Nedelec elements for the Darcy velocity, and piecewise constant approximation for the pressure. Mass lumping technique is introduced for the Raviart-Thomas-Nedelec elements in order to eliminate the Darcy velocity and, therefore, reduce the computational cost. We show convergence of the discrete scheme which is implicit in time and use these types of elements in space with and without mass lumping. Finally, numerical experiments illustrate the convergence of the method and show its effectiveness to avoid spurious pressure oscillations when mass lumping for the Raviart-Thomas-Nedelec elements is used.

preprint2015arXiv

A Finite Element Framework for Some Mimetic Finite Difference Discretizations

In this work we derive equivalence relations between mimetic finite difference schemes on simplicial grids and modified Nédélec-Raviart-Thomas finite element methods for model problems in $\mathbf{H}(\operatorname{\mathbf{curl}})$ and $H(\operatorname{div})$. This provides a simple and transparent way to analyze such mimetic finite difference discretizations using the well-known results from finite element theory. The finite element framework that we develop is also crucial for the design of efficient multigrid methods for mimetic finite difference discretizations, since it allows us to use canonical inter-grid transfer operators arising from the finite element framework. We provide special Local Fourier Analysis and numerical results to demonstrate the efficiency of such multigrid methods.

preprint2015arXiv

On a local Fourier analysis for overlapping block smoothers on triangular grids

A general local Fourier analysis for overlapping block smoothers on triangular grids is presented. This analysis is explained in a general form for its application to problems with different discretizations. This tool is demonstrated for two different problems: a stabilized linear finite element discretization of Stokes equations and an edge-based discretization of the curl-curl operator by lowest-order Nédélec finite element method. In this latter, special Fourier modes have to be considered in order to perform the analysis. Numerical results comparing two- and three-grid convergence factors predicted by the local Fourier analysis to real asymptotic convergence factors are presented to confirm the predictions of the analysis and show their usefulness.

preprint2014arXiv

Local Fourier Analysis of Multigrid Methods with Polynomial Smoothers and Aggressive coarsening

We focus on the study of multigrid methods with aggressive coarsening and polynomial smoothers for the solution of the linear systems corresponding to finite difference/element discretizations of the Laplace equation. Using local Fourier analysis we determine automatically the optimal values for the parameters involved in defining the polynomial smoothers and achieve fast convergence of cycles with aggressive coarsening. We also present numerical tests supporting the theoretical results and the heuristic ideas. The methods we introduce are highly parallelizable and efficient multigrid algorithms on structured and semi-structured grids in two and three spatial dimensions.