Yes, now I see this. So, maybe mathematically I have setup the problem
incorrectly. I would like to allow the concentrations at the interface to
increase if the rate of back diffusion becomes greater than the rate of decay.
The rate of decay at the interface is dc/dt = -kt.
________________________________
From: Jonathan Guyer <[email protected]>
To: FIPY list <[email protected]>
Sent: Monday, May 13, 2013 2:16 PM
Subject: Re: Transient Diffusion with Time Varying BC
On May 13, 2013, at 3:58 PM, Chuck Holbert <[email protected]> wrote:
> To estimate back diffusion and an accumulation at the interface (left-most
> boundary), don't I want the precursor of the exponential to change at every
> timestep as phi changes?
>
> Forgive me if this should be evident, but want to make sure I understand.
I don't know that it's evident either way, but I don't think you want the
precursor changing as a function of time other than when the decay rate
changes. If nothing else, I note that you don't change the precursor during the
"INITIAL PERIOD FROM 1970 TO 2013", it's just 250. for the entire time period.
In subsequent periods, the reservoir to the left has already decayed away to
some lower value and you want to make it now decay at a lower rate; if you have
the precursor changing at every timestep, I think you're double counting the
decay.
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
