unimainzgeo / lamem Goto Github PK

Phase Transitions don't work if the involved phases have phase==0.
This can be reproduced with t16_PhaseTransitions/Plume_PhaseTransitions.dat by making the air phase 3 instead of 0, and the mantle phase 0 instead of 3.

This results in:

whereas it should look like:

A test file is here: [Plume_PhaseTransitions_SwappedPhases.dat.zip]
(https://github.com/UniMainzGeo/LaMEM/files/14157744/Plume_PhaseTransitions_SwappedPhases.dat.zip)

Hello!

I'm trying to use LaMEM to run a 3D model which requires the existence of a plume boundary. I've used these boundaries before with success, but after updating the code recently I've been met with this error:

[124]PETSC ERROR: ------------------------------------------------------------------------
[124]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[124]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[124]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind and https://petsc.org/release/faq/
[124]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[124]PETSC ERROR: to get more information on the crash.
[124]PETSC ERROR: Run with -malloc_debug to check if memory corruption is causing the crash.
--------------------------------------------------------------------------

The plume inflow was defined like this:

# Set a plume in the plateau region:
        Plume_InflowBoundary 	= 	1	
        Plume_Type                     =      Inflow_Type
	Plume_Dimension           =	3D	
	Plume_Phase 		        =	7	
	Plume_Depth		        = 	0
	Plume_Temperature 	=	1700	
	Plume_Inflow_Velocity 	= 	15.0
	Plume_VelocityType        =	Gaussian
	Plume_Center 		        =	-200.0 -100.0
	Plume_Radius 		        =	50		

The only change from my previously used plumes are the coordinates and the number of the inflow phase. It was attempting to run on 128 cores.

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:

Nico

Hi,

I'm using LaMEM to run a 2D subduction model with phase transition from basalt to eclogite. I set an internal velocity 5cm/yr on the horizontal plate to initializing subduction. The viscosity field seems normal, but the pressure field looks very strange, very different from the example in LaMEM. Besides, I use the eclogite phase transition in LaMEM. However, even the temeprature of basalt parts is higher than 1000 degree, there is no phase transition. I use 32 cpu cores.

I use 20 SNES iterations, and the SNES residual is ~100.

The following are Petsc options and result pictures.

`
-snes_npicard 10
-snes_monitor
-snes_atol 1e-4
-snes_rtol 1e-5
-snes_stol 1e-16
-snes_max_it 20
-snes_max_funcs 50000
-snes_max_linear_solve_fail 10000
-snes_type newtonls

-snes_ksp_ew
-snes_ksp_ew_rtolmax 1e-2

-snes_PicardSwitchToNewton_rtol 1e-4
-snes_NewtonSwitchToPicard_it 15
-snes_NewtonSwitchToPicard_rtol 1e-3

-snes_linesearch_monitor
-snes_linesearch_type cp
-sens_linesearch_maxstep 1.0

-js_ksp_type fgmres
-js_ksp_max_it 25
-js_ksp_converged_reason
-js_ksp_monitor
-js_ksp_rtol 1e-4
-js_ksp_atol 1e-4`

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!

For instance: (I cannot copy the code, it put everything as a block)

Nico

I am currently working with LaMEM and encountered an issue while compiling using the "make mode=opt all" command. The error message I received is as follows:

In file paraViewOutBin.cpp, line 378:
Error: macro "SETERRQ" is passed 4 arguments but takes only 3
  378 |   SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_USER, "Too many phase aggregates specified! Max allowed: %lld", (LLD)_max_num_phase_agg_);

Could you kindly provide guidance on how to resolve this compilation error? I appreciate your assistance in this matter.