The simulation of the Navier-Stokes equations requires accurate and robust numerical methods that
take into account diffusion operators, gradient and convection terms. Operational approaches have
shown their effectiveness on simplexes. However, in some models or codes
(TrioCF, Flica5), it may be useful to improve the accuracy of solutions locally using an
error estimator or to take into account general meshes. We are here interested in staggered schemes.
This means that the pressure is calculated at the centre of the mesh and the velocities on the edges
(or faces) of the mesh. This results in methods that are naturally accurate at low Mach numbers .
New schemes have recently been presented in this context and have shown their
robustness and accuracy. However, these discretisations can be very costly in terms of memory and
computation time compared with MAC schemes on regular meshes
We are interested in the "gradient" type methods. Some of them are based on a
variational formulation with pressure unknowns at the mesh centres and velocity vector unknowns on
the edges (or faces) of the cells. This approach has been shown to be effective, particularly in terms of
robustness. It should also be noted that an algorithm with the same degrees of freedom as the
MAC methods has been proposed and gives promising results.
The idea would therefore be to combine these two approaches, namely the "gradient" method with the same degrees of freedom as MAC methods. Initially, the focus will be on recovering MAC schemes on regular meshes. Fundamental
questions need to be examined in the case of general meshes: stability, consistency, conditioning of
the system to be inverted, numerical locking. An attempt may also be made to recover the gains in
accuracy using the methods presented in for discretising pressure gradients.
During the course of the thesis, time will be taken to settle the basic problems of this method (first and
second years), both on the theoretical aspects and on the computer implementation. It may be carried
out in the Castem, TrioCFD, Trust or POLYMAC development environments. The focus will be on
application cases that are representative of the community.