[aspect-devel] boundary velocity function expressions
John Naliboff
jbnaliboff at ucdavis.edu
Tue Oct 11 10:28:23 PDT 2016
Hi Payman,
A tip for locating files in unix … you can use the “find” command to locate files based on all or part of a file name. The commands below were entered from the main ASPECT folder:
jbnaliboff >> find . -name "lithosphere_boundary_indicator"
./tests/lithosphere_boundary_indicator
This indicates that the relevant file is in the “tests” folder.
The “grep” command is also quite useful for determining what files contain a certain string pattern. For example, if you want to find files that contain the string “two_merged_boxes” in the “source” directory:
jbnaliboff >> grep -r "two_merged_boxes" source/
source//boundary_composition/two_merged_boxes.cc:#include <aspect/boundary_composition/two_merged_boxes.h>
source//boundary_composition/two_merged_boxes.cc:#include <aspect/geometry_model/two_merged_boxes.h>
source//boundary_temperature/two_merged_boxes.cc:#include <aspect/boundary_temperature/two_merged_boxes.h>
source//boundary_temperature/two_merged_boxes.cc:#include <aspect/geometry_model/two_merged_boxes.h>
source//geometry_model/two_merged_boxes.cc:#include <aspect/geometry_model/two_merged_boxes.h>
jbnaliboff >>
The output above shows that the string “two_merged_boxes” appears 5 times in the file two_merged_boxes.cc.
I highly recommend doing some online tutorials on the “find” and “grep” commands, which can make your daily workflow significantly more efficient.
Cheers,
John
*************************************************
John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis
> On Oct 11, 2016, at 8:41 AM, Payman Janbakhsh <payman.janbakhsh at mail.utoronto.ca> wrote:
>
> Found it,
> Many Thanks
>
> payman
>
> From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org] On Behalf Of Glerum, A.C. (Anne)
> Sent: October 11, 2016 11:31 AM
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] boundary velocity function expressions
>
> Hi Payman,
>
> The input files are in the “tests” directory. So they’re only small models, but they should show you how to use to geometry.
>
> Cheers,
>
> Anne
>
> On 11 Oct 2016, at 16:53, Payman Janbakhsh <payman.janbakhsh at mail.utoronto.ca <mailto:payman.janbakhsh at mail.utoronto.ca>> wrote:
>
> Hi Anne,
> Can you point out where “lithosphere_boundary_indicator.prm” or other related examples are located? Couldn’t find it in cookbook or cookbook/future directories.
>
> Thanks
>
> payman
>
> From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org <mailto:aspect-devel-bounces at geodynamics.org>] On Behalf Of Glerum, A.C. (Anne)
> Sent: October 11, 2016 4:14 AM
> To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
> Subject: Re: [aspect-devel] boundary velocity function expressions
>
> Hi Mohamed,
>
> I don't think anything has changed. The plugin file is called two_merged_boxes.cc, but to use the plugin, in the input file select "box with lithosphere boundary indicators".
>
> Payman, there are two tests using this geometry model ("lithosphere_boundary_indicator.prm" for the 2d case), in case you'd like an example.
>
> Cheers,
> Anne
>
>
> On 11 Oct 2016, at 08:32, Mohamed Gouiza <M.Gouiza at leeds.ac.uk <mailto:M.Gouiza at leeds.ac.uk>> wrote:
>
> Hi Anne,
>
> Is the two_merged_boxes geometry the new name for the box_with_lithosphere geometry?
> If it is the case, then the manual needs to be updated
>
> Mohamed
>
> From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org <mailto:aspect-devel-bounces at geodynamics.org>] On Behalf Of John Naliboff
> Sent: 11 October 2016 07:14
> To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
> Subject: Re: [aspect-devel] boundary velocity function expressions
>
> Hi Anne,
>
> Thanks for pointing that out!
>
> Payman - from the file two_merged_boxes.cc documentation section:
> "An example application of this geometry is to prescribe a velocity on "
> "the lithospheric plates, but use open boundary conditions underneath. "
>
> Seems like this is exactly what you are looking for.
>
> Cheers,
> John
>
> *************************************************
> John Naliboff
> Assistant Project Scientist, CIG
> Earth & Planetary Sciences Dept., UC Davis
>
>
>
>
>
>
> On Oct 10, 2016, at 10:37 PM, Glerum, A.C. (Anne) <A.C.Glerum at uu.nl <mailto:A.C.Glerum at uu.nl>> wrote:
>
> Hi Payman and John,
>
> It's possible to prescribe two different velocity boundary conditions on the vertical boundaries of a box domain when you use the geometry model "two merged boxes". For example, a prescribed inflow of the lithosphere and free slip in the mantle.
>
> Cheers,
> Anne
>
>
> On 11 Oct 2016, at 06:55, John Naliboff <jbnaliboff at ucdavis.edu <mailto:jbnaliboff at ucdavis.edu>> wrote:
>
> Hi Payman,
>
> Glad the input file helped.
>
> However, how can I prescribe velocity to a segment of a boundary and let the rest of that boundary be tangential?
>
> In effect, you would like part of a vertical wall (right boundary) to be free-slip, while the other part has the vertical velocity prescribed? Hmm, I’m not sure ASPECT is setup to do a ‘mixed’ velocity boundary. In fact, I’m having a hard time thinking of a study that does this.
>
> In the case you describe below, do you need to prescribed both the vertical velocity and not just the horizontal velocity?
>
> Cheers,
> John
>
> *************************************************
> John Naliboff
> Assistant Project Scientist, CIG
> Earth & Planetary Sciences Dept., UC Davis
>
>
>
>
>
>
> On Oct 10, 2016, at 7:22 PM, Payman Janbakhsh <payman.janbakhsh at mail.utoronto.ca <mailto:payman.janbakhsh at mail.utoronto.ca>> wrote:
>
> Thank you John,
> Actually the continental extension prm file that you sent to Mohamed helped me to figure out the velocity prescription to different boundaries using functions.
>
> However, how can I prescribe velocity to a segment of a boundary and let the rest of that boundary be tangential?
> For instance, considering a right boundary, a mantle lithospheric slab having its own westward prescribed velocity, but the rest of that boundary (mantle) be tangential.
>
> Thanks again
>
> Payman
>
> From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org <mailto:aspect-devel-bounces at geodynamics.org>] On Behalf Of John Naliboff
> Sent: October 10, 2016 12:25 PM
> To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
> Subject: Re: [aspect-devel] boundary velocity function expressions
>
> Hi Payman,
>
> Specifying the velocity on the boundaries works a bit differently than what you have outlined below.
>
> In 2D you need to specify velocity components for x and y, which you do via the function. You can control how the velocity on different boundaries (ex: left & right) by including conditional statements in your functions based on the position (defined by 'x' and 'y' in 2D).
>
> You can also specify that only certain components of the function apply to a certain boundary (i.e. only x velocity from functions is applied to left/right boundary).
>
> If there is a specific setup you are thinking of, can you sent over a basic parameter file that we can look over?
>
> Cheers,
> John
>
>
>
>
> *************************************************
> John Naliboff
> Assistant Project Scientist, CIG
> Earth & Planetary Sciences Dept., UC Davis
> On 10/09/2016 03:08 PM, Payman Janbakhsh wrote:
> Hi,
> I’m trying to implement the following for a 2D problem:
>
> Set Prescribed velocity boundary indicators = left : function , right : function
>
> But get error when trying to give Function expressions in: set Function expression =
> I have tried separating 2 velocity components with semicolon and separate function expressions with comma .
> Example:
> Left boundary function right boundary function
> Set Function expression = if(z>=1000 , 1, 0) ; 0 , if(z>=1000, -1,0) ; 0
>
> I’ve tried many ways including below but all give errors or nonsensical velocity Glyphs:
> Set Function expression = left: if(z>=1000 , 1, 0) ; 0 , right: if(z>=1000, -1,0) ; 0
> Set Function expression = left = if(z>=1000 , 1, 0) ; 0 , right= if(z>=1000, -1,0) ; 0
> Set Function expression = ( if(z>=1000 , 1, 0) ; 0 ) , ( if(z>=1000, -1,0) ; 0 )
> …
>
> How can I separate multiple function expressions each for its corresponding boundary indicator ?
>
> Thanks
>
> Payman
>
>
>
>
>
>
> _______________________________________________
> 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 <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 <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 <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 <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 <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 <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>
> _______________________________________________
> 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/20161011/7f5d31d9/attachment-0001.html>
More information about the Aspect-devel
mailing list