[aspect-devel] Temperature solve

[Timo Heister]

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=c08Btfq4m9QEScXN3ZQwLZzzWQE7S8CYq1IYuzKV_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/