[aspect-devel] Temperature solve

Thu Oct 6 13:53:02 PDT 2016

Thank you Timo!

Best regards,
Lev Karatun.

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

> The problem is that you describe no temperature boundary:
>   set Fixed temperature boundary indicators   = #bottom, top, left, right
>
> I created an issue here https://urldefense.proofpoint.
> com/v2/url?u=https-3A__github.com_geodynamics_aspect_issues_
> 1243&d=CwIFaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=
> c08Btfq4m9QEScXN3ZQwLZzzWQE7S8CYq1IYuzKV_Zk&m=y_
> TFlIfMfInZtERUh97h257MkQg48ek7TZU444VFaRY&s=dwWO6KB1n46LOVu8zqAWt9qqIjX_
> Qk6ScjvzhpfmDAE&e=
>
> On Thu, Oct 6, 2016 at 3:11 PM, Lev Karatun <lev.karatun at gmail.com> wrote:
> > Hi Timo, John,
> >
> > I attached the source and prm files showing the issue.
> > If it's easier, the git branch is here:
> > https://urldefense.proofpoint.com/v2/url?u=https-3A__github.
> com_lkaratun_aspect_tree_hillslope-5Fdiffusion&d=CwIFaQ&c=Ngd-
> ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=c08Btfq4m9QEScXN3ZQwLZzzWQE7S8
> CYq1IYuzKV_Zk&m=y_TFlIfMfInZtERUh97h257MkQg48ek7TZU444VFaRY&s=z_
> rYokhlsBQMbaPfMLnYItZGevDWNQLf60RW9klkjGc&e=
> > 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
> >> 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=y_TFlIfMfInZtERUh97h257MkQg48ek7TZU444VFaRY&s=
> jvcJxvQb9bZXHWkieCPN_c3MDsBjr1PmbURXZrFv-uc&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=uJNnn3tyIUu6qXvrGpmsCu_f4LEBMiHMdXQDNjFA8ss&s=
> OOPvMf1lRNvKLAvagZx_1kvv-oYVn-nq7PEiGtNddWc&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/49a5b620/attachment.html>