Précédent : Les modèles de transition de
Remonter : Problèmes à frontières librescontrôle
Suivant : Contrôle optimal et équation
d'Hamilton-Jacobi
Participants : Olivier Coulaud, Philippe Laurençot, Bruno Pinçon
Il s'agit ici d'un investissement récent en collaboration avec
les chercheurs du laboratoire de Chimie Théorique de l'Université
de Nancy I. Il relève de la modélisation en génie
moléculaire. Dans un premier temps, nous nous sommes concentrés
sur la simulation de molécules plongées dans un solvant
et, plus particulièrement, sur le calcul du potentiel électrique.
Les calculs doivent conduire à déterminer la déformation de la
molécule et le potentiel intermoléculaire à la surface de la
molécule. La situation est fortement changée par rapport au vide
et soulève des difficultés nouvelles. Deux d'entre elles ont pu
être résolues :
- L'utilisation d'un programme existant (dû à M. Sanner)
nous a permis de réaliser la triangulation de ``grosses
molécules'' ; à partir de ces triangulations que nous
améliorons un peu en les régularisant, nous créons une
triangulation P2 ``courbe''.
- Nous avions initialement programmé deux méthodes de résolution
de léquation intégrale décrivant le potentiel électrique. Notre
appoche donne maintenant les résultats attendus :
convergence plus rapide que la méthode élémentaire et les
éléments courbes nous permettent de jouer sur le maillage sans
tronquer la forme de la molécule (même s'il reste beaucoup à
faire pour aboutir à un véritable maillage adaptatif). Un rapport
[77] sera soumis
prochainement pour publication.
Les collaborations avec le Laboratoire de Chimie Théorique de Nancy doivent encore s'intensifier. Des contacts ont été aussi établis avec M. Karplus de Strasbourg (membre du Conseil Scientifique de l'Inria).