[CIG-SHORT] pylith info

[Brad Aagaard]

Luigi-

You can perform the calculation in one simulation as long as you can 
setup your boundary conditions for the two calculations simultaneously.
If you have zero displacement boundary conditions during your 
gravitational calculation and start with zero displacement on the 
boundaries in your quasi-static problem, them you can combine them into 
one calculation. If you want to start with nonzero displacements on the 
boundaries in your quasi-static problem, them you would need to separate 
the problem into two simulations.

Right now, we have not setup PyLith to easily use the output from one 
simulation as an initial condition or input in another simulation. This 
is something we plan to do in the next year. The current workaround is 
to create a spatial database from the output of the first simulation to 
use as input in the second simulation.

If the strain is concentrated near the boundary, make sure the solution 
is converging. In your .cfg file, use

[pylithapp.petsc]
ksp_monitor = true
ksp_view = true
ksp_converged_reason = true

and look at the stdout to see whether the solver reports that the 
solution converged or diverged.

See examples/3d/hex8/pylithapp.cfg for an example. You may need to 
increase the maximum number of iterations (ksp_max_it).

Brad

On 03/25/2011 10:16 AM, luigi vadacca wrote:
> Dear all,
> I calculated the gravity of a block of dimensions 130 km x 130 km x 20 km (with
> elastic rheology). Now I want to use the values obtained as an initial condition
> for a second model. The latter will consider static friction with quasi-static
> dirichlet boundary conditions. Is this possible?
> Otherwise, can I perform the calculation in one step? Thanks in advance.
>
> Regards
>
> Luigi
>
> p.s: I have already tried to calculate a model with static Friction and
> quasi-static dirichlet boundary conditions , but the strain remained
> concentrated only near the boundary of the block.
>
> p.s.2: fault info:
> min depth: 0 km
> max depth: 10 km
> strike (deg): 140
> dip (deg): 25
> rake: 270
>
>
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short