Résolution des équations intégrales des fluides à potentiels intermoléculaires anisotropes par l'algorithme Général de Minimisation du RESte

Publisher Website
Google Scholar
Abstract
Après avoir rappelé les principales méthodes de résolution numérique des équations intégrales donnant les fonctions de distribution de paires on propose d'appliquer le nouvel algorithme Général de Minimisation du RESte étendu aux systèmes d'équations Non Linéaires (GMRESNL). Cet algorithme, de type Newton-Raphson, est général, stable, efficace et très simple à programmer. Il est utilisé pour résoudre l'approximation Relative Percus-Yevick (RPY) dans le cas d'un fluide pur dense de sphères dures à grands moments dipolaires où les méthodes itératives traditionnelles de type Picard ne convergent pas. Les résultats obtenus sont comparés aux données Monte-Carlo et à celles d'autres théories d'équations intégrales. L'approximation RPY apparaît comme bien moins précise que la théorie Relative des chaînes hypertressées (RHNC) pour ce fluide à potentiel intermoléculaire à relativement longue portée. On analyse enfin les propriétés de convergence de l'algorithme GMRESNL dans le cas des équations intégrales de la théorie des liquides