<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Keep in mind that Dirichlet BC would not exactly be the same as a volumetric expansion or contraction, since you are fixing the absolute displacement in the interior of your domain.<div><br></div><div>Depending on the geometry of your problem, you could still use Neumann BC, or you might even consider having material in the reservoir and applying an initial strain.</div><div><br></div><div>Cheers,</div><div>Charles</div><div><br></div><div><br><div><div>On 14/07/2013, at 4:19 AM, Matthew Knepley wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">On Sat, Jul 13, 2013 at 11:17 AM, Brad Aagaard <span dir="ltr"><<a href="mailto:baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Scott,<br>
<br>
Specifying a Dirichlet boundary condition using an orientation other than the coordinate axes is not implemented. I think you are the first person to request this feature. We do have a DirichletBoundary setup that associates a boundary mesh with a Dirichlet boundary condition so the orientation information could be setup. We will add this to a list of features to consider adding in a future release (right now we have several features to add that are higher priority).<br>
</blockquote><div><br></div><div style="">Right now, since you know the normal, you could specify a Cartesian Dirichlet BC as a stopgap.</div><div style=""><br></div><div style=""> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Regards,<br>
Brad<div class="HOEnZb"><div class="h5"><br>
<br>
On 7/13/13 9:01 AM, Scott Henderson wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I'm wondering if it's possible to specify a 'normal displacement' boundary condition along a curved surface in PyLith. For example 'displacement-normal', similar to 'traction-normal' for the Neumann condition.<br>
<br>
This would be useful to represent fixed volumetric expansion or contraction of a subsurface reservoir (e.g. hydrothermal or magmatic).<br>
<br>
Thanks,<br>
Scott<br>
______________________________<u></u>_________________<br>
CIG-SHORT mailing list<br>
<a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-<u></u>bin/mailman/listinfo/cig-short</a><br>
<br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
CIG-SHORT mailing list<br>
<a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-<u></u>bin/mailman/listinfo/cig-short</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>
_______________________________________________<br>CIG-SHORT mailing list<br><a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" style="font: normal normal normal 12px/normal Helvetica; ">Charles A. Williams</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Scientist<br>GNS Science<br>1 Fairway Drive, Avalon<br>PO Box 30368<br>Lower Hutt 5040<br>New Zealand<br>ph (office): 0064-4570-4566<br>fax (office): 0064-4570-4600</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="mailto:C.Williams@gns.cri.nz">C.Williams@gns.cri.nz</a></div></div></div></span></span>
</div>
<br></div></body></html>