Great. I find that what I suggested is consistent with
http://www.ctcms.nist.gov/fipy/documentation/numerical/discret.html#linear-equations
and
with eq. 1 in "numerical schemes".
On 10 July 2013 22:45, Jonathan Guyer wrote:
>
> On Jul 10, 2013, at 3:14 PM, Noam Yorav-Raphael
> wrote:
>
> > In
On Jul 10, 2013, at 3:14 PM, Noam Yorav-Raphael wrote:
> In the "numerical schemes" page, in the first paragraph:
>
> aA-Ff>0 should be aA+Ff>0
> Pf=-Ff/Df should be Pf=Ff/Df
>
> At least that's what makes sense to me.
I'll let Wheeler address that. I know we had some sign changes in how we
In the "numerical schemes" page, in the first paragraph:
aA-Ff>0 should be aA+Ff>0
Pf=-Ff/Df should be Pf=Ff/Df
At least that's what makes sense to me.
I really tried to file a ticket, with no success. Perhaps, if you want to
encourage contributions, you can consider to migrate to github - the t
Thanks! Now I managed to get the expected result.
It took me some time to realize that since ρ is time-dependent I had to use
hasOld=True and updateOld(). I wonder if it's possible for fipy to detect
that the TransientTerm coefficient has changed and issue an informative
warning.
Thanks again for
On Wed, Jul 10, 2013 at 9:02 AM, Noam Yorav-Raphael wrote:
> I'm sorry, it turns out I was wrong about the equation. The correct
> equation is ∂(e/ρ)/∂t = 1/3 ∂²e/∂m², which is solvable by using
> TransientTerm(coeff=1/rho).
> Thanks!
>
I hope that things work out for you and thanks for the kind
I'm sorry, it turns out I was wrong about the equation. The correct
equation is ∂(e/ρ)/∂t = 1/3 ∂²e/∂m², which is solvable by using
TransientTerm(coeff=1/rho).
Thanks!
On 9 July 2013 17:27, Noam Yorav-Raphael wrote:
> Hello,
>
> I'm trying to simulate radiation diffusion in an expanding matter.