[CIG-SHORT] 2D Simulations

James Moore james.moore at earth.ox.ac.uk
Wed Sep 14 09:42:11 PDT 2011


I'm trying to approximate an infinitely long strike slip fault. which is
essentially a 2D anti-plane problem. Hence my thoughts regarding using
periodic domains to make the 3D case resemble a stretched 2D case with no
end effects.

 

James

 

 

---------------

D.Phil Research Student,

COMET+ Research Group <http://comet.nerc.ac.uk/> 

 

Department of Earth Sciences,

University of Oxford,

South Parks Road,

Oxford,

OX1 3AN,

United Kingdom

 

From: knepley at gmail.com [mailto:knepley at gmail.com] On Behalf Of Matthew
Knepley
Sent: 14 September 2011 16:56
To: James Moore
Cc: Brad Aagaard; cig-short at geodynamics.org
Subject: Re: [CIG-SHORT] 2D Simulations

 

On Wed, Sep 14, 2011 at 10:50 AM, James Moore <james.moore at earth.ox.ac.uk>
wrote:

Dear Brad (or anyone else who can help me),

Further to my previous question, is it possible to implement periodic
boundary conditions for a 3 D problem, thus negating edge effects?

 

We do not support periodic domains. That seems like a strange setup to me

when your unknowns are displacements. Wouldn't a stress-free condition be

more what you are looking for?

 

  Thanks,

 

     Matt

 

James

---------------
D.Phil Research Student,
COMET+ Research Group

Department of Earth Sciences,
University of Oxford,
South Parks Road,
Oxford,
OX1 3AN,
United Kingdom

-----Original Message-----
From: cig-short-bounces at geodynamics.org
[mailto:cig-short-bounces at geodynamics.org] On Behalf Of Brad Aagaard
Sent: 17 August 2011 18:11
To: cig-short at geodynamics.org
Subject: Re: [CIG-SHORT] 2D Simulations

James-

PyLith does not support 2-D out-of-plane deformation. We have not
implemented this capability because it involves writing extra code to handle
this special case. The workaround is to construct a 3-D model and apply
Dirichlet boundary conditions to enforce the desired symmetry and extract
the solution along the centerline. This is more computational intensive but
works reasonably well.

Brad


On 8/17/11 7:39 AM, James Moore wrote:
> Hello all,
>
> I have recently come back to Pylith after a year or so of not having
> used it and wanted to run some toy problems in 2D. The manual has
> examples of plane strain problems being solved, however I was interested
> in running the anti-plane problem, corresponding to shear out of the
> plane of a 2 dimensional slice through a strike slip fault. I can run
> the 3D equivalent, and tried applying 3D boundary conditions to a 2D
> mesh, but it didn't seem to like it very much.
>
> Any suggestions?
>
> Many thanks,
>
> James Moore
>
> ---------------
>
> D.Phil Research Student,
>
> COMET+ Research Group <http://comet.nerc.ac.uk/>
>
> Department of Earth Sciences,
>
> University of Oxford,
>
> South Parks Road,
>
> Oxford,
>
> OX1 3AN,
>
> United Kingdom
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short

_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org
http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short

_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org
http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short





 

-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-short/attachments/20110914/6e13ad52/attachment-0001.htm 


More information about the CIG-SHORT mailing list