<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1" style="word-wrap:break-word">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Thanks for your suggestions,<br>
<br>
I tried the pseudo-dirichlet you suggested Brad, but I did not have good results. I think the best way of dealing with this would be to probably be to do transfer the information from the 2D mesh in Pylith to the 3D mesh for the flow simulation. The 3D mesh
would then be generated by extruding the 2D mesh so they have exactly the same in-plane geometry.
<br>
<br>
I haven't written the flow code or the connecting piece myself and I am cooperating with a colleague for that part. I'm wondering if it would be hard to connect the information in the different meshes. Charles, do you have any suggestions on how you have done
this? Is there some pattern in the node numbers you exploited when feeding your 3D pressures into the 2D Pylith mesh? Or did you look at the x-z coordinates for nodes in the 3D mesh and find the corresponding node with same x-z coordinates in the 2D mesh?<br>
<br>
I will discuss this option with my colleague and if it turns out it is not viable I may reopen this thread to find out why the pseudo-dirichlet does not seem to work properly.<br>
<br>
cheers,<br>
Lucas<br>
<br>
<br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF129751"><font color="#000000" face="Tahoma" size="2"><b>From:</b> cig-short-bounces@geodynamics.org [cig-short-bounces@geodynamics.org] on behalf of Charles Williams [willic3@gmail.com]<br>
<b>Sent:</b> Friday, April 05, 2013 18:14<br>
<b>To:</b> cig-short@geodynamics.org<br>
<b>Subject:</b> Re: [CIG-SHORT] pseudo plane-strain in Pylith<br>
</font><br>
</div>
<div></div>
<div>Hi Lucas,
<div><br>
</div>
<div>I've done a similar model, but going the other way (feeding fluid pressures from a flow code into the fault-normal stress for PyLith. To do this, I used a 2D mesh for PyLith and a 3D mesh (1 element thick) for the flow code. If you are careful when designing
the meshes, you can have elements that are identical, except that the quadrilaterals are extruded to form hexahedra. The flow code provided identical fluid pressures on the front and back faces, so I was able to select one to use. Hopefully, you can do something
similar for your problem.</div>
<div><br>
</div>
<div>Cheers,</div>
<div>Charles</div>
<div><br>
</div>
<div><br>
<div>
<div>On 6/04/2013, at 9:07 AM, Lucas Abraham Willemsen wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite"><span class="Apple-style-span" style="border-collapse:separate; font-family:Helvetica; 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; font-size:medium">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Hello,<br>
<br>
Currently I am trying to investigate the geomechanical response of a reservoir with a single fault when CO2 is injected. The fluid-flow part of the problem is performed by a code which is coupled to Pylith. The results from Pylith are fed into the flow code,
which marches time forward and the new fluid pressures are then fed back into Pylith, etc. <span class="Apple-converted-space"> </span><br>
<br>
The CO2 sequestration geometry I am considering is plane-strain, but the complication is that the flow-code (currently experimental) only works with hexahedrals. So I turned the 2D problem into a 3D-problem that is 1-element thick in the out-of-plane direction
(y-direction). To enforce plane strain conditions I put zero displacement in the y-direction on the front and back plane. But the problem is that the fault intersects these front and back planes. This is a natural consequence of turning the 2D geometry with
fault into a 1-element thick 3D geometry with fault.<span class="Apple-converted-space"> </span><br>
<br>
All of the input parameters are independent of 'y' since I am trying to model a plane strain problem. So there should be no reason for slip in the y-direction on any of the fault nodes (intuitively speaking). But Pylith does not allow me to explicitly constrain
the y-displacements of the fault where they intersect the front and back plane. Therefore, the only way to run the model is by excluding the nodes on the fault from the front- and back- plane nodesets. But this way the boundary condition on these fault nodes
is effectively zero normal stress, and elastic extrusion of material takes places. (Applying Neumann B.C. with normal stress to these fault nodes on the front and back plane to somehow reduce the extrusion does not seem to affect the fault nodes either. It
would not be a great solution anyway).<br>
<br>
Is there a way to deal with faults in this 'pseudo plane-strain' geometry I described?<br>
<br>
sincerely,<br>
Lucas<br>
</div>
_______________________________________________<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-bin/mailman/listinfo/cig-short</a></div>
</span></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; 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">
<div style="word-wrap:break-word">
<div>
<div style="margin-top:0px; margin-right:0px; margin-bottom:0px; margin-left:0px">
<font style="font:normal normal normal 12px/normal Helvetica" face="Helvetica" size="3">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" target="_blank">C.Williams@gns.cri.nz</a></div>
</div>
</div>
</span></span></div>
<br>
</div>
</div>
</div>
</div>
</body>
</html>