Hello all,
I successfully installed Pylith 2.22. But I had a little problem with the simulation .
1、 I want to simulate the deformation of the interior of the lithosphere under gravity. I take the 2D model, and I take the quadrilateral mesh.
2、 I set up two lithospheres, plateau 、transition zone and basin, and set their brittle layer as elastic material and ductile layer as viscoelastic material through rheological strength.
3、 For the transition zone, the density and viscosity are given as a linear vary along the transverse direction and depth. For plateaus and basins, only vary in depth are considered .
4、 No faults are defined within the lithosphere.
5、 I set the sides and bottom of the model to free slip boundaries .
6、 In order to test the correctness of the simulation, I have set the mesh to be rough. In previous simulations, it might have taken only a few minutes to compute, but now it takes so long to execute that it seems impossible to converge .
The petsc message is as follows ：
[pylithapp.petsc]
Preconditioner settings.
pc_type = asm
Convergence parameters.
ksp_rtol = 1.0e-9
ksp_atol = 1.0e-9
Linear solver monitoring options.
ksp_monitor = true
#ksp_view = true
ksp_converged_reason = true
ksp_error_if_not_converged = true
Nonlinear solver monitoring options.
snes_rtol = 1.0e-9
snes_atol = 1.0e-9
snes_max_it = 500
snes_monitor = true
#snes_view = true
snes_converged_reason = true
snes_linesearch_monitor = true
snes_error_if_not_converged = true
The error message is as follows ：
/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/problems/Implicit.py:207:step
– implicit(info)
– Solving equations.
0 SNES Function norm 3.504722374800e-09
0 KSP Residual norm 9.999723312736e-10
Linear solve converged due to CONVERGED_ATOL iterations 0
Line search: Initial direction and size is 0
Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 0
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR:
[0]PETSC ERROR: SNESSolve has not converged
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.10.2, Jul, 01, 2019
[0]PETSC ERROR: /home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/bin/mpinemesis on a named ubuntu by wuhang Sun Jun 18 20:33:56 2023
[0]PETSC ERROR: Configure options --prefix=/home/brad/pylith-binary/dist --with-c2html=0 --with-x=0 --with-clanguage=C --with-mpicompilers=1 --with-shared-libraries=1 --with-64-bit-points=1 --with-large-file-io=1 --download-chaco=1 --download-ml=1 --download-f2cblaslapack=1 --with-hwloc=0 --with-ssl=0 --with-x=0 --with-c2html=0 --with-lgrind=0 --with-hdf5=1 --with-hdf5-dir=/home/brad/pylith-binary/dist --with-zlib=1 --LIBS=-lz --with-debugging=0 --with-fc=0 CPPFLAGS=“-I/home/brad/pylith-binary/dist/include -I/home/brad/pylith-binary/dist/include " LDFLAGS=”-L/home/brad/pylith-binary/dist/lib -L/home/brad/pylith-binary/dist/lib64 -L/home/brad/pylith-binary/dist/lib -L/home/brad/pylith-binary/dist/lib64 " CFLAGS=“-g -O2” CXXFLAGS=“-g -O2 -DMPICH_IGNORE_CXX_SEEK” FCFLAGS= PETSC_DIR=/home/brad/pylith-binary/build/petsc-pylith PETSC_ARCH=arch-pylith
[0]PETSC ERROR: #1 SNESSolve() line 4408 in /home/brad/pylith-binary/build/petsc-pylith/src/snes/interface/snes.c
[0]PETSC ERROR: #2 void pylith::problems::SolverNonlinear::solve(pylith::topology::Field*, pylith::topology::Jacobian*, const pylith::topology::Field&)() line 152 in …/…/…/pylith-2.2.2/libsrc/pylith/problems/SolverNonlinear.cc
Fatal error. Calling MPI_Abort() to abort PyLith application.
Traceback (most recent call last):
File “/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/apps/PetscApplication.py”, line 74, in onComputeNodes
self.main(*args, **kwds)
File “/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/apps/PyLithApp.py”, line 138, in main
self.problem.run(self)
File “/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/problems/TimeDependent.py”, line 203, in run
self.formulation.step(t, dt)
File “/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/problems/Implicit.py”, line 212, in step
self.solver.solve(dispIncr, self.jacobian, residual)
File “/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/problems/problems.py”, line 186, in solve
def solve(self, *args): return _problems.SolverNonlinear_solve(self, *args)
RuntimeError: Error detected while in PETSc function.
application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/bin/nemesis: mpirun: exit 255
/home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/bin/pylith: /home/wuhang/Desktop/pylith/pylith/pylith-2.2.2-linux-x86_64/bin/nemesis: exit 1
Does anyone know what caused this and how I can fix it?
Thank you
gravity_infstrain.cfg (14.7 KB)
pylithapp.cfg (24.0 KB)