[Getdp] quelques points à éclaircir

Christophe Geuzaine geuzaine at acm.caltech.edu
Wed Feb 5 19:31:53 CET 2003


Ali ABAKAR wrote:
> 
> Bonjour Patrick et merci pour tes réponses, nous avons bien avancé
> mais Il reste encore quelques points à éclaircir.
> Pour ce qui est de notre couplage circuit, ce n'est pas un vrai couplage
> avec
> circuit extérieur. Il s'agit en fait de faire une résolution avec une
> condition limite
> ensuite, le résultat (courant d'électrode) obtenu permet de recalculer
> les nouvelles
>  conditions limites et aisi de suite ... jusqu'à ce que le courant
> d'électrode ne varie plus.

Comme la grandeur que tu veux imposer est une grandeur globale, tu
pourrais simplement lancer getdp a partir d'un programme externe (e.g.
un shell script) qui modifierait a chaque "iteration" un fichier
contenant les courants d'electrode.

Le script pourrait ressembler ceci (si tu remplaces l'initialisation
explicite des courants dans 'four_3_electrodes.pro' par l'inclusion d'un
fichier, e.g 'Include "param.pro";') :

I1 = I2 = I3 = 0

while (1) {

  creer 'param.pro' contenant 'I1 = $I1;', etc.

  lancer 'getdp four_3_electrodes.pro -solve EleKin_vf -pos I'  

  parser le fichier de post-processing afin de recuperer les nouvelles
  valeurs de I1, I2, etc.

  comparer les nouvelles valeurs avec les anciennes

}

Il est egalement possible de lancer une resolution en pre-processing (a
l'interieur meme de getdp, cf. "AssignFromResolution") pour calculer les
valeurs a affecter a une contrainte, mais ce mecanisme est plutot
destine a etre effectue une seule fois pour calculer des grandeurs
locales (i.e. qu'il serait difficile/impossible de traiter en
"externe").

> J'ai vu qu'il y a une commande printf (dans GetDP), peux tu me dire
> comment l'utiliser pour
> afficher ou imprimer dans un fichier qques variables?
> Pour la commande printf de gmsh est il possible de diriger le print sur
> un fichier de sortie?

Il n'y a malheureusement actuellement aucune option pour specifier un
fichier de sortie.

Christophe

-- 
Christophe Geuzaine

Tel: (626) 395-4552    http://www.geuz.org
Fax: (626) 578-0124    mailto:geuzaine at acm.caltech.edu