<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi C.S.,<br>
<br>
Here is the stress vs distance plot at step 0. The shear stress (xy
component) at x=z=0.5 is plotted. The shear stress decreases toward 0
as y approaches 0 and 1.<br>
<img alt="" src="cid:part1.02070702.09030402@geodynamics.org"
height="452" width="600"><br>
The shear stress is not exactly zero in the plot. The slope of the
curve changes near the boundary, which is caused by interpolation. The
stress is computed on the elements and then interpolated to the nodes.
The interpolation is less accurate at the boundary nodes. <br>
<br>
If you really care about boundary stress, you need to apply additional
correction, like the correction at line 439 in Topo_gravity.c<br>
tpg[snode] = -2 * SZZ[node] + SZZ[node - 1];<br>
Here, tpg is the non-dimensional dynamic topography at surface, which
is the same as -SZZ at the boundary node. The correction is to add
(SZZ[node0] - SZZ[node1]) to SZZ[node0], where node0 is the boundary
node and node1 is the neighboring node in the normal direction, is this
case, right below node0.<br>
<br>
After the correction, the plot becomes:<br>
<img alt="" src="cid:part2.05050809.01000509@geodynamics.org"
height="452" width="600"><br>
Now the boundary stress is indeed very close to zero (1e-3).<br>
<br>
Cheers,<br>
<br>
Eh<br>
<br>
<br>
Shijie Zhong wrote:
<blockquote cite="mid20090109085738.AGM25737@riddler.int.colorado.edu"
type="cite">
<pre wrap="">Dear Mr. Natarajan,
Eh is probably out of town at the moment, and I would be happy to answer your question.
Can you please provide more information? 3D cartesian or regional spherical models? Generally speaking, shear stresses are non-zero numerically at the free-slip boundaries, but if you plot them against the distance from the boundaries, you should see that the stresses decrease towards zero rapidly as the distance decreases.
To reduce the email traffic, you can send your follow-up questions to me and Eh.
Cheers,
Shijie
Shijie Zhong
Department of Physics
University of Colorado at Boulder
Boulder, CO 80309
Tel: 303-735-5095; Fax: 303-492-7935
Web: <a class="moz-txt-link-freetext" href="http://anquetil.colorado.edu/szhong">http://anquetil.colorado.edu/szhong</a>
---- Original message ----
</pre>
<blockquote type="cite">
<pre wrap="">Date: Thu, 8 Jan 2009 23:03:50 -0800 (PST)
From: "C.S Natarajan" <a class="moz-txt-link-rfc2396E" href="mailto:natarajan@geology.ucdavis.edu"><natarajan@geology.ucdavis.edu></a>
Subject: [CIG-MC] Zero shear stresses with free-slip boundary conditions
To: <a class="moz-txt-link-abbreviated" href="mailto:cig-mc@geodynamics.org">cig-mc@geodynamics.org</a>
Hello everybody,
I was wondering if anybody has outputted the shear stresses from
Citcom-CU and actually seen that the free-slip boundary condition
gives 0 shear stresses. I was trying to run the sample input file
provided from the distribution version and output the stresses. I am
getting non-zero shear stresses on the vertical faces. I'm wondering
what I might be doing wrong. The only code I had to add was a routine
to output stresses (calculated in Topo_gravity.c) in Output.c
Thanks in advance,
C.S.N
_______________________________________________
CIG-MC mailing list
<a class="moz-txt-link-abbreviated" href="mailto:CIG-MC@geodynamics.org">CIG-MC@geodynamics.org</a>
<a class="moz-txt-link-freetext" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc</a>
</pre>
</blockquote>
<pre wrap=""><!---->_______________________________________________
CIG-MC mailing list
<a class="moz-txt-link-abbreviated" href="mailto:CIG-MC@geodynamics.org">CIG-MC@geodynamics.org</a>
<a class="moz-txt-link-freetext" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc</a>
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Eh Tan
Staff Scientist
Computational Infrastructure for Geodynamics
2750 E. Washington Blvd. Suite 210
Pasadena, CA 91107
(626) 395-1693
<a class="moz-txt-link-freetext" href="http://www.geodynamics.org">http://www.geodynamics.org</a>
</pre>
</body>
</html>