search

UniMainzGeo / LaMEM

LaMEM (Lithosphere and Mantle Evolution Model) - a 3D parallel code to simulate geodynamic and geomechanical processes
MIT License
50 stars 16 forks source link

Phase change problem when using phase diagram #13

Closed NicolasRiel closed 8 months ago

NicolasRiel commented 8 months ago

Hello,

When trying to use phase change with phase diagram ("rho_ph" option in the Material structure) we get the following error message. Without the phase diagram option "rho_ph" everything works fine.

 --------------------------------------------------------------------------
 ============================== INITIAL GUESS =============================
 --------------------------------------------------------------------------
 --------------------------------------------------------------------------
 **************   NONLINEAR SOLVER FAILED TO CONVERGE!   ****************** 
 --------------------------------------------------------------------------
 SNES Divergence Reason  : residual norm is NAN
 [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
 [0]PETSC ERROR: Severe divergence reason detected (see above)
 [0]PETSC ERROR: WARNING! There are option(s) set that were not used! Could be the program crashed before they were used or a spelling mistake, etc!

This error message is obtained for instance using this in the input.dat file:

Untitled-1

Nico

NicolasRiel commented 8 months ago

On the other hand, there is NaN in the phase diagram look-up-table. I might be the guilty one here...

boriskaus commented 8 months ago

Hmm, that is no good. We have no check in LaMEM for NaN‘s…

NicolasRiel commented 8 months ago

The issue has been fixed, it was coming from default "mfmax" not set to 1.0