[CIG-LONG] CIG-LONG Digest, Vol 43, Issue 14

[Karen Paczkowski]

Hi Walter,

I tried adding a slab velocity into the program you sent me and have  
been having issues.  I have attached an xml file with the simplest  
case I tried.  It is just the hydrostatically balanced box you sent  
me, but instead of vx = vy = 0 on the left hand side, it has vx=0 and  
vy = -1.  I expected the down going velocity on the left side would  
cause mass to be pulled into the simulation box on the right hand side  
(like in a cornerflow problem), but instead I just get a velocity  
field that is straight down everywhere.  I feel like it may be  
creating mass at the top boundary, but I'm not sure.  If I set vy=0 on  
the top of the box I get the results I expect, but I would prefer to  
not do that since I would like to track the finite deformation of the  
top interface.

Do you know how I can fix this?

Thanks,
Karen

-------------- next part --------------
A non-text attachment was scrubbed...
Name: Slab.xml
Type: text/xml
Size: 20063 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20100617/c6af6c4d/attachment-0001.bin 
-------------- next part --------------




On Jun 17, 2010, at 3:36 PM, cig-long-request at geodynamics.org wrote:

> Send CIG-LONG mailing list submissions to
> 	cig-long at geodynamics.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> 	http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
> or, via email, send a message with subject or body 'help' to
> 	cig-long-request at geodynamics.org
>
> You can reach the person managing the list at
> 	cig-long-owner at geodynamics.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of CIG-LONG digest..."
>
>
> Today's Topics:
>
>   1. Re: Brittle and Ductile deformation (Walter Landry)
>   2. Fw: Modifed HydrostaticTerm.c => VTK pressure output	problem
>      (Walter Landry)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 17 Jun 2010 12:28:19 -0700 (PDT)
> From: Walter Landry <walter at geodynamics.org>
> Subject: Re: [CIG-LONG] Brittle and Ductile deformation
> To: misl5790 at uni.sydney.edu.au
> Cc: cig-long at geodynamics.org
> Message-ID:
> 	<20100617.122819.506052730756394805.walter at geodynamics.org>
> Content-Type: Text/Plain; charset=us-ascii
>
> "Md. Tariqul Islam" <misl5790 at uni.sydney.edu.au> wrote:
>> Hi Walter,
>> Thanks for your previous help. Now I would like to share some  
>> issues with you.
>> One is the brittle deformation. As you mentioned in the user manual  
>> that if
>> the crust has both yielding and viscus rheology it should deform  
>> half the
>> crust ; but for my case it deformed until Moho and the surface did  
>> not
>> subsided enough indicating it could be totally brittle ( attached a  
>> strain
>> rate invariant image).
>
> I do not quite see the problem.  Were you expecting the faults to not
> reach down to the Moho?  You would have to fiddle with the parameters
> a bit to get that to work.
>
>> The other issue ( this is different model) is the change of  
>> viscosity - the
>> mantle is showing high viscosity though I used low viscosity ( the  
>> input file
>> and image for viscosity field is attached).
>> As I am using non Newtonian viscosity it is very sensitive, can you  
>> suggest me
>> how it can be consistent as it is crushing around thirty steps.
>
> You have set the maxViscosity of the mantle to 1.  My guess is that
> your coefficients are off.  For example, your T_0 looks rather high
> (62545).  I posted an input file with a working NonNewtonian Rheology
> at
>
>  http://www.geodynamics.org/pipermail/cig-long/2010-May/000376.html
>
> It is for the as yet unreleased version of Gale, but you should be
> able to use the same numbers for scaling.
>
>> The other thing is I am using the direct solver, do you think that
>> this is making any difference?
>
> A direct solver works better than anything else, so you do not need to
> change that.
>
> Cheers,
> Walter Landry
> walter at geodynamics.org
>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 17 Jun 2010 12:34:48 -0700 (PDT)
> From: Walter Landry <walter at geodynamics.org>
> Subject: [CIG-LONG] Fw: Modifed HydrostaticTerm.c => VTK pressure
> 	output	problem
> To: cig-long at geodynamics.org
> Message-ID: <20100617.123448.85384807014128993.walter at geodynamics.org>
> Content-Type: text/plain; charset="us-ascii"
>
> Forwarded.
>
> Walter Landry
> walter at geodynamics.org
>
> -------------- next part --------------
> An embedded message was scrubbed...
> From: John Naliboff <jbnaliboff at ucdavis.edu>
> Subject: Modifed HydrostaticTerm.c => VTK pressure output problem
> Date: Thu, 10 Jun 2010 16:37:57 -0700
> Size: 512616
> Url: http://geodynamics.org/pipermail/cig-long/attachments/20100617/372e7fbf/attachment.eml
>
> ------------------------------
>
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>
>
> End of CIG-LONG Digest, Vol 43, Issue 14
> ****************************************