search

geodynamics / pylith

PyLith is a finite element code for the solution of dynamic and quasi-static tectonic deformation problems.
Other
153 stars 96 forks source link

Periodic Boundary Implementations in Pylith #73

Open kimtaeho07 opened 5 years ago

kimtaeho07 commented 5 years ago

What is the best way to address periodic boundary conditions in Pylith? It seems that Issue #59 was opened to address the issue and that a direct implementation may still be in development. Any updates?

baagaard-usgs commented 5 years ago

PyLith does not currently include periodic boundary conditions nor is it on our current list of features to add. Periodic boundary conditions is quite different than infinite boundary conditions. Can you describe the context in which you want to use periodic boundary conditions?

kimtaeho07 commented 5 years ago

I plan to create some toy models with simple geometries and meshing to compare to other codes that I have using spectral boundary methods with periodic boundary conditions. Perhaps absorbing boundary conditions will be the best alternative. What do you mean by infinite boundary conditions?

knepley commented 5 years ago

Spectral code have to have periodic conditions, but its not very realistic. Infinite boundaries are also called radiation conditions. It just means that nothing reflects from the boundary.