Search

A posteriori error estimator is derived for an elliptic interface problem in the fictitious domain formulation with distributed Lagrange multiplier considering a discontinuous Lagrange multiplier finite element space. A posteriori error estimation plays a pivotal role in assessing the accuracy and reliability of computational solutions across various domains of science and engineering. This study delves into the theoretical underpinnings and computational considerations of a residual-based estimator. Theoretically, the estimator is studied for cases with constant coefficients which jump across an interface as well as generalized scenarios with smooth coefficients that jump across an interface. Theoretical findings demonstrate the reliability and efficiency of the proposed estimators under all considered cases. Numerical experiments are conducted to validate the theoretical results, incorporating various immersed geometries and instances of high coefficients jumps at the interface. Leveraging an adaptive algorithm, the estimator identifies regions with singularities and applies refinement accordingly. Results substantiate the theoretical findings, highlighting the reliability and efficiency of the estimators. Furthermore, numerical solutions exhibit optimal convergence properties, demonstrating resilience against geometric singularities or coefficients jumps., Comment: 32 pages, 36 figures

We consider a fictitious domain formulation for fluid-structure interaction problems based on a distributed Lagrange multiplier to couple the fluid and solid behaviors. How to deal with the coupling term is crucial since the construction of the associated finite element matrix requires the integration of functions defined over non-matching grids: the exact computation can be performed by intersecting the involved meshes, whereas an approximate coupling matrix can be evaluated on the original meshes by introducing a quadrature error. The purpose of this paper is twofold: we prove that the discrete problem is well-posed also when the coupling term is constructed in approximate way and we discuss quadrature error estimates over non-matching grids., Comment: 27 pages, 6 figures

This paper focuses on the numerical solution of elliptic partial differential equations (PDEs) with Dirichlet and mixed boundary conditions, specifically addressing the challenges arising from irregular domains. Both finite element method (FEM) and finite difference method (FDM), face difficulties in dealing with arbitrary domains. The paper introduces a novel nodal symmetric ghost {method based on a variational formulation}, which combines the advantages of FEM and FDM. The method employs bilinear finite elements on a structured mesh and provides a detailed implementation description. A rigorous a priori convergence rate analysis is also presented. The convergence rates are validated with many numerical experiments, in both one and two space dimensions.

In this paper we introduce an abstract setting for the convergence analysis of the virtual element approximation of an acoustic vibration problem. We discuss the effect of the stabilization parameters and remark that in some cases it is possible to achieve optimal convergence without the need of any stabilization. This statement is rigorously proved for lowest order triangular element and supported by several numerical experiments., Comment: 36 pages, 7 figures, 14 tables

We discuss the approximation of the eigensolutions associated with the Maxwell eigenvalues problem in the framework of least-squares finite elements. We write the Maxwell curl curl equation as a system of two first order equation and design a novel least-squares formulation whose minimum is attained at the solution of the system. The eigensolution are then approximated by considering the eigenmodes of the underlying solution operator. We study the convergence of the finite element approximation and we show several numerical tests confirming the good behavior of the method. It turns out that nodal elements can be successfully employed for the approximation of our problem also in presence of singular solutions.

We perform a systematic comparison of various numerical schemes for the approximation of interface problems. We consider unfitted approaches in view of their application to possibly moving configurations. Particular attention is paid to the implementation aspects and to the analysis of the costs related to the different phases of the simulations., Comment: 25 pages, 15 figures, 15 tables

We consider nodal-based Lagrangian interpolations for the finite element approximation of the Maxwell eigenvalue problem. The first approach introduced is a standard Galerkin method on Powell-Sabin meshes, which has recently been shown to yield convergent approximations in two dimensions, whereas the other two are stabilized formulations that can be motivated by a variational multiscale approach. For the latter, a mixed formulation equivalent to the original problem is used, in which the operator has a saddle point structure. The Lagrange multiplier introduced to enforce the divergence constraint vanishes in an appropriate functional setting. The first stabilized method we consider consists of an augmented formulation with the introduction of a mesh dependent term that can be regarded as the Laplacian of the multiplier of the divergence constraint. The second formulation is based on orthogonal projections, which can be recast as a residual based stabilization technique. We rely on the classical spectral theory to analyze the approximating methods for the eigenproblem. The stability and convergence aspects are inherited from the associated source problems. We investigate the numerical performance of the proposed formulations and provide some convergence results validating the theoretical ones for several benchmark tests, including ones with smooth and singular solutions.

We use a Gaussian Process Regression (GPR) strategy that was recently developed [3,16,17] to analyze different types of curves that are commonly encountered in parametric eigenvalue problems. We employ an offline-online decomposition method. In the offline phase, we generate the basis of the reduced space by applying the proper orthogonal decomposition (POD) method on a collection of pre-computed, full-order snapshots at a chosen set of parameters. Then, we generate our GPR model using four different Mat\'{e}rn covariance functions. In the online phase, we use this model to predict both eigenvalues and eigenvectors at new parameters. We then illustrate how the choice of each covariance function influences the performance of GPR. Furthermore, we discuss the connection between Gaussian Process Regression and spline methods and compare the performance of the GPR method against linear and cubic spline methods. We show that GPR outperforms other methods for functions with a certain regularity.

In this article we apply reduced order techniques for the approximation of parametric eigenvalue problems. The effect of the choice of sampling points is investigated. Here we use the standard proper orthogonal decomposition technique to obtain the basis of the reduced space and Galerking orthogonal technique is used to get the reduced problem. We present some numerical results and observe that the use of sparse sampling is a good idea for sampling if the dimension of parameter space is high., Comment: 6 pages

A finite element discretization is developed for the Cai-Hu model, describing the formation of biological networks. The model consists of a non linear elliptic equation for the pressure $p$ and a non linear reaction-diffusion equation for the conductivity tensor $\mathbb{C}$. The problem requires high resolution due to the presence of multiple scales, the stiffness in all its components and the non linearities. We propose a low order finite element discretization in space coupled with a semi-implicit time advancing scheme. The code is {verified} with several numerical tests performed with various choices for the parameters involved in the system. In absence of the exact solution, we apply Richardson extrapolation technique to estimate the order of the method., Comment: 11 pages, 3 figures, 18 plots, 2 tables

In this paper we discuss reduced order models for the approximation of parametric eigenvalue problems. In particular, we are interested in the presence of intersections or clusters of eigenvalues. The singularities originating by these phenomena make it hard a straightforward generalization of well known strategies normally used for standards PDEs. We investigate how the known results extend (or not) to higher order frequencies.

We present and analyze a parallel solver for the solution of fluid structure interaction problems described by a fictitious domain approach. In particular, the fluid is modeled by the non-stationary incompressible Navier-Stokes equations, while the solid evolution is represented by the elasticity equations. The parallel implementation is based on the PETSc library and the solver has been tested in terms of robustness with respect to mesh refinement and weak scalability by running simulations on a Linux cluster., Comment: Contribution to the 5th African Conference on Computational Mechanics

In this article, we propose a data-driven reduced basis (RB) method for the approximation of parametric eigenvalue problems. The method is based on the offline and online paradigms. In the offline stage, we generate snapshots and construct the basis of the reduced space, using a POD approach. Gaussian process regressions (GPR) are used for approximating the eigenvalues and projection coefficients of the eigenvectors in the reduced space. All the GPR corresponding to the eigenvalues and projection coefficients are trained in the offline stage, using the data generated in the offline stage. The output corresponding to new parameters can be obtained in the online stage using the trained GPR. The proposed algorithm is used to solve affine and non-affine parameter-dependent eigenvalue problems. The numerical results demonstrate the robustness of the proposed non-intrusive method.

The aim of this work is to present a parallel solver for a formulation of fluid-structure interaction (FSI) problems which makes use of a distributed Lagrange multiplier in the spirit of the fictitious domain method. The fluid subproblem, consisting of the non-stationary Stokes equations, is discretized in space by $\mathcal{Q}_2$-$\mathcal{P}_1$ finite elements, whereas the structure subproblem, consisting of the linear or finite incompressible elasticity equations, is discretized in space by $\mathcal{Q}_1$ finite elements. A first order semi-implicit finite difference scheme is employed for time discretization. The resulting linear system at each time step is solved by a parallel GMRES solver, accelerated by block diagonal or triangular preconditioners. The parallel implementation is based on the PETSc library. Several numerical tests have been performed on Linux clusters to investigate the effectiveness of the proposed FSI solver., Comment: 29 pages, 8 figures, 14 tables

In this paper, new unfitted mixed finite elements are presented for elliptic interface problems with jump coefficients. Our model is based on a fictitious domain formulation with distributed Lagrange multiplier. The relevance of our investigations is better seen when applied to the framework of fluid structure interaction problems. Two finite elements schemes with piecewise constant Lagrange multiplier are proposed and their stability is proved theoretically. Numerical results compare the performance of those elements, confirming the theoretical proofs and verifying that the schemes converge with optimal rate., Comment: 28 pages, 24 figures. Numer Methods Partial Differential Eq. 2023

We compare the solutions of two systems of partial differential equations (PDE), seen as two different interpretations of the same model that describes formation of complex biological networks. Both approaches take into account the time evolution of the medium flowing through the network, and we compute the solution of an elliptic-parabolic PDE system for the conductivity vector $m$, the conductivity tensor $\mathbb{C}$ and the pressure $p$. We use finite differences schemes in a uniform Cartesian grid in the spatially two-dimensional setting to solve the two systems, where the parabolic equation is solved by a semi-implicit scheme in time. Since the conductivity vector and tensor appear also in the Poisson equation for the pressure $p$, the elliptic equation depends implicitly on time. For this reason we compute the solution of three linear systems in the case of the conductivity vector $m\in\mathbb{R}^2$, and four linear systems in the case of the symmetric conductivity tensor $\mathbb{C}\in\mathbb{R}^{2\times 2}$, at each time step. To accelerate the simulations, we make use of the Alternating Direction Implicit (ADI) method. The role of the parameters is important for obtaining detailed solutions. We provide numerous tests with various values of the parameters involved, to see the differences in the solutions of the two systems., Comment: 22 pages, 8 figures, 6 tables

In this paper we introduce an algorithm based on a sparse grid adaptive refinement, for the approximation of the eigensolutions to parametric problems arising from elliptic partial differential equations. In particular, we are interested in detecting the crossing of the hypersurfaces describing the eigenvalues as a function of the parameters. The a priori matching is followed by an a posteriori verification, driven by a suitably defined error indicator. At a given refinement level, a sparse grid approach is adopted for the construction of the grid of the next level, by using the marking given by the a posteriori indicator. Various numerical tests confirm the good performance of the scheme.

In this paper a novel numerical approximation of parametric eigenvalue problems is presented. We motivate our study with the analysis of a POD reduced order model for a simple one dimensional example. In particular, we introduce a new algorithm capable to track the matching of eigenvalues when the parameters vary.

We study a non standard mixed formulation of the Poisson problem, sometimes known as dual mixed formulation. For reasons related to the equilibration of the flux, we use finite elements that are conforming in H(div) for the approximation of the gradients, even if the formulation would allow for discontinuous finite elements. The scheme is not uniformly inf-sup stable, but we can show existence and uniqueness of the solution, as well as optimal error estimates for the gradient variable when suitable regularity assumptions are made. Several additional remarks complete the paper, shedding some light on the sources of instability for mixed formulations.

We study a recent formulation for fluid-structure interaction problems based on the use of a distributed Lagrange multiplier in the spirit of the fictitious domain approach. In this paper, we focus our attention on a crucial computational aspect regarding the interface matrix for the finite element discretization: it involves integration of functions supported on two different meshes. Several numerical tests show that accurate computation of the interface matrix has to be performed in order to ensure the optimal convergence of the method., Comment: 30 pages, 18 figures

We prove convergence of the Maxwell eigenvalue problem using quadratic or higher Lagrange finite elements on Worsey-Farin splits in three dimensions. To do this, we construct two Fortin-like operators to prove uniform convergence of the corresponding source problem. We present numerical experiments to illustrate the theoretical results.

In this paper we consider a reduced order method for the approximation of the eigensolutions of the Laplace problem with Dirichlet boundary condition. We use a time continuation technique that consists in the introduction of a fictitious time parameter. We use a POD approach and we present some theoretical results showing how to choose the optimal dimension of the POD basis. The results of our computations, related to the first eigenvalue, confirm the optimal behavior of our approximate solution.

Stokes flows are a type of fluid flow where convective forces are small in comparison with viscous forces, and momentum transport is entirely due to viscous diffusion. Besides being routinely used as benchmark test cases in numerical fluid dynamics, Stokes flows are relevant in several applications in science and engineering including porous media flow, biological flows, microfluidics, microrobotics, and hydrodynamic lubrication. The present study concerns the discretization of the equations of motion of Stokes flows in three dimensions utilizing the MINI mixed finite element, focusing on the superconvergence of the method which was investigated with numerical experiments using five purpose-made benchmark test cases with analytical solution. Despite the fact that the MINI element is only linearly convergent according to standard mixed finite element theory, a recent theoretical development proves that, for structured meshes in two dimensions, the pressure superconverges with order 1.5, as well as the linear part of the computed velocity with respect to the piecewise-linear nodal interpolation of the exact velocity. The numerical experiments documented herein suggest a more general validity of the superconvergence in pressure, possibly to unstructured tetrahedral meshes and even up to quadratic convergence which was observed with one test problem, thereby indicating that there is scope to further extend the available theoretical results on convergence.

We continue the investigation on the spectrum of operators arising from the discretization of partial differential equations. In this paper we consider a three field formulation recently introduced for the finite element least-squares approximation of linear elasticity. We discuss in particular the distribution of the discrete eigenvalues in the complex plane and how they approximate the positive real eigenvalues of the continuous problem. The dependence of the spectrum on the Lam\'e parameters is considered as well and its behavior when approaching the incompressible limit., Comment: 11 pages, 12 figures

In this paper we provide some more details on the numerical analysis and we present some enlightening numerical results related to the spectrum of a finite element least-squares approximation of the linear elasticity formulation introduced recently. We show that, although the formulation is robust in the incompressible limit for the source problem, its spectrum is strongly dependent on the Lam\'e parameters and on the underlying mesh., Comment: 21 pages, 14 figures

In this paper we describe a computational model for the simulation of fluid-structure interaction problems based on a fictitious domain approach. We summarize the results presented over the last years when our research evolved from the Finite Element Immersed Boundary Method (FE-IBM) to the actual Finite Element Distributed Lagrange Multiplier method (FE-DLM). We recall the well-posedness of our formulation at the continuous level in a simplified setting. We describe various time semi-discretizations that provide unconditionally stable schemes. Finally we report the stability analysis for the finite element space discretization where some improvements and generalizations of the previous results are obtained., Comment: This paper is dedicated to the memory of Prof. Claudio Baiocchi

We discuss the approximation of eigenvalue problems associated with elliptic partial differential equations using the virtual element method. After recalling the abstract theory, we present a model problem, describing in detail the features of the scheme, and highligting the effects of the stabilizing parameters. We conlcude the discussion with a survey of several application examples.

In this paper, the discontinuous Petrov--Galerkin approximation of the Laplace eigenvalue problem is discussed. We consider in particular the primal and ultra weak formulations of the problem and prove the convergence together with a priori error estimates. Moreover, we propose two possible error estimators and perform the corresponding a posteriori error analysis. The theoretical results are confirmed numerically and it is shown that the error estimators can be used to design an optimally convergent adaptive scheme.

In this paper we consider the linearized version of a system of partial differential equations arising from a fluid-structure interaction model. We prove the existence and the uniqueness of the solution under natural regularity assumptions.

The scaled boundary finite element method (SBFEM) is a relatively recent boundary element method that allows the approximation of solutions to PDEs without the need of a fundamental solution. A theoretical framework for the convergence analysis of SBFEM is proposed here. This is achieved by defining a space of semi-discrete functions and constructing an interpolation operator onto this space. We prove error estimates for this interpolation operator and show that optimal convergence to the solution can be obtained in SBFEM. These theoretical results are backed by a numerical example., Comment: 15 pages, 3 figures

We consider finite element approximations of the Maxwell eigenvalue problem in two dimensions. We prove, in certain settings, convergence of the discrete eigenvalues using Lagrange finite elements. In particular, we prove convergence in three scenarios: piecewise linear elements on Powell--Sabin triangulations, piecewise quadratic elements on Clough--Tocher triangulations, and piecewise quartics (and higher) elements on general shape-regular triangulations. We provide numerical experiments that support the theoretical results. The computations also show that, on general triangulations, the eigenvalue approximations are very sensitive to nearly singular vertices, i.e., vertices that fall on exactly two "almost" straight lines.

In this paper we study the approximation of eigenvalues arising from the mixed Hellinger--Reissner elasticity problem by using the simple finite element using partial relaxation of $C^0$ vertex continuity of stresses introduced recently by Jun Hu and Rui Ma. We prove that the method converge when a residual type error estimator is considered and that the estimator decays optimally with respect to the number of degrees of freedom.

We study the approximation of the spectrum of least-squares operators arising from linear elasticity. We consider a two-field (stress/displacement) and a three-field (stress/displacement/vorticity) formulation; other formulations might be analyzed with similar techniques. We prove a priori estimates and we confirm the theoretical results with simple two-dimensional numerical experiments., Comment: Minor changes after referee reports

In this paper we discuss spectral properties of operators associated with the least-squares finite element approximation of elliptic partial differential equations. The convergence of the discrete eigenvalues and eigenfunctions towards the corresponding continuous eigenmodes is studied and analyzed with the help of appropriate $L^2$ error estimates. A priori and a posteriori estimates are proved.

We discuss the solution of eigenvalue problems associated with partial differential equations that can be written in the generalized form $\m{A}x=\lambda\m{B}x$, where the matrices $\m{A}$ and/or $\m{B}$ may depend on a scalar parameter. Parameter dependent matrices occur frequently when stabilized formulations are used for the numerical approximation of partial differential equations. With the help of classical numerical examples we show that the presence of one (or both) parameters can produce unexpected results., Comment: v2 contains minor descriptive modifications with respect to v1

In this paper we review the hypercircle method of Prager and Synge. This theory inspired several studies and induced an active research in the area of a posteriori error analysis. In particular, we review the Braess--Sch\"oberl error estimator in the context of the Poisson problem. We discuss adaptive finite element schemes based on two variants of the estimator and we prove the convergence and optimality of the resulting algorithms.

We consider a recently introduced formulation for fluid-structure interaction problems which makes use of a distributed Lagrange multiplier in the spirit of the fictitious domain method. In this paper we focus on time integration methods of second order based on backward differentiation formulae and on the Crank-Nicolson method. We show the stability properties of the resulting method; numerical tests confirm the theoretical results.

This paper derives a posteriori error estimates for the mixed numerical approximation of the Laplace eigenvalue problem with homogeneous Dirichlet boundary conditions. In particular, the resulting error estimator constitutes an upper bound for the error and is shown to be local efficient. Therefore, we present a reconstruction in the standard $H^1_0$-conforming space for the primal variable of the mixed Laplace eigenvalue problem. This reconstruction is performed locally on a set of vertex patches.

Super-convergence of order 1.5 in pressure and velocity has been experimentally investigated for the two-dimensional Stokes problem discretised with the MINI mixed finite element. Even though the classic mixed finite element theory for the MINI element guarantees linear convergence for the total error, recent theoretical results indicate that super-convergence of order 1.5 in pressure and of the linear part of the computed velocity to the piecewise linear nodal interpolation of the exact velocity is in fact possible with structured, three-directional triangular meshes. The numerical experiments presented here suggest a more general validity of super-convergence of order 1.5, possibly to automatically generated and unstructured triangulations. In addition, the approximating properties of the complete computed velocity have been compared with the approximating properties of the piecewise-linear part of the computed velocity, finding that the former is generally closer to the exact velocity, whereas the latter conserves mass better.

In this paper we prove the optimal convergence of a standard adaptive scheme based on edge finite elements for the approximation of the solutions of the eigenvalue problem associated with Maxwell's equations. The proof uses the known equivalence of the problem of interest with a mixed eigenvalue problem.

We present a distributed Lagrange multiplier formulation of the Finite Element Immersed Boundary Method to couple incompressible fluids with compressible solids. This is a generalization of the formulation presented in Heltai and Costanzo (2012), that offers a cleaner variational formulation, thanks to the introduction of distributed Lagrange multipliers, that acts as intermediary between the fluid and solid equations, keeping the two formulation mostly separated. Stability estimates and a brief numerical validation are presented., Comment: Contribution to a conference in honor of Piero Colli Franzone