[aspect-devel] Checkerboard at large scale

[Lev Karatun]

Hi everyone,

just wanted to let you know, it seems I did have a slight imbalance in the
boundary velocity, and after fixing it (and making net volume flux equal to
zero) the solver successfully converged with constant mesh. It still chokes
when I enable AMR, so I'll keep the free surface enabled as you suggested.

Really appreciate the help!

Best regards,
Lev Karatun.

2018-04-08 0:34 GMT-04:00 Wolfgang Bangerth <bangerth at colostate.edu>:

> On 04/07/2018 02:09 PM, Juliane Dannberg wrote:
>
>>
>> If that is not the case, and all boundaries are closed, you will not get
>> a reasonable solution (or the model will just not converge). John,
>> Wolfgang, please correct me if you think that works differently.
>>
>
> Correct. For incompressible models, there must be a zero net in/outflux or
> either the solution you get is likely wrong, or the solver will simply not
> converge to high accuracy.
>
>
> If, on the other hand, your model is _compressible_, this would work
>> differently (so it would be good if you could tell us if it is).
>>
>
> Indeed. If the model is compressible, then a net influx will simply
> increase the pressure to the point where the density increases (due to the
> higher pressure) to accommodate the net influx.
>
> Best
>  W.
>
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth          email:                 bangerth at colostate.edu
>                            www: http://www.math.colostate.edu/~bangerth/
>
> _______________________________________________
> 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/20180411/76544288/attachment.html>