Formulation implicite ou explicite d'un problème de dynamique

En simulation numérique, un problème dépendant du temps peut être formulé de manière implicite ou explicite. Un problème dépendant du temps décrit une situation qui évolue ; le système est modélisé à différents instants t discrets appelés « pas de temps ». La méthode explicite consiste à déterminer la solution à t + Δt en fonction de la valeur de la fonction en t. Si la fonction à évaluer s'appelle y(t), alors le problème se formule de la manière suivante : y(t + Δt) = F(y(t)). La méthode d'Euler est une méthode explicite. La méthode implicite consiste à déterminer la solution à t + Δt en résolvant une équation prenant en compte la valeur de la fonction en t et en t + Δt. Le problème se formule de la manière suivante : G(y(t), y(t + Δt)) = 0. Les méthodes de Runge-Kutta sont des méthodes dites implicites-explicites (ou « imex ») car une partie est résolue par une méthode implicite et l'autre par une méthode explicite. La formulation implicite est la formulation la plus simple mais elle est limitée aux problèmes quasi statiques. La formulation explicite permet de modéliser plus finement les phénomènes mais est très gourmande en ressources (nombre d'opérations, durée du calcul, mémoire nécessaire). Elle est donc utilisée pour des simulations de phénomènes de courte durée et ne pouvant être simulés avec la formulation implicite : il s'agit essentiellement de problèmes de dynamique rapide, chocs, propagation d'ondes. En termes de logiciels, on parle de solveur implicite ou de solveur explicite. En formulation implicite, le phénomène est représenté par l'équation où est le vecteur position ; est le vecteur vitesse et le vecteur accélération ; m est la masse ; c est un facteur d'amortissement ; k est la raideur ; F est la force extérieure s'exerçant sur le système ; t est le temps. Pour être plus précis, dans la méthode des éléments finis, cette équation décrit le comportement des éléments, les différents termes de l'équation sont donc des matrices. À chaque pas de temps, le solveur cherche une solution stationnaire à cette équation, cela représente donc une situation d'équilibre.