Dr. Gatica-Simpertigue, Luis

The numerical approximation of hyperelasticity must address nonlinear constitutive laws, geometric nonlinearities associated with large strains and deformations, the imposition of the incompressibility of the solid, and the solution of large linear systems arising from the discretisation of 3D problems in complex geometries. We adapt the three-field formulation for nearly incompressible hyperelasticity introduced in Chavan et al. (2007) to the fully incompressible case. The mixed formulation is of Hu–Washizu type and it differs from other approaches in that we use the Kirchhoff stress, displacement, and pressure as principal unknowns. We also discuss the solvability of the linearised problem restricted to neo-Hookean materials, illustrating the interplay between the coupling blocks. We construct a family of mixed finite element schemes (with different polynomial degrees) for simplicial meshes and verify its error decay through computational tests. We also propose a new augmented Lagrangian preconditioner that improves convergence properties of iterative solvers. The numerical performance of the family of mixed methods is assessed with benchmark solutions, and the applicability of the formulation is further tested in a model of cardiac biomechanics using orthotropic strain energy densities. The proposed methods are advantageous in terms of physical fidelity (as the Kirchhoff stress can be approximated with arbitrary accuracy and no locking is observed) and convergence (the discretisation and the preconditioners are robust and computationally efficient, and they compare favourably at least with respect to classical displacement–pressure schemes).

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.

We propose and analyze an augmented mixed formulation for the time-dependent Brinkman–Forchheimer equations written in terms of vorticity, velocity and pressure. The weak formulation is based on the introduction of suitable least squares terms arising from the incompressibility condition and the constitutive equation relating the vorticity and velocity. We establish existence and uniqueness of a solution to the weak formulation, and derive the corresponding stability bounds, employing classical results on nonlinear monotone operators. We then propose a semidiscrete continuous-in-time approximation based on stable Stokes elements for the velocity and pressure, and continuous or discontinuous piecewise polynomial spaces for the vorticity. In addition, by means of the backward Euler time discretization, we introduce a fully discrete finite element scheme. We prove well-posedness and derive the stability bounds for both schemes, and establish the corresponding error estimates. We provide several numerical results verifying the theoretical rates of convergence and illustrating the performance and flexibility of the method for a range of domain configurations and model parameters.