Section: New Results

Multirate methods

Participants : Alain Dervieux, Bruno Koobus, Emmanuelle Itam, Stephen Wornom.

This study is performed in collaboration with IMAG-Montpellier II. It addresses an important complexity issue in unsteady mesh adaptation and takes place in the work done in the ANR Maidesc. Unsteady high-Reynolds computations are strongly penalized by the very small time-step imposed by accuracy requirements on regions involving small space-time scales. Unfortunately, this is also true for sophisticated unsteady mesh adaptive calculations. This small time-step is an important computational penalty for mesh adaptive methods of AMR type. This is also the case for the Unsteady Fixed-Point mesh-adaptive methods developed by Ecuador in cooperation with the Gamma3 team of Inria-Saclay. In the latter method, the loss of efficiency is even more crucial when the anisotropic mesh is locally strongly streched. This loss is evaluated as limiting the numerical convergence order for discontinuities to 8/5 instead of second-order convergence. An obvious remedy is to design time-consistent methods using different time steps on different parts of the mesh, as far as they are efficient and not too complex. The family of time-advancing methods in which unsteady phenomena are computed with different time steps in different regions is referred to as the multirate methods. In our cooperation with university of Montpellier, a novel multirate method using cell agglomeration has been designed and developed in our AIRONUM CFD platform. A series of large-scale test cases show that the new method is much more efficient than an explicit method, while retaining a similar time accuracy over the whole computational domain. The comparison with an implicit scheme shows that the implicit scheme is in some cases one order less accurate due to higher time steps and higher dissipation. A communication has been presented at ECCOMAS [17] and an article is submitted to a journal.