Dear FiPy developers,
I am using FiPy to solve a coupled system of two partial
differential equations governing the
evolution of two surface densities, including convection,
diffusion and source terms. To check if
everything was going well, I checked if mass
was conserved in the system by checking if
the inflowing mass equals the outflowing mass.
This was only the case if I set the desired residual low enough in:
res = 1e10
nres = 0
while res>1e-2:
res = Eq.sweep(dt=timestepduration)
nres += 1
But I have no clue how to determine the desired residual other than
by just checking if mass is
conserved in the solution.
Could anybody explain to me what the residuals actually mean?
If one is working with very large numbers in FiPy, can the desired
residuals then also be larger
than if one is working with small numbers,
and get the same accuracy?
Thanks a lot in advance!
Best regards,
Djoeke
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
