90 likes | 159 Vues
Detailed improvements to hydrogrid code for saltwater simulation, focusing on time/memory efficiency, code modification, parallelism, and visualization. Includes case tests, meshing, and finite elements.
E N D
Le code d´eau saléeréunion hydrogridle 11/10/2004 E. Canot J. Erhel C. De Dieuleveult
Sommaire • Présentation générale • Gain de temps et de mémoire • Modification du code • Parallélisme/MUMPS • Réécriture du code • Modification du code • CORBA • OUTILS • Cas tests • Mailleur & Visualisation
Présentation générale • Écoulement + Transport • Eléments Finis, schéma point fixe • Maillage 2D rectangulaire et régulier • ρ variable • εconstante • S=0 • Pas de terme source et de gradient de concentration pour le transport.
Gain de temps et de Mémoire • Modification du code • Allocation dynamique des tableaux • Calcul de nom à l´initialisation • Revenir sur le calcul de TP puis de P • Modification du couplage écoulement et transport • MUMPS/Parallélisme
Gain de temps et de mémoireMUMPS / Parallélisme • Modification du parallélisme, MUMPS dans le transport, analyse à l´initialisation • test de parallélisme • Sur le code entier ( + METIS, renumérotation) • Mauvais résultats sur les machines du projet, paraci, idefix • Transport : • Mauvais sur les paraci (TPC/IP) • Bon résultat à Lyon (Myrinet) • Tester le code sur le réseau Myrinet des paraci avec MPICH-GM • Distribution des données (METIS)
Réécriture du code • Modification du code de base • caldis →integrate_quad (de Trace-1.2) • rmprecgc • calmtc • mat, matc à changer • CORBA • Création d´une première version • Validation du code (performance) • Parallélisation de la visualisation…
Réécriture du code - CORBA gestionnaire 1 Gestion boucle sur les pas de temps boucle sur les itérations ecoulement → contrôle transport → contrôle convergence? OUI- affichage? OUI- log → appel NON- divergence? OUI fin fin 4 2 Pas de temps suivant On recommence le même pas de temps ecoulement transport Contrôle convergence? modifications variables calcul transport → donnee calcul donnée convergence Donnee Affichage Contrôle convergence? modifications variables calcul ecoulement → donnee calcul donnée convergence Donnee Affichage 5 3 6 log Appel ecoulement → affichage transport → affichage affichage et sauvegarde 7 8
Outils • Cas tests Issus de l´article de P.Ackerer, A. Younes et R. Mose : « Modeling Variable Density Flow and Solute Transport in Porous Medium : 1. Numerical Model and Verification. » • Paramètres d´entrées et de sorties • Unités SI
Outils • Mailleur • Version fortran ou Matlab • Adapté à de gros maillages • CL linéaire en charge • Possibilité de partitionnement • Visualisation • Sauvegarde : Images GIF ou EPS • Dégradés de couleurs et lignes de niveaux • Stretch sur la colormap