[aspect-devel] Questions about surface process
Lev Karatun
lev.karatun at gmail.com
Sat Dec 10 10:53:00 PST 2016
Hi Yimin,
surface processes is still work in progress. I'm not modifying any of the
mesh refinement related functions, only solving the 1- or 2-D diffusion
equation on the surface nodes displacement and adjusting their velocities
accordingly.
Ian Rose would probably be the best person to ask the questions about
partitioning since he's the one who implemented free surface in Aspect.
Best regards,
Lev Karatun.
2016-12-10 0:03 GMT-05:00 Yimin Jin <choby1234567 at gmail.com>:
> Sorry for the unedited subject last time. Shall I send this mail to this
> address or directly to the person that answered me?
>
> Contents of the previous mail I sent:
>
> Hi John,
>
> Thanks for answering. It is inspiring to know that surface erosion models
> have been implemented in ASPECT, but I still have some doubts: 1. The
> adaptive refinements make the partition changes as the timestep goes on, so
> if we use the function GridGenerator::extract_boundary_mesh() to extract
> the surface mesh as an object of
> parallel::distributed::Triangulation<dim-1,dim>, how can we keep the
> partition of the surface mesh the same as the volume mesh? Or shall we
> partition the two meshes independently and transport velocity data between
> processes? 2. By saying 'I’m not sure what the status is of further
> development pushing it to the main branch', did you mean that it is
> difficult to merge this module to the program structure, or the surface
> process module is not suitable to be added to ASPECT?
>
> Cincerely,
> Yimin Jin
>
> 2016-12-10 11:14 GMT+08:00 <aspect-devel-request at geodynamics.org>:
>
>> Send Aspect-devel mailing list submissions to
>> aspect-devel at geodynamics.org
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect
>> -devel
>> or, via email, send a message with subject or body 'help' to
>> aspect-devel-request at geodynamics.org
>>
>> You can reach the person managing the list at
>> aspect-devel-owner at geodynamics.org
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of Aspect-devel digest..."
>>
>>
>> Today's Topics:
>>
>> 1. Re: Van Keken et al. 2008 subduction benchmark (Wolfgang Bangerth)
>> 2. Re: Aspect-devel Digest, Vol 61, Issue 4 (Yimin Jin)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Fri, 9 Dec 2016 14:42:06 -0700
>> From: Wolfgang Bangerth <bangerth at tamu.edu>
>> To: maxwellr at gmail.com, aspect-devel at geodynamics.org, Scott King
>> <sdk at vt.edu>
>> Subject: Re: [aspect-devel] Van Keken et al. 2008 subduction benchmark
>> Message-ID: <e5858bf2-4b77-cafd-087c-d500324a45a9 at tamu.edu>
>> Content-Type: text/plain; charset=utf-8; format=flowed
>>
>> On 12/09/2016 12:12 PM, Max Rudolph wrote:
>> >
>> > I did not try the modification to allow a linear ramp. The problem that
>> I
>> > encountered was not associated with the pressure field in the wedge. It
>> was
>> > actually along the subduction interface where the downgoing plate meets
>> the
>> > over-riding plate. There is a discontinuity in the velocity field
>> there. The
>> > velocity field is locally divergence-free but it may not be
>> divergence-free in
>> > the integral sense when interpolated using the finite element basis
>> functions.
>>
>> That's one possible issue. When you say that it's locally divergence free,
>> does this also hold at the points along the interface?
>>
>>
>> > One possible work-around would be to prescribe zero dynamic pressure
>> > throughout the rigid downgoing and over-riding plate (the correct
>> solution).
>> > Doing so would mean prescribing all of the degrees of freedom for every
>> node
>> > in the prescribed velocity region. I think that this would be a very
>> > straightforward modification of the prescribed velocity plugin. Maybe
>> we can
>> > chat for a few minutes about this at the ASPECT dinner?
>>
>> Sure. Yours sounds like a reasonable solution. I don't recall who wrote
>> this
>> plugin, so I don't know how difficult it would be to implement what you
>> suggest; but it can't be very difficult at all.
>>
>> Best
>> W.
>>
>>
>> --
>> ------------------------------------------------------------------------
>> Wolfgang Bangerth email: bangerth at colostate.edu
>> www: http://www.math.colostate.edu/~bangerth/
>>
>>
>>
>> ------------------------------
>>
>> Message: 2
>> Date: Sat, 10 Dec 2016 10:53:42 +0800
>> From: Yimin Jin <choby1234567 at gmail.com>
>> To: aspect-devel at geodynamics.org
>> Subject: Re: [aspect-devel] Aspect-devel Digest, Vol 61, Issue 4
>> Message-ID:
>> <CALhZPeCxzDFjkH2+478GDZ+ou=nphObb3iY7D4-dzkEm_01HRQ at mail.gm
>> ail.com>
>> Content-Type: text/plain; charset="utf-8"
>>
>> Hi John,
>>
>> Thanks for answering. It is inspiring to know that surface erosion models
>> have been implemented in ASPECT, but I still have some doubts: 1. The
>> adaptive refinements make the partition changes as the timestep goes on,
>> so
>> if we use the function GridGenerator::extract_boundary_mesh() to extract
>> the surface mesh as an object of
>> parallel::distributed::Triangulation<dim-1,dim>, how can we keep the
>> partition of the surface mesh the same as the volume mesh? Or shall we
>> partition the two meshes independently and transport velocity data between
>> processes? 2. By saying 'I’m not sure what the status is of further
>> development pushing it to the main branch', did you mean that it is
>> difficult to merge this module to the program structure, or the surface
>> process module is not suitable to be added to ASPECT?
>>
>> Cincerely,
>> Yimin Jin
>>
>> 2016-12-10 5:39 GMT+08:00 <aspect-devel-request at geodynamics.org>:
>>
>> > Send Aspect-devel mailing list submissions to
>> > aspect-devel at geodynamics.org
>> >
>> > To subscribe or unsubscribe via the World Wide Web, visit
>> > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspec
>> t-devel
>> > or, via email, send a message with subject or body 'help' to
>> > aspect-devel-request at geodynamics.org
>> >
>> > You can reach the person managing the list at
>> > aspect-devel-owner at geodynamics.org
>> >
>> > When replying, please edit your Subject line so it is more specific
>> > than "Re: Contents of Aspect-devel digest..."
>> >
>> >
>> > Today's Topics:
>> >
>> > 1. Re: Questions about surface process (John Naliboff)
>> > 2. Re: Van Keken et al. 2008 subduction benchmark (Max Rudolph)
>> > 3. Re: Van Keken et al. 2008 subduction benchmark
>> > (Thieulot, C.A.P. (Cedric))
>> >
>> >
>> > ----------------------------------------------------------------------
>> >
>> > Message: 1
>> > Date: Fri, 9 Dec 2016 11:03:02 -0800
>> > From: John Naliboff <jbnaliboff at ucdavis.edu>
>> > To: aspect-devel at geodynamics.org
>> > Subject: Re: [aspect-devel] Questions about surface process
>> > Message-ID: <822D9590-B1B1-4FDB-96E4-E2BDF4B862E0 at ucdavis.edu>
>> > Content-Type: text/plain; charset="utf-8"
>> >
>> > Hi Yimin,
>> >
>> > I believe a number of people implemented a preliminary" diffusive"
>> surface
>> > erosion model this summer. However, I’m not sure what the status is of
>> > further development pushing it to the main branch. Hopefully someone
>> can
>> > chime in on this.
>> >
>> > Cheers,
>> > John
>> >
>> > *************************************************
>> > John Naliboff
>> > Assistant Project Scientist, CIG
>> > Earth & Planetary Sciences Dept., UC Davis
>> >
>> >
>> >
>> >
>> >
>> >
>> > > On Dec 9, 2016, at 2:30 AM, Yimin Jin <choby1234567 at gmail.com> wrote:
>> > >
>> > > Dear developers of ASPECT:
>> > >
>> > > I am a student of University of Chinese Academy of Sciences and I have
>> > been using ASPECT to do simulations on crack propagation of small-scaled
>> > geological bodies in the passing year. I have found that the existing
>> > velocity pojection strategies in the free-surface module introduces too
>> > much error to the result (the 'normal' strategy leads to severe
>> distortion
>> > of cells, and the 'vertical' strategy ignores horizontal displacements),
>> > while the sticky-air method will double the computational complexity of
>> the
>> > problem (my model can 'grow' tall so the air layer must be thick
>> enough).
>> > Therefore I intended to apply a diffusive surface process on the
>> > free-surface, but when I looked into this problem I found it quite
>> > difficult because the partition changes each time we execute refinement
>> and
>> > coarsening. The finest way seems to be transporting faces between
>> processes
>> > using functions provided by p4est (I don't know if there are any), then
>> it
>> > will be unavoidable to read the codes of p4est. Could you please tell me
>> > whether you have started to do such things or not? If it has been
>> > implemented and will soon be released, then I will be happily looking
>> > forward to the next version of ASPECT; otherwise, I would like to make
>> my
>> > effort in this part of work.
>> > >
>> > > Cincerely,
>> > > Yimin Jin
>> > >
>> > > Key Laboratory of Computational Geodynamics, UCAS
>> > > _______________________________________________
>> > > 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/20161209/d7d9703e/attachment-0001.html>
>> >
>> > ------------------------------
>> >
>> > Message: 2
>> > Date: Fri, 9 Dec 2016 11:12:13 -0800
>> > From: Max Rudolph <maxwellr at gmail.com>
>> > To: Scott King <sdk at vt.edu>
>> > Cc: aspect-devel <aspect-devel at geodynamics.org>
>> > Subject: Re: [aspect-devel] Van Keken et al. 2008 subduction benchmark
>> > Message-ID:
>> > <CAFNmARsEzxYJb6LbOW6s4EEsntcpCJFVgcCASsqCzBk9r8D5xw at mail.
>> > gmail.com>
>> > Content-Type: text/plain; charset="utf-8"
>> >
>> > Scott and Wolfgang,
>> > I did not try the modification to allow a linear ramp. The problem that
>> I
>> > encountered was not associated with the pressure field in the wedge. It
>> was
>> > actually along the subduction interface where the downgoing plate meets
>> the
>> > over-riding plate. There is a discontinuity in the velocity field there.
>> > The velocity field is locally divergence-free but it may not be
>> > divergence-free in the integral sense when interpolated using the finite
>> > element basis functions. One possible work-around would be to prescribe
>> > zero dynamic pressure throughout the rigid downgoing and over-riding
>> plate
>> > (the correct solution). Doing so would mean prescribing all of the
>> degrees
>> > of freedom for every node in the prescribed velocity region. I think
>> that
>> > this would be a very straightforward modification of the prescribed
>> > velocity plugin. Maybe we can chat for a few minutes about this at the
>> > ASPECT dinner?
>> >
>> > Max
>> >
>> > On Fri, Dec 9, 2016 at 1:22 AM, Scott King <sdk at vt.edu> wrote:
>> >
>> > > Max;
>> > >
>> > > Do you use a linear ramp something from the point where the slab meets
>> > the
>> > > over-riding plate down into the wedge? This *smooths* the
>> discontinuity
>> > > (at least as far as the wedge sees it). With this the pressure
>> fields in
>> > > the wedge seems to become more smooth for many codes. I always had
>> to
>> > you
>> > > a much larger "ramp" (25,-50 km) because my elements were not refined
>> as
>> > > much as Herr. Dr. Prof. van Keken. This puts the discontinuity into
>> the
>> > > part of the domain where the velocities are entirely prescribed. It
>> > seems
>> > > to be ok for some codes because they don't "solve" any equation when
>> all
>> > > the velocities are prescribed. Years of arguing and twiddling.
>> > >
>> > > Scott
>> > >
>> > > On Dec 8, 2016, at 10:25 PM, Wolfgang Bangerth <bangerth at tamu.edu>
>> > wrote:
>> > >
>> > > > On 10/27/2016 10:13 AM, Max Rudolph wrote:
>> > > >> I was wondering if anyone has been successful in setting up the Van
>> > > Keken et
>> > > >> al. 2008 subduction benchmark using ASPECT? I have been using this
>> as
>> > a
>> > > >> starting point to set up a more complicated corner flow model and
>> have
>> > > not yet
>> > > >> been successful.
>> > > >>
>> > > >> This benchmark solution requires a rigid overriding plate and
>> > kinematic
>> > > >> subducting plate. There is a discontinuity along the subduction
>> > > interface. I
>> > > >> tried imposing kinematic boundary conditions on the left and
>> bottom of
>> > > the box
>> > > >> and using Jonathan Perry-Houts' prescribed velocity plugin
>> > > >> (cookbooks/prescribed_velocity), but the Stokes solver fails to
>> > > converge.
>> > > >> Looking at the pressure field, it's not surprising to see very
>> large
>> > > >> oscillations in this region where the velocity field is
>> discontinuous.
>> > > It
>> > > >> occurred to me that using the locally conservative discretization
>> > might
>> > > help.
>> > > >> However, when the locally conservative discretization is enabled,
>> the
>> > > internal
>> > > >> prescribed velocities are no longer enforced anywhere. Is it
>> obvious
>> > > why this
>> > > >> might be the case?
>> > > >
>> > > > Max -- I haven't tried implementing the benchmark in ASPECT, but
>> recall
>> > > trying to do so in an early version of step-32 when I visited Scott in
>> > > Blacksburg for the very first time many many years ago. I remember it
>> to
>> > be
>> > > very difficult because the benchmark definition really doesn't fit the
>> > > finite element method.
>> > > >
>> > > > One of the issues with the use of the prescribed velocity is that
>> the
>> > > prescribed velocity needs to be divergence free on each cell. Is it,
>> in
>> > > your case? Even on cells where you only prescribe values at some
>> nodes?
>> > > >
>> > > > 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/20161209/294bcf0f/attachment-0001.html>
>> >
>> > ------------------------------
>> >
>> > Message: 3
>> > Date: Fri, 9 Dec 2016 21:32:48 +0000
>> > From: "Thieulot, C.A.P. (Cedric)" <c.thieulot at uu.nl>
>> > To: "maxwellr at gmail.com" <maxwellr at gmail.com>,
>> > "aspect-devel at geodynamics.org" <aspect-devel at geodynamics.org>
>> > Subject: Re: [aspect-devel] Van Keken et al. 2008 subduction benchmark
>> > Message-ID: <FB2A0542-6101-4613-B2BA-AB5BB4EED818 at uu.nl>
>> > Content-Type: text/plain; charset="windows-1252"
>> >
>> > hi,
>> > I believe that Scott’s advice still holds: if you look at subsection 2.3
>> > of the paper, you’ll see that along the interface
>> > between downgoing plate and overriding place there is no velocity
>> > discontinuity.
>> > Ce/
>> >
>> >
>> >
>> >
>> > [cid:83FFAEC2-5940-4970-B26E-88C7176B394B at lan]
>> >
>> > On 09 Dec 2016, at 20:12, Max Rudolph <maxwellr at gmail.com<mailto:max
>> > wellr at gmail.com>> wrote:
>> >
>> > Scott and Wolfgang,
>> > I did not try the modification to allow a linear ramp. The problem that
>> I
>> > encountered was not associated with the pressure field in the wedge. It
>> was
>> > actually along the subduction interface where the downgoing plate meets
>> the
>> > over-riding plate. There is a discontinuity in the velocity field there.
>> > The velocity field is locally divergence-free but it may not be
>> > divergence-free in the integral sense when interpolated using the finite
>> > element basis functions. One possible work-around would be to prescribe
>> > zero dynamic pressure throughout the rigid downgoing and over-riding
>> plate
>> > (the correct solution). Doing so would mean prescribing all of the
>> degrees
>> > of freedom for every node in the prescribed velocity region. I think
>> that
>> > this would be a very straightforward modification of the prescribed
>> > velocity plugin. Maybe we can chat for a few minutes about this at the
>> > ASPECT dinner?
>> >
>> > Max
>> >
>> > On Fri, Dec 9, 2016 at 1:22 AM, Scott King <sdk at vt.edu<mailto:
>> sdk at vt.edu>>
>> > wrote:
>> > Max;
>> >
>> > Do you use a linear ramp something from the point where the slab meets
>> the
>> > over-riding plate down into the wedge? This *smooths* the
>> discontinuity
>> > (at least as far as the wedge sees it). With this the pressure fields
>> in
>> > the wedge seems to become more smooth for many codes. I always had to
>> you
>> > a much larger "ramp" (25,-50 km) because my elements were not refined as
>> > much as Herr. Dr. Prof. van Keken. This puts the discontinuity into
>> the
>> > part of the domain where the velocities are entirely prescribed. It
>> seems
>> > to be ok for some codes because they don't "solve" any equation when all
>> > the velocities are prescribed. Years of arguing and twiddling.
>> >
>> > Scott
>> >
>> > On Dec 8, 2016, at 10:25 PM, Wolfgang Bangerth <bangerth at tamu.edu
>> <mailto:
>> > bangerth at tamu.edu>> wrote:
>> >
>> > > On 10/27/2016 10:13 AM, Max Rudolph wrote:
>> > >> I was wondering if anyone has been successful in setting up the Van
>> > Keken et
>> > >> al. 2008 subduction benchmark using ASPECT? I have been using this
>> as a
>> > >> starting point to set up a more complicated corner flow model and
>> have
>> > not yet
>> > >> been successful.
>> > >>
>> > >> This benchmark solution requires a rigid overriding plate and
>> kinematic
>> > >> subducting plate. There is a discontinuity along the subduction
>> > interface. I
>> > >> tried imposing kinematic boundary conditions on the left and bottom
>> of
>> > the box
>> > >> and using Jonathan Perry-Houts' prescribed velocity plugin
>> > >> (cookbooks/prescribed_velocity), but the Stokes solver fails to
>> > converge.
>> > >> Looking at the pressure field, it's not surprising to see very large
>> > >> oscillations in this region where the velocity field is
>> discontinuous.
>> > It
>> > >> occurred to me that using the locally conservative discretization
>> might
>> > help.
>> > >> However, when the locally conservative discretization is enabled, the
>> > internal
>> > >> prescribed velocities are no longer enforced anywhere. Is it obvious
>> > why this
>> > >> might be the case?
>> > >
>> > > Max -- I haven't tried implementing the benchmark in ASPECT, but
>> recall
>> > trying to do so in an early version of step-32 when I visited Scott in
>> > Blacksburg for the very first time many many years ago. I remember it
>> to be
>> > very difficult because the benchmark definition really doesn't fit the
>> > finite element method.
>> > >
>> > > One of the issues with the use of the prescribed velocity is that the
>> > prescribed velocity needs to be divergence free on each cell. Is it, in
>> > your case? Even on cells where you only prescribe values at some nodes?
>> > >
>> > > Best
>> > > W.
>> > >
>> > > --
>> > > ------------------------------------------------------------
>> ------------
>> > > Wolfgang Bangerth email:
>> bangerth at colostate.edu
>> > <mailto:bangerth at colostate.edu>
>> > > www: http://www.math.colostate.edu/
>> ~bangerth/
>> > >
>> > > _______________________________________________
>> > > 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
>> >
>> >
>> > _______________________________________________
>> > 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
>> >
>> > -------------- next part --------------
>> > An HTML attachment was scrubbed...
>> > URL: <http://lists.geodynamics.org/pipermail/aspect-devel/
>> > attachments/20161209/38d13197/attachment.html>
>> > -------------- next part --------------
>> > A non-text attachment was scrubbed...
>> > Name: Screen Shot 2016-12-09 at 22.30.30.png
>> > Type: image/png
>> > Size: 184354 bytes
>> > Desc: Screen Shot 2016-12-09 at 22.30.30.png
>> > URL: <http://lists.geodynamics.org/pipermail/aspect-devel/
>> > attachments/20161209/38d13197/attachment.png>
>> >
>> > ------------------------------
>> >
>> > Subject: Digest Footer
>> >
>> > _______________________________________________
>> > Aspect-devel mailing list
>> > Aspect-devel at geodynamics.org
>> > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>> >
>> > ------------------------------
>> >
>> > End of Aspect-devel Digest, Vol 61, Issue 4
>> > *******************************************
>> >
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachm
>> ents/20161210/5dca7ea5/attachment.html>
>>
>> ------------------------------
>>
>> Subject: Digest Footer
>>
>> _______________________________________________
>> Aspect-devel mailing list
>> Aspect-devel at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>>
>> ------------------------------
>>
>> End of Aspect-devel Digest, Vol 61, Issue 5
>> *******************************************
>>
>
>
> _______________________________________________
> 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/20161210/411efd38/attachment-0001.html>
More information about the Aspect-devel
mailing list