[aspect-devel] Temperature solve

[Timo Heister]

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