Étude de méthodes d'intégration temporelle basées sur des schémas de corrections différées d'ordres élevés et A-stables

Dans le domaine de la simulation numérique en mécanique des fluides, deux aspects sont à prendre en compte : la discrétisation spatiale et la discrétisation temporelle. Ces deux aspects vont de pair et influent tous deux sur la précision et le coût d'un calcul. Dans notre travail, nous nous intéressons à la discrétisation temporelle et plus particulièrement aux méthodes BDF – Backward Differentiations Formulae – car elles permettent d'effectuer des calculs stables dans des cas de problèmes raides dont font notamment partie les problèmes de Navier-Stokes. Mais ces méthodes sont stables pour des ordres de convergence de un et deux seulement. Pour augmenter l'ordre de convergence, nous construisons une méthode de correction différée, c'est-à-dire que nous modifions le schéma de discrétisation avec une correction. Cette correction est basée sur l'erreur de troncature propre à chaque BDF. Nous construisons ainsi une méthode que nous nommons DC5/BDF1 qui est un schéma BDF1 corrigé successivement quatre fois afin d'obtenir une méthode stable d'ordre cinq. De même, nous construisons une méthode DC5/BDF2, également à l'ordre cinq, mais basée cette fois sur le schéma BDF2. De plus nous construisons ces méthodes pour qu'elles soient fonctionnelles, c'est-à-dire stables et avec le bon taux de convergence, à pas de temps alterné. En effet le but est de faire un premier pas en direction d'un intégrateur permettant de fonctionner avec des pas de temps adaptatifs. Nous testons ensuite nos méthodes pour des EDO simples et des systèmes d'EDO. Nous vérifions que les taux de convergence sont bien atteints. Nous discutons de la raideur des équations ainsi que du pas de temps utilisé pour les simulations. Ensuite des tests sur des EDP en une dimension sont effectués pour comprendre l'influence des conditions limites sur nos intégrateurs. Nous remarquons alors que des pertes d'ordre peuvent se produire sur des corrections successives, et cela en fonction des conditions limites et en utilisant des pas de temps alternés. Pour remédier à ce problème, nous développons donc une méthode pour imposer différemment les conditions limites : la correction de Verwer. Précédemment nous utilisions des conditions de Dirichlet, le principe est maintenant d'imposer faiblement une équation différentielle du premier ordre aux bords de notre domaine. Pour le faire directement au sein des formes faibles, nous utilisons un multiplicateur de Lagrange. Cette correction permet de remédier au problème de perte d'ordre et fournit également des informations supplémentaires sur la simulation via le multiplicateur de Lagrange. Dans un problème de thermique, le multiplicateur correspond au flux de chaleur au bord. Dans un problème d'écoulement isotherme, le multiplicateur nous donne la résultante des forces aux frontières. Cette information est très utile dans le cadre d'une application aux Interactions Fluide-Structure, car cela évite le difficile calcul des forces en post-traitement. Nous fournissons également des résultats en deux dimensions pour des problèmes de diffusion thermique et d'écoulements isothermes. Cela nous permet de discuter de l'effet de la correction de Verwer. Ces résultats nous permettent de valider nos méthodes. Enfin nous effectuons une simulation sur des tourbillons de Von Karman. Pour ce cas nous devons combiner nos intégrateurs DC/BDF1 et DC/BDF2 afin d'initialiser correctement la méthode DC/BDF2. Les résultats obtenus comparés à ceux de la littérature nous permettent d'apporter un élément de vérification à nos méthodes.