<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Louis, <br>
<br>
I didn't get completely your message.. <br>
<br>
Is there a way or not in the current version of gale to get Q1P0
and/or Q2P1 working?<br>
Or do you recommend using Q1P0 with UW? <br>
if so, it would be very usefull to provide users with example files
with the source of gale.<br>
<br>
Laetitia <br>
<br>
On 9/24/11 1:28 PM, Louis Moresi wrote:
<blockquote
cite="mid:CACRTXu1Rpa=kG3OqLY574KWyqi9H7+J594w9yug+h-U3MzKJeA@mail.gmail.com"
type="cite">For your information cig-long readers, Walter, Louise
and I had a conversation around this issue two weeks ago.
<div>
<div><br>
</div>
<div>I don't think there are any supporters of Q1Q1 for this
particular problem left any more. Underworld defaults to Q1P0
but provides recovered deviatoric stress and pressure
quantities which mitigate the checkerboard problem. This can
be configured quickly and is at least well-understood and
tested.</div>
<div><br>
</div>
<div>Q2P1 would be useful to have next. It is implemented in
underworld (and therefore gale) but there are some
shortcomings / bugs left to iron out with the current
integration scheme for PIC and (my guess is that) there are
some outstanding issues with multigrid and some of the
processing plugins. I am not sure about the implementation of
the stress / pressure recovery for these elements although
they are likely to be less important.<br>
<div><br>
</div>
<div>Walter was very keen to get this particular element
working in gale before he steps back from active
development.</div>
<div><br>
</div>
<div>Louis</div>
<div><br>
<br>
<div class="gmail_quote">On Sat, Sep 24, 2011 at 20:16,
laetitia le pourhiet <span dir="ltr"><<a
moz-do-not-send="true"
href="mailto:laetitia.le_pourhiet@upmc.fr">laetitia.le_pourhiet@upmc.fr</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div bgcolor="#FFFFFF" text="#000000"> Hi Brian and the
others,<br>
<br>
There is not much you can do about that I gave up
trying to do such models with Gale.<br>
It is now established Q1Q1 elements are not suited for
this kind of problems.<br>
What you can do is to<br>
1) Remove the free surface/sky to work with delta rho
instead of rho<br>
2) Increase the resolution.<br>
<br>
Basically you should only use these elements for
extended thin shit problems. <br>
<br>
The old Q1P0 had other problems, but they could be
used to model drip/subduction problems.<br>
<br>
If CIG still wants to spend some time/money on Gale<br>
It would be usefull to have a switch between elements,
like in UW so that depending on the model set up we
can choose the best suited element.<br>
<br>
hope this was helpfull<br>
<br>
Laetitia
<div>
<div class="h5"><br>
<br>
On 9/24/11 12:44 AM, Brian Wilson wrote: </div>
</div>
<blockquote type="cite">
<div>
<div class="h5">Hi all,<br>
<br>
I'm working on a model of lithospheric drip and
despite having zero flux on the boundaries, I am
getting downward velocities with no compensating
upwelling, which would seem to violate the
divergence free condition. Has anybody else had
a similar problem and/or figured out a solution?
I've attached a copy of my input file if anyone
wants to look at it; everything is scaled as
follows: distance: km, mass: 10^15 kg, time:
m.y.<br>
<br>
Thanks,<br>
Brian<br>
<br>
<fieldset></fieldset>
<br>
</div>
</div>
<div class="im">
<pre>_______________________________________________
CIG-LONG mailing list
<a moz-do-not-send="true" href="mailto:CIG-LONG@geodynamics.org" target="_blank">CIG-LONG@geodynamics.org</a>
<a moz-do-not-send="true" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long</a>
</pre>
</div>
</blockquote>
<font color="#888888"> <br>
<br>
<div>-- <br>
<b>Dr Laetitia Le Pourhiet</b><\n> MCF -
Institut des Sciences de la Terre Paris -
UPMC<\n> tel : <a moz-do-not-send="true"
href="tel:%2B33%201%2044%2027%2058%2083"
value="+33144275883" target="_blank">+33 1 44 27
58 83</a><\n> web : <a
moz-do-not-send="true"
href="http://perso.numericable.com/laetitia.le.pourhiet/index.html"
target="_blank">http://perso.numericable.com/laetitia.le.pourhiet/index.html</a><\n></div>
</font></div>
<br>
_______________________________________________<br>
CIG-LONG mailing list<br>
<a moz-do-not-send="true"
href="mailto:CIG-LONG@geodynamics.org">CIG-LONG@geodynamics.org</a><br>
<a moz-do-not-send="true"
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>
<br clear="all">
<div><br>
</div>
-- <br>
<div><font face="Verdana, Arial, Helvetica, sans-serif"
size="1"><span style="font-size:9px">
<div style="color:rgb(0, 0, 0);font-family:Verdana,
Arial, Helvetica,
sans-serif;font-size:80%;margin-top:8px;margin-right:8px;margin-bottom:8px;margin-left:8px;background-color:rgb(255,
255,
255);padding-top:0.3em;padding-right:0.3em;padding-bottom:0.3em;padding-left:0.3em">
<hr>
<table
style="font-size:medium;border-top-width:1px;border-right-width:1px;border-bottom-width:1px;border-left-width:1px;border-top-style:solid;border-right-style:solid;border-bottom-style:solid;border-left-style:solid;border-top-color:rgb(255,
255, 255);border-right-color:rgb(255, 255,
255);border-bottom-color:rgb(255, 255,
255);border-left-color:rgb(255, 255,
255);width:541px">
<tbody>
<tr>
<td style="color:rgb(0, 0,
0);font-family:Verdana, Arial, Helvetica,
sans-serif;font-size:medium;margin-top:8px;margin-right:8px;margin-bottom:8px;margin-left:8px;border-top-width:0px;border-right-width:0px;border-bottom-width:0px;border-left-width:0px;border-top-style:dashed;border-right-style:dashed;border-bottom-style:dashed;border-left-style:dashed;border-top-color:rgb(187,
187, 187);border-right-color:rgb(187, 187,
187);border-bottom-color:rgb(187, 187,
187);border-left-color:rgb(187, 187,
187);border-style:initial;border-color:initial">
<a moz-do-not-send="true"
href="http://picasaweb.google.com/lh/photo/gQ7_To7dX8sI0sEe2ty0AaO6j9vpN3ZPghWjdy2iDOg?feat=embedwebsite"
target="_blank"><img
moz-do-not-send="true"
src="http://lh4.ggpht.com/_LUlLz13t5HQ/TFJicxcKudI/AAAAAAAABKE/kwovRKlv36Q/s800/Gunnedah3D-LowRes.00001.png"
alt=""
style="border-top-width:0px;border-right-width:0px;border-bottom-width:0px;border-left-width:0px;border-style:initial;border-color:initial"
height="100"></a></td>
<td style="color:rgb(0, 0,
0);font-family:arial,
sans-serif;font-size:11px;margin-top:8px;margin-right:8px;margin-bottom:8px;margin-left:8px;border-top-width:0px;border-right-width:0px;border-bottom-width:0px;border-left-width:0px;border-top-style:dashed;border-right-style:dashed;border-bottom-style:dashed;border-left-style:dashed;border-top-color:rgb(187,
187, 187);border-right-color:rgb(187, 187,
187);border-bottom-color:rgb(187, 187,
187);border-left-color:rgb(187, 187,
187);text-align:right;border-style:initial;border-color:initial">
Professor LouisMoresi<br>
<a moz-do-not-send="true"
href="mailto:louis.moresi@monash.edu"
target="_blank">louis.moresi@monash.edu<br>
</a><a moz-do-not-send="true" title="Go to
my public google profile"
href="http://www.google.com/profiles/louis.moresi"
target="_blank">GoogleProfile</a> | <a
moz-do-not-send="true" title="Link to my
public google calendar"
href="http://www.moresi.info/pages/lm-calendar"
target="_blank">Busy or not ?</a><br>
<a moz-do-not-send="true"
href="http://www.underworldproject.org/"
target="_blank">Underworld Geodynamics
Software</a><br>
Mobile phone: +61 4 2850 1907 ☎</td>
</tr>
</tbody>
</table>
<hr></div>
</span></font></div>
<br>
</div>
</div>
</div>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<b>Dr Laetitia Le Pourhiet</b><\n>
MCF - Institut des Sciences de la Terre Paris - UPMC<\n>
tel : +33 1 44 27 58 83<\n>
web :
<a class="moz-txt-link-freetext" href="http://perso.numericable.com/laetitia.le.pourhiet/index.html">http://perso.numericable.com/laetitia.le.pourhiet/index.html</a><\n></div>
</body>
</html>