does not satisfy the zero-flux boundary
>>>>>>>>> condition on species 3. In the absence of a potential gradient, that
>>>>>>>>> species would need to also have zero gradient (like species 2 in the
>>>>>>>>> curr
>>>>> use a boundary condition as that would be more applicable in higher
>>>>>>>> dimensions.
>>>>>>>>
>>>>>>>> ☘
>>>>>>>>
>>>>>>>>
>>>>>>>&
t;>>> Scott -
>>>>>>>>
>>>>>>>> Eq2 appears to have a typo:
>>>>>>>>
>>>>>>>> -Eq2 = DiffusionTerm(coeff=z2*C1, var=Phi) + DiffusionTerm(coeff=1.0,
>>>>>>>> var
t;
>>>>>>> -Eq2 = DiffusionTerm(coeff=z2*C1, var=Phi) + DiffusionTerm(coeff=1.0,
>>>>>>> var=C2)
>>>>>>> +Eq2 = DiffusionTerm(coeff=z2*C2, var=Phi) + DiffusionTerm(coeff=1.0,
>>>>>>> var=C2)
>>>
ffusionTerm(coeff=z3*C3, var=Phi) + (C3.faceGrad).divergence
>>>>>> +Eq3 = DiffusionTerm(coeff=z3*C3, var=Phi) - DiffusionTerm(coeff=-z1/z3,
>>>>>> var=C1) - DiffusionTerm(coeff=-z2/z3, var=C2)
>>>>>>
>>>>>> [Note that
ht and it's
>>>>> unusual to care.]
>>>>>
>>>>> With these changes, the solution converges in one sweep, but does not
>>>>> agree with the Analytical Solution you give from West.
>>>>>
>>>>> I don't u
ou don't need Eq3 (it'll get taken care of automatically) and
> you need some form of Poisson's equation.
>
> - Jon
>
>> On Jun 11, 2019, at 4:48 PM, Scott Calabrese Barton wrote:
>>
>> I’ve had some trouble using fipy for the 1D, steady-state ternary
>> el
or the 1D, steady-state ternary electrolyte
> problem, which involves coupled, nonlinear migration and diffusion.
>
> A text version of the model is listed below. A more complete development is
> at
> https://gcc01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fbit.
I’ve had some trouble using fipy for the 1D, steady-state ternary electrolyte
problem, which involves coupled, nonlinear migration and diffusion.
A text version of the model is listed below. A more complete development is at
http://bit.ly/2wJDXob
from fipy import *
# Parameters
nx=50
C20=4.0