Section: New Results

Numerical Methods

Effective and energy-preserving time discretization for a general nonlinear poromechanical formulation

Participants : Bruno Burtschell, Dominique Chapelle [correspondant] , Philippe Moireau.

In this work, we consider a general nonlinear poromechanical model, formulated based on fundamental thermodynamics principle, suitable for representing the coupling of rapid internal fluid flows with large deformations of the solid, and compatible with a wide class of constitutive behavior. The objective of the present work is to propose for this model a time discretization scheme of the partitioned type, to allow the use of existing time schemes - and possibly separate solvers - for each component of the model, i.e. for the fluid and the solid. To that purpose, we adapt and extend an earlier proposed approach devised for fluid-structure interaction in an Arbitrary Lagrangian-Eulerian framework. We then establish an energy estimate for the resulting time scheme, in a form that is consistent with the underlying energy principle in the poromechanical formulation, up to some numerical dissipation effects and some perturbations that we have carefully identified and assessed. In addition, we provide some numerical illustrations of our numerical strategy with test problems that present typical features of large strains and rapid fluid flows, and also a case of singular transition related to total drainage. An example of challenging application envisioned for this model and associated numerical coupling scheme concerns the perfusion of the heart. This work has resulted in the publication [15].

Delayed feedback control method for calculating space-time periodic solutions of viscoelastic problems

Participants : Ustim Khristenko, Patrick Le Tallec.

We are interested in fast techniques for calculating a periodic solution to viscoelastic evolution problems with a space-time periodic condition. In order to avoid the invertion of very large matrices, such a solution is often computed as an asymptotic limit of the initial value problem with arbitrary initial data. We have developed a control method, accelerating the convergence to the periodic state. The main idea is to modify our problem by introducing a feedback control term, based on a periodicity error.

First, an abstract evolution problem has been studied. From the analytic solution of the modified (controlled) problem, an efficient control has been constructed, optimizing the spectrum of the problem. The proposed control term can be mechanically interpreted, and its efficiency increases with the relaxation time.

In order to confirm numerically the theoretical results, a finite element simulation has been carried out on a full 3D model for a steady rolling of a viscoelastic tyre with periodic sculpture. It has demonstrated that the controlled solution converges indeed faster than the non-controlled one, and that the efficiency of the method increases with the problem's relaxation time, that is when the memory of the underlying problem is large.

Construction and analysis of an adapted spectral finite element method to convective acoustic equations

Participant : Sébastien Imperiale [correspondant] .

This work addresses the construction of a non spurious mixed spectral finite element (FE) method to problems in the field of computational aeroacous-tics. Based on a computational scheme for the conservation equations of linear acoustics, the extension towards convected wave propagation is investigated. In aeroacoustic applications, the mean flow effects can have a significant impact on the generated sound field even for smaller Mach numbers. For those convective terms, the initial spectral FE discretization leads to non-physical, spurious solutions. Therefore, a regularization procedure is proposed and qualitatively investigated by means of discrete eigenvalues analysis of the discrete operator in space. A study of convergence and an application of the proposed scheme to simulate the flow induced sound generation in the process of human phonation underlines stability and validity. This work has resulted in the publication [19].

Space/time convergence analysis of a class of conservative schemes for linear wave equations

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

This work concerns the space/time convergence analysis of conservative two-steps time discretizations for linear wave equations. Explicit and implicit, second and fourth order schemes are considered, while the space discretization is given and satisfies minimal hypotheses. The convergence analysis is done using energy techniques and holds if the time step is upper-bounded by a quantity depending on space discretization parameters. In addition to showing the convergence for recently introduced fourth order schemes, the novelty of this work consists in the independency of the convergence estimates with respect to the difference between the time step and its greatest admissible value. This work has resulted in the publication [16].