Calcul erreur juste après la fin des itérations

[pleblanc86]

Bonjour à tous,

Récemment, nous avons au sein du laboratoire implémenté dans le Code Saturne un module de Conduction Ionique permettant de simuler le développement de la Double Couche Électrique à l'interface solide/liquide. (J'avais d’ailleurs effectué un stage de prise en main dans vos locaux en Novembre dernier).

Après pas mal de codage et de débogage le code "calcul" mais lorsque il termine les itérations il nous sort un message d'erreur :

Signal SIGFPE (exception en virgule flottante) intercepté !

Code: Select all

Pile d'appels :
   1: 0x7f57d1f1e4a0 <+0x364a0>                       (libc.so.6)
   2: 0x7f57d2d0210e <+0x44710e>                      (libsaturne.so.0)
   3: 0x7f57d2d05ddb <fvm_convert_array+0x1d9>        (libsaturne.so.0)
   4: 0x7f57d2911402 <+0x56402>                       (libsaturne.so.0)
   5: 0x7f57d29129ea <fvm_to_ensight_export_field+0x5ca> (libsaturne.so.0)
   6: 0x7f57d291659d <fvm_writer_export_field+0xff>   (libsaturne.so.0)
   7: 0x7f57d2a51d3e <cs_post_write_var+0x604>        (libsaturne.so.0)
   8: 0x7f57d2a4f70a <+0x19470a>                      (libsaturne.so.0)
   9: 0x7f57d2a52e4c <cs_post_write_vars+0x5e8>       (libsaturne.so.0)
  10: 0x7f57d2a5467a <pstvar_+0x102>                  (libsaturne.so.0)
  11: 0x7f57d2948efb <caltri_+0x88ef>                 (libsaturne.so.0)
  12: 0x7f57d290b865 <cs_run+0xda9>                   (libsaturne.so.0)
  13: 0x7f57d290bab4 <main+0x15d>                     (libsaturne.so.0)
  14: 0x7f57d1f0976d <__libc_start_main+0xed>         (libc.so.6)
  15: 0x404109     <>                               (cs_solver)
Fin de la pile

On pense que cela vient du post-processing mais voilà déjà un petit moment quand bloque dessus.

De plus dans le listing il marque certains warnings qui ne sont pas trop claire pour moi :

Code: Select all

 [CSDS]           : CV-DIF-TS    1 IT - RES=  0.13953E-05 NORM=  0.45533E+18
@
@ @@ WARNING: [CSDS]   CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE [CSBL]
      ----------------

 [CSBL]           : CV-DIF-TS    1 IT - RES=  0.15872E+13 NORM=  0.21602E+29
@
@ @@ WARNING: [CSBL]   CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE POT_EL_R
      ----------------

 POT_EL_R         : CV-DIF-TS    1 IT - RES=  0.62110E-05 NORM=  0.48742E+03
@
@ @@ WARNING: POT_EL_R CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

En cherchant un peu j'ai vu que c'est le fichier codits.f90 qui généré les messages, il s'occupe de la résolution de l'équation d'advection diffusion avec les termes sources. J'ai un peu de mal à voir comment je vais pouvoir fixer ce problème.

Voilà donc pourquoi je me tourne vers vous, en espérant pouvoir obtenir un début de piste et/ou une orientation qui me permettra de résoudre le problème. Je joint mon listing file à la fin du post

En vous remerciant d’avance
Cordialement
Paul Leblanc

[pleblanc86]

Re,

Nous venons de corriger l'erreur (mauvaise lecture du fichier de donnée), le code fini bien le calcul et s'arrête correctement. Néanmoins nous ne comprenons pas pourquoi on a toujours un Warning, voir exemple extrait du listing :

Code: Select all

   ** SOLVING VARIABLE [CSBL]
      ----------------

 [CSBL]           : CV-DIF-TS    1 IT - RES=  0.79589E-13 NORM=  0.10884E+05
@
@ @@ WARNING: [CSBL]   CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE POT_EL_R
      ----------------

@
@ @@ WARNING: POT_EL_R CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

Merci d'avance pour vos éclaircissements
Cordialement
P.L

[Jacques Fontaine]

This warning appears only when iwarni(ivar) > 0 (verbosity lvl). It indicates that the iterative process on RHS does not converge. This warning is not very important for steady calculations and should disappear when the flow is converged.

you can refer to the theory doc. for more details.
Regards

[pleblanc86]

Merci de votre réponse