Hi Everyone,

I  remind you what my problem  is about,  there is a diffusing field c, a
growing interface (defined by the level set parameter phi) and I wanted an
absorbing boundary at the interface,

With the valuable help of Daniel and Jonathan, I modified this file
(MetalIonSourcevariable.py) like this

To do this replace line 88 in
<http://matforge.org/fipy/browser/trunk/fipy/models/levelSet/electroChem/metalIonSourceVariable.py#L88>
with "return self.distanceVar._cellInterfaceFlag * 1e+20".  That will
make c=0 on the interface
(phi=0).


I want to understand this trick, because I am now interested in putting an
appropriate source term in this problem.
If I understand correctly the return value of this function is going to be
used later on, in a source term for the diffusing field. So why is it an
infinitely large value, if the boundary is absorbing ? and not a -infinity
term ?

also if I want to add a positive source term, creation of a diffusive
species at the interface of a given finite value c1, how can I combinate
that with this infinity trick ??

thanks for your help !!

Julien

Reply via email to