!------------------------------------------------------------------------------- ! Code_Saturne version 3.0.0 ! -------------------------- ! This file is part of Code_Saturne, a general-purpose CFD tool. ! ! Copyright (C) 1998-2013 EDF S.A. ! ! This program is free software; you can redistribute it and/or modify it under ! the terms of the GNU General Public License as published by the Free Software ! Foundation; either version 2 of the License, or (at your option) any later ! version. ! ! This program is distributed in the hope that it will be useful, but WITHOUT ! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS ! FOR A PARTICULAR PURPOSE. See the GNU General Public License for more ! details. ! ! You should have received a copy of the GNU General Public License along with ! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin ! Street, Fifth Floor, Boston, MA 02110-1301, USA. !------------------------------------------------------------------------------- subroutine uselrc & !================ ( nvar , nscal , & dt , rtpa , rtp , propce , propfa , propfb ) !=============================================================================== ! FONCTION : ! ---------- ! ROUTINE PHYSIQUE PARTICULIERE POUR LE MODULE ELECTRIQUE ! CALCULS DU COEFFICIENT DE RECALAGE ! POUR LES VARIABLES ELECTIQUES ! RECALAGE DES VARIABLES ELECTRIQUES ! EN FONCTION DE CE COEFFICIENT !------------------------------------------------------------------------------- !ARGU ARGUMENTS !__________________.____._____.________________________________________________. ! name !type!mode ! role ! !__________________!____!_____!________________________________________________! ! nvar ! i ! <-- ! total number of variables ! ! nscal ! i ! <-- ! total number of scalars ! ! itypsm ! te ! <-- ! type de source de masse pour les ! ! (ncesmp,nvar) ! ! ! variables (cf. ustsma) ! ! dt(ncelet) ! ra ! <-- ! time step (per cell) ! ! rtp, rtpa ! ra ! <-- ! calculated variables at cell centers ! ! (ncelet, *) ! ! ! (at current and previous time steps) ! ! propce(ncelet, *)! ra ! <-- ! physical properties at cell centers ! ! propfa(nfac, *) ! ra ! <-- ! physical properties at interior face centers ! ! propfb(nfabor, *)! ra ! <-- ! physical properties at boundary face centers ! ! smacel ! tr ! <-- ! valeur des variables associee a la ! ! (ncesmp,* ) ! ! ! source de masse ! ! ! ! ! pour ivar=ipr, smacel=flux de masse ! !__________________!____!_____!________________________________________________! ! Type: i (integer), r (real), s (string), a (array), l (logical), ! and composite types (ex: ra real array) ! mode: <-- input, --> output, <-> modifies data, --- work array !=============================================================================== !=============================================================================== ! Module files !=============================================================================== use paramx use numvar use entsor use optcal use cstphy use cstnum use parall use period use ppppar use ppthch use ppincl use elincl use mesh !=============================================================================== implicit none ! Arguments integer nvar , nscal double precision dt(ncelet), rtp(ncelet,*), rtpa(ncelet,*) double precision propce(ncelet,*) double precision propfa(nfac,*), propfb(nfabor,*) ! Local variables integer iel , ifac , iutile integer ipcefj , ipcdc1 , ipcdc2 , ipcdc3 , ipcsig integer ipdcrp , idimve , jaiex integer dtimes double precision somje , coepoa , coefav , coepot double precision emax , aiex , amex double precision rayo , econs , z1 , z2 , posi double precision dtj , dtjm , delhsh , cdtj , cpmx double precision xelec , yelec , zelec, diff double precision, allocatable, dimension(:) :: w1 !=============================================================================== !=============================================================================== !=============================================================================== ! 1. INITIALISATION !=============================================================================== !=============================================================================== ! 2. ARC ELECTRIQUE !=============================================================================== if ( ippmod(ielarc).ge.1 ) then ! 2.1 : exemple : cas avec claquage ! ======================================= ! Ceci est un cas particulier et doit etre adapte en fonction ! du cas et du maillage (intervenir aussi dans uselcl) ! Utilisation d'une rampe d'intensite ! ----------------------------------- dtimes = ntcabs if ( ntcabs.le.200 ) then couimp = 5.d0 + (100.d0-5.d0) * ntcabs/200 else couimp = 100.d0 endif ! Calcul de l'integrale sur le Volume de J.E ! ----------------------------------- ! (c'est forcement positif ou nul) ipcefj = ipproc(iefjou) somje = 0.d0 do iel = 1, ncel somje = somje+propce(iel,ipcefj)*volume(iel) enddo if(irangp.ge.0) then call parsom (somje) endif if (somje .ne. 0) then coepot = couimp*dpot/max(somje,epzero) endif write(nfecra,1001) couimp,dpot,somje ! Calcul de l'intensite du courant d'arc ! -------------------------------------- ! Calcul de l'integrale de J sur une surface plane ! perpendiculaire a l'axe de l'arc ! ATTENTION : changer la valeur des tests sur CDGFAC(3,IFAC) ! en fonction du maillage ipcdc3 = ipproc(idjr(3)) elcou = 0.d0 ! do ifac = 1, nfac if( surfac(1,ifac).lt.1.d-8 .and. surfac(1,ifac).gt.-1.d-8 & .and. surfac(2,ifac).lt.1.d-8 .and. surfac(2,ifac).gt.-1.d-8 & .and. cdgfac(3,ifac).lt.5.d-3 .and. cdgfac(3,ifac) .gt.-5.d-3) then iel = ifacel(1,ifac) elcou = elcou + propce(iel,ipcdc3) * surfac(3,ifac) endif enddo !elcou = abs(elcou) if(irangp.ge.0) then call parsom (elcou) endif if ( abs(elcou).ge.1.d-06 ) then elcou=abs(elcou) else elcou=0.d0 endif if(elcou.ne.0.d0) coepoa = couimp/elcou coepot = coepoa WRITE(NFECRA,*) ' ELCOU = ',ELCOU, elcou, dtimes dtj = 1.d15 dtjm =dtj delhsh = 0.d0 cdtj= 2.0d2 do iel = 1, ncel if(propce(iel,ipproc(irom)).ne.0.d0) & delhsh = propce(iel,ipcefj) * dt(iel) & /propce(iel,ipproc(irom)) if(delhsh.ne.0.d0) then dtjm= rtp(iel,isca(iscalt))/delhsh else dtjm= dtj endif dtjm=abs(dtjm) dtj =min(dtj,dtjm) enddo if(irangp.ge.0) then call parmin (dtj) endif cpmx= sqrt(cdtj*dtj) coepot=cpmx if(ntcabs.gt.3) then if(coepoa.ge.1.05d0) then coepot=cpmx else coepot=coepoa endif endif write(nfecra,1008)cpmx,coepoa,coepot write(nfecra,1009)elcou,dpot*coepot ! RECALAGE DES VARIABLES ELECTRIQUES ! ---------------------------------- ! Valeur de DPOT ! -------------- dpot = dpot*coepot ! Potentiel Electrique (on pourrait eviter ; c'est pour le post) ! -------------------- do iel = 1, ncel rtp(iel,isca(ipotr)) = rtp(iel,isca(ipotr))*coepot enddo ! Densite de courant (sert pour A et pour jXB) ! ------------------ if(ippmod(ielarc).ge.1 ) then do idimve = 1, ndimve do iel = 1, ncel ipdcrp = ipproc(idjr(idimve)) propce(iel,ipdcrp) = propce(iel,ipdcrp) * coepot enddo enddo endif ! Effet Joule (sert pour H au pas de temps suivant) ! ----------- ipcefj = ipproc(iefjou) do iel = 1, ncel propce(iel,ipcefj) = propce(iel,ipcefj)*coepot**2 enddo endif !-------- ! FORMATS !-------- 1001 format(/, ' Courant impose= ',E14.5, /, & ' Dpot= ',E14.5,/, & ' Somje= ',E14.5) 1008 format(/,' Cpmx = ',E14.5,/, & ' COEPOA = ',E14.5,/, & ' COEPOT = ',E14.5) 1009 format(/,' Courant calcule = ',E14.5,/, & ' Dpot recale = ',E14.5) !---- ! FIN !---- return end subroutine uselrc