Problem source terme in eltssc.f90 (particular physical)

pleblanc86 · Post by **pleblanc86** » Wed May 21, 2014 11:31 am

Hello everybody,

We recently implemented a new mode in the electrical module to simulate the double electrical layer at the interface solid/liquid. But we have the question about writing our source terms in eltssc file.

Before, We did that with CRVEXP and CRVIMP but now we use smbrs, rovsdt. I think best you give an example of our source term.

Considering the term source of $n^L$

$k_r^L n_p^L n_n^L \rho_{mass}^2 - k_d^L n^L \rho_{mass}$

We can determine the explicit term equal to

tsexp = $volume \times k_r^L n_p^L n_n^L \rho_{mass}^2$

We can determine the implicit term equal to

tsimp = $volume \times k_d^L \rho_{mass}$

And so, I'm writing un eltssc.f90 :

$smbrs(iel) = smbrs(iel) + tsexp(iel) + tsimp(iel)*n^L rovdsdt(iel) = rovdsdt(iel) + max(-timp,zero)$

My question is very simple, Does my term source is correct

Thank you for your assistance
Best Regards
P.L.