<font color="#330099"><font size="2"><font face="georgia,serif">Sorry, I should have forwarded this to the list !<br></font></font></font><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Louis Moresi</b> <span dir="ltr"><<a href="mailto:louis.moresi@monash.edu">louis.moresi@monash.edu</a>></span><br>
Date: 17 June 2010 13:38<br>Subject: Re: [CIG-LONG] Extension - Compressible Layer in Gale<br>To: "<a href="mailto:atyl0789@uni.sydney.edu.au">atyl0789@uni.sydney.edu.au</a>" <<a href="mailto:atyl0789@uni.sydney.edu.au">atyl0789@uni.sydney.edu.au</a>><br>
<br><br><font color="#330099"><font size="2"><font face="georgia,serif">The question of strain in the "air" layer has been mentioned before on this list and the conclusion is that, since the</font></font></font><div>
<font size="2"></font><font color="#330099" face="georgia, serif">layer is there for numerical convenience, the behaviour is not especially important. If you don't plot this strain / strain rate </font></div>
<div><font color="#330099" face="georgia, serif">then the model looks pretty good. <br></font><div><div><font color="#330099"><font size="2"><font face="georgia,serif"><br></font></font></font></div>
<div><font color="#330099"><font size="2"><font face="georgia,serif">As with any modelling, determining the applicability / suitability of the tool with which you have been provided is</font></font></font></div><div><font color="#330099"><font size="2"><font face="georgia,serif">your task — like a pilot checking the plane before flying it, I suggest you run a number of tests before you make</font></font></font></div>
<div><font color="#330099"><font size="2"><font face="georgia,serif">any bold scientific conclusions. Does the air layer matter ? A more relevant question is how much / when does the air </font></font></font></div><div><font color="#330099"><font size="2"><font face="georgia,serif">influence the model. This is something you can determine empirically. </font></font></font></div>
<div><font color="#330099"><font size="2"><font face="georgia,serif"><br></font></font></font></div><div><font color="#330099"><font size="2"><font face="georgia,serif">I'll leave the question of the hydrostatic correction to Walter or somebody else to deal with.</font></font></font></div>
<div><font color="#330099"><font size="2"><font face="georgia,serif"><br></font></font></font></div><div><font color="#330099"><font size="2"><font face="georgia,serif">L.<br></font></font></font><br><div class="gmail_quote">
<div><div></div><div class="h5">
On 17 June 2010 12:59, <a href="mailto:atyl0789@uni.sydney.edu.au" target="_blank">atyl0789@uni.sydney.edu.au</a> <span dir="ltr"><<a href="mailto:atyl0789@uni.sydney.edu.au" target="_blank">atyl0789@uni.sydney.edu.au</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div class="h5">
G'day all,<br>
<br>
My name is Aaron Tyler and I am currently undertaking my honours year in Geophysics at the University of Sydney, with a project focused on the modelling of tectonic inversion using the Gale software.<br>
<br>
As a first step I have been attempting to create a simple extensional model consisting of 8km of brittle crust overlying 2km of weak viscous material, extending at a constant rate of 1cm/year. I have been so far using the Gale_1_4_1 binary and have encountered a number of problems in getting this model to run. Using either iterative or direct solver, Gale generally crashes after around 50-60 timesteps coinciding with the development of high strain and the first glimpses of faulting. I initially increased resolution to 100m (x and y) in order to determine whether this was the issue, however the same problem was encountered.<br>
<br>
I read that John Naliboff encountered a similar sounding problem in the cig-long archive (March 2010), and following Guillaume and Walter's advice I attempted to implement a compressible air layer above the crust and added a basal velocity. The compressible 'air' layer was 2km thick, had a density of zero, a viscosity of 1e19, and compressibility factor of 10. The model seems to run smoothly, however I observe some extreme deformation of the air layer as extension proceeds, and I am unsure whether this is expected and/or whether it affects the deformation of the layers beneath it. I am also unsure whether this would begin to cause problems if I were to run the model longer. Finally, I am unsure how to use the hydrostatic function with the addition of this compressible layer. Is it possible to make it so that it affects only the lithosphere? I have attached an image of what my model output looks like after 152 timesteps (arbitrary, its just when i stopped the model).<br>
<br>
I have corresponded with John and he has told me that he was unable to get the compressible layer working with the hydrostatic term, and ultimately ended up using an a weak incompressible layer with a density of 1000, simulating a body of water above the crust. Is this a more effective approach for me to follow?<br>
<br>
I am unsure how to proceed at the moment, so any advice on how best to tackle my problems would be fantastic :) I have attached my initial input file without the compressible layer (1012), and an input file with the compressible layer (2001). They are relatively well commented so it should hopefully be clear what I have done. Thanks everyone!<br>
<br>
Cheers,<br>
Aaron<br></div></div>_______________________________________________<br>
CIG-LONG mailing list<br>
<a href="mailto:CIG-LONG@geodynamics.org" target="_blank">CIG-LONG@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long</a><br>
<br></blockquote></div><br></div></div></div>
</div><br>