Cerfacs Enter the world of high performance ...

The 24 March 2021 at 15h00

PhD defense of Adèle VEILLEUX : ” Extension of the spectral difference method to simplex cells and hybrid grids”

Marie LABADENS |  Cerfacs, Toulouse (France) - Salle de réunion Administration en visioconférence |  

Webex link:

access code: 128 243 1066
password: rXN5DAMGj63

Please switch off your camera and micro during the defense except for the jury’s members

This thesis examines the extension of the Spectral Difference (SD) method on unstructured hybrid grids involving simplex cells (triangles, tetrahedra) and prismatic elements. The Spectral Difference method is part of high-order spectral discontinuous numerical methods. These methods rely on piecewise continuous polynomial approximation to obtain high-order accuracy with a good parallel efficiency. The standard SD scheme is first presented in the one-dimensional case and then for tensor-product elements (quadrangles and hexahedra). The treatment of simplex cells using Raviart-Thomas elements is detailed for triangles (in 2D) and tetrahedra (in 3D), followed by the implementation for prismatic elements. 

The linear stability of the Spectral Difference method using Raviart-Thomas elements (SDRT) is studied on triangles and tetrahedra. The SDRT scheme stability is strongly dependent on the interior flux points location. On triangles, the SDRT implementation based on interior flux points located at Williams-Shunn-Jameson quadrature points is found stable up to the fourth-order of accuracy but shown as spatially unstable for higher orders. Nevertheless, it is shown that this implementation can be stabilized for fifth- and sixth-order schemes using suitable temporal integration schemes. This approach being submitted to strict conditions, an optimization of the interior flux points location is conducted to determine spatially stable SDRT formulations for orders higher than four. The optimization process leads to spatially stable schemes up to the sixth-order of accuracy. Finally, the stability analysis on tetrahedra proves that the SDRT scheme based on the interior flux points located at Shunn-Ham quadrature points is stable up to the third-order. 

The SD/SDRT numerical method is validated on several academic cases for first and second-order Partial Differential Equations (linear advection equation, Euler equations, Navier-Stokes equations). Both proposed implementations (based either on Williams-Shunn-Jameson quadrature points or optimization points) are used. Numerical experiments involve grids composed of quadratic triangles, linear tetrahedral elements as well as 2D hybrid meshes.

Keywords : High-order discontinuous method, Spectral Difference, Hybrid grid, Stability analysis, Coupled time-space analysis


Chunlei LIANGClarkson University, Postdam-New York, USAReferee
Georg MAYVon Karman Institute, Sint-Genesius-Rode, BelgiumReferee
Hélène BARUCQINRIA Toulouse, FranceMember
Eusebio VALEROUniversidad Politécnica de Madrid, SpainMember
Hugues DENIAUONERA Toulouse, FranceInvited member
Jean-Luc ESTIVALEZESIMFT Toulouse, FranceInvited member
Guillaume PUIGTONERA Toulouse, FranceAdvisor
Guillaume DAVILLERCERFACS Toulouse, FranceCo Advisor


Thierry Poinsot officially entered the French Academy of Sciences

CERFACS |  8 November 2021

Thierry Poinsot officially entered the French Academy of Sciences on October 12. See presentation here :Read more

The AVBP code from CERFACS at the heart of for PRACE projects from the 23rd call

CERFACS |  30 September 2021

Cerfacs is involved in three PRACE projects of the 23rd call for which hour allocation runs from 01/10/2021 to 30/09/2022. Researchers from ECL/LMFA UMR5509 (Ecole Centrale de Lyon) and IMFT (UMR 5502) laboratories have earned projects entirely based on the use of the LES solver developed by Cerfacs AVBP and involve the support of experts from the CFD and COOP teams underling the importance and effectiveness of collaborations between French labs and Cerfacs. Alexis Giauque from ECL/LMFA UMR5509 (Ecole Centrale de Lyon) has obtained not only one but two PRACE projects! The first project LESFAN (RA0101, 30 000 000 CPU hours on Irene/Rome TGCC) is based on the use of AVBP in the turbomachinery version to study the generation of noise by a fan of a real airplane engine. The second, PRACE-EDGES (RA0101, 40 000 000 CPU hours on Irene/Rome TGCC) focuses on LES modeling of dens gas in complex geometries. To do so, the LMFA Team has developed advanced thermodynamic closures in AVBP allowing the direct simulation of such flows. Laurent Selle from IMFT (UMR 5502) has received CPU hours for the GASTON project (RA0061, 30 000 000 CPU hours on Marenostrum BSC) which aims to study the structure of hydrogen flames in porous materials. For this, IMFT and Cerfacs will perform coupled simulations considering the reactive flow with AVBP as well as the conduction in the porous medium with AVTP which is known to play an central role in the flame stabilization process. Carlos Perez Arroyo from IMFT (UMR 5502) received 16 Mh CPU hours on Joliot-Curie Skylake partition to support his project WONDER.Read more