Hi, We are solving a system of 5 coupled non-linear PDEs.
As shown in this plot of residuals vs. sweep count https://imperialcollegelondon.box.com/s/9davbq2gq5eani98xuuj2cw9tmz4mbu3 , our residuals die down very slowly, i.e. the solver bottoms out. The drop in all the residuals is linear at first, and then asymptotically bottoms out to a value. How do we get our residuals to drop faster, i.e. with lesser sweeps and faster convergence ? I tried changing solvers and tolerances, but curiously enough the results remain identical. Any pointers on this will be much appreciated. Krishna
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]