Dr. Gatica-Simpertigue, Luis

2016, Dr. Gatica-Simpertigue, Luis, Gatica, Gabriel, Sequeira, Filánder

In this paper we present the apriori and aposteriori error analyses of a non-standard mixed finite element method for the linear elasticity problem with non-homogeneous Dirichlet boundary conditions. More precisely, the approach introduced here is based on a simplified interpretation of the pseudostress–displacement formulation originally proposed in Arnold and Falk (1988), which does not require symmetric tensor spaces in the finite element discretization. In addition, physical quantities such as the stress, the strain tensor of small deformations, and the rotation, are computed through a simple postprocessing in terms of the pseudostress variable. Furthermore, we also introduce a second elementby-element postprocessing formula for the stress, which yields an optimally convergent approximation of this unknown with respect to the broken H(div)-norm. We apply the classical Babuška–Brezzi theory to prove that the corresponding continuous and discrete schemes are well-posed. In particular, Raviart–Thomas spaces of order k ≥ 0 for the pseudostress and piece wise polynomials of degree≤ k for the displacement can be utilized. Moreover, were mark that in the 3D case the number of unknowns behaves approximately as 9 times the number of elements (tetrahedra) of the triangulation when k = 0. This factor increases to 12.5 when one uses the classical PEERS. Next, we derive a reliable and efficient residual-based a posteriori error estimator for the mixed finite element scheme. Finally, several numerical results illustrating the performance of the method, confirming the theoretical properties of the estimator, and showing the expected behaviour of the associated adaptive algorithm, are provided.

2023, Dr. Caucao-Paillán, Sergio, Dr. Gatica-Simpertigue, Luis, Gatica, Gabriel

We propose and analyze a new mixed finite element method for the nonlinear problem given by the stationary convective Brinkman–Forchheimer equations. In addition to the original fluid variables, the pseudostress is introduced as an auxiliary unknown, and then the incompressibility condition is used to eliminate the pressure, which is computed afterwards by a postprocessing formula depending on the aforementioned tensor and the velocity. As a consequence, we obtain a mixed variational formulation consisting of a nonlinear perturbation of, in turn, a perturbed saddle point problem in a Banach spaces framework. In this way, and differently from the techniques previously developed for this model, no augmentation procedure needs to be incorporated into the formulation nor into the solvability analysis. The resulting non-augmented scheme is then written equivalently as a fixed-point equation, so that recently established solvability results for perturbed saddle-point problems in Banach spaces, along with the well-known Banach–Nečas–Babuška and Banach theorems, are applied to prove the well-posedness of the continuous and discrete systems. The finite element discretization involves Raviart–Thomas elements of order for the pseudostress tensor and discontinuous piecewise polynomial elements of degree for the velocity. Stability, convergence, and optimal a priori error estimates for the associated Galerkin scheme are obtained. Numerical examples confirm the theoretical rates of convergence and illustrate the performance and flexibility of the method. In particular, the case of flow through a 2D porous media with fracture networks is considered.

2015, Dr. Gatica-Simpertigue, Luis, Gatica, Gabriel, Sequeira, Filánder

We present a non-standard mixed finite element method for the linear elasticity problem in Rn with non-homogeneous Dirichlet boundary conditions. More precisely, our approach his based on a simplified interpretation of the pseudo stress–displacement formulation originally proposed in Arnold and Falk (1988), which does not require symmetric tensor spaces in the finite element discretization. We apply the classical Babuˇ ska–Brezzi theory to prove that the corresponding continuous and discrete schemes are well-posed. In particular, Raviart–Thomas spaces of orderk≥0 for the pseudo stress and piece wise polynomials of degree ≤k for the displacement can be utilized. In addition, complementing the results in the afore mentioned reference, we introduceanewpostprocessingformulaforthestressrecoveringtheoptimally convergent approximation of the broken H(div)-norm. Numerical results confirm our theoretical findings.