La réponse courte est: cela nécessite un travail spécifique pour différentes équations, mais il existe certaines techniques générales qui suggèrent comment le faire. Essentiellement, étant donné une PDE d'évolution de premier ordre
ut= A u + B u
A , B
A u + B u = 0.
UNEBUNEB
Une approche que j'aime est l'utilisation de solveurs Riemann à ondes f comme proposé par Bale et. Al. . L'idée est de discrétiser les termes convectifs avec une méthode de type Godunov, mais de soustraire la contribution des autres termes à l' intérieur du solveur de Riemann. Dans le cas de l'état stationnaire, aucune onde n'est générée. Cependant, cela nécessite que les termes convectifs et sources soient calculés exactement (afin d'annuler exactement). C'est possible pour les équations en eau peu profonde, mais plus difficile pour de nombreux autres systèmes.
David Ketcheson
la source