[aspect-devel] boundary velocity function expressions

[John Naliboff]

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
> 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/20161010/d5a0105f/attachment.html>