112 research outputs found

    A Hybrid High-Order Method for a Class of Strongly Nonlinear Elliptic Boundary Value Problems

    Full text link
    In this article, we design and analyze a Hybrid High-Order (HHO) finite element approximation for a class of strongly nonlinear boundary value problems. We consider an HHO discretization for a suitable linearized problem and show its well-posedness using the Gardings type inequality. The essential ingredients for the HHO approximation involve local reconstruction and high-order stabilization. We establish the existence of a unique solution for the HHO approximation using the Brouwer fixed point theorem and contraction principle. We derive an optimal order a priori error estimate in the discrete energy norm. Numerical experiments are performed to illustrate the convergence histories.Comment: arXiv admin note: substantial text overlap with arXiv:2110.1557

    Superconvergent P1 honeycomb virtual elements and lifted P3 solutions

    Full text link
    When solving the Poisson equation on honeycomb hexagonal grids, we show that the P1P_1 virtual element is three-order superconvergent in H1H^1-norm, and two-order superconvergent in L2L^2 and L∞L^\infty norms. We define a local post-process which lifts the superconvergent P1P_1 solution to a P3P_3 solution of the optimal-order approximation. The theory is confirmed by a numerical test

    Mixed finite element approximation of porous media flows

    Get PDF
    The reliable simulation of flow in fractured porous media is a key aspect in the decision making process of stakeholders within politics and the geosciences, for example when assessing the suitability of burial sites for storage of high–level radioactive waste. This thesis aims to tackle the challenge that is the accurate simulation of these flows and does so via three computational developments. That is, suitable models for porous media flow with fractures; obtaining rigorous and reliable estimates of errors generated through these models; and the accurate simulation of the times–of–flight for particles transported by groundwater within the porous medium. Firstly, an expansion procedure for fractures in porous media is developed so that physical fluid laws are still retained when tracking particles across fracture–bulk interfaces. Moreover, the second contribution of this work is the utilisation of the dual–weighted–residual method to define suitable elementwise indicators for generic quantities of interest. The third contribution of this thesis is the attainment of accurate simulations of travel times for particles in porous media, achieved through linearising the functional representing the time–of–flight; in practice, numerical examples, including one inspired by the Sellafield site in Cumbria, UK, validate the performance of the proposed error estimator, and hence are useful in the safety assessment of storage facilities intended for radioactive waste

    A family of stabilizer-free virtual elements on triangular meshes

    Full text link
    A family of stabilizer-free PkP_k virtual elements are constructed on triangular meshes. When choosing an accurate and proper interpolation, the stabilizer of the virtual elements can be dropped while the quasi-optimality is kept. The interpolating space here is the space of continuous PkP_k polynomials on the Hsieh-Clough-Tocher macro-triangle, where the macro-triangle is defined by connecting three vertices of a triangle with its barycenter. We show that such an interpolation preserves PkP_k polynomials locally and enforces the coerciveness of the resulting bilinear form. Consequently the stabilizer-free virtual element solutions converge at the optimal order. Numerical tests are provided to confirm the theory and to be compared with existing virtual elements

    Scalable Recovery-based Adaptation on Quadtree Meshes for Advection-Diffusion-Reaction Problems

    Full text link
    We propose a mesh adaptation procedure for Cartesian quadtree meshes, to discretize scalar advection-diffusion-reaction problems. The adaptation process is driven by a recovery-based a posteriori estimator for the L2(Ω)L^2(\Omega)-norm of the discretization error, based on suitable higher order approximations of both the solution and the associated gradient. In particular, a metric-based approach exploits the information furnished by the estimator to iteratively predict the new adapted mesh. The new mesh adaptation algorithm is successfully assessed on different configurations, and turns out to perform well also when dealing with discontinuities in the data as well as in the presence of internal layers not aligned with the Cartesian directions. A cross-comparison with a standard estimate--mark--refine approach and with other adaptive strategies available in the literature shows the remarkable accuracy and parallel scalability of the proposed approach

    A posteriori error estimation and modeling of unsaturated flow in fractured porous media

    Get PDF
    This doctoral thesis focuses on three topics: (1) modeling of unsaturated flow in fractured porous media, (2) a posteriori error estimation for mixed-dimensional elliptic equations, and (3) contributions to open-source software for complex multiphysics processes in porous media. In our first contribution, following a Discrete-Fracture Matrix (DFM) approach, we propose a model where Richards' equation governs the water flow in the matrix, whereas fractures are represented as lower-dimensional open channels, naturally providing a capillary barrier to the water flow. Therefore, water in the matrix is only allowed to imbibe the fracture if the capillary barrier is overcome. When this occurs, we assume that the water inside the fracture flows downwards without resistance and, therefore, is instantaneously at hydrostatic equilibrium. This assumption can be justifiable for fractures with sufficiently large apertures, where capillary forces play no role. Mathematically, our model can be classified as a coupled PDE-ODE system of equations with variational inequalities, in which each fracture is considered a potential seepage face. Our second contribution deals with error estimation for mixed-dimensional (mD) elliptic equations, which, in particular, model single-phase flow in fractured porous media. Here, based on the theory of functional a posteriori error estimates, we derive guaranteed upper bounds for the mD primal and mD dual variables, and two-sided bounds for the mD primal-dual pair. Moreover, we improve the standard results of the functional approach by proposing four ways of estimating the residual errors based on the conservation properties of the approximations, that is, (1) no conservation, (2) subdomain conservation, (3) local conservation, and (4) pointwise conservation. This results in sharper and fully-computable bounds when mass is conserved either locally or exactly. To our knowledge, to date, no error estimates have been available for fracture networks, including fracture intersections and floating subdomains. Our last contribution is related to the development of open-source software. First, we present the implementation of a new multipoint finite-volume-based module for unsaturated poroelasticity, compatible with the Matlab Reservoir Simulation Toolbox (MRST). Second, we present a new Python-based simulation framework for multiphysics processes in fractured porous media, named PorePy. PorePy, by design, is particularly well-suited for handling mixed-dimensional geometries, and thus optimal for DFM models. The first two contributions discussed above were implemented in PorePy.Denne avhandlingen tar for seg tre emner: (1) modellering av flyt i umettet porøst medium med sprekker, (2) a posteriori feilestimater for blandet-dimensjonale elliptiske ligninger, og (3) bidrag til åpen kildekode for komplekse multifysikk-prosesser i porøse medier. I det første bidraget anvender vi en Discrete-Fracture Matrix (DFM) (Diskret-Sprekk Matrise) metode til å sette opp en modell hvor Richard's ligning modellerer vann-flyt i matrisen, og sprekkene representeres som lavere-dimensjonale åpne kanaler, som naturlig virker som kapillærbarrierer til vann-flyten. Derfor vil vann i matrisen kun få tilgang til sprekken når kapillærbarrieren blir brutt. Når det inntreffer, antar vi at vannet i sprekken flyter nedover uten motstand, og at hydrostatisk ekvilibrium derfor inntreffer øyeblikkelig. Slike antakelser kan rettferdiggjøres for sprekker med tilstrekkelig stor apertur (åpning), hvor kapillærkrefter ikke har noen innvirkning. Fra et matematisk standpunkt kan modellen klassifiseres som en sammenkoblet PDE-ODE med variasjonelle ulikheter hvor hver sprekk behandles som en filtreringsfase. Det andre bidraget tar for seg feilestimater for blandet-dimensjonale elliptiske ligninger, som modellerer en-fase flyt i porøse medier med sprekker. Her anvender vi teorien for "funksjonal a posteriori feilestimater" til å finne øvre skranker for primær og dual variablene, samt øvre og nedre skranker for primær-dual paret. Dessuten viser vi at vi kan forbedre standardresultatene fra "funksjonal a posteriori feilestimater" ved å foreslå fire måte å estimere residualfeilen basert på bevaringsegenskapene til diskretiseringen. De fire forskjellige bevaringsegenskapene er; ingen bevaringsegenskap, under- domene bevaring, lokal bevaring og punktvis bevaring. Dette fører til skarpere skranker som er mulige å beregne når masse er bevart enten lokalt, eller eksakt. Vi kjenner ikke til andre tilgjengelige feilestimater for sprekknettverk som inkluderer snitt av sprekker og sprekkrender som ligger innenfor domenets rand. Det siste bidraget omhandler utvikling av åpen kildekode. Først presenterer vi imple- menteringen av en multipunktfluks-basert modul for flyt i umettet deformerbart porøst medium som er kompatibelt med "Matlab Reservoir Simulation Toolbox"(MRST). I tillegg presenterer vi et nytt Python-basert rammeverk for simulering av multifysikkprosesser i porøse medier med sprekker, som heter PorePy. Dette rammeverket er designet for å håndtere geometrier med blandede dimensjoner og er derfor optimalt for DFM modeller. De to første bidragene i avhandlingen (nevnt over) er implementert i PorePy.Doktorgradsavhandlin

    A C0 Finite Element Method For The Biharmonic Problem In A Polygonal Domain

    Get PDF
    This dissertation studies the biharmonic equation with Dirichlet boundary conditions in a polygonal domain. The biharmonic problem appears in various real-world applications, for example in plate problems, human face recognition, radar imaging, and hydrodynamics problems. There are three classical approaches to discretizing the biharmonic equation in the literature: conforming finite element methods, nonconforming finite element methods, and mixed finite element methods. We propose a mixed finite element method that effectively decouples the fourth-order problem into a system of one steady-state Stokes equation and one Poisson equation. As a generalization to the above-decoupled formulation, we propose another decoupled formulation using a system of two Poison equations and one steady-state Stokes equation. We solve Poisson equations using linear and quadratic Lagrange\u27s elements and the Stokes equation using Hood-Taylor elements and Mini finite elements. It is shown that the solution of each system is equivalent to that of the original fourth-order problem on both convex and non-convex polygonal domains. Two finite element algorithms are, in turn, proposed to solve the decoupled systems. Solving this problem in a non-convex domain is challenging due to the singularity occurring near re-entrant corners. We introduce a weighted Sobolev space and a graded mesh refine Algorithm to attack the singularity near re-entrant corners. We show the regularity results of each decoupled system in both Sobolev space and weighted Sobolev space. We derive the H1H^1 and L2L^2 error estimates for the numerical solutions on quasi-uniform and graded meshes. We present various numerical test results to justify the theoretical findings. Given the availability of fast Poisson solvers and Stokes solvers, our Algorithm is a relatively easy and cost-effective alternative to existing algorithms for solving the biharmonic equation

    Unified discontinuous Galerkin scheme for a large class of elliptic equations

    Get PDF
    We present a discontinuous Galerkin internal-penalty scheme that is applicable to a large class of linear and nonlinear elliptic partial differential equations. The unified scheme can accommodate all second-order elliptic equations that can be formulated in first-order flux form, encompassing problems in linear elasticity, general relativity, and hydrodynamics, including problems formulated on a curved manifold. It allows for a wide range of linear and nonlinear boundary conditions, and accommodates curved and nonconforming meshes. Our generalized internal-penalty numerical flux and our Schur-complement strategy of eliminating auxiliary degrees of freedom make the scheme compact without requiring equation-specific modifications. We demonstrate the accuracy of the scheme for a suite of numerical test problems. The scheme is implemented in the open-source SpECTRE numerical relativity code

    Estimation d’erreur a posteriori pour l’approximation de problèmes Laplaciens fractionnaires et applications en poro-élasticité

    Get PDF
    This manuscript is concerned with a posteriori error estimation for the finite element discretization of standard and fractional partial differential equations as well as an application of fractional calculus to the modeling of the human meniscus by poro-elasticity equations. In the introduction, we give an overview of the literature of a posteriori error estimation in finite element methods and of adaptive refine- ment methods. We emphasize the state–of–the–art of the Bank–Weiser a posteriori error estimation method and of the adaptive refinement methods convergence results. Then, we move to fractional partial differential equations. We give some of the most common discretization methods of fractional Laplacian operator based equations. We review some results of a priori error estimation for the finite element discretization of these equations and give the state–of–the–art of a posteriori error estimation. Finally, we review the literature on the use of the Caputo’s fractional derivative in applications, focusing on anomalous diffusion and poro-elasticity applications. The rest of the manuscript is organized as follow. Chapter 1 is concerned with a proof of the reliability of the Bank–Weiser estimator for three–dimensional problems, extending a result from the literature. In Chapter 2 we present a numerical study of the Bank–Weiser estimator, provide a novel implementation of the estimator in the FEniCS finite element software and apply it to a variety of elliptic equations as well as goal-oriented error estimation. In Chapter 3 we derive a novel a posteriori estimator for the L2 error induced by the finite element discretization of fractional Laplacian operator based equations. In Chapter 4 we present new theoretical results on the convergence of a rational approximation method with consequences on the approximation of fractional norms as well as a priori error estimation results for the finite element discretization of fractional equations. Finally, in Chapter 5 we provide an application of fractional calculus to the study of the human meniscus via poro-elasticity equations.Ce manuscrit traite d’estimation d’erreur a posteriori pour la discrétisation d’équations aux dérivées partielles standard et fractionnaires par les méthodes éléments finis ainsi que de l’application de l’analyse fractionnaire à la modélisation du ménisque humain par les équations de poro-élasticité. Dans l’introduction, nous donnons un aperçu de la littérature sur l’estimation d’erreur a posteriori pour les méth- odes éléments finis et des méthodes de raffinement adaptatif. Nous insistons particulièrement sur l’état de l’art de la méthode d’estimation d’erreur a posteriori de Bank-Weiser et sur les résultats de convergence des méthodes adaptatives. Ensuite, nous nous intéressons aux équations aux dérivées partielles fractionnaires. Nous présentons certaines méthodes de discrétisation d’équations basées sur l’opérateur Laplacien fractionnaire et donnons l’état de l’art sur l’estimation d’erreur a posteriori. Finalement, nous donnons un aperçu de la littérature concernant les applications de la dérivée fractionnaire au sens de Caputo en nous concentrant sur le phénomène de diffusion anormale et les applications en poro-élasticité
    • …
    corecore