[aspect-devel] mesh refinement strategy question

Menno Fraters menno.fraters at outlook.com
Mon Jun 11 03:40:17 PDT 2018


Hey Magali,


Just as a note, the mail was send between the mail of Felipe and Rene, but delivered very late because of the large image required approval of someone.


What I forgot to mention in my mail is that you set a minimum an maximum refinement level in the plugin, so it can be easily combined with the gradient approach if needed. For me there where some other reasons to go with this approach (such as the mesh dependence of plasticity and the unpredictability of the relative errors during a changing simulation), but like you said there are many ways to tackle this problem :)


Cheers,


Menno

________________________________
From: Aspect-devel <aspect-devel-bounces at geodynamics.org> on behalf of Magali Billen <mibillen at ucdavis.edu>
Sent: Monday, June 11, 2018 9:16:37 AM
To: aspect-devel at geodynamics.org
Subject: Re: [aspect-devel] mesh refinement strategy question

Hello Menno,
Than you for this suggestion also and sharing the plug-in.  Once the simulation gets going, the mesh refinement
will depend on multiple parameters (temperature, strain-rate, composition) and so having it depend on the gradients of these parameters is, I think, the right thing to do. But, it is helpful to have other options to think about if this doesn’t work out the way I anticipate.

-Magali

On Jun 7, 2018, at 11:28 PM, Menno Fraters <menno.fraters at outlook.com<mailto:menno.fraters at outlook.com>> wrote:

Hey Magali,

This is a problem I also struggled quite a bit with during my master thesis, and in the end I arrived at the solution which works really well for my applications (both in 2d and 3d) and I think it will also well for your application since it looks very similar.

The problem here is that most mesh refinement strategies in aspect are gradient based with an error indicator or location based. The trick here is to realise that you probably only want to refine in cold area's, because that is where the slab currently is. So what you want to do is refined based on absolute temperature or isotherms. I have written a plugin for this quite a while ago, but never got around to fully merging it into the main branch. You can find it at: https://github.com/MFraters/aspect/tree/MFraters-isoterms-plugin, or in the pull requesthttps://github.com/geodynamics/aspect/pull/295. I have attached a figure which shows how the grid follows the evolution of the slab really nicely based on this approach.

If you have questions on this, or need help with it, please let me know.

Cheers,

Menno
________________________________
From: Aspect-devel <aspect-devel-bounces at geodynamics.org<mailto:aspect-devel-bounces at geodynamics.org>> on behalf of FELIPE ORELLANA ROVIROSA <f_orellana at berkeley.edu<mailto:f_orellana at berkeley.edu>>
Sent: Thursday, June 7, 2018 9:20:00 PM
To: aspect-devel at geodynamics.org<mailto:aspect-devel at geodynamics.org>
Subject: Re: [aspect-devel] mesh refinement strategy question


 Hi Magali,

      Nice project.

      I am no expert on this, but I have some experience. I hope other modelers can speak out.

     Time steps between mesh refinement = 1  !!! watch out with this, it may be appropriate for some running time, but not always.

     For the strategy, if you use Temperature, that should help the appropriate solving of the energy equation. More so, if your slab is defined as  thermal lithosphere, a proper temp-based refinement should reveal the slab properly, with lots of mesh points defining its shape.

    The sole temp strategy may fail to properly solve the momentum equation, for which you should turn to the strain rate strategy.

    These days I am using:

       set Strategy = strain rate, temperature
       set Refinement criteria merge operation = plus

but in my experience it 's kind of an art..sometimes I do find problems, even with all of these precautions, and I need to play with refinement/coarsening fractions.

   I guess you should start with a purely temp-based refinement, maybe at the very start, or maybe even for some run-geologic time. Then, I would surely add an extra criteria, something that takes care of the velocity field (e.g. the strain rate strategy).

cheers,
Felipe


On Thu, Jun 7, 2018 at 8:48 AM, Magali Billen <mibillen at ucdavis.edu<mailto:mibillen at ucdavis.edu>> wrote:
Hello,

I am starting a project for subduction modeling in Aspect and I have a  question about how certain mesh refinement strategies work together. Specifically, I need to define a “proto-slab” based on temperature,
 which is needed to start the subduction process. However, once the slab starts sinking its location in the model domain depends on the dynamics and it can move 1000s of kilometers from its starting location.

So, I have first set this up using the  Strategy = minimum refinement function (see below for the
full subsection definition). And, I can get a nice looking initial slab (see below).  However, after the initialization,
I would like the mesh refinement to depend on things like the temperature, (and later composition, and strain-rate).  It does not work to just use Strategy = temperature, because this just refines on the temperature field defined on the coarsest mesh (see picture).

So, what is not clear to me is whether there is a way to use one Strategy during the initial adaptive refinement
and a different (set) of criteria after that point, or to combine them such that the minimum refinement function
is no longer have much impact after the first step.

Has anyone else had a problem like this?  Any ideas?
Instead of thinking about refining, should I start with a really fine mesh everywhere and then use the
Strategy = temperature to coarsen the mesh where refinement isn’t needed. It doesn’t seem like it set up
to work that way?

- Magali
Here’s the prm file for the nice looking slab using the function to refine based on location:


Here’s just the mesh  refinement section:

subsection Mesh refinement
  set Initial global refinement                = 4
  set Minimum refinement level                 = 4
  set Initial adaptive refinement              = 5
  set Time steps between mesh refinement       = 1
  set Strategy = minimum refinement function
  subsection Minimum refinement function
    set Coordinate system   = cartesian
    set Variable names      = x,y
    set Function constants = YMAX=1e6,LITH=1.5e5,SM=2.5e5,XTR=2.2e6,RC=4e5
    set Function expression = if( (YMAX-y<=LITH)||((YMAX-y<=SM)&&(x<=XTR+3*RC)&&(x>=XTR-0.5*RC)),9,0)
  end
end

And, here are some screen-shots

Full domain




Zoom-in on slab



This is what happens if I just use Strategy = temperature


And if I combine them  (this is just the 0th time-step), I need to set up some other things before
I can run forward in time.


____________________________________________________________
Professor of Geophysics
Earth & Planetary Sciences Dept., UC Davis
Davis, CA 95616
2129 Earth & Physical Sciences Bldg.
Office Phone: (530) 752-4169
http://magalibillen.faculty.ucdavis.edu<http://magalibillen.faculty.ucdavis.edu/>

Currently on Sabbatical at Munich University (LMU)
Department of Geophysics (PST + 9 hr)

Avoid implicit bias - check before you submit:
http://www.tomforth.co.uk/genderbias/
___________________________________________________________


_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

<ma12-m4-25nli-cfls-tstp100-grid.png>_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

____________________________________________________________
Professor of Geophysics
Earth & Planetary Sciences Dept., UC Davis
Davis, CA 95616
2129 Earth & Physical Sciences Bldg.
Office Phone: (530) 752-4169
http://magalibillen.faculty.ucdavis.edu

Currently on Sabbatical at Munich University (LMU)
Department of Geophysics (PST + 9 hr)

Avoid implicit bias - check before you submit:
http://www.tomforth.co.uk/genderbias/
___________________________________________________________

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180611/c7b28083/attachment-0001.html>


More information about the Aspect-devel mailing list