Précédent : GENIE phase I - Remonter : Optimisation de forme et différentiation Suivant : Perspectives -

Contrat relais INRIA -Dassault -

Dans le cadre de ce contrat relais, nous nous sommes engagés à réaliser une première version d'une maquette de logiciel d'optimisation 3D avec un gradient d'ordre deux. Pour ce faire, il est nécessaire de réécrire le noyau du programme de résolution (calculs des flux) de la même manière qu'en 2D. Il est également nécessaire de calculer les dérivées de ces flux par rapport aux valeurs des inconnues sur les points du maillage.

C'est ici que les différences sensibles apparaissent entre le 2D et le 3D. En trois dimensions d'espace, il est hors de question de stocker une matrice adjointe sur un maillage raisonnable, même en ne stockant que les termes non nuls. Il faut donc modifier notre façon de procéder et faire générer par ODYSSÉE, non plus , mais , où nous avons noté le flux élémentaire sur un segment du maillage, et W les inconnues dont dépend . Cette modification nécessite la réécriture d'un solveur linéaire sans stockage de type Jacobi multipas ; ces deux tâ ches sont en cours de réalisation. D'autre part, une version parallèle du programme 3D est en cours de mise au point avec S. Lanteri ; elle fonctionne avec un gradient d'ordre 1 et une résolution d'état d'ordre 2.