[aspect-devel] Temperature solve

[John Naliboff]

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
> 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/43c845b1/attachment.html>