[CIG-SHORT] Applying Gravitational Forces with Pylith

[Niloufar Abolfathian]

Hi,

I followed the 2D gravity example and I could run my models with gravity
and add initial stresses. But in models including finite strain, the
nonlinear solver could not converge. I have tried running the model with
gravitation on and off, applying depth dependent initial stress and no
initial stress, applying multigrid preconditioner and also not applying the
preconditioner and also modifying the linear and nonlinear tolerances, but
could not converge the model. I have attached the setting when I am using
the gravity on, applying initial stress and using preconditioning, I am
wondering why it is always giving the error after 6000 years. I really
appreciate all your help solving these issues.

Thanks,
Niloufar



On Tue, Aug 14, 2018 at 8:55 PM, Niloufar Abolfathian <
niloufar.abolfathian at gmail.com> wrote:

> Hi,
>
> I followed the 2D gravity example and I could run my models with gravity
> and add initial stresses. But in models including finite strain, the
> nonlinear solver could not converge. I have tried running the model with
> gravitation on and off, applying depth dependent initial stress and no
> initial stress, applying multigrid preconditioner and also not applying the
> preconditioner and also modifying the linear and nonlinear tolerances, but
> could not converge the model. I have attached the setting when setting when
> I am using the gravity on, applying initial stress and using
> preconditioning, I am wondering why it is always giving the error after
> 6000 years. I really appreciate all your help solving these issues.
>
> Thanks,
> Niloufar
>
>
>
> On Thu, Aug 2, 2018 at 8:43 AM, Brad Aagaard <baagaard at usgs.gov> wrote:
>
>> Is this the gravity only problem to compute the initial stresses? If so,
>> then you should remove the fault.
>>
>> On 08/01/2018 08:59 PM, Niloufar Abolfathian wrote:
>>
>>> Thanks, Charles.
>>>
>>> My error is due to the convergence of the nonlinear solver. I have made
>>> some modifications but not sure if it is a good idea to decrease my
>>> residual norms.
>>>
>>> Best,
>>> Niloufar
>>>
>>>
>>> This is the setting I am using:
>>>
>>>
>>> [pylithapp.timedependent]
>>> formulation =pylith.problems.ImplicitLgDeform
>>>
>>>
>>> [pylithapp.timedependent.formulation]
>>> output =[domain,subdomain]
>>> output.subdomain =pylith.meshio.OutputSolnSubset
>>> solver =pylith.problems.SolverNonlinear
>>>
>>> [pylithapp.timedependent.formulation.time_step]
>>> total_time =50000.0*year
>>> dt =2.0*year
>>>
>>> [pylithapp.timedependent.interfaces]
>>> fault =pylith.faults.FaultCohesiveDyn
>>>
>>> [pylithapp.timedependent.interfaces.fault]
>>> label =fault
>>> edge =fault_edge
>>> zero_tolerance =1.0e-11
>>>
>>> # Use the slip-weakening friction model.
>>> friction.db_properties =spatialdata.spatialdb.UniformDB
>>> friction.db_properties.label =Slip weakening
>>> friction.db_properties.values =[static-coefficient,dynamic-c
>>> oefficient,slip-weakening-parameter,cohesion]
>>>
>>> friction.db_properties.data =[0.6,0.4,0.2*m,0.0*Pa]
>>>
>>> # ----------------------------------------------------------------------
>>> # PETSc settings
>>> # ----------------------------------------------------------------------
>>> [pylithapp.petsc]
>>> malloc_dump =
>>>
>>> # Preconditioner settings.
>>> pc_type =asm
>>> sub_pc_factor_shift_type =nonzero
>>>
>>
>> When you encounter convergence errors, a good strategy is to revert to
>> the direct solver. Note that this only works when you don't have a fault.
>>
>> pc_type = lu
>>
>>
>> For most real applications (after you have debugged the setup)l, we
>> recommend using the multigrid preconditioner. See
>> solver_algebraicmultigrid.cfg (no fault) and solver_fieldsplit.cfg )with
>> fault) in examples/3d/subduction.
>>
>>
>>> # Friction sensitivity solve used to compute the increment in slip
>>> # associated with changes in the Lagrange multiplier imposed by the
>>> # fault constitutive model.
>>> friction_pc_type =asm
>>> friction_sub_pc_factor_shift_type =nonzero
>>> friction_ksp_max_it =50
>>> friction_ksp_gmres_restart =30
>>>
>>> # Uncomment to view details of friction sensitivity solve.
>>> #friction_ksp_monitor = true
>>> #friction_ksp_view = true
>>> friction_ksp_converged_reason =true
>>>
>>> # Reduce convergence tolerances.
>>> ksp_rtol =1.0e-16
>>> ksp_atol =1.0e-12
>>> ksp_max_it =4000
>>> ksp_gmres_restart =100
>>>
>>> snes_rtol =1.0e-9
>>> snes_atol =1.0e-8
>>> snes_max_it =4000
>>>
>>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180814/59637bcd/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: m1st24gb.cfg
Type: application/octet-stream
Size: 11111 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180814/59637bcd/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pylithapp.cfg
Type: application/octet-stream
Size: 3494 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180814/59637bcd/attachment-0003.obj>