Page 1 of 1

Inlets problems on electrical arc model

Posted: Wed Aug 03, 2011 8:44 am
by Christophe Rehmet
Hello everybody,
 
I have a problem when I introduce 3 differents Inlets in arc elect module. My boundaries conditions works in simple hydrodynamics models but not in arc elect study.

@@ ATTENTION : ARRET LORS DE LA VERIFICATION DES COND. LIM.
@    =========                                              
@    CONDITIONS AUX LIMITES INCORRECTES OU INCOMPLETES       @                                                           
@    Au moins une face de bord declaree en entree           
@      (ou sortie) a vitesse imposee avec un flux rentrant  
@      pour laquelle la valeur d'une variable n'a pas ete 
@      specifiee (condition de Dirichlet).                  
@    Le calcul ne sera pas execute                           @                                                           
@    Verifier les conditions aux limites dans l'Interface
@    ou dans le sous-programme utilisateur correspondant.    @                                                            @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                           

Premire face avec une mauvaise définition de conditions aux limites
  (sur 2560)
  a une condition aux limites de type 2, un centre (-0.018, -0.0123563, -0.000485482)
/usr/local/saturne/ncs-2.0.0-rc2/src/base/cs_boundary_conditions.c:429: Erreur fatale.


Moreover when i change the direction of the 2 inlets to outlet, it works.  I don't know what for Dirichlet conditions I need to introduce in my code.

Thank.

Re: Inlets problems on electrical arc model

Posted: Wed Aug 03, 2011 12:00 pm
by Mickael Hassanaly
Dear Christophe,
 
Could you send us your mesh and your simulation set up (f90 routines + *.xml)?
 
Best regards
 
Mickaël Hassanaly
 

Re: Inlets problems on electrical arc model

Posted: Sat Aug 06, 2011 9:11 am
by Christophe Rehmet
Dear Mickaël,

I found my mystake, I forgot to introduce at all Inlets:

 if ( ippmod(ielarc).ge.1 ) then
    do idimve = 1, ndimve
      do iel = 1, ncel
        rtp(iel,isca(ipotva(idimve))) = 0.d0
      enddo
    enddo
  endif
 
Best regards.
 
Christophe