Electrical DC arc simulation problem

[husanhao]

Dear all,

I would like to know if I want to model the DC electrical arc between two vertical electrodes in Code Saturne 3.0.0, do I need to consider the subroutine uselrc.f90? Please refer to the picture to see the detail about my model.

When I try to use uselrc.f90 subroutine in my model, I always got the error "error (divergence) solving for Enthalpy", so I think probability I did something wrong in this subroutine.

The only I changed in uselrc.f90 is the current and the calculation of the intensity of the arc current. I changed this based on my geometry, I set these two parts as below:

a. Current (I believe that this part should be ok)

if ( ntcabs.le.200 ) then
couimp = 5.d0 + (100.d0-5.d0) * ntcabs/200
else
couimp = 100.d0
endif

b. Calculation of the intensity of the arc current

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

I have attached my subroutine 'uselrc.f90'.
Can anyone please help me to figure out this problem, thanks for your help, I really appreciate it !!!

Best Regards
Eric