New Software and Platforms
Bilateral Contracts and Grants with Industry
Partnerships and Cooperations
Bibliography
 PDF e-Pub

## Section: New Results

### Numerical Methods

#### Numerical analysis for an energy-stable total discretization of a poromechanics model with inf-sup stability

Participants : Dominique Chapelle [correspondant] , Philippe Moireau.

In this joint work with Bruno Burtschell [16], we consider a previously proposed general nonlinear poromechanical formulation, and we derive a linearized version of this model. For this linearized model, we obtain an existence result and we propose a complete discretization strategy–in time and space–with a special concern for issues associated with incompressible or nearly-incompressible behavior. We provide a detailed mathematical analysis of this strategy, the main result being an error estimate uniform with respect to the compressibility parameter. We then illustrate our approach with detailed simulation results and we numerically investigate the importance of the assumptions made in the analysis, including the fulfillment of specific inf-sup conditions.

#### Conservative and entropy controlled remap for multi-material ALE simulations

Participant : Patrick Le Tallec.

#### Multipatch isogeometric analysis for complex structures

Participant : Patrick Le Tallec.

This work – done in collaboration with Nicolas Adam (École Polytechnique and PSA) and Malek Zarroug (PSA) – introduces, analyzes and validates isogeometric mortar methods for the solution of thick shells problems which are set on a multipatch geometry. It concerns industrial parts of complex geometries for which the effects of transverse shear cannot be neglected. For this purpose, Reissner-Mindlin model was retained and rotational degrees of freedom (DOF) of the normal are taken into account. A particular attention is devoted to the introduction of a proper formulation of the coupling conditions at patches interfaces, with a particular interest on augmented lagrangian formulations, to the choice and validation of mortar spaces, and to the derivation of adequate integration rules. The relevance of the proposed approach is assessed numerically on various significative examples of industrial relevance. This work has been submitted for publication in an international journal.

#### Mathematical and numerical study of transient wave scattering by obstacles with the Arlequin Method

Participant : Sébastien Imperiale.

In this work [14] we extend the Arlequin method to overlapping domain decomposition technique for transient wave equation scattering by obstacles. The main contribution of this work is to construct and analyze from the continuous level up to the fully discrete level some variants of the Arlequin method. The constructed discretizations allow to solve wave propagation problems while using non-conforming and overlapping meshes for the background propagating medium and the surrounding of the obstacle respectively. Hence we obtain a flexible and stable method in terms of the space discretization – an inf-sup condition is proven – while the stability of the time discretization is ensured by energy identities.

#### Construction and analysis of fourth-order, energy consistent, family of explicit time discretizations for dissipative linear wave equations

Participants : Juliette Chabassier [MAGIQUE-3D] , Julien Diaz [MAGIQUE-3D] , Sébastien Imperiale [correspondant] .

This work and the corresponding article [19], deal with the construction of a family of fourth order, energy consistent, explicit time discretizations for dissipative linear wave equations. The schemes are obtained by replacing the inversion of a matrix, that comes naturally after using the technique of the Modified Equation on the second order Leap Frog scheme applied to dissipative linear wave equations, by explicit approximations of its inverse. The stability of the schemes are studied using an energy analysis and a convergence analysis is carried out. Numerical results in 1D illustrate the space/time convergence properties of the schemes and their efficiency is compared to more classical time discretizations.

#### Energy decay and stability of a perfectly matched layer For the wave equation

Participants : Sébastien Imperiale [correspondant] , Maryna Kachanovska [POEMS] .

We follow a previous work where PML formulations was proposed for the wave equation in its standard second-order form. In the present work [15], energy decay and ${L}^{2}$ stability bounds in two and three space dimensions are rigorously proved both for continuous and discrete formulations with constant damping coefficients. Numerical results validate the theory.

#### A high-order spectral element fast Fourier transform for the poisson equation

Participants : Federica Caforio, Sébastien Imperiale [correspondant] .

The aim of this work [17] is to propose a novel, fast solver for the Poisson problem discretised with High-Order Spectral Element Methods (HO-SEM) in canonical geometries (rectangle in 2D, rectangular parallelepiped in 3D). This method is based on the use of the Discrete Fourier Transform to reduce the problem to the inversion of the symbol of the operator in the frequency space. The proposed solver is endowed with several properties. First, it preserves the efficiency of the standard FFT algorithm; then, the matrix storage is drastically reduced (in particular, it is independent of the space dimension); a pseudo-explicit Singular Value Decomposition (SVD) is used for the inversion of the symbols; finally, it can be extended to non-periodic boundary conditions. Furthermore, due to the underlying HO-SEM discretisation, the multi-dimensional symbol of the operator can be efficiently computed from the one-dimensional symbol by tensorisation.

#### Thermodynamic properties of muscle contraction models and associated discrete-time principles

Participants : François Kimmig, Dominique Chapelle [correspondant] , Philippe Moireau.

#### Mechanical and imaging models-based image registration

Participants : Radomir Chabiniok, Martin Genet [correspondant] .

#### Validation of finite element image registration-based cardiac strain estimation from magnetic resonance images

Participants : Martin Genet [correspondant] , Philippe Moireau.