Numerical discretization of a Darcy-Forchheimer problem coupled with a singular heat equation

In Lipschitz domains, we study a Darcy-Forchheimer problem coupled with a singular heat equation by a nonlinear forcing term depending on the temperature. By singular we mean that the heat source corresponds to a Dirac measure. We establish the existence of solutions for a model that allows a diffusion coefficient in the heat equation depending on the temperature. For such a model, we also propose a finite element discretization scheme and provide an a priori convergence analysis. In the case that the aforementioned diffusion coefficient is constant, we devise an a posteriori error estimator and investigate reliability and efficiency properties. We conclude by devising an adaptive loop based on the proposed error estimator and presenting numerical experiments.Comment: arXiv admin note: text overlap with arXiv:2208.1288

An adaptive stabilized finite element method for the Darcy's equations with pressure dependent viscosities

This work aims to introduce and analyze an adaptive stabilized finite element method to solve a nonlinear Darcy equation with a pressure-dependent viscosity and mixed boundary conditions. We stated the discrete problem's well-posedness and optimal error estimates, in natural norms, under standard assumptions. Next, we introduce and analyze a residual-based a posteriori error estimator for the stabilized scheme. Finally, we present some two- and three-dimensional numerical examples which confirm our theoretical results

Finite element methods for parameter dependent problems

This thesis develops finite element methods for parameter dependent equations. The interest lies in cases where the nature of the problem and the numerical methods used change with the parameters. The studied examples are the reaction-diffusion problem, the Robin boundary condition, which generalizes the Dirichlet and Neumann conditions, and the Brinkman equation, which generalizes the Stokes and the Darcy equations. The developed methods depend continuously on the problem parameters and work even for the limiting values. A posteriori estimates are derived for all the proposed methods, taking into account the parameters. Significant parts of this work are the implementation of the proposed methods and the numerical verfication of the developed theory

A mixed finite element method for Darcy’s equations with pressure dependent porosity

In this work we develop the a priori and a posteriori error analyses of a mixed finite element method for Darcy’s equations with porosity depending exponentially on the pressure. A simple change of variable for this unknown allows to transform the original nonlinear problem into a linear one whose dual-mixed variational formulation falls into the frameworks of the generalized linear saddle point problems and the fixed point equations satisfied by an affine mapping. According to the latter, we are able to show the well-posedness of both the continuous and discrete schemes, as well as the associated Cea estimate, by simply applying a suitable combination of the classical Babuska-Brezzi theory and the Banach fixed point Theorem. In particular, given any integer k ≥ 0, the stability of the Galerkin scheme is guaranteed by employing Raviart-Thomas elements of order k for the velocity, piecewise polynomials of degree k for the pressure, and continuous piecewise polynomials of degree k+1 for an additional Lagrange multiplier given by the trace of the pressure on the Neumann boundary. Note that the two ways of writing the continuous formulation suggest accordingly two different methods for solving the discrete schemes. Next, we derive a reliable and efficient residualbased a posteriori error estimator for this problem. The global inf-sup condition satisfied by the continuous formulation, Helmholtz decompositions, and the local approximation properties of the Raviart-Thomas and Cl´ement interpolation operators are the main tools for proving the reliability. In turn, inverse and discrete inequalities, and the localization technique based on triangle-bubble and edge-bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of both methods, confirming the aforementioned properties of the estimator, and showing the behaviour of the associated adaptive algorithm, are reported.Centro de Investigación en Ingeniería Matemática (CI2MA), Universidad de ConcepciónUniversity of LausanneMinistry of Education, Youth and Sports of the Czech Republi