[aspect-devel] ASCII files - Cut a 2D spherical model in slices
Eleonora Ficini
eleonora.ficini at uniroma1.it
Fri Jan 19 07:37:11 PST 2018
Dear Wolfgang (and dear all),
I'm sorry if it was a bit unspecific :) Trying to explain better what I
would need:
I would like, first of all, to look (in numbers) how viscosity, velocities,
density, temperature etc. evolves with depth in a model.
I am also sorry because, trying to learn ASPECT everyday, a lot of
questions come to my mind (and very often they are really basic).
So, for example, I now have another one. I don't know if it's correct to
write a new question within an old one, but I would like to avoid a
spamming effect for you all! Please, tell me if it's not correct.
The new question is: I implemented GPlates velocities as boundary
conditions for my 2D spherical model.The slice at which my first model was
cut was the geographic equator (the default Point one and Point two section
of the cookbooks...). Since I would like to cut my model using different
coordinates, I tried to change them several times (paying attention to
change them to colatitude and longitude -in radians-), but what I got was
always a cut slice through the geographic equator.
I visualize it with Paraview, using the coastline file.
What am I missing?
Thank you in advance and sorry if I have a lot of questions!
With my best regards,
Eleonora Ficini
2018-01-18 21:00 GMT+01: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: ASPECT Data to ASCII file - aspect-devel Digest, Vol 74,
> Issue 8 (Wolfgang Bangerth)
> 2. Re: Flow through free surface (Lev Karatun)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 17 Jan 2018 15:36:00 -0700
> From: Wolfgang Bangerth <bangerth at colostate.edu>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] ASPECT Data to ASCII file - aspect-devel
> Digest, Vol 74, Issue 8
> Message-ID: <66865122-ea85-1b5b-d114-80f17255314a at colostate.edu>
> Content-Type: text/plain; charset=utf-8; format=flowed
>
>
> > I would like to have an XYZ ascii file to statistically visualize
> > numbers associated with my results and, thus, to have a comprehensive
> > view of them.
>
> That's a bit non-specific :-) Concretely *what* numbers do you need?
> Pressures, velocities, temperatures? And what kind of statistical
> analysis do you want to do with them beyond visualizing them?
>
> Best
> W.
>
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bangerth at colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 17 Jan 2018 20:50:54 -0500
> From: Lev Karatun <lev.karatun at gmail.com>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] Flow through free surface
> Message-ID:
> <CAPGyz275af_WCT1stXicDOMaQcXhTfxCMWKtBGb8M
> MAtg+m+8w at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Wolfgang,
>
> * You want the boundary to be free and move according to whatever flow you
> > get below
> >
> My understanding is that this is the current implementation of free
> surface?
>
> For any more complicated procedure, it would be useful to first define
> > exactly what equations you want to solve, and then (and *only* then)
> think
> > about how you want to solve them.
> >
> What I was thinking of doing is either something as simple as eroding all
> material that is above a specified elevation (in which case I don't even
> need to solve a system of equations -- correct me if I'm wrong please), or
> solving a diffusion equation (which (I think) I know how to solve).
>
> Essentially, right now I know how to adjust the velocity of the surface
> vertices according to the solution of the diffusion equation. But I don't
> know how to superimpose the new boundary on top of the advected material.
> Is it a realistic thing to do? How would I go about it?
>
> Best regards,
> Lev Karatun.
>
> 2018-01-17 16:03 GMT-05:00 Wolfgang Bangerth <bangerth at colostate.edu>:
>
> > On 01/16/2018 01:53 PM, Lev Karatun wrote:
> >
> >>
> >> I was wondering if you could give me an idea of how hard it would be to
> >> allow material flow through the free surface? Right now I'm just able to
> >> dampen the elevation change by applying the diffusion equation, but is
> it
> >> possible to actually take material out of the system? Would the hardest
> >> part be figuring out if a point is inside a curved domain? Or is there
> an
> >> inherent limitation in Aspect that makes it not feasible in principle?
> >> To rephrase what I mean: is it possible to let the material move at the
> >> surface with Stokes velocity, and then draw the new top boundary on top
> of
> >> it, getting rid of material that was eroded away (falls outside of the
> new
> >> model boundaries), and filling in the now less deep trenches with
> crustal
> >> material?
> >>
> >
> > Mathematically, the free boundary is of course determined just so that no
> > flux goes across it.
> >
> > I think what you want is essentially a two-step process:
> > * You want the boundary to be free and move according to whatever flow
> you
> > get below
> > * Then you want to erode away or deposit material on top of it.
> >
> > The easiest way to implement may be exactly to do this: do one; then do
> > the other. For any more complicated procedure, it would be useful to
> first
> > define exactly what equations you want to solve, and then (and *only*
> then)
> > think about how you want to solve them.
> >
> > 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/20180117/39e7d2e4/attachment-0001.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 74, Issue 11
> ********************************************
>
--
Eleonora Ficini, PhD student
Department of Earth Sciences
Sapienza University of Rome
P.le A. Moro, 5 - 00185, Rome, Italy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180119/d37b99ae/attachment.html>
More information about the Aspect-devel
mailing list