[aspect-devel] Composition refinement strategy

[John Naliboff]

Hi Lev,

The AMR should indeed be refining according to the compositional gradients. If you want to refine in regions without compositional gradients or another property (e.g. viscosity, velocity, etc), you can use the “minimum refinement function” option to initially refine based on position.

Cheers,
John

*************************************************
John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis






> On May 17, 2017, at 12:20 PM, Lev Karatun <lev.karatun at gmail.com> wrote:
> 
> Hi,
> 
> I was using the composition refinement strategy for my models, and it's not working I expected it to. I have 5 c.f. in my model - upper and lower crust, weak zone, lithospheric and sub-lithospheric mantle. Lithoispheric mantle is shown on screenshot "composition_setup" (attached). The refinement scaling factors are as follows (in order): 
> Compositional field scaling factors = 1,1,2,1,0
> So upper, lower crust and the lithosphere have the same weight. However, when I look at the error indicator (and the resulting mesh), it looks like the refinement plugin prioritizes the boundaries of the c.f. over the interior (so basically it's looking at gradients, not c.f. value itself) -- see screenshot composition_error. 
> 
> Can someone please have a look and tell me if it's just me or if the plugin is indeed not working properly?
> 
> Thanks in advance!
> 
> Best regards,
> Lev Karatun.
> <composition_error.png><composition_setup.png>_______________________________________________
> 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/20170517/653e2173/attachment.html>