Intégrateurs exponentiels modifiés pour la simulation des vagues non linéaires

Pour réaliser des simulations précises aux temps longs pour des vagues non linéaires, il faut faire appel à des algorithmes d'évolution temporelle précis. En particulier, la combinaison d'un pas de temps adaptatif avec un facteur intégrant est connue pour être très efficace. Nous proposons...

Full description

Bibliographic Details
Main Author: Eichwald, Brice
Language:FRE
Published: Université Nice Sophia Antipolis 2013
Subjects:
Online Access:http://tel.archives-ouvertes.fr/tel-00873578
http://tel.archives-ouvertes.fr/docs/00/87/35/78/PDF/2013NICE4048.pdf
Description
Summary:Pour réaliser des simulations précises aux temps longs pour des vagues non linéaires, il faut faire appel à des algorithmes d'évolution temporelle précis. En particulier, la combinaison d'un pas de temps adaptatif avec un facteur intégrant est connue pour être très efficace. Nous proposons une modification de cette technique. Le principe consiste à soustraire un certain polynôme à une EDP. Puis, comme pour le facteur intégrant, nous faisons un changement de variable pour retirer la partie linéaire. Mais nous espérons retirer quelque chose de plus afin de rendre l'EDP moins raide pour les calculs numériques. Le polynôme choisi est une expansion de Taylor autour du temps initial de la solution. Afin de calculer les différentes dérivées nécessaires, nous utilisons le Dense Output qui donne la possibilité d'approximer les dérivées de la solution à tout temps. Une fois le facteur intégrant modifié appliqué, nous faisons appel à une avance temporelle classique afin de résoudre l'équation d'évolution. Il a été considéré plusieurs schémas de Runge-Kutta avec pas de temps adaptatif. Nous avons tiré avantage des méthodes emboîtées, afin de ne pas calculer de nouvelles fonctions et perdre du temps de calcul, en utilisant uniquement des données déjà calculées durant l'évolution temporelle. Les résultats numériques montrent que l'efficacité de notre méthode varie selon les cas. Par exemple, nous avons vérifié que plus le profil de l'onde est pentue, plus notre méthode est efficace. Pour le modèle de vagues non linéaires le plus compliqué à notre disposition, le modèle HOS, nous avons pu réduire le nombre de pas de temps de calcul jusqu'à près de 30 % avec un schéma de Runge-Kutta de Dormand-Prince et jusqu'à plus de 99 % pour un schéma de Bogacki-Shampine.