[aspect-devel] Aspect hangs after several time steps

[Timo Heister]

> Did you activate the compositional fields Timo?

No, I just ran the .prm he sent.

On Tue, Feb 23, 2016 at 6:18 PM, Wolfgang Bangerth <bangerth at tamu.edu> wrote:
> On 02/23/2016 04:58 PM, Rene Gassmoeller wrote:
>>
>>
>> The problem is that the 'Simpler' material model is not using
>> compositional
>> fields, and does not consider that the model might still be solved with
>> passive compositional fields. Therefore the 'reaction_terms' variable,
>> which
>> describes the change of the compositional fields over one timestep remains
>> unfilled and contains the uninitialized value 'NaN'.
>
>
> Ah, wow, good detective work. I think that the material models should
> actually set these reaction terms to zero, since that's the material
> behavior they describe. This is likely an issue for not just the Simpler
> material model.
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel



-- 
Timo Heister
http://www.math.clemson.edu/~heister/