[aspect-devel] Temperature solve

[Lev Karatun]

Hi Timo, John,

I attached the source and prm files showing the issue.
If it's easier, the git branch is here:
https://github.com/lkaratun/aspect/tree/hillslope_diffusion
Commit a7a68c1 (the latest one as of right now)

Best regards,
Lev Karatun.

2016-10-06 10:02 GMT-04:00 Timo Heister <heister at clemson.edu>:

> You could be hitting a bug where we are setting the linear solver
> tolerance to a value that is too small. Can you try to create a simple
> .prm that shows the issue?
>
> On Thu, Oct 6, 2016 at 9:33 AM, John Naliboff <jbnaliboff at ucdavis.edu>
> wrote:
> > Hi Lev,
> >
> > Yes, something is definitely off.  What are your temperature BC and what
> > kind non-linear solver are you using?
> >
> > The issues may be completely unrelated, but in other codes I have had
> > situations arise where the temperature solver had trouble when surface
> > erosion laws were in place.  Never figured out the issue as it only
> occurred
> > for a very small subset of models.
> >
> >
> > Cheers,
> > John
> >
> > *************************************************
> > John Naliboff
> > Assistant Project Scientist, CIG
> > Earth & Planetary Sciences Dept., UC Davis
> >
> >
> >
> > On Oct 6, 2016, at 12:56 AM, Lev Karatun <lev.karatun at gmail.com> wrote:
> >
> > Hi,
> >
> > I was testing one of the models for the hillslope diffusion algorithm
> with
> > constant temperature (1) througout the entire domain, but at some point
> > during the simulation (about when the surface gets flat and velocities
> are
> > basically zero) it starts taking increasingly bigger number of iterations
> > for temperature solver to converge, failing eventually. I was wondering
> if
> > this is the correct behaviour? Isn't the temerature solver supposed to
> > converge right away if there is no temperature gradient? Or am I missing
> > something?
> >
> > I'm not sure if it's relevant, but when I vizualized temperature in
> > paraview, the scale ranged from 1 to 2, although all the points had
> > temperature 1 (none with 2).
> >
> > Best regards,
> > Lev Karatun.
> > _______________________________________________
> > Aspect-devel mailing list
> > Aspect-devel at geodynamics.org
> > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.
> geodynamics.org_cgi-2Dbin_mailman_listinfo_aspect-2Ddevel&d=CwIFaQ&c=Ngd-
> ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=c08Btfq4m9QEScXN3ZQwLZzzWQE7S8
> CYq1IYuzKV_Zk&m=HNLCa6l0Ml2HYYJsrhmwxHL1YqE4D0LeVZyUjalr6BM&s=
> gXWLQmtAIoa8olhfUTpFolIRLj3c2caoDQE8krNVoKM&e=
> >
> >
> >
> > _______________________________________________
> > Aspect-devel mailing list
> > Aspect-devel at geodynamics.org
> > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.
> geodynamics.org_cgi-2Dbin_mailman_listinfo_aspect-2Ddevel&d=CwIGaQ&c=Ngd-
> ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=R5lvg9JC99XvuTgScgbY_
> QFS80R7PEA2q0EPwDy7VQw&m=piOX_59vSmVFG7de-1ArbqczU5LnD6246cYs-S0qXZk&s=
> BMkFdhy5uqtk_S_jX60s2X5XdXcyfw39FsKPObganx8&e=
>
>
>
> --
> Timo Heister
> http://www.math.clemson.edu/~heister/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161006/91cacb8a/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: free_surface.h
Type: text/x-chdr
Size: 8108 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161006/91cacb8a/attachment-0001.h>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: free_surface.cc
Type: application/octet-stream
Size: 50014 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161006/91cacb8a/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 11.prm
Type: application/octet-stream
Size: 6614 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161006/91cacb8a/attachment-0003.obj>