Fleurianne Bertrand

Fleurianne Bertrand, Jakub Wiktor Both, Tugay Dağlı

We study the fully mixed formulation of the Biot equations, which is characterized by a symmetric coupling between flow and deformation. This structure enables the use of stable mixed finite elements for each subproblem without a strong compatibility condition across the two subphysics. To exploit this flexibility while preserving the conservation structure of both subproblems, we consider fully mixed finite element methods in which the symmetry of the elastic stress tensor is enforced weakly. The resulting mixed formulation exhibits a saddle-point structure whose stability is determined by suitable inf--sup conditions. Inf--sup stability is established for several families of discrete spaces of arbitrary order, leading to optimal a priori error estimates. Iterative splitting strategies following the classical fixed-stress split with additional tuning are specifically investigated for the fully mixed formulation, with proof of convergence and rates depending on the coupling strength. Contrary to previous analyses on coupled problems with a symmetric structure, we theoretically prove the efficacy of negative stabilization, consistent with Schur-complement ideas. Numerical results based on analytical solutions and the classical Mandel problem support the theory.

Fleurianne Bertrand, Marcel Moldenhauer, Gerhard Starke

A stress equilibration procedure for hyperelastic material models is proposed andanalyzed in this paper. Based on the displacement-pressure approximation computed with a stable finite element pair, it constructs, in a vertex-patch-wise manner, an $H(div)$-conforming approximation to the first Piola-Kirchhoff stress. This is done in such a way that its associated Cauchy stress is weakly symmetric in the sense that its anti-symmetric part is zero tested against continuous piecewise linear functions. Our main result is the identification of the subspace of test functions perpendicular to the range of the local equilibration system on each patch which turn out to be rigid body modes associated with the current configuration. Momentum balance properties are investigated analytically and numerically and the resulting stress reconstruction is shown to provide improved results for surface traction forces by computational experiments.

Fleurianne Bertrand

A first-order system least squares formulation for the sea-ice dynamics is presented. In addition to the displacement field, the stress tensor is used as a variable. As finite element spaces, standard conforming piecewise polynomials for the displacement approximation are combined with Raviart-Thomas elements for the rows in the stress tensor. Computational results for a test problem illustrate the least-squares approach.

Fleurianne Bertrand, Marcel Moldenhauer, Gerhard Starke

The nonconforming triangular piecewise quadratic finite element space by Fortin and Soulie can be used for the displacement approximation and its combination with discontinuous piecewise linear pressure elements is known to constitute a stable combination for incompressible linear elasticity computations. In this contribution, we extend the stress reconstruction procedure and resulting guaranteed a posteriori error estimator developed by Ainsworth, Allendes, Barrenechea and Rankin \cite{AinAllBarRan:12} and by Kim \cite{Kim:12a} to linear elasticity. In order to get a guaranteed reliability bound with respect to the energy norm involving only known constants, two modifications are carried out: (i) the stress reconstruction in next-to-lowest order Raviart-Thomas spaces is modified in such way that its anti-symmetric part vanishes in average on each element; (ii) the auxiliary conforming approximation is constructed under the constraint that its divergence coincides with the one for the nonconforming approximation. An important aspect of our construction is that all results hold uniformly in the incompressible limit. Local efficiency is also shown and the effectiveness is illustrated by adaptive computations involving different Lamé parameters including the incompressible limit case.